Remote Sensing Monitoring and Analysis of Spatiotemporal Changes in China’s Anthropogenic Carbon Emissions Basedon XCO2 Data

The monitoring and analysis of the spatiotemporal distribution of anthropogenic carbon emissions is an important part of realizing China’s regional “dual carbon” goals; that is, the aim is for carbon emissions to peak in 2030 an to achieve carbon neutrality by 2060, as well as achieving sustainable development of the ecological environment. The column-averaged CO2 dry air mole fraction (XCO2) of greenhouse gas remote sensing satellites has been widely used to monitor anthropogenic carbon emissions. However, selecting a reasonable background region to eliminate the influence of uncertainty factors is still an important challenge to monitor anthropogenic carbon emissions by using XCO2. Aiming at the problems of the imprecise selection of background regions, this study proposes to enhance the anthropogenic carbon emission signal in the XCO2 by using the regional comparison method based on the idea of zoning. First, this study determines the background region based on the Open-Data Inventory for Anthropogenic Carbon dioxide (ODIAC) dataset and potential temperature data. Second, the average value of the XCO2 in the background area was extracted and taken as the XCO2 background. On this basis, the XCO2 anomaly (XCO2ano) was obtained by regional comparison method. Finally, the spatiotemporal variation characteristics and trends of XCO2ano were analyzed, and the correlations between the number of residential areas and fossil fuel emissions were calculated. The results of the satellite observation data experiments over China from 2010 to 2020 show that the XCO2ano and anthropogenic carbon emissions have similar spatial distribution patterns. The XCO2ano in China changed significantly and was in a positive growth trend as a whole. The XCO2ano values have a certain positive correlation with the number of residential areas and observations of fossil fuel emissions. The purpose of this research is to enhance the anthropogenic carbon emission signals in satellite observation XCO2 data by combining ODIAC data and potential temperature data, achieve the remote sensing monitoring and analysis of spatiotemporal changes in anthropogenic carbon emissions over China, and provide technical support for the policies and paths of regional carbon emission reductions and ecological environmental protection.


Introduction
Carbon dioxide (CO 2 ) can remain in the air for hundreds to thousands of years and is a major greenhouse gas [1]. The continuous increase in atmospheric CO 2 caused by human activities since the Industrial Revolution is one of the main reasons for global warming [2][3][4]. With the rapid development of the global social economy, especially industry, a large amount of CO 2 produced by energy consumption is absorbed by terrestrial and Related studies have widely used satellite XCO 2 data products to monitor and evaluate anthropogenic carbon emissions. Some of these studies select the median or average value of the XCO 2 in the study area as the background to enhance the anthropogenic carbon emission signal and achieve the monitoring of anthropogenic carbon emissions [29][30][31][32]. In other studies, the median XCO 2 in the less anthropogenically affected area around the study area was selected as the background [33][34][35][36]. In addition, there are related studies that will build an anthropogenic carbon emission model based on the estimation of anomalous XCO 2 (XCO 2ano ) obtained from satellite observations in order to estimate and monitor anthropogenic carbon emissions [30,37]. The above research results show the feasibility of remote sensing monitoring of anthropogenic carbon emissions, but the selection of background needs to be further discussed. Although most of the relevant studies have shown that the anthropogenic carbon emission signal in the atmosphere is susceptible to interference, less consideration is given to selecting a background area in combination with atmospheric transmission and other factors to eliminate background interference.
In the current study, the determination of background region of XCO 2 is the key content of extracting anthropogenic carbon emission signal from XCO 2 . However, few previous studies combined multisource data to eliminate the influence of XCO 2 background. In addition, the fine analysis of the spatiotemporal characteristics of the XCO 2ano was insufficient, and the relevant influencing factors affecting the anthropogenic carbon emission signal were not explored. Therefore, there are three main objectives of the study: (1) to combine multisource data to reasonably select the background area to obtain the XCO 2ano that can enhance the anthropogenic carbon emission signal; (2) to finely analyze the spatial and temporal distribution characteristics of XCO 2ano ; and (3) to investigate the influencing factors that affect the extraction of anthropogenic carbon emission signal.
In view of the above problems and deficiencies, this study combined multisource data to select the background region, obtained the monthly XCO 2ano results in China from 2010 to 2020 by the selected region comparison method, and analyzed the spatiotemporal characteristics and trends of the XCO 2ano . The feasibility of using satellite observation data to monitor and evaluate anthropogenic carbon emissions was explored, and it provided technical support for the subsequent monitoring and evaluation of carbon emission reduction and regional sustainable development. This study can provide methodological support for targeted spatially differentiated carbon reduction measures.

