R. M. Law

Towards robust regional estimates of CO2 sources and sinks using atmospheric transport models

Information about regional carbon sources and sinks can be derived from variations in observed atmospheric CO2 concentrations via inverse modelling with atmospheric tracer transport models. A consensus has not yet been reached regarding the sizexa0and distribution of regional carbon fluxes obtained using this approach, partly owing to the use of several different atmospheric transport models. Here we report estimates of surface–atmosphere CO2 fluxes from an intercomparison of atmospheric CO2 inversion models (the TransCom 3 project), which includes 16 transport models and model variants. We find an uptake of CO2 in the southern extratropical ocean less than that estimated from ocean measurements, a result that is not sensitive to transport models or methodological approaches. We also find a northern land carbon sink that is distributed relatively evenly among the continents of the Northern Hemisphere, but these results show some sensitivity to transport differences among models, especially in how they respond to seasonal terrestrial exchange of CO2. Overall, carbon fluxes integrated over latitudinal zones are strongly constrained by observations in the middle to high latitudes. Further significant constraints to our understanding of regional carbon fluxes will therefore require improvements in transport models and expansion of the CO2 observation network within the tropics.

TransCom 3 inversion intercomparison: Impact of transport model errors on the interannual variability of regional CO2 fluxes, 1988–2003

Monthly CO2 fluxes are estimated across 1988–2003 for 22 emission regions using data from 78 CO2 measurement sites. The same inversion (method, priors, data) is performed with 13 different atmospheric transport models, and the spread in the results is taken as a measure of transport model error. Interannual variability (IAV) in the winds is not modeled, so any IAV in the measurements is attributed to IAV in the fluxes. When both this transport error and the random estimation errors are considered, the flux IAV obtained is statistically significant at P ≤ 0.05 when the fluxes are grouped into land and ocean components for three broad latitude bands, but is much less so when grouped into continents and basins. The transport errors have the largest impact in the extratropical northern latitudes. A third of the 22 emission regions have significant IAV, including the Tropical East Pacific (with physically plausible uptake/release across the 1997–2000 El Nino/La Nina) and Tropical Asia (with strong release in 1997/1998 coinciding with large-scale fires there). Most of the global IAV is attributed robustly to the tropical/southern land biosphere, including both the large release during the 1997/1998 El Nino and the post-Pinatubo uptake.

Transcom 3 inversion intercomparison: Model mean results for the estimation of seasonal carbon sources and sinks

[1] The TransCom 3 experiment was begun to explore the estimation of carbon sources and sinks via the inversion of simulated tracer transport. We build upon previous TransCom work by presenting the seasonal inverse results which provide estimates of carbon flux for 11 land and 11 ocean regions using 12 atmospheric transport models. The monthly fluxes represent the mean seasonal cycle for the 1992 to 1996 time period. The spread among the model results is larger than the average of their estimated flux uncertainty in the northern extratropics and vice versa in the tropical regions. In the northern land regions, the model spread is largest during the growing season. Compared to a seasonally balanced biosphere prior flux generated by the CASA model, we find significant changes to the carbon exchange in the European region with greater growing season net uptake which persists into the fall months. Both Boreal North America and Boreal Asia show lessened net uptake at the onset of the growing season with Boreal Asia also exhibiting greater peak growing season net uptake. Temperate Asia shows a dramatic springward shift in the peak timing of growing season net uptake relative to the neutral CASA flux while Temperate North America exhibits a broad flattening of the seasonal cycle. In most of the ocean regions, the inverse fluxes exhibit much greater seasonality than that implied by the DpCO2 derived fluxes though this may be due, in part, to misallocation of adjacent land flux. In the Southern Ocean, the austral spring and fall exhibits much less carbon uptake than implied by DpCO2 derived fluxes. Sensitivity testing indicates that the inverse estimates are not overly influenced by the prior flux choices. Considerable agreement exists between the model mean, annual mean results of this study and that of the previously published TransCom annual mean inversion. The differences that do exist are in poorly constrained regions and tend to exhibit compensatory fluxes in order to match the global mass constraint. The differences between the estimated fluxes and the prior model over the northern land regions could be due to the prior model respiration response to temperature. Significant phase differences, such as that in the Temperate Asia region, may be due to the limited observations for that region. Finally, differences in the boreal land regions between the prior model and the estimated fluxes may be a reflection of the timing of spring thaw and an imbalance in respiration versus photosynthesis. INDEX TERMS: 0322 Atmospheric Composition and Structure: Constituent sources and sinks; 1615 Global Change: Biogeochemical processes (4805); 0315 Atmospheric Composition and Structure: Biosphere/atmosphere interactions; KEYWORDS: carbon transport, inversion

TransCom model simulations of hourly atmospheric CO2: Analysis of synoptic-scale variations for the period 2002-2003

The ability to reliably estimate CO2 fluxes from current in situ atmospheric CO2 measurements and future satellite CO2 measurements is dependent on transport model performance at synoptic and shorter timescales. The TransCom continuous experiment was designed to evaluate the performance of forward transport model simulations at hourly, daily, and synoptic timescales, and we focus on the latter two in this paper. Twenty-five transport models or model variants submitted hourly time series of nine predetermined tracers (seven for CO2) at 280 locations. We extracted synoptic-scale variability from daily averaged CO2 time series using a digital filter and analyzed the results by comparing them to atmospheric measurements at 35 locations. The correlations between modeled and observed synoptic CO2 variabilities were almost always largest with zero time lag and statistically significant for most models and most locations. Generally, the model results using diurnally varying land fluxes were closer to the observations compared to those obtained using monthly mean or daily average fluxes, and winter was often better simulated than summer. Model results at higher spatial resolution compared better with observations, mostly because these models were able to sample closer to the measurement site location. The amplitude and correlation of model-data variability is strongly model and season dependent. Overall similarity in modeled synoptic CO2 variability suggests that the first-order transport mechanisms are fairly well parameterized in the models, and no clear distinction was found between the meteorological analyses in capturing the synoptic-scale dynamics.

