Detection of Anthropogenic CO2 Emission Signatures with TanSat CO2 and with Copernicus Sentinel-5 Precursor (S5P) NO2 Measurements: First Results

Dongxu YANG, Janne HAKKARAINEN, Yi LIU, Iolanda IALONGO,Zhaonan CAI, and Johanna TAMMINEN

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

2Finnish Meteorological Institute, Helsinki FI-00560, Finland

ABSTRACT China’s first carbon dioxide (CO2) measurement satellite mission, TanSat, was launched in December 2016.This paper introduces the first attempt to detect anthropogenic CO2 emission signatures using CO2 observations from TanSat and NO2 measurements from the TROPOspheric Monitoring Instrument (TROPOMI) onboard the Copernicus Sentinel-5 Precursor (S5P) satellite.We focus our analysis on two selected cases in Tangshan, China and Tokyo, Japan.We found that the TanSat XCO2 measurements have the capability to capture the anthropogenic variations in the plume and have spatial patterns similar to that of the TROPOMI NO2 observations.The linear fit between TanSat XCO2 and TROPOMI NO2 indicates the CO2-to-NO2 ratio of 0.8 × 10-16 ppm (molec cm-2)-1 in Tangshan and 2.3 × 10-16 ppm (molec cm-2)-1 in Tokyo.Our results align with the CO2-to-NOx emission ratios obtained from the EDGAR v6 emission inventory.

Key words: TanSat, CO2, Remote sensing, city carbon emission, climate change

1.Introduction

By absorbing infrared radiation, greenhouse gases trap heat in the Earth’s atmosphere, which would otherwise escape into space, and greenhouse gas concentrations in the atmosphere have increased since the onset of the industrial revolution.In the 2015 United Nations Climate Change Conference held in Paris, participants agreed to reduce greenhouse gas emissions to prevent an average global surface temperature increase of more than 1.5°C.Hence, mitigating or slowing down global warming represents a challenge faced by the global population in the 21st century.Concentrations of carbon dioxide (CO2)have risen by more than 40% due to anthropogenic activities such as fossil fuel combustion and land-use change.The emissions related to the combustion of fossil fuels are particularly localized, with urban areas being the dominant contributor responsible for more than 70% of global emissions.

Even though global economies are struggling due, in part, to the recovery from the impact of COVID-19, in 2021, fossilfuel-related CO2emissions have shown an increase of 5% with an uncertainty of ±5% (Friedlingstein et al., 2022).Large uncertainty remains in fossil-fuel emission estimates when focused on specific countries (Andres et al., 2014), and uncertainties are even larger than expected for cities (Ciais et al., 2014; Gately et al., 2015; Miller and Michalak, 2017; Gately and Hutyra, 2017; Oda et al., 2018; Gurney et al., 2019).For example, one investigation indicates a 20.4% variation in the fossilfuel emission estimates of Los Angeles, while variations exceeding 200% were found in Beijing (Han et al., 2020).The uncertainty mainly comes from two parts, the remaining bias in the emission inventory statistics, for example, less than optimal estimations of energy consumption and emission factors (Liu et al., 2015), and the inability to adequately capture the rapidly increasing magnitude of fossil-fuel consumption related to economic growth, which is not updated in a timely manner.

To better understand the uncertainty of carbon emissions (sources) and sinks, satellite missions dedicated to atmospheric greenhouse measurements have been developed in the last decade.China’s Global Carbon Dioxide Monitoring Scientific Experimental Satellite (TanSat) was launched in December 2016 after Japan’s Greenhouse gases Observing SATellite(GOSAT, 2009) (Kuze et al., 2009) and NASA’s Orbiting Carbon Observatory-2 (OCO-2, 2014) (Crisp et al., 2017) mission.Similar to the OCO-2 mission, TanSat has a hyperspectral grating spectrometer covering the near-infrared (NIR) and short-wave infrared (SWIR) regions that are sensitive to CO2absorption (Liu and Yang, 2016).The nadir footprint size of TanSat measurement is ～2 × 2 km2, and it has nine footprints across the track in a frame, which means TanSat has the potential to capture anthropogenic emission enhancements in a city and downwind of a power station.

