Articles | Volume 22, issue 22
Research article
16 Nov 2022
Research article |  | 16 Nov 2022

Towards sector-based attribution using intra-city variations in satellite-based emission ratios between CO2 and CO

Dien Wu, Junjie Liu, Paul O. Wennberg, Paul I. Palmer, Robert R. Nelson, Matthäus Kiel, and Annmarie Eldering

Carbon dioxide (CO2) and air pollutants such as carbon monoxide (CO) are co-emitted by many combustion sources. Previous efforts have combined satellite-based observations of multiple tracers to calculate their emission ratio (ER) for inferring combustion efficiency at the regional to city scale. Very few studies have focused on combustion efficiency at the sub-city scale or related it to emission sectors using space-based observations. Several factors are important for interpreting and deriving spatially resolved ERs from asynchronous satellite measurements, including (1) variations in meteorological conditions given the mismatch in satellite overpass times, (2) differences in vertical sensitivity of the retrievals (i.e., averaging kernel profiles), (3) interferences from the biosphere and biomass burning, and (4) the mismatch in the daytime variations of CO and CO2 emissions. In this study, we extended an established emission estimate approach to arrive at spatially resolved ERs based on retrieved column-averaged CO2 (XCO2) from the Snapshot Area Mapping (SAM) mode of the Orbiting Carbon Observatory-3 (OCO-3) and column-averaged CO from the TROPOspheric Monitoring Instrument (TROPOMI).

To evaluate the influences of the confounding factors listed above and further attribute intra-urban variations in ERs to certain sectors, we leveraged a Lagrangian atmospheric transport model with an urban land cover classification dataset and reported ERCO values from the sounding level to the overpass and city level. We found that the differences in overpass times and averaging kernels between OCO and TROPOMI strongly affect the estimated spatially resolved ERCO. Specifically, a time difference of >3 h typically led to dramatic changes in wind directions and urban plume shapes, thereby making the calculation of accurate sounding-specific ERCO difficult. After removing such cases from consideration and applying a simple plume shift method when necessary to account for changes in wind direction and speed, we discovered significant contrasts in combustion efficiencies between (1) two megacities versus two industry-oriented cities and (2) different regions within a city, based on six nearly coincident overpasses per city. Results suggest that the ERCO impacted by heavy industry in Los Angeles is slightly lower than the overall city-wide value (<10 ppb-CO/ppm-CO2). In contrast, the ERCO related to heavy industry in Shanghai is much higher than Shanghai's city mean and more aligned with the city means of two selected industry-oriented cities in China (approaching 20 ppb-CO/ppm-CO2). Although investigations based on a larger number of satellite overpasses are needed, our unique approach (i.e., without using sector-specific information from emission inventories) provides new insights into assessing combustion efficiency within a city from future satellite missions, such as those that will map column CO2 and CO concentrations simultaneously with high spatiotemporal resolutions.

1 Introduction

Home to more than half of the total global population, urban areas have been expanding, especially in Asia and Africa, which had urbanization rates of 1.3 % and 1.1 % yr−1, respectively, between 2015 and 2020 (United Nations et al.2020). Urban regions are also responsible for a significant amount of anthropogenic emissions of greenhouse gases (GHG) and air pollutants into the atmosphere, including carbon dioxide (CO2), methane, carbon monoxide (CO), and nitrogen oxides (Duncan et al.2016; Lin et al.2018; Super et al.2017; Plant et al.2019). Satellite observations have become indispensable for monitoring the abundances of several atmospheric species in a globally consistent manner (Yokota et al.2009; Crisp et al.2012; Veefkind et al.2012). For example, carbon-monitoring satellites such as Orbiting Carbon Observatory-2 (OCO-2, Crisp et al.2012) have made the quantification of city-scale CO2 emissions and emission trends possible (e.g., Hedelius et al.2018; Ye et al.2020; Wu et al.2020; Shekhar et al.2020; Lei et al.2021). Quantifying the spatial gradient of atmospheric concentrations and relating this gradient to emissions within the city domain has become the next critical yet challenging task. Understanding such spatial heterogeneity in emissions and the environmental consequences of it can support better decisions in urban planning and the pinpointing of hotspots for emission mitigation.

Given the co-benefit of GHG reduction and improved air quality at various scales (Zhang et al.2017), controlling the consumption of fossil fuels altogether is the key. The efficiencies associated with various combustion activities are linked to their underlying processes and conditions (e.g., oxygen-to-fuel ratio and temperature). For example, the amount of CO2 emitted from coal-fired power plants varies with thermal and pressure conditions, the type of fuel consumed, the technology deployed, and the service duration of power plants (Yuan and Smith2011). A modern power generation plant that uses scrubbing technology is often regarded as a “clean” emitter, leading to minimal CO and NOx enhancement (Lindenmaier et al.2014). The commonly used approach when estimating combustion efficiency is to combine atmospheric observations of multiple trace gases and report the ratio of the total or excess measured concentrations (above a defined background value) of the tracers (Silva and Arellano2017; Reuter et al.2019; Park et al.2021). Such a tracer-to-tracer ratio calculation has the benefit that errors in describing the atmospheric transport that carries tracers to the measurement site can be canceled. A few notable studies have further utilized derived emission ratios (ERs) from ground or airborne measurements to infer sector-specific emission signals (Wennberg et al.2012; Lindenmaier et al.2014; Nathan et al.2018; Tang et al.2020).

CO and NOx often serve as tracers for anthropogenic CO2 as they arise from similar sources (e.g., Palmer et al.2006; Wunch et al.2009; Hedelius et al.2018). Analyzing remotely sensed NOx plumes with relatively short lifetimes can help identify local fossil fuel CO2 (FFCO2) sources that would otherwise be difficult to detect (Reuter et al.2019; Fujinawa et al.2021). At the same time, such reactivity requires that chemical transformations are accurately accounted for and complicates the interpretation of emission signals or ERs from NOx observations (Lama et al.2020; Hakkarainen et al.2021). Given its much longer lifetime, CO is much easier to interpret and more likely to be found during incomplete combustion. The emission ratio of CO to CO2 (ERCO) is usually estimated from sparse ground-based measurements within a city (Bares et al.2018; Chandra et al.2016; Lindenmaier et al.2014) and from satellites at the city scale (Park et al.2021; Silva and Arellano2017). Sector-specific activities and ERCO values such as those from the traffic sector have been analyzed in limited but valuable tunnel studies (Ammoura et al.2014; Bradley et al.2000; Popa et al.2014).

Figure 1ERCO [ppb ppm−1] values associated with specific processes (a) and ERCO values integrated over the entire city/region (b), as summarized from previous studies. The x axis indicates the times at which these estimates were made, except for Akagi et al. (2011); in the latter case, 2011 was chosen for the x axis since the paper was published in 2011. The error bars represent the uncertainties in the estimated ERs and the shaded rectangles indicate ranges of ERs over multiple years. Paper citations are omitted from the figure but are included in Appendix A. ERs related to the biomass burning and shipping sectors are derived using EFCO and EFCO2. The range of overpass-specific ERCO estimates for Shanghai, LA, Baotou, and Zibo derived from our study are also shown in the figure as a dashed black box.


We performed a literature search for ERCO values derived from observations (Appendix A) and the values are summarized in Fig. 1. The combustion efficiency fluctuates (1) over time (e.g., Turnbull et al.2011b), likely due to technological improvements, and (2) between sub-sectors, e.g., gasoline vs. diesel vehicles or moving vs. congestion traffic (Westerdahl et al.2009; Popa et al.2014). Despite differences in measurement platforms and analytical approaches, the observed urban-integrated ERCO values, especially those in Europe and the United States, are well constrained within the range of 4 to 15 ppb ppm−1 (Fig. 1b). ERCO values for biomass burning and shipping sectors are estimated based on fuel-specific emission factors, i.e., ERCO (= EFCO/EFCO2) with proper unit conversions, where each emission factor EFX indicates the emission of gas X per kg of fuel burned.

When estimating fossil fuel emissions from a bottom-up perspective, most inventories rely on activity data and may involve prior knowledge of emission factors (Gurney et al.2019; Solazzo et al.2021). One notable example is Hestia, a high-resolution inventory for the US that estimates CO2 emissions of non-point sources based on CO emissions from the National Emission Inventory and EFs and carefully evaluates the adopted EFs (Gurney et al.2019). However, when constructing emission inventories across regions and nations, the large variability in ERs across combustion processes, sectors, years, and regions (as seen in Fig. 1a) makes the choice of EFs extremely challenging. Accurate bottom-up emission estimates require accurate activity data and EFX values that naturally vary with the combustion conditions (e.g., temperature, fuel load, oxygen level) and are generally not well known, especially over data-scarce regions. To our knowledge, only a few global inventories, such as the Emissions Database for Global Atmospheric Research (EDGAR, Solazzo et al.2021), offer global anthropogenic CO and CO2 emissions. Considering the challenge involved in approximating ERs, certain knowledge derived from atmospheric observations may (1) complement inventory-based ERs (e.g., the CO:NOx ratio in Lama et al.2020) and (2) facilitate emission constraints for a desired gas with relatively large uncertainties (Wunch et al.2009; Palmer et al.2006; Wang et al.2009; Brioude et al.2012; Nathan et al.2018). Such prior achievements motivate us to examine ERs using satellite observations of multiple tracers.

Most existing studies have focused on quantifying an integrated ER for a whole city or region. We take a step forward, zooming into an urban area and leveraging spatially resolved satellite observations. Intra-city variations in the satellite-based concentration of a specific air pollutant such as NOx have been analyzed and linked to societal inequalities regarding income and educational attainment (Demetillo et al.2021; Kerr et al.2021). Yet, no one has attempted to study the intra-urban gradient in combustion efficiency from space and to relate this gradient to a specific combustion sector. This is now possible by virtue of Orbiting Carbon Observatory-3 (OCO-3, mounted on the International Space Station), which can sample a city landscape during the Snapshot Area Mapping (SAM) mode (Eldering et al.2019; Taylor et al.2020; Kiel et al.2021). In an effort to arrive at spatially varying ERs from sensors with asynchronous orbits, we must account for several factors that have not been thoroughly investigated. These include (a) variations in meteorological conditions and emission patterns during different overpass times, (b) discrepancies in the vertical sensitivities of the retrievals (i.e., averaging kernel (AK) profiles), and (c) interference from non-anthropogenic sources and sinks, especially from the biosphere.

In this study, we explore the spatial distribution of ERCO within four urban areas, mainly using XCO2 observations from OCO-3 and XCO observations from the TROPOspheric Monitoring Instrument onboard the Sentinel-5 Precursor (TROPOMI; Veefkind et al.2012). To avoid relying on prior sector-specific information on ERCO from emission inventories, we adopt the urban land cover data from the high-resolution World Urban Database and Access Portal Tools (WUDAPT; Ching et al.2018). WUDAPT offers the so-called local climate zone (LCZ) that considers the building structure/spacing along with the vegetation coverage (Stewart and Oke2012), which shed light on the urban infrastructure.

Our work seeks to answer the following two questions:

  1. Is it possible to accurately quantify the spatially resolved ERCO from asynchronous satellite measurements?

  2. Can the combustion efficiency for a given sector be determined without using prior emission inventories?

In Sect. 2, we describe the satellite data and methodology used for obtaining emissions, the ERCO, and associated uncertainties. In Sect. 3, we show intra-city variations in ERCO (including the ERCO tied to heavy industry in a megacity) and how multiple factors may interfere when deriving ERCO. In Sect. 4, we discuss the implications and limitations of this analysis.

2 Data and methodology

We target two types of cities: (1) industry- and energy-oriented cities (Baotou, China and Zibo, China) and (2) megacities with more diverse emission sectors (Shanghai, China and Los Angeles, USA). The four cities are selected considering the amount and quality of XCO2 data from OCO-3 SAMs and TROPOMI XCO data. The two industry- and energy-oriented cities in China are selected given their large amount of metal production plants for aluminum or iron and steel (Global Energy Infrastructure Emission Database, GID; Wang et al.2019) and surrounding coal-fired power plants (Global Energy Monitor, GEM; and the Global Power Plant Dataset, Byers et al.2018) that support the nearby industries.

Our goal is to calculate ERCO from every satellite sounding within an urban plume, which is a downwind area affected by urban emissions (Sect. 2.2). Sounding-dependent ERCO values are calculated as ratios of CO emissions over CO2 emissions (Eq. 3) that are estimated from satellite-derived fossil fuel (FF) enhancements and further refined with the “scaling factor” in Eq. (1). This scaling factor accounts for several mismatches between OCO-2/3 and TROPOMI (Sect. 2.1) and is obtained from an atmospheric transport model (Sect. 2.2.1). Since we do not differentiate emission signals due to biofuel and fossil fuel combustion, the term “FF enhancement” is used to refer to the column enhancement induced by any anthropogenic combustion processes in the target city. The determination of FF enhancements requires estimates of the background values (Sect. 2.2.2) and “second-order” correction terms for biogenic and pyrogenic sources (Sect. 2.2.3). Sounding-specific ERs and uncertainties (Sect. 2.2.4) are aggregated to yield an ER per overpass and per city. Lastly, we illustrate how the ERCO values associated with heavy industry in Los Angeles and Shanghai can be extracted with the assistance of WUDAPT (Sect. 2.3).

