Enhancing scientific transparency in national CO2 emissions reports … – Nature.com

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 13, Article number: 15427 (2023) Cite this article
Metrics details
Biennial Update Reports (BURs) are essential requirements from the United Nations Framework Convention on Climate Change (UNFCCC). However, many non-Annex I countries have not submitted these reports due to difficulties in compiling the inventories. We developed a satellite-based method for the top-down inverse estimation of CO₂ emissions using partial-column data in the lower troposphere obtained by the Greenhouse Gases Observing Satellite, adopted to validate the Mongolian 2^nd BUR (BUR2) for the energy sector in 2018. The estimated CO₂ emissions were only 1.5% higher than those reported in the BUR2; these were also very close (4.2% smaller) to estimates from the Emission Database for Global Atmospheric Research. Mongolia is the first country to introduce an independent inverse estimate in its BUR, thereby increasing scientific transparency. Our method could be applied into other countries and could be incorporated into UNFCCC reporting guidelines, significantly improving global CO₂ emission estimates.
Anthropogenic CO₂ is the most significant contributor to climate change. An increase in global temperature of more than 1.5–2 °C will occur during the twenty-first century unless a considerable reduction in CO₂ and other greenhouse gas (GHG) emissions occurs in the coming decades¹. National GHG inventories are an essential input to the Global Stocktake (GST), which was established in 2021 by the 26^th Conference of the Parties. The GST aims to assess the collective progress towards achieving the aims of the Paris Agreement and its long-term goals. All parties are required to submit a GHG inventory as part of the Biennial Transparency Report (BTR) under the Enhanced Transparency Framework^2,3. Countries under the United Nations Framework Convention on Climate Change (UNFCCC) calculate CO₂ emissions based on a bottom-up approach according to standard guidelines developed by the Intergovernmental Panel on Climate Change (IPCC)⁴. IPCC guidelines for National Greenhouse Gas inventories provide comprehensive information for each GHG emission and sink, including standard calculation methods along with parameters for standardized emission and absorption factors. Moreover, individual countries determine their own calculation methods, considering the actual domestic situation, data availability, and scientific knowledge of each emission and sink. However, such bottom-up approaches can lead to significant uncertainties in relation to missing social and economic information and inaccurate emission factors. In 2022 the National Academies of Sciences, Engineering, and Medicine published an assessment of the current capabilities of the inverse top-down, bottom-up, and hybrid approaches against the six evaluation pillars, including transparency, usability and timeliness, evaluation and validation, completeness, inclusivity, and communication⁵. However, important considerations regarding cost-effectiveness and less labour requirements were omitted.
Since satellite data and the inversion method are features of the top-down approach, cost-effective and less labour-intensive methods are advantages of the top-down approach compared with the bottom-up approach.
A comparison of the inverse top-down and bottom-up approaches indicated that transparency are difficult to assess. Nevertheless, the former has an overwhelming advantage because it is cost-effectiveness, less labour-intensive, and has favourable evaluation and validation⁵ features.
Approximately half of non-Annex I countries have not yet submitted Biennial Update Reports (BURs) to the UNFCCC due to knowledge and resource (financial and labour) gaps^6,7, and practical solutions for shrinking these gaps are necessary. Specifically, we believe it is possible to shrink these knowledge and resource gaps by using a cost-effective, less labour-intensive, transparent, top-down approach with a great capacity for evaluation and validation, and making it possible for all non-Annex I countries to submit BTRs to the UNFCCC.
In practice, hybrid approaches have a strong usability that optimizes the integration of bottom-up and top-down approaches to provide users with the best available comprehensive GHG emissions information⁵. The 2019 Refinement of the IPCC 2006 report describes an independent approach based on atmospheric observations and inverse analysis (top-down approach) and recommends that bottom-up emission inventories be scientifically validated⁸. A successful GST is based on a submission of validated GHG inventory from all non-Annex I countries; the critical challenge is to accelerate these independent validation efforts, close the knowledge and resource gaps so that all participating countries can submit reports, and improve the transparency of emission inventories in each country. In this direction, an ambitious project was carried out in Mongolia, which entailed the use of position satellite-based inverse analysis results in validating the BUR. Providing information on knowledge and resources in the BTR will be a significant challenge. Here, considering Mongolia as a case study, we propose an a posteriori method that uses satellite data to verify CO₂ emissions estimates contained in the BUR.
Mongolia plans to submit its 2^nd BUR (BUR2) with data through 2019. We propose including our estimates, obtained with a top-down approach, in BUR2. Mongolia has massive coal reserves⁹; thus, its primary energy source is coal¹⁰. In addition, Mongolia’s national CO₂ emissions for 2018¹¹ were the highest in the energy, contributing to 97.8% of the total, followed by those relating to industrial processes and product use with 2.2%. CO₂ emissions from the energy sector amount to 20.3 Tg. Furthermore, urban areas are responsible for approximately 70% of global fossil fuel-related emissions¹² and play an essential role in mitigation strategies under the Paris Agreement’s action plan. Consistently, emissions from the capital Ulaanbaatar account for approximately 69% of anthropogenic CO₂ emissions in Mongolia (Supplementary Table S1). Therefore, we define Ulaanbaatar as our target city for obtaining CO₂ emission estimates using the top-down approach.
Approximately 60% of CO₂ emissions in Ulaanbaatar originate from coal-fired power plants, 29% from coal stoves in “ger” (traditional transportable dwelling) districts, and 9% from automobiles¹³. Coal-fired power plants supply electricity and hot water to the urban area in the centre of the capital. Climate change has caused frequent mass livestock death events (called “dzud”) owing to summer droughts and heavy winter snowfalls (approximately 30% of all livestock deaths in Mongolia in 2010 were caused by dzud¹⁴). Dzud forced small herding families who lost their livestock to migrate with their gers and household goods around the capital city of Ulaanbaatar. Due to the Mongolian tradition of pastoralism, the government allowed people to occupy a specific portion of land anywhere in the country, resulting in the formation of sprawling ger settlements surrounding Ulaanbaatar (currently, migration is restricted by the government). In addition, the severe weather in winter, with temperature as low as − 40 °C, and the city’s location in a topographic depression that creates an atmospheric inversion layer, results in the trapping smoke from ger districts below the inversion layer, causing the worst air pollution in the world^15,16,17, especially in the winter season (Fig. 1).
Copyright: 2017 Reuters/B. Rentsendorj < https://static.reuters.com/resources/r/?m=02&d=20170207&t=2&i=1171734569&r=LYNXMPED160TP&w=1920 > (Thomson Reuters permitted the use of this photograph).
Typical winter scenery in Ulaanbaatar, Mongolia.
Our approach, targeting 2018, relies on a Bayesian framework and uses data from the Greenhouse Gases Observation Satellite (GOSAT) and a high spatial resolution regional atmospheric transport model^{18,19,20,21,22}. A high-resolution regional atmospheric transport model is capable of capturing fine-scale variability in XCO2 distributions caused by transport and emission processes at urban scales¹⁸. Second, using the lower-layer concentrations from the GOSAT-based 2-layer analysis product by the Japan Aerospace Exploration Agency/Earth Observation Research Center (JAXA/EORC), the approach allows the quantification of an increase in observed CO₂ concentrations within the city^23,24.
We introduced independent science-based techniques to validate GHG emissions inventory reporting and contribute to scientific transparency²⁵ in Mongolia. This approach could be applied to other countries and regions, integrated into the Paris Agreement’s enhanced transparency framework, and may be proposed as a milestone in the UNFCCC’s inventory assessment for the GST.
Three series of satellites observe the XCO₂ concentration: the GOSAT series, the TanSat series, and the Orbiting Carbon Observatory (OCO) series. Among them, only the GOSAT series can observe the vertical distribution of the XCO₂ concentration by sensing shortwave infrared (SWIR) and thermal infrared spectral band data based on the Fourier transform method. The lower and upper partial columns in the troposphere are identified according to the pressure ranges of 0.6–1 and 0.2–0.6 P_surf, respectively, where P_surf is the pressure at the ground^23,24. XCO₂LT products are advantageous because they are not affected by thin clouds or aerosols²⁴. For instance, XCO₂LT data for Station 6 in the city centre of Ulaanbaatar are retrieved more than twice as frequently as National Institute for Environmental Studies (NIES) GOSAT v02.95-02.97 products. Here, we present the analysis results for anthropogenic CO₂ emissions considering the atmospheric environment in the city of Ulaanbaatar using GOSAT-EORC-Daily-Partial-Column-GHG data from JAXA/EORC²⁶.
A comparison between in-situ CO₂ observations and the results of the WRF-Chem simulation is shown in Supplementary Fig. S1, Supplementary Table S2, and Supplementary Note 1. Figure 2 displays an observation time series of GOSAT in terms of the difference between the averaged XCO₂ concentration in the lower troposphere (i.e., 0.6–1 P_surf in Ulaanbaatar; XCO₂LT) and that in the upper troposphere (i.e., 0.2–0.6 P_surf in Ulaanbaatar; XCO₂UT), which is defined as the XCO₂LT enhancement in Ulaanbaatar. The XCO₂UT concentration is calculated by taking monthly average area-averaged value data as a reference; therefore, the XCO₂LT enhancement is not affected by annual CO₂ increases or seasonal changes. Furthermore, in Ulaanbaatar, the maximum mixed-layer height, i.e., the height of a capping temperature inversion or statically stable layer of air, is approximately 2000 m on a typical clear summer day. In winter, a thick temperature inversion layer (> 500 m) covers a weak and thin mixed layer (< 300 m)²⁷ because radiative cooling at night lowers the near-surface temperature, which is stratified in the absence of convection. Therefore, the actual boundary layer lies well below 4 km in altitude²⁴. In addition, CO₂ emissions from the chimneys of thermal power plants and the ger districts move and diffuse horizontally in the layer below 4 km (Fig. 1, Supplementary Fig. S2). As CO₂ emissions from Ulaanbaatar remain confined to the lower troposphere, a typical condition in a high latitude basin area, XCO₂LT enhancement data can be considered representative of the CO₂ emissions from the city²⁴ (Supplementary Note 2).
Area-averaged enhancement of XCO₂LT in Ulaanbaatar. The enhancement (daily averaged) was calculated by subtracting the observed XCO₂LT from the calculated monthly area-averaged XCO₂UT.
In this study, the model domain is set as slightly larger than the city size (50 km × 50 km and 0.6–1 P_surf; Supplementary Note 3). Fourteen target observed data from GOSAT for XCO₂LT with 10.5 km resolution cover the model domain almost completely. Thus, we assume that such data completely capture the CO₂ emissions trapped in the urban dome. This is especially true in winter, when CO₂ emissions are more easily quantifiable than in summer owing to the thick inversion layer confining the gas to lower altitudes, which contributes to the amount of available observed data from GOSAT. Therefore, there are more available data for winter than for summer.
Figure 2 shows the XCO₂LT enhancement in Ulaanbaatar. GOSAT can observe not only points along the orbit but also target points. The trend between the summer of 2017 and winter of 2018, when GOSAT began targeted observations in Ulaanbaatar, shows that CO₂ concentrations in the lower troposphere were consistently higher than those in the upper troposphere. This confirms that the XCO₂LT enhancement is closely related to the emissions from the city, thereby reflecting its CO₂ emissions. Furthermore, the XCO₂LT enhancement, which became increasingly prominent from fall to winter/spring, shows a seasonal change with a decrease towards the beginning of summer, consistent with a trend of increased emissions in winter and decreased emissions in summer.
Figure 3 shows the results of the inverse analysis of CO₂ emissions in the energy sector in Ulaanbaatar using GOSAT satellite data with a top-down approach. We estimated the a posteriori CO₂ emissions in the model domain surrounding Ulaanbaatar in 2018 using Green’s function between the CO₂ concentration estimated by the atmospheric WRF-Chem transport model (see the comparison between modelled and observed XCO₂ values in Supplementary Fig. S3) and 24 h-averaged XCO₂LT values from GOSAT. In 2018, GOSAT conducted targeted observations in the model domain at 14 target observation points with a spatial resolution of 10.5 km. We only considered XCO₂LT concentration data at times in which GOSAT successfully observed more than 10 target observation points. Finally, we obtained data for 23 days that we could use for the inverse analysis.
Comparison of Mongolia’s bottom-up and top-down energy sector CO₂ emissions in 2018. Mongolia’s CO₂ emissions (in Tg yr⁻¹) for BUR, top-down estimates, and EDGAR. The vertical lines in BUR and top-down columns correspond to 95% confidence intervals.
Green’s function requires a priori errors of GOSAT observations and CO₂ emissions data as inputs. An observation uncertainty for XCO₂LT concentration and the prior emission uncertainty is set to 2 ppmv and 400 ton h⁻¹ (see Method). Additionally, a sensitivity analysis was performed on additional scenarios derived from different combinations with the observation and the prior emission uncertainties (Supplementary Table S3, Supplementary Note 4, and Supplementary Figs. S4–S15). Figure 4 shows the monthly-averaged CO₂ emissions obtained from a posteriori CO₂ emission estimates. The average percentage of a posteriori CO₂ emissions to a priori CO₂ emissions (the modification rate) is only 1.1% throughout the year (Supplementary Table S3). The posterior uncertainty of the scaling factor ranges from 0.12 to 0.34 (Supplementary Table S4), which is comparable to that evaluated in large cities such as Riyadh (0.17–0.24), Cairo (0.10–0.25), and Los Angeles (0.11–0.16)¹⁸. We derived a scaling factor (± 1σ uncertainty) for all observation days throughout the year (Supplementary Table S4), which was approximately 3%. Our result is consistent with the study of Ye et al. (2020)¹⁸ in that utilizing satellite observation over a longer time period potentially obtains emission estimates with lower uncertainty. This uncertainty underscores the value of integrating atmospheric constraints into our model, leading to more stable and reliable emissions estimates.
Inversion results of monthly-averaged CO₂ emissions in Ulaanbaatar using GOSAT observation data in 2018. Each month other than June, September, November, and December features several data points. The vertical lines show the posterior uncertainties.
The a posteriori CO2 emissions estimates are only 0.2% higher than the value expected from a global emission inventory, the Emission Database for Global Atmospheric Research (EDGAR)²⁸. In countries with limited data sources, such as Mongolia, there are cases in which BURs and National Inventory Reports use the same data sources as the global database to produce GHG inventories based on a bottom-up approach. Notably, although emissions in this study were estimated with a particular focus on the energy sector, a priori CO₂ emissions and EDGAR values were already close to each other. Moreover, a posteriori CO₂ emissions based on Green’s function reduce the gap to the EDGAR values even further.
According to Mongolia’s GHG inventory report, the CO₂ emissions in 2014 were as follows: the energy sector accounted for 97.4% of total CO₂ emissions, while industrial processes and product use sector, agriculture, forestry and other land use sector, and waste sector accounted for the remainder. There are no CO₂ emission sources other than the energy sector in Ulaanbaatar. In the energy sector in Ulaanbaatar, the largest emissions are attributed to fuel combustion activities, accounting for 99.99%. Fuel combustion includes energy industries (electricity, heat, and other energy generation), manufacturing and construction industries, residential and commercial activities, and transport. Consequently, the top-down CO₂ emission estimates consist of the energy sector. According to the percentage of CO₂ emissions in Ulaanbaatar to those in the whole of Mongolia (Supplementary Table S1), we assume that this percentage did not change between 2014 and 2018 because economic growth in Mongolia during this period was minimal, partly owing to a conflict with the Mongolian government over taxation and financial agreements related to a foreign mining development project, which led to a pause in the project²⁹. Thus, the economic structure of the capital city and rest of the country remained unchanged. Therefore, we consider that CO₂ emissions from Ulaanbaatar account for 69% of CO₂ emissions in Mongolia in our study.
The top-down CO₂ emissions (20.6 ± 0.3 Tg yr⁻¹) (all reported ranges are 95 percent confidence intervals) and energy sector CO₂ emissions (20.3 ± 1.4 Tg yr⁻¹) in BUR2 for 2018 differed by only 1.5% (0.3 Tg yr⁻¹) (Fig. 3). The top-down CO₂ emissions and energy sector emissions from EDGAR v6.0 differ by 4.2% (21.5 Tg yr⁻¹). EDGAR v6.0’s energy sector CO₂ emissions estimates were 6% higher than Mongolia’s energy sector CO₂ emissions reported in BUR2. Therefore, the top-down approach presented in this study could reduce the difference between the estimate from the global CO₂ emissions database and the value calculated in BUR2 in Mongolia.
We showed good agreement with the CO₂ emissions of the energy sector computed in the upcoming BUR2 (with a discrepancy of only 1.5%). In addition, Mongolia’s a posteriori CO₂ emission estimate was 4.2% smaller than that obtained from EDGAR v6.0²⁸. The low uncertainty in our study can be primarily ascribed to the topographic features in Ulaanbaatar, as well as utilizing lower troposphere data with longer time from satellite observations, and applying detailed emission sources to the chimneys of thermal power plants and the ger districts in the atmospheric transport model.
GOSAT satellites make observations covering an area with a 10.5 km diameter in a 160 km-wide grid. Therefore, there are many unobserved areas in Mongolia, and interpolation in those spaces becomes a problem. In this study, we conducted 14 intensive observations (with a diameter of 10.5 km) to cover the entire city of Ulaanbaatar, Mongolia. The use of an observation method that captures all CO₂ emissions from Ulaanbaatar improves the accuracy of the observations. As such, the inverse analysis makes it possible to estimate the amount of CO₂ emissions from Ulaanbaatar.
Regarding global fossil fuel emissions, one report states that urban areas account for approximately 70% of global fossil fuel emissions¹² and play an essential role in mitigation strategies under the Paris Agreement’s action plan. Ulaanbaatar is the largest city in Mongolia, and its fossil fuel-derived CO₂ emissions account for 69% of the country’s total, which is most of the country’s emissions. Notably, to prevent air pollution, population influx is restricted. Therefore, it is reasonable to assume that the emissions from Ulaanbaatar and Mongolia will remain unchanged¹².
We establish that the inverse analysis accurately estimates 69% of fossil fuel-derived CO₂ emissions in Mongolia. The remaining 31% can, thus, be obtained according to the ratio between the city and nationwide emissions rather than by inverse analysis. This is acceptable for Mongolia as it is a non-Annex I country with limited financial and human resources. “Hybrid Approaches,” combining the latest inverse analysis methods and all other applicable data analysis-based methods, have recently been recommended as a powerful means of reflecting estimated GHG emissions in the policy-making process more quickly⁵.
The Government of Mongolia decided to include the results of the National Emission Inventory 2018 by the “Hybrid Approaches” method described in Chapter 3, 3.1 Inventory overview, Additional information/best practice of the BUR2 to be submitted to the UNFCCC.
This approach is applicable in most, if not all, non-Annex I countries.
The ratio between the city and nationwide emissions of fossil fuel-derived CO₂ emissions could influence technological innovation towards decarbonization in each country in the future. To apply this ratio to countries other than Mongolia, it is necessary to upgrade the latest inverse analysis to a more general emission estimation method for a vast area that comprises not only urban areas but also an entire country. To this end, further development of future GHG observation satellite technology and improvement of the analysis system are necessary.
JAXA/NIES/MOEJ plans to launch the Next-Generation Greenhouse Gases Observing Satellite (GOSAT-GW) in 2024, it will make it possible to make planar observations over a swath of 911 km or more with a resolution of 10 km. Using these next-generation satellites, all areas, including many urban areas and areas with other land uses, will be covered with a 10 km resolution mesh. In addition, targeted observations such as those conducted in Ulaanbaatar in this research will be carried out in all other areas. With the advent of GOSAT-GW, it will become possible to uniformly apply the top-down method using satellite data to all land use areas. To incorporate GOSAT-GW, we are expanding and improving all analysis systems. We are preparing to apply the top-down method using satellite data in a unified manner to countries other than Mongolia.
Our approach is to first clarify the CO₂ emissions from fossil fuel-based energy sources, which are Mongolia’s most important emission sources. However, there are various sinks and emission sources in Mongolia and other regions. Livestock is increasing yearly in Mongolia. Methane emissions from livestock accounted for 47.9% and 48.9% of the total emissions inventories in Mongolia in 2014 and 2018, respectively.
Observations covering a wide area can effectively detect methane emissions and CO₂ absorption and emission from wide-ranging emission sources, such as grasslands and forest areas. Additionally, using GOSAT-GW in these areas, it is expected that a top-down method using satellites will be applicable to estimate the amount of CO₂ absorption from grassland/forest areas and methane emissions from grazing livestock.
Two approaches can be used for estimating GHG emissions: a bottom-up approach relying on forward analysis of global databases²⁸ and top-down approach based on inverse analysis of observational data²⁹.
Top-down approaches rely on inverse models such as Lagrangian particle models³¹ and Bayesian models based on Green’s function³². These models, however, have to contain and rely on a CO₂ atmospheric transport model. These model can be a global transport model, such as the Goddard Earth Observing System model coupled to chemistry (GEOS-Chem)³³, or a regional-scale model based on WRF-Chem³⁴ (i.e., a regional model that relies on appropriately downscaled data extracted from a global transport model).
Ganesan et al.³⁰ and Miller et al.³⁵ reported examples of inverse analyses based on global transport models for India and China, respectively. In both study cases, methane emissions were analyzed by an inverse model using satellite data elaborated with GEOS-Chem and compared with BURs. In terms of CO₂ inversion analysis, there have been studies such as Janarganan et al. (2020)³⁶ who used a global model, and Zhang et al. (2021)³⁷ who utilized a regional model. However, for most of non-Annex I countries like Mongolia, where observational data for CO₂ emissions are limited, application of these analysis methods proves challenging. This is a common issue among non-Annex I countries in the UNFCCC, highlighting the need for developing a method that begins with the analysis at city level, where knowledge and observational data is relatively abundant such as air pollution monitoring data, regional weather model, and meteorological knowledge, which extends them to a nationwide estimation. Several examples of inverse analysis of CO₂ emissions based on regional models targeting cities have been reported in the literature^{18,19,38,39,40}. Pillai et al.³², for instance, used a 100-km square grid to estimate CO₂ emissions in Berlin, Germany. Therefore, in this study, we first focus on a top-down approach to estimate CO2 emissions at the city level, specifically in Ulaanbaatar, Mongolia. This city-level data is then used as a basis for the nationwide CO₂ emission estimation, creating a comprehensive picture of Mongolia’s CO2 emissions.
The process of estimating CO₂ emissions consisted of three sequential steps (Supplementary Fig. S16): forward analysis using the regional chemical transport model WRF-Chem, inverse analysis using Green’s function, and comparison with a global database such as EDGAR.
The spatial distribution of a priori CO₂ emissions in Ulaanbaatar was estimated based on CO₂ emission data from thermal power plants (four sites), automobiles, household stoves in gers, heat-only boilers, and coal-fired water heaters, which accounted for almost all the energy sector-related emissions in Mongolia in 2015. These data were estimated by the Japanese Ministry of Environment in the “Project on Development of Innovative Green Technology and MRV Method for JCM in Mongolia”, which was based on the same methodology as the Mongolian GHG emission estimate. Comparing the CO₂ inventory for the energy sector in Mongolia from EDGAR with Mongolia’s gross domestic product (GDP) for the same period, a constant ratio can be evaluated, suggesting that CO₂ emissions in the energy sector are strongly related to the economic growth in Mongolia, and that no changes in economic structure occurred during the monitoring period. We calculated the a priori total CO₂ emissions from Ulaanbaatar in 2018 according to this constant ratio with GDP⁴¹ (Supplementary Fig. S17). We calculated the hourly input data of a prior CO₂ emissions for the forward model by differentiating observed data on hourly CO concentration changes provided by the Information and Research Institute of Meteorology, Hydrology and Environment in Mongolia (Supplementary Note 5, Supplementary Fig. S18).
We used the GHG Tracer developed by the Max Planck Institute and the National Oceanic and Atmospheric Administration as an atmospheric transport model for CO₂. Consequently, we were able to sequentially analyse CO₂ concentrations with a horizontal resolution of 9 km and 34 vertical levels. The modelled XCO₂ was in good agreement with the GOSAT XCO₂ data, with a correlation coefficient of 0.89. Supplementary Table S5 provides details on WRF-Chem physical conditions and main input parameters, including the modelling of the CO₂ flux from vegetation (Vegetation Photosynthesis and Respiration Model, VPRM)^42,43, and to Supplementary Fig. S3 shows the simulation results of XCO₂.
Ye et al.¹⁸ estimated the error variance by combining the measurement error and forward model error variances. The former accounts for GOSAT measurements, whereas the latter considers transport model errors in wind speed, wind direction, and boundary conditions. WRF-Chem transport model uncertainties range within 0.31–3.03 ppm in Riyadh, Cairo, and Los Angeles¹⁸. In our model, we considered uncertainties in CO₂ boundary conditions and the prior emissions.
Regarding boundary conditions, meteorological data and background CO₂ concentration data—taken as boundary conditions in WRF-Chem—were obtained, from the final analysis data from the National Center for Environmental Prediction, USA⁴⁴ and JENA CarboScope inversion system JENA_s04oc_v4.3 (shortened as JENA_s04oc)⁴⁵ data from the Max Planck Institute, Germany. Although the resolution of JENA_s04oc was lower in our region (4° × 5°), the CO₂ distribution was optimized using a Bayesian inversion framework^46,47. The monthly average difference with the Atmospheric Infrared Sounder product of NASA’s AQUA satellite was approximately 1.1 ppmv in terms of column CO₂ concentration in the Mongolian region (44–49 N, 100–110 E)⁴⁸. In situ CO₂ measurements by aircraft over Europe showed that JENA_s04oc products had a lower bias than those from the Copernicus Atmosphere Monitoring Service (CAMS)⁴⁹, 0.8 (1.3) μmol mol^-1 vs. 3.7 (1.5) μmol mol⁻¹⁵⁰, where the standard uncertainty in the final digits is given in brackets. This suggests that CAMS products are preferred for setting lateral boundary conditions for regional modelling. Consequently, we opted for JENA_s04oc products because the bias of CAMS in Mongolia was poorly constrained.
Regarding CO₂ emissions, we quantified the discrepancy in CO₂ concentrations between observed and model-predicted values (Supplementary Fig. S3).
We conducted the inverse analysis through a Bayesian inversion based on the synthesis inversion method (Green’s function method)³² and constructed Green’s function to modify CO₂ emissions to minimize the difference between the WRF-Chem simulated estimate and GOSAT observations by the following cost function J. This resulted in an a posteriori estimate of CO₂ emissions.
where c_obs is the observed CO₂ mole fraction, c_fwd is the simulated CO₂ mole fraction, H is Jacobi’s matrix (number of emission sources × number of observations per optimization period; e.g., 1 day), R is the c_obs—c_fwd error covariance matrix, x is the modified amount of CO₂ at the observation point, x₀ is the prior emission estimate of x, and P₀ is the error covariance matrix of x₀.
In (1), the measurement error covariance matrix (diagonal matrix), R depends on the XCO₂ random errors³² using GOSAT and TCCON observation data collected independently⁵¹. As a reference for R, Yoshida et al.⁵¹ showed that the standard deviation of the differences between XCO₂ (SWIR L2 v02.xx product) and TCCON data are evaluated to be 2.1 ppm. Wang et al.⁵² also stated that using GOSAT data (ACOS dataset), setting 2 times (land) and 1.25 times (sea) of the presented retrieval error would give realistic errors. In (1), P0 represents the error covariance of the a priori CO₂ emission and is defined as a diagonal matrix when the CO₂ emission data as referred to by Pillai et al.³² are independently generated.
In this case, the optimal solution (a posteriori estimate of CO₂ emissions) is given by:
Jacobi’s matrix H³², which corresponds to Green’s function, is a function of the impact of CO₂ emissions from a given point on Earth’s surface on the spatiotemporal variation in atmospheric concentrations. Therefore, forward modelling results based on an atmospheric transport model inputting pulsed emission fluxes or a Lagrangian analysis tracking the trajectories of CO₂ gas can be used to populate H. The period during which emissions from the ground surface remain within the targeting area varies depending on the size of the targeting area and meteorological conditions (such as the wind speed). Therefore, we must formulate H for the period during which CO₂ gas concentrations are recorded after the gas is emitted from within the model domain (Supplementary Fig. S19).
In R, Kuze et al.²⁴ used a standard deviation of 2.09 ppm for both XCO₂LT and XCO₂UT because the value was validated by the TCCON network. Therefore, we also set a standard deviation of 2 ppm for our a priori error in the XCO₂LT concentration data, which is also close to the value provided by the NIES⁵¹.
For P_0, Pillai et al.³² gave a 40% a priori emission error in CO₂ emissions based on the difference between their detailed inventory data and the EDGAR inventory. Ye et al.¹⁸ chose an uncertainty of 20% for cities with low statistical uncertainty in CO₂ emissions (Los Angeles) and 40% for cities with high statistical uncertainty (Riyadh, Cairo). In Ulaanbaatar, the statistical uncertainty is assumed to be high; however, the city is smaller than Riyadh and Cairo, and the number and location of primary sources of CO₂ are known. Moreover, comparing the a priori CO₂ emission data in this study with EDGAR energy sector data, the difference is as small as 1.3%. Therefore, we set the a priori emission uncertainty to 400 ton h⁻¹ (26%).
In this study, Eqs. (1) and (2) were applied to a 50 km × 50 km model domain (refer to the Results, Supplementary Note 3, and Supplementary Fig. S20) in and around Ulaanbaatar in 2018 to calculate a posteriori estimates of GHG emissions, and XCO₂LT concentrations were calculated from GOSAT/GOSAT-2 EORC Daily Partial Column GHG data^24,25. CO₂ emission values in Green’s function corresponded to hourly-averaged emissions within the model domain over 24 h from the observation time of GOSAT. We calculated the a posteriori emissions only if 70% of the GOSAT observation points (10 out of 14 points) were available. Finally, we calculated the averaged ratio of a posteriori to a priori CO₂ emissions within the model domain over the entire year.
After solving the cost function J, the following two formulas¹⁸ can be used to calculate an optimized estimate of the scaling factor and the a posteriori error variance:
where (widehat{lambda }=widehat{x}), y_m = H, S_o = R, and σ_a=({{varvec{P}}}_{0}) based on (1) and (2).
Here, we describe the method to calculate a posteriori estimates of CO₂ emissions for Ulaanbaatar and Mongolia.
First, we calculated the average ratio of a posteriori to a priori CO₂ emissions at each GOSAT observation point. Second, we multiplied this ratio by the a priori estimate of energy sector CO₂ emissions for Ulaanbaatar and Mongolia, thereby obtaining the a posteriori CO₂ emissions of the energy sector. Finally, we calculated the 95% confidence interval (across all observation points during the whole of 2018) of the difference between a priori and a posteriori emissions at each observation point.
We calculated the a posteriori estimates of CO₂ emissions in Mongolia using the following procedure. First, we calculated the a priori estimate of CO₂ emissions in 2018 in Mongolia by multiplying the CO₂ emissions in 2014 from the 1^st BUR in Mongolia by the 2018 GDP/2014 GDP³⁹ ratio. Second, we calculated the a posteriori estimate of CO₂ emissions in 2018 in Mongolia by multiplying the averaged ratio of a posteriori to a priori CO₂ emissions within the model domain in Sect. 2 by the a priori estimate of CO₂ in 2018 for Mongolia. Third, we calculated the difference between top-down and bottom-up CO₂ emissions by calculating the difference between the a posteriori estimate of CO₂ emissions and energy sector CO₂ emissions in BUR2 for 2018. We used ArcGIS 10.2 (https://www.esri.com/en-us/arcgis/about-arcgis/overview) to extract subsets for Ulaanbaatar and Mongolia from the global CO₂ emissions dataset in EDGAR v6.0. The latter contains energy sector CO₂ emissions data for each administrative district. Finally, we calculated the difference between each subset and the a posteriori estimate of CO₂ emissions.
GOSAT XCO2LT data were obtained from https://www.eorc.jaxa.jp/GOSAT/CO2_monitor/index_Ver.K.html. Observed CO and CO₂ data are available from the corresponding authors on reasonable request.
Masson-Delmotte, V. et al. IPCC, 2021: Summary for policymakers. In Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (eds Masson-Delmotte, V. et al.) 2–32 (Cambridge University Press, 2021).
Google Scholar 
UNFCCC. Proposal by the president Conference of the Parties (COP), Adoption of the Paris Agreement. http://unfccc.int/resource/docs/2015/cop21/eng/l09.pdf (2015).
Tsukui, A., Murun, T. & Umemiya, C. Event: Round table for “from satellite observation to national GHG inventories: inputs to the Global Stocktake of the Paris Agreement”. https://www.iges.or.jp/en/pub/gst-satellite-roundtable-summary-report/en (2022).
IPCC 2006. 2006 IPCC Guidelines for National Greenhouse Gas Inventories (eds Eggleston, H. S. et al.) (IGES, 2006).
National Academies of Sciences, Engineering, and Medicine. Greenhouse Gas Emissions Information for Decision Making: A Framework Going Forward (National Academies Press, 2022).
Google Scholar 
Hattori, T. & Umemiya, C. IGES Biennial Update Report (BUR) Database. https://pub.iges.or.jp/pub/iges-bur-database (2022)
Deng, Z. et al. Comparing national greenhouse gas budgets reported in UNFCCC inventories against atmospheric inversions. Earth Syst. Sci. Data 14, 1639–1675 (2022).
Article  ADS  Google Scholar 
Calvo Buendia, E. et al. 2019 Refinement to the 2006 IPCC Guidelines for National Greenhouse Gas Inventories (2019).
BP. p.l.c. Statistical Review of World Energy – All Data, 1965–2021. https://www.bp.com/en/global/corporate/energy-economics/statistical-review-of-world-energy/using-the-review/definitions-and-explanatory-notes.html#accordion_Tonnes (2022).
National Statistics Office of Mongolia. BALANCE OF COAL, by national level, by year. https://1212.mn/tables.aspx?TBL_ID=DT_NSO_1100_010V1 (2022).
Climate Change Research and Cooperation Centre, Ministry of Environment and Tourism, Mongolia: Mongoila’s Second Biennial Update Report under the United Nation Framework Convention on Climate Change, to be submitted in October 2023.
Seto, K. C. et al. Human settlements, infrastructure, and spatial planning. In Clim. Change: Mitigation of Climate Change. Contribution of Working Group III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press (2014).
The Climate Change Project Implementing Unit of Environment and Climate Fund under the Ministry of Environment and Tourism: Mongolia’s Initial Biennial Update Report under the United Nations Framework Convention on Climate Change. (2017).
Sternberg, T. Investigating the presumed causal links between drought and dzud in Mongolia. Nat. Hazards 92, 27–43 (2018).
Article  Google Scholar 
Guttikunda, S. K. et al. Particulate pollution in Ulaanbaatar, Mongolia. Air Qual. Atmos. Health 6, 589–601 (2013).
Article  CAS  Google Scholar 
Kohn, M. World’s Worst Air Has Mongolians Seeing Red, Planning Action. https://www.bloomberg.com/news/articles/2016-12-22/world-s-worst-air-has-mongolians-seeing-red-planning-protest?leadSource=uverify%20wall (2016).
Allen, R. W. et al. An assessment of air pollution and its attributable mortality in Ulaanbaatar, Mongolia. Air Qual. Atmos. Health 6, 137–150 (2013).
Article  CAS  PubMed  Google Scholar 
Ye, X. et al. Constraining fossil fuel CO₂ emissions from urban area using OCO-2 observations of total column CO₂. J. Geophys. Res. Atmos. 125, 1–29 (2020).
Article  Google Scholar 
Lei, R. et al. Fossil fuel CO₂ emissions over metropolitan areas from space: A multi-model analysis of OCO-2 data over Lahore, Pakistan. Remote Sens. Environ. 264, 1–11 (2021).
Article  Google Scholar 
Lespinas, F. et al. The potential of a constellation of low earth orbit satellite imagers to monitor worldwide fossil fuel CO₂ emissions from large cities and point sources. Carbon Balance Manag. 15, 1 (2020).
Article  Google Scholar 
Zheng, T., Nassar, R. & Baxter, M. Estimating power plant CO₂ emission using OCO-2 XCO₂ and high resolution WRF-Chem simulations. Environ. Res. Lett. 14, 085001 (2019).
Article  ADS  CAS  Google Scholar 
Broquet, G. et al. The potential of satellite Spectro-imagery for monitoring CO₂ emissions from large cities. Atmos. Meas. Tech. 11, 681–708 (2018).
Article  CAS  Google Scholar 
Kikuchi, N., Yoshida, Y., Uchino, O., Morino, I. & Yokota, T. An advanced retrieval algorithm for greenhouse gases using polarization information measured by GOSAT TANSO-FTS SWIR I: Simulation study. J. Geophys. Res. Atmos. 121, 13–129 (2016).
Article  Google Scholar 
Kuze, A. et al. Examining partial-column density retrieval of lower-tropospheric CO₂ from GOSAT target observations over global megacities. Remote Sens. Environ. 273, 112966 (2022).
Article  ADS  Google Scholar 
Aganaba-Jeanty, T. & Huggins, A. Satellite measurement of GHG emissions: Prospects for enhancing transparency and answerability under international law. Transnat. Environ. Law 8, 303–326 (2019).
Article  Google Scholar 
Japan Aerospace Exploration Agency Earth Observation Research Center (JAXA/EORC). GOSAT/GOSAT-2 EORC Daily Partial Column GHGs. https://www.eorc.jaxa.jp/GOSAT/GPCG/index_GOSAT.html (2022).
Wang, M., Kai, K., Sugimoto, N. & Enkhmaa, S. Meteorological factors affecting winter particulate air pollution in Ulaanbaatar from 2008 to 2016. Asian J. Atmos. Environ. 12, 244–254 (2018).
Article  CAS  Google Scholar 
Crippa, M. et al. EDGAR v6.0 greenhouse gas emissions. European Commission, Joint Research Centre (JRC). http://data.europa.eu/89h/97a67d67-c62e-4826-b873-9d972c4f670b (2021).
Batjargal, D. & Doojav, G. O. Measuring Sudden Stops in Mongolia. https://mpra.ub.uni-muenchen.de/105115/1/MPRA_paper_105115.pdf (2020).
Ganesan, A. L. et al. Atmospheric observations show accurate reporting and little growth in India’s methane emissions. Nat. Commun. 8, 836 (2017).
Article  ADS  PubMed  PubMed Central  Google Scholar 
Lauvaux, T. et al. High-resolution atmospheric inversion of urban CO₂ emissions during the dormant season of the Indianapolis flux experiment (INFLUX). J. Geophys. Res. Atmos. 121, 5213–5236 (2016).
Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 
Pillai, D. et al. Tracking city CO₂ emissions from space using a high-resolution inverse modelling approach: A case study for Berlin. Germany. Atmos. Chem. Phys. 16, 9591–9610 (2016).
Article  ADS  CAS  Google Scholar 
Henze, D. K., Hakami, A. & Seinfeld, J. H. Development of the adjoint of GEOS-Chem. Atmos. Chem. Phys. 7, 2413–2433 (2007).
Article  ADS  CAS  Google Scholar 
Beck, V. et al. WRF-Chem simulations in the Amazon region during wet and dry season transitions: Evaluation of methane models and wetland inundation maps. Atmos. Chem. Phys. 13, 7961–7982 (2013).
Article  ADS  Google Scholar 
Miller, S. M. et al. China’s coal mine methane regulations have not curbed growing emissions. Nat. Commun. 10(1), 303 (2019).
Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 
Janardanan, J. et al. Country-scale analysis of methane emissions with a high-resolution inverse model using GOSAT and surface observations. Remote Sens. 12, 375 (2020).
Article  ADS  Google Scholar 
Zhang, Q. et al. CO₂ flux over the contiguous United States in 2016 inverted by WRF-Chem/DART from OCO-2 XCO₂ retrievals. Remote Sens. 13, 2996 (2021).
Article  ADS  Google Scholar 
Pitt, J. R. et al. New York City greenhouse gas emissions estimated with inverse modeling of aircraft measurements. Elementa 10, 1–13 (2022).
Google Scholar 
Lian, J. et al. Assessing the effectiveness of an urban CO₂ monitoring network over the Paris region through the COVID-19 lockdown natural experiment. Environ. Sci. Technol. 56, 2153–2162 (2022).
Article  ADS  CAS  PubMed  Google Scholar 
Brioude, J. et al. Top-down estimate of surface flux in the Los Angeles Basin using a mesoscale inverse modeling technique: Assessing anthropogenic emissions of CO, NOx and CO₂ and their impacts. Atmos. Chem. Phys. 13, 3661–3677 (2013).
Article  ADS  Google Scholar 
World Bank. World Bank national accounts data, and OECD National Accounts data files (constant LCU). https://data.worldbank.org/indicator/NY.GDP.MKTP.KN?locations=MN (2022).
Mahadevan, P. et al. A satellite-based biosphere parameterization for net ecosystem CO₂ exchange: Vegetation Photosynthesis and Respiration Model (VPRM). Glob. Biogeochem. Cycles 22, GB2005 (2008).
Article  ADS  Google Scholar 
Ahmadov, R. et al. Comparing high resolution WRF-VPRM simulations and two global CO₂ transport models with coastal tower measurements of CO₂. Biogeosciences 6, 807–817 (2009).
Article  ADS  CAS  Google Scholar 
National Center for Atmospheric Research. NCEP FNL operational model global tropospheric analyses, continuing from July 1999. https://rda.ucar.edu/datasets/ds083.2/ (2022).
Max Planck Institute for Biogeochemistry. Atmospheric CO₂ Inversion s04oc_v4.3 Atmospheric CO₂ fields. https://www.bgc-jena.mpg.de/CarboScope/?ID=s (2022).
Rödenbeck, C. Estimating CO₂ Sources and Sinks from Atmospheric Mixing Ratio Measurements Using a Global Inversion of Atmospheric Transport. Technical Report No. 6 (Max-Planck-Institut für Biogeochemie, 2005).
Rödenbeck, C., Houweling, S., Gloor, M. & Heimann, M. CO₂ flux history 1982–2001 inferred from atmospheric data using a global inversion of atmospheric transport. Atmos. Chem. Phys. 3, 1919–1964 (2003).
Article  ADS  Google Scholar 
Tiwari, Y. K. et al. Comparing CO₂ retrieved from Atmospheric infrared Sounder with model predictions: Implications for constraining surface fluxes and lower-to-upper troposphere transport. J. Geophys. Res. 111, D17106 (2006).
Article  ADS  Google Scholar 
Inness, A. et al. The CAMS reanalysis of atmospheric composition. Atmos. Chem. Phys. 19, 3515–3556 (2019).
Article  ADS  CAS  Google Scholar 
Gałkowski, M. et al. In situ observations of greenhouse gases over Europe during the CoMet 1.0 campaign aboard the HALO aircraft. Atmos. Meas. Tech. 14, 1525–1544 (2021).
Article  Google Scholar 
Yoshida, Y. et al. Improvement of the retrieval algorithm for GOSAT SWIR XCO₂ and XCH₄ and their validation using TCCON data. Atmos. Meas. Tech. 6, 1533–1547 (2013).
Article  Google Scholar 
Wang, J. S. et al. A global synthesis inversion analysis of recent variability in CO₂ fluxes using GOSAT and in situ observations. Atmos. Chem. Phys. 18, 11097–11124 (2018).
Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 
Download references
The authors acknowledge support from the Ministry of Environment of Japan through the “Project on Development of Innovative Green Technology and MRV Method for JCM in Mongolia”, “Project on Assessment of the accuracy of GHG emission estimates using GOSAT-series satellite data in Mongolia”, and “Project on Upgrading Techniques for Estimating Greenhouse Gas Emissions using Satellite Observation data in Mongolia”. The authors thank the National Institute for Environmental Studies (NIES) and Japan Aerospace Exploration Agency (JAXA) for conducting the GOSAT target observations in Ulaanbaatar, Mongolia. The authors also thank Tomohiro ODA for invaluable scientific feedback on the manuscript.
Research and Development Initiative, Chuo University, Tokyo, 1128551, Japan
Masataka Watanabe, Akihiro Oba & Yoko Saito
Information and Research Institute of Meteorology, Hydrology and Environment, Ulaanbaatar, 15160, Mongolia
Gomboluudev Purevjav, Batjargal Gankhuyag & Munkhbat Byamba-Ochir
Climate Change Research and Cooperation Centre, Ministry of Environment and Tourism, Mongolia, Ulaanbaatar, 14191, Mongolia
Batjargal Zamba
Graduate School of Science and Engineering, Chuo University, Tokyo, 1128551, Japan
Tomohiro Shishime
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
M.W. developed the concept; M.W., A.O., and Y.S. performed the data analysis and wrote the manuscript; A.O., G.P., and B.G. developed the atmospheric model. M.B. and Y.S. collected and analyzed the atmospheric samples. B.Z. provided the national CO₂ emissions data. T.S. was responsible for research ethics and project management. All authors reviewed the final manuscript.
Correspondence to Masataka Watanabe.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Watanabe, M., Oba, A., Saito, Y. et al. Enhancing scientific transparency in national CO₂ emissions reports via satellite-based a posteriori estimates. Sci Rep 13, 15427 (2023). https://doi.org/10.1038/s41598-023-42664-3
Download citation
Received: 28 January 2023
Accepted: 13 September 2023
Published: 18 September 2023
DOI: https://doi.org/10.1038/s41598-023-42664-3
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
© 2023 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

source