Data Sources
The experimental data in this study include the monthly Mapping-XCO 2 dataset, annual ODIAC dataset and potential temperature data from 2010 to 2020.
(1) Mapping-XCO 2 dataset This study used the monthly Global Land Mapping-XCO 2 (Mapping-XCO 2 ) dataset from January 2010 to December 2020 which has a spatial resolution of 1 • latitude × 1 • longitude. The Mapping-XCO 2 dataset was generated by applying weighted spatiotemporal kriging interpolation methods to XCO 2 obtained from the GOSAT satellite (January 2010 to August 2014) and the OCO-2 satellite (September 2014 to December 2020) [38,39]. The XCO 2 products include the ACOS Level 2 Lite data product (v9r) from the GOSAT satellite and the OCO-2 Level 2 Lite data product (v10r) from the OCO-2 satellite. The generation of the Mapping-XCO 2 dataset consists of three steps: (1) Adjusting the a priori CO 2 profiles of satellite XCO 2 retrievals. The prior CO 2 profile affects the inversion of XCO 2 data. (2) Adjusting the observing time of satellite XCO 2 data. Observation time is the local overpass times of the satellite. (3) Unifying spatiotemporal scales of satellite observations. Spatiotemporal scale is the data time interval and spatial resolution. He et al. [38] showed that the estimation uncertainty for the Mapping-XCO 2 dataset is small, cross-validation shows that the exact weighted spatiotemporal kriging interpolation method used has good reliability, and the Mapping-XCO 2 dataset has high accuracy. (2) ODIAC dataset The ODIAC dataset is a global anthropogenic carbon emissions product that estimates carbon emissions from fossil fuel combustion based on national-level fossil fuel carbon emissions estimates, fossil fuel consumption statistics, satellite-observed night-time light data, and fossil fuel point source information [40,41]. This dataset can effectively reflect the spatiotemporal distribution of anthropogenic carbon emissions. The latest dataset, version ODIAC2020b, provides the monthly spatial distribution of anthropogenic carbon emissions from 2000 to 2019. This study selects a dataset with a spatial resolution of 1 • latitude × 1 • longitude. In this study, monthly ODIAC datasets from 2000 to 2019 were downloaded from the National Institute of Environmental Research Center for Global Environmental Studies.
(3) Potential temperature data The potential temperature data used in this study came from the National Center for Environmental Prediction/National Center for Atmospheric Research. Using the monthly potential temperature data at an atmospheric pressure level of 1000 mb for reanalysis, the average potential temperature data from January 2010 to December 2020 were obtained, and the average potential temperature contours were regenerated. Potential temperature is a dynamic tracer of stable air mass transport [42] and is not affected by the physical lifting or sinking associated with flow over obstacles or large-scale atmospheric turbulence [43]. The latitudinal and zonal spatial distribution patterns of XCO 2 have a high degree of similarity with the distribution of potential temperature contours. Using the potential temperature contours to divide the study area into different potential temperature zones can eliminate the influence of atmospheric transport on the extraction of anthropogenic carbon emission signals to a certain extent. This study divides the study area into five zones with 10K intervals, area I, area II, area III, area IV, and area V. The specific regions are shown in Figure 1.
(2) ODIAC dataset The ODIAC dataset is a global anthropogenic carbon emissions product mates carbon emissions from fossil fuel combustion based on national-level fossi bon emissions estimates, fossil fuel consumption statistics, satellite-observed n light data, and fossil fuel point source information [40,41]. This dataset can effec flect the spatiotemporal distribution of anthropogenic carbon emissions. The taset, version ODIAC2020b, provides the monthly spatial distribution of anthr carbon emissions from 2000 to 2019. This study selects a dataset with a spatial r of 1° latitude × 1° longitude. In this study, monthly ODIAC datasets from 200 were downloaded from the National Institute of Environmental Research C Global Environmental Studies.
(3) Potential temperature data The potential temperature data used in this study came from the National C Environmental Prediction/National Center for Atmospheric Research. Using the potential temperature data at an atmospheric pressure level of 1000 mb for reana average potential temperature data from January 2010 to December 2020 were and the average potential temperature contours were regenerated. Potential tem is a dynamic tracer of stable air mass transport [42] and is not affected by the lifting or sinking associated with flow over obstacles or large-scale atmospher lence [43]. The latitudinal and zonal spatial distribution patterns of XCO2 have a gree of similarity with the distribution of potential temperature contours. Usin tential temperature contours to divide the study area into different potential tem zones can eliminate the influence of atmospheric transport on the extraction of a genic carbon emission signals to a certain extent. This study divides the study five zones with 10K intervals, area I, area II, area III, area IV, and area V. The sp gions are shown in Figure 1.  The China vector map data used in this research comes from the National Geographic Information Resource Catalog Service System, which belongs to the data of China's 1:1 million national basic geographic databases. This China vector map has been reviewed by the Ministry of Natural Resources of China, and the research area of this study is the range of the China vector map.
The land use data are obtained from the European Space Agency (ESA) Climate Change Initiative (CCI) project. The land use data include data on 22 types of land, including urban land, water bodies and grassland, with a resolution of 300 m. We integrated land cover types into six categories: urban, cropland, vegetation, bare areas, Permanent snow, and ice and water (Table 1). The residential area data used in this research comes from the National Geographic Information Resource Catalog Service System, which belongs to the data of China's 1:1 million public basic geographic information data (2021). Residential areas are places where people gather and settle down. The main elements of residential areas in this study include houses, sheds, cave dwellings, yurts, grazing spots, and other residential buildings.

