No abstract

[1] Monthly CO 2 fluxes are estimated across 1988-2003 for 22 emission regions using data from 78 CO 2 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 Niño/La Niña) 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 Niño and the post-Pinatubo uptake.

Abstract. Carbon storage by many terrestrial ecosystems can be limited by nutrients, predominantly nitrogen (N) and phosphorus (P), in addition to other environmental constraints, water, light and temperature. However the spatial distribution and the extent of both N and P limitation at the global scale have not been quantified. Here we have developed a global model of carbon (C), nitrogen (N) and phosphorus (P) cycles for the terrestrial biosphere. Model estimates of steady state C and N pool sizes and major fluxes between plant, litter and soil pools, under present climate conditions, agree well with various independent estimates. The total amount of C in the terrestrial biosphere is 2767 Gt C, and the C fractions in plant, litter and soil organic matter are 19%, 4% and 77%. The total amount of N is 135 Gt N, with about 94% stored in the soil, 5% in the plant live biomass, and 1% in litter. We found that the estimates of total soil P and its partitioning into different pools in soil are quite sensitive to biochemical P mineralization. The total amount of P (plant biomass, litter and soil) excluding occluded P in soil is 17 Gt P in the terrestrial biosphere, 33% of which is stored in the soil organic matter if biochemical P mineralization is modelled, or 31 Gt P with 67% in soil organic matter otherwise.This model was used to derive the global distribution and uncertainty of N or P limitation on the productivity of terrestrial ecosystems at steady state under present conditions. Our model estimates that the net primary productivity of most tropical evergreen broadleaf forests and tropical savannahs Correspondence to: Y. P. Wang (Yingping.wang@csiro.au) is reduced by about 20% on average by P limitation, and most of the remaining biomes are N limited; N limitation is strongest in high latitude deciduous needle leaf forests, and reduces its net primary productivity by up to 40% under present conditions.

