Philippe Peylin

Recent patterns and mechanisms of carbon exchange by terrestrial ecosystems

Knowledge of carbon exchange between the atmosphere, land and the oceans is important, given that the terrestrial and marine environments are currently absorbing about half of the carbon dioxide that is emitted by fossil-fuel combustion. This carbon uptake is therefore limiting the extent of atmospheric and climatic change, but its long-term nature remains uncertain. Here we provide an overview of the current state of knowledge of global and regional patterns of carbon exchange by terrestrial ecosystems. Atmospheric carbon dioxide and oxygen data confirm that the terrestrial biosphere was largely neutral with respect to net carbon exchange during the 1980s, but became a net carbon sink in the 1990s. This recent sink can be largely attributed to northern extratropical areas, and is roughly split between North America and Eurasia. Tropical land areas, however, were approximately in balance with respect to carbon exchange, implying a carbon sink that offset emissions due to tropical deforestation. The evolution of the terrestrial carbon sink is largely the result of changes in land use over time, such as regrowth on abandoned agricultural land and fire prevention, in addition to responses to environmental changes, such as longer growing seasons, and fertilization by carbon dioxide and nitrogen. Nevertheless, there remain considerable uncertainties as to the magnitude of the sink in different regions and the contribution of different processes.

Net carbon dioxide losses of northern ecosystems in response to autumn warming

The carbon balance of terrestrial ecosystems is particularly sensitive to climatic changes in autumn and spring, with spring and autumn temperatures over northern latitudes having risen by about 1.1 °C and 0.8 °C, respectively, over the past two decades. A simultaneous greening trend has also been observed, characterized by a longer growing season and greater photosynthetic activity. These observations have led to speculation that spring and autumn warming could enhance carbon sequestration and extend the period of net carbon uptake in the future. Here we analyse interannual variations in atmospheric carbon dioxide concentration data and ecosystem carbon dioxide fluxes. We find that atmospheric records from the past 20 years show a trend towards an earlier autumn-to-winter carbon dioxide build-up, suggesting a shorter net carbon uptake period. This trend cannot be explained by changes in atmospheric transport alone and, together with the ecosystem flux data, suggest increasing carbon losses in autumn. We use a process-based terrestrial biosphere model and satellite vegetation greenness index observations to investigate further the observed seasonal response of northern ecosystems to autumnal warming. We find that both photosynthesis and respiration increase during autumn warming, but the increase in respiration is greater. In contrast, warming increases photosynthesis more than respiration in spring. Our simulations and observations indicate that northern terrestrial ecosystems may currently lose carbon dioxide in response to autumn warming, with a sensitivity of about 0.2 PgC °C-1, offsetting 90% of the increased carbon dioxide uptake during spring. If future autumn warming occurs at a faster rate than in spring, the ability of northern ecosystems to sequester carbon may be diminished earlier than previously suggested.

The carbon balance of terrestrial ecosystems in China

Global terrestrial ecosystems absorbed carbon at a rate of 1–4 Pg yr-1 during the 1980s and 1990s, offsetting 10–60 per cent of the fossil-fuel emissions. The regional patterns and causes of terrestrial carbon sources and sinks, however, remain uncertain. With increasing scientific and political interest in regional aspects of the global carbon cycle, there is a strong impetus to better understand the carbon balance of China. This is not only because China is the world’s most populous country and the largest emitter of fossil-fuel CO2 into the atmosphere, but also because it has experienced regionally distinct land-use histories and climate trends, which together control the carbon budget of its ecosystems. Here we analyse the current terrestrial carbon balance of China and its driving mechanisms during the 1980s and 1990s using three different methods: biomass and soil carbon inventories extrapolated by satellite greenness measurements, ecosystem models and atmospheric inversions. The three methods produce similar estimates of a net carbon sink in the range of 0.19–0.26 Pg carbon (PgC) per year, which is smaller than that in the conterminous United States but comparable to that in geographic Europe. We find that northeast China is a net source of CO2 to the atmosphere owing to overharvesting and degradation of forests. By contrast, southern China accounts for more than 65 per cent of the carbon sink, which can be attributed to regional climate change, large-scale plantation programmes active since the 1980s and shrub recovery. Shrub recovery is identified as the most uncertain factor contributing to the carbon sink. Our data and model results together indicate that China’s terrestrial ecosystems absorbed 28–37 per cent of its cumulated fossil carbon emissions during the 1980s and 1990s.

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.