Research Methods
In this study, based on Mapping-XCO 2 , potential temperature data and ODIAC data, a calculation method of XCO 2ano was designed and constructed, and the spatiotemporal variation characteristics and trends of XCO 2ano results were analyzed. To analyze the potential of the XCO 2ano to monitor anthropogenic carbon emissions, a correlation analysis of the XCO 2ano with fossil fuel emissions and residential area is finally carried out. Correlation analysis was used to quantitatively analyze the ability of the XCO 2ano to monitor anthropogenic carbon emissions. The research process is shown in Figure 2.
Remote Sens. 2023, 15, 3207 6 of 2 variation characteristics and trends of XCO2ano results were analyzed. To analyze the po tential of the XCO2ano to monitor anthropogenic carbon emissions, a correlation analysi of the XCO2ano with fossil fuel emissions and residential area is finally carried out. Corre lation analysis was used to quantitatively analyze the ability of the XCO2ano to monito anthropogenic carbon emissions. The research process is shown in Figure 2. (1) XCO2ano calculation In this study, the time series of average XCO2 in emission zones and background zones in China was analyzed from 2010 to 2020 ( Figure 3). The atmospheric CO2 has stron seasonal variation. The linear fitting degree of atmospheric CO2 and time variables in th emission area and background area is relatively high, and the goodness of fit R 2 is greate than 0.9, indicating that atmospheric CO2 increased during the study period. From Figur 3, it can be found that the atmospheric CO2 has a stable periodic seasonal variation pattern This seasonal variation represents a strong background signal of atmospheric CO2, whic seriously affects the ability of satellites to observe anthropogenic carbon emissions. I summer, XCO2 in the background area is more than that in the emission area. This is be cause the background region contains the Tibetan Plateau, where XCO2 mainly come from the upper troposphere, whose seasonal signal is weaker than that of the lower terrai troposphere. As can be seen from Figure 3b, XCO2 in the background area in summer i (1) XCO 2ano calculation In this study, the time series of average XCO 2 in emission zones and background zones in China was analyzed from 2010 to 2020 ( Figure 3). The atmospheric CO 2 has strong seasonal variation. The linear fitting degree of atmospheric CO 2 and time variables in the emission area and background area is relatively high, and the goodness of fit R 2 is greater than 0.9, indicating that atmospheric CO 2 increased during the study period. From Figure 3, it can be found that the atmospheric CO 2 has a stable periodic seasonal variation pattern. This seasonal variation represents a strong background signal of atmospheric CO 2 , which seriously affects the ability of satellites to observe anthropogenic carbon emissions. In summer, XCO 2 in the background area is more than that in the emission area. This is Remote Sens. 2023, 15, 3207 7 of 25 because the background region contains the Tibetan Plateau, where XCO 2 mainly comes from the upper troposphere, whose seasonal signal is weaker than that of the lower terrain troposphere. As can be seen from Figure 3b, XCO 2 in the background area in summer is similar to or lower than the emission area after removing the area with an altitude higher than 3000 m over the Tibetan Plateau.
To weaken the background signal of CO 2 and enhance the anthropogenic carbon emission signal, this study selects the regional comparison method to obtain the XCO 2ano in the anthropogenic carbon emission area. Most current studies use regional comparison method to remove the influence of XCO 2 background and enhance anthropogenic emission signals [35]. The regional comparison method uses the difference in XCO 2 between the anthropogenic carbon emission area and the background area as the XCO 2ano . The key step of the regional comparison method is to select the background area.  To weaken the background signal of CO2 and enhance the anthropogenic carbon emission signal, this study selects the regional comparison method to obtain the XCO2ano in the anthropogenic carbon emission area. Most current studies use regional comparison method to remove the influence of XCO2 background and enhance anthropogenic emission signals [35]. The regional comparison method uses the difference in XCO2 between the anthropogenic carbon emission area and the background area as the XCO2ano. The key step of the regional comparison method is to select the background area.
This study combines potential temperature data and ODIAC data to determine the background area. First, the ODIAC data is used to identify the background areas by screening out the areas without anthropogenic carbon emissions, and the remaining area is the anthropogenic carbon emission area. According to the monthly ODIAC data, the areas without anthropogenic carbon emissions are consistent in each month from 2010 to 2019, which are mainly distributed in the less-traveled areas of Tibet, Xinjiang, Inner Mongolia, and Northeast China. Second, China is divided into 5 regions by using potential temperature data. In this study, the XCO2ano in different potential temperature regions was calculated according to the background regions of different potential temperature regions. In simulations with zonally uniform surface fluxes, XCO2 is tightly correlated with potential temperature [43]. However, the influence of topographic and climatic factors may lead to differences in atmospheric transport modes in some regions within the same potential temperature isoline. Therefore, it is difficult to completely eliminate the influence of atmospheric transport. The average potential temperature data of atmospheric altitude of 1000 mb was selected for analysis in this study. However, there are differences in atmospheric altitude in different latitudes and seasons, so the partitioning result cannot completely eliminate the influence of atmospheric transmission.
Monthly XCO2ano is calculated for each partition based on their respective monthly background region average XCO2 for 2010-2020. The formula of the area comparison method is as follows: where XCO is the XCO2 in the anthropogenic carbon emission area and XCO is the XCO2 in the background area. This study combines potential temperature data and ODIAC data to determine the background area. First, the ODIAC data is used to identify the background areas by screening out the areas without anthropogenic carbon emissions, and the remaining area is the anthropogenic carbon emission area. According to the monthly ODIAC data, the areas without anthropogenic carbon emissions are consistent in each month from 2010 to 2019, which are mainly distributed in the less-traveled areas of Tibet, Xinjiang, Inner Mongolia, and Northeast China. Second, China is divided into 5 regions by using potential temperature data. In this study, the XCO 2ano in different potential temperature regions was calculated according to the background regions of different potential temperature regions. In simulations with zonally uniform surface fluxes, XCO 2 is tightly correlated with potential temperature [43]. However, the influence of topographic and climatic factors may lead to differences in atmospheric transport modes in some regions within the same potential temperature isoline. Therefore, it is difficult to completely eliminate the influence of atmospheric transport. The average potential temperature data of atmospheric altitude of 1000 mb was selected for analysis in this study. However, there are differences in atmospheric altitude in different latitudes and seasons, so the partitioning result cannot completely eliminate the influence of atmospheric transmission.
Monthly XCO 2ano is calculated for each partition based on their respective monthly background region average XCO 2 for 2010-2020. The formula of the area comparison method is as follows: where XCO 2 emi is the XCO 2 in the anthropogenic carbon emission area and XCO 2 bck is the XCO 2 in the background area.
(2) Analysis of the change trend of the XCO 2ano In this study, the coefficient of variation (CV) was chosen to reflect the variation range of the XCO 2ano . The coefficient of variation was the absolute value reflecting the degree of dispersion of the data, which eliminated the influence of the data measurement scale and could more objectively reflect the variation in the XCO 2ano at different locations. The formula for the coefficient of variation is as follows: where σ is the standard deviation, and µ is the mean.
To address the deficiency that the CV in reflecting the changing trend, this study chooses the skewness coefficient (SKEW) to represent the changing trend of the XCO 2ano . The SKEW describes the asymmetry of the distribution in terms of a characteristic number of the degree of deviation. In this study, if the SKEW is greater than 0, the high value of the XCO 2ano is distributed in the first row of the time series. If it is less than 0, the high value of the XCO 2ano is distributed at the end of the time series. The smaller the value is, the closer the distribution of the high XCO 2ano value is to the end of the time series, which reflects the growth trend. The SKEW formula is as follows: where µ is the mean, σ is the standard deviation and M 0 is the mode.

