Greenhouse Gases Observing SATellite XCO2 Research Articles

Improved Consistency of Satellite XCO2 Retrievals Based on Machine Learning

AbstractQuantifying atmospheric CO2 over long periods from space is crucial in understanding the carbon cycle's response to climate change. However, a single satellite offers limited spatiotemporal coverage, making comprehensive monitoring challenging. Moreover, biases among various satellite retrievals hinder their direct integration. This study proposed a machine learning framework for fusing the column‐averaged dry‐air mole fraction of CO2 (XCO2) retrievals from Greenhouse Gases Observing Satellite (GOSAT) and OCO‐2 satellites. The best model (R2 = 0.85) presented improved consistency of GOSAT retrievals by reducing 71.5% of the average monthly bias while using OCO‐2 retrievals as a benchmark, indicating the fusion data set's potential to enhance observation coverage. Incorporating the adjusted GOSAT XCO2 retrievals into the OCO‐2 data set added an average of 84.7 thousand observations annually, enhancing the yearly temporal coverage by 53.6% (from 14 to 21.5 days per grid). This method can be adapted to other satellites, maximizing satellite resources for a more robust carbon flux inversion.

Influences of aerosols and thin cirrus clouds on GOSAT XCO2 and XCH4 using Total Carbon Column Observing Network, sky radiometer, and lidar data

ABSTRACT The column-averaged, dry air mole fractions of CO2 and CH4 (XCO2 and XCH4, respectively) were retrieved from short-wavelength infrared (SWIR) spectra observed by the Greenhouse gases Observing SATellite (GOSAT). Continuous measurements of SWIR spectra have been made via GOSAT since 2009, but there has been insufficient investigation of the effects of cirrus clouds and aerosols on the observations. In this work, we investigated the influences of aerosols and cirrus clouds on the differences between GOSAT observations and Total Carbon Column Observing Network (TCCON) data for XCO2 and XCH4 (ΔXCO2 and ΔXCH4) at three sites: Tsukuba and Saga in Japan, and Lauder in New Zealand. We used aerosol optical thickness (AOT), Angstrom exponents (AEs), and single scattering albedo (SSA), all obtained from sky radiometer observations, as well as vertical profiles of aerosols and thin cirrus clouds from lidar observations. Matchups were performed within ±0.1° latitude/longitude rectangular areas of each TCCON site, and within 30 min of the GOSAT overpass time. The results show a negative slope between ΔXCO2 and AOT at 500 nm determined from sky radiometer data at Tsukuba and Saga. The GOSAT XCO2 values tended to be lowered in the presence of cirrus clouds and dense boundary-layer aerosols. Moreover, a significant negative ΔXCO2 was observed at times of large AOTs that resulted from dust-like events. At Lauder, ΔXCO2 was negatively correlated with the AOT at 500 nm, although the AOT at this site was generally small. The mean ± standard deviation for ΔXCO2 and ΔXCH4 at Lauder are −0.80 ± 1.83 (ppm) and −5.27 ± 10.79 (ppb) with correlation coefficients r between GOSAT and TCCON of 0.94 and 0.83, respectively. Both ΔXCO2 and ΔXCH4 were significantly and negatively correlated with the AOT during Sep-Oct-Nov. In addition, stratospheric aerosols caused large negative biases of ΔXCO2 and ΔXCH4 at Lauder despite the small stratospheric aerosol optical depth at that site.

An 11-year record of XCO<sub>2</sub> estimates derived from GOSAT measurements using the NASA ACOS version 9 retrieval algorithm

