The First Global Carbon Dioxide Flux Map Derived from TanSat Measurements

Dongxu YANG ,Yi LIU ,Liang FENG ,Jing WANG* ,Lu YAO ,Zhaonan CAI ,Sihong ZHU,Naimeng LU,and Daren LYU

1Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029,China

2National Centre for Earth Observation,University of Edinburgh,Edinburgh,EH9 3FF,UK

3University of Edinburgh,Edinburgh,EH9 3FF,UK

4National Satellite Meteorological Center,China Meteorological Administration Beijing 100081,China

ABSTRACT Space-borne measurements of atmospheric greenhouse gas concentrations provide global observation constraints for top-down estimates of surface carbon flux.Here,the first estimates of the global distribution of carbon surface fluxes inferred from dry-air CO2 column (XCO2) measurements by the Chinese Global Carbon Dioxide Monitoring Scientific Experimental Satellite (TanSat) are presented.An ensemble transform Kalman filter (ETKF) data assimilation system coupled with the GEOS-Chem global chemistry transport model is used to optimally fit model simulations with the TanSat XCO2 observations,which were retrieved using the Institute of Atmospheric Physics Carbon dioxide retrieval Algorithm for Satellite remote sensing (IAPCAS).High posterior error reduction (30%–50%) compared with a priori fluxes indicates that assimilating satellite XCO2 measurements provides highly effective constraints on global carbon flux estimation.Their impacts are also highlighted by significant spatiotemporal shifts in flux patterns over regions critical to the global carbon budget,such as tropical South America and China.An integrated global land carbon net flux of 6.71±0.76 Gt C yr−1 over 12 months (May 2017–April 2018) is estimated from the TanSat XCO2 data,which is generally consistent with other inversions based on satellite data,such as the JAXA GOSAT and NASA OCO-2 XCO2 retrievals.However,discrepancies were found in some regional flux estimates,particularly over the Southern Hemisphere,where there may still be uncorrected bias between satellite measurements due to the lack of independent reference observations.The results of this study provide the groundwork for further studies using current or future TanSat XCO2 data together with other surfacebased and space-borne measurements to quantify biosphere–atmosphere carbon exchange.

Key words:TanSat,carbon flux,CO2,flux inversion

1.Introduction

Top-down carbon flux inversion by assimilating atmospheric carbon dioxide (CO2) measurements is a useful tool for evaluating the global carbon budget and was used in the 2019 Refinement to the 2006 Intergovernmental Panel on Climate Change (IPCC) Guidelines for National Greenhouse Gas Inventories.Ground-based measurements provide highly accurate continuous data and have improved our understanding of global carbon flux.However,the sparseness and spatial inhomogeneity of the existing ground-based network limits our ability to infer consistent global-and regional-scale CO2sources and sinks (Scholes et al.,2009),especially for tropical and high-latitude regions.

To improve observation coverage,tailor-made satellites have been developed to provide atmospheric greenhouse gas(GHG) measurements at unprecedented precision.The Japanese GHG monitoring satellite mission,Greenhouse Gases Observing Satellite (GOSAT),was launched in 2009 (Kuze et al.,2009),and the U.S.satellite mission,Orbiting Carbon Observatory 2 (OCO-2),was launched in 2014 (Crisp et al.,2017).

Satellite measurements are advantageous for constraining top-down flux inversions because they provide continuous global observation coverage.Although satellite measurements are not as accurate as ground-based measurements,the increased spatial coverage can provide additional information not available from sparse surface networks (Buchwitz et al.,2007),and such data have been used in surface carbon flux inversion studies (Peylin et al.,2013;Saeki et al.,2013;Jiang et al.,2016;Chevallier et al.,2019;Wang et al.,2020).The CO2column-averaged dry-air mole fractions (XCO2) measured by GOSAT and OCO-2 have been applied in regional carbon flux optimization studies (Basu et al.,2013;Maksyutov et al.,2013;Feng et al.,2017;Chevallier et al.,2019;Wang et al.,2020).

Recent studies have shown that the accuracy of inversion results inferred from the GOSAT and OCO-2 retrievals is comparable to traditional inversions using accurate but sparse surface networks,indicating that satellite GHG observations provide valuable complementary data for global carbon budget studies (Chevallier et al.,2019;Wang et al.,2020).Assimilation of satellite GHG observations has significantly changed or raised questions about our understanding of spatiotemporal patterns of surface CO2fluxes.For example,GOSAT inversion data suggest high European biospheric uptake,nearly twice that suggested by in situ-only studies (Houweling et al.,2015;Feng et al.,2016;Reuter et al.,2017),which would modify the global carbon flux distribution (Houweling et al.,2015).Assimilated GOSAT and OCO-2 XCO2measurements have also revealed unexpectedly high net emissions from tropical Africa,which are thought to be mainly caused by substantial land-use change,resulting in the release of carbon from large soil organic carbon stores (Palmer et al.,2019).