Evaluation of terrestrial carbon cycle models for their response to climate variability and to CO2 trends

The purpose of this study was to evaluate 10 process-based terrestrial biosphere models that were used for the IPCC fifth Assessment Report. The simulated gross primary productivity (GPP) is compared with flux-tower-based estimates by Jung et al. [Journal of Geophysical Research 116 (2011) G00J07] (JU11). The net primary productivity (NPP) apparent sensitivity to climate variability and atmospheric CO2 trends is diagnosed from each model output, using statistical functions. The temperature sensitivity is compared against ecosystem field warming experiments results. The CO2 sensitivity of NPP is compared to the results from four Free-Air CO2 Enrichment (FACE) experiments. The simulated global net biome productivity (NBP) is compared with the residual land sink (RLS) of the global carbon budget from Friedlingstein et al. [Nature Geoscience 3 (2010) 811] (FR10). We found that models produce a higher GPP (133 ± 15 Pg C yr(-1) ) than JU11 (118 ± 6 Pg C yr(-1) ). In response to rising atmospheric CO2 concentration, modeled NPP increases on average by 16% (5-20%) per 100 ppm, a slightly larger apparent sensitivity of NPP to CO2 than that measured at the FACE experiment locations (13% per 100 ppm). Global NBP differs markedly among individual models, although the mean value of 2.0 ± 0.8 Pg C yr(-1) is remarkably close to the mean value of RLS (2.1 ± 1.2 Pg C yr(-1) ). The interannual variability in modeled NBP is significantly correlated with that of RLS for the period 1980-2009. Both model-to-model and interannual variation in model GPP is larger than that in model NBP due to the strong coupling causing a positive correlation between ecosystem respiration and GPP in the model. The average linear regression slope of global NBP vs. temperature across the 10 models is -3.0 ± 1.5 Pg C yr(-1) °C(-1) , within the uncertainty of what derived from RLS (-3.9 ± 1.1 Pg C yr(-1) °C(-1) ). However, 9 of 10 models overestimate the regression slope of NBP vs. precipitation, compared with the slope of the observed RLS vs. precipitation. With most models lacking processes that control GPP and NBP in addition to CO2 and climate, the agreement between modeled and observation-based GPP and NBP can be fortuitous. Carbon-nitrogen interactions (only separable in one model) significantly influence the simulated response of carbon cycle to temperature and atmospheric CO2 concentration, suggesting that nutrients limitations should be included in the next generation of terrestrial biosphere models.

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

Carbon residence time dominates uncertainty in terrestrial vegetation responses to future climate and atmospheric CO2

Future climate change and increasing atmospheric CO2 are expected to cause major changes in vegetation structure and function over large fractions of the global land surface. Seven global vegetation models are used to analyze possible responses to future climate simulated by a range of general circulation models run under all four representative concentration pathway scenarios of changing concentrations of greenhouse gases. All 110 simulations predict an increase in global vegetation carbon to 2100, but with substantial variation between vegetation models. For example, at 4 °C of global land surface warming (510–758 ppm of CO2), vegetation carbon increases by 52–477 Pg C (224 Pg C mean), mainly due to CO2 fertilization of photosynthesis. Simulations agree on large regional increases across much of the boreal forest, western Amazonia, central Africa, western China, and southeast Asia, with reductions across southwestern North America, central South America, southern Mediterranean areas, southwestern Africa, and southwestern Australia. Four vegetation models display discontinuities across 4 °C of warming, indicating global thresholds in the balance of positive and negative influences on productivity and biomass. In contrast to previous global vegetation model studies, we emphasize the importance of uncertainties in projected changes in carbon residence times. We find, when all seven models are considered for one representative concentration pathway × general circulation model combination, such uncertainties explain 30% more variation in modeled vegetation carbon change than responses of net primary productivity alone, increasing to 151% for non-HYBRID4 models. A change in research priorities away from production and toward structural dynamics and demographic processes is recommended.