Abstract. The Thermal And Near infrared Sensor for carbon Observation – Fourier Transform Spectrometer (TANSO-FTS) on the Japanese Greenhouse gases Observing SATellite (GOSAT) has been returning data since April 2009. The version 9 (v9) Atmospheric Carbon Observations from Space (ACOS) Level 2 Full Physics (L2FP) retrieval algorithm (Kiel et al., 2019) was used to derive estimates of carbon dioxide (CO2) dry air mole fraction (XCO2) from the TANSO-FTS measurements collected over its first 11 years of operation. The bias correction and quality filtering of the L2FP XCO2 product were evaluated using estimates derived from the Total Carbon Column Observing Network (TCCON) as well as values simulated from a suite of global atmospheric inversion systems (models) which do not assimilate satellite-derived CO2. In addition, the v9 ACOS GOSAT XCO2 results were compared with collocated XCO2 estimates derived from NASA's Orbiting Carbon Observatory-2 (OCO-2), using the version 10 (v10) ACOS L2FP algorithm. These tests indicate that the v9 ACOS GOSAT XCO2 product has improved throughput, scatter, and bias, when compared to the earlier v7.3 ACOS GOSAT product, which extended through mid 2016. Of the 37 million soundings collected by GOSAT through June 2020, approximately 20 % were selected for processing by the v9 L2FP algorithm after screening for clouds and other artifacts. After post-processing, 5.4 % of the soundings (2×106 out of 37×106) were assigned a “good” XCO2 quality flag, as compared to 3.9 % in v7.3 (<1 ×106 out of 24×106). After quality filtering and bias correction, the differences in XCO2 between ACOS GOSAT v9 and both TCCON and models have a scatter (1σ) of approximately 1 ppm for ocean-glint observations and 1 to 1.5 ppm for land observations. Global mean biases against TCCON and models are less than approximately 0.2 ppm. Seasonal mean biases relative to the v10 OCO-2 XCO2 product are of the order of 0.1 ppm for observations over land. However, for ocean-glint observations, seasonal mean biases relative to OCO-2 range from 0.2 to 0.6 ppm, with substantial variation in time and latitude. The ACOS GOSAT v9 XCO2 data are available on the NASA Goddard Earth Science Data and Information Services Center (GES-DISC) in both the per-orbit full format (https://doi.org/10.5067/OSGTIL9OV0PN, OCO-2 Science Team et al., 2019b) and in the per-day lite format (https://doi.org/10.5067/VWSABTO7ZII4, OCO-2 Science Team et al., 2019a). In addition, a new set of monthly super-lite files, containing only the most essential variables for each satellite observation, has been generated to provide entry level users with a light-weight satellite product for initial exploration (CaltechDATA, https://doi.org/10.22002/D1.2178, Eldering, 2021). The v9 ACOS Data User's Guide (DUG) describes best-use practices for the GOSAT data (O'Dell et al., 2020). The GOSAT v9 data set should be especially useful for studies of carbon cycle phenomena that span a full decade or more and may serve as a useful complement to the shorter OCO-2 v10 data set, which begins in September 2014.

Comparative Evaluation of Top-Down GOSAT XCO2 vs. Bottom-Up National Reports in the European Countries

Submitting national inventory reports (NIRs) on emissions of greenhouse gases (GHGs) is obligatory for parties of the United Nations Framework Convention on Climate Change (UNFCCC). The NIR forms the basis for monitoring individual countries’ progress on mitigating climate change. Countries prepare NIRs using the default bottom–up methodology of the Intergovernmental Panel on Climate Change (IPCC), as approved by the Kyoto protocol. We provide tangible evidence of the discrepancy between official bottom–up NIR reporting (unit: tons) versus top–down XCO2 reporting (unit: ppm) within the European continent, as measured by the Greenhouse Gases Observing Satellite (GOSAT). Bottom–up NIR (annual growth rate of CO2 emission from 2010 to 2016: −1.55%) does not show meaningful correlation (geographically weighted regression coefficient = −0.001, R2 = 0.024) to top–down GOSAT XCO2 (annual growth rate: 0.59%) in the European countries. The top five countries within the European continent on carbon emissions in NIR do not match the top five countries on GOSAT XCO2 concentrations. NIR exhibits anthropogenic carbon-generating activity within country boundaries, whereas satellite signals reveal the trans-boundary movement of natural and anthropogenic carbon. Although bottom–up NIR reporting has already gained worldwide recognition as a method to track national follow-up for treaty obligations, the single approach based on bottom–up did not present background atmospheric CO2 density derived from the air mass movement between the countries. In conclusion, we suggest an integrated measuring, reporting, and verification (MRV) approach using top–down observation in combination with bottom–up NIR that can provide sufficient countrywide objective evidence for national follow-up activities.

Evaluation of earth system model and atmospheric inversion using total column CO2 observations from GOSAT and OCO-2

The measurements of one of the major greenhouse gases, carbon dioxide (CO2), are being made using dedicated satellite remote sensing since the launch of the greenhouse gases observing satellite (GOSAT) by a three-way partnership between the Japan Aerospace Exploration Agency (JAXA), the Ministry of Environment (MoE) and the National Institute for Environmental Studies (NIES), and the National Aeronautics and Space Administration (NASA) Orbiting Carbon Observatory-2 (OCO-2). In the past 10 years, estimation of CO2 fluxes from land and ocean using the earth system models (ESMs) and inverse modelling of in situ atmospheric CO2 data have also made significant progress. We attempt, for the first time, to evaluate the CO2 fluxes simulated by an earth system model (MIROC-ES2L) and the fluxes estimated by an inverse model (MIROC4-Inv) using in situ data by comparing with GOSAT and OCO-2 observations. Both MIROC-ES2L and MIROC4-Inv fluxes are used in the MIROC4-atmospheric chemistry transport model (referred to as ACTM_ES2LF and ACTM_InvF, respectively) for calculating total column CO2 mole fraction (XCO2) that are sampled at the time and location of the satellite measurements. Both the ACTM simulations agreed well with the GOSAT and OCO-2 satellite observations, within 2 ppm for the spatial maps and time evolutions of the zonal mean distributions. Our results suggest that the inverse model using in situ data is more consistent with the OCO-2 retrievals, compared with those of the GOSAT XCO2 data due to the higher accuracy of the former. This suggests that the MIROC4-Inv fluxes are of sufficient quality to evaluate MIROC-ES2L simulated fluxes. The ACTM_ES2LF simulation shows a slightly weaker seasonal cycle for the meridional profiles of CO2 fluxes, compared with that from the ACTM_InvF. This difference is revealed by greater XCO2 differences for ACTM_ES2LF vs GOSAT, compared with those of ACTM_InvF vs GOSAT. Using remote sensing–based global products of leaf area index (LAI) and gross primary productivity (GPP) over land, we show a weaker sensitivity of MIROC-ES2L biospheric activities to the weather and climate in the tropical regions. Our results clearly suggest the usefulness of XCO2 measurements by satellite remote sensing for evaluation of large-scale ESMs, which so far remained untested by the sparse in situ data.

Specific patterns of XCO2 observed by GOSAT during 2009–2016 and assessed with model simulations over China

Spatiotemporal patterns of column-averaged dry air mole fraction of CO2 (XCO2) have not been well characterized on a regional scale due to limitations in data availability and precision. This paper addresses these issues by examining such patterns in China using the long-term mapping XCO2 dataset (2009–2016) derived from the Greenhouse gases Observing SATellite (GOSAT). XCO2 simulations are also constructed using the high-resolution nested-grid GEOS-Chem model. The following results are found: Firstly, the correlation coefficient between the anthropogenic emissions and XCO2 spatial distribution is nearly zero in summer but up to 0.32 in autumn. Secondly, on average, XCO2 increases by 2.08 ppm every year from 2010 to 2015, with a sharp increase of 2.6 ppm in 2013. Lastly, in the analysis of three typical regions, the GOSAT XCO2 time series is in better agreement with the GEOS-Chem simulation of XCO2 in the Taklimakan Desert region (the least difference with bias 0.65±0.78 ppm), compared with the northern urban agglomeration region (−1.3±1.2 ppm) and the northeastern forest region (−1.4±1.4 ppm). The results are likely attributable to uncertainty in both the satellite-retrieved XCO2 data and the model simulation data.

Identifying local anthropogenic CO2 emissions with satellite retrievals: a case study in South Korea

ABSTRACTWe used multiyear Greenhouse Gases Observing Satellite (GOSAT) dry air, column-integrated CO2 (XCO2) retrievals (2010–2013) to evaluate urban and local-scale CO2 emissions over East Asia and examined whether GOSAT XCO2 captures the impact of strong local CO2 emissions over South Korea, an East Asian downwind region with high atmospheric aerosol loading and strong summer monsoons. We chose a region in western Mongolia (upwind region) as the XCO2 background, and estimated XCO2 enhancements in South Korea to gauge local and regional emissions. We found that the cold season (November–February) was better suited for estimating XCO2 enhancements of local emissions due to the summer monsoon and stronger transboundary impacts in other seasons. In particular, we focused on three local GOSAT XCO2 footprints (about 10.5 km in diameter) in South Korea: the Seoul Metropolitan Area (SMA), the Gwangyang Steelworks and Hadong power plants (GYG), and the Samcheonpo power plants (SCH). The range of XCO2 enhancement was 7.3–10.7 ppm (14.1–21.3 mg m−3 in standard temperature and pressure (STP)). By estimating other important contributions to XCO2 enhancements such as the XCO2 latitudinal gradients and Chinese fossil fuel combustions, we estimated the net enhancements caused mainly by local CO2 emissions in the range of 4.2–7.6 ppm (8.1–14.7 mg m−3 in STP) These high enhancements imply that large point source contributions are an important factor in determining these enhancements, even if contributions are also made by broader-scale emissions. Additionally, differences in net XCO2 enhancements and trends between GYG (+ 4.2 ppm (+ 8.2 mg m−3 in STP), – 0.2 ppm year−1 (–0.4 mg m−3 year−1 in STP)) and SCH (+ 7.6 ppm (+ 14.9 mg m−3 in STP), + 1.3 ppm year−1 (+ 2.6 mg· m−3 year−1 in STP)) indicate that these closely located footprints (approximately 26 km apart) are separable. These results will be useful in evaluating and reducing uncertainties in regional and local anthropogenic greenhouse gas (GHG) emissions over East Asia.

Regional uncertainty of GOSAT XCO<sub>2</sub> retrievals in China: quantification and attribution

Abstract. The regional uncertainty of the column-averaged dry air mole fraction of CO2 (XCO2) retrieved using different algorithms from the Greenhouse gases Observing SATellite (GOSAT) and its attribution are still not well understood. This paper investigates the regional performance of XCO2 within a latitude band of 37–42° N segmented into 8 cells in a grid of 5° from west to east (80–120° E) in China, where typical land surface types and geographic conditions exist. The former includes desert, grassland and built-up areas mixed with cropland; and the latter includes anthropogenic emissions that change from small to large from west to east, including those from the megacity of Beijing. For these specific cells, we evaluate the regional uncertainty of GOSAT XCO2 retrievals by quantifying and attributing the consistency of XCO2 retrievals from four algorithms (ACOS, NIES, OCFP and SRFP) by intercomparison. These retrievals are then specifically compared with simulated XCO2 from the high-resolution nested model in East Asia of the Goddard Earth Observing System 3-D chemical transport model (GEOS-Chem). We also introduce the anthropogenic CO2 emissions data generated from the investigation of surface emitting point sources that was conducted by the Ministry of Environmental Protection of China to GEOS-Chem simulations of XCO2 over the Chinese mainland. The results indicate that (1) regionally, the four algorithms demonstrate smaller absolute biases of 0.7–1.1 ppm in eastern cells, which are covered by built-up areas mixed with cropland with intensive anthropogenic emissions, than those in the western desert cells (1.0–1.6 ppm) with a high-brightness surface from the pairwise comparison results of XCO2 retrievals. (2) Compared with XCO2 simulated by GEOS-Chem (GEOS-XCO2), the XCO2 values from ACOS and SRFP have better agreement, while values from OCFP are the least consistent with GEOS-XCO2. (3) Viewing attributions of XCO2 in the spatio-temporal pattern, ACOS and SRFP demonstrate similar patterns, while OCFP is largely different from the others. In conclusion, the discrepancy in the four algorithms is the smallest in eastern cells in the study area, where the megacity of Beijing is located and where there are strong anthropogenic CO2 emissions, which implies that XCO2 from satellite observations could be reliably applied in the assessment of atmospheric CO2 enhancements induced by anthropogenic CO2 emissions. The large inconsistency among the four algorithms presented in western deserts which displays a high albedo and dust aerosols, moreover, demonstrates that further improvement is still necessary in such regions, even though many algorithms have endeavored to minimize the effects of aerosols scattering and surface albedo.

Lower-tropospheric CO<sub>2</sub> from near-infrared ACOS-GOSAT observations

Abstract. We present two new products from near-infrared Greenhouse Gases Observing Satellite (GOSAT) observations: lowermost tropospheric (LMT, from 0 to 2.5 km) and upper tropospheric–stratospheric (U, above 2.5 km) carbon dioxide partial column mixing ratios. We compare these new products to aircraft profiles and remote surface flask measurements and find that the seasonal and year-to-year variations in the new partial column mixing ratios significantly improve upon the Atmospheric CO2 Observations from Space (ACOS) and GOSAT (ACOS-GOSAT) initial guess and/or a priori, with distinct patterns in the LMT and U seasonal cycles that match validation data. For land monthly averages, we find errors of 1.9, 0.7, and 0.8 ppm for retrieved GOSAT LMT, U, and XCO2; for ocean monthly averages, we find errors of 0.7, 0.5, and 0.5 ppm for retrieved GOSAT LMT, U, and XCO2. In the southern hemispheric biomass burning season, the new partial columns show similar patterns to MODIS fire maps and MOPITT multispectral CO for both vertical levels, despite a flat ACOS-GOSAT prior, and a CO–CO2 emission factor comparable to published values. The difference of LMT and U, useful for evaluation of model transport error, has also been validated with a monthly average error of 0.8 (1.4) ppm for ocean (land). LMT is more locally influenced than U, meaning that local fluxes can now be better separated from CO2 transported from far away.

Study of the footprints of short-term variation in XCO<sub>2</sub> observed by TCCON sites using NIES and FLEXPART atmospheric transport models

Abstract. The Total Carbon Column Observing Network (TCCON) is a network of ground-based Fourier transform spectrometers (FTSs) that record near-infrared (NIR) spectra of the sun. From these spectra, accurate and precise observations of CO2 column-averaged dry-air mole fractions (denoted XCO2) are retrieved. TCCON FTS observations have previously been used to validate satellite estimations of XCO2; however, our knowledge of the short-term spatial and temporal variations in XCO2 surrounding the TCCON sites is limited. In this work, we use the National Institute for Environmental Studies (NIES) Eulerian three-dimensional transport model and the FLEXPART (FLEXible PARTicle dispersion model) Lagrangian particle dispersion model (LPDM) to determine the footprints of short-term variations in XCO2 observed by operational, past, future and possible TCCON sites. We propose a footprint-based method for the collocation of satellite and TCCON XCO2 observations and estimate the performance of the method using the NIES model and five GOSAT (Greenhouse Gases Observing Satellite) XCO2 product data sets. Comparison of the proposed approach with a standard geographic method shows a higher number of collocation points and an average bias reduction up to 0.15 ppm for a subset of 16 stations for the period from January 2010 to January 2014. Case studies of the Darwin and Reunion Island sites reveal that when the footprint area is rather curved, non-uniform and significantly different from a geographical rectangular area, the differences between these approaches are more noticeable. This emphasises that the collocation is sensitive to local meteorological conditions and flux distributions.

Interpreting Temporal Changes of Atmospheric CO2Over Fire Affected Regions Based on GOSAT Observations

The carbon dioxide (CO2) emissions released from biomass burning significantly affect the temporal variations of atmospheric CO2 concentrations. Based on a long-term (July 2009–June 2015) retrieved data sets by the greenhouse gases observing satellite (GOSAT), the seasonal cycle and interannual variations of column-averaged volume mixing ratios of atmospheric carbon dioxide (XCO2) in four fire affected continental regions were analyzed. The results showed that Northern Africa (NA) had the largest seasonal variations after removing its regional trend of XCO2 with peak-to-peak amplitude of 6.2 ppm within the year, higher than central South America (CSA) (2.4 ppm), Southern Africa (SA) (3.8 ppm), and Australia (1.7 ppm). The detrended regional XCO2 ( $\triangle $ XCO2) was found to be positively correlated with the fire CO2 emissions during the fire activity period but with different seasonal variabilities. NA recorded the largest change of seasonal variations of $\triangle $ XCO2 with a total of 12.8 ppm during fire seasons, higher than CSA, SA, and Australia with 5.4, 6.7, and 2.2 ppm, respectively. During the fire episode, the positive $\triangle $ XCO2 was noticed during June–November in CSA, December to next June in NA, and May–November in SA. The Pearson correlation coefficients between the variations of $\triangle $ XCO2 and fire CO2 emissions achieved the best correlations in SA ( $R = 0.77$ and $p ). This letter revealed that fire CO2 emissions and GOSAT XCO2 presented consistent seasonal variations.

Inter-annual variability of summertime CO2 exchange in Northern Eurasia inferred from GOSAT XCO2

Northern Eurasia is one of the largest terrestrial carbon reservoirs on the Earth’s surface. However, since the coverage of surface CO2 observations is still limited, the response to the climate variability remains uncertain. We estimated monthly CO2 fluxes for three sub-regions in Northern Eurasia (north of ∼60°N), Northeastern Europe, Western Siberia and Eastern Siberia, using CO2 retrievals from the Japanese Greenhouse Gases Observing SATellite (GOSAT). The variations of estimated CO2 fluxes were examined in terms of the regional climate variability, for the three consecutive growing seasons of 2009–2011. The CO2 fluxes estimated using GOSAT data are highly correlated with the surface temperature anomalies in July and August (r > 0.8) while no correlation is found in the CO2 fluxes estimated only using surface observations. The estimated fluxes from GOSAT data exhibit high negative correlations with one-month lagged positive precipitation anomalies in late summer (r > −0.7) through surface temperature and the Normalized Difference Vegetation Index (NDVI). The results indicate that GOSAT data reflects the changes in terrestrial biospheric processes responding to climate anomalies. In 2010, a large part of Eurasia experienced an extremely hot and dry summer, while cold and wet weather conditions were recorded in Western Siberia. The CO2 fluxes estimated from GOSAT data showed a reduction of net CO2 uptake in Northeastern Europe and Eastern Siberia, but the enhancement of net CO2 uptake in Western Siberia. These opposite sub-regional flux anomalies can be explained by the different climate anomalies on a sub-regional scale in Northern Eurasia. Thus, this study demonstrates that space-based observations by GOSAT compensate for the lack of ground-based observational coverage so as to better capture the inter-annually varying atmosphere-terrestrial biosphere CO2 exchange on a regional scale.

Ability of the 4-D-Var analysis of the GOSAT BESD XCO<sub>2</sub> retrievals to characterize atmospheric CO<sub>2</sub> at large and synoptic scales

Abstract. This study presents results from the European Centre for Medium-Range Weather Forecasts (ECMWF) carbon dioxide (CO2) analysis system where the atmospheric CO2 is controlled through the assimilation of column-averaged dry-air mole fractions of CO2 (XCO2) from the Greenhouse gases Observing Satellite (GOSAT). The analysis is compared to a free-run simulation (without assimilation of XCO2), and they are both evaluated against XCO2 data from the Total Carbon Column Observing Network (TCCON). We show that the assimilation of the GOSAT XCO2 product from the Bremen Optimal Estimation Differential Optical Absorption Spectroscopy (BESD) algorithm during the year 2013 provides XCO2 fields with an improved mean absolute error of 0.6 parts per million (ppm) and an improved station-to-station bias deviation of 0.7 ppm compared to the free run (1.1 and 1.4 ppm, respectively) and an improved estimated precision of 1 ppm compared to the GOSAT BESD data (3.3 ppm). We also show that the analysis has skill for synoptic situations in the vicinity of frontal systems, where the GOSAT retrievals are sparse due to cloud contamination. We finally computed the 10-day forecast from each analysis at 00:00 UTC, and we demonstrate that the CO2 forecast shows synoptic skill for the largest-scale weather patterns (of the order of 1000 km) even up to day 5 compared to its own analysis.

Estimates of European uptake of CO<sub>2</sub> inferred from GOSAT X<sub>CO<sub>2</sub></sub> retrievals: sensitivity to measurement bias inside and outside Europe

Abstract. Estimates of the natural CO2 flux over Europe inferred from in situ measurements of atmospheric CO2 mole fraction have been used previously to check top-down flux estimates inferred from space-borne dry-air CO2 column (XCO2) retrievals. Several recent studies have shown that CO2 fluxes inferred from XCO2 data from the Japanese Greenhouse gases Observing SATellite (GOSAT) and the Scanning Imaging Absorption Spectrometer for Atmospheric CHartographY (SCIAMACHY) have larger seasonal amplitudes and a more negative annual net CO2 balance than those inferred from the in situ data. The cause of this elevated European uptake of CO2 is still unclear, but some recent studies have suggested that this is a genuine scientific phenomenon. Here, we put forward an alternative hypothesis and show that realistic levels of bias in GOSAT data can result in an erroneous estimate of elevated uptake over Europe. We use a global flux inversion system to examine the relationship between measurement biases and estimates of CO2 uptake from Europe. We establish a reference in situ inversion that uses an Ensemble Kalman Filter (EnKF) to assimilate conventional surface mole fraction observations and XCO2 retrievals from the surface-based Total Carbon Column Observing Network (TCCON). We use the same EnKF system to assimilate two independent versions of GOSAT XCO2 data. We find that the GOSAT-inferred European terrestrial biosphere uptake peaks during the summer, similar to the reference inversion, but the net annual flux is 1.40 ± 0.19 GtC a−1 compared to a value of 0.58 ± 0.14 GtC a−1 for our control inversion that uses only in situ data. To reconcile these two estimates, we perform a series of numerical experiments that assimilate observations with added biases or assimilate synthetic observations for which part or all of the GOSAT XCO2 data are replaced with model data. We find that for our global flux inversions, a large portion (60–90 %) of the elevated European uptake inferred from GOSAT data in 2010 is due to retrievals outside the immediate European region, while the remainder can largely be explained by a sub-ppm retrieval bias over Europe. We use a data assimilation approach to estimate monthly GOSAT XCO2 biases from the joint assimilation of in situ observations and GOSAT XCO2 retrievals. The inferred biases represent an estimate of systematic differences between GOSAT XCO2 retrievals and the inversion system at regional or sub-regional scales. We find that a monthly varying bias of up to 0.5 ppm can explain an overestimate of the annual sink of up to 0.20 GtC a−1. Our results highlight the sensitivity of CO2 flux estimates to regional observation biases, which have not been fully characterized by the current observation network. Without further dedicated measurements we cannot prove or disprove that European ecosystems are taking up a larger-than-expected amount of CO2. More robust inversion systems are also needed to infer consistent fluxes from multiple observation types.

CO2 annual and semiannual cycles from multiple satellite retrievals and models

AbstractSatellite CO2 retrievals from the Greenhouse gases Observing SATellite (GOSAT), Atmospheric Infrared Sounder (AIRS), and Tropospheric Emission Spectrometer (TES) and in situ measurements from the National Oceanic and Atmospheric Administration ‐ Earth System Research Laboratory (NOAA‐ESRL) Surface CO2 and Total Carbon Column Observing Network (TCCON) are utilized to explore the CO2 variability at different altitudes. A multiple regression method is used to calculate the CO2 annual cycle and semiannual cycle amplitudes from different data sets. The CO2 annual cycle and semiannual cycle amplitudes for GOSAT XCO2 and TCCON XCO2 are consistent but smaller than those seen in the NOAA‐ESRL surface data. The CO2 annual and semiannual cycles are smallest in the AIRS midtropospheric CO2 compared with other data sets in the Northern Hemisphere. The amplitudes for the CO2 annual cycle and semiannual cycle from GOSAT, TES, and AIRS CO2 are small and comparable to each other in the Southern Hemisphere. Similar regression analysis is applied to the Model for OZone And Related chemical Tracers‐2 and CarbonTracker model CO2. The convolved model CO2 annual cycle and semiannual cycle amplitudes are similar to those from the satellite CO2 retrievals, although the models tend to underestimate the CO2 seasonal cycle amplitudes in the Northern Hemisphere midlatitudes and underestimate the CO2 semiannual cycle amplitudes in the high latitudes. These results can be used to better understand the vertical structures for the CO2 annual cycle and semiannual cycle and help identify deficiencies in the models, which are very important for the carbon budget study.

Does GOSAT capture the true seasonal cycle of carbon dioxide?

Abstract. 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.

Source‐receptor relationships of column‐average CO2 and implications for the impact of observations on flux inversions

AbstractSource‐receptor relationships are the fundamental quantities used in atmosphere CO2 flux inversions. In this study, we systematically investigate the global source‐receptor relationships of column CO2 (XCO2) in 12 continental‐scale receptors in terms of transport and local versus nonlocal flux contributions using the GEOS‐Chem adjoint model. Using simulated Greenhouse gases Observing Satellite (GOSAT) XCO2, we quantify the impact of inclusion (add‐on) or exclusion of observations (data‐denial) within a receptor region on flux inversion. We discuss the connections between the observation impact and the underlying source‐receptor relationships. The strong sensitivity of XCO2 to nonlocal fluxes makes the XCO2 observations have strong impact on nonlocal flux estimation. On an annual mean, the impact of GOSAT XCO2 over Europe on North America flux estimation is 10%–39% of the full observation impact. Because the Southern Hemisphere midlatitude XCO2 are most sensitive to local and tropical fluxes, the mean impact of GOSAT XCO2 over Southern South America on Northern South America (N‐S‐America) flux estimation is 30%–59% of the full observation impact. Because of the strong sensitivity of the XCO2 over N‐S‐America to Central Africa (C‐Africa) fluxes, the mean impact on C‐Africa flux estimation is between 14% and 31% of the full observation impact in spite of the sparse observation coverage. The results also show that XCO2 have similar sensitivity to local and nonlocal fluxes occurring 3 months before, which indicates that XCO2 observations cannot differentiate any fluxes occurring beyond 3 months.

Inferring regional sources and sinks of atmospheric CO<sub>2</sub> from GOSAT XCO<sub>2</sub> data

Abstract. We have examined the utility of retrieved column-averaged, dry-air mole fractions of CO2 (XCO2) from the Greenhouse Gases Observing Satellite (GOSAT) for quantifying monthly, regional flux estimates of CO2, using the GEOS-Chem four-dimensional variational (4D-Var) data assimilation system. We focused on assessing the potential impact of biases in the GOSAT CO2 data on the regional flux estimates. Using different screening and bias correction approaches, we selected three different subsets of the GOSAT XCO2 data for the 4D-Var inversion analyses, and found that the inferred global fluxes were consistent across the three XCO2 inversions. However, the GOSAT observational coverage was a challenge for the regional flux estimates. In the northern extratropics, the inversions were more sensitive to North American fluxes than to European and Asian fluxes due to the lack of observations over Eurasia in winter and over eastern and southern Asia in summer. The regional flux estimates were also sensitive to the treatment of the residual bias in the GOSAT XCO2 data. The largest differences obtained were for temperate North America and temperate South America, for which the largest spread between the inversions was 1.02 and 0.96 Pg C, respectively. In the case of temperate North America, one inversion suggested a strong source, whereas the second and third XCO2 inversions produced a weak and strong sink, respectively. Despite the discrepancies in the regional flux estimates between the three XCO2 inversions, the a posteriori CO2 distributions were in good agreement (with a mean difference between the three inversions of typically less than 0.5 ppm) with independent data from the Total Carbon Column Observing Network (TCCON), the surface flask network, and from the HIAPER Pole-to-Pole Observations (HIPPO) aircraft campaign. The discrepancy in the regional flux estimates from the different inversions, despite the agreement of the global flux estimates suggests the need for additional work to determine the minimum spatial scales at which we can reliably quantify the fluxes using GOSAT XCO2. The fact that the a posteriori CO2 from the different inversions were in good agreement with the independent data although the regional flux estimates differed significantly, suggests that innovative ways of exploiting existing data sets, and possibly additional observations, are needed to better evaluate the inferred regional flux estimates.

The seasonal variation of the CO2 flux over Tropical Asia estimated from GOSAT, CONTRAIL, and IASI

AbstractWe estimate the CO2 flux over Tropical Asia in 2009, 2010, and 2011 using Greenhouse Gases Observing Satellite (GOSAT) total column CO2(XCO2) and in situ measurements of CO2. Compared to flux estimates from assimilating surface measurements of CO2, GOSAT XCO2 estimates a more dynamic seasonal cycle and a large source in March–May 2010. The more dynamic seasonal cycle is consistent with earlier work by Patra et al. (2011), and the enhanced 2010 source is supported by independent upper air CO2 measurements from the Comprehensive Observation Network for Trace gases by Airliner (CONTRAIL) project. Using Infrared Atmospheric Sounding Interferometer (IASI) measurements of total column CO (XCO), we show that biomass burning CO2 can explain neither the dynamic seasonal cycle nor the 2010 source. We conclude that both features must come from the terrestrial biosphere. In particular, the 2010 source points to biosphere response to above‐average temperatures that year.

Validation of XCO<sub>2</sub> derived from SWIR spectra of GOSAT TANSO-FTS with aircraft measurement data

Abstract. Column-averaged dry air mole fractions of carbon dioxide (XCO2) retrieved from Greenhouse gases Observing SATellite (GOSAT) Short-Wavelength InfraRed (SWIR) observations were validated with aircraft measurements by the Comprehensive Observation Network for TRace gases by AIrLiner (CONTRAIL) project, the National Oceanic and Atmospheric Administration (NOAA), the US Department of Energy (DOE), the National Institute for Environmental Studies (NIES), the HIAPER Pole-to-Pole Observations (HIPPO) program, and the GOSAT validation aircraft observation campaign over Japan. To calculate XCO2 based on aircraft measurements (aircraft-based XCO2), tower measurements and model outputs were used for additional information near the surface and above the tropopause, respectively. Before validation, we investigated the impacts of GOSAT SWIR column averaging kernels (CAKs) and the shape of a priori profiles on the aircraft-based XCO2 calculation. The differences between aircraft-based XCO2 with and without the application of GOSAT CAK were evaluated to be less than ±0.4 ppm at most, and less than ±0.1 ppm on average. Therefore, we concluded that the GOSAT CAK produces only a minor effect on the aircraft-based XCO2 calculation in terms of the overall uncertainty of GOSAT XCO2. We compared GOSAT data retrieved within ±2 or ±5° latitude/longitude boxes centered at each aircraft measurement site to aircraft-based data measured on a GOSAT overpass day. The results indicated that GOSAT XCO2 over land regions agreed with aircraft-based XCO2, except that the former is biased by −0.68 ppm (−0.99 ppm) with a standard deviation of 2.56 ppm (2.51 ppm), whereas the averages of the differences between the GOSAT XCO2 over ocean and the aircraft-based XCO2 were −1.82 ppm (−2.27 ppm) with a standard deviation of 1.04 ppm (1.79 ppm) for ±2° (±5°) boxes.