The Chinese Global Carbon Dioxide Monitoring Scientific Experimental Satellite (TanSat),funded by the Ministry of Science and Technology of China,the Chinese Academy of Sciences,and the China Meteorological Administration,was launched in December 2016 (Liu and Yang,2016;Ran et al.,2019).The first global XCO2map measured by TanSat was reported in a previous study (Yang et al.,2018),and Total Column Carbon Observing Network validation shows 2.2 ppm accuracy for version 1 of the TanSat L2 data product (Liu et al.,2018).The accuracy and precision of TanSat XCO2retrievals were further improved using a wavelength dependence gain factor to correct the spectrum continuum feature (Yang et al.,2020).A new version of TanSat XCO2was recently released to the public (Yang et al.,2021),which provides an opportunity to improve our knowledge of global carbon flux from flux inversions based on TanSat measurements.

In this study,we introduce the first estimate of global carbon flux based on TanSat global measurements.The following section shows the TanSat measurements and methods used in the flux inversion.Section 3 introduces the main results of the optimized carbon flux.We then close with our conclusions.

2.Method

2.1.TanSat measurement

TanSat flies in a sun-synchronous low Earth orbit (LEO) with an equator crossing time around 1330 local time and operates in three observation modes including nadir (ND),glint (GL),and target modes.ND mode is used when TanSat flies over land surfaces,and GL mode is activated over the ocean to increase the incident signal and ensure the signal-to-noise ratios meet the requirements.The atmospheric carbon dioxide grating spectrometer on board TanSat provides hyperspectral measurements of the O2A band (0.76 μm) and CO2weak (1.61 μm) and strong (2.04 μm) bands.The ND footprint size of TanSat measurement is～2 × 2 km with nine footprints across the track in a frame,while the total field-of-view width is～18 km.The TanSat repeat measurement cycle is 16 days,and the ND and GL measurements are staggered.In ND mode,TanSat has the same ground track interval as OCO-2.

XCO2retrieval is performed using the Institute of Atmospheric Physics carbon dioxide retrieval algorithm for satellite remote sensing (Liu et al.,2013;Yang et al.,2015).The version 1 TanSat XCO2retrievals only used the CO2weak band due to calibration issues on the O2A band (which can cause critically biased estimates of surface pressure).Therefore,the flux inversion precision and accuracy were not adequate for studying global surface carbon flux.Further study on solar calibration measurements revealed a remaining feature on the solar spectrum fitting residual.A Fourier series model was applied as a gain factor to the continuum in the retrieval,significantly improving the fitting residual,especially on the O2A band.The O2A band and CO2weak band measurements have been used together in new retrievals that correct the parameters of water vapor,surface pressure,temperature,aerosol,cirrus,and instrument model.A genetic algorithm was used in post screening,and then a multiple linear regression model was applied to correct bias (Yang et al.,2020).In this study,we used a 15-month (March 2017–May 2018) XCO2retrieval from TanSat ND measurements.

Only “good” retrievals have been provided in the TanSat v2 product used in flux inversion.We used single sounding uncertainty (posterior error) to construct weighted time–space average data (Crowell et al.,2019).First,we calculated a 1-second span average:

whereXCO2,iis the ith sounding value in a 1-second span with a total of N soundings,and σiis the corresponding posterior error.The summary includes all soundings in a 1-second span,and the uncertainty for the 1-second average is:

where we assumed soundings in a 1-second span were highly correlated.Because the posterior error mainly considered the theoretical error from the measurement noise level,theestimation was sometimes lower than real measurement uncertainty.Therefore,we also considered the standard deviation;the uncertainty of a 1-second span,,was the maximum of either.This differs from the 1-second spans of OCO-2 because TanSat L2 sounding values have not been corrected for the bias of small areas.Thus,small systematic uncertainty,such as floor error,was not considered.

In this study,we used a 5-second span average for flux inversion.Thevalues were constructed using the methods for 1-second spans of sounding measurements.

2.2.Carbon flux inversion system

We used an ensemble transform Kalman filter (ETKF) data assimilation system (Feng et al.,2009,2011,2016) to estimate global and regional carbon flux.The carbon flux and CO2concentration were optimized by TanSat XCO2measurements via:

wherexaandxfare the a posteriori and a priori fluxes,respectively.Hoois the observation operator that describes the relationship between the state vector and the observations.yobsis XCO2measured by TanSat,and K is the Kalman gain matrix,which determines the adjustment to the a priori based on the difference between the model and observations and their uncertainties as follows:

where R is the observation error covariance matrix,namely,a diagonal matrix representing measurement errors,including instrument error,retrieval error,model error,and representation error (Peylin et al.,2002).Pfis the a priori error covariance matrix.The Jacobian of the observation operator,H,mapsPfinto observation space.In our carbon flux inversion system,we used the chemistry transport model (CTM) GEOS-Chem (v9-02) in forward simulationand performed an ensemble run to establish the connection between surface fluxes and atmospheric CO2concentrations (H).

We used sequentially assimilated TanSat 5-second averaged measurements day by day in each assimilation step (1 month) to optimize the a priori estimation of surface CO2fluxes.The GEOS-5 meteorological analyses,provided by the NASA Goddard Global Modelling and Assimilation Office,were used to drive the GEOS-Chem run.The model was run at a horizontal resolution of 4° (latitude) × 5° (longitude) with 47 vertical levels,which spanned from the surface to the mesosphere,typically with 35 levels in the troposphere.Monthly inventories have been used as the a priori flux for GEOS-Chem runs,including ocean flux (Takahashi et al.,2009),biomass burning fluxes (van der Werf et al.,2010),3-hourly terrestrial biosphere fluxes (Olsen,2004),and fossil fuel emissions (Oda and Maksyutov,2011).For the terrestrial biosphere,we used a 3-hour flux to better represent and optimize the uptakes and emissions from photosynthesis and respiration processes,respectively (Olsen,2004).The inversion for surface flux spanned 792 regions globally,including 475 land regions and 317 ocean regions,based on a TransCom 3 study (Gurney et al.,2002).

3.Results

3.1.Annually integrated global carbon net flux

In this study,only the ND model land observations were applied to constrain the flux optimization.Due to the low signal-to-noise ratio,the ND measurements over the ocean are not provided in the L2 XCO2data product,and no further information is provided by TanSat ocean measurements.Therefore,a very strong prior constraint was applied to the ocean CO2flux.The inversion result is shown as carbon flux based on CO2flux using the mass ratio between CO2and carbon.

Error reduction for global carbon sinks indicates how well the measurements optimized the estimations.The error reduction in July was greater than in other months (Fig.1);for example,it was >50% in northern Asia,Europe,and America,but lower in Australia,southern Asia,and India.South America and Africa showed 20%–40% error reduction.TanSat measurements improved results in South Africa mostly in April and October,in North and South America mostly in July,and in northern Eurasia in all three months.

Fig.1.The reductions in uncertainty for (a) July 2017,(b) October 2017,(c) January 2018,and (d) April 2018.

The global annual carbon sinks (Fig.2) constrained by TanSat measurements differed significantly from a priori fluxes.For example,the land sink increased in most of Northern Asia,Europe,and the Americas,but decreased in the middle of Africa and India.Similar results have been reported in GOSAT and OCO-2 CO2flux inversion studies (Wang et al.,2019,2020).The land sink also increased in southwest China,which is consistent with carbon flux estimations based on additional ground-based in-situ measurements from Asia (Wang et al.,2020).The optimization was not significant over the ocean because no direct measurements were used to constrain the oceanic carbon flux.

