A city-level comparison of fossil-fuel and industry processes-induced CO2 emissions over the Beijing-Tianjin-Hebei region from eight emission inventories
Carbon Balance and Management volume 15, Article number: 25 (2020)
Quantifying CO2 emissions from cities is of great importance because cities contribute more than 70% of the global total CO2 emissions. As the largest urbanized megalopolis region in northern China, the Beijing-Tianjin-Hebei (Jing-Jin-Ji, JJJ) region (population: 112.7 million) is under considerable pressure to reduce carbon emissions. Despite the several emission inventories covering the JJJ region, a comprehensive evaluation of the CO2 emissions at the prefectural city scale in JJJ is still limited, and this information is crucial to implementing mitigation strategies.
Here, we collected and analyzed 8 published emission inventories to assess the emissions and uncertainty at the JJJ city level. The results showed that a large discrepancy existed in the JJJ emissions among downscaled country-level emission inventories, with total emissions ranging from 657 to 1132 Mt CO2 (or 849 ± 214 for mean ± standard deviation (SD)) in 2012, while emission estimates based on provincial-level data estimated emissions to be 1038 and 1056 Mt. Compared to the mean emissions of city-data-based inventories (989 Mt), provincial-data-based inventories were 6% higher, and national-data-based inventories were 14% lower. Emissions from national-data-based inventories were 53–75% lower in the high-emitting industrial cities of Tangshan and Handan, while they were 47–160% higher in Beijing and Tianjin than those from city-data-based inventories. Spatially, the emissions pattern was consistent with the distribution of urban areas, and urban emissions in Beijing contributed 50–70% of the total emissions. Higher emissions from Beijing and Tianjin resulted in lower estimates of prefectural cities in Hebei for some national inventories.
National-level data-based emission inventories produce large differences in JJJ prefectural city-level emission estimates. The city-level statistics data-based inventories produced more consistent estimates. The consistent spatial distribution patterns recognized by these inventories (such as high emissions in southern Beijing, central Tianjin and Tangshan) potentially indicate areas with robust emission estimates. This result could be useful in the efficient deployment of monitoring instruments, and if proven by such measurements, it will increase our confidence in inventories and provide support for policy makers trying to reduce emissions in key regions.
Cities play a significant role in global greenhouse gas emissions, especially in urban areas, which are responsible for 67–76% of the global CO2 emissions and energy consumption . Cities have become the critical and basic units for implementing emissions mitigation policies [2,3,4,5]. However, city-level mitigation actions remain daunting challenges [1, 6, 7]. City carbon emissions are influenced by the physical environment, economic development, urbanized density, industry structure, and energy use patterns specific to each city [1, 8]. Cities with a heavy industry, high traffic load, and high population density more easily have high emissions . Discrepancies in the emissions and emission-socioeconomic characteristics among different cities require the development of corresponding policies . Moreover, how to deploy observational instruments to form an efficient network is rather challenging [9,10,11], especially when there is no robust understanding of emissions. With emissions patterns identified by multiple inventories, such deployments would have more scientific guidelines. Furthermore, city strategies that reduce carbon emissions are expected to achieve emissions mitigation and meet a city’s economic growth goals . Therefore, an accurate understanding of city-level CO2 emissions is of great importance in developing and implementing mitigation strategies.
With rapid economic development and urbanization in China, cities account for 85% of China’s CO2 emissions . Emission error could be much larger at subnational levels . However, most of the existing studies focus on national [15,16,17,18,19], provincial [4, 16, 20, 21], or sectoral CO2 emissions inventories in China [22,23,24,25]. For example, using nine inventories, Han et al.  estimated that the national total fossil fuel and industrial process-related CO2 emissions were 9.8 (9.2–10.4) Gt CO2 in 2016, and the emissions estimated from provincial-data-based inventories were more consistent than those from spatial disaggregation of national energy statistics. A few efforts have been made to estimate city-level emissions, but these efforts have mainly focused on megacities or provincial capital cities due to the limited energy data [26,27,28,29]. To date, the studies by Shan et al.  and Cai et al.  included 182 and 305 Chinese cities, respectively, and still had gaps in city coverage. Zheng et al.  estimated all the cities’ CO2 emissions intensities in mainland China, yet the latest year included in the study was 2013, and the study lacked temporal dynamics. Thus, a more comprehensive assessment of city-level CO2 emissions in China is critical for understanding the role of cities in carbon emissions.
The Beijing-Tianjin-Hebei (Jing-Jin-Ji, JJJ) region is the largest urbanized megalopolis region in northern China, covering an area of 218,000 km2 and home to 112.70 million people [32,33,34]. Cities in the JJJ region include the municipalities Beijing and Tianjin and eleven prefecture cities in Hebei Province. In 2018, the total energy consumption in the JJJ region accounted for more than 10% of China’s total energy consumption . Moreover, coal is the primary energy source in this region [33, 34]. The JJJ region is under considerable pressure to reduce CO2 as well as air pollutant emissions [31, 35]. Beijing and Tianjin have committed to peak their CO2 emissions by 2020 and approximately 2025, respectively, in the 13th Five-Year Plan . Hebei Province is experiencing rapid industrial and urban development , contributing greatly to the national CO2 emissions . However, there is a wide range of CO2 emission estimates in the JJJ region, especially in the areas with high emissions. The emissions from Beijing and Tianjin estimated by Cai et al.  and Wang  were 41–57% higher than those estimated by Shan et al. . In addition, based on the results of Mi et al.  and Cai et al. , Tangshan’s carbon emissions differed by 50%. To our knowledge, there is not a comprehensive assessment on prefectural city level CO2 emissions in the JJJ region. Cai et al.  reported the provincial emissions, but not covered the prefectural cities in JJJ. Therefore, an accurate estimation of the CO2 emissions in the JJJ region is of great significance in terms of providing accurate information for developing mitigation policies.
The assessment of existing CO2 emissions inventories is urgently needed, yet direct observations of CO2 emissions at the city scale is limited, especially in developing countries [10, 11]. Here, we conducted a comprehensive analysis of 8 state-of-the-art inventories and presented the temporal dynamics, spatial distributions, and urban and non-urban fractions of 13 cities. We recognized the similarities and differences in emissions and thus improved the understanding of current inventories; this research provides useful information for policy making related to reducing city emissions and monitoring CO2.
Data and methods
We used annual CO2 emissions data from eight emission inventories, including the China High Resolution Emission Database (CHRED); China Emission Accounts and Datasets (CEADs); Multi-resolution Emission Inventory for China (MEIC), version 1.3; the Nanjing University CO2 emission inventory (NJU); the Peking University CO2 emission inventory (PKU), version 2 (PKU-CO2-V2); the Open-source Data Inventory for Anthropogenic CO2, version 2018 (ODIAC2018); the Emissions Database for Global Atmospheric Research, version 5.0 (EDGARv5.0); and the Fossil Fuel Data Assimilation System, version 2.2 (FFDAS v2.2). Below, the eight inventories were categorized into two based on the emission calculation methodologies (Table 1).
Emission inventories based on city-level data
We used two inventories that are based on city-level data. CHRED was constructed by enterprise-level point-source data and China’s city and provincial statistics, and then carbon emissions were allocated to 10 km resolution using proxy data in 2007 and 2012 [15, 30]. The CEADs inventory provides total CO2 emission estimates at the provincial and city levels from 2000 to 2016 based on apparent energy consumption data and local optimized emission factors [4, 17, 49].
Emission inventories based on provincial or national level data
We used two inventories that are based on province-level data. China’s carbon emissions from MEIC are developed by a technology-based methodology based on provincial energy consumption, combustion/industrial/control technologies and emission factor databases covering 2000 to 2016 at 0.25, 0.5, and 1 degree spatial resolutions by Tsinghua University [18, 31, 50]. NJU calculated China’s CO2 emissions using provincial energy statistics and spatially distributed emissions based on the location of large point sources (power plants and cement plants) and various proxy data, i.e., using GDP for industry-related emissions and population for transportation and other emissions, at 0.25 degree resolution from 2000 to 2016 [5, 18].
We used 4 inventories that are based on national-level emission estimates. These subnational emission distributions are largely achieved by emission downscaling. PKU distributes national or subnational fuel data with various proxies (e.g., power plants as point sources, night-time light to distribute national gas flaring and population for others) based on the subnational disaggregation method at 0.1 degree resolution from 1960 to 2014 . The year 2018 version of the ODIAC emissions data product (ODIAC2018) distributes national emissions into 1 km and 1 degree grids from 2000 to 2017 based on spatial proxies, such as geographical locations of power plant emissions, satellite observations of nightlights to distribute nonpoint emissions, and aircraft and ship fleet tracks [52, 53]. EDGAR provides carbon emissions on the 0.1 degree grid from 1970 to 2018 based on national emissions by a variety of spatial proxy data, including power plants from CARMA3.0 for point sources, road network and different weighting factors for line sources and population for residential and commercial emissions [54, 55]. FFDAS quantifies carbon emissions using a data assimilation technique incorporating remote-sensing and national fuel accounts and power plants at a 0.1 degree resolution from 1997 to 2015 [42, 56]. All data sets used here have yearly data. For more details, please refer to Han et al. .
These inventories were first extracted by a JJJ mask (in shapefile format) from the National Geomatics Center of China using ArcGIS 10.02 software (ESRI, 2012). The emissions from urban and non-urban areas were separated by using an urban mask from the European Space Agency (ESA) Climate Change Initiative (CCI) land cover maps with a 300 m resolution (https://www.esa-landcover-cci.org), and the urban area here mainly refers to impermeable surfaces, with high coherence (stability) and bright reflections, maintained in time and under varying angles detected by satellite. This means there is a density of human structures such as houses, commercial buildings, roads, bridges, and railways, while a city refers to prefectural level zones with a legal definition that defines a physical geographic boundary. ArcGIS was used to obtain prefectural cities’ total emissions in urban and non-urban areas using city mask data from the National Geomatics Center of China. The data pre-process procedures are as follows: (1) Convert the grid data to polygon data, which will keep the original value of the grid; (2) Split the polygon data using the prefectural city boundaries, and this step produces the real areas of a polygon within the city boundary; (3) Calculate the actual areas of each polygon using “Calculate geometry” and multiply the areas with emissions intensity (e.g. kg CO2/km2); (4) Finally, sum all the values within the boundary of a city. Emission intensity was calculated as the CO2 emissions divided by the GDP, data which were derived from the National Bureau of Statistics of the People's Republic of China (NBS) and websites of prefectural cities’ statistics bureaus in Hebei. Linear regressions were conducted using the Python scipy package between inventories with CEADs or CHRED.
City-level CO2 emissions
In 2012, the emissions of 13 cities varied widely from 13 to 282 Mt CO2 (or 72 ± 45, mean ± SD) (Fig. 1). Tangshan, Tianjin, Handan and Beijing are high-emission areas. The total emissions estimated from provincial-data-based inventories (i.e., MEIC and NJU) were 6% higher than those from city-data-based inventories (i.e., CHRED and CEADs) but were 14% lower from downscaled national-level emissions (i.e., PKU, ODIAC, EDGAR, and FFDAS) in 2012. There was a great discrepancy in national-data-based inventories with a range from 657 to 1132 Mt CO2 (or 849 ± 214, mean ± SD). The CO2 emissions from EDGAR and PKU were 28% and 34% lower than the average from city-data-based inventories. The 13 cities have substantial differences in natural resources and socioeconomic conditions. A city with high productivity, rapid economic growth, and a large population tends to have high carbon emissions. Emissions from metropolises (e.g., Beijing and Tianjin) with advanced economics and high urbanization rates greatly contributed to the JJJ total emissions (23–43%, Fig. 1). In Hebei Province, high emissions were also located in the provincial capital and in industrial cities, such as Shijiazhuang, Tangshan, Handan, and Baoding, which accounted for 57%–68% of the total emissions in Hebei. The 7 remaining cities accounted for 32–43% of the total emissions in Hebei.
Since CHRED and CEADs use city-level statistical data and these two datasets are more consistent with each other than with the other datasets, they are used as references for national and provincial data-based inventories. The emissions of Beijing from EDGAR, MEIC, and NJU are comparable to those from CHRED, with differences ranging from 1 to 8 (1–7%) Mt CO2 in 2012. However, compared to the values estimated by CHRED, the emissions were 115% and 32% higher when estimated by ODIAC and FFDAS, respectively, and 31% lower when estimated by PKU. These trends were also found in Tianjin, and the emissions from ODIAC were equal to 227 Mt CO2, a value that was 47% higher than that from CHRED but 53% lower than that from PKU. Emissions in prefectural cities with coal mines and heavy-intensity industries also showed large differences, such as Tangshan and Handan. Tangshan’s emissions varied largely across these datasets. CHRED and CEADs both produced larger estimates (182 and 282 Mt), while estimates from PKU, EDGAR, and FFDAS were only 69–87 Mt, or 52–62% lower, than that by CHRED. Handan’s emissions were up to 121 and 130 Mt based on CHRED and CEADs, respectively, but ranged from 61 to 89 Mt (26–53% lower than CHRED and CEADs) for the other inventories. It should also be noted that CHRED and CEADs produced large differences in some cities, such as Tangshan, Shijiazhuang, and Cangzhou, and these differences highlighted future directions for city-level inventories in JJJ. According to the plan of “Collaborative Development of Beijing, Tianjin and Hebei Province”, Zhangjiakou and Chengde belong to ecological conservation areas; thus, heavy industry is not recommended, and more efforts should be concentrated in high-emission cities such as Tangshan, Shijiazhuang and Handan.
Spatial pattern of CO2 emissions
The uneven spatial distribution of CO2 emissions reflects the highly diverse conditions in population, economic development, and natural environment in the JJJ region. The spatial distributions showed reasonably good agreement in patterns such as the high emissions areas of Beijing-Tianjin-Tangshan, although they varied in detail (Fig. 2). Low emissions (e.g., < 100 ton CO2 km−2) from different datasets were mainly located in the northwestern part of the region. This result is because cities in the north, such as Zhangjiakou and Chengde, have a low population density and have lower economic development . High-emission areas (> 10,000 t CO2 km−2) are clustered in urban centers in the south and east. Hotspots of CO2 emissions (> 50,000 t CO2 km−2) recognized by most inventories are located in the urban areas of Beijing, Tianjin, and Tangshan.
Specifically, the spatial distributions of CO2 emissions from EDGAR and FFDAS are consistent with the pattern from CHRED. However, the total amount of emissions in the region from EDGAR is 23% lower than that from CHRED. This difference is caused by the discrepancies estimated in Handan, Tianjin, and Tangshan from EDGAR, which is altogether approximately 204 Mt CO2 lower than that from CHRED. Although the difference in total emissions between FFDAS and CHRED is less than 3%, the emissions in Tangshan and Beijing from FFDAS are 62% lower and 32% higher than those from CHRED. ODIAC distributes the non-power plant portion of national emissions based on satellite nighttime light [52, 53]. The pattern of carbon emissions in ODIAC is highly consistent with the spatial distribution of urban areas. Low emissions ranging from 50 to 500 t CO2 km−2 are not found in ODIAC. This result is probably due to the low emissions mostly located in the non-urban areas, which are mainly covered by vegetation and do not emit strong night light (Fig. 2). ODIAC and FFDAS have similar patterns, shown in Fig. 2, which can be because they both used night-time light to distribute nonpoint emissions. Compared with CHRED, the carbon emissions for the remaining inventories all showed higher emissions in the southeastern areas, which are mainly covered by croplands, with the results of MEIC, NJU, PKU and ODIAC being more notable and the results of EDGAR and FFDAS being relatively weaker. For the spatial allocation of emissions, CHRED, EDGAR and FFDAS give more emissions to high emitting grids, with > 50,000 t CO2/km2 grids contributing more than 50% of the total emissions (Fig. 2c), while for NJU and PKU, the same level of high emitting grids contributed less than 10%.
It should be noted that the urban area contributed largely to the total emissions because of its large energy consumption and population . Here, according to the urban extent from the ESA CCI land cover maps with a 300 m resolution (https://www.esa-landcover-cci.org) (Fig. 2), we extracted urban emissions from the seven gridded inventories, namely, CHRED, NJU, MEIC, PKU, ODIAC, EDGAR and FFDAS. As illustrated in Fig. 3, the proportion of urban emissions shows obvious differences among datasets, ranging from 17% (NJU) to 50% (ODIAC). As the largest contributors, the proportions of urban emissions in Beijing and Tianjin varied between 37% (NJU) –73% (ODIAC) and 23% (NJU) –63% (ODIAC), respectively. ODIAC tended to overestimate emissions in urban areas, resulting in the strongest urban-non-urban emissions gradients, especially in Tianjin (63% from urban) and Hengshui (79% from urban). This result is due mainly to the use of the nightlight proxy [14, 58]. The urban emissions from FFDAS also show a similar pattern with those from ODIAC, except for Zhangjiakou and Hengshui. In addition, in high-emitting cities in Hebei, such as Tangshan, all inventories identified urban emissions accounting for a smaller fraction (20%–40%) than those in Beijing and Tianjin.
Temporal variation in city-level CO2 emissions and emissions intensity
The interannual variations in carbon emissions from existing inventories all showed an increasing trend during 2000–2012, ranging from 355 ± 58 Mt to 915 ± 178 Mt CO2, and then tended to level off or showed a slight downward trend afterwards (Fig. 4). Total emissions of the JJJ region from ODIAC, CEADs and EDGAR increased faster than others, with average growth rates of 9.7%, 9.2% and 9.1%, respectively, during 2000–2012. However, the regional total emissions from FFDAS, PKU, MEIC, and NJU showed a relatively small trend, with an annual growth rate of 6.0–8.4%. PKU and EDGAR tended to underestimate emissions compared to CHRED and CEADs, especially in high-emitting cities. The interannual variation at the city scale from MEIC was consistent with that from NJU, both of which were based on provincial statistical data.
Specifically, for Beijing, the CO2 emissions from MEIC, NJU, CEADs, PKU, and CHRED tended to be stable since 2007, which was consistent with the studies by Li et al.  and Shan et al. . This result may be because Beijing decreased coal use by 43% and increased natural gas consumption by 144% from 2007 to 2014 [4, 16]. However, the annual emissions in Beijing from ODIAC fluctuated greatly, ranging from 73 Mt CO2 in 2000 to 263 Mt CO2 in 2017. PKU is a downscaled inventory; however, it uses provincial consumption fractions to rescale the IEA (International Energy Agency) total fossil fuel consumption when distributing emissions to grids, and thus, it can capture the Beijing decreasing trend after 2007, while other national downscaled inventories cannot. The emissions in Tianjin and Tangshan experienced rapid growth and then outpaced the emissions of other cities in recent years. The time series of emissions in Tianjin from ODIAC, MEIC, NJU, CEADs, and CHRED all showed an apparent growing trend during 2009–2014, which mainly resulted from the increment of coal use (22%) and crude oil (90%) . The decline between 2005 and 2010 for CEADs was because of the change in statistical methods and dimensions, and blast furnace gas was lacking from the energy balance sheet during this period. For other periods, high emissions and large fluctuations in Tangshan’s emissions from CEADs were partially caused by its heavy industrial system [12, 39], whose iron and steel production accounted for 60–70% of Hebei’s total production.
The emissions intensity (CO2 emissions per unit of GDP) in the JJJ region showed a decreasing trend since 2000 (Fig. 4). Interannual changes in emissions intensity among datasets were consistent with the total emissions. Moreover, the decoupled relationship between an increase in GDP and a decrease in emissions intensity indicated that a reduction in carbon emissions intensity could be achieved while also maintaining economic growth. For example, the GDP in Beijing and Tianjin significantly increased by average annual rates of 14% and 16%, respectively, during 2000–2016, but emissions intensity declined by 8–11% annually from PKU, MEIC, NJU, and CEADs and by 5–7% from ODIAC, EDGAR, and FFDAS.
Comparison of city-level CO2 emissions from different inventories to CHRED and CEADs
CHRED and CEADs were used as references to evaluate other datasets because they are both based on city-level energy statistics and consistent with each other. Furthermore, CHRED included unique comprehensive point sources (over 1.5 million enterprises) , and CEADs used measured local emission factors from 602 coal samples and 4243 coal mines . As expected, the city-level CO2 emissions showed the best correlation between CEADs and CHRED, with a correlation coefficient (R) as high as 0.9 (Figs. 5a and 6a) and with a slope close to 1 and a smaller intercept than the others. The emissions from MEIC and NJU were both highly correlated with those from CEADs and CHRED, with R values ranging from 0.8 to 0.9 and slopes ranging from 0.5 to 0.7. These results were probably due to MEIC and NJU using provincial energy statistics as input data. CHRED and CEADs also used provincial inventories to calculate some city-level CO2 emissions because of the lack of consistent and accurate energy statistics in certain cities [12, 30].
PKU was strongly correlated with CEADs and CHRED, with R values both equal to 0.8. The slopes were only 0.2 and 0.4 for the correlation relationships of PKU and CEADs and of PKU and CHRED, respectively (Figs. 5 and 6). The lower emissions estimated from PKU were mainly due to the rescaled energy data for China’s total by IEA statistics. Emissions from EDGAR had a relatively weak relationship with those of CEADs and CHRED, with R values of 0.6 and 0.8, respectively. ODIAC and FFDAS showed the lowest relationships with CEADs and CHRED, with R values of 0.5 and 0.7, respectively. These results were probably due to ODIAC, EDGAR and FFDAS all disaggregating emissions based on national energy statistics.
Differences in prefectural-city-, provincial- and national-statistical data-based inventories
The differences between provincial- and national-statistical data-based emissions have been well discussed in previous studies [19, 62, 63]. However, the differences between provincial and prefectural city statistical data-based emissions are poorly understood. Thus, it is difficult to conclude which is closer to the truth. Of these inventories, CEADs provided both prefectural-city- and provincial-statistical data-based estimates. We found that prefectural-city-statistical data-based estimates were 60% higher than provincial-statistical data-based estimates (Fig. 7) from 2000–2010 and gradually decreased to approximately zero after 2012, i.e., the sum of prefectural-statistical data became more consistent with provincial data in recent years due to the wider coverage of large-volume industries in local statistical authorities, indicating the improvement of statistical data between these two levels.
At the prefectural city level for JJJ and other cities, the estimated differences generally increased from less than 10% to more than 300% when more inventories were included (Table 2) [5, 64,65,66]. Generally, the prefectural-city-statistical data-based inventories of CHRED and CEADs were more consistent than the others. For Hebei prefectural cities (such as Tangshan and Handan), the spread was relatively large among the three methods (i.e., the national, provincial and prefectural statistics data-based inventories), with differences ranging from 69 to 282 Mt CO2 for Tangshan and from 61 to 130 Mt CO2 for Handan in 2012 (Table 2). For Tangshan, the national inventories (ODIAC, EDGAR, PKU and FFDAS) were 69–117 Mt CO2, which were all smaller than the prefectural-city level and the provincial estimates of 138–282 Mt CO2 (MEIC, CHRED, MEIC and NJU). CEADs estimated Tangshan’s emissions to be much higher than the other inventories due to the large consumption of coke, e.g., the amount for 2015 was 31,986,000 tons, which was 41.4% that of Hebei Province . The estimates for provincial-statistic data-based inventories were lower in Tangshan and Handan and thus tended to calculate higher emissions for other cities such as Baoding, Langfang and Hengshui due to the constraint of provincial total emissions. Moreover, the estimates for Beijing and Tianjin were more consistent for provincial estimates than the national estimates.
Urban and non-urban CO2 emissions and implications for carbon monitoring instrument distribution
The urban extent of the ESACCI-LC product was 17%, 14%, and 5% of the land area for Beijing, Tianjin and Hebei, respectively, and the corresponding direct CO2 emissions were 64%, 37%, and 20%, respectively, of the total emissions for Beijing, Tianjin and Hebei for the mean of all inventories, which indicates that high emissions are spatially located more often in the urban areas of Beijing and Tianjin and are more diffusive in non-urban areas in Hebei. Similarly, an urban CO2 emissions study conducted by Cai et al.  showed that the urban extents were 17%, 17%, and 4% of land area in Beijing, Tianjin, and Hebei, which contributed 84%, 60%, and 41% of direct CO2 emissions, respectively. The differences between these two studies were mainly due to the inventories and urban land masks used. The ESA 300 m data showed more urban details, while Cai et al.  used a homemade dataset based mainly on the county/district and town/township GIS data, which are more continuous in space, to depict the urban areas. However, our study and Cai et al.  are largely consistent in their urban extent estimates for Beijing, Tianjin and Hebei and were approximately 20% lower than the CO2 emissions for urban areas calculated by Cai et al. . Specifically, at the prefectural city level, emissions from national-data-based inventories were 53–75% lower in the industrial cities of Tangshan and Handan and 47–160% higher in Beijing and Tianjin than those from city-data-based inventories, and this difference was consistent with Gately and Hutyra . The implication to the community and society is that national-based inventories have more biases in prefectural city levels than in local data-based inventories, and thus, we should be cautious when using national-based inventories for city-level use, such as in emissions evaluation, modeling research and policy making.
The high-emission areas recognized by most of the present inventories have significant implications for monitoring instrument deployment [7, 10, 11, 67], although this type of comparison does not allow us to further discuss the accuracy of emission estimates. Thus, this information needs to be supplemented by objective physical measurements to validate the accuracy (e.g., [14, 53]). The areas with good agreement among inventories can be key areas for observation deployment, and in return, we can use these measurements to validate the accuracy of inventories. For example, third-party monitoring of CO2 emissions using high-density low-cost sensor networks is becoming possible due to the development of nondispersive infrared (NDIR) technology. Martin et al.  investigated a low-cost NDIR sensor and compared it with the standard instrument Los Gatos and found that the accuracy could reach 2–5 ppm after environmental factor corrections. More than 300 of such CO2 sensors were deployed in a network in Switzerland, and they were able to resolve CO2 changes and differences with magnitudes larger than ~ 20 ppm . Bao et al.  proved such sensors to be promising in high-emission areas near Shijiazhuang, Hebei Province. However, how to deploy such nodes to form an efficient network has been a challenge . With emissions patterns identified by most inventories, such deployment would have more scientific references. For example, in Beijing and Tangshan, emissions from urban areas consisted of 50–70% and 20–30%, respectively, which indicated that more sensors should be deployed in the urban areas of Beijing, while more nodes should be deployed in the non-urban areas of Tangshan. Moreover, atmospheric CO2 measurements need to consider the combined effects of biospheric and anthropogenic signals since the JJJ region is surrounded by vegetation, especially in the western and northern parts, and atmospheric transport also plays a significant role in the measurements.
The identified high-emission areas might also indicate potential target areas for emission reductions for policy-makers if such high-emission areas are also confirmed by instruments. Beijing and Tianjin have committed to peak their CO2 emissions by 2020 and approximately 2025, respectively, in the 13th Five-Year Plan . More specifically, local governments of Beijing have proposed a clear peak time of total and per capita CO2 emissions in 2020 in its "13th Five-Year Plan" for energy conservation and consumption reduction and climate change; thus, the city must promote the revolution in energy production and consumption, improve energy efficiency, and accelerate the construction of low-carbon cities . High-emission areas (e.g., > 50,000 ton emissions in Fig. 2) need to allocate high priority to emissions control. The joint emissions control of JJJ also needs to identify the high-emission areas in Fig. 2.
Point source contributions
Point emissions consist of a large proportion of total emissions [15, 33]. For the JJJ regional total, CHRED estimated the highest point emissions proportion of 78.1%, including 12,991 industrial key emission sources and 187 industrial process sites . Other inventories produced much smaller proportions ranging from 19% for NJU to 43% for EDGAR (Fig. 8), and MEIC and EDGAR both had more than 100 large emission grids associated with power plants. Point sources for CHRED consisted of power plants, industries and industrial processes, while other inventories (e.g., PKU, ODIAC) mostly used the CARMA dataset, which included only power plants. NJU included power plants ranking with the top 80% in electricity production and cement production that exceeded 1 Mt yr−1 , and there were only 42 power plants and 23 cement plants and thus fewer high-emission grids (N = 37, Fig. 8), which may be the reason why NJU point emission fractions were much smaller than those of other inventories. Since the gridded maps are mostly produced based on point sources, line sources (transportation emissions) and finally area sources, and area sources are distributed using the total estimate minus the point and line sources and proxies (e.g., night light, population, GDP), point source numbers, geolocations and emission magnitudes determined a very large degree of the gridded products (e.g., ). To improve the accuracy in mapping emissions and mitigate the errors in emission estimates, it is important to include reliable information, preferably reported information, on large point sources as much as possible. Thus, it is preferable for the community to share the point source information and for the information to be accurate to improve the understanding of the point source emissions (e.g., ). This is also important for emission monitoring purposes.
Here, we conducted a comprehensive analysis of city-level fossil fuel and industrial process-related CO2 emissions from cities in the Beijing-Tianjin-Hebei region. We showed their temporal dynamics, spatial distributions, and urban and non-urban emissions fractions. We recognized the similarities and differences in emissions and thus improved the understanding of current inventories and provided useful information for policy making in terms of reducing city emissions and monitoring CO2. The results showed that compared to city-data-based inventories, provincial-data-based inventories were 6% higher, and national-data-based inventories were 14% lower in 2012. Compared with city-data-based inventories, the lower estimates (53–75%) in the industrial cities of Tangshan and Handan resulted in higher estimates (47–160%) in Beijing and Tianjin for national-data-based inventories. Due to the more complete data of industrial enterprises above the state designated scale, the differences between city statistical data-based estimates and provincial statistical data-based estimates decreased from 60% to approximately zero from 2000 to 2012, indicating the improvement in local statistical authorities. Spatially, all datasets agreed with high emissions in the triangular spatial distribution pattern of Beijing-Tianjin-Tangshan and low emissions in the northern parts of Zhangjiakou and Chengde. The implications of the consistent spatial distribution patterns recognized by these inventories provide useful information for the efficient deployment of monitoring instruments, and in return, the independent measurements from these areas will increase our confidence in inventories and thus provide support for policy makers in joint emissions reductions.
Availability of data and materials
The data sets of ODIAC, EDGAR, PKU, CEADs and FFDAS are freely available from http://db.cger.nies.go.jp/dataset/ODIAC/DL_odiac2018.html, https://edgar.jrc.ec.europa.eu/overview.php?v=50_GHG, http://inventory.pku.edu.cn/download/download.htmlhttp://www.ceads.net/, and http://ffdas.rc.nau.edu/Data.html, respectively. CHRED, MEIC and NJU are available from the data developers upon request.
Open-source Data Inventory for Anthropogenic CO2
Emissions Database for Global Atmospheric Research
Peking University CO2 emission inventory
Fossil Fuel Data Assimilation System
China Emission Accounts and Datasets
China High Resolution Emission Database
Multi-resolution Emission Inventory for China
Nanjing University CO2 emission inventory
CARbon Monitoring for Action
The First China Pollution Source Census
China Power Emissions Database
Commission for Environmental Cooperation
China Cement Almanac
China Cement Industry Enterprise Indirectory
Gross domestic product
European Space Agency
Climate Change Initiative
National Bureau of Statistics
IPCC. IPCC, 2014: Summary for Policymakers. In: Climate Change 2014: Mitigation of Climate Change. Contribution of Working Group III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [Edenhofer O, Pichs-Madruga R, Sokona Y, Farahani E, Kadner S, Seyboth K, Adler A, Baum I, Brunner S, Eickemeier P, Kriemann B, Savolainen J, Schlömer S, von Stechow C, Zwickel T, Minx JC (eds.)]. Cambridge University Press, Cambridge, United Kingdom and New York. 2014.
Li X, Yang L, Zheng H, Shan Y, Zhang Z, Song M, et al. City-level water-energy nexus in Beijing-Tianjin-Hebei region. Appl Energy. 2019;235:827–34.
Cai B, Lu J, Wang J, Dong H, Liu X, Chen Y, et al. A benchmark city-level carbon dioxide emission inventory for China in 2005. Appl Energy. 2019;233–234:659–73.
Shan Y, Liu J, Liu Z, Xu X, Shao S, Wang P, et al. New provincial CO2 emission inventories in China based on apparent energy consumption data and updated emission factors. Appl Energy. 2016;184:742–50.
Wang H, Lu X, Deng Y, Sun Y, Nielsen CP, Liu Y, et al. China’s CO2 peak before 2030 implied from characteristics and growth of cities. Nat Sustain. 2019;2(8):748–54.
Kort EA, Frankenberg C, Miller CE, Oda T. Space-based observations of megacity carbon dioxide. Geophys Res Lett. 2012;39(17).
Zheng B, Chevallier F, Ciais P, Broquet G, Wang Y, Lian J, et al. Observing carbon dioxide emissions over China’s cities with the Orbiting Carbon Observatory-2. Atmos Chem Phys Discuss. 2020;2020:1–17.
Kennedy CA, Ibrahim N, Hoornweg D. Low-carbon infrastructure strategies for cities. Nat Climate Change. 2014;4(5):343–6.
Wu L, Broquet G, Ciais P, Bellassen V, Vogel F, Chevallier F, et al. What would dense atmospheric observation networks bring to the quantification of city CO2 emissions? Atmos Chem Phys. 2016;16(12):7743–71.
Shusterman AA, Teige VE, Turner AJ, Newman C, Kim J, Cohen RC. The BErkeley atmospheric CO2 observation network: initial evaluation. Atmos Chem Phys. 2016;16(21):13449–63.
Mueller M, Graf P, Meyer J, Pentina A, Dominik B, Perez-Cruz F, et al. Integration and calibration of NDIR CO2 low-cost sensors, and their operation in a sensor network covering Switzerland. Atmos Meas Tech Discuss. 2019;2019:1–25.
Shan Y, Guan D, Hubacek K, Zheng B, Davis SJ, Jia L, et al. City-level climate change mitigation in China. Sci Adv. 2018;4(6):0390.
Shan Y, Guan D, Liu J, Mi Z, Liu Z, Liu J, et al. Methodology and applications of city level CO2 emission accounts in China. J Cleaner Prod. 2017;161:1215–25.
Oda T, Bun R, Kinakh V, Topylko P, Halushchak M, Marland G, et al. Errors and uncertainties in a gridded carbon dioxide emissions inventory. Mitigation Adaptation Strategies Global Change. 2019;24(6):1007–50.
Wang J, Cai B, Zhang L, Cao D, Liu L, Zhou Y, et al. High resolution carbon dioxide emission gridded data for China derived from point sources. Environ Sci Technol. 2014;48(12):7085–93.
Shan Y, Guan D, Zheng H, Ou J, Li Y, Meng J, et al. Data Descriptor: China CO2 emission accounts 1997–2015. Sci Data. 2017;5:170201. https://doi.org/10.1038/sdata.2017.201.
Guan D, Meng J, Reiner DM, Zhang N, Shan Y, Mi Z, et al. Structural decline in China’s CO 2 emissions through transitions in industry and energy systems. Nat Geosci. 2018;11(8):551–5.
Liu M, Wang H, Oda T, Zhao Y, Yang X, Zang R, et al. Refined estimate of China’s CO 2 emissions in spatiotemporal distributions. Atmos Chem Phys. 2013;13(21):10873–82.
Han P, Zeng N, Oda T, Lin X, Crippa M, Guan D, et al. Evaluating China’s fossil-fuel CO2 emissions from a comprehensive dataset of nine inventories. Atmos Chem Phys Discuss. 2020;2020:1–21.
Wang K, Zhang X, Wei Y-M, Yu S. Regional allocation of CO2 emissions allowance over provinces in China by 2020. Energy Policy. 2013;54:214–29.
Zheng B, Zhang Q, Tong D, Chen C, Hong C, Li M, et al. Resolution dependence of uncertainties in gridded emission inventories: a case study in Hebei, China. Atmos Chem Phys. 2017;17(2):921–33.
Liu F, Zhang Q, Tong D, Zheng B, Li M, Huo H, et al. High-resolution inventory of technologies, activities, and emissions of coal-fired power plants in China from 1990 to 2010. Atmos Chem Phys. 2015;15(23):13299–317.
Zheng B, Huo H, Zhang Q, Yao Z, Wang X, Yang X, et al. High-resolution mapping of vehicle emissions in China in 2008. Atmospheric Chemistry & Physics. 2014;14(18).
Shan Y, Liu Z, Guan D. CO2 emissions from China’s lime industry. Appl Energy. 2016;166:245–52.
Guo B, Geng Y, Franke B, Hao H, Liu Y, Chiu A. Uncovering China’s transport CO2 emission patterns at the regional level. Energy Policy. 2014;74:134–46.
Dhakal S. Urban energy use and carbon emissions from cities in China and policy implications. Energy Policy. 2009;37(11):4208–19.
Fang C, Wang S, Li G. Changing urban forms and carbon dioxide emissions in China: a case study of 30 provincial capital cities. Appl Energy. 2015;158:519–31.
Mi Z, Zhang Y, Guan D, Shan Y, Liu Z, Cong R, et al. Consumption-based emission accounting for Chinese cities. Appl Energy. 2016;184:1073–81.
Tan X, Dong L, Chen D, Gu B, Zeng Y. China’s regional CO2 emissions reduction potential: a study of Chongqing city. Appl Energy. 2016;162:1345–54.
Cai B, Cui C, Zhang D, Cao L, Wu P, Pang L, et al. China city-level greenhouse gas emissions inventory in 2015 and uncertainty analysis. Appl Energy. 2019;253:113579.
Zheng B, Tong D, Li M, Liu F, Hong C, Geng G, et al. Trends in China’s anthropogenic emissions since 2010 as the consequence of clean air actions. Atmos Chem Phys. 2018;18:14095–111. https://doi.org/10.5194/acp-18-2018.
Chen W, Wu S, Lei Y, Li S. Interprovincial transfer of embodied energy between the Jing-Jin-Ji area and other provinces in China: a quantification using interprovincial input-output model. Sci Total Environ. 2017;584:990–1003.
Cai B, Li W, Dhakal S, Wang J. Source data supported high resolution carbon emissions inventory for urban areas of the Beijing-Tianjin-Hebei region: spatial patterns, decomposition and policy implications. J Environ Manag. 2018;206:786–99.
NBS NBoSotPsRoC. China Statistical Yearbook 2019. China Statistics Press. 2019.
Meng W, Zhong Q, Chen Y, Shen H, Yun X, Smith KR, et al. Energy and air pollution benefits of household fuel policies in northern China. Proc Natl Acad Sci. 2019;116(34):16773.
NDRC. The 13th five-year plan for energy saving and emissions reduction of the People's Repubic of China (2016–2020). 2016. http://www.ndrc.gov.cn/zcfb/zcfbqt/201701/t20170105_834500.html.
Wang Z, Yang L. Delinking indicators on regional industry development and carbon emissions: Beijing–Tianjin–Hebei economic band case. Ecol Indicators. 2015;48:41–8.
Wang H, Zhang R, Liu M, Bi J. The carbon emissions of Chinese cities. Atmos Chem Phys. 2012;12:6197–206. https://doi.org/10.5194/acp-12-6197-2012.
Mi Z, Zheng J, Meng J, Zheng H, Li X, Coffman DM, et al. Carbon emissions of cities from a consumption-based perspective. Appl Energy. 2019;235:509–18.
Oda T, Maksyutov S, Andres RJ. The Open-source Data Inventory for Anthropogenic CO2, version 2016 (ODIAC2016): a global monthly fossil fuel CO2 gridded emissions data product for tracer transport simulations and surface flux inversions. Earth Syst Sci Data. 2018a;10:87–107. https://doi.org/10.5194/essd-10-87-2018.
Janssens-Maenhout G, Crippa M, Guizzardi D, Muntean M, Schaaf E, Dentener F, Bergamaschi P, Pagliari V, Olivier JGJ, Peters JAHW, van Aardenne JA, Monni S, Doering U, Petrescu AMR. EDGAR v4.3.2 Global Atlas of the three major Greenhouse Gas Emissions for the period 1970–2012. Earth Syst Sci Data Discuss. 2017. https://doi.org/10.5194/essd-2017-79.
Asefi-Najafabady S, Rayner P, Gurney K, McRobert A, Song Y, Coltin K, et al. A multiyear, global gridded fossil fuel CO2 emission data product: evaluation and analysis of results. J Geophys Res Atmospheres. 2014;119(17):10–213.
Cai B, Liang S, Zhou J, Wang J, Cao L, Qu S, et al. China high resolution emission database (CHRED) with point emission sources, gridded emission data, and supplementary socioeconomic data. Resources Conserv Recycling. 2018;129:232–9.
Wang J, Bofeng C, Lixiao Z, Dong C, Lancui L, Ying Z, et al. High resolution carbon dioxide emission gridded data for China derived from point sources. Environ Sci Technol. 2014;48(12):7085–93.
Zheng B, Tong D, Li M, Liu F, Hong C, Geng G, Li H, Li X, Peng L, Qi J, Yan L, Zhang Y, Zhao H, Zheng Y, He K, Zhang Q. Trends in China’s anthropogenic emissions since 2010 as the consequence of clean air actions. Atmos Chem Phys. 2018;18:14095–111. https://doi.org/10.5194/acp-18-2018.
Liu Z, Guan D, Wei W, Davis SJ, Ciais P, Bai J, et al. Reduced carbon emission estimates from fossil fuel combustion and cement production in China. Nature. 2015a;524:335.
Liu M, Wang H, Wang H, Oda T, Zhao Y, Yang X, Zang R, Zang B, Bi J, Chen J. Refined estimate of China’s CO2 emissions in spatiotemporal distributions. Atmos Chem Phys. 2013;13:10873–82. https://doi.org/10.5194/acp-13-2013.
Shan Y, Guan D, Zheng H, Ou J, Li Y, Meng J, et al. China CO2 emission accounts 1997–2015. Sci Data. 2018;5:170201.
Liu Z, Guan D, Wei W, Davis SJ, Ciais P, Bai J, et al. Reduced carbon emission estimates from fossil fuel combustion and cement production in China. Nature. 2015b;524(7565):335–46.
Li M, Zhang Q, Kurokawa J-I, Woo J-H, He K, Lu Z, et al. MIX: a mosaic Asian anthropogenic emission inventory under the international collaboration framework of the MICS-Asia and HTAP. Atmospheric Chem Phys. 2017;17:935–63.
Wang R, Tao S, Ciais P, Shen H, Huang Y, Chen H, et al. High-resolution mapping of combustion processes and implications for CO 2 emissions. Atmos Chem Phys. 2013;13(10):5189–203.
Oda T, Maksyutov S. A very high-resolution (1 km× 1 km) global fossil fuel CO 2 emission inventory derived using a point source database and satellite observations of nighttime lights. Atmos Chem Phys. 2011;11(2):543–56.
Oda T, Maksyutov S, Andres RJ. The Open-source Data Inventory for Anthropogenic CO 2, version 2016 (ODIAC2016): a global monthly fossil fuel CO 2 gridded emissions data product for tracer transport simulations and surface flux inversions. Earth Syst Sci Data. 2018b;10(1):87–107.
Crippa M, Oreggioni G, Guizzardi D, Muntean M, Schaaf E, Lo Vullo E, et al. Fossil CO2 and GHG emissions of all world countries - 2019 Report, EUR 29849 EN, Publications Office of the European Union, Luxembourg, ISBN 978-92-76-11100-9. 2019. https://doi.org/10.2760/687800, JRC117610.
Janssens-Maenhout G, Crippa M, Guizzardi D, Muntean M, Schaaf E, Dentener F, et al. EDGAR v4.3.2 Global Atlas of the three major Greenhouse Gas Emissions for the period 1970–2012. Earth Syst Sci Data. 2019;11(3):959–1002.
Rayner P, Raupach M, Paget M, Peylin P, Koffi E. A new global gridded data set of CO2 emissions from fossil fuel combustion: methodology and evaluation. J Geophys Res Atmospheres. 2010;115:D19306. https://doi.org/10.1029/2009JD013439.
Ramaswami A, Tong K, Fang A, Lal RM, Nagpure AS, Li Y, et al. Urban cross-sector actions for carbon mitigation with local health co-benefits in China. Nat Climate Change. 2017;7(10):736–42.
Gurney KR, Liang J, Patarasuk R, Song Y, Huang J, Roest G. The vulcan version 3.0 high-resolution fossil fuel Co2 emissions for the United States. Earth Syst Sci Data Discuss. 2019;2019:1–42.
Li JS, Zhou HW, Meng J, Yang Q, Chen B, Zhang YY. Carbon emissions and their drivers for a typical urban economy from multiple perspectives: A case analysis for Beijing city. Appl Energy. 2018;226:1076–86.
Shan Y, Guan D, Liu J, Liu Z, Liu J, Schroeder H, et al. CO2 emissions inventory of Chinese cities. Atmos Chem Phys Discuss. 2016;2016:1–26.
NBS NBoSotPsRoC. China Statistical Yearbook 2016. China Statistics Press. 2016.
Guan D, Liu Z, Geng Y, Lindner S, Hubacek K. The gigatonne gap in China’s carbon dioxide inventories. Nat Climate Change. 2012;2(9):672–6.
Hong C, Zhang Q, He K, Guan D, Li M, Liu F, et al. Variations of China’s emission estimates: response to uncertainties in energy statistics. Atmos Chem Phys. 2017;17(2):1227–39.
Chen J, Zhao F, Zeng N, Oda T. Comparing a global high-resolution downscaled fossil fuel CO2 emission dataset to local inventory-based estimates over 14 global cities. Carbon Balance Manag. 2019;15(1):9.
Gurney KR, Liang J, O’Keeffe D, Patarasuk R, Hutchins M, Huang J, et al. Comparison of global downscaled versus bottom-up fossil fuel CO2 emissions at the urban scale in four US urban areas. J Geophys Res Atmospheres. 2019. https://doi.org/10.1029/2018JD028859.
Gately CK, Hutyra LR. Large uncertainties in urban-scale carbon emissions. J Geophys Res Atmospheres. 2017;122(20):11–242.
Wang Y, Ciais P, Broquet G, Bréon FM, Oda T, Lespinas F, et al. A global map of emission clumps for future monitoring of fossil fuel CO2 emissions from space. Earth Syst Sci Data. 2019;11(2):687–703.
Martin CR, Ning Z, Karion A, Dickerson RR, Weber KJ. Evaluation and environmental correction of ambient CO2 measurements from a low-cost NDIR sensor. Atmospheric Measurement Tech. 2017;10(7):2383–95.
Bao Z, Han P, Zeng N, Liu D, Wang Y, Tang G, et al. Observation and modeling of vertical carbon dioxide distribution in a heavily polluted suburban environment. Atmospheric Oceanic Sci Lett. 2020. https://doi.org/10.1080/16742834.2020.1746627.
The People's Government of Beijing Municipality P. Beijing's "13th Five-Year Plan" for energy conservation and consumption reduction and climate change. 2016. http://www.beijing.gov.cn/zfxxgk/110001/szfwj/2016-08/07/content_c7607556c0e74fe58c1c85a5d25183b6.shtml.
This work was supported by the National Key R&D Program of China (No. 2017YFB0504000). Tomohiro Oda is supported by the NASA Carbon Cycle Science program (grant no. NNX14AM76G). We thank Dr. Dabo Gun’s team, Dr. Qiang Zhang, Dr. Shu Tao and Dr. Bofeng Cai for kindly providing CEADs, MEIC, PKU-CO2 and CHRED data.
Authors are listed alphabetically after Qixiang Cai.
This research was supported by the National Key R&D Program of China (No. 2017YFB0504000) and the NASA Carbon Cycle Science program (No. NNX14AM76G).
The authors declare that they have no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Han, P., Zeng, N., Oda, T. et al. A city-level comparison of fossil-fuel and industry processes-induced CO2 emissions over the Beijing-Tianjin-Hebei region from eight emission inventories. Carbon Balance Manage 15, 25 (2020). https://doi.org/10.1186/s13021-020-00163-2