2.1 Satellite observations and data pre-processing

We evaluate all coincident OCO-3 SAM and TROPOMI overpass observations, but only select those with relatively small differences in overpass times. Considering the limited number of coincidences between sensors, two non-SAM overpasses from OCO-3 and one OCO-2 overpass are added to the analysis. As a result, six OCO–TROPOMI coincidences with high data quality from October 2019 to June 2021 are integrated into the final result for every city. Two of the total of 24 overpasses fall within the Northern Hemisphere summer months (both in June).

2.1.1 OCO-2/3 XCO2

The column-averaged dry-air mole fraction of CO2 (XCO2) is retrieved from the reflected sunlight over two CO2 bands centered on 1.6 and 2.0 µm and the oxygen A band for obtaining the surface pressure (Eldering et al.2019; Taylor et al.2020). In addition to the standard nadir, glint, and target modes, OCO-3 collects several adjacent swaths of XCO2 observations over a spatial area of approximately 80 km by 80 km during its SAM mode, e.g., four individual swaths in an overpass over LA on 24 February 2020 (Fig. 2a). Similar to OCO-2, each satellite swath comprises eight spatial footprints/soundings, and each sounding has an area of  1.6 × 2.2 km2 at nadir (Fig. 2a). Our analysis only uses screened OCO-3 B10r/B10p4r data (Eldering2021) with an XCO2 quality flag of zero (QF = 0). It is worth highlighting that the B10r/B10p4r product is superior to the Early version of OCO-3 (Taylor et al.2020); it has improved geo-location, advanced radiometric calibration, improved quality filters, and customized post-processing bias correction. As OCO-3 is mounted on the International Space Station, which is in a precessing orbit, its overpass time varies (for example, from 07:00 to 15:00 LT (local time)) for the overpasses we examine, unlike OCO-2.


The TROPOMI column density of CO molecules [mole cm−2] is retrieved via the measured radiation from shortwave infrared wavebands centered at  2.3 µm (Veefkind et al.2012). We select soundings with a quality assurance of ≥0.5 as recommended by the TROPOMI README document (Landgraf et al.2020) and convert the vertical column density to the total column-averaged dry-air mole fraction of CO [XCO in ppb] by calculating the dry-air column density [mole cm−2] using the retrieved surface pressure and total column water vapor. TROPOMI CO is retrieved from a larger pixel area of  7 × 7 km2 at nadir, which reduces to 5.5 × 7 km2 after 6 October 2019 (Fig. 2c). The overpass time of TROPOMI is  13:30 LT for an equatorial overpass in nadir measurements, with a time span of 1–2 h for soundings on the edge of the wide swath (i.e.,  2600 km).

Figure 2Spatial maps of FF XCO2 enhancements with corrections for non-FF sources at the native OCO-3 scale (a, ppm) and aggregated at the TROPOMI scale (b, ppm) along with FF XCO enhancements (c, ppb) over Los Angeles on 24 February 2020. Power stations with different primary fuel types are displayed as different white symbols based on the Global Power Plant Dataset (Byers et al.2018). The overall X-STILT column footprint [ppm (µmol m−2 s−1)−1] from all soundings is drawn in light gray (see explanations for the column footprint in Sect. 2.2.1). The underlying hybrid maps were created using the ggmap library in R with the hybrid view from Google Maps over LA (copyright: map data © 2021, imagery © 2021 TerraMetrics).

2.1.3 Differences between the two sensors/species

Four mismatches between OCO-3 XCO2 and TROPOMI XCO that pose challenges when extracting FF enhancements and ERCO from atmospheric observations are accounted for in this analysis:

  1. Satellite pixel area. XCO2 enhancements from multiple OCO-2/3 soundings falling within a given TROPOMI polygon are grouped and averaged (Fig. 2a vs. b). For simplicity, the centered lat/long coordinate of an OCO pixel is used to determine its corresponding TROPOMI polygon. The retrieval uncertainty tied to each OCO sounding is also aggregated according to the TROPOMI sampling, contributing to the total observational uncertainty (Sect. 2.2.4).

  2. Averaging kernel profile. Within the planetary boundary layer, where most emissions occur, TROPOMI XCO retrieval is affected by cloud height/fractions, which yields a lower-than-unity AK (Supplement Fig. S1). The OCO-2/3 XCO2 retrieved under cloudy conditions is typically omitted from Lite files and when QF = 0 is applied; thus, its AK normally approaches 1 near the surface for cloud-free scenes. The mismatch in AK between sensors must be accounted for as it can affect the interpretation of ERs. In this work, we account for AKs within an atmospheric transport model (Sect. 2.2.1).

  3. Overpass times, meteorological conditions, and emission variations. As a result of the overpass time difference between sensors, variations in meteorological conditions (e.g., wind direction and speed) can lead to changes in the urban plume shapes detected by the two sensors as they pass by. We deal with changes in wind speed and wind direction separately. The former is resolved by using the “scaling factor” inferred from an atmospheric transport model and the latter undergoes manual evaluations (Sect. 3.1). Also, the CO and CO2 emissions themselves can vary over the course of a day, driven by, e.g., the road transportation and residential sectors. Given the overpass time difference between sensors, it is likely that such a mismatch in the timing of CO versus CO2 emissions may affect the observed ERCO.

  4. Non-fossil-fuel sources/sinks. Not accounting for the influences from the biosphere and biomass burning may bias ERCO. Given our definition of the “local background”, the contrast in non-FF concentration anomalies between the urban and the background regions needs to be included (for more explanation, see Sect. 2.2.3).

2.2 Estimates of Egas, ERCO, and uncertainties

Previous studies (Mitchell et al.2018; Wu et al.2020; Lin et al.2021) proposed an approach for calculating an overall CO2 or CH4 flux using atmospheric measurements and an atmospheric transport model without relying on prior information from emission inventories. Here we briefly describe this approach to obtaining the overall emission of either CO2 (Eq. 1) or CO (Eq. 2) for a single sounding S, modified from Wu et al. (2020):


All the X terms in the numerator contribute to the estimate of the FF column enhancement (Xff, Fig. 2). Xff at a downwind satellite sounding S is the net result of FF sources over the source region (x,y). To describe the source region and attribute it to each satellite sounding, we adopt the column version of the Stochastic Time-Inverted Lagrangian Transport (X-STILT) model (Lin et al.2003; Fasoli et al.2018; Wu et al.2018). This model helps to provide a “scaling factor” 〈XFgas that accounts for the sounding-specific AK profile and meteorology (Sect.  2.2.1). Xbg denotes the local background values from satellite observations uncontaminated by the emission from the city, which are constant for a group of background observations. The background region is usually chosen to be a rural region outside the urban plume while considering the wind direction (Sect. 2.2.2). From a Lagrangian viewpoint, the air parcels arriving at an urban sounding might be traced back to different origins from the air parcels arriving at a rural sounding, meaning that observations at the two soundings may be influenced differently by the surrounding biosphere. Hence, two background correction δ terms are attached to account for the urban–background gradient in concentration anomalies due to net ecosystem exchange (NEE) and biomass burning (Sect. 2.2.3).

For a given sounding, the estimated flux Egas, with units of µmol m−2 s−1, represents the average emission over the corresponding source region of that sounding, which should not be confused with the direct emission at that sounding location. The ERCO for a given sounding S is then derived from Eqs. (1) and (2) as follows:

(3) ER CO , s = E CO , s E CO 2 , s = X ffCO , s X ffCO 2 , s XF CO 2 , s XF CO , s = X ffCO , s X ffCO 2 , s γ foot , s ,

where XffCO,sXffCO2,s is the observed enhancement ratio and γfoot,s measures how enhancement ratios with no consideration of AKs and meteorology differ from emission ratios. We simply use ppb-CO/ppm-CO2 as the units of ERCO (i.e., the same as mmol-CO/mol-CO2).

2.2.1 The X-STILT model

The X-STILT model is adopted in this study (1) to provide the scaling factor 〈XFgas that resolves differences in AKs and changes in wind speeds, (2) to identify an overpass-specific urban plume for determining background regions (Sect. 2.2.2), and (3) to estimate the sounding-specific biogenic and pyrogenic anomalies for background corrections (Sect. 2.2.3).

STILT releases an ensemble of air parcels from target observations (known as the “receptor”) and tracks the movement of those air parcels backward in time. The source region corresponding to each sounding is inferred from the “source–receptor relation” or the STILT “footprint” (Lin et al.2003; Fasoli et al.2018). The STILT footprint [ppm (µmol m−2 s−1)−1] describes the change in atmospheric concentration [ppm] at a downwind location due to possible upwind sources/sinks [µmol m−2 s−1]. The magnitude of the STILT footprint tends to be higher close to the target observation or under steadier wind conditions; thus, air parcels within the boundary layer can interact more closely with fluxes from the surface.

To accommodate the use of satellite-based column data, X-STILT incorporates retrieval-specific AK and pressure weighting profiles into the footprint calculation (Wu et al.2018) such that influences on air parcels originating from various altitudes of an atmospheric column are weighted by the sensor/species/sounding-specific vertical profile (Fig. S2). The “column footprint” (XFgas) measures the sensitivity of the total column concentration to upwind fluxes from the perspective of a specific satellite sensor. For instance, XFgas for TROPOMI XCO differs from XFgas for OCO-2/3 XCO2, even for concurrent observations. Since the airflow arriving at each satellite observation is unique, the magnitude and spatial distribution of XFgas vary across soundings (Fig. S3). By taking an average of these sounding-dependent column footprints, as shown in Fig. S3, we can identify the source region for all soundings in a SAM (light gray area in Fig. 2b, c). In this work, we only traced air parcels back for 12 h to calculate column footprints, which is sufficient to capture the near-field influence from the target city and better aligned with the local background region outside the city (Sect. 2.2.2).

In short, the spatial summation of column footprints 〈XFgas is regarded as a scaling factor to address the sounding-specific meteorological conditions and AK profile. The term γfoot derived from Eq. (3) reveals the difference between a simple enhancement ratio and a more robust, model-corrected emission ratio.

2.2.2 Background definition

Defining accurate background levels to extract urban FF enhancements has always been a challenge in top-down analyses, especially when dealing with column measurements with small signal-to-noise ratios. Wu et al. (2018) compared several approaches to determine a localized XCO2 background for extracting urban signals from OCO-2, including approaches that (1) solely use satellite observations with statistics (e.g., daily median); (2) solely use an atmospheric transport model (e.g., the “curtain method” based on global concentration fields); and (3) combine observations and transport information from models. Here, we expand the third approach to arrive at localized swath-dependent background values. The broader spatial coverage compared to the narrow swath of OCO-2 and multiple swaths stretching out of the city domain of OCO-3 SAMs help improve such background determination by introducing spatial variations in the background. Accurately describing latitudinal or spatial gradients in the background XCO2 has been emphasized recently (Ye et al.2020; Schuh et al.2021).

The process of background determination used in this work involves the first step of identifying the urban plume and differentiating soundings as being within or outside of the plume. To outline the urban plume shape at the overpass time, we utilize the forward mode of STILT with the inclusion of wind uncertainty in atmospheric dispersion. Specifically, 1000 air parcels are released continuously from a rectangle representing the city domain (dashed black box in Fig. 3) every 30 min starting 10 h ahead of the overpass time. All air parcels are allowed to travel forward in time for 12 h from their initial release times. A random wind component typifying model-data wind errors is added to the parcel dispersion (Lin and Gerbig2005). We subset the air parcels only during the overpass time and apply a two-dimensional kernel density estimate (KDE) based on the parcels' spatial distributions (blue to purple contours in Fig. 3). KDE is carried out using the kde2d function provided by the MASS library in R (Venables and Ripley2002). These normalized KDE contours indicate the likelihood and shape of an urban plume when the satellite scans through. The extent of the urban plume is finalized using a normalized KDE contour of 0.15 (black curve in Fig. 3), which is appropriate to include soundings with a possible influence from the target city and to exclude observations elevated by another city (e.g., the red polygons centered at  32 N and 120 E in Fig. 3c). This procedure is carried out separately for OCO-2/3 and TROPOMI to reveal the impact of changing meteorology on urban plumes at different overpass times (see Sect. 3.1). It is worth stressing that only enhancements within the urban plume are used for ERCO estimates.

Figure 3Demonstrations of background determination from OCO-2 XCO2 on 4 February 2020 (a) and OCO-3 SAM XCO2 [ppm] and TROPOMI XCO [ppb] on 20 February 2020 (b, c) over Shanghai. The model-based urban plume (solid black curve) is determined by the normalized 2-D kernel density of the air parcel distribution during the overpass time (blue-purple contours). The “background soundings” outside the urban plume are highlighted with black outlines, while other soundings are outlined in white. For example, OCO-2 observations to the north outside the plume ( 121.1 E, 31.9 N, a) and OCO-3 and TROPOMI soundings to the southeast outside the plume (b, c) are used to estimate background values and background uncertainties. The underlying hybrid maps were created using the ggmap library in R with the hybrid view from Google Maps over Shanghai (copyright: map data © 2021, imagery © 2021 TerraMetrics).