Precision requirements for space-based XCO 2 data

Precision requirements are determined for space-based column-averaged CO_2 dry air mole fraction (X_(CO)_2) data. These requirements result from an assessment of spatial and temporal gradients in (X_(CO)_2) the relationship between (X_(CO)_2) precision and surface CO_2 flux uncertainties inferred from inversions of the (X_(CO)_2) data, and the effects of (X_(CO)_2) biases on the fidelity of CO_2 flux inversions. Observational system simulation experiments and synthesis inversion modeling demonstrate that the Orbiting Carbon Observatory mission design and sampling strategy provide the means to achieve these (X_(CO)_2) data precision requirements.

TransCom model simulations of hourly atmospheric CO2 : experimental overview and diurnal cycle results for 2002

[1] A forward atmospheric transport modeling experiment has been coordinated by the TransCom group to investigate synoptic and diurnal variations in CO2. Model simulations were run for biospheric, fossil, and air-sea exchange of CO2 and for SF6 and radon for 2000-2003. Twenty-five models or model variants participated in the comparison. Hourly concentration time series were submitted for 280 sites along with vertical profiles, fluxes, and meteorological variables at 100 sites. The submitted results have been analyzed for diurnal variations and are compared with observed CO2 in 2002. Mean summer diurnal cycles vary widely in amplitude across models. The choice of sampling location and model level account for part of the spread suggesting that representation errors in these types of models are potentially large. Despite the model spread, most models simulate the relative variation in diurnal amplitude between sites reasonably well. The modeled diurnal amplitude only shows a weak relationship with vertical resolution across models; differences in near-surface transport simulation appear to play a major role. Examples are also presented where there is evidence that the models show useful skill in simulating seasonal and synoptic changes in diurnal amplitude.

Interannual variability of the global carbon cycle (1992–2005) inferred by inversion of atmospheric CO2 and δ13CO2 measurements

[1]xa0We present estimates of the surface sources and sinks of CO2 for 1992–2005 deduced from atmospheric inversions. We use atmospheric CO2 records from 67 sites and 10 δ13CO2 records. We use two atmospheric models to increase the robustness of the results. The results suggest that interannual variability is dominated by the tropical land. Statistically significant variability in the tropical Pacific supports recent ocean modeling studies in that region. The northern land also shows significant variability. In particular, there is a large positive anomaly in 2003 in north Asia, which we associate with anomalous biomass burning. Results using δ13CO2 and CO2 are statistically consistent with those using only CO2, suggesting that it is valid to use both types of data together. An objective analysis of residuals suggests that our treatment of uncertainties in CO2 is conservative, while those for δ13CO2 are optimistic, highlighting problems in our simple isotope model. Finally, δ13CO2 measurements offer a good constraint to nearby land regions, suggesting an ongoing value in these measurements for studies of interannual variability.

The relationship between tropical CO2 fluxes and the El Niño‐Southern Oscillation

This paper summarizes some features of the interannual variability of tropical CO2 sources during 1980–95. Sources are derived from inversion of atmospheric concentration and isotopic data using three different techniques and two different transport models. We show that the tropical source is significantly correlated with the SOI. Composite CO2 sources for ENSO events show an initial negative anomaly followed by a positive anomaly. We tentatively attribute the negative anomaly to the ocean and the positive anomaly to a terrestrial response.

Comment on “Saturation of the Southern Ocean CO2 Sink Due to Recent Climate Change”

Unlike Le Quéré et al. (Reports, 22 June 2007, p. 1735), we do not find a saturating Southern Ocean carbon sink due to recent climate change. In our ocean model, observed wind forcing causes reduced carbon uptake, but heat and freshwater flux forcing cause increased uptake. Our inversions of atmospheric carbon dioxide show that the Southern Ocean sink trend is dependent on network choice.

Global observations of the carbon budget 3. Initial assessment of the impact of satellite orbit, scan geometry, and cloud on measuring CO2 from space

[1]xa0In this third of three companion papers we assess the utility of column-integrated measurements of CO2 in constraining surface source estimates. The measurements have error characteristics defined by OBrien and Rayner [2002]. There we described a retrieval algorithm which appeared capable of unbiased retrievals. An assessment of the sensitivity to dominant sources of error in these measurements suggested a precision about 0.5% of the background mixing ratio. We include this information along with the probability of clear sky and the sampling pattern imposed by a particular orbit and scan geometry. We carry out a synthesis inversion to recover specified surface sources. The uncertainty of the inferred sources quantifies the strength of the constraint offered by such measurements. The difference between input and retrieved sources demonstrates biases in the whole procedure. We show that source uncertainties are low in the presence of a sunlit surface but rise dramatically at high latitudes in winter. Provided we use sufficient spatial resolution for the sources we estimate, the inversion can also produce relatively accurate estimates. The inversion procedure is sensitive to biases caused by poor sampling of the diurnal cycle. Errors that isotropically affect all measured mixing ratios do not cause biases in estimated sources.