Our TanSat inversion indicated a land net flux integrated over 12 months (May 2017–April 2018) of 6.71 Gt C yr−1with an uncertainty of 0.76 Gt C yr−1.Table 1shows the land carbon net flux over the same period constrained by the new L2 OCO-2 XCO2 data product (v10) as a reference.The annul integrated OCO-2 net carbon flux was 5.13±0.59 Gt C yr−1,which is 23.5% lower than the TanSat results.Several studies have reported different carbon flux estimations based on satellite measurements for 2009 onwards.The GOSAT measurement was the first spaced-based GHG measurement used in global CO2flux estimation (Basu et al.,2013).GOSAT measurements indicated land net carbon fluxes estimated by the carbon tracker (Peters et al.,2007;https://www.esrl.noaa.gov/gmd/ccgg/carbontracker/) of 7.86±1.26 Gt C yr−1in 2017 and 7.61±1.78 Gt C yr−1in 2018.The OCO-2 model intercomparison project showed 10 inversion model parallel studies with results ranging from 5.75 to 8.97 Gt C yr−1from June 2017 to May 2018.

Table 1.The a priori and a posteriori fluxes estimated by TanSat (v2) and OCO-2 (v10) for 11 TransCom regions and continents.

3.2.Seasonal and regional CO2 flux

Seasonal net carbon flux was grouped every three months and denoted following the Northern Hemisphere seasons,namely,summer (June,July,August),autumn (September,October,November),winter (December,January,February),and spring (March,April,May).The distribution of global carbon flux differed significantly from a priori flux after assimilation of TanSat measurements (Fig.3).For example,the measurements (satellite data) significantly changed carbon fluxes for South America for the whole year (Fig.4).The South America tropical region showed the largest impacts from TanSat measurements,with the seasonal phase shifted significantly from a priori fluxes.We found that the satellite did not provide dense measurements over tropical South America,but even a small amount of TanSat data can still significantly change the carbon flux feature.TanSat measurements increased carbon sinks in spring and summer around Eurasia,including in Europe and boreal and temperate Eurasia,compared with a priori fluxes.North Africa also showed significant increases in the carbon sink in winter.In addition,the seasonal variation of India decreased during spring and summer.

We also observed changes in China’s carbon flux in autumn,winter,and spring;in general,the carbon sink increased compared to a priori fluxes,especially in southwest China in the summer and autumn.Similar results have been reported by studies of vegetation trend measurements (Chen et al.,2019) and flux estimations (Wang et al.,2020).

Figure 5shows regional carbon flux from OCO-2 and TanSat measurements,as well as a priori fluxes.OCO-2 and TanSat have similar optimization effects on a priori fluxes in most regions.The biggest discrepancy was in the Southern Hemisphere,where TanSat measurements and a priori fluxes indicated a source,whereas OCO-2 measurements indicated a sink;the difference in outcomes was driven in large part by results from tropical and temperate South America.This contradiction has also been reported when comparing results from GOSAT and OCO-2 measurements (Wang et al.,2019).

Fig.2.The annual carbon sink optimized by GOSAT measurements (b) from a priori estimates (a).The color bar indicates the carbon sink.

Fig.3.The a priori (upper row) and a posteriori (lower row) seasonal carbon sinks constrained by TanSat measurements.The columns indicate summer,autumn,winter,and spring from left to right.

Fig.3.(Continued).

Fig.3.(Continued).

4.Conclusions and outlook

We described the first attempt at estimating global carbon flux using TanSat XCO2measurements.We used an ETKF data assimilation system coupled with the GEOS-Chem CTM in a top-down carbon flux inversion.The error reduction compared to a priori estimates demonstrated the impact of assimilating global satellite XCO2measurements,which resulted in significant changes in spatiotemporal flux patterns throughout the year from May 2017 to April 2018.Our TanSat inversion indicated an annual integrated global land carbon net flux of 6.71±0.76 Gt C yr−1,compared to annual totals of 5.75–8.97 Gt C yr−1for inversions based on OCO-2 XCO2data.

These carbon flux data are publicly accessible worldwide on the China GEO TanSat data service archive,the Cooperation on the Analysis of Carbon Satellites Data (CASA;www.chinageoss.org/tansat),which is hosted by the International Reanalysis Cooperation on Carbon Satellites Data (IRCSD).

TanSat provides preliminary global measurements that extend the ground-based network from local measurements to global estimates.However,the coverage,repeat,pixel size,and accuracy of measurements are still inadequate to meet the Global Stocktake of the Paris Agreement and inventory validation requirements of IPCC standards.For example,GOSAT measurements have indicated that flux estimation is sensitive to measurement coverage (Deng et al.,2014).XCO2measurements are also sensitive to regional flux inversion (Feng et al.,2009;Palmer et al.,2011).The future European Copernicus Carbon Dioxide Monitoring satellite mission (Kuhlmann et al.,2019) will improve measurements of global coverage and investigations of anthropogenic emissions.The plan for the next TanSat mission has been discussed and is in the design phase.One possibility is a multi-satellite LEO constellation.Several satellites would capture global background measurements with large pixel size and a wide swath,while others would focus on target measurements for city and point sources with a small pixel size.This would capture carbon sinks and sources at multiple scales to elucidate the influences of natural processes and anthropogenic activities on the carbon cycle.

Fig.4.The a priori and a posteriori CO2 flux trends for 11 TransCom land regions shown with measurement quantities.The names at the top of each subplot indicate regions.

Column CO2measurements are still inadequate for separating local carbon exchanges from long-range transport processes because the signal is mixed into the column value (Keppel-Aleks et al.,2011).The signal of local flux is critical to city and point-source research.Ground-based measurements near the surface will be helpful for future investigations.Therefore,coordinating satellites and surface measurement networks will significantly improve the study of small-scale anthropogenic emission sources.

Fig.5.The carbon flux estimated by TanSat and OCO-2 measurements and the a priori flux in 11 TransCom land regions,South Land (SL),North Land (NL),and globally.

Acknowledgements.This work is supported by the National Key R&D Program of China (Grant No.2016YFA0600203),the Key Research Program of the Chinese Academy of Sciences (ZDRW-ZS-2019-1),the National Key R&D Program of China (Grant No.2017YFB0504000),and the Youth Program of the National Natural Science Foundation of China (Grant No.41905029).Liang FENG is supported by the UK NERC National Centre for Earth Observation (NCEO).The TanSat L1B data service is provided by IRCSD and CASA (131211KYSB20180002).We also thank the FENGYUN Satellite Data Center of the National Satellite Meteorological Center,who provided the TanSat L1B data service.The authors thank the TanSat mission and highly appreciate the support from everyone involved.