Next, the background value is calculated as the median value of observed Xgas per swath over the background region. For example, the background region is the area to the east outside the urban plume since southeasterly wind dominates (Fig. 3b, c). Background values vary with swaths if an OCO-3 SAM is examined. We choose the median instead of the mean to minimize the impact of any “outliers” that may be from a second FF source (other than our target cities) in the background region. Background uncertainty is estimated as a component of the total observed uncertainty (Sect. 2.2.4).

2.2.3 Background correction terms for non-FF sources/sinks

The swath-dependent local background approach described above explicitly assumes equal contributions from non-FF sources and sinks for soundings in the background versus soundings in the urban plume, which may not always be the case. We then correct for the spatial gradient in contributions from biogenic and pyrogenic fluxes.

As proposed in Wu et al. (2021), rather than absolute biogenic concentration anomalies, it is the contrast in these anomalies between the background versus the urban plume that is required, considering our localized background definition. Specifically, hourly X-STILT column footprints are convolved, respectively, with hourly mean NEE from a biospheric model representation and daily mean wildfire emissions from the Global Fire Assimilation System (GFAS, Kaiser et al.2012) to estimate the sounding-specific absolute column anomalies Xbio and Xbb. The Solar-Induced Fluorescence (SIF) for Modeling Urban biogenic Fluxes (SMUrF, Wu et al.2021) model estimates gross primary production (GPP) from a contiguous SIF product (CSIF, trained based on OCO-2 SIF, Zhang et al.2018) and respiration based on modeled SIF-based GPP and air and soil temperatures.

Next, the urban–background gradient in such anomalies is calculated as the difference between sounding-specific anomalies and the mean anomaly within the background region:

(4) δ X bioCO 2 ( s ) = X bioCO 2 ( s ) - X bioCO 2 ( s bg ) ,

where s or sbg represents all the soundings or select soundings in the background region, respectively. Let us imagine a summer day at noon in the Northern Hemisphere. The urban core is normally associated with a weaker biospheric uptake than the surrounding rural region. Biogenic signals XbioCO2(s) for soundings in the city are less negative than the mean biogenic signal over the rural background XbioCO2(sbg). Hence, the urban–background biogenic gradient δXbioCO2(s) is normally positive and subtracted from the total column (Eq. 1). The estimated Xbio values and their urban–background gradient δXbio are shown in Sect. 3.1.

Flux exchanges from the ocean and chemical transformations (e.g., the CO sink from the hydroxyl radical (OH) and the source from the oxidation of volatile organic compounds, VOCs) are not considered. The average lifetime of CO against OH ranges from a few weeks to several months depending on the season – much longer than the few-hours timescale we care about. Yet, CO can be generated from the oxidation of CH4 and non-methane VOCs at various rates, which is discussed in Sect. 4.3.

2.2.4 Uncertainty sources

The uncertainty related to emissions should contain uncertainties from (1) the atmospheric transport (i.e., column footprints), (2) observations, and (3) non-FF sources and sinks, according to Eqs. (1) or (2). We neglect uncertainties from column footprints, assuming that no transport bias exists during either the OCO or the TROPOMI overpass time. The urban–background gradient in non-FF fluxes remains very small compared to FF enhancements (Sect. 3.1).

We estimate the uncertainties of observed FF enhancements following Eq. (5). As previously described, observations from a few screened OCO soundings ( 5 to 28 OCO soundings, depending on the TROPOMI footprint size) are averaged to arrive at a mean XCO2 at the TROPOMI scale. Due to this averaging/binning process, the XCO2 uncertainty due to binning is considered using the standard deviation of XCO2 observations (σε,bin2 in Eq. 5) within a TROPOMI polygon.

(5) σ ε , obs 2 = σ ε , bin 2 + σ ε , bg 2 + σ ε , retrv 2 .

Here, σε,bin2 is not required for estimating the XCO uncertainty. The background uncertainty (σε,bg2) contains both the retrieval error and the variability of column observations (as standard deviations) within background regions.

The retrieval uncertainty (σε,retrv2) of XCO is available for each TROPOMI sounding, whereas that of XCO2 is reported for individual OCO-2/3 soundings (as read from Level 2 Lite files), which need to be aggregated at the TROPOMI scale. Due to possible correlations in retrieval errors between nearby OCO soundings, we estimate the error correlation length scale (Lx) using exponential variograms, as demonstrated in Fig. S4. Within a TROPOMI polygon that contains N (i.e., the number of) OCO soundings, an error variance–covariance matrix with dimensions of N×N is constructed with diagonal elements filled with OCO sounding-specific retrieval error variances. Then, Lx is used to form the normalized covariance matrix, i.e., exp(-D(Si,Sj)Lx), where D(Si,Sj) denotes the distance between each two OCO soundings (1i<jN). Lastly, the sum of all elements in the error covariance matrix (both variance and covariance elements) is divided by N2 to obtain one σε,retrv2 per TROPOMI grid. As a result, the overall uncertainty of FF enhancement per sounding is often dominated by the background error component.

2.3 Identifying the ERCO for heavy industry within a city

A key objective of this study is to explain the intra-city variability of ERCO by exploring sector-specific or sector-dominant combustion activities. While certain combustion processes and sectors tend to have higher ERs than others, sectorally dependent ERs are variable within and across cities. The ERs derived from atmospheric observations comprise a mixed effect of different activities in the city. Previous attempts include reducing the number of sectors and relying on prior sector-specific ERs via a (joint) Bayesian inversion (Brioude et al.2012; Nathan et al.2018).

Figure 4Maps of 120 m Local Climate Zone (LCZ) from WUDAPT (a and d) along with magnified images (b and e) and interpolated areal coverage of the heavy industry [%] at 1 km around Shanghai and Los Angeles (c, f). LCZ classifications centered on Wuxi and Shanghai are combined. Based on the Global Power Plant Dataset (Byers et al.2018), power stations are drawn as white symbols. The dashed light gray or white rectangles in the maps indicate the magnified region.

Here, we propose a novel approach to identifying ERs associated with heavy industry in a city. Instead of relying on prior emission inventories that can sometimes be erroneous regarding the magnitudes and the locations of sector-specific activities (see discussions in Sect. 4.4), we utilized an urban land cover classification dataset, WUDAPT, that provides Local Climate Zone (LCZ) classifications at a grid spacing of 120 m (Ching et al.2018). As shown in Fig. 4a, d, LCZ categories include street canyons (e.g., compact/open/lightweight, high/mid/low rise), building spacing (e.g., sparsely built, heavy industry), and tree spacing (e.g., dense/scattered trees, low plants, rocks, etc.). Each LCZ is unique in its thermal, radiative, and metabolic properties. For instance, the compact high-rise (LCZ 1) and heavy industry (LCZ 10) categories have the highest anthropogenic heat outputs of 50–300 and >300 W m−2, respectively (Stewart and Oke2012). Heavy industry is defined as low-rise and mid-rise industrial structures (towers, tanks, stacks) and mostly paved or hard-packed metal with steel and concrete construction materials and few or no trees in WUDAPT (Ching et al.2018), which differs from the industry-relevant sectors defined by the Intergovernmental Panel on Climate Change (e.g., as used in EDGAR). We clarify that we are not trying to tackle individual industrial processes, which is much more difficult. As of this analysis, LCZ maps are only available for a limited number of cities, including Shanghai and LA, but they have recently been generalized to the entire globe (Demuzere et al.2022a).

To relate ERCO to heavy industry, the percentage of heavy industry is first interpolated using 1 km grid spacing from WUDAPT LZC maps (%, Fig. 4c, f). The industrial coverage map is then convolved with the X-STILT column footprint (Fig. S3) to quantify the industrial influence on each TROPOMI polygon Pind(x,y), which is defined as the column footprint-normalized industry fraction (Fig. S5). For example, soundings in the city center farther away from the heavy industry in LA are related to smaller influences. Lastly, we sum Pind(x,y) across the space to arrive at Pind, which serves as a metric of how much the observation at a given sounding is affected by heavy industry. Specifically, soundings with Pind larger than the 75th or 90th percentile are marked as locations that are “impacted” or “strongly impacted” by heavy industry within the city. Sensitivity and significance analyses are conducted and presented in Sect. 3.2.2; these test if: (1) the results are subject to the percentile threshold when defining industry-dominated soundings; and (2) ERs over industry-dominated soundings are statistically significantly different from ERs for the remaining soundings.

Figure 5Examples of modeled urban plumes during OCO-3 (red curve) and TROPOMI (blue curve) overpass times (in UTC). The likelihood of these meteorology-only urban plumes (no emission involved) is quantified by the normalized KDE binned in 10 intervals of the modeled air parcel distribution (yellow-green-purple contours). Three types of overpasses are shown, as follows: (a, b) “good cases” with almost identical urban plumes at two times, e.g., Baotou on 8 October 2019 with a Δt of 8 min; (c, d) cases to be used with caution, where the urban plume shifts from one time to another and requires a simple plume rotation, e.g., Shanghai on 20 February 2020 with a Δt of 1.5 h; (e, f) outliers where two urban plumes change significantly, e.g., for Los Angeles on 3 March 2020 with a Δt of over 4 h. The underlying hybrid maps were created using the ggmap library in R with the hybrid view from Google Maps over LA (copyright: map data © 2021, imagery © 2021 TerraMetrics).

3 Results

ERCO values and uncertainties are reported at multiple spatial scales, from the spatially resolved sounding level (Eq. 3) to the overall overpass and city level. Again, only ERs at soundings within the urban plume are selected. Overpasses with too few valid soundings in a plume area are also removed from the results. Before presenting ERs at different spatial scales, we assess factors that may influence the derived ERCO.

3.1 Interference factors that modify ERCO

We examine impacts on ERCO from the following interference factors: (a) differences in AKs between OCO-2/3 XCO2 and TROPOMI XCO; (b) shifts in wind fields between two overpass times; (c) the urban–background contrast in biogenic and pyrogenic contributions; and (d) temporal variation in the emissions themselves. In summary, we find that differences in AKs and wind directions between sensors can significantly affect the spatially resolved ERCO. For the final 24 overpasses we selected, temporal variations in the emission pattern and urban–background gradients in biogenic/pyrogenic contributions play minor roles in overpass- or city-level ERs.

Recall that sounding-specific AKs and wind speeds were considered in the sounding-specific column footprint using X-STILT (Sect. 2.2), and γfoot = XFCO2XFCO measures the overall contributions from AKs and wind speeds to the spatially resolved ERCO (Supplement Fig. S6c). For instance, the mean γfoot spans from 1.20 to 1.57 over LA and from 1.02 to 1.38 over Shanghai across different overpasses (printed in Fig. 7a, c). γfoot is generally larger than 1 because AKs of TROPOMI XCO near the surface are smaller than surface AKs for OCO-2/3. Simply using enhancement ratios without accounting for mismatches in AKs and wind speeds between sensors will likely lead to an underestimation of emission ratios (Eq. 3). On average, the overpass-level ERCO can be  20 % higher than enhancement ratios across our 24 overpasses.

The second factor is the change in wind directions between two overpass times, which is evaluated using the same algorithm as the urban plume detection in Sect. 2.2.2. Again, colored contours and curves in Fig. 5 indicate neither the intensities of concentrations nor flux fields (as no prior emissions are used) but rather the likelihood of urban plumes, determined by atmospheric dispersion with random wind uncertainties. Matching between OCO-3 soundings and TROPOMI polygons as described earlier would be fine for concurrent observations (Fig. 5a, b), but this becomes problematic if Δt becomes large (“outliers” with significant changes in plumes in Fig. 5e, f). The cases between the good cases and the outliers are to be used with caution (Fig. 5c, d). By comparing the overlap of plumes at the two times, we shifted OCO-3 soundings to better align with TROPOMI polygons. For example, on 20 February 2020, because the modeled plume at the OCO-3 overpass time (06:06 UTC) appears northward compared to the plume at the TROPOMI overpass time (04:44 UTC), OCO-3 soundings were shifted southward by zero to two grids, depending on their longitudinal coordinates (Fig. S7). In other words, by shifting the FF XCO2 enhancements, we better align the urban plume at the OCO-3 time with the plume at the TROPOMI time. Every OCO–TROPOMI coincidence is manually examined and assigned to one of the three categories, which are further summarized in Sect. 4.1. Outliers are removed from this analysis, since no simple wind or plume rotation would improve their ERCO estimates.

Besides changes in wind directions, the CO and CO2 emissions themselves can vary across daytime hours, likely driven by the road transportation and residential sectors. As a result, variations in the derived ERCO across multiple overpasses may reflect not only the variation in combustion efficiencies but also the mismatch in the emission timing. LA may be a city with more distinct daytime changes in emissions compared to industry-centered cities. Fortunately, based on a supplementary sensitivity analysis using measurements from the Total Carbon Column Observing Network in Pasadena (TCCON, Wennberg et al.2017), observed ERCO values appear to be less variable when limiting satellite overpasses to those with a smaller time difference (Fig. S8). Future geostationary satellite monitoring of NOx (e.g., TEMPO, Chance et al.2022) may provide better guidance regarding the hourly pattern in urban emissions, especially from the traffic sector, which show more daytime fluctuations, as discovered using surface monitoring networks (e.g., over Chicago; de Foy2018).