[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 DpCO 2 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 DpCO 2 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.

Abstract. Atmospheric CO2 inversions estimate surface carbon fluxes from an optimal fit to atmospheric CO2 measurements, usually including prior constraints on the flux estimates. Eleven sets of carbon flux estimates are compared, generated by different inversions systems that vary in their inversions methods, choice of atmospheric data, transport model and prior information. The inversions were run for at least 5 yr in the period between 1990 and 2010. Mean fluxes for 2001–2004, seasonal cycles, interannual variability and trends are compared for the tropics and northern and southern extra-tropics, and separately for land and ocean. Some continental/basin-scale subdivisions are also considered where the atmospheric network is denser. Four-year mean fluxes are reasonably consistent across inversions at global/latitudinal scale, with a large total (land plus ocean) carbon uptake in the north (−3.4 Pg C yr−1 (±0.5 Pg C yr−1 standard deviation), with slightly more uptake over land than over ocean), a significant although more variable source over the tropics (1.6 ± 0.9 Pg C yr−1) and a compensatory sink of similar magnitude in the south (−1.4 ± 0.5 Pg C yr−1) corresponding mainly to an ocean sink. Largest differences across inversions occur in the balance between tropical land sources and southern land sinks. Interannual variability (IAV) in carbon fluxes is larger for land than ocean regions (standard deviation around 1.06 versus 0.33 Pg C yr−1 for the 1996–2007 period), with much higher consistency among the inversions for the land. While the tropical land explains most of the IAV (standard deviation ~ 0.65 Pg C yr−1), the northern and southern land also contribute (standard deviation ~ 0.39 Pg C yr−1). Most inversions tend to indicate an increase of the northern land carbon uptake from late 1990s to 2008 (around 0.1 Pg C yr−1, predominantly in North Asia. The mean seasonal cycle appears to be well constrained by the atmospheric data over the northern land (at the continental scale), but still highly dependent on the prior flux seasonality over the ocean. Finally we provide recommendations to interpret the regional fluxes, along with the uncertainty estimates.

Abstract. Results from the fully and biogeochemically coupled simulations in which CO2 increases at a rate of 1 % yr−1 (1pctCO2) from its preindustrial value are analyzed to quantify the magnitude of carbon–concentration and carbon–climate feedback parameters which measure the response of ocean and terrestrial carbon pools to changes in atmospheric CO2 concentration and the resulting change in global climate, respectively. The results are based on 11 comprehensive Earth system models from the most recent (sixth) Coupled Model Intercomparison Project (CMIP6) and compared with eight models from the fifth CMIP (CMIP5). The strength of the carbon–concentration feedback is of comparable magnitudes over land (mean ± standard deviation = 0.97 ± 0.40 PgC ppm−1) and ocean (0.79 ± 0.07 PgC ppm−1), while the carbon–climate feedback over land (−45.1 ± 50.6 PgC ∘C−1) is about 3 times larger than over ocean (−17.2 ± 5.0 PgC ∘C−1). The strength of both feedbacks is an order of magnitude more uncertain over land than over ocean as has been seen in existing studies. These values and their spread from 11 CMIP6 models have not changed significantly compared to CMIP5 models. The absolute values of feedback parameters are lower for land with models that include a representation of nitrogen cycle. The transient climate response to cumulative emissions (TCRE) from the 11 CMIP6 models considered here is 1.77 ± 0.37 ∘C EgC−1 and is similar to that found in CMIP5 models (1.63 ± 0.48 ∘C EgC−1) but with somewhat reduced model spread. The expressions for feedback parameters based on the fully and biogeochemically coupled configurations of the 1pctCO2 simulation are simplified when the small temperature change in the biogeochemically coupled simulation is ignored. Decomposition of the terms of these simplified expressions for the feedback parameters is used to gain insight into the reasons for differing responses among ocean and land carbon cycle models.

Abstract. The scientific understanding of the Earth's climate system, including the central question of how the climate system is likely to respond to human-induced perturbations, is comprehensively captured in GCMs and Earth System Models (ESM). Diagnosing the simulated climate response, and comparing responses across different models, is crucially dependent on transparent assumptions of how the GCM/ESM has been driven -especially because the implementation can involve subjective decisions and may differ between modelling groups performing the same experiment. This paper outlines the climate forcings and setup ofCorrespondence to: C. D. Jones (chris.d.jones@metoffice.gov.uk) the Met Office Hadley Centre ESM, HadGEM2-ES for the CMIP5 set of centennial experiments. We document the prescribed greenhouse gas concentrations, aerosol precursors, stratospheric and tropospheric ozone assumptions, as well as implementation of land-use change and natural forcings for the HadGEM2-ES historical and future experiments following the Representative Concentration Pathways. In addition, we provide details of how HadGEM2-ES ensemble members were initialised from the control run and how the palaeoclimate and AMIP experiments, as well as the "emissiondriven" RCP experiments were performed.

Abstract. A chemistry-transport model (CTM) intercomparison experiment (TransCom-CH4) has been designed to investigate the roles of surface emissions, transport and chemical loss in simulating the global methane distribution. Model simulations were conducted using twelve models and four model variants and results were archived for the period of 1990–2007. All but one model transports were driven by reanalysis products from 3 different meteorological agencies. The transport and removal of CH4 in six different emission scenarios were simulated, with net global emissions of 513 ± 9 and 514 ± 14 Tg CH4 yr−1 for the 1990s and 2000s, respectively. Additionally, sulfur hexafluoride (SF6) was simulated to check the interhemispheric transport, radon (222Rn) to check the subgrid scale transport, and methyl chloroform (CH3CCl3) to check the chemical removal by the tropospheric hydroxyl radical (OH). The results are compared to monthly or annual mean time series of CH4, SF6 and CH3CCl3 measurements from 8 selected background sites, and to satellite observations of CH4 in the upper troposphere and stratosphere. Most models adequately capture the vertical gradients in the stratosphere, the average long-term trends, seasonal cycles, interannual variations (IAVs) and interhemispheric (IH) gradients at the surface sites for SF6, CH3CCl3 and CH4. The vertical gradients of all tracers between the surface and the upper troposphere are consistent within the models, revealing vertical transport differences between models. An average IH exchange time of 1.39 ± 0.18 yr is derived from SF6 time series. Sensitivity simulations suggest that the estimated trends in exchange time, over the period of 1996–2007, are caused by a change of SF6 emissions towards the tropics. Using six sets of emission scenarios, we show that the decadal average CH4 growth rate likely reached equilibrium in the early 2000s due to the flattening of anthropogenic emission growth since the late 1990s. Up to 60% of the IAVs in the observed CH4 concentrations can be explained by accounting for the IAVs in emissions, from biomass burning and wetlands, as well as meteorology in the forward models. The modeled CH4 budget is shown to depend strongly on the troposphere-stratosphere exchange rate and thus on the model's vertical grid structure and circulation in the lower stratosphere. The 15-model median CH4 and CH3CCl3 atmospheric lifetimes are estimated to be 9.99 ± 0.08 and 4.61 ± 0.13 yr, respectively, with little IAV due to transport and temperature.