The CEOS Chair commissioned the Atmospheric Composition Virtual Constellation (AC-VC) to define a global architecture for monitoring atmospheric greenhouse gas concentrations and their fluxes from space (Crisp et al., 2018).Quantifying anthropogenic CO2emissions is one of the most important requirements of the global greenhouse gas monitoring system that increases the verification support (MVS) capacity to support the Paris Agreement.Carbon dioxide measurements alone cannot provide enough information to screen out the atmospheric background, which includes CO2turbulence from the natural carbon cycle; therefore, it is useful to coordinately use measurements of short-lived anthropogenic tracers, e.g., NO2and CO.In this study, we introduce the first attempts to use TanSat CO2measurements to detect anthropogenic CO2emission signatures together with the Copernicus Sentinel-5 Precursor (S5P) NO2measurements.

2.Satellite data

2.1.TanSat XCO2 retrieval

TanSat XCO2data are retrieved by the Institute of Atmospheric Physics Carbon dioxide retrieval Algorithm for Satellite remote sensing (IAPCAS) (Yang et al., 2018, 2020, 2021).The standard TanSat XCO2retrieval has been introduced and released to the public (http://www.chinageoss.cn/tansat/index.html).The cross footprint bias significantly impacts research that compares city-based emissions to global carbon flux inversions.The bias correction method used in the TanSat v2 data product cannot meet the requirements of this study due to the bias correction reference to the global absolute bias.The relative biases between a footprint with a scale of tens of kilometers to a footprint with a scale of hundreds of kilometers have not been efficiently resolved.In this study, we introduce a new bias correction strategy.To reduce a cross-footprint bias, we first apply a dynamic bias-related parameter selection to a cloud-screened data product.The bias has been separated into absolute bias and relative bias components.The former part has been corrected to a TCCON measurement which uses a coupling method similar to TanSat v2 processing.For each footprint, the relative bias correction parameters are selected from a rank of bias relativity, including solar zenith angle, surface pressure, albedo, aerosol optical depth, cloud optical depth, spectrum correction coefficient, etc.Additionally, in the relative bias correction, we apply a 0.1° moving window cross-footprint bias correction for each footprint.This special bias-corrected XCO2data product is only used in this study, and the bias correction method will be considered in the next version of TanSat XCO2data production [Fig.S1 in the electronic supplementary material (ESM)].

2.2.TROPOMI/Sentinel-5P NO2 data

The TROPOspheric Monitoring Instrument (TROPOMI) is the sole satellite instrument onboard the Copernicus Sentinel-5 Precursor (S5P) satellite.It flies on a sun-synchronous polar orbit with an equatorial overpass time at 1330 LST(LST = UTC + 8), close to the TanSat mission equatorial overpass time.One of the advantages of the TROPOMI tropospheric NO2observations compared to its predecessors is the high spatial resolution (currently 5.5 km by 3.5 km at nadir) which enables the detection of NO2plumes.As nitrogen oxides (NOx= NO + NO2) are often co-emitted with CO2, this data can be used to facilitate the identification and analysis of the CO2plumes.For example, TROPOMI NO2data were used to derive the plume width of the CO2plume using the cross-sectional flux method (Reuter et al., 2019) and to analyze the NOxto-CO2emission ratio together with OCO-2 XCO2observations of the Matimba power station in South Africa (Hakkarainen et al., 2021).Previously, tropospheric NO2observations had been used, for example, to characterize global XCO2anomalies(Hakkarainen et al., 2016).

3.Results

Here, we analyze two TanSat overpasses over anthropogenic emission plumes.In the first case (6 May 2018), TanSat captures an XCO2enhancement near Tangshan, China, while in the second case (29 March 2018), it detects an emission plume originating from Tokyo, Japan.

Tangshan:Figure 1 illustrates the case of Tangshan on 6 May 2018.The map in Fig.1 (left panel) shows the TROPOMI tropospheric NO2columns overlapped by the TanSat XCO2observations.The arrows indicate the ERA5 wind fields.For illustrative purposes, we removed the local background from the TanSat XCO2observations.The background iscalculated from the median of the XCO2values north of the white dashed line in Fig.1, and we thus consider this to be the local XCO2anomaly.With respect to this background, the largest enhancements near Tangshan are about 4-5 ppm.The colocated TROPOMI tropospheric NO2columns indicate values of 3.5 × 1016molec cm-2.Figure 1 (right, top, and bottom panels, respectively) shows the NO2and XCO2observations (blue dots) along the TanSat track, while the black dots indicate the running means.We observe two enhancements at about 39.5°N and 40.1°N from both datasets, corresponding to one emission area near Tangshan and another about 60 km north).We calculated the linear fit between TROPOMI NO2and TanSat XCO2observations to transform NO2tropospheric columns into XCO2(orange dots in Fig.1; top right panel).The slope of the linear fit indicates the CO2-to-NO2ratio (Hakkarainen et al., 2019, 2021) for this case is about 0.8 × 10-16ppm(molec cm-2)-1, roughly corresponding to the ratio of the CO2and NO2peak values near Tangshan.The correlation coefficient of the linear fit is 0.54.

Fig.1.The left panel illustrates TROPOMI NO2 observations overlapped by TanSat XCO2 observation near Tangshan(China) on 6 May 2018.ERA5 wind fields are shown as arrows.The XCO2 background is calculated as a median, north of the white dashed line.The right panels (top and bottom) show the XCO2 and NO2 observations (respectively) along the TanSat track.The black dots indicate the running mean.Orange colors indicate the XCO2 converted from NO2 tropospheric columns via a linear fit.

Tokyo:Figure 2 shows the case of the Tokyo plume on 29 March 2018.We again calculate the background north of the white dashed line and consider local XCO2anomalies.We note that the enhancements are spread over a wide city area.The highest running mean values are between 2-3 ppm.The highest TROPOMI NO2enhancements co-located with TanSat are about 1 × 1016molec cm-2.In addition to the main enhancement east of Tokyo, another smaller enhancement is visible,e.g., near latitude 36.5°N.The linear fit between TanSat XCO2and TROPOMI NO2indicates that the CO2-to-NO2ratio is about 2.3 × 10-16ppm (molec cm-2)-1, which roughly corresponds to the ratio of the peak values.This ratio is about 2-3 times higher than the ratio observed near Tangshan.The correlation coefficient of the linear fit is 0.47.The main reason for the higher CO2-to-NO2ratio is the higher CO2-to-NOxemission ratio in Tokyo with respect to Tangshan.For example, considering the pixels ±0.5° around Tangshan and Tokyo in the EDGAR v6 emission inventory (Crippa et al., 2021) gives an emission ratio of about 400 for Tangshan and 1000 for Tokyo.In addition, the higher ratio in Tokyo could be partly due to the distance from the source, as the NO2lifetime is relatively short (about 4 hours).However, in both cases, the linear fit is calculated over a large area.

4.Outlook

In this study, we introduced the first attempt to use TanSat measurements to detect anthropogenic CO2emission signatures.Carbon dioxide enhancements of 2-3 ppm in Tokyo and 3-4 ppm in Tangshan have been observed by TanSat nadir measurements, corresponding to enhancements in the S5P NO2measurement.To accurately measure anthropogenic CO2signals, the coverage, footprint size, and repeat cycle of satellite observations need to be improved.The next generation of China’s Global Carbon Dioxide Monitoring Satellite mission, TanSat-2, is now in the design phase.The target measurement will focus on cities with an 800-1000 km swath to record the gradient of XCO2from its central to rural areas using an imaging process and a 500 m footprint size to improve the emission estimation accuracy.TanSat-2 will not be a single satellite but a constellation distributed into at least two orbits in the morning and afternoon to cover a city or point source twice daily.Theco-located measurements of NO2help constrain anthropogenic CO2emissions from the total CO2budget.Hence, a NO2-detecting instrument will be on board the TanSat-2 satellite to identify the emission plume and help constrain anthropogenic emissions in the atmospheric inversion.

Fig.2.Same as Fig.1, but for Tokyo on 29 March 2018.

Acknowledgements.This work is supported by the National Key Research And Development Plan (2019YFE0127500), International Partnership Program of the Chinese Academy of Sciences (060GJHZ2022070MI).The authors thank the Finland-China mobility cooperation project funded by the Academy of Finland (No.348596) and the Key Research Program of the Chinese Academy of Sciences (ZDRWZS-2019-1).Financial support for the Academy of Finland (No.336798) is kindly acknowledged.The authors thank the TanSat mission,and the support from everyone that worked with TanSat mission is highly appreciated.

Electronic supplementary material:Supplementary material is available in the online version of this article at https://doi.org/10.1007/s00376-022-2237-5.