CO2 surface fluxes at grid point scale estimated from a global 21 year reanalysis of atmospheric measurements

This paper documents a global Bayesian variational inversion of CO2 surface fluxes during the period 1988-2008. Weekly fluxes are estimated on a 3.75 degrees x 2.5 degrees (longitude-latitude) grid throughout the 21 years. The assimilated observations include 128 station records from three large data sets of surface CO2 mixing ratio measurements. A Monte Carlo approach rigorously quantifies the theoretical uncertainty of the inverted fluxes at various space and time scales, which is particularly important for proper interpretation of the inverted fluxes. Fluxes are evaluated indirectly against two independent CO2 vertical profile data sets constructed from aircraft measurements in the boundary layer and in the free troposphere. The skill of the inversion is evaluated by the improvement brought over a simple benchmark flux estimation based on the observed atmospheric growth rate. Our error analysis indicates that the carbon budget from the inversion should be more accurate than the a priori carbon budget by 20% to 60% for terrestrial fluxes aggregated at the scale of subcontinental regions in the Northern Hemisphere and over a year, but the inversion cannot clearly distinguish between the regional carbon budgets within a continent. On the basis of the independent observations, the inversion is seen to improve the fluxes compared to the benchmark: the atmospheric simulation of CO2 with the Bayesian inversion method is better by about 1 ppm than the benchmark in the free troposphere, despite possible systematic transport errors. The inversion achieves this improvement by changing the regional fluxes over land at the seasonal and at the interannual time scales. (Less)

Inferring CO2 sources and sinks from satellite observations: Method and application to TOVS data

Properly handling satellite data to constrain the inversion of CO2 sources and sinks at the Earth surface is a challenge motivated by the limitations of the current surface observation network. In this paper we present a Bayesian inference scheme to tackle this issue. It is based on the same theoretical principles as most inversions of the flask network but uses a variational formulation rather than a pure matrix-based one in order to cope with the large amount of satellite data. The minimization algorithm iteratively computes the optimum solution to the inference problem as well as an estimation of its error characteristics and some quantitative measures of the observation information content. A global climate model, guided by analyzed winds, provides information about the atmospheric transport to the inversion scheme. A surface flux climatology regularizes the inference problem. This new system has been applied to 1 years worth of retrievals of vertically integrated CO2 concentrations from the Television Infrared Observation Satellite Operational Vertical Sounder (TOVS). Consistent with a recent study that identified regional biases in the TOVS retrievals, the inferred fluxes are not useful for biogeochemical analyses. In addition to the detrimental impact of these biases, we find a sensitivity of the results to the formulation of the prior uncertainty and to the accuracy of the transport model. Notwithstanding these difficulties, four-dimensional inversion schemes of the type presented here could form the basis of multisensor data assimilation systems for the estimation of the surface fluxes of key atmospheric compounds.

Evidence for a weakening relationship between interannual temperature variability and northern vegetation activity.

Satellite-derived Normalized Difference Vegetation Index (NDVI), a proxy of vegetation productivity, is known to be correlated with temperature in northern ecosystems. This relationship, however, may change over time following alternations in other environmental factors. Here we show that above 30°N, the strength of the relationship between the interannual variability of growing season NDVI and temperature (partial correlation coefficient RNDVI-GT) declined substantially between 1982 and 2011. This decrease in RNDVI-GT is mainly observed in temperate and arctic ecosystems, and is also partly reproduced by process-based ecosystem model results. In the temperate ecosystem, the decrease in RNDVI-GT coincides with an increase in drought. In the arctic ecosystem, it may be related to a nonlinear response of photosynthesis to temperature, increase of hot extreme days and shrub expansion over grass-dominated tundra. Our results caution the use of results from interannual time scales to constrain the decadal response of plants to ongoing warming.