The seasonal cycle accounts for a dominant mode of total column CO2 (XCO2) annual variability and is connected to CO2 uptake and release; it thus represents an important quantity to test the accuracy ...of the measurements from space. We quantitatively evaluate the XCO2 seasonal cycle of the Greenhouse Gases Observing Satellite (GOSAT) observations from the Atmospheric CO2 Observations from Space (ACOS) retrieval system and compare average regional seasonal cycle features to those directly measured by the Total Carbon Column Observing Network (TCCON). We analyse the mean seasonal cycle amplitude, dates of maximum and minimum XCO2, as well as the regional growth rates in XCO2 through the fitted trend over several years. We find that GOSAT/ACOS captures the seasonal cycle amplitude within 1.0 ppm accuracy compared to TCCON, except in Europe, where the difference exceeds 1.0 ppm at two sites, and the amplitude captured by GOSAT/ACOS is generally shallower compared to TCCON. This bias over Europe is not as large for the other GOSAT retrieval algorithms (NIES v02.21, RemoTeC v2.35, UoL v5.1, and NIES PPDF-S v.02.11), although they have significant biases at other sites. We find that the ACOS bias correction partially explains the shallow amplitude over Europe. The impact of the co-location method and aerosol changes in the ACOS algorithm were also tested and found to be few tenths of a ppm and mostly non-systematic. We find generally good agreement in the date of minimum XCO2 between ACOS and TCCON, but ACOS generally infers a date of maximum XCO2 2–3 weeks later than TCCON. We further analyse the latitudinal dependence of the seasonal cycle amplitude throughout the Northern Hemisphere and compare the dependence to that predicted by current optimized models that assimilate in situ measurements of CO2. In the zonal averages, models are consistent with the GOSAT amplitude to within 1.4 ppm, depending on the model and latitude. We also show that the seasonal cycle of XCO2 depends on longitude especially at the mid-latitudes: the amplitude of GOSAT XCO2 doubles from western USA to East Asia at 45–50° N, which is only partially shown by the models. In general, we find that model-to-model differences can be larger than GOSAT-to-model differences. These results suggest that GOSAT/ACOS retrievals of the XCO2 seasonal cycle may be sufficiently accurate to evaluate land surface models in regions with significant discrepancies between the models.
The recent increase of atmospheric methane is investigated by using two atmospheric inversions to quantify the distribution of sources and sinks for the 2006-2008 period, and a process-based model of ...methane emissions by natural wetland ecosystems. Methane emissions derived from the two inversions are consistent at a global scale: emissions are decreased in 2006 (-7 Tg) and increased in 2007 (+21 Tg) and 2008 (+18 Tg), as compared to the 1999-2006 period. The agreement on the latitudinal partition of the flux anomalies for the two inversions is fair in 2006, good in 2007, and not good in 2008. In 2007, a positive anomaly of tropical emissions is found to be the main contributor to the global emission anomalies (~60-80%) for both inversions, with a dominant share attributed to natural wetlands (~2/3), and a significant contribution from high latitudes (~25%). The wetland ecosystem model produces smaller and more balanced positive emission anomalies between the tropics and the high latitudes for 2006, 2007 and 2008, mainly due to precipitation changes during these years. At a global scale, the agreement between the ecosystem model and the inversions is good in 2008 but not satisfying in 2006 and 2007. Tropical South America and Boreal Eurasia appear to be major contributors to variations in methane emissions consistently in the inversions and the ecosystem model. Finally, changes in OH radicals during 2006-2008 are found to be less than 1% in inversions, with only a small impact on the inferred methane emissions.
With the densification of surface observing networks and the development of remote sensing of greenhouse gases from space, estimations of methane (CH4) sources and sinks by inverse modeling are ...gaining additional constraining data but facing new challenges. The chemical transport model (CTM) linking the flux space to methane mixing ratio space must be able to represent these different types of atmospheric constraints for providing consistent flux estimations. Here we quantify the impact of sub-grid-scale physical parameterization errors on the global methane budget inferred by inverse modeling. We use the same inversion setup but different physical parameterizations within one CTM. Two different schemes for vertical diffusion, two others for deep convection, and one additional for thermals in the planetary boundary layer (PBL) are tested. Different atmospheric methane data sets are used as constraints (surface observations or satellite retrievals). At the global scale, methane emissions differ, on average, from 4.1 Tg CH4 per year due to the use of different sub-grid-scale parameterizations. Inversions using satellite total-column mixing ratios retrieved by GOSAT are less impacted, at the global scale, by errors in physical parameterizations. Focusing on large-scale atmospheric transport, we show that inversions using the deep convection scheme of Emanuel (1991) derive smaller interhemispheric gradients in methane emissions, indicating a slower interhemispheric exchange. At regional scale, the use of different sub-grid-scale parameterizations induces uncertainties ranging from 1.2 % (2.7 %) to 9.4 % (14.2 %) of methane emissions when using only surface measurements from a background (or an extended) surface network. Moreover, spatial distribution of methane emissions at regional scale can be very different, depending on both the physical parameterizations used for the modeling of the atmospheric transport and the observation data sets used to constrain the inverse system. When using only satellite data from GOSAT, we show that the small biases found in inversions using a coarser version of the transport model are actually masking a poor representation of the stratosphere-troposphere methane gradient in the model. Improving the stratosphere-troposphere gradient reveals a larger bias in GOSAT CH4 satellite data, which largely amplifies inconsistencies between the surface and satellite inversions. A simple bias correction is proposed. The results of this work provide the level of confidence one can have for recent methane inversions relative to physical parameterizations included in CTMs.
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° × 2.5° (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.
The land and ocean absorb on average just over half of the anthropogenic emissions of carbon dioxide (CO2) every year. These CO2 "sinks" are modulated by climate change and variability. Here we use a ...suite of nine dynamic global vegetation models (DGVMs) and four ocean biogeochemical general circulation models (OBGCMs) to estimate trends driven by global and regional climate and atmospheric CO2 in land and oceanic CO2 exchanges with the atmosphere over the period 1990–2009, to attribute these trends to underlying processes in the models, and to quantify the uncertainty and level of inter-model agreement. The models were forced with reconstructed climate fields and observed global atmospheric CO2; land use and land cover changes are not included for the DGVMs. Over the period 1990–2009, the DGVMs simulate a mean global land carbon sink of −2.4 ± 0.7 Pg C yr−1 with a small significant trend of −0.06 ± 0.03 Pg C yr−2 (increasing sink). Over the more limited period 1990–2004, the ocean models simulate a mean ocean sink of −2.2 ± 0.2 Pg C yr−1 with a trend in the net C uptake that is indistinguishable from zero (−0.01 ± 0.02 Pg C yr−2). The two ocean models that extended the simulations until 2009 suggest a slightly stronger, but still small, trend of −0.02 ± 0.01 Pg C yr−2. Trends from land and ocean models compare favourably to the land greenness trends from remote sensing, atmospheric inversion results, and the residual land sink required to close the global carbon budget. Trends in the land sink are driven by increasing net primary production (NPP), whose statistically significant trend of 0.22 ± 0.08 Pg C yr−2 exceeds a significant trend in heterotrophic respiration of 0.16 ± 0.05 Pg C yr−2 – primarily as a consequence of widespread CO2 fertilisation of plant production. Most of the land-based trend in simulated net carbon uptake originates from natural ecosystems in the tropics (−0.04 ± 0.01 Pg C yr−2), with almost no trend over the northern land region, where recent warming and reduced rainfall offsets the positive impact of elevated atmospheric CO2 and changes in growing season length on carbon storage. The small uptake trend in the ocean models emerges because climate variability and change, and in particular increasing sea surface temperatures, tend to counter\\-act the trend in ocean uptake driven by the increase in atmospheric CO2. Large uncertainty remains in the magnitude and sign of modelled carbon trends in several regions, as well as regarding the influence of land use and land cover changes on regional trends.
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 year's 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.
The Monitoring Atmospheric Composition and Climate Interim Implementation (MACC-II) delayed-mode (DM) system has been producing an atmospheric methane (CH4) analysis 6 months behind real time since ...June 2009. This analysis used to rely on the assimilation of the CH4 product from the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) instrument onboard Envisat. Recently the Laboratoire de Météorologie Dynamique (LMD) CH4 products from the Infrared Atmospheric Sounding Interferometer (IASI) and the SRON Netherlands Institute for Space Research CH4 products from the Thermal And Near-infrared Sensor for carbon Observation (TANSO) were added to the DM system. With the loss of Envisat in April 2012, the DM system now has to rely on the assimilation of methane data from TANSO and IASI. This paper documents the impact of this change in the observing system on the methane tropospheric analysis. It is based on four experiments: one free run and three analyses from respectively the assimilation of SCIAMACHY, TANSO and a combination of TANSO and IASI CH4 products in the MACC-II system. The period between December 2010 and April 2012 is studied. The SCIAMACHY experiment globally underestimates the tropospheric methane by 35 part per billion (ppb) compared to the HIAPER Pole-to-Pole Observations (HIPPO) data and by 28 ppb compared the Total Carbon Column Observing Network (TCCON) data, while the free run presents an underestimation of 5 ppb and 1 ppb against the same HIPPO and TCCON data, respectively. The assimilated TANSO product changed in October 2011 from version v.1 to version v.2.0. The analysis of version v.1 globally underestimates the tropospheric methane by 18 ppb compared to the HIPPO data and by 15 ppb compared to the TCCON data. In contrast, the analysis of version v.2.0 globally overestimates the column by 3 ppb. When the high density IASI data are added in the tropical region between 30° N and 30° S, their impact is mainly positive but more pronounced and effective when combined with version v.2.0 of the TANSO products. The resulting analysis globally underestimates the column-averaged dry-air mole fractions of methane (xCH4) just under 1 ppb on average compared to the TCCON data, whereas in the tropics it overestimates xCH4 by about 3 ppb. The random error is estimated to be less than 7 ppb when compared to TCCON data.
For the first time, carbon monoxide (CO) and formaldehyde (HCHO) satellite retrievals are used together with methane (CH4) and methyl choloroform (CH3CCl3 or MCF) surface measurements in an advanced ...inversion system. The CO and HCHO are respectively from the MOPITT and OMI instruments. The multi-species and multi-satellite dataset inversion is done for the 2005–2010 period. The robustness of our results is evaluated by comparing our posterior-modeled concentrations with several sets of independent measurements of atmospheric mixing ratios. The inversion leads to significant changes from the prior to the posterior, in terms of magnitude and seasonality of the CO and CH4 surface fluxes and of the HCHO production by non-methane volatile organic compounds (NMVOC). The latter is significantly decreased, indicating an overestimation of the biogenic NMVOC emissions, such as isoprene, in the GEIA inventory. CO and CH4 surface emissions are increased by the inversion, from 1037 to 1394 TgCO and from 489 to 529 TgCH4 on average for the 2005–2010 period. CH4 emissions present significant interannual variability and a joint CO-CH4 fluxes analysis reveals that tropical biomass burning probably played a role in the recent increase of atmospheric methane.
Over the last few years, solar‐induced chlorophyll fluorescence (SIF) observations from space have emerged as a promising resource for evaluating the spatio‐temporal distribution of gross primary ...productivity (GPP) simulated by global terrestrial biosphere models. SIF can be used to improve GPP simulations by optimizing critical model parameters through statistical Bayesian data assimilation techniques. A prerequisite is the availability of a functional link between GPP and SIF in terrestrial biosphere models. Here we present the development of a mechanistic SIF observation operator in the ORCHIDEE (Organizing Carbon and Hydrology In Dynamic Ecosystems) terrestrial biosphere model. It simulates the regulation of photosystem II fluorescence quantum yield at the leaf level thanks to a novel parameterization of non‐photochemical quenching as a function of temperature, photosynthetically active radiation, and normalized quantum yield of photochemistry. It emulates the radiative transfer of chlorophyll fluorescence to the top of the canopy using a parametric simplification of the SCOPE (Soil Canopy Observation Photosynthesis Energy) model. We assimilate two years of monthly OCO‐2 (Orbiting Carbon Observatory‐2) SIF product at 0.5° (2015–2016) to optimize ORCHIDEE photosynthesis and phenological parameters over an ensemble of grid points for all plant functional types. The impact on the simulated GPP is considerable with a large decrease of the global scale budget by 28 GtC/year over the period 1990–2009. The optimized GPP budget (134/136 GtC/year over 1990–2009/2001–2009) remarkably agrees with independent GPP estimates, FLUXSAT (137 GtC/year over 2001–2009) in particular and FLUXCOM (121 GtC/year over 1990–2009). Our results also suggest a biome dependency of the SIF‐GPP relationship that needs to be improved for some plant functional types.
Key Points
We developed a process‐based SIF observation operator in a terrestrial biosphere model
We assimilated monthly OCO‐2 SIF products to optimize model photosynthesis and phenology‐related parameters
The optimized GPP is considerably reduced with spatio‐temporal patterns in closer agreement with independent products
Satellite retrievals of methane weighted atmospheric columns are assimilated within a Bayesian inversion system to infer the global and regional methane emissions and sinks for the period August 2009 ...to July 2010. Inversions are independently computed from three different space-borne observing systems and one surface observing system under several hypotheses for prior-flux and observation errors. Posterior methane emissions are compared and evaluated against surface mole fraction observations via a chemistry-transport model. Apart from SCIAMACHY (SCanning Imaging Absorption spectroMeter for Atmospheric CartograpHY), the simulations agree fairly well with the surface mole fractions. The most consistent configurations of this study using TANSO-FTS (Thermal And Near infrared Sensor for carbon Observation – Fourier Transform Spectrometer), IASI (Infrared Atmospheric Sounding Interferometer) or surface measurements induce posterior methane global emissions of, respectively, 565 ± 21 Tg yr−1, 549 ± 36 Tg yr−1 and 538 ± 15 Tg yr−1 over the one-year period August 2009–July 2010. This consistency between the satellite retrievals (apart from SCIAMACHY) and independent surface measurements is promising for future improvement of CH4 emission estimates by atmospheric inversions.