1 Introduction
On November 17, 2023, the ERA5 (ECMWF Reanalysis v5) data (
Hersbach et al., 2020) recorded a temporary exceedance of the 2 °C global warming, underscoring the urgent need to address climate change. Over recent decades, the detrimental effects of climate change have become increasingly evident, with a marked rise in the frequency and intensity of extreme weather events, including downpours, droughts, and wildfires (
Swain et al., 2020;
Newman and Noy, 2023;
Ombadi et al., 2023;
Yuan et al., 2023). In response to this critical global challenge, most nations have committed to reducing greenhouse gas emissions through Nationally Determined Contributions (NDCs) as stipulated by the Paris Agreement. These commitments, updated every five years since 2020, constitute a crucial part of the global effort to combat climate change (
Röser et al., 2020). To timely evaluate collective progress toward limiting global warming to below 2 °C (ideally 1.5 °C) and motivate parties to ratchet up their NDCs, the Global Stocktake (GST) was initiated in 2023 and will recur every five years (
Editorial, 2023). Rigorous and comprehensive implementation of NDCs could potentially constrain global warming to just under 2 °C, which emphasizes the necessity for intensified and unequivocal actions from all participating entities (
Meinshausen et al., 2022).
To bolster climate ambitions, the development of anthropogenic carbon dioxides (CO
2) Emissions Monitoring & Verification Support (MVS) capacity is essential (
European et al., 2017). Anthropogenic CO
2 emissions are primarily categorized into two sources: fossil fuel combustion and land use, accounting for 88.1% and 11.9% of emissions, respectively, according to the Global Carbon Budget 2023 (
Friedlingstein et al., 2023). Given the dominance of fossil fuel CO
2 emissions (FFCO
2), the quest for a thorough comprehensive understanding of FFCO
2 emissions has spurred numerous research endeavors (
Jones et al., 2023;
Liu et al., 2023). The 2006 IPCC Guidelines for National Greenhouse Gas Inventories (
IPCC, 2006) identified bottom-up inventory compilation as the primary method for countries to report their estimated emissions. However, the bottom-up method entails substantial labor and time investment in acquiring detailed data on activity levels and emission factors. These requirements pose challenges for many countries, particularly those grappling with the development of reliable statistical databases. In contrast, advancements in satellite remote sensing technologies have provided indispensable tools for monitoring changes in atmospheric composition concentrations (
Yang et al., 2013). Concurrently, inversion algorithms for deriving emissions from satellite-observed column concentrations have flourished. In this context, the top-down emission inversion was proposed as an alternative to the bottom-up approach in the 2019 Refinement to 2006 IPCC Guidelines (
Eduardo Calvo Buendia et al., 2019). The top-down inversed emissions have been currently recognized and adopted as a viable method for international climate commitments.
Researchers have dedicated considerable effort to developing top-down FFCO
2 emission inversions directly from CO
2 satellite observation, broadly categorized into data-driven (e.g., cross-sectional flux method (CSF)) and model-driven approaches (e.g., based on Lagrangian or Eulerian models) (
Nehrkorn et al., 2010;
Nassar et al., 2017). CO
2-based inversions have been widely adopted to infer FFCO
2 emissions from point and city scales, demonstrating their feasibility. However, persistent challenges remain despite these advances. For instance, natural fluxes exhibit substantial spatial and seasonal fluctuations, whereas FFCO
2 emissions remain relatively stable, complicating the detection of CO
2 enhancement signals from human activities (
Buchwitz et al., 2021;
Byrne et al., 2023). To address these limitations, integrating co-emitted short-lived proxies such as nitrogen dioxides (NO
2) into CO
2 emission inversions has recently emerged as an effective approach (
Berezin et al., 2013;
Reuter et al., 2014).
Amidst the rapid advancements in FFCO2 emission inversion methodologies, there is a critical need for a comprehensive and systematic review to consolidate existing knowledge and guide future research endeavors. This review aims to synthesize the current research progress and challenges in FFCO2 emission inversions. It is noteworthy that, unless otherwise specified, the CO2 emissions in this study refer to FFCO2 emissions. This paper is organized as follows: Section 2 provides an overview of CO2 emission inversion methods from CO2 satellite observations including both advantages and limitations; Section 3 summarizes CO2 emission inversion methods integrating NO2 retrievals, detailing the principles and applications; Section 4 examines the inherent limitations and uncertainties within current NO2-based CO2 emission inversion frameworks; Section 5 outlines potential future research directions for both technology and methodology for this field; and finally, Section 6 offers conclusions.
2 Utilizing satellite CO2 data to monitor CO2 emissions
2.1 CO2 monitoring satellites
Over the past two decades, substantial progress has been made in satellite monitoring of CO
2 column concentrations, providing a robust data foundation for CO
2 emission inversion (Tab.1). The European Space Agency (ESA) launched the Environmental Satellite (Envisat) in 2002, equipped with the Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY), enabling the monitoring of various atmospheric compositions (
Burrows et al., 1995;
Bovensmann et al., 1999). Until 2009, Envisat remained the sole satellite capable of measuring column-averaged dry-air mole fractions of carbon dioxide (XCO
2), providing long-term data on the global distribution of atmospheric CO
2 concentrations (
Schneising et al., 2011). In 2009, the Japan Aerospace Exploration Agency launched the Greenhouse gases Observing SATellite (GOSAT), the first satellite dedicated to monitoring CO
2 and methane (CH
4) concentration (
Kuze et al., 2009). Building on the success of GOSAT, Japan launched the GOSAT-2 in 2018, featuring the advanced Fourier Transform Spectrometer-2 (TANSO-FTS-2) and the Cloud and Aerosol Imager-2 (TANSO-CAI-2), substantially expanding coverage of the shortwave infrared spectrum. These advancements are particularly beneficial for reducing interference in CO
2 inversion caused by clouds and aerosols, thereby improving observation accuracy (
Suto et al., 2021).
Since 2014, satellites have been substantially advanced in achieving finer resolution and reducing uncertainty, thereby enhancing their ability to detect CO
2 signals. The Orbital Carbon Observing-2 (OCO-2), launched by the National Aeronautics and Space Administration (NASA) in 2014, utilizes three high-resolution grating spectrometers to achieve precise CO
2 measurements (
Crisp et al., 2017). It plays a pivotal role in assessing the gross primary productivity and carbon budgets of terrestrial ecosystems by revealing seasonal and regional variations in CO
2 concentrations and Solar-induced chlorophyll fluorescence (SIF) (
Sun et al., 2017). OCO-2 can provide unique insights into the complex interactions between climate phenomena (e.g., El Niño events) and natural disasters (e.g., wildfires) within the global carbon cycle (
Chatterjee et al., 2017;
Liu et al., 2017). Furthermore, it demonstrates the potential to characterize CO
2 emission hotspots from both anthropogenic and natural sources (
Eldering et al., 2017;
Heymann et al., 2017;
Schwandner et al., 2017). Despite these technological advancements, uncertainty persists among different satellites, such as GOSAT and OCO-2, which differ in spatial coverage, data amount, and retrieval biases (
Wang et al., 2019). To enhance the monitoring capabilities of CO
2 in the atmosphere, NASA launched and hosted the Orbiting Carbon Observatory-3 (OCO-3) sensor on the International Space Station in 2019. Its Snapshot Area Map (SAM) and target mode measurements can scan larger contiguous areas and capture plume structures over city areas more effectively compared to OCO-2 (
Eldering et al., 2019).
The Chinese CO
2 satellite TanSat was launched in 2016. The CO
2 retrieval data was first utilized to map global CO
2 distribution (
Yang et al., 2018). Liu et al. (
2018) pioneered efforts to validate TanSat retrieval against eight ground-based measurements from the Total Carbon Column Observing Network (TCCON), demonstrating a precision of 2.11 ppm. Subsequent endeavors have persistently aimed at refining TanSat data precision through spectral calibration and algorithm optimization, achieving retrievals with a precision of 1.47 ppm, comparable to OCO-2 (
Yang et al., 2020;
Hong et al., 2022). Building upon these advancements,
Hong et al. (
2024) further improved TanSat XCO
2 retrievals by implementing spectral recalibration, spectral window optimization, and explicit radiative transfer simulation. This work extended the utility of TanSat data to capture CO
2 enhancements over the ocean and quantify emission flux.
2.2 Data- and model-driven approaches
The initial mission of CO
2 satellites is to assess carbon budgets on a large spatial scale, while advancements in spatial resolution and signal-to-noise ratios have demonstrated their potential for quantifying FFCO
2 local emissions (
Wu et al., 2020;
Nassar et al., 2021). Currently, there are two main methods for estimating FFCO
2 emissions from CO
2 satellites: data-driven and model-driven approaches (Tab.2).
Data-driven methods typically involve processing satellite XCO
2 along with local wind field information, which is suitable for estimating emissions from large point sources (e.g., power plants) and isolated cities assuming steady-state conditions (
Velazco et al., 2011;
Nassar et al., 2022). These methods primarily include four categories (Tab.2), generally starting with the detection of CO
2 plumes from satellite imagery to determine regions for emission quantification and attribution. Plumes can be identified using segmentation algorithms (
Dumont Le Brazidec et al., 2023) or thresholding based on signal-to-noise ratio (
Kuhlmann et al., 2021). Plume positioning can also be achieved using known source locations and wind fields (
Lin et al., 2023). Quantifying XCO
2 enhancements above the background entails carefully defining and deducting the background levels. Various approaches have been employed for this purpose, such as smoothing the background field with low-pass filters (
Thoning et al., 1989), estimating the background based on the upwind direction of point source emissions (
Kort et al., 2012), and incorporating background term into a linear polynomial (
Reuter et al., 2019). The identified XCO
2 enhancements above the background levels are finally combined with local wind field data to estimate emissions related to these enhancements. Data Science and Artificial intelligence (AI) techniques, such as clustering analysis and convolutional neural network (CNN), have been increasingly integrated into data-driven methods to extract emission information from concentration distributions and enhance the efficiency of data processing (
Schuit et al., 2023).
Model-driven methods generally involve the utilization of Eulerian (e.g., Weather Research and Forecasting-Chem (WRF-Chem)) (
Lei et al., 2021) and Lagrangian models (e.g., Column-Stochastic Time-Inverted Lagrangian Transport model (X-STILT)) (
Wu et al., 2020), which are commonly employed for city or regional-scale CO
2 emission inversion. Eulerian models utilize fixed grid cells and a known prior emission inventory to simulate XCO
2 by forwardly transporting emitted CO
2 molecules through the simulation domain. Lagrangian models trace the backward trajectory of air parcels from a given receptor location (e.g., the position of observed XCO
2 enhancement from satellite) to determine the source emission area or footprint contributing to the concentration at the receptor. These models simulate the transport of CO
2 and establish the relationship between CO
2 emissions and satellite-observed XCO
2 concentrations.
Pillai et al. (2016) utilized the Bayesian inversion approach, which minimizes the cost function accounting for model- and observation-specific errors, yielding a posterior estimate of CO
2 emissions.
Zheng et al. (2019) applied a maximum likelihood estimation approach to directly estimate CO
2 emissions from XCO
2 observations without prior emission information, by calculating the Jacobian matrix.
2.3 Application in point source and city emissions
Satellite XCO
2 images have been utilized to estimate CO
2 emissions from both point sources and city scales. Regarding point sources,
Bovensmann et al. (2010) proposed a conceptual technique (GP method) to quantify CO
2 emissions from power plants based on the Observing System Simulation Experiment (OSSE).
Nassar et al. (2017) quantified CO
2 emissions from 7 large and middle-sized power plants based on OCO-2 CO
2 observations and subsequently improved their methods to capture emissions from more point sources (
Nassar et al., 2021). Utilizing a high-resolution WRF-Chem model with a spatial resolution of 1.33 km × 1.33 km,
Zheng et al. (2019) successfully captured emissions from 7 point sources identified by
Nassar et al. (2017), without prior information on emissions from these point sources. Considering the difficulty of accurately detecting XCO
2 enhancements,
Reuter et al. (2019) identified XCO
2 enhancements using NO
2 tropospheric vertical column densities (TVCDs) from TROPOspheric Monitoring Instrument (TROPOMI) as a proxy, and then estimated CO
2 emissions from point sources using the CSF method.
Cusworth et al. (2023) detected plumes worldwide from OCO-3 and estimated emissions from the world’s largest coal-fired power plants using the IME method. The Divergence (Div) method has demonstrated potential for estimating CO
2 point source emissions based on Copernicus Anthropogenic Carbon Dioxide Monitoring (CO2M) observations in the future (
Hakkarainen et al., 2022).
For city-scale CO
2 emission inversion,
Nehrkorn et al. (2010) integrated the Weather Research and Forecasting (WRF) model with the Stochastic Time-Inverted Lagrangian Transport (STILT) model, a Lagrangian Particle Dispersion Model and demonstrated the model’s capability in capturing regional CO
2 fluxes.
Wu et al. (2018) developed a modified version of the STILT model, “X-STILT”, which was used to estimate per capita CO
2 emissions from OCO-2 XCO
2 for 20 cities across multiple continents (
Wu et al., 2020).
Zheng et al. (2020a) used the light cross-sectional flux (LCSF) method focusing on the XCO
2 enhancements detected in satellite images, capturing emissions from 46 cities and industrial areas in China.
Chevallier et al. (2020) extended the approach of
Zheng et al. (2020a) to the global context and investigated the changes in global CO
2 emissions during the Coronavirus disease 2019 (COVID-19) pandemic.
Ye et al. (2020) optimized the ability of FFCO
2 emission estimation based on OCO-2 satellite observations in urban areas by combining a high-resolution chemical transport model (CTM) and multiple OCO-2 tracks.
Lei et al. (2021) evaluated the ability of three emission estimation methods (WRF-Chem, X-STILT, CSF method) to detect trends in urban fossil fuel CO
2 emissions and found that Open-Data Inventory for Anthropogenic Carbon dioxide (ODIAC) may have underestimated Lahore’s FFCO
2 emissions.
2.4 Challenges in monitoring large-scale emissions
Although CO
2 satellites hold promise in quantifying point sources and urban CO
2 emissions, substantial challenges persist in their regional and global scale application due to the satellite’s capability limitations, such as narrow swath width, low revisit frequency, and cloud interference (
Zheng et al., 2019;
Santaren et al., 2021).
The limited sampling coverage hinders the identification of CO
2 enhancements over broader coverage.
Chevallier et al. (2022) merged a 7-year dataset from OCO-2 with data from the OCO-3 satellite. Their findings only cover about one-third of global CO
2 emissions. There is a critical need for the development of CO
2 detection satellites with wider swath coverage and higher data density to enhance global monitoring capabilities. Additionally, current satellites only capture CO
2 concentrations during their overpass times, which may not necessarily align with actual XCO
2 enhancements considering daily variation. This discrepancy complicates CO
2 emission inversion (
Ye et al., 2020). Moreover, the emission inversion method also has potential limitations. Data-driven methods typically assume stable atmospheric conditions and relatively flat terrain, making them more suitable for point sources or urban scales, while the expensive computational costs associated with the model-driven methods may limit their application.
The insensitivity of CO
2 concentration to anthropogenic emission variations and the cloud’s impact on CO
2 sensing coverage further complicates the CO
2-based FFCO
2 emission inversion.
Weir et al. (2021) pointed out that the OCO-2 satellite has limited capacity to distinguish between changes in anthropogenic emissions and biosphere changes. The minor CO
2 fluctuations resulting from human activities (several ppm) are notably smaller than the background field (~400 ppm), typically falling within the margin of satellites’ observation error (
Eldering et al., 2017). During the COVID-19 pandemic, a valuable opportunity to investigate the response of atmospheric CO
2 concentrations to human activity changes, the decline in CO
2 observed from satellites was not as pronounced (
Buchwitz et al., 2021) compared to the decline in CO
2 emissions (
Liu et al., 2022).
Chevallier et al. (2020) also illustrated the difficulty in identifying large-scale decreases in CO
2 concentrations due to FFCO
2 emissions during the COVID-19 pandemic. Some key regions such as Hubei Province in China and the Bay of Bengal in India are shrouded in clouds, hindering effective CO
2 concentration observation and emission monitoring. Over sunny areas, yearly variations in cloud cover can still obscure the detection of CO
2 emission changes (
Chevallier et al., 2020).
3 Satellite NO2-based CO2 emission monitoring
3.1 NO2 satellites and NOx emission inversion
NO
2 satellites have undergone rapid and substantial advancements in recent decades (Tab.3). The Global Ozone Monitoring Experiment (GOME), launched in 1995, was a pioneer in global NO
2 observations (
Burrows et al., 1999). Successive satellites, including GOME-2A, GOME-2B, and GOME-2C, have since been launched, toward enhanced spatial resolution, advanced data processing algorithms, and broader spatial coverage (
Hassinen et al., 2016). The SCIAMACHY, launched in 2002 and initially designed for ozone monitoring, has also been utilized to detect NO
2 columns. In 2014, NASA and the Royal Netherlands Meteorological Institute (KNMI) collaborated to launch the Ozone Monitoring Instrument (OMI), which improved spatial resolution to 13 km × 24 km (
Levelt et al., 2006). This technological progress considerably enhanced the capability to resolve NO
2 distribution at fine scales like cities. The launch of TROPOMI in 2017 further advanced this ability, offering a resolution of 7 km × 3.5 km (5.5 km × 3.5 km since 2019) (
Veefkind et al., 2012). In May 2018, China launched the Environmental Trace Gases Monitoring Instrument (EMI) onboard the GaoFen-5 satellite.
Zhang et al. (2020) optimized the retrieval algorithm, demonstrating that EMI NO
2 vertical column densities (VCDs) perform comparably to those of OMI and TROPOMI. This satellite has also shown potential in quantifying point source CH
4 emissions (
He et al., 2024b). These satellites operate in sun-synchronous polar orbits, crossing the equator at approximately the same local times each day to deliver regular imagery.
Geostationary satellites, positioned above the equator, offer extensive temporal coverage, which enables monitoring of daily atmospheric concentration variations. Geostationary satellites continuously observe the same region of the Earth’s surface, covering approximately one-third of it. For instance, the Geostationary Environment Monitoring Spectrometer (GEMS) (
Kim et al., 2020b), launched by the Republic of Korea in 2020, has been monitoring atmospheric compositions over East Asia. NASA’s Tropospheric Emissions: Monitoring of Pollution (TEMPO), launched in 2023, has been providing hourly observations of atmospheric concentrations in North America (
Zoogman et al., 2017). In the near future, the ESA plans to launch Sentinel-4 to monitor air quality across Europe (
Ingmann et al., 2012).
The progress in remote sensing of NO
2 tends to have outpaced that of CO
2, featuring increased satellite revisit frequencies, finer spatial resolutions, broader coverage, and improved signal-to-noise ratios in column retrievals (
Cooper et al., 2022;
MacDonald et al., 2023). These advantages render satellite NO
2 retrievals highly effective for constraining nitrogen oxides (NO
x, sum of NO and NO
2) emission inversions. Research and methodologies (Fig.1) in this field have proliferated over the past two decades. The Mass Balance (MB) method, proposed by
Martin et al. (2003), employs a three-dimensional CTM to establish a linear localized response function between NO
2 TVCDs and surface NO
x emissions by grid. It then infers posterior NO
x emissions based on satellite observations and prior emissions. This method simplifies atmospheric transport and nonlinear chemistry to some extent and has proven feasible for inferring daily emission changes on a regional scale (
Zhao and Wang, 2009;
Gu et al., 2014). However, it is not appropriate for high spatial resolution due to the inter-grid transport and nonlinear chemistry of NO
2 in the atmosphere. The four-dimensional variational (4D-Var) (
Stavrakou et al., 2008) and the Ensemble Kalman Filter (EnKF) method (
Miyazaki et al., 2017) utilize data assimilation techniques to minimize the gap between observations and model simulations by adjusting the CTM initial conditions and emissions input. These methods account for atmospheric transport and nonlinear chemistry but require multiple iterations of simulations with high computational costs.
To avoid the limitations associated with CTM, researchers have proposed several simplified models to expedite NO
x emission estimation, such as the Exponential Modified Gaussian (EMG), Div, and Peking University High-resolution Lifetime-Emission-Transport (PHLET) methods. These approaches integrate satellite-monitored NO
2 column concentrations with wind fields to update NO
x emissions at high spatial resolutions on a monthly or annual basis. In the EMG method, NO
x emissions are estimated by integrating wind speed with NO
2 line density based on different wind directions, assuming that downwind NO
2 TVCDs from a point source follow a Gaussian distribution. This method has been widely used to quantify emissions from isolated point sources like power plants (
Beirle et al., 2011;
Liu et al., 2016;
Goldberg et al., 2019b). The Div method defines the divergence of flux (
F) as the difference between emissions (
E) and sinks (
S) (▿
F =
E −
S), creating a vector field of flux by multiplying the wind field with NO
2 columns. The finite difference method is then used to calculate the divergence of this flux field, with peak fitting employed to infer annual emissions for each source from the divergence field (
Beirle et al., 2019;
de Foy and Schauer, 2022;
Goldberg et al., 2022). The PHLET method estimates NO
x emissions at a high resolution of 0.05° × 0.05° using the adjoint method. It assumes that the NO
2 TVCDs remain in a steady-state for a period, considering the influence of local nonlinear chemistry on transport and lifetime (
Kong et al., 2019).
3.2 Co-emitted CO2 and NOx from combustion sources
During fossil fuel combustion, CO
2 and NO
x exhibit co-emission characteristics temporally and spatially (Fig.2). CO
2 is generated from the combustion of fossil fuels, while NO
x arises from the oxidation of nitrogen within the fuel and the decomposition of atmospheric N
2 at elevated temperatures including thermal, prompt, and fuel NO
x (
Correa, 1993;
Graven et al., 2009). This congruence in emission processes inspires a concept and presents an opportunity to advance the simultaneous control of greenhouse gases and air pollutants, thereby addressing the dual challenge of climate change mitigation and air quality improvement (
Shi et al., 2022). The co-inversion and co-monitoring of NO
x and FFCO
2 emissions are not only feasible but also imperative.
The CO
2-to-NO
x emission ratios (ERs) vary substantially across different source sectors. Thermal power generation, characterized by high coal consumption and the implementation of flue gas denitrification technologies, tends to have the highest ERs. The transport sector, predominantly reliant on oil consumption, has the lowest ERs (
Zheng et al., 2018;
Zhao et al., 2024). These ERs are influenced by factors including fuel characteristics, oxidation ratio, operating condition, and post-treatment device (
Ammoura et al., 2014;
Silva and Arellano, 2017;
Bares et al., 2018;
Fu et al., 2022). CO
2 emission factors differ markedly across fuel types, declining from coal (97.3–96.8 gCO
2/MJ) to oil (59.1–73.9 gCO
2/MJ), and natural gas (50 gCO
2/MJ) (
IPCC, 2006). NO
x emission factors from combustion processes are more influenced by temperature and post-treatment technologies, which vary by source. For instance, in thermal power plants, NO
x emissions increase with temperature due to the contribution of thermal NO
x (
Thompson et al., 1972;
Abel et al., 2017). In the transport sector, diesel-powered vehicles emit more NO
x at lower ambient temperatures due to the decreased efficiency of post-combustion NO
x control technologies in cold conditions (
Grange et al., 2019). Higher CO
2-to-NO
x ERs generally indicate stricter emission controls and improved oxidation ratios (
Lei et al., 2022), and they can serve as indicators of national economic development (
Miyazaki and Bowman, 2023).
The complexities associated with CO
2-to-NO
x ERs pose challenges in data acquisition. Direct
in situ measurement of CO
2 and NO
x emissions across all processes is impractical due to the time-consuming and labor-intensive nature of the task. Recent advancements in remote sensing satellites have enabled the use of observed CO
2 and NO
2 concentrations (mostly spaceborne monitor) as a novel approach to characterize oxidation ratios or emission properties for point sources and broader scales (
Reuter et al., 2014;
Silva and Arellano, 2017;
Park et al., 2021).
Ammoura et al. (2014) employed
in situ monitored CO
2 and NO
2 concentration ratios to characterize their co-emission from traffic in a tunnel study.
Lindenmaier et al. (2014) utilized multiscale CO
2 and NO
2 concentrations to differentiate emission factors from two coal-fired power plants with varying scrubbing technologies.
MacDonald et al. (2023) analyzed the enhancement ratios of OCO-2 CO
2 and TROPOMI NO
2 concentrations in 27 large urban areas worldwide, utilizing this data to identify potential biases contributing to discrepancies between observations and inventories.
3.3 Framework inferring CO2 emissions from NO2 retrievals
The co-emission characteristics of FFCO
2 and NO
x has sparked the growing research interest in introducing NO
2 as a proxy to infer FFCO
2 emissions with several advantages. NO
2 anomalies induced by emission sources are more distinguishable from the background compared to CO
2 anomalies, showing magnitudes higher (
Hakkarainen et al., 2016;
Liu et al., 2024). The short atmospheric lifetime of NO
2, lasting several hours, renders its concentration highly responsive to emission sources and their changes (
Lange et al., 2022;
Li et al., 2023a). Conversely, CO
2 possesses a much longer lifetime of hundreds of years (
Eby et al., 2009), resulting in a high background concentration (hundreds of ppm) and only a marginal increase of several ppm due to proximate emission sources (
Nassar et al., 2017;
Reuter et al., 2019). Consequently, the response of CO
2 concentration to local and regional anthropogenic emissions is insensitive (
Chevallier et al., 2020). Furthermore, NO
2 predominantly originates from anthropogenic sources, whereas CO
2 stems from both anthropogenic and biogenic sources (
Ciais et al., 2014;
Zheng et al., 2023). This makes NO
2 a valuable proxy for disentangling anthropogenic CO
2 emissions. Study indicates that approximately 92% of CO
2 emissions can be detected through satellite NO
2 plumes, highlighting the potential of NO
2 as an effective proxy for FFCO
2 emissions (
Finch et al., 2022).
There are two primary methodologies for integrating satellite NO
2 data into CO
2 emission estimation (Fig.3). For the first method, researchers utilize NO
2 observations as locators, constraints, or validators for deciphering CO
2 emission plumes. Satellite NO
2 column observations serve to pinpoint fossil fuel emission sources and provide constraints on plume shapes for CO
2 plume analysis (
Fujinawa et al., 2021;
Fuentes Andrade et al., 2024). By initially detecting and determining the plume shape with easily monitored NO
2, this approach avoids the need for extensive scanning of the CO
2 swath, which has a low signal-to-noise ratio, thereby enhancing computational accuracy and efficiency (
Reuter et al., 2019). Likewise, NO
2 observations can serve as validators for evaluating CO
2 enhancements resulting from nearby sources (
Kiel et al., 2021).
For the second method, studies employ the ERs-based inversion approaches, wherein NO
2 observation-constrained NO
x emissions are first estimated. These NO
x emissions are subsequently converted to CO
2 emissions using CO
2-to-NO
x ERs. A crucial component in the first step involves establishing the relationship between NO
2 concentration and NO
x emissions, which can be achieved through various approaches (detailed discussion in Section 3.1). For instance,
Berezin et al. (2013) utilized the trend of GOME and SCIAMACHY NO
2 column data to derive long-term NO
x emission trends from 1996 to 2008 based on an assumed linear relationship.
Zheng et al. (2020b) calculated the gird-by-grid responses of NO
2 concentration changes to NO
x emission variations by introducing a 40% emission perturbation in GEOS-Chem simulations.
Zhang et al. (2023) employed a superposition column model to fit the TROPOMI NO
2 line density over a domain to estimate daily NO
x emissions.
Obtaining CO
2-to-NO
x ERs is relatively challenging, with most studies directly adopting regional emission ratios from existing inventories, such as the Regional Emission Inventory in Asia (REAS) (
Berezin et al., 2013), the continuous emissions monitoring system (CEMS) (
Liu et al., 2020), and the Air Benefit and Attainment and Cost Assessment System Emission Inventory (ABACAS) (
Zhang et al., 2023). Efforts are underway to improve these ratios. Some studies have attempted to use observed CO
2-to-NO
2 concentration ratios as reported ERs constraints, assuming that observed atmospheric enhancement ratios partly reflect the ERs (
Reuter et al., 2014;
Wren et al., 2023). Additionally, Zheng et al. (
2020b) dynamically updated sectoral ERs at a daily and grid-scale based on the Multi-resolution Emission Inventory for China (MEIC), considering the spatiotemporal variations in sectoral contribution to emissions.
3.4 Application spanning from point source to country
The integration of NO
2 into CO
2 emission inversion methods has been widely applied across various spatial scales, ranging from point sources and urban regions to national levels. At the point sources level, researchers typically employ plume-based methods to estimate CO
2 emissions from power and industrial facilities (
Konovalov et al., 2016;
Fujinawa et al., 2021). Notably,
Reuter et al. (2019) utilized co-located NO
2 observations from TROPOMI to constrain the shape of CO
2 plumes observed by OCO-2, enabling the calculation of cross-sectional fluxes and CO
2 emissions for 6 hotspot areas.
Liu et al. (2020) estimated CO
2 emissions from 21 individual power plants by integrating EMG-updated NO
x emissions and CO
2-to-NO
x ERs. Additionally,
Hakkarainen et al. (2023) applied the CSF method to estimate CO
2 emissions from 6 power plants based on OCO-3 CO
2 plumes, leveraging TROPOMI NO
2 data for shape constraints. It is noteworthy that with the rapid development of computing power, deep learning algorithms, such as CNN, are being used to calculate local CO
2 emissions from satellite images (
Dumont Le Brazidec et al., 2024).
At the city scale, researchers begin by estimating NO
x emissions using the plume method before converting them to CO
2 emissions based on the reported ERs.
Goldberg et al. (2019a) utilized the EMG method to update NO
x emissions for 8 cities using OMI observations, subsequently converting them to CO
2 emissions with CO
2-to-NO
x ERs from inventories provided by the US Environmental Protection Agency.
Zhang et al. (2023) improved the superposition column model to estimate daily NO
x emissions in Wuhan, utilizing reported CO
2-to-NO
x ERs for CO
2 emission estimation. The MB method applies to city-scale emission inversion, as well.
Yang et al. (2023) established empirical relationships between OCO-3 CO
2 SAM and TROPOMI NO
2 in three urban cities, applying these relationships to NO
2 fields to infer NO
2-derived CO
2 fields for CO
2 emission estimation using a simplified mass balance approach.
At the national level, the ERs-based inversion approaches are prevalent.
Berezin et al. (2013) quantified multiannual changes in FFCO
2 emissions from China using GOME and SCIAMACHY NO
2 retrievals in conjunction with ERs from inventories. In an early application of data assimilation,
Konovalov et al. (2016) optimized NO
x emissions in Europe by minimizing the bias between modeled and observed NO
2 column concentrations, subsequently converting them to CO
2 emissions based on reported ERs.
Zheng et al. (2020b) and
Li et al. (2023b) integrated satellite constraints and inventory sectoral profiles to infer sectoral NO
x emissions using the MB method, further translating them into sectoral CO
2 emissions in China. Miyazaki and Bowman (2023) employed a state augmentation technique (a kind of data assimilation), to infer NO
x emissions first, then converted them to CO
2 emissions over multiple countries using reported ERs.
Feng et al. (2024) utilized a data assimilation system (RAPAS) to deduce NO
x emissions and thereafter projected CO
2 emissions, with surface observation of NO
2 concentrations.
4 Structural and data uncertainties
4.1 Structural uncertainties and impacts
Considerable progress has been made in estimating CO
2 emissions using NO
2 as a proxy, yet several structural uncertainties remain (Fig.4). Establishing the relationship between atmospheric species concentration and emissions is fundamental to all inversion frameworks, which also represents the primary source of uncertainty. Current methodologies inherently involve certain assumptions. For example, plume-based methods presuppose that the shapes and spatial coverage of CO
2 and NO
2 plumes are similar or even identical (
Reuter et al., 2019). However, this assumption does not always hold due to their different chemical reactivity and diffusion characteristics (
Massman, 1998;
Squadrito and Pryor, 1998). Another assumption is that the distribution of concentrations (e.g., plume) reflects the magnitude of emissions nearby as an aggregated source. This assumption is highly sensitive to wind fields and the spatial coverage of the plume (
Kim et al., 2020a;
Kurchaba et al., 2024), while not applicable at larger scales like region or country. In the case of the ERs-based inversion approaches, the uncertainty sources vary according to their underlying principles. For instance, CTM-associated inversions are subject to systematic errors intrinsic to the physical and chemical mechanism inside CTM (
Shah et al., 2023), while MB methods are affected by the ‘smearing effect’ due to horizontal transport of the species (
Palmer et al., 2003). Although eliminating biases in the relationship between concentration and emissions is rather challenging, researchers have employed various strategies to approximate this relationship, such as using relative changes to reduce the systematic errors’ impact.
Satellite remote sensing, typically aboard polar-orbiting satellites and used as observational constraints, samples only at specific times of the day. For example, TROPOMI provides global daily NO
2 column concentrations with a local time of approximately 13:30 (
van Geffen et al., 2022), while OCO-3 monitors CO
2 at varying local overpass times (
Taylor et al., 2023). Current studies utilizing the ERs-based inversion approaches generally rely on concentrations at a specific moment or averaged over a period to represent average concentrations. However, NO
2 and CO
2 column concentrations typically decrease around noon and peak at night (
Olsen and Randerson, 2004;
Li et al., 2021), while the FFCO
2 emissions peak around noon but decline at night (
Nassar et al., 2013). This asynchrony in their diurnal variation weakens the temporal representativeness of remote sensing data. To tackle this issue, some studies combine hourly-monitored concentrations with CTM simulations to restore diurnal changes in remote sensing data (
Boersma et al., 2009). Another challenge arises when jointly using multiple satellites, such as OCO-3 and TROPOMI, where the research assumes that the captured plumes are temporally consistent. However, the sampling times of different satellites often do not align. Consequently, this assumption leads to data incompatibility issues, particularly in the plume-based methods (
Hakkarainen et al., 2023). To mitigate this time incompatibility, researchers try to limit the time difference within a certain range when incorporating different satellite data for usage.
Distinguishing anthropogenic emissions from biogenic sources and attributing them to specific sectors remains challenging given the mixed nature of all sourced NO
2 or CO
2 concentrations in the atmosphere. The plume-based method is susceptible to natural flux fluctuations (
Attermeyer et al., 2021), while the sectoral attribution of anthropogenic emissions through the ERs-based inversion approaches is challenging (
Qu et al., 2022). Although NO
2 aids in disentangling anthropogenic emissions, the impact of biogenic sources, such as soil emissions, becomes noticeable in summer (
Lu et al., 2021). Currently, the plume-based methods are mainly conducted near known emission sources, and studies with the ERs-based inversion approaches usually exclude natural backgrounds by setting a concentration threshold, or exclusively perturbing anthropogenic emissions while fixing biogenic sources. However, neither approach can reliably estimate FFCO
2 emissions in any industrialized region. Future advancements in finer and more accurate emission inventory compilation and satellite observation technology are expected to enhance source identification and attribution.
4.2 Data uncertainties and management
The uncertainty of observational constraints, particularly those from remote sensing, ranks foremost among data uncertainties in inversion methods. Satellite monitoring translates spectral signals into concentrations, a process sensitive to retrieval parameters including viewing geometry, cloud fraction, aerosol presence, surface albedo, and snow content (
Lin et al., 2014;
Lamsal et al., 2021). Additionally, the spatiotemporal coverage provided by satellites is often insufficient for emission inversion on a daily scale or fine spatial resolution (e.g., kilometer level). These retrieval errors, especially random ones, and poor coverage can be mitigated by using averaged data over extended periods. For instance,
Beirle et al. (2023) targeted point sources by calculating the divergence of NO
2 flux averaged over three years. AI approaches have also proven effective in filling in missing satellite retrieval data (
Wei et al., 2022).
The CO
2-to-NO
x ERs are key parameters in the ERs-based inversion approaches, which are primarily derived from available emission inventories. Therefore, the accuracy and timeliness of converted CO
2 emissions are highly contingent on the inventory, whose emission factors are mainly derived from widely adopted references (e.g., IPCC guidelines) or empirical extrapolation.
Konovalov et al. (2016) demonstrated that different ERs settings could result in emission biases of up to 24.9%, comparable to or even larger than the uncertainty of emission inventory (
Solazzo et al., 2021). Currently, comprehensive industry-wide emission factors are still lacking, despite efforts to measure in situ sector-specific emission factors (
Park et al., 2011;
Shen et al., 2014). Moreover, ERs can vary over time due to changes in emission controls, energy structures, and sectoral contributions (
Luderer et al., 2022;
Li et al., 2023c). As obtaining comprehensive in situ measured data are impractical, researchers continue to optimize CO
2-to-NO
x ERs given their critical role in CO
2 conversion. For example, simultaneous monitoring of CO
2 and NO
2 concentrations has been used as constraints of ERs (
Wren et al., 2023). Future advancements in satellite technology capable of simultaneously monitoring NO
2 and CO
2 will further strengthen these constraints.
Most inversion methods rely on a specific prior, usually a bottom-up inventory, serving as the base for posterior emission updates. However, compiling bottom-up inventories introduces uncertainties stemming from activity data and emission factors (
Solazzo et al., 2021). The mapping of emission inventory with spatial proxies introduces additional uncertainties in emissions due to their spatial mismatch (
Zheng et al., 2017). These uncertainties can propagate into the NO
2-based CO
2 emission inversion process. In the plume-based methods, the prior acts as an anchor for emission attribution, implying that uncertainty in the emission location within the prior can lead to misinterpretation of emission estimates. In contrast, in the ERs-based inversion approaches, the prior establishes a baseline that influences the posterior estimates by introducing inaccurate ERs (
Solazzo et al., 2021). Data assimilation requires grid-level uncertainty (
Choulga et al., 2021), yet most inventories only report uncertainty ranges for total emission estimates. Efforts have been made to optimize priors and develop more accurate inventories. For instance,
Cooper et al. (2017) and
Li et al. (2019) iteratively updated priors for subsequent posteriors in conventional MB methods, resulting in notable improvements. Researchers are also working on establishing plant-specific emission inventories (
Lei et al., 2023;
Xu et al., 2023). With the development of multisource data fusion, developing accurate and timely updated inventories is becoming increasingly feasible, thereby reducing the uncertainties associated with their use as priors.
5 Emerging trends and future avenues
5.1 Next-generation satellites
The aforementioned uncertainties have proposed several requirements for the next-generation satellites (Fig.5). These satellites are expected to possess several improved capabilities, including simultaneous monitoring of CO
2 and NO
2 concentration, hourly remote sensing, high signal-to-noise ratios and spatial resolution, and wide swath coverage. The co-monitoring of multiple species addresses spatiotemporal discrepancies in joint satellite usage, offering better mutual verification and constraints for ERs (
Bernd et al., 2021). Hourly observations, such as those from geostationary satellites, capture diurnal profiles, addressing the insufficient representativeness of specific time sampling (
Penn and Holloway, 2020). High signal-to-noise ratios and spatial resolution enable finer mapping of concentrations, facilitating the identification of emission-induced enhancements from background levels (
Broquet et al., 2018). Wide swath capabilities extend spatial coverage to encompass both emission sources and background areas, allowing for the capture of entire plumes (
Pillai et al., 2016).
While achieving high spatial and temporal resolutions simultaneously with satellites seems unlikely at present, several upcoming satellites are equipped with some of these characteristics. For instance, the ESA and the European Commission have committed to launching a CO
2 satellite imagery, the CO2M, in 2026. This satellite will monitor atmospheric CO
2 and NO
2 concentrations at a high spatial resolution (4 km × 4 km) along a 250 km wide swath and achieve global coverage within 5 days (
Bernd et al., 2019). GHGSat, a company specializing in commercial greenhouse gas satellites (
Jervis et al., 2021), launched GHGSat-C10 in 2023, which is capable of detecting individual CO
2 sources with a resolution of up to 25 m, such as refineries, steel mills, and cement plants. Likewise, advancements in NO
2 satellite technology are promising, particularly with geostationary satellites. TEMPO, launched in 2023, can provide hourly observational data for North America at a resolution of 2 km × 4.5 km (
Zoogman et al., 2017). Additionally, the ESA plans to launch Sentinel-4 in 2024 to monitor hourly air quality across Europe (
Ingmann et al., 2012). Hourly monitoring of NO
2 column concentrations enables the inversion estimate of hourly NO
x emissions, representing more accurately diurnal fluctuations in NO
2 concentrations and lifetimes. This enhanced monitoring capability improves the establishment of the concentration-emission relationship, thereby reducing uncertainties in emission inversion (
Xu et al., 2024). Nevertheless, the global application of geostationary satellites is hindered by their limited coverage of specific regions. The deployment of these satellites is expected to reduce existing uncertainties in NO
2-based CO
2 emission inversion systems.
The satellites mentioned above, as well as those currently available, are passive remote sensors that measure the sunlight spectrum. Besides, active sensors first illuminate the target and then capture the reflected light. This self-emission characteristic enables them to provide high-precision data even under complex conditions. Presently, active remote sensors have been solely employed in aerosol detection (e.g., CALIPSO satellite) within atmospheric composition studies (
Winker et al., 2009). Future advancements are anticipated to extend the application of active remote sensing satellites to monitor additional atmospheric components, such as NO
2 and CO
2.
Advancements have also been made in satellite constellations (
Ulybyshev, 2008), which refers to a network of multiple satellites working together in space with specific configurations and orbital arrangements to achieve global coverage, real-time monitoring, and high-resolution observation of the Earth’s surface (
Singh et al., 2020). Currently, constellations support research and applications in numerous fields, including global navigation, telecommunications, environmental monitoring, and natural disaster response (
Maral et al., 1991;
Hao et al., 2018;
del Portillo et al., 2019). Several studies have investigated the optimal number of satellites required for effective global coverage in a constellation, specifically examining whether three or four satellites suffice (
Zhu and Gao, 2017). Additionally, research has explored the potential of satellite constellations for CO
2 emission inversion (
Lespinas et al., 2020). As satellite constellations continue to mature with wider adoption, their benefits of extensive coverage, near-real-time data retrieval, and high-resolution capabilities are expected to drive further advancements in the field of CO
2 emission inversion.
5.2 Development of new inversion methods
In addition to refining observational techniques, the development of novel inversion methods is imperative, albeit challenging. Key avenues for advancement are elucidated as follows:
1) Establishing a concentration-emission relationship at a granular level (e.g., hundred meters, point source) in tandem with the advancement of satellite technologies. Finer resolution approximates the localized relationship between concentration and emissions, aiding in the attribution to specific sources and thus targeting those with substantial potential for emission reduction. For example, incorporating inter-grid transport and nonlinear responses in the MB method can overcome resolution limitations, facilitating subsequent emission attribution (
Turner et al., 2012). High resolution also benefits the capture of plumes downwind from sources and the extraction of emission magnitude therein (
He et al., 2024a). Another inquiry arises: How can we precisely assess the reliability and accuracy of estimated emissions? As high-resolution emission inversion techniques continue to advance, ensuring the reliability and accuracy of estimated emissions becomes increasingly crucial. One effective approach is to refine plant-level emission monitoring (regarded as a “truth”) to verify inversion results (
Tang et al., 2023). This method also provides foundational data for delineating emission reduction responsibilities and evaluating the effectiveness of emissions reduction efforts.
2) Amalgamating diverse methodologies that harbor complementary strengths to achieve spatiotemporal refinement. Current techniques often require a trade-off between refinement across different dimensions, typically sacrificing one to enhance another. For instance, achieving high spatial resolution (kilometer scale) often necessitates a reduction in temporal resolution (monthly or annual average). These methods have independently advanced various dimensions, such as finer temporal resolution (days or even hours), spatial resolution (kilometers), and detailed sectoral attribution. The integration of strengths from various approaches can enable simultaneous multi-dimensional refinement.
Moon et al. (2024) combined IFDMB and 4D-Var as a hybrid inversion system to estimate NO
x emissions in South Korea, wherein IFDMB constraining spatial distribution and 4D-Var estimating temporal variations on an hourly scale.
Tibrewal et al. (2024) integrated an ensemble of bottom-up inventories and inversions from the Global Carbon Project to resolve CH
4 emissions from coal, oil, and natural gas exploitation for major producing nations. The synergistic use of multiple tools can overcome the limitations of individual methods, enabling simultaneous advancements across multiple dimensions.
3) Integrating emerging analytical paradigms, Data Science and AI, to enhance the efficiency and accuracy of processing extensive multi-source datasets. With the surge in satellite retrievals and the growing computational burden of model simulations, coupled with the demand for near-real-time inversion, conventional data processing methods are becoming inadequate. In this context, leveraging new analytical methods to improve efficiency is becoming imperative (
Vance et al., 2024). Data Science plays a pivotal role in processing big data by applying advanced statistical techniques and algorithms (
zu Castell et al., 2022). For instance, data fusion combines information from multiple sources, enhancing the accuracy and completeness of datasets. This method offers special advantages in NO
2-based CO
2 emission inversion due to the multi-source datasets involved, including CO
2 and NO
2 satellites, emission inventories, and model simulations. Besides, AI substantially improves efficiency and automation capabilities. Using AI to train on a subset of data allows for the extraction of effective information and internal connections within big data, which can then be quickly applied to newly generated data. A sophisticated AI can supersede manual plume identification and the modeling of concentration-emission relationships.
Finch et al. (2022) utilized a CNN to automatically identify NO
2 plumes observed by TROPOMI.
Joyce et al. (2023) implemented a deep neural network to effectively identify and quantify CH
4 emissions from point sources using satellite imagery at 30 m resolution.
Dumont Le Brazidec et al. (2024) demonstrated that the CNN can automatically learn to identify the XCO
2 plume and assess emissions based on the plume concentrations. Overall, the adaptability, learning, and scalability of AI techniques demonstrate strong potential for reliable and efficient big-data-based emission inversion, such as NO
2-based CO
2 emission inversion involving multiple species.
5.3 Data stewardship and sharing with caution
Despite extensive efforts and opportunities to enhance the accuracy of inferred emissions, clear and detailed caveats must accompany the released datasets to users. The compilation of a rigorous dataset for universal use necessitates publishers to include a statement of uncertainty and scope of application, similar to the IPCC report or the Product User Manual of satellite data. This attachment should provide comprehensive information to enable users to accurately interpret and utilize the dataset, thus mitigating the risk of misusage. The principles in data stewardship and sharing include transparency and clarity, delineation of application scope, uncertainty statements, provision of examples and usage recommendations, and ongoing updates and communication.
Firstly, clear and understandable documentation should be published to elucidate the data’s sources, processing methodologies, and potential influencing factors. This ensures users grasp the background and limitations of the data. Secondly, defining the scope of application for the data and elucidating it in the documentation aids users in understanding their specific utility and limitations. Thirdly, accompanying data releases with uncertainty statements is crucial, detailing the data’s accuracy and reliability. This may involve confidence intervals, potential biases or error ranges, and uncertainties associated with the data generation process. Incorporating attributes quantifying the data quality, such as the quality_flag in TROPOMI NO
2 retrievals (
van Geffen et al., 2020), enhances the practicability. Fourthly, offering examples and usage recommendations assists users in correctly interpreting and utilizing the data. This could include data visualization examples, data processing advice, practices tailored to specific domains, and processing scripts. Finally, ensuring regular updates to data products and maintaining open communication channels with users to promptly address any issues or feedback fosters trust and enhances user satisfaction with the data products.
Some data stewardship and sharing have met the aforementioned standards. For instance, when KNMI released the TROPOMI L2 NO2 satellite data product, they included comprehensive documentation such as data definitions and product user manuals. Additionally, several bottom-up inventory datasets are continuously updated post-release, maintaining user communication and improving emission data. Notable examples include the MEIC, Carbon Monitor, Emissions Database for Global Atmospheric Research (EDGAR), and the Community Emissions Data System (CEDS). Due to their later development, inversion data currently lacks a mature framework for data stewardship and sharing. However, these frameworks can draw valuable insights from the established practices of emission inventory data.
6 Conclusions
To facilitate the attainment of NDCs, the GST has been implemented as a “ratchet mechanism” every five years since 2023, emphasizing the necessity for establishing anthropogenic CO2 MVS capacity, particularly FFCO2. Research endeavors aimed at compiling more accurate, detailed, and timely FFCO2 emission inventories are flourishing. Despite progress in direct CO2 emission inversion based on CO2 observations, challenges persist due to the long lifetime of atmospheric CO2, which obscures concentration enhancements above the background. CO2 emission inversion solely reliant on CO2 concentration does not work for larger spatial scales, such as regional and national levels. In contrast, the short lifespan of co-emitted NOx, with enhancements substantially higher than background levels, presents a more viable proxy. Therefore, integrating NO2 observation has emerged as a promising strategy for CO2 emission inversion. Two prevalent methodological systems include: 1) using NO2 observations to locate, constrain, and validate CO2 plumes from nearby sources, particularly for point sources (plume-based method); and 2) inferring NOx emissions from NO2 observations, followed by their conversion to CO2 emissions using CO2-to-NOx ERs (ERs-based inversion approach), which is often applied at the national scale.
Despite the advantages of using NO2 as a proxy for inferring FFCO2 emissions, inherent uncertainties persist. Structural uncertainties include 1) the establishment of relationships between concentrations and emissions at fine scales; 2) the representativeness of satellite retrievals sampled at specific times; and 3) the disentanglement of anthropogenic emissions from biogenic sources. Additionally, data uncertainties stem from satellite retrievals, CO2-to-NOx ERs, and priors. To date, numerous studies have been dedicated to mitigating these challenges, yet they remain formidable. Future efforts hold substantial promise, including the deployment of next-generation satellites and the development of advanced inversion systems. Achieving these objectives will necessitate collaboration among diverse research communities, including remote sensing, emission inventory, transport model improvement, and atmospheric inversion algorithm development. This interdisciplinary cooperation is vital for enhancing the accuracy and reliability of CO2 emission estimates, supporting the GST ratchet mechanism, and ultimately contributing to the Paris Agreement’s objective of limiting global temperature rise to within 2 °C (ideally 1.5 °C).