Characteristics of Spatiotemporal Distribution in XCO 2ano
In this study, the spatial distribution of monthly average XCO 2ano values from 2010 to 2020 was calculated (Figure 4), and the results showed that the overall spatial distribution pattern of XCO 2ano values in China was similar to the spatial distribution pattern of fossil fuel emissions from ODIAC data. The high-value areas of China's XCO 2ano are mainly distributed in the Yangtze River Delta, the Pearl River Delta, and the Beijing-Tianjin-Hebei urban agglomeration. The areas around these three urban agglomerations also have high XCO 2ano values. The low-value areas of China's XCO 2ano are mainly distributed in western China. However, there is a high XCO 2ano in the economically underdeveloped Xinjiang. Some studies believe that Xinjiang's anthropogenic carbon emissions should be less, and the results are more uncertain due to the small number of satellite observations [31,37]. In fact, statistics show that energy consumption in Xinjiang is at a relatively high level [44,45]. Therefore, XCO 2ano high in Xinjiang is consistent with actual anthropogenic carbon emissions. The spatial distribution of monthly average XCO 2ano values from 2010 to 2020 was calculated with the background without the high altitudes is shown in Appendix A Figure A1.
The results show that there are obvious differences between XCO 2ano and ODIAC data in the Yangtze River Delta. Combined with land use data, it can be found that the differences are located in the southern area near urban land, mainly because local atmospheric transport may cause anthropogenic carbon emissions in the Yangtze River Delta region to be transported to the south, resulting in differences. There are also obvious differences between XCO 2ano and ODIAC data in the Pearl River Delta region and Northeast China. In addition to the diffusion of anthropogenic carbon emissions due to atmospheric transport, the rich vegetation resources around the Pearl River Delta and northeastern China will absorb a large amount of anthropogenic carbon emissions, resulting in differences.
To analyze the inter-annual variability of XCO 2ano in China, this study selected the XCO 2ano in 2010, 2015, and 2020 for spatial variability analysis ( Figure 5). The results show that the XCO 2ano in 2020 is significantly higher than that of the other two years. Although there were spatial differences in the XCO 2ano in the three years, the high values of XCO 2ano in the three years were mainly distributed in the southeast coastal areas, and  To analyze the inter-annual variability of XCO2ano in China, this study selected the XCO2ano in 2010, 2015, and 2020 for spatial variability analysis ( Figure 5). The results show that the XCO2ano in 2020 is significantly higher than that of the other two years. Although there were spatial differences in the XCO2ano in the three years, the high values of XCO2ano in the three years were mainly distributed in the southeast coastal areas, and the low values of XCO2ano were mainly distributed in the west and northeast regions. The spatial distribution results of XCO2ano for the remaining years in China from 2010 to 2020 are shown in Appendix A Figure A2. Although there are obvious changes in XCO2ano in China from 2010 to 2020, XCO2ano in coastal areas of China has always maintained a high level. The spatial distribution results of XCO2ano were calculated with the background without the high altitudes from 2010 to 2020 are shown in Appendix A Figure A3. In addition, to analyze the seasonality of the XCO 2ano in China, this study plotted the spatial distribution of the average XCO 2ano of four seasons (Spring: March-May; Summer: June to August; Autumn: September-November; Winter: December to February) from 2010 to 2020 ( Figure 6). The results show that the XCO 2ano in China has obvious seasonality. The XCO 2ano in winter is higher than that in the other three seasons, plant photosynthesis is the weakest in winter, and heating in winter increases the use of fossil fuels so that there is a higher XCO 2ano in winter. In summer, plant photosynthesis is the strongest, and the terrestrial ecosystem absorbs a large amount of the anthropogenic carbon emissions, so the XCO 2ano in summer is lower than that in the other three seasons, and the XCO 2ano in spring and autumn is between summer and winter. Although there were differences in the spatial distribution of XCO 2ano in the four seasons, they all had similar spatial distribution patterns. Due to the influence of vegetation photosynthesis and other factors, the seasonal variation of XCO 2ano cannot show the same seasonal variation of anthropogenic carbon emissions. If carbon sink information is subsequently added, XCO 2ano can be better linked to anthropogenic carbon emissions. The spatial distribution of the average XCO 2ano of four seasons was calculated with the background without the high altitudes from 2010 to 2020 are shown in Appendix A Figure A4. In addition, to analyze the seasonality of the XCO2ano in China, this study plotted the spatial distribution of the average XCO2ano of four seasons (Spring: March-May; Summer: June to August; Autumn: September-November; Winter: December to February) from 2010 to 2020 ( Figure 6). The results show that the XCO2ano in China has obvious seasonality. The XCO2ano in winter is higher than that in the other three seasons, plant photosynthesis is the weakest in winter, and heating in winter increases the use of fossil fuels so that there is a higher XCO2ano in winter. In summer, plant photosynthesis is the strongest, and the terrestrial ecosystem absorbs a large amount of the anthropogenic carbon emissions, so the XCO2ano in summer is lower than that in the other three seasons, and the XCO2ano in spring and autumn is between summer and winter. Although there were differences in the spatial distribution of XCO2ano in the four seasons, they all had similar spatial distribution patterns. Due to the influence of vegetation photosynthesis and other factors, the seasonal variation of XCO2ano cannot show the same seasonal variation of anthropogenic carbon emissions. If carbon sink information is subsequently added, XCO2ano can be better linked to anthropogenic carbon emissions. The spatial distribution of the average XCO2ano of four seasons was calculated with the background without the high altitudes To better analyze the seasonal variation in the XCO 2ano in China, this study calculated the monthly average monthly value of the XCO 2ano in different potential temperature zones (Figure 7). The results showed that the XCO 2ano in different potential temperature zones showed obvious seasonal changes, with the highest XCO 2ano in winter and the lowest XCO 2ano in summer (which even contained negative values). The negative value exists more in summer because vegetation photosynthesis strongly absorbs anthropogenic carbon emissions in summer. According to formula 1, the CO 2 absorbed by vegetation photosynthesis in the emission area is higher than that in the background area, which may lead to negative XCO 2ano . If carbon sink information can be added later, we can better eliminate the effect of vegetation photosynthesis. However, XCO 2ano in some years of area I showed an opposite time trend, such as 2012, 2013, and 2019, which showed negative peaks in winter. To better analyze the seasonal variation in the XCO2ano in China, this study calculated the monthly average monthly value of the XCO2ano in different potential temperature zones (Figure 7). The results showed that the XCO2ano in different potential temperature zones showed obvious seasonal changes, with the highest XCO2ano in winter and the lowest XCO2ano in summer (which even contained negative values). The negative value exists more in summer because vegetation photosynthesis strongly absorbs anthropogenic carbon emissions in summer. According to formula 1, the CO2 absorbed by vegetation photosynthesis in the emission area is higher than that in the background area, which may lead to negative XCO2ano. If carbon sink information can be added later, we can better eliminate the effect of vegetation photosynthesis. However, XCO2ano in some years of area I showed an opposite time trend, such as 2012, 2013, and 2019, which showed negative peaks in winter.