The last factor is the urban–background contrast in contributions from non-FF sources and sinks. The biogenic XCO2 anomaly modeled using SMUrF and X-STILT ranges from −0.7 to 0.3 ppm per OCO-3 sounding, depending on the hour of the day (i.e., the solar zenith angle), season, and wind direction (Fig. S9). As explained in Sect. 2.2.2, urban–background gradients in these biogenic anomalies (i.e., δXbio) were used to correct the constant localized background Xbg (Eq. 1). Take the two overpasses with the largest urban–background contrast as examples: as biospheric uptake is normally weaker in urban areas than in surrounding rural areas (often used as background regions), the urban–rural gradient for locations in the plume region becomes more positive (Fig. S10b). Nonetheless, even for the one summertime SAM over Zibo on 21 June 2020, the sounding-level δXbio ranges from 0 to 0.4 ppm, which remains small compared to the sounding-level FFCO2 enhancements of 2 to 7 ppm (Fig. S11a, b). For most other overpasses, δXbio aggregated according to the TROPOMI sampling stays low, with an absolute value of <0.3 ppm (as printed in each panel of Figs. 6 and 7). Even with a bias in the resultant δXbio resulting from an incorrect prior NEE, the effect on the derived FF enhancements and ERCO would be small.

Although LA is surrounded by occasional intense wildfire activities, the column anomalies due to biomass burning that are suggested by the coupling of GFAS and X-STILT are minimal for the dates we examined. Yet, since wildfire-related ERCO values are usually higher than the FF-related ERCO (Fig. 1), properly accounting for pyrogenic contributions and gradients between urban and surrounding rural areas is important for cities in mountainous and forested areas during fire seasons. For instance, Crounse et al. (2009) leveraged aircraft measurements of HCN and C2H2 over Mexico City as indicators to disentangle the CO signals due to biomass burning and urban emissions, respectively.

3.2 Intra-city ERCO variations and signals from heavy industry

Observed enhancements are the net consequence of associated sources/sinks from source regions. That is, a high atmospheric content of CO2 or CO at the sounding location does not necessarily indicate a high emission rate at this location (Kiel et al.2021). Our derived emissions and ERs, although reported for each sounding, indicate the overall emission and combustion efficiency over its source region.

In the following subsections, we present ERs for each sounding and the aggregate for each overpass and city. Since the aggregation of sounding-level ERs to a single value per overpass or city is sensitive to the method/statistic adopted, we bootstrapped ECO and ECO2 based on their sounding-specific values and uncertainties to generate a linear regression fit per bootstrap loop (light gray lines in Fig. 6). Specifically, 1000 random sets of ECO and ECO2 were generated according to assumed normal distributions, where sounding-level emission estimates provided mean statistics with observational uncertainties as standard deviations. We used the standardized major axis (SMA) solution for linear regression to minimize deviations of data points from the regression line for both axes. Eventually, we obtained 1000 bootstrapped slopes and selected slopes with positive values, which yielded the overpass-level ERCO and uncertainty (e.g., dashed colored lines and text in Fig. 6). Also, sounding-level ERCO values from all overpasses are presented in histograms and generally follow a log-normal distribution (Fig. 7b, d).

Table 1A summary of the total power generation capacity (from the Global Power Plant Dataset of the World Resources Institute, Byers et al.2018) and information on heavy industry, including the annual crude iron capacity (GID, Wang et al.2019). Power plants are selected from a 0.5× 0.5 region around each city; percentages indicate the partitioning by fuel type.

Download Print Version | Download XLSX

Figure 6Scatter plots of CO and CO2 fluxes [µmol m−2 s−1] and their uncertainties (error bars as dotted lines) for Baotou (a) and Zibo (b). Linear regressions are applied to the data from each overpass (dashed colored lines) and from all overpasses (dashed black lines). Only bootstrapped regression lines with positive slopes were chosen from the Monte Carlo experiment (dashed gray lines;  98.4 % and 93.3 % of the total 6000 bootstrapped lines for Baotou and Zibo, respectively). The TROPOMI overpass time (in UTC), the total TROPOMI sounding number, the discrepancy in overpass times (Δt, min), the impact of the AK and wind conditions between sensors (γfoot, unitless), and the range of the urban–rural biogenic gradient (δXbio, ppm) are labeled in each panel. XCO2 values on 6 December and 28 December 2019 over Zibo came from the non-SAM nadir OCO-3 observations.


3.2.1 Baotou and Zibo

Combustion efficiencies are generally poor for the two industry- and energy-oriented cities. The overpass-specific ERs span from 9.3 ± 1.2 to 24.6 ± 3.8 mmol mol−1, with an integrated city-level estimate of 17.3 ± 0.5 mmol mol−1 for Baotou (Fig. 6a). According to GID, the Baotou Iron and Steel Group is located within the city and contributes to an annual capacity of crude iron of 12 619 kt yr−1, with estimated CO2 emissions of 20 462 kt yr−1 (Table 1). The slightly lower ERCO and FF enhancements in February 2021 coincide with the timing of the Spring Festival in 2021 ( 12 February). Standard deviations of bootstrapped slopes are higher for overpasses with fewer high-quality satellite soundings, e.g., 3.8 mmol mol−1 for overpasses with seven available TROPOMI polygons in the urban plume on 31 May 2020. Utilizing the bootstrap method helps to account for the impact of the sounding number on the overall city-level ER estimate.

Figure 7(a, c) Same as Fig. 6 but for LA and Shanghai. Only bootstrapped regression lines with positive slopes are presented as light gray lines ( 94.3 % and 88.4 % of the 1000 bootstrapped lines per overpass for LA and Shanghai, respectively). (b, d) Histogram of ERCO for all soundings (black bars). Soundings impacted or strongly impacted by heavy industry are defined as having Pind larger than their 75th or 90th percentile (blue or red bars); the corresponding median ERs are also shown (vertical dashed lines). The industrial impact is quantified using column footprints from X-STILT (to account for the atmospheric transport) together with localized information from the urban land cover data WUDAPT.


Zibo, along with the nearby county-level city of Zouping, accounted for over one-eighth of the total coal consumption of Shandong Province in 2017. The coal-fired power plants in the area contribute to a total power generation capacity of 9720 MW (Table 1), which is likely to support local metal industries, especially the producers of electrolytic aluminum (they are the world's top producers). The maximum XCO2 enhancement per OCO-3 sounding can even reach up to 10 ppm for a few overpasses (not shown). Interestingly, the ERCO for Zibo first declined from 10.1 ± 1.1 to 6.1 ± 0.6 mmol mol−1 during February 2020 and then gradually increased back to 18.2 ± 1.1 mmol mol−1 by June 2020 (Fig. 6b). Such temporal variations in ERCO agree nicely with the timing of the initial phase of the COVID-19 lockdown in China (i.e., February to May 2020) (e.g., Laughner et al.2021). We suspect that changes in ERCO could be driven by the partial shutdown and reopening of the multiple coal-fired power plants and metal industries in the area.

3.2.2 Los Angeles and Shanghai

Although OCO-3 has sampled the Los Angeles Basin dozens of times to date, many overpasses did not pass the quality check (i.e., QF) and were removed from the final result due to their noticeable shifts in urban plumes between two overpass times (e.g., 3 March, 15 April, and 5 May 2020 for LA; discussed in Sect. 4.1). The overpass-level ER ranges from 7.4 ± 0.8 to 11.7 ± 1.5 mmol-CO/mol-CO2, with a city-level value of either 9.6 ± 0.5 mmol mol−1 (obtained using the regression approach; Fig. 7a) or 9.7 mmol mol−1 (obtained using the histogram approach; Fig. 7b). Our space-based ERCO estimates over LA fall within the range of 7.1 to 12.4 mmol mol−1 reported from prior studies (Wennberg et al.2012; Brioude et al.2013; Hedelius et al.2016; Silva and Arellano2017). Small discrepancies in ERCO between studies may be attributed to discrepancies in the times of interest, sampling strategies, and techniques used for ERCO calculations (e.g., background definition).

In contrast to LA, where urban plumes are usually well constrained by the basin, wind speeds and directions vary across different overpasses over Shanghai – i.e., there is a southeasterly wind on 4 February and 20 February 2020, a southwesterly wind on 24 February 2020 and 19 February 2021, and a northerly wind on 23 April and 30 December 2020. Such changes in the wind regime between overpasses over Shanghai suggest that soundings from an individual overpass may reflect emission patterns over different source regions, which emphasizes the importance of integrating atmospheric transport when interpreting temporal variations in observation-based ERs. In other words, one cannot simply use all the soundings over a city to calculate ERs; it is necessary to select those soundings that are affected by emissions from that city. The overpass-specific ER ranges from 4.2 ± 1.2 to 17.1 ± 6.2 mmol mol−1, with a city-level average of 10.2 ± 0.4 mmol mol−1 based on the linear regression approach (Fig. 7c) or 12.9 mmol mol−1 using the histogram approach (Fig. 7d).

Now we focus on the distribution of sounding-level ERs for these two megacities (Fig. 7b, d) to see if ERs associated with a part of a city (i.e., the heavy industry region) can be revealed. As described earlier, to locate the soundings affected or strongly affected by the heavy industry in a city while accounting for the overpass-specific meteorology, we coupled the LZC-based industrial coverage (Fig. 4c, f) with X-STILT column footprints and quantified the industrial influence, Pind, at each sounding location.

Industrial regions within the LA Basin are concentrated to the south, near the Port of LA; and to the west of downtown, near Los Angeles Airport and the Chevron Refinery in El Segundo (Fig. 4e). The distribution of ERs for industry-dominated soundings tends to shift slightly towards the lower end (blue or red bars in Fig. 7b) compared to the distribution for all soundings (gray bars in Fig. 7b). For example, ERs of >15 ppb ppm−1 are less frequently found for industry-dominated soundings (red bars in Fig. 7b). The industry-oriented soundings generally have slightly lower CO but higher CO2 enhancements (Fig. 2b, c) compared to other soundings within the basin, resulting in slightly lower ERCO. No iron and steel facilities or coal-fired power plants are found over the heavy industry area in LA according to GID and GEM. We hypothesize that the slight shift of ERCO towards the lower end may be explained by the heavy-duty diesel engines and natural gas power plants in the Port of LA versus the predominately gasoline vehicles across the city, because the ERCO for heavy-duty diesel vehicles and non-coal-fired power plants is generally lower than that for light-duty gasoline vehicles. For example, a field campaign in 2007 in Beijing that split observations into daytime versus nighttime observations suggested that the ER linked to nighttime diesel transportation is much lower than that for the gasoline sub-sector (Westerdahl et al.2009, Fig. 1a). Similar to LA, a higher fuel efficiency was found over the ship channel of Houston (ER of  4 ppb ppm−1) compared to downtown Houston (ER of  10 ppb ppm−1) (Brioude et al.2012, Fig. 1b). Unfortunately, only two good SAMs near Houston from late 2019 to June 2021 are available, but future work can further validate the urban-industry contrast in ERs from space.

In Shanghai, the heavy industry is concentrated to the north of the city center (Fig. 4a). Interestingly, in contrast to LA, ERs affected by heavy industry are skewed towards the higher end, with medians of 16.8 or 18.8 ppb ppm−1 (blue or red bars in Fig. 7d) compared to the city-level median of 12.9 ppb ppm−1 (black bars in Fig. 7d). CO and CO2 enhancements and ERCO are all higher for industry-oriented soundings than for all soundings combined. Such spatial divergence in enhancements and ERs between heavy industry and the entire city may be attributed to substantial CO emissions from iron and steel production. Schneising et al. (2019) also found that many hotspots with high TROPOMI CO enhancements in China and India are tied to the iron and steel industries. During their production processes, iron ores are reduced to crude iron and steel, with CO involved. According to a plant-level estimate in 2019 from GID, Baoshan Iron & Steel Co., Ltd., located to the north of downtown Shanghai, has an annual crude steel capacity of 25 099 kt yr−1 (Table 1) and a total CO2 emission of 32 148 kt yr−1 for all coke, sinter, iron, and crude steel combined.

To validate the robustness of such ER shifts related to heavy industry, we tested the use of different percentile thresholds other than the 75th and 95th percentiles to determine industry-dominated soundings (Sect. 2.3). The above statements on industry-impacted ERCO hold if using alternative thresholds, including the 50th, 60th, and 80th percentiles. An additional Welch two-sample t-test confirms that ERs from industry-dominated soundings significantly differ from the remaining soundings that are less affected by heavy industry. When the adopted percentile threshold increases from the 50th to the 95th, the divergence in ERs between industrial and non-industrial soundings becomes more apparent, and the p value for the statistical significance of this difference becomes smaller (p values are <0.05 for all thresholds). In addition, the average number of OCO-3 soundings in a TROPOMI polygon is roughly the same for industry-affected soundings versus the rest (e.g., 11.8 vs. 10.3 for LA and 7.3 vs. 8.7 for Shanghai).

We acknowledge that although many iron/steel plants may aim at combusting as much CO as possible before releasing CO into the atmosphere, the indispensable role that CO plays in the iron/steel industry makes it unique when assessing its ERCO and combustion efficiency for various industrial processes. Furthermore, it is difficult to separate combustion signals of individual sectors from observations without prior assumptions of sector-specific contributions, since the atmospheric concentration at a given location arises from various underlying combustion processes spread over the source region. Even using additional co-emitted species, it would be risky to assume that a co-emitted species (e.g., CO or NOx) comes solely from one individual emission sector.

4 Discussion

This study is one of the first to analyze intra-city variations of emission ratios between CO and CO2 using two asynchronous satellite sensors. We describe complications induced by discrepancies between satellite sensors and retrievals and demonstrate methods to mitigate these complications by accounting for sounding-specific averaging kernels, atmospheric transport, and urban–background contrast in the contribution from non-FF sources/sinks using an atmospheric transport model.

4.1 Influences from non-FF components and atmospheric transport

Pyrogenic anomalies are minimal for the overpasses we examined but should be considered for certain cities (e.g., during dry seasons over Mexico City, Crounse et al.2009), considering the high ERCO from forest wildfires of 35 to 80 ppb ppm−1 (Fig. 1a). Most overpasses we analyzed fall within the dormant seasons. For the three overpasses during the growing season, biogenic anomalies modeled using the SMUrF model for a given OCO-2/3 sounding may reach up to 0.5 ppm (Fig. S9). Even though modeled NEE and resultant biogenic contributions/gradient can be uncertain, we stress again that it is the urban–background biogenic contrast (δXbio in Eq. 1) that is important for estimating FFCO2 enhancements, given our setup for a local background value. Satellite missions such as TROPOMI and the upcoming Geostationary Carbon Cycle Observatory (GeoCarb) will provide solar-induced fluorescence (SIF), which may help improve spatially explicit SIF-based GPP and NEE estimates (Turner et al.2020; Wu et al.2021), specifically by reducing the dependence on other remote sensing products and the assumption of model parameters for each plant functional type.

The biggest challenge affecting the robust estimation of spatially resolved ERCO is the shift in wind directions between two overpass times. Substantial changes in wind directions and urban plumes (e.g., the “outliers” in Fig. 5e, f) were mostly found for overpasses with an absolute time difference |Δt| of >2 h (implied by the bars labeled with an asterisk in Fig. 8). If TROPOMI pixel sizes are relatively large (i.e., non-nadir observations) or the wind is steadier, this |Δt| constraint may be relaxed as long as emissions for a specific city are less driven by sectors with noticeable diurnal cycles (e.g., road transportation). For instance, on 31 May 2020, TROPOMI polygon sizes for the industry-dominated city Baotou are sufficiently large compared to the shift in urban plumes, despite its |Δt| of 3 h (Fig. 8). In addition, we manually re-positioned the OCO-3 soundings to TROPOMI polygons for a few cases (the bars with nonzero numbers on top in Fig. 8) using the simple wind/plume shift demonstrated in Sect. 3.1. Fortunately, future geostationary satellites will be capable of mapping XCO and XCO2 simultaneously at a higher temporal frequency, which will eliminate this issue.

Figure 8A summary of the wind directional shift between OCO-2/3 and TROPOMI overpass times. The y axis denotes the absolute time difference (|Δt| in hours) per overpass. The color of the bar represents the instantaneous modeled surface wind speed [m s−1] from OCO-2/3 Lite files. Bars labeled with an asterisk indicate that urban plumes between two overpass times differ so much that they cannot be brought into agreement with a simple plume rotation. The number on each bar (e.g., 0–1) denotes the number of TROPOMI polygons that had to be shifted to align the urban plumes at the two times. For example, 0–1 means that TROPOMI polygons over certain locations are shifted by one grid. Bars labeled with a zero on top indicate that a manual plume shift was not required. Bars with white outlines indicate that the sampled TROPOMI soundings on that date were non-nadir with larger pixel sizes.


4.2 ERs for an individual sounding, overpass, city, and heavy industry region within a city

Contrary to previous work relying on inventory-based sector-based ERs, we attribute the intra-urban gradient to heavy industry using an urban land cover classification dataset. Such high-resolution localized maps help identify the observations strongly influenced by heavy industry. Based on a limited sample size, the heavy industry within the Greater Shanghai area is tied to an ERCO that is higher than the city average, reflecting poorer combustion efficiency (Fig. 7d). Industry- and energy-centered cities such as Baotou and Zibo are less efficient in their combustion activities. In particular, the industry-dominated ER over Shanghai (18.8 mmol mol−1, as indicated by the dashed red line in Fig. 7d) aligns better with the overall city-scale ER over Baotou of 17.3 mmol mol−1 (Fig. 6). The previously reported urban-integrated ERCO values are mostly constrained within the range of 5 to 20 ppb ppm−1, with a few exceptions in East Asia before 2010 that are over 30 ppb ppm−1 (Fig. 1). Our city-level estimates from space agree well with the range of previously reported values.

A city-scale ERCO derived from spatially explicit ERs can be influenced by (1) the adopted statistic, (2) overpass dates and overpass-specific wind conditions, and (3) estimated uncertainties in ERCO. For example, the overall ER derived from all soundings within the urban plume differs from that derived from a selection of soundings. Even though we started with all quantified OCO-2/3 observations in a SAM, only those located within the urban plumes (black curve in Fig. 3) can be used to estimate ERs, as this allows ERs from overpasses with different meteorological conditions to be compared in an unbiased way. The mean or median value of the sounding-level ERCO (e.g., 13.4 or 9.6 ppb ppm−1 for LA in Fig. S12) differs slightly from the city average obtained using the regression slope method when observational uncertainties are taken into account (e.g., 9.6 ppb ppm−1 for LA in Fig. 7a). Apart from these bulk quantities, the distribution of ERCO in the linear space is negatively skewed and roughly follows the log-normal distribution (Fig. S12), where a few observations with higher ERCO values are influenced by point sources with poorer combustion efficiency. More observations with finer satellite pixels across the city would improve the robustness of both the spatial distribution and bulk estimates of ERs.

4.3 Limitation

The main limitation of this work is the relatively low sample size, which is largely constrained by the requirement for small differences in overpass times. When more satellite data or upcoming data from geostationary satellites become accessible, intra-city ERs can be used to more robustly assess the temporal variation in sector-oriented combustion efficiency, including across seasons or times (e.g., business-as-usual scenarios versus pandemic-disturbed time frames). Beyond the sheer number of soundings, uncertainty arises when aggregating CO2 enhancements from the finer-resolution OCO-3 grid to the TROPOMI sampling. The centered lat/long coordinates of OCO-2/3 soundings are chosen to determine the corresponding TROPOMI polygon, while very few OCO soundings may be located right on the boundaries of TROPOMI polygons. Nevertheless, we find no significant bias associated with the number of OCO soundings per TROPOMI polygon in the heavy industry analysis.

Another factor that we did not explicitly account for is the secondary CO production from both anthropogenic and biogenic VOCs (AVOCs, BVOCs). Under a cascade of reactions in favorable conditions, VOCs emitted from the upwind source location are oxidized to CO at various rates, which may result in higher CO at the downwind sounding location and divergence between enhancement ratios and emission ratios. As BVOCs are usually associated with shorter lifetimes compared to many AVOCs (e.g., Surl et al.2018), we discuss BVOCs and AVOCs separately. BVOCs can contribute significantly to the total CO source at the regional scale, especially during growing seasons (e.g., Miller et al.2008; Hudman et al.2008; Gonzalez et al.2021). However, since BVOCs, like biogenic CO2, come mainly from rural areas outside the city, subtracting the localized CO background using CO observations outside the urban plume minimizes the impact of BVOCs on the derivation of CO enhancements. The lifetimes of AVOCs are long enough, except for a few species, including alkenes (Surl et al.2018). Without a good observational constraint on the VOC composition and group-specific emissions for different cities around the globe over the years, it would be challenging to accurately quantify the impact on atmospheric XCO and ERCO of AVOCs emitted from urban areas or specifically from industrial areas. More future efforts regarding urban VOCs may include (1) identifying good proxies that can be measured from space and can capture the bulk AVOC characteristics (e.g., formaldehyde, Zhu et al.2014) and (2) interpreting such observations, e.g., by utilizing chemical transport models for source attribution (Gonzalez et al.2021). Note that the noise/uncertainty in current daily TROPOMI formaldehyde observations may be too large for daily resolved analyses.

4.4 Implications for inventory evaluation

This work provides insights into estimating emission ratios from future satellite sensors, as ERs help pinpoint hotspots with poor combustion efficiency, which inform sub-city emission/pollution control efforts.

Satellite-based ER estimates help in the evaluation of sector-specific emission factors and source locations adopted in bottom-up emission inventories (e.g., Silva and Arellano2017). Substantial contrast in both the magnitudes and spatial distribution of enhancement ratios can be found between the observations and forward simulations (using X-STILT column footprints and sectoral emissions from EDGARv5). Taking Shanghai as an example, simulated enhancement ratios using prior emissions appear to be much higher (>50 ppb ppm−1, Fig. S13b) than observed ratios (mostly <30 ppb ppm−1, Fig. S13a). Regarding the spatial distribution, simulated enhancement ratios using total FF emissions mimic simulated enhancement ratios using only industry-related emissions (Fig. S13b versus c), and both simulated ratios differ substantially from observed enhancement ratios. Such model–data mismatches may result from inaccurate activity data and emission factors of EDGAR as well as atmospheric transport uncertainties. This preliminary analysis illustrates that satellite observations of trace gases could be used to evaluate emission factors adopted in bottom-up emission inventories. More sophisticated approaches, such as flux inversions (Hedelius et al.2018; Brioude et al.2011, 2012; Palmer et al.2006), may better constrain sector-specific CO and CO2 emissions from inventories.

Spatial proxies, including nightlight data from the Black Marble (, last access: 1 July 2022) and sophisticated urban land cover datasets, can support not only the development of emission inventories but also sector-orientated evaluations with atmospheric observations of CO2 and co-emitted pollutants. This work demonstrates the benefit of using high-resolution urban land cover classifications to provide independent information about locations of various anthropogenic activities, building structures, and vegetation coverage.

5 Conclusions

We investigated fossil fuel combustion efficiency by quantifying the emission ratios of CO and CO2 across Los Angeles, Shanghai, Baotou, and Zibo (Zouping) using nearly coincident observations of TROPOMI XCO and OCO-2/3 XCO2. Multiple swaths of observations collected by OCO-3 SAMs cover a much broader area relative to the OCO-2 swaths, facilitating the determination of background values and the separation of emission signals from different parts of a city. We incorporated spatial gradients in background values by calculating the background per swath and correcting for the urban–background gradient due to non-anthropogenic sources and sinks. Sensor-specific averaging kernel profiles and meteorological conditions were accounted for using an atmospheric transport model (X-STILT). The ratio between XCO and XCO2 enhancements without considering such sensor-specific factors is normally lower than the emission ratio. Cases with severe asynchronicity, specifically those with overpass time differences of over 3 h, correspond to significant changes in urban plumes. Properly accounting for the overpass-specific meteorological conditions or source–receptor relationship and identifying only the soundings influenced by urban emissions is critical when estimating ERs for cities, and is realized using an atmospheric transport model. Our model approach is then used to identify soundings strongly affected by heavy industry. As a result, the overall city-level ERCO for Shanghai (10.2 ± 0.4 mmol mol−1) is slightly larger than that for Los Angeles (9.6 ± 0.5 mmol mol−1). Industry-related ERCO values for Shanghai are much larger than its city-level average, whereas industry-related ERCO values for LA are slightly lower than its city-level average. ERs tied to heavy-industry regions in Shanghai (18.8 mmol mol−1) are approximately equal to the city-level ERCO for the industry-orientated city of Baotou (17.3 ± 0.5 mmol mol−1). High ERs highlight the poor combustion efficiency tied to certain industrial activities, e.g., metal production (Table 1).

Future satellites (e.g., GeoCarb, TEMPO, CO2M) will provide better spatial and temporal coverage of XCO2 and relevant co-located tracer observations, making it possible to monitor and verify temporal trends and variations in the combustion efficiency over hotspots within an urban area, which will provide significant guidance for urban planning and emission control.

Code and data availability

OCO-3 L2 B10r XCO2 data and TROPOMI XCO data were accessed at (Eldering2021) and (ESA2018), respectively. X-STILT code has been modified to work with TROPOMI data archived on GitHub branch at (last access: 1 July 2022) (DOI:, Wu et al.2019). We kindly ask users to follow the code policy for utilizing and acknowledging the X-STILT code for interpreting TROPOMI column data. Hourly NEE fluxes from SMUrF are archived in the Oak Ridge National Lab DAAC at (Wu and Lin2021). The urban land cover classification from WUDAPT can be downloaded from (DOI:, Stewart and Oke2012; Demuzere et al.2022b).


The supplement related to this article is available online at:

Author contributions

DW designed and carried out this analysis. JL, POW, and PIP supervised this study. RRN, MK, and AE provided the bias-corrected B10 data for the OCO-3 SAMs used in this work. All authors participated in the interpretation of the results and in the writing and editing of the paper.

Competing interests

The contact author has declared that none of the authors has any competing interests.


Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


The computations presented here were conducted in the Resnick High-Performance Computing Center, a facility supported by the Resnick Sustainability Institute at the California Institute of Technology. The first author appreciates the discussion with Joshua Laughner, Eric Kort, Tomohiro Oda, and John Lin. We thank Julia Marshall and the second anonymous referee for their careful reading of our submitted manuscript and for their constructive suggestions that have helped improve our study.

Financial support

The production of the OCO-3 science data products used in this paper was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (prime contract number 80NM0018D0004). The research effort was funded by the Jet Propulsion Laboratory Research and Technology Development project R.21.023.106. The analysis was supported by the W. M. Keck Institute for Space Studies and by the National Aeronautics and Space Administration (grant no. 80NSSC21k1064).

Review statement

This paper was edited by Jason West and reviewed by Stijn Dellaert and Julia Marshall.


Akagi, S. K., Yokelson, R. J., Wiedinmyer, C., Alvarado, M. J., Reid, J. S., Karl, T., Crounse, J. D., and Wennberg, P. O.: Emission factors for open and domestic biomass burning for use in atmospheric models, Atmos. Chem. Phys., 11, 4039–4072,, 2011. a, b

Ammoura, L., Xueref-Remy, I., Gros, V., Baudic, A., Bonsang, B., Petit, J.-E., Perrussel, O., Bonnaire, N., Sciare, J., and Chevallier, F.: Atmospheric measurements of ratios between CO2 and co-emitted species from traffic: a tunnel study in the Paris megacity, Atmos. Chem. Phys., 14, 12871–12882,, 2014. a, b

Ammoura, L., Xueref-Remy, I., Vogel, F., Gros, V., Baudic, A., Bonsang, B., Delmotte, M., Té, Y., and Chevallier, F.: Exploiting stagnant conditions to derive robust emission ratio estimates for CO2, CO and volatile organic compounds in Paris, Atmos. Chem. Phys., 16, 15653–15664,, 2016. a

Bares, R., Lin, J. C., Hoch, S. W., Baasandorj, M., Mendoza, D. L., Fasoli, B., Mitchell, L., and Stephens, B. B.: The wintertime covariation of CO2 and criteria pollutants in an urban valley of the Western United States, J. Geophys. Res.-Atmos., 123, 2684–2703, 2018. a, b

Bishop, G. A. and Stedman, D. H.: A decade of on-road emissions measurements, Environ. Sci. Technol., 42, 1651–1656, 2008. a

Bradley, K. S., Brooks, K. B., Hubbard, L. K., Popp, P. J., and Stedman, D. H.: Motor vehicle fleet emissions by OP-FTIR, Environ. Sci. Technol., 34, 897–899,, 2000. a, b

Brioude, J., Kim, S. W., Angevine, W. M., Frost, G. J., Lee, S. H., McKeen, S. A., Trainer, M., Fehsenfeld, F. C., Holloway, J. S., Ryerson, T. B., Williams, E. J., Petron, G., and Fast, J. D.: Top-down estimate of anthropogenic emission inventories and their interannual variability in Houston using a mesoscale inverse modeling technique, J. Geophys. Res.-Atmos., 116, D20305,, 2011. a

Brioude, J., Petron, G., Frost, G. J., Ahmadov, R., Angevine, W. M., Hsie, E. Y., Kim, S. W., Lee, S. H., McKeen, S. A., Trainer, M., Fehsenfeld, F. C., Holloway, J. S., Peischl, J., Ryerson, T. B., and Gurney, K. R.: A new inversion method to calculate emission inventories without a prior at mesoscale: Application to the anthropogenic CO2 emission from Houston, Texas, J. Geophys. Res.-Atmos., 117, D05312,, 2012. a, b, c, d

Brioude, J., Angevine, W. M., Ahmadov, R., Kim, S.-W., Evan, S., McKeen, S. A., Hsie, E.-Y., Frost, G. J., Neuman, J. A., Pollack, I. B., Peischl, J., Ryerson, T. B., Holloway, J., Brown, S. S., Nowak, J. B., Roberts, J. M., Wofsy, S. C., Santoni, G. W., Oda, T., and Trainer, M.: Top-down estimate of surface flux in the Los Angeles Basin using a mesoscale inverse modeling technique: assessing anthropogenic emissions of CO, NOx and CO2 and their impacts, Atmos. Chem. Phys., 13, 3661–3677,, 2013. a, b

Byers, L., Friedrich, J., Hennig, R., Kressig, A., Li, X., McCormick, C., and Malaguzzi Valeri, L.: A Global Database of Power Plants, World Resources Institute, Washington, DC, (last access: 1 July 2022), 2018. a, b, c, d

Chance, K., Liu, X., Miller, C. C., González Abad, G., Huang, G., Nowlan, C., Souri, A., Suleiman, R., Sun, K., Wang, H., Zhu, L., Zoogman, P., Al-Saadi, J., Antuña-Marrero, J. C., Carr, J., Chatfield, R., Chin, M., Cohen, R., Edwards, D., Fishman, J., Flittner, D., Geddes, J., Grutter, M., Herman, J. R., Jacob, D. J., Janz, S., Joiner, J., Kim, J., Krotkov, N. A., Lefer, B., Martin, R. V., Mayol-Bracero, O. L., Naeger, A., Newchurch, M., Pfister, G. G., Pickering, K., Pierce, R. B., Rivera Cárdenas, C., Saiz-Lopez, A., Simpson, W., Spinei, E., Spurr, R. J. D., Szykman, J. J., Torres, O., and Wang, J.: TEMPO Green Paper: Chemistry, Physics, and Meteorology Experiments with the Tropospheric Emissions: Monitoring of Pollution Instrument, in: Sensors, Systems, and Next-Generation Satellites XXIII, edited by: Neeck, S. P., Kimura, T., and Martimort, P., p. 10, SPIE, Strasbourg, France,, 2019. a

Chandra, N., Lal, S., Venkataramani, S., Patra, P. K., and Sheel, V.: Temporal variations of atmospheric CO2 and CO at Ahmedabad in western India, Atmos. Chem. Phys., 16, 6153–6173,, 2016. a

Che, K., Liu, Y., Cai, Z., Yang, D., Wang, H., Ji, D., Yang, Y., and Wang, P.: Characterization of Regional Combustion Efficiency using ΔXCO: ΔXCO2 Observed by a Portable Fourier-Transform Spectrometer at an Urban Site in Beijing, Adv. Atmos. Sci., 39, 1299–1315, 2022. a

Ching, J., Mills, G., Bechtel, B., See, L., Feddema, J., Wang, X., Ren, C., Brousse, O., Martilli, A., Neophytou, M., Mouzourides, P., Stewart, I., Hanna, A., Ng, E., Foley, M., Alexander, P., Aliaga, D., Niyogi, D., Shreevastava, A., Bhalachandran, P., Masson, V., Hidalgo, J., Fung, J., Andrade, M., Baklanov, A., Dai, W., Milcinski, G., Demuzere, M., Brunsell, N., Pesaresi, M., Miao, S., Mu, Q., Chen, F., and Theeuwes, N.: WUDAPT: An Urban Weather, Climate, and Environmental Modeling Infrastructure for the Anthropocene, B. Am. Meteorol. Soc., 99, 1907–1924,, 2018. a, b, c

Crisp, D., Fisher, B. M., O'Dell, C., Frankenberg, C., Basilio, R., Bösch, H., Brown, L. R., Castano, R., Connor, B., Deutscher, N. M., Eldering, A., Griffith, D., Gunson, M., Kuze, A., Mandrake, L., McDuffie, J., Messerschmidt, J., Miller, C. E., Morino, I., Natraj, V., Notholt, J., O'Brien, D. M., Oyafuso, F., Polonsky, I., Robinson, J., Salawitch, R., Sherlock, V., Smyth, M., Suto, H., Taylor, T. E., Thompson, D. R., Wennberg, P. O., Wunch, D., and Yung, Y. L.: The ACOS CO2 retrieval algorithm – Part II: Global XCO2 data characterization, Atmos. Meas. Tech., 5, 687–707,, 2012. a, b

Crounse, J. D., DeCarlo, P. F., Blake, D. R., Emmons, L. K., Campos, T. L., Apel, E. C., Clarke, A. D., Weinheimer, A. J., McCabe, D. C., Yokelson, R. J., Jimenez, J. L., and Wennberg, P. O.: Biomass burning and urban air pollution over the Central Mexican Plateau, Atmos. Chem. Phys., 9, 4929–4944,, 2009. a, b

de Foy, B.: City-level variations in NOx emissions derived from hourly monitoring data in Chicago, Atmos. Environ., 176, 128–139, 2018. a

Demetillo, M. A. G., Harkins, C., McDonald, B. C., Chodrow, P. S., Sun, K., and Pusede, S. E.: Space-Based Observational Constraints on NO2 Air Pollution Inequality From Diesel Traffic in Major US Cities, Geophys. Res. Lett., 48, e2021GL094333,, 2021. a

Demuzere, M., Kittner, J., Martilli, A., Mills, G., Moede, C., Stewart, I. D., van Vliet, J., and Bechtel, B.: A global map of local climate zones to support earth system modelling and urban-scale environmental science, Earth Syst. Sci. Data, 14, 3835–3873,, 2022a. a

Demuzere, M., Kittner, J., Martilli, A., Mills, G., Moede, C., Stewart, I. D., van Vliet, J., and Bechtel, B.: Global map of Local Climate Zones (1.0.0), Zenodo [data set],, 2022b. a

Djuricin, S., Pataki, D. E., and Xu, X.: A comparison of tracer methods for quantifying CO2 sources in an urban region, J. Geophys. Res.-Atmos., 115, D11303,, 2010. a

Duncan, B. N., Lamsal, L. N., Thompson, A. M., Yoshida, Y., Lu, Z., Streets, D. G., Hurwitz, M. M., and Pickering, K. E.: A space-based, high-resolution view of notable changes in urban NOx pollution around the world (2005–2014), J. Geophys. Res.-Atmos., 121, 976–996, 2016. a

Eldering, A.: OCO-3 B10 QTS Evaluation XCO2 Lite Files, Caltech Data [data set],, 2021. a, b

Eldering, A., Taylor, T. E., O'Dell, C. W., and Pavlick, R.: The OCO-3 mission: measurement objectives and expected performance based on 1 year of simulated data, Atmos. Meas. Tech., 12, 2341–2370,, 2019. a, b

ESA: Copernicus Sentinel-5P, TROPOMI Level 2 Carbon Monoxide total column products, Version 01, European Space Agency [data set],, 2018. a

Famulari, D., Nemitz, E., Di Marco, C., Phillips, G. J., Thomas, R., House, E., and Fowler, D.: Eddy-covariance measurements of nitrous oxide fluxes above a city, Agr. Forest Meteorol., 150, 786–793, 2010. a

Fasoli, B., Lin, J. C., Bowling, D. R., Mitchell, L., and Mendoza, D.: Simulating atmospheric tracer concentrations for spatially distributed receptors: updates to the Stochastic Time-Inverted Lagrangian Transport model's R interface (STILT-R version 2), Geosci. Model Dev., 11, 2813–2824,, 2018. a, b

Fujinawa, T., Kuze, A., Suto, H., Shiomi, K., Kanaya, Y., Kawashima, T., Kataoka, F., Mori, S., Eskes, H., and Tanimoto, H.: First concurrent observations of NO2 and CO2 from power plant plumes by airborne remote sensing, Geophys. Res. Lett., 48, e2021GL092685,, 2021. a

Ghasemifard, H., Vogel, F. R., Yuan, Y., Luepke, M., Chen, J., Ries, L., Leuchner, M., Schunk, C., Noreen Vardag, S., and Menzel, A.: Pollution events at the high-altitude mountain site Zugspitze-Schneefernerhaus (2670 m a.s.l.), Germany, Atmosphere, 10, 330,, 2019. a

Gonzalez, A., Millet, D. B., Yu, X., Wells, K. C., Griffis, T. J., Baier, B. C., Campbell, P. C., Choi, Y., DiGangi, J. P., Gvakharia, A., Halliday, H. S., Kort, E. A., McKain, K., Nowak, J. B., and Plant, G.: Fossil versus nonfossil CO sources in the US: New airborne constraints from ACT-America and GEM, Geophys. Res. Lett., 48, e2021GL093361,, 2021. a, b

Gurney, K. R., Patarasuk, R., Liang, J., Song, Y., O'Keeffe, D., Rao, P., Whetstone, J. R., Duren, R. M., Eldering, A., and Miller, C.: The Hestia fossil fuel CO2 emissions data product for the Los Angeles megacity (Hestia-LA), Earth Syst. Sci. Data, 11, 1309–1335,, 2019. a, b

Hakkarainen, J., Szelag, M. E., Ialongo, I., Retscher, C., Oda, T., and Crisp, D.: Analyzing nitrogen oxides to carbon dioxide emission ratios from space: A case study of Matimba Power Station in South Africa, Atmospheric Environment: X, 10, 100110,, 2021. a

Han, S., Kondo, Y., Oshima, N., Takegawa, N., Miyazaki, Y., Hu, M., Lin, P., Deng, Z., Zhao, Y., Sugimoto, N., and Wu, Y.: Temporal variations of elemental carbon in Beijing, J. Geophys. Res.-Atmos., 114, D23202,, 2009. a

Harrison, R. M., Dall'Osto, M., Beddows, D. C. S., Thorpe, A. J., Bloss, W. J., Allan, J. D., Coe, H., Dorsey, J. R., Gallagher, M., Martin, C., Whitehead, J., Williams, P. I., Jones, R. L., Langridge, J. M., Benton, A. K., Ball, S. M., Langford, B., Hewitt, C. N., Davison, B., Martin, D., Petersson, K. F., Henshaw, S. J., White, I. R., Shallcross, D. E., Barlow, J. F., Dunbar, T., Davies, F., Nemitz, E., Phillips, G. J., Helfter, C., Di Marco, C. F., and Smith, S.: Atmospheric chemistry and physics in the atmosphere of a developed megacity (London): an overview of the REPARTEE experiment and its conclusions, Atmos. Chem. Phys., 12, 3065–3114,, 2012. a

Haszpra, L., Ferenczi, Z., and Barcza, Z.: Estimation of greenhouse gas emission factors based on observed covariance of CO2, CH4, N2O and CO mole fractions, Environmental Sciences Europe, 31, 95,, 2019. a

Hedelius, J. K., Viatte, C., Wunch, D., Roehl, C. M., Toon, G. C., Chen, J., Jones, T., Wofsy, S. C., Franklin, J. E., Parker, H., Dubey, M. K., and Wennberg, P. O.: Assessment of errors and biases in retrievals of XCO2, XCH4, XCO, and XN2O from a 0.5 cm−1 resolution solar-viewing spectrometer, Atmos. Meas. Tech., 9, 3527–3546,, 2016. a, b

Hedelius, J. K., Liu, J., Oda, T., Maksyutov, S., Roehl, C. M., Iraci, L. T., Podolske, J. R., Hillyard, P. W., Liang, J., Gurney, K. R., Wunch, D., and Wennberg, P. O.: Southern California megacity CO2, CH4, and CO flux estimates using ground- and space-based remote sensing and a Lagrangian model, Atmos. Chem. Phys., 18, 16271–16291,, 2018. a, b, c

Huang, X., Wang, T., Talbot, R., Xie, M., Mao, H., Li, S., Zhuang, B., Yang, X., Fu, C., and Zhu, J.: Temporal characteristics of atmospheric CO2 in urban Nanjing, China, Atmos. Res., 153, 437–450, 2015. a

Hudman, R. C., Murray, L. T., Jacob, D. J., Millet, D., Turquety, S., Wu, S., Blake, D., Goldstein, A., Holloway, J., and Sachse, G. W.: Biogenic versus anthropogenic sources of CO in the United States, Geophys. Res. Lett., 35, L04801,, 2008. a

Kaiser, J. W., Heil, A., Andreae, M. O., Benedetti, A., Chubarova, N., Jones, L., Morcrette, J.-J., Razinger, M., Schultz, M. G., Suttie, M., and van der Werf, G. R.: Biomass burning emissions estimated with a global fire assimilation system based on observed fire radiative power, Biogeosciences, 9, 527–554,, 2012. a

Kerr, G. H., Goldberg, D. L., and Anenberg, S. C.: COVID-19 pandemic reveals persistent disparities in nitrogen dioxide pollution, P. Natl. Acad. Sci., 118, e2022409118,, 2021. a

Kiel, M., Eldering, A., Roten, D. D., Lin, J. C., Feng, S., Lei, R., Lauvaux, T., Oda, T., Roehl, C. M., Blavier, J.-F., et al.: Urban-focused satellite CO2 observations from the Orbiting Carbon Observatory-3: A first look at the Los Angeles megacity, Remote Sens. Environ., 258, 112314,, 2021. a, b

Lama, S., Houweling, S., Boersma, K. F., Eskes, H., Aben, I., Denier van der Gon, H. A. C., Krol, M. C., Dolman, H., Borsdorff, T., and Lorente, A.: Quantifying burning efficiency in megacities using the NO2/CO ratio from the Tropospheric Monitoring Instrument (TROPOMI), Atmos. Chem. Phys., 20, 10295–10310,, 2020. a, b

Landgraf, J., Borsdorff, T., Langerock, B., and Keppens, A.: S5P Mission Performance Centre Carbon Monoxide [L2 CO] Readme V1.4, Tech. Rep., Netherlands Institute for Space Research (SRON),, last access: 16 September 2020. a

Laughner, J. L., Neu, J. L., Schimel, D., Wennberg, P. O., Barsanti, K., Bowman, K. W., Chatterjee, A., Croes, B. E., Fitzmaurice, H. L., Henze, D. K., Kim, J., Kort, E. A., Liu, Z., Miyazaki, K., Turner, A. J., Anenberg, S., Avise, J., Cao, H., Crisp, D., de Gouw, J., Eldering, A., Fyfe, J. C., Goldberg, D. L., Gurney, K. R., Hasheminassab, S., Hopkins, F., Ivey, C. E., Jones, D. B. A., Liu, J., Lovenduski, N. S., Martin, R. V., McKinley, G. A., Ott, L., Poulter, B., Ru, M., Sander, S. P., Swart, N., Yung, Y. L., and Zeng, Z. C.: Societal shifts due to COVID-19 reveal large-scale complexities and feedbacks between atmospheric chemistry and climate change, P. Natl. Acad. Sci. USA, 118, e210948118,, 2021. a

Lei, R., Feng, S., Danjou, A., Broquet, G., Wu, D., Lin, J. C., O'Dell, C. W., and Lauvaux, T.: Fossil fuel CO2 emissions over metropolitan areas from space: A multi-model analysis of OCO-2 data over Lahore, Pakistan, Remote Sens. Environ., 264, 112625,, 2021. a

Lin, J. and Gerbig, C.: Accounting for the effect of transport errors on tracer inversions, Geophys. Res. Lett., 32, L01802,, 2005. a

Lin, J., Gerbig, C., Wofsy, S., Andrews, A., Daube, B., Davis, K., and Grainger, C.: A near-field tool for simulating the upstream influence of atmospheric observations: The Stochastic Time-Inverted Lagrangian Transport (STILT) model, J. Geophys. Res.-Atmos., 108, 4493,, 2003. a, b

Lin, J. C., Mitchell, L., Crosman, E., Mendoza, D. L., Buchert, M., Bares, R., Fasoli, B., Bowling, D. R., Pataki, D., Catharine, D., Strong, C., Gurney, K. R., Patarasuk, R., Baasandorj, M., Jacques, A., Hoch, S., Horel, J., and Ehleringer, J.: CO2 and carbon emissions from cities linkages to air quality, socioeconomic activity, and stakeholders in the Salt Lake city urban area, B. Am. Meteorol. Soc., 99, 2325–2339,, 2018. a

Lin, J. C., Bares, R., Fasoli, B., Garcia, M., Crosman, E., and Lyman, S.: Declining methane emissions and steady, high leakage rates observed over multiple years in a western US oil/gas production basin, Scientific Reports, 11, 22291,, 2021. a

Lindenmaier, R., Dubey, M. K., Henderson, B. G., Butterfield, Z. T., Herman, J. R., Rahn, T., and Lee, S. H.: Multiscale observations of CO2, 13CO2, and pollutants at Four Corners for emission verification and attribution, P. Natl. Acad. Sci. USA, 111, 8386–8391,, 2014. a, b, c

Lopez, M., Schmidt, M., Delmotte, M., Colomb, A., Gros, V., Janssen, C., Lehman, S. J., Mondelain, D., Perrussel, O., Ramonet, M., Xueref-Remy, I., and Bousquet, P.: CO, NOx and 13CO2 as tracers for fossil fuel CO2: results from a pilot study in Paris during winter 2010, Atmos. Chem. Phys., 13, 7343–7358,, 2013. a

Makarova, M. V., Alberti, C., Ionov, D. V., Hase, F., Foka, S. C., Blumenstock, T., Warneke, T., Virolainen, Y. A., Kostsov, V. S., Frey, M., Poberovskii, A. V., Timofeyev, Y. M., Paramonova, N. N., Volkova, K. A., Zaitsev, N. A., Biryukov, E. Y., Osipov, S. I., Makarov, B. K., Polyakov, A. V., Ivakhov, V. M., Imhasin, H. Kh., and Mikhailov, E. F.: Emission Monitoring Mobile Experiment (EMME): an overview and first results of the St. Petersburg megacity campaign 2019, Atmos. Meas. Tech., 14, 1047–1073,, 2021. a

Miller, S. M., Matross, D. M., Andrews, A. E., Millet, D. B., Longo, M., Gottlieb, E. W., Hirsch, A. I., Gerbig, C., Lin, J. C., Daube, B. C., Hudman, R. C., Dias, P. L. S., Chow, V. Y., and Wofsy, S. C.: Sources of carbon monoxide and formaldehyde in North America determined from high-resolution atmospheric data, Atmos. Chem. Phys., 8, 7673–7696,, 2008. a

Mitchell, L. E., Lin, J. C., Bowling, D. R., Pataki, D. E., Strong, C., Schauer, A. J., Bares, R., Bush, S. E., Stephens, B. B., Mendoza, D., et al.: Long-term urban carbon dioxide observations reveal spatial and temporal dynamics related to urban characteristics and growth, P. Natl. Acad. Sci. USA, 115, 2912–2917, 2018. a

Moldanová, J., Fridell, E., Popovicheva, O., Demirdjian, B., Tishkova, V., Faccinetto, A., and Focsa, C.: Characterisation of particulate matter and gaseous emissions from a large ship diesel engine, Atmos. Environ., 43, 2632–2641, 2009. a

Nathan, B., Lauvaux, T., Turnbull, J., and Gurney, K.: Investigations into the use of multi-species measurements for source apportionment of the Indianapolis fossil fuel CO2 signal, Elementa, 6, 21,, 2018. a, b, c

O'Shea, S. J., Allen, G., Fleming, Z. L., Bauguitte, S. J.-B., Percival, C. J., Gallagher, M. W., Lee, J., Helfter, C., and Nemitz, E.: Area fluxes of carbon dioxide, methane, and carbon monoxide derived from airborne measurements around Greater London: A case study during summer 2012, J. Geophys. Res.-Atmos., 119, 4940–4952, 2014. a

Palmer, P. I., Suntharalingam, P., Jones, D. B., Jacob, D. J., Streets, D. G., Fu, Q., Vay, S. A., and Sachse, G. W.: Using CO2: CO correlations to improve inverse analyses of carbon fluxes, J. Geophys. Res.-Atmos., 111, D12318,, 2006. a, b, c

Park, H., Jeong, S., Park, H., Labzovskii, L. D., and Bowman, K. W.: An assessment of emission characteristics of Northern Hemisphere cities using spaceborne observations of CO2, CO, and NO2, Remote Sens. Environ., 254, 112246,, 2021. a, b

Pitt, J. R., Allen, G., Bauguitte, S. J.-B., Gallagher, M. W., Lee, J. D., Drysdale, W., Nelson, B., Manning, A. J., and Palmer, P. I.: Assessing London CO2, CH4 and CO emissions using aircraft measurements and dispersion modelling, Atmos. Chem. Phys., 19, 8931–8945,, 2019. a

Plant, G., Kort, E. A., Floerchinger, C., Gvakharia, A., Vimont, I., and Sweeney, C.: Large fugitive methane emissions from urban centers along the US East Coast, Geophys. Res. Lett., 46, 8500–8507, 2019. a

Popa, M. E., Vollmer, M. K., Jordan, A., Brand, W. A., Pathirana, S. L., Rothe, M., and Röckmann, T.: Vehicle emissions of greenhouse gases and related tracers from a tunnel study: CO:CO2, N2O:CO2, CH4:CO2, O2:CO2 ratios, and the stable isotopes 13C and 18O in CO2 and CO, Atmos. Chem. Phys., 14, 2105–2123,, 2014. a, b, c

Reuter, M., Buchwitz, M., Schneising, O., Krautwurst, S., O'Dell, C. W., Richter, A., Bovensmann, H., and Burrows, J. P.: Towards monitoring localized CO2 emissions from space: co-located regional CO2 and NO2 enhancements observed by the OCO-2 and S5P satellites, Atmos. Chem. Phys., 19, 9371–9383,, 2019. a, b

Schneising, O., Buchwitz, M., Reuter, M., Bovensmann, H., Burrows, J. P., Borsdorff, T., Deutscher, N. M., Feist, D. G., Griffith, D. W. T., Hase, F., Hermans, C., Iraci, L. T., Kivi, R., Landgraf, J., Morino, I., Notholt, J., Petri, C., Pollard, D. F., Roche, S., Shiomi, K., Strong, K., Sussmann, R., Velazco, V. A., Warneke, T., and Wunch, D.: A scientific algorithm to simultaneously retrieve carbon monoxide and methane from TROPOMI onboard Sentinel-5 Precursor, Atmos. Meas. Tech., 12, 6771–6802,, 2019. a

Schuh, A. E., Otte, M., Lauvaux, T., and Oda, T.: Far-field biogenic and anthropogenic emissions as a dominant source of variability in local urban carbon budgets: A global high-resolution model study with implications for satellite remote sensing, Remote Sens. Environ., 262, 112473, 2021. a

Shekhar, A., Chen, J., Paetzold, J. C., Dietrich, F., Zhao, X., Bhattacharjee, S., Ruisinger, V., and Wofsy, S. C.: Anthropogenic CO2 emissions assessment of Nile Delta using XCO2 and SIF data from OCO-2 satellite, Environ. Res. Lett., 15, 095010,, 2020. a

Silva, S. and Arellano, A.: Characterizing Regional-Scale Combustion Using Satellite Retrievals of CO, NO2 and CO2, Remote Sensing, 9, 744,, 2017. a, b, c, d, e

Silva, S. J., Arellano, A. F., and Worden, H. M.: Toward anthropogenic combustion emission constraints from space-based analysis of urban CO2/CO sensitivity, Geophys. Res. Lett., 40, 4971–4976,, 2013. a

Sim, S., Jeong, S., Park, H., Park, C., Kwak, K. H., Lee, S. B., Kim, C. H., Lee, S., Chang, J. S., Kang, H., and Woo, J. H.: Co-benefit potential of urban CO2 and air quality monitoring: A study on the first mobile campaign and building monitoring experiments in Seoul during the winter, Atmos. Pollut. Res., 11, 1963–1970,, 2020. a

Solazzo, E., Crippa, M., Guizzardi, D., Muntean, M., Choulga, M., and Janssens-Maenhout, G.: Uncertainties in the Emissions Database for Global Atmospheric Research (EDGAR) emission inventory of greenhouse gases, Atmos. Chem. Phys., 21, 5655–5683,, 2021. a, b

Stewart, I. D. and Oke, T. R.: Local climate zones for urban temperature studies, B. Am. Meteorol. Soc., 93, 1879–1900, 2012. a, b, c

Super, I., van der Gon, H. A. D., Visschedijk, A. J., Moerman, M. M., Chen, H., van der Molen, M. K., and Peters, W.: Interpreting continuous in-situ observations of carbon dioxide and carbon monoxide in the urban port area of Rotterdam, Atmos. Pollut. Res., 8, 174–187,, 2017. a, b

Surl, L., Palmer, P. I., and González Abad, G.: Which processes drive observed variations of HCHO columns over India?, Atmos. Chem. Phys., 18, 4549–4566,, 2018. a, b

Tang, W., Arellano, A. F., DiGangi, J. P., Choi, Y., Diskin, G. S., Agustí-Panareda, A., Parrington, M., Massart, S., Gaubert, B., Lee, Y., Kim, D., Jung, J., Hong, J., Hong, J.-W., Kanaya, Y., Lee, M., Stauffer, R. M., Thompson, A. M., Flynn, J. H., and Woo, J.-H.: Evaluating high-resolution forecasts of atmospheric CO and CO2 from a global prediction system during KORUS-AQ field campaign, Atmos. Chem. Phys., 18, 11007–11030,, 2018. a

Tang, W., Gaubert, B., Emmons, L., Choi, Y., DiGangi, J. P., Diskin, G. S., Xu, X., He, C., Worden, H., Tilmes, S., Buchholz, R., Halliday, H. S., and Arellano, A. F.: On the relationship between tropospheric CO and CO2 during KORUS-AQ and its role in constraining anthropogenic CO2, Atmos. Chem. Phys. Discuss. [preprint],, 2020. a

Taylor, T. E., Eldering, A., Merrelli, A., Kiel, M., Somkuti, P., Cheng, C., Rosenberg, R., Fisher, B., Crisp, D., Basilio, R., et al.: OCO-3 early mission operations and initial (vEarly) XCO2 and SIF retrievals, Remote Sens. Environ., 251, 112032,, 2020. a, b, c

Turnbull, J. C., Karion, A., Fischer, M. L., Faloona, I., Guilderson, T., Lehman, S. J., Miller, B. R., Miller, J. B., Montzka, S., Sherwood, T., Saripalli, S., Sweeney, C., and Tans, P. P.: Assessment of fossil fuel carbon dioxide and other anthropogenic trace gas emissions from airborne measurements over Sacramento, California in spring 2009, Atmos. Chem. Phys., 11, 705–721,, 2011a. a

Turnbull, J. C., Tans, P. P., Lehman, S. J., Baker, D., Conway, T. J., Chung, Y. S., Gregg, J., Miller, J. B., Southon, J. R., and Zhou, L. X.: Atmospheric observations of carbon monoxide and fossil fuel CO2 emissions from East Asia, J. Geophys. Res.-Atmos., 116, D24306,, 2011b. a, b

Turnbull, J. C., Sweeney, C., Karion, A., Newberger, T., Lehman, S. J., Tans, P. P., Davis, K. J., Lauvaux, T., Miles, N. L., Richardson, S. J., Cambaliza, M. O., Shepson, P. B., Gurney, K., Patarasuk, R., and Razlivanov, I.: Toward quantification and source sector identification of fossil fuel CO2 emissions from an urban area: Results from the INFLUX experiment, J. Geophys. Res.-Atmos., 120, 292–312,, 2015. a

Turner, A. J., Köhler, P., Magney, T. S., Frankenberg, C., Fung, I., and Cohen, R. C.: A double peak in the seasonality of California's photosynthesis as observed from space, Biogeosciences, 17, 405–422,, 2020. a

United Nations, Department of Economic and Social Affairs, and Population Division: World Urbanization Prospects: The 2018 Revision (ST/ESA/SER.A/420), United Nations, New York, 2019. a

Veefkind, J. P., Aben, I., McMullan, K., Förster, H., de Vries, J., Otter, G., Claas, J., Eskes, H. J., de Haan, J. F., Kleipool, Q., van Weele, M., Hasekamp, O., Hoogeveen, R., Landgraf, J., Snel, R., Tol, P., Ingmann, P., Voors, R., Kruizinga, B., Vink, R., Visser, H., and Levelt, P. F.: TROPOMI on the ESA Sentinel-5 Precursor: A GMES mission for global observations of the atmospheric composition for climate, air quality and ozone layer applications, Remote Sens. Environ., 120, 70–83,, 2012. a, b, c

Venables, W. N. and Ripley, B. D.: Modern Applied Statistics with S, fourth edn., Springer, New York, ISBN 0-387-95457-0, (last access: 1 July 2022​​​​​​​), 2002. a

Vollmer, M. K., Juergens, N., Steinbacher, M., Reimann, S., Weilenmann, M., and Buchmann, B.: Road vehicle emissions of molecular hydrogen (H2) from a tunnel study, Atmos. Environ., 41, 8355–8369, 2007. a

Wang, H., Jacob, D. J., Kopacz, M., Jones, D. B. A., Suntharalingam, P., Fisher, J. A., Nassar, R., Pawson, S., and Nielsen, J. E.: Error correlation between CO2 and CO as constraint for CO2 flux inversions using satellite data, Atmos. Chem. Phys., 9, 7313–7323,, 2009. a

Wang, X., Lei, Y., Yan, L., Liu, T., Zhang, Q., and He, K.: A unit-based emission inventory of SO2, NOx and PM for the Chinese iron and steel industry from 2010 to 2015, Sci. Total Environ., 676, 18–30, 2019. a, b

Wang, Y., Munger, J. W., Xu, S., McElroy, M. B., Hao, J., Nielsen, C. P., and Ma, H.: CO2 and its correlation with CO at a rural site near Beijing: implications for combustion efficiency in China, Atmos. Chem. Phys., 10, 8881–8897,, 2010. a

Wennberg, P. O., Mui, W., Wunch, D., Kort, E. A., Blake, D. R., Atlas, E. L., Santoni, G. W., Wofsy, S. C., Diskin, G. S., Jeong, S., and Fischer, M. L.: On the sources of methane to the Los Angeles atmosphere, Environ. Sci. Technol., 46, 9282–9289,, 2012. a, b, c, d

Wennberg, P. O., Wunch, D., Roehl, C., Blavier, J.-F., Toon, G. C., and Allen, N.: TCCON data from Caltech (US), Release GGG2020R0, TCCON Data Archive [data set],, 2017. a

Westerdahl, D., Wang, X., Pan, X., and Zhang, K. M.: Characterization of on-road vehicle emission factors and microenvironmental air quality in Beijing, China, Atmos. Environ., 43, 697–705, 2009. a, b, c

Williams, E., Lerner, B., Murphy, P., Herndon, S., and Zahniser, M.: Emissions of NOx, SO2, CO, and HCHO from commercial marine shipping during Texas Air Quality Study (TexAQS) 2006, J. Geophys. Res.-Atmos., 114, D21306,, 2009. a

Wu, D. and Lin, J. C.: Urban Biogenic CO2 fluxes: GPP, Reco and NEE Estimates from SMUrF, 2010–2019, ORNL DAAC, Oak Ridge, Tennessee, USA [data set],, 2021. a

Wu, D., Lin, J. C., Fasoli, B., Oda, T., Ye, X., Lauvaux, T., Yang, E. G., and Kort, E. A.: A Lagrangian approach towards extracting signals of urban CO2 emissions from satellite observations of atmospheric column CO2 (XCO2): X-Stochastic Time-Inverted Lagrangian Transport model (“X-STILT v1”), Geosci. Model Dev., 11, 4843–4871,, 2018. a, b, c

Wu, D., Fasoli, B., and Lin, J. C.: uataq/X-STILT: X-STILT (v1.4.1), Zenodo [data set],, 2019. a

Wu, D., Lin, J. C., Oda, T., and Kort, E. A.: Space-based quantification of per capita CO2 emissions from cities, Environ. Res. Lett., 15, 035004,, 2020. a, b, c

Wu, D., Lin, J. C., Duarte, H. F., Yadav, V., Parazoo, N. C., Oda, T., and Kort, E. A.: A model for urban biogenic CO2 fluxes: Solar-Induced Fluorescence for Modeling Urban biogenic Fluxes (SMUrF v1), Geosci. Model Dev., 14, 3633–3661,, 2021. a, b, c

Wunch, D., Wennberg, P., Toon, G., Keppel-Aleks, G., and Yavin, Y.: Emissions of greenhouse gases from a North American megacity, Geophys. Res. Lett., 36, L15810,, 2009. a, b

Xia, L., Zhang, G., Liu, L., Li, B., Zhan, M., Kong, P., and Wang, H.: Atmospheric CO2 and CO at Jingdezhen station in central China: Understanding the regional transport and combustion efficiency, Atmos. Environ., 222, 117104,, 2020. a

Ye, X., Lauvaux, T., Kort, E. A., Oda, T., Feng, S., Lin, J. C., Yang, E. G., and Wu, D.: Constraining Fossil Fuel CO2 Emissions From Urban Area Using OCO-2 Observations of Total Column CO2, J. Geophys. Res.-Atmos., 125, e2019JD030528,, 2020.  a, b

Yokota, T., Yoshida, Y., Eguchi, N., Ota, Y., Tanaka, T., Watanabe, H., and Maksyutov, S.: Global concentrations of CO2 and CH4 retrieved from GOSAT: First preliminary results, Sola, 5, 160–163, 2009. a

Yuan, L. and Smith, A. C.: CO and CO2 emissions from spontaneous heating of coal under different ventilation rates, Int. J. Coal Geol., 88, 24–30, 2011. a

Zhang, F., Chen, Y., Tian, C., Lou, D., Li, J., Zhang, G., and Matthias, V.: Emission factors for gaseous and particulate pollutants from offshore diesel engine vessels in China, Atmos. Chem. Phys., 16, 6319–6334,, 2016. a

Zhang, Y., Smith, S. J., Bowden, J. H., Adelman, Z., and West, J. J.: Co-benefits of global, domestic, and sectoral greenhouse gas mitigation for US air quality and human health in 2050, Environ. Res. Lett., 12, 114033,, 2017. a

Zhang, Y., Joiner, J., Alemohammad, S. H., Zhou, S., and Gentine, P.: A global spatially contiguous solar-induced fluorescence (CSIF) dataset using neural networks, Biogeosciences, 15, 5779–5800,, 2018. a

Zhu, L., Jacob, D. J., Mickley, L. J., Marais, E. A., Cohan, D. S., Yoshida, Y., Duncan, B. N., Abad, G. G., and Chance, K. V.: Anthropogenic emissions of highly reactive volatile organic compounds in eastern Texas inferred from oversampling of satellite (OMI) measurements of HCHO columns, Environ. Res. Lett., 9, 114004,, 2014. a

Short summary
Prior studies have derived the combustion efficiency for a region/city using observed CO2 and CO. We further zoomed into the urban domain and accounted for factors affecting the calculation of spatially resolved combustion efficiency from two satellites. The intra-city variability in combustion efficiency was linked to heavy industry within Shanghai and LA without relying on emission inventories. Such an approach can be applied when analyzing data from future geostationary satellites.
Final-revised paper