Spatiotemporal Variation of the XCO2ano
To analyze the spatiotemporal trends of the XCO2ano in China, this study calculated the coefficient of variation of the XCO2ano in China from 2010 to 2020 (Figure 8). The CV can reflect the magnitude of the change in XCO2ano. The results showed that the of XCO2ano in China changed greatly from 2010 to 2020, with an average CV of 36.16%. The average CV of area III and area IV was larger, 37.74% and 38.85%, respectively, while the average CV of area I was smaller, only 30.10%, and the average coefficients of variation of area II and area V were 35.06% and 34.48%, respectively. The average CV of the five regions was all higher than 30%, indicating that China's XCO2ano changed greatly during 2010-2020. The areas with large changes in XCO2ano in China are mainly distributed in the Yangtze River Delta urban agglomeration and western China located in area III and area IV, and the XCO2ano in northeastern China located in area I and Yunnan Province in area V has small changes. XCO2ano contains anthropogenic carbon emission signals, which are gener-

Spatiotemporal Variation of the XCO 2ano
To analyze the spatiotemporal trends of the XCO 2ano in China, this study calculated the coefficient of variation of the XCO 2ano in China from 2010 to 2020 (Figure 8). The CV can reflect the magnitude of the change in XCO 2ano . The results showed that the of XCO 2ano in China changed greatly from 2010 to 2020, with an average CV of 36.16%. The average CV of area III and area IV was larger, 37.74% and 38.85%, respectively, while the average CV of area I was smaller, only 30.10%, and the average coefficients of variation of area II and area V were 35.06% and 34.48%, respectively. The average CV of the five regions was all higher than 30%, indicating that China's XCO 2ano changed greatly during 2010-2020. The areas with large changes in XCO 2ano in China are mainly distributed in the Yangtze River Delta urban agglomeration and western China located in area III and area IV, and the XCO 2ano in northeastern China located in area I and Yunnan Province in area V has small changes. XCO 2ano contains anthropogenic carbon emission signals, which are generated by human production and life. Therefore, the change degree of XCO 2ano has a certain relationship with social and economic development. We can find that regions III and IV with high average CV include the Pearl River Delta and the Yangtze River Delta, which are the most economically active regions in China, while region I with low average CV is mainly located in northeast China, which has slow economic development. The coefficient of variation of the XCO 2ano from 2010 to 2020 was calculated with the background without the high altitudes are shown in Appendix A Figure A5.

Spatiotemporal Variation of the XCO2ano
To analyze the spatiotemporal trends of the XCO2ano in China, this study calculated the coefficient of variation of the XCO2ano in China from 2010 to 2020 (Figure 8). The CV can reflect the magnitude of the change in XCO2ano. The results showed that the of XCO2ano in China changed greatly from 2010 to 2020, with an average CV of 36.16%. The average CV of area III and area IV was larger, 37.74% and 38.85%, respectively, while the average CV of area I was smaller, only 30.10%, and the average coefficients of variation of area II and area V were 35.06% and 34.48%, respectively. The average CV of the five regions was all higher than 30%, indicating that China's XCO2ano changed greatly during 2010-2020. The areas with large changes in XCO2ano in China are mainly distributed in the Yangtze River Delta urban agglomeration and western China located in area III and area IV, and the XCO2ano in northeastern China located in area I and Yunnan Province in area V has small changes. XCO2ano contains anthropogenic carbon emission signals, which are generated by human production and life. Therefore, the change degree of XCO2ano has a certain relationship with social and economic development. We can find that regions III and IV with high average CV include the Pearl River Delta and the Yangtze River Delta, which are the most economically active regions in China, while region I with low average CV is mainly located in northeast China, which has slow economic development. The coefficient of variation of the XCO2ano from 2010 to 2020 was calculated with the background without the high altitudes are shown in Appendix A Figure A5.  In addition, this study calculated the SKEW of the XCO 2ano in China from 2010 to 2020 (Figure 9). SKEW can reflect the change direction of the XCO 2ano . The results show that the SKEW of the XCO 2ano in most areas of China from 2010 to 2020 is less than 0, and the average SKEW is −0.26, indicating that the high value of the XCO 2ano in most areas of China is biased toward the end of the time series, and the overall XCO 2ano in China is increasing. The average SKEWs of area I and area II are smaller, −0.53 and −0.52, respectively. The average SKEWs of area III, area IV, and area V are similar, −0.12, −0.11, and −0.19, respectively. The variation coefficient of the XCO 2ano in Northeast China located in area I is relatively small, indicating that although the variation in the XCO 2ano in Northeast China is small, it is in a significant increasing trend. In addition, there are many cases of positive SKEWs in China, indicating that in 2010-2020, the high value of the XCO 2ano in China is more inclined towards the start of the time series. For the region with significant growth of XCO 2ano , it is necessary to consider improving the regional low-carbon production capacity, adjusting the energy structure and improving energy utilization efficiency to control carbon emissions. The SKEW of the XCO 2ano from 2010 to 2020 was calculated with the background without the high altitudes are shown in Appendix A Figure A6.
In addition, the CV and SKEW of XCO 2ano in China also have seasonality. In this study, the average CV and average SKEW of the XCO 2ano in China in the four seasons were calculated ( Table 2). The results show that the CV in winter is the smallest, and the CV in summer is significantly larger than that in the other three seasons. Summer is the season where photosynthesis in plants is most obvious, and it is also a season with large changes in meteorological conditions. Therefore, there will be obvious differences in the meteorological conditions in the three months of summer, so that not only the internal summer XCO 2ano is different but also different months in summer in the same year, resulting in an excessively large CV in summer. Compared with summer, the effect of photosynthesis in winter is weaker, and the meteorological conditions are relatively more stable, so XCO 2ano in winter can better reflect anthropogenic carbon emissions. The SKEWs of the four seasons are all less than 0, among which the SKEW is the smallest in summer, indicating that the XCO 2ano does not only change greatly in summer, but also exhibits an increasing trend. Among them, the anthropogenic carbon emissions increased the fastest in summer and the slowest in spring. The average CV and average SKEW of the XCO 2ano in China in the four seasons were calculated with the background without the high altitudes and are shown in Appendix A Table A1. In addition, this study calculated the SKEW of the XCO2ano in China from 2010 to 2020 (Figure 9). SKEW can reflect the change direction of the XCO2ano. The results show that the SKEW of the XCO2ano in most areas of China from 2010 to 2020 is less than 0, and the average SKEW is −0.26, indicating that the high value of the XCO2ano in most areas of China is biased toward the end of the time series, and the overall XCO2ano in China is increasing. The average SKEWs of area I and area II are smaller, −0.53 and −0.52, respectively. The average SKEWs of area III, area IV, and area V are similar, −0.12, −0.11, and −0.19, respectively. The variation coefficient of the XCO2ano in Northeast China located in area I is relatively small, indicating that although the variation in the XCO2ano in Northeast China is small, it is in a significant increasing trend. In addition, there are many cases of positive SKEWs in China, indicating that in 2010-2020, the high value of the XCO2ano in China is more inclined towards the start of the time series. For the region with significant growth of XCO2ano, it is necessary to consider improving the regional low-carbon production capacity, adjusting the energy structure and improving energy utilization efficiency to control carbon emissions. The SKEW of the XCO2ano from 2010 to 2020 was calculated with the background without the high altitudes are shown in Appendix A Figure A6. In addition, the CV and SKEW of XCO2ano in China also have seasonality. In this study, the average CV and average SKEW of the XCO2ano in China in the four seasons were calculated ( Table 2). The results show that the CV in winter is the smallest, and the CV in summer is significantly larger than that in the other three seasons. Summer is the season where photosynthesis in plants is most obvious, and it is also a season with large changes in meteorological conditions. Therefore, there will be obvious differences in the meteorological conditions in the three months of summer, so that not only the internal summer XCO2ano is different but also different months in summer in the same year, resulting in an excessively large CV in summer. Compared with summer, the effect of photosynthesis in winter is weaker, and the meteorological conditions are relatively more stable, so XCO2ano in winter can better reflect anthropogenic carbon emissions. The SKEWs of the four seasons are all less than 0, among which the SKEW is the smallest in summer, indicating that the XCO2ano does not only change greatly in summer, but also exhibits an increasing trend. Among them, the anthropogenic carbon emissions increased the fastest in summer and the slowest in spring. The average CV and average SKEW of the XCO2ano in China in the four seasons were calculated with the background without the high altitudes and are shown in Appendix A Table A1.

Correlation Analysis of XCO 2ano
In this study, the regional comparison method was used to calculate the XCO 2ano in the emission area. The anthropogenic carbon emissions in the emission area have a cumulative effect on the XCO 2 , so XCO 2ano has a positive correlation with the number of residential areas ( Figure 10). The results show that there is a certain positive correlation between the XCO 2ano in China and the number of residential areas but the correlation is low. In addition, this study analyzes the correlation between the XCO 2ano and the number of residential areas in different potential temperature zones. Among them, in Area III, the two have the strongest positive correlation. Area III includes not only economically developed Yangtze River Delta and Beijing-Tianjin-Hebei urban agglomerations but also economically underdeveloped western China. The correlation analysis of average XCO 2ano and the number of residential areas for the other four areas is in Appendix A Figure A7.
The purpose of calculating the XCO 2ano through satellite observations is to monitor and track anthropogenic carbon emissions. In this study, the correlation between the average XCO 2ano and fossil fuel emissions was analyzed. The results show that there is only a weak positive correlation between the two (Figure 11). In Area V, the two have the strongest positive correlation, but the correlation coefficient r is only 0.51, mainly because there is a nonlinear relationship between the XCO 2ano and fossil fuel emissions, and the XCO 2ano is also affected by terrestrial ecosystems and atmospheric transport [35]. Therefore, the XCO 2ano does not only reflect anthropogenic carbon emissions, but is also affected by plant photosynthesis and wind transport. There are many cases where the XCO 2ano overestimates fossil fuel emissions. Atmospheric transport will transport some of the fossil fuel emissions to the surrounding areas, enhancing the XCO 2 observed by satellites in the surrounding areas. The underestimation is mainly due to the absorption of some anthropogenic carbon emissions by plant photosynthesis. The correlation analysis of average XCO 2ano and fossil fuel emissions for the other four areas is in Appendix A Figure A8.
In this study, the regional comparison method was used to calculate the XCO2ano in the emission area. The anthropogenic carbon emissions in the emission area have a cumulative effect on the XCO2, so XCO2ano has a positive correlation with the number of residential areas (Figure 10). The results show that there is a certain positive correlation between the XCO2ano in China and the number of residential areas but the correlation is low. In addition, this study analyzes the correlation between the XCO2ano and the number of residential areas in different potential temperature zones. Among them, in Area III, the two have the strongest positive correlation. Area III includes not only economically developed Yangtze River Delta and Beijing-Tianjin-Hebei urban agglomerations but also economically underdeveloped western China. The correlation analysis of average XCO2ano and the number of residential areas for the other four areas is in Appendix A Figure A7. The purpose of calculating the XCO2ano through satellite observations is to monitor and track anthropogenic carbon emissions. In this study, the correlation between the average XCO2ano and fossil fuel emissions was analyzed. The results show that there is only a weak positive correlation between the two (Figure 11). In Area V, the two have the strongest positive correlation, but the correlation coefficient r is only 0.51, mainly because there is a nonlinear relationship between the XCO2ano and fossil fuel emissions, and the XCO2ano is also affected by terrestrial ecosystems and atmospheric transport [35]. Therefore, the XCO2ano does not only reflect anthropogenic carbon emissions, but is also affected by plant photosynthesis and wind transport. There are many cases where the XCO2ano overestimates fossil fuel emissions. Atmospheric transport will transport some of the fossil fuel emissions to the surrounding areas, enhancing the XCO2 observed by satellites in the surrounding areas. The underestimation is mainly due to the absorption of some anthropogenic carbon emissions by plant photosynthesis. The correlation analysis of average XCO2ano and fossil fuel emissions for the other four areas is in Appendix A Figure A8.

Selection of Background Area
The XCO2 in the background area of different potential temperature partitions has some differences. According to the XCO2ano calculation method of Formula (1), the XCO2 in the background area will directly affect the XCO2ano result. The average XCO2ano in different regions is shown in Table 3. The results show that the average XCO2ano in Area I and Area II is lower, and the XCO2ano in the other three regions is significantly higher. Among them, Area IV has the highest XCO2ano. Area IV includes the Pearl River Delta urban agglomeration and the Yangtze River Delta urban agglomeration with active economic activities in China. Urban development needs to be driven by energy consumption, which produces anthropogenic carbon emissions that lead to an increase in atmospheric CO2. Relevant studies have shown that cities are the promoters of the increase in atmospheric CO2, and their anthropogenic carbon emissions also promote the growth of XCO2ano [35]. Partitioning according to potential temperature can reduce the influence of atmospheric transport, and there are obvious differences in the background XCO2 different par-

Selection of Background Area
The XCO 2 in the background area of different potential temperature partitions has some differences. According to the XCO 2ano calculation method of Formula (1), the XCO 2 in the background area will directly affect the XCO 2ano result. The average XCO 2ano in different regions is shown in Table 3. The results show that the average XCO 2ano in Area I and Area II is lower, and the XCO 2ano in the other three regions is significantly higher. Among them, Area IV has the highest XCO 2ano . Area IV includes the Pearl River Delta urban agglomeration and the Yangtze River Delta urban agglomeration with active economic activities in China. Urban development needs to be driven by energy consumption, which produces anthropogenic carbon emissions that lead to an increase in atmospheric CO 2 . Relevant studies have shown that cities are the promoters of the increase in atmospheric CO 2 , and their anthropogenic carbon emissions also promote the growth of XCO 2ano [35].
Partitioning according to potential temperature can reduce the influence of atmospheric transport, and there are obvious differences in the background XCO 2 different partitions, so it is necessary to calculate the XCO 2ano by partition. Some related studies have used the idea of zoning and selected the median or average value of XCO 2 in different zones as the background [31,35], while this study chose the average value as the background. To quantitatively analyze the impact of different background area selection methods on the capturing of anthropogenic carbon emission signals, this study selected the median XCO 2 in China as the background to calculate the XCO 2ano in the anthropogenic carbon emission area (Figure 12a). The results show that the XCO 2ano result cannot better reflect the spatial distribution pattern of China's anthropogenic carbon emissions. In addition, the median XCO 2 value in five zones was extracted according to the background area determined in this study and the XCO 2ano was calculated (Figure 12b). The results show that the XCO 2ano results are similar to the spatial distribution pattern of this study.
To quantify the potential of selecting median and average XCO 2 as background value to monitor anthropogenic carbon emissions, this study also performed a correlation analysis between the above XCO 2ano results and fossil fuel emissions ( Figure 13). The results show that the correlations between the XCO 2ano and fossil fuel emissions using the Chinese median XCO 2 and the regional median XCO 2ano as the background are 0.40 and 0.42, respectively. This is lower than the correlation coefficient between the XCO 2ano and fossil fuel emissions as determined in this study (Figure 11a). Therefore, choosing the average XCO 2 can improve the monitoring of anthropogenic carbon emissions. Most previous studies chose the median XCO 2 for subsequent XCO 2ano calculation. The selection of median can remove the influence of outliers, but for the small area of this study, it may cause information omission.

Uncertainty Factor Analysis of the XCO 2ano
Compared to bottom-up anthropogenic emissions inventory data, satellite observations are susceptible to meteorological, biological, and atmospheric conditions. The area explored in this study is large, and the wind field in a small area still interferes with the anthropogenic carbon emission signal in the XCO 2ano . Wind farms can diffuse CO 2 from anthropogenic sources to surrounding areas. XCO 2 over high altitude points has a less pronounced seasonal cycle. As shown in Figure 3b, the problem of higher XCO 2 values in the background area than in the emission area in summer largely disappears after removing the high altitude background points. (except for years 2010 and 2012). At present, WRF modelled wind field/wind measurements from LiDAR equipment has been widely used, among which Doppler wind LiDAR is a relatively new technology to acquire wind measurement [46]. In the future, relevant technologies can be considered to eliminate the influence of wind field. In addition, the XCO 2ano also includes the effect of biological sinks. Using only deserts or bare land as background to compare with vegetated emission areas causes a strong low bias in the emission effect because of removal of CO 2 by photosynthesis. We will work on this in a subsequent study. This is especially true in summer when the photosynthesis of plants absorbs a large amount of CO 2 and makes the correlation between XCO 2ano and fossil fuel emissions weaker. Subsequent auxiliary data related to CO 2 absorption and emissions should be added to enhance anthropogenic carbon emissions signals, such as net primary productivity and nighttime light images, to enhance the ability of satellite observations to assess and monitor anthropogenic carbon emissions. At the same time, satellites are sensitive to clouds and aerosols, resulting in relatively little available data in some areas [47]. Thus, satellite observations will have greater uncertainty [48].

Discussion on the Accurate Monitoring of Anthropogenic Carbon Emission
This study calculated the XCO 2ano in China based on satellite observation data, but there are some remaining shortcomings. The goal of the study is to monitor and evaluate anthropogenic carbon emissions. The results show that the XCO 2ano in this study has only a weak correlation with the fossil fuel emissions of the ODIAC data. The issue of how to solve the influence of biological and atmospheric factors is the main focus of subsequent research. Future research will consider adding data related to carbon absorption and carbon emissions such as vegetation index and solar-induced chlorophyll fluorescence. In addition, regional carbon neutrality is also an ecological issue of social concern, and subsequent studies will analyze the current status of regional carbon neutrality. To quantify the potential of selecting median and average XCO2 as background value to monitor anthropogenic carbon emissions, this study also performed a correlation analysis between the above XCO2ano results and fossil fuel emissions ( Figure 13). The results show that the correlations between the XCO2ano and fossil fuel emissions using the Chinese median XCO2 and the regional median XCO2ano as the background are 0.40 and 0.42, respectively. This is lower than the correlation coefficient between the XCO2ano and fossil fuel emissions as determined in this study (Figure 11a). Therefore, choosing the average XCO2 can improve the monitoring of anthropogenic carbon emissions. Most previous studies chose the median XCO2 for subsequent XCO2ano calculation. The selection of median can remove the influence of outliers, but for the small area of this study, it may cause information omission.
(a) (b) Figure 13. Correlation analysis of XCO2ano and fossil fuel emissions. (a,b) are the correlations between the XCO2ano and fossil fuel emissions using the Chinese median XCO2 and the regional median XCO2ano as the background value, respectively. To quantify the potential of selecting median and average XCO2 as background value to monitor anthropogenic carbon emissions, this study also performed a correlation analysis between the above XCO2ano results and fossil fuel emissions ( Figure 13). The results show that the correlations between the XCO2ano and fossil fuel emissions using the Chinese median XCO2 and the regional median XCO2ano as the background are 0.40 and 0.42, respectively. This is lower than the correlation coefficient between the XCO2ano and fossil fuel emissions as determined in this study (Figure 11a). Therefore, choosing the average XCO2 can improve the monitoring of anthropogenic carbon emissions. Most previous studies chose the median XCO2 for subsequent XCO2ano calculation. The selection of median can remove the influence of outliers, but for the small area of this study, it may cause information omission.
(a) (b) Figure 13. Correlation analysis of XCO2ano and fossil fuel emissions. (a,b) are the correlations between the XCO2ano and fossil fuel emissions using the Chinese median XCO2 and the regional median XCO2ano as the background value, respectively.

Uncertainty Factor Analysis of the XCO2ano
Compared to bottom-up anthropogenic emissions inventory data, satellite observations are susceptible to meteorological, biological, and atmospheric conditions. The area explored in this study is large, and the wind field in a small area still interferes with the Figure 13. Correlation analysis of XCO 2ano and fossil fuel emissions. (a,b) are the correlations between the XCO 2ano and fossil fuel emissions using the Chinese median XCO 2 and the regional median XCO 2ano as the background value, respectively.

Conclusions
This study aims to eliminate the influence of uncertainty factors in XCO 2 as far as possible to enhance the anthropogenic carbon emission signal. Based on the regional comparison method based on the idea of zoning, this study uses the potential temperature data and the ODIAC data set to effectively enhance the anthropogenic carbon emission signal in XCO 2 , and strengthen the remote sensing monitoring ability of anthropogenic carbon emission spatiotemporal changes. In this study, 2010-2020 monthly Mapping-XCO 2 data were used to calculate the spatiotemporal distributions of XCO 2ano in China, and the spatiotemporal changes in the results were analyzed. The XCO 2ano in China has obvious spatiotemporal differences. In addition, there are obvious seasonal variations, with the highest XCO 2ano in winter and the lowest in summer. The variation range of XCO 2ano in China is large, and showed an increasing trend, and the variation of XCO 2ano also has seasonality. The XCO 2ano has a high similarity with the spatial distribution of fossil fuel carbon emissions, which provides remote sensing observation means for anthropogenic carbon emissions monitoring. Compared with previous studies, the regional comparison method based on the idea of zoning can better enhance the anthropogenic carbon emission signal in XCO 2 , and using the average background regional as the background can better monitor anthropogenic carbon emission.
In order to explore the feasibility of satellite observation for assessing and monitoring anthropogenic carbon emissions, this study combined ODIAC data and potential temperature data to select a regionalized background region, and designed a regional comparison method based on the idea of zoning to eliminate background CO 2 and enhance anthropogenic carbon emission signals in XCO 2 . Relevant research results can provide a policy reference for China's "dual carbon" strategy. The anthropogenic carbon emission in the atmospheric carbon is very low, and the fluctuation caused by it is difficult to accurately measure. With the development of remote sensing technology, it is hoped that carbon monitoring satellites can provide higherprecision XCO 2 data. In the future, we will add CO 2 -related data for background region selection and supplement carbon neutrality research.

Data Availability Statement:
No new data were created or analyzed in this study. Data sharing is not applicable to this article.               There is a high correlation between the FFCO2 emissions and the average XCO2ano in areas II (b), III (c), and IV (d), while the lowest correlation between the two is found in area I (a). Figure A8. The correlation analysis between XCO 2ano and fossil fuel emissions in different area. There is a high correlation between the FFCO 2 emissions and the average XCO 2ano in areas II (b), III (c), and IV (d), while the lowest correlation between the two is found in area I (a).