Development of a unit-based industrial emission inventory in the Beijing-Tianjin-Hebei region and resulting improvement in air quality modeling

The Beijing-Tianjin-Hebei (BTH) region is a metropolitan area with the most severe fine particle (PM2.5) pollution in China. Accurate emission inventory plays an important role in air pollution control policy making. In this study, we develop a unit-based emission inventory for industrial sectors in the BTH region, including power plants, industrial boilers, and steel, non-ferrous metal, coking, cement, glass, brick, lime, 15 ceramics, refinery, and chemical industries, based on detailed information for each enterprise, such as location, annual production, production technology/process and air pollution control facilities. In the BTH region, the emissions of sulfur dioxide (SO2), nitrogen oxide (NOx), particulate matter with diameter less than 10 m (PM10), PM2.5, black carbon (BC), organic carbon (OC), and non-methane volatile organic compounds (NMVOCs) from industrial sectors are 869 kt, 1164 kt, 910 kt, 622 kt, 71 kt, 63 kt and 1390 20 kt in 2014, respectively, accounting for 61%, 55%, 62%, 56%, 58%, 22% and 36%, respectively, of the total emissions. Compared with the traditional proxy-based emission inventory, much less emissions in the high-resolution unit-based inventory are allocated to the urban center because of the accurate positioning of industrial enterprises. We apply the Community Multi-scale Air Quality (CMAQ) model simulation to evaluate the unit-based inventory. The simulation results show that the unit-based emission 25 inventory gives better performance of both PM2.5 and gaseous pollutants than the proxy-based emission Atmos. Chem. Phys. Discuss., https://doi.org/10.5194/acp-2018-858 Manuscript under review for journal Atmos. Chem. Phys. Discussion started: 8 October 2018 c © Author(s) 2018. CC BY 4.0 License.


Introduction
The Beijing-Tianjin-Hebei (BTH) region is the political, economic, and cultural center of China.According to China National Environmental Monitoring Centre (2018), in 2017, the annual average concentrations of PM 2.5 in Beijing, Tianjin,and Hebei were 65.6,63.8,and 57.1 µg m −3 , respectively, ranking them second, third and sixth among all provinces.The severe PM 2.5 pollution in the BTH region is largely attributed to the substantial emissions of air pollutants (B.Zhao et al., 2017).An accurate emission inventory, in terms of both emission rates and spatial distribution, is imperative for an adequate understanding of the sources and the formation mechanism of serious air pollution in this area.
The spatial distribution is one of the most uncertain components of emission inventories considering the diverse source categories and complex emission characteristics.The traditional method of spatial allocation is to distribute the emissions by administrative region into grids based on spatial proxies such as population, gross domestic product (GDP), road map, land use data, and nighttime lights (Geng et al., 2017;Oda and Maksyutov, 2011;Streets et al., 2003).The results may deviate significantly from the actual spatial distributions of many sources (Zhou and Gurney, 2011), especially the power and industrial sources, which contribute over 50 % of the total PM 2.5 emissions in China (Zhao et al., 2013a).Due to the stricter air quality regulations and the higher land prices in urban areas, people tend to build factories in suburban areas where the population density and GDP are lower.Zheng et al. (2017) studied the influence of the resolution of gridded emission inventories and found that there were large biases when the inventories were distributed to very fine resolutions following the traditional proxy-based allocation method.The emission inventory could be significantly improved with detailed information regarding point sources such as power plants, steel plants, and cement plants.The high spatial resolution of the inventory may subsequently improve the air quality modeling results and enable a better source apportionment of air pollution (Y.Zhao et al., 2017).
A number of studies have developed the emission inventory in the BTH region (Li et al., 2017;Wang et al., 2014), whereas others have provided emission estimates for this region as part of national or larger-scale emission inventories (Ohara et al., 2007;Stohl et al., 2015).However, only limited studies have estimated the emissions from individual point sources (i.e., a unit-based emission inventory).Zhao et al. (2008), Chen et al. (2014), and Liu et al. (2015) established unit-based emission inventories of coal-fired power plants in China.K. Wang et al. (2016) and Wu et al. (2015) developed an emission inventory for the steel industry.Lei et al. (2011) and Chen et al. (2015) established an emission inventory for the cement industry in China.Qi et al. (2017) established an emission inventory in the BTH region in which power and major industrial sources were treated as point sources.These studies usually focused on one or several major industries, and did not cover all industrial sectors in the BTH region.Moreover, these previous studies seldom validated the unit-based emission inventory or evaluated the improvement it brings to air quality simulation.
In this study, we developed a unit-based emission inventory of industrial sectors for the BTH region.A three-domain nested simulation by the WRF-CMAQ (Weather Research and Forecasting-Community Multi-scale Air Quality) model was applied to evaluate the emission inventory.In order to study the influence of the point sources, we compared the simulation results of this emission inventory with those of a traditional proxy-based emission inventory.

High-resolution emission inventory for the BTH region
A unit-based method is applied to quantify the emissions from industrial sectors such as power plants, industrial boilers, iron and steel production, non-ferrous metal smelters, coking plants, cement, glass, brick, lime, ceramics, refineries, and chemical industries in 2014.The product yields used for estimating emissions of each sector are shown in Table S4 in the Supplement.The pollutant emissions from each industrial enterprise are calculated from activity level (energy consumption for power plants and industrial boilers, and product yield for other sectors), the emission factor, and the removal efficiency of control technology, as shown in the following equation: where E i,j is emissions of pollutant i from industrial enterprise j , A j is the activity level of industrial enterprise j , EF i,j is the uncontrolled emission factor of pollutant i from industrial enterprise j , and η i,j is the removal efficiency of pollutant i by control technology in enterprise j .η i,j is determined by the production process and control technology of the industrial enterprise.The EF i,j values, which depend on the production process of the industrial enterprise, are calculated according to the sulfur and ash contents of fuels, e.g., coal, used in each province (for PM and SO 2 ), or obtained from our previous study (Zhao et al., 2013b), for other pollutants.Some industrial sources involve multiple production processes, such as iron and steel production and cement production.Using cement production as an example, emissions are calculated using the following equation: where E i,j values are the emissions of pollutant i from industrial enterprise j , AK j,m is the amount of clinker produced by the clinker burning process m of the enterprise j , EF i,m is the uncontrolled emission factor for pollutant i from the clinker burning process m, η i,j,m is the removal efficiency of pollutant i from the clinker burning process m in enterprise j , AC j is the amount of cement produced by enterprise j , ef i values are the uncontrolled emission factors from the clinker processing stage (ef i = 0 if i is not particulate matter), η i,j is the removal efficiency of pollutant i in enterprise j .η i,j,m and η i,j both depend on the control technology of the industrial enterprise.
The production processes represented by the first and second terms of Eq. ( 2) are frequently performed in different enterprises.For example, for cement production, clinker may be produced in one enterprise and subsequently processed in another enterprise, which is very common.For each enterprise, we calculate the emission of each production process.Specifically, the total emission of enterprise j is the sum of the emissions of all of the production processes in that enterprise.If processes are divided between multiple enterprises, the emission will be considered in the calculation of the emission of each individual enterprise.
In this study, we collected detailed information for all power and industrial sources except industrial boilers, including latitude/longitude, annual product, production technology/process, and pollution control facilities from a compilation of power industry statistics (China Electricity Council, 2015b) Plume rise is caused by the buoyancy effect and momentum rise (Briggs, 1982).Therefore, stack information, including stack height, flue gas temperature, chimney diameter, and flue gas velocity, is essential for plume rise calculation.For power plants, we get the stack height from the "Compilation of Power Industry Statistics" (China Electricity Council, 2015b).For the stack height of cement factories, we refer to the emission standard of air pollutants for the cement industry (Ministry of Environmental Protection of China, 2013).For the stack height of glass, brick, lime, and ceramic industries, we refer to the emission standard of air pollutants for industrial kilns and furnaces (Ministry of Environmental Protection of China, 1997).For the stack height of non-ferrous metal smelters, coking plants, refineries, and chemical industries, as well as the flue gas temperature, chimney diameter, and flue gas velocity for all industrial sectors, we refer to the national information platform of pollutant discharge permits (http://114.251.10.126/permitExt/outside/default.jsp, last access: 11 March 2019), where very detailed information can be found regarding plants with pollutant discharge permits.For sources without pollutant discharge permits, we use the parameters of plants with a similar production output or coal consumption.Individual information regarding the stacks is applied to each production process.The locations of different processes from the same enterprise are usually assumed to be the same.
The emission inventory for other sources, including residential sources, transportation, solvent use, and open burning, is developed based on the "top-down method" following our previous work (Fu et al., 2013;Wang et al., 2014;Zhao et al., 2013b).The method is the same as Eq. ( 1) except that the emissions are calculated for an individual prefecture-level city rather than individual enterprises.The activity data and technology distribution for each sector are derived based on the statistical yearbooks (Beijing Municipal Bureau of Statistics, 2015;Hebei Municipal Bureau of Statistics, 2015;National Bureau of Statistics, 2015h, g, f, e, i, j, a, b, c, d;Tianjin Municipal Bureau of Statistics, 2015), a wide variety of Chinese technology reports (China Electricity Council, 2015a;National Bureau of Statistics, 2012), and an energy demand modeling approach.Figure S1 shows the energy consumption in the BTH region in 2014.We compared the sum of the energy consumption for each plant with the energy statistics.The sum of individual plants accounts for over 90 % of the energy consumption or product yield reported in the statistics.For the plants not included in the preceding data sources, we calculate the emission using the top-down method.The emission factors are obtained from Zhao et al. (2013b).The speciation of PM 2.5 in both are and point sources is from Fu et al. (2013), whereas the speciation of NMVOCs is updated by Wu et al. (2017).The application rates of removal technologies are obtained from the evolution of emission standards and a variety of technical reports (Chinese State Council, 2013).

Air quality model configuration
In this work, we use CMAQ version 5.0.2 to simulate the concentration of pollutants.A three-domain nested simulation is established as shown in Fig. 1a.The first domain covers almost the entire area of China, Korea, Japan, and parts of India and Southeast Asia with a horizontal grid resolution of 36 km × 36 km.The second domain covers eastern China with a resolution of 12 km × 12 km.The third domain with a horizontal resolution of 4 km × 4 km focuses on the BTH region.The observational sites in the BTH region are marked in Fig. 1b.All of the grids are divided into 14 layers vertically from the surface to an altitude of about 19 km above the ground, and the thickness of the first layer is about 40 m.
In order to minimize the influence of the initial conditions, we choose a 5-day spin-up period.The Carbon Bond 05 (CB05) and AERO6 (Sarwar et al., 2011) are chosen as the gas-phase and aerosol chemical mechanisms, respectively.The simulation periods are January and July of 2014, representing winter and summer, respectively.
We use the Weather Research and Forecasting (WRF) model version 3.7.1 (Skamarock et al., 2008) to simulate the meteorological fields.The physics options for the WRF simulation are the Kain-Fritsch cumulus scheme (Kain, 2004), the Morrison double-moment scheme for cloud microphysics (Morrison et al., 2005), the Pleim-Xiu land surface model (Xiu and Pleim, 2001), the Pleim-Xiu surface layer scheme (Pleim, 2006), the Asymmetric Convective Model (ACM2; Pleim) boundary layer parameterization (Pleim, 2007), and the Rapid Radiative Transfer Model for GCMs radiation scheme (Mlawer et al., 1997).The meteorological initial and boundary conditions are generated from the Final Operational Global Analysis data (ds083.2) of the National Center for Environmental Prediction (NCEP) at 1.0 • × 1.0 • and 6 h resolutions.Default profile data are used for chemical initial and boundary conditions.The Meteorology Chemistry Interface Processor (MCIP) version 4.1 is applied to process the meteorological data into the format required by CMAQ.The simulated wind speed, wind direction, temperature, and humidity agree well with the observation data from the National Climate Data Center (NCDC), as detailed in the Supplement.
In order to evaluate the high-resolution emission inventory with the unit-based industrial sources, we develop a traditional proxy-based emission inventory with the same amount of emissions and compare the simulation results of these two emission inventories.In the proxy-based emission inventory, all sectors are allocated as area sources using spatial proxies such as population, GDP, road map, and land use data.The proxies used for each sector are described in detail in Table S2.In order to separate the influences of the horizontal and vertical distributions of the emissions, we developed another unit-based inventory with emission heights the same as the proxy-based inventory; we call this inventory the "hypounit-based inventory".The anthropogenic emission inventories for other provinces in China were developed in our previous studies (Wang et al., 2014;Zhao et al., 2018).The emissions outside China are obtained from the MIX emission inventory (Li et al., 2017) for 2010, which is the most current year available.In the simulation with the unit-based inventory, plume rise is calculated using the built-in algorithm in CMAQ.Meteorological data are used to calculate plume rise for all point sources.Then, the plume is distributed into the vertical layers that the plume intersects based on the pressure in each layer.

Air pollutant emissions in the BTH region
In the BTH region, the emissions of sulfur dioxide (SO 2 ), nitrogen oxide (NO x ), PM 10 , PM 2.5 , black carbon (BC), organic carbon (OC), non-methane volatile organic compounds (NMVOCs) and ammonia (NH 3 ) were 1417, 2100, 1479, 1106, 213, 289, 2381, and 712 kt in 2014, respectively. Figure , respectively. Figure 2 shows the sectoral emissions for major pollutants in the BTH region by city.Figure S2 shows the NMVOCs speciation by sector.The emission estimates are compared with previous studies in Fig. S3. Figure 3 shows the locations and emissions of power and industrial sources.
Power plants account for 13 %, 16 %, and 4 % of the total SO 2 , NO x , and PM 2.5 emissions, respectively, whereas the contributions to NMVOC and NH 3 emissions are neg- The emissions from industrial boilers account for 27 %, 19 %, 8 %, 1 %, and < 1 % of the total SO 2 , NO x , PM 2.5 , NMVOCs, and NH 3 emissions, respectively.As shown in Fig. 3, there are many industrial boilers in the BTH region.Industrial boilers are one of the most important emission sources for SO 2 and NO x .
The emissions from cement contribute 6 %, 9 %, and 10 % of the total SO 2 , NO x , and PM 2.5 emissions, respectively, whereas the contributions to NMVOC and NH 3 emissions are negligible (< 1 %).Most cement plants are located in southern and eastern Hebei.
The emissions from steel production represent 8 %, 3 %, and 22 % of the total SO 2 , NO x , and PM 2.5 emissions, respectively, whereas the contributions to NMVOC and NH 3 emissions are negligible (< 1 %).Tangshan has the largest number of steel plants in the BTH region, and steel production accounts for over half of the PM 2.5 emissions in Tangshan.
Besides the aforementioned sectors, 8 %, 8 %, 13 %, 36 %, and < 1 % of the total respective SO 2 , NO x , PM 2.5 , NMVOCs, and NH 3 emissions come from other industrial processes (chemistry, coking plants, nonferrous metal smelting, brick, ceramics, lime, glass, and refineries).Industrial processes are the most important emission source for NMVOCs, accounting for nearly half of the emissions in Tianjin and Shijiazhuang.
Considering the large contribution of industrial sources to the total emissions, the application of a unit-based method results in remarkable changes in the spatial distribution of air pollutant emissions.The emission rates of PM 2.5 , NO x , and SO 2 in the proxy-based and the unit-based inventories and their differences are shown in Fig. 4. In the unit-based emission inventory, the emissions are lower than those in the proxy-based emission inventory in urban centers in the BTH region.Instead, a large amount of the emissions are concentrated in certain points in suburban areas, where large plants are located.

Evaluation of the unit-based emission inventory
In order to study the accuracy of the unit-based inventory, the simulation results of SO 2 , NO 2 , O 3 , and PM 2.5 with the unitbased inventory are compared with the observational data from the China National Environmental Monitoring Centre.The observations are available for 80 sites located in 13 cities in the BTH region, including 70 sites in urban areas and 10 sites in suburban areas.The accurate locations of urban and suburban sites in Beijing are shown in Figs.S5-S6.The analysis of the results is shown in Table 1.We use the normalized mean bias (NMB), the normalized mean error (NME), the mean fractional bias (MFB), and the mean fractional error (MFE) (U.S. EPA, 2007) to quantitatively evaluate the model performance.
SO 2 and NO 2 are precursors of PM 2.5 , so we first compare the simulation results of gaseous pollutants with observations.For NO 2 , the results with the proxy-based inventory overestimate the observations by 22 %, whereas results with the unit-based inventory overestimate the observations by 9 % in January.Similarly, in July, the simulated NO 2 concentrations show an overestimation in simulations with both inventories, but the overestimation is less with the unit-based inventory.The simulation results of SO 2 are similar to those of NO 2 .However, the overestimation is higher with both inventories, and the differences between the concentrations with the two inventories are larger.The overestimation of SO 2 concentrations may be due to the lack of several SO 2 reaction mechanisms in CMAQ, such as heterogeneous reactions of SO 2 on the surface of dust particles (Fu et al., 2016), the oxidation of SO 2 by NO x in aerosol liquid water (Cheng et al., 2016;G. Wang et al., 2016), and the effects of SO 2 and NH 3 on secondary organic aerosol formation (Chu et al., 2016).It may also be due to uncertainty in the emission inventory, especially the uncertainty regarding the removal efficiencies of SO 2 control facilities.The biased spatial distribution of SO 2 emissions from residential combustion may also contribute to the overestimation.A large fraction of residential combustion takes place in rural areas.However, in this work the emissions from residential combustion are allocated by GDP and population, which leads to an overestimation of SO 2 emissions in urban areas and hence an overestimation of the SO 2 concentration.
For O 3 , the simulation results in January with the proxybased inventory underestimate the observations by 21 %, whereas the results with the unit-based inventory underestimate the observations by only 5 %.The simulation results in July follow the same trend.China is experiencing increasingly severe O 3 pollution (Li et al., 2019), which usually occurs in summer.Therefore, we analyze two extra indices of O 3 , 1 h-peak O 3 and daily maximum 8 h averaged (MDA8) O 3 concentrations in July, which are shown in Table 2.The results of 1 h-peak O 3 and MDA8 O 3 concentrations are similar to those of the monthly average O 3 concentration.The concentration with the unit-based inventory is slightly higher than that with the proxy-based inventory and is also closer to the observation.The reason for the changes in the O 3 concentrations will be discussed later.
The simulated PM 2.5 concentrations with the unit-based inventory are lower than those with the proxy-based inventory in both winter and summer.In January, the simulated PM 2.5 concentrations with the proxy-based inventory overestimate the observed values by 25 %, whereas the overestimation is only 7 % with the unit-based inventory.In July, the simulated PM 2.5 concentrations with both inventories are 17 % and 30 % lower than the observations, respectively.An overall underestimation is expected because the default CMAQ model significantly underestimates the concentrations of secondary organic aerosols (Zhao et al., 2016) and the fugitive dust emission is not included in the emission inventory.According to Boylan and Russell (2006), the simulation results of PM are acceptable when the mean fractional bias (MFB) is less than or equal to ±60 % and the mean fractional error (MFE) is less than 75 %; furthermore, a model performance goal is met when the MFB is less than ±30 % and the MFE is less than 50 %.The statistical indices of the simulation results for PM 2.5 with both inventories and both months are within the performance goal values, which means that the simulation results are relatively accurate.
Figure 5 further shows the spatial distribution of SO 2 , NO 2 , O 3 , 1 h-peak O 3 , MDA8 O 3 , and PM 2.5 concentrations with the proxy-based inventory, and the differences between the other two simulations and the proxy-based inventory.For SO 2 , NO 2 , and PM 2.5 , the concentrations in urban areas are generally higher with the proxy-based inventory than those with the unit-based inventory, especially in winter.In January, large differences with respect to simulated concentrations with the two inventories are found in urban Tianjin, Tangshan, Baoding, and Shijiazhuang, where a large amount of industrial emissions are allocated in the proxy-based inventory due to the large population densities in these areas.The simulations for July follow the same pattern but the concentrations and the differences between the concentrations with the two inventories are lower than those for January.In some areas where many factories are located, such as the northern part of Xingtai city, the concentration   with the unit-based inventory is higher because of the high emission intensity.There are two reasons for the difference between the results with the two inventories.The first reason is the spatial distribution.With detailed information regarding the industrial sectors, more emissions are allocated to certain locations in suburban/rural areas in the unit-based emission inventory.From "Diff1" (hypo-unit-based minus proxy-based), we can see that the improved horizontal distribution of the unit-based emission inventory significantly decreases the PM 2.5 , SO 2 , and NO 2 concentrations in most urban centers, and significantly increases the concentrations in a large fraction of suburban and rural areas, especially the areas where large industrial plants are located.The other reason is the vertical distribution.Plume rise is calculated in the simulation with the unit-based inventory, which causes the difference in emissions in the vertical layers.Therefore, the higher the pollutants are emitted, the lower the ground concentration becomes.From the differences between Diff1 and Diff2 we can see that plume rise leads to lower concentrations over the whole region.
For O 3 , the difference in the concentrations is evident but is opposite to that for PM 2.5 .This is because the urban cen- ters of Beijing/Tianjin are located within the VOC-control chemical regime (Liu et al., 2010).Hence, the emissions of NO x in the surface layer are lower in the unit-based inventory than in the proxy-based inventory, which leads to higher O 3 concentrations in urban areas.
The spatial distribution of the concentrations of these pollutants are significantly heterogeneous.The NME and MFE of most pollutants averaged over a 2-month period are lower with the unit-based inventory than with the proxy-based inventory; this means that the spatial distribution with the unitbased inventory more with the observations than that of the proxy-based inventory.For SO 2 , NO 2 , and PM 2.5 , peak concentrations usually occur in the urban center whereas the opposite is noted for O 3 .We apply the "concentration gradient" metric, which is defined as the ratio of the urban monthly mean concentrations to the suburban concentrations, to quan-  titatively characterize the heterogeneous spatial distributions.We calculate the concentration gradients for Beijing and Tianjin (Fig. 6), as both urban and suburban observational sites are found in these two cities.The concentration gradients of NO 2 and SO 2 between the urban and suburban areas are closer to the observations in the simulation with the unitbased inventory than those with the proxy-based inventory (Fig. 6).The simulated O 3 concentration gradients with the unit-based inventory, the proxy-based inventory, and the observations are 0.7, 0.5, and 0.9 in January and 0.9, 0.8, and 1.1 in July, respectively.As for the 1h-peak O 3 and MDA8 O 3 in July, the simulated results with the unit-based inventory are also closer to the observations.As previously stated, this is explained by the VOC-limited photochemical regime and lower NO x emissions in the unit-based inventory over urban areas.As for PM 2.5 , the concentration gradients for simulations with the unit-based inventory, the proxy-based inventory, and the observations in Beijing are 1.6, 2.1, and 1.5 in January and 1.3, 1.5, and 1.3 in July, respectively.The results imply that the unit-based emission inventory better reproduces the distributions of pollutant emissions between the urban and suburban areas.
To further elucidate the reasons for the difference between the PM 2.5 concentrations with two emission inventories, we examine the simulation results of different chemical components, including sulfate (SO 2− 4 ), nitrate (NO − 3 ), ammonium (NH + 4 ), elemental carbon (EC), and organic carbon (OC), as shown in Fig. 7 and Table 2.The concentrations of EC and OC in the simulation with the unit-based inventory are generally lower than those with the proxy-based inventory in both January and July, especially in urban Beijing, Baoding, and Shijiazhuang.This pattern is similar to that of PM 2.5 .In some cities such as Xingtai, the concentrations of EC and OC in the simulation with the unit-based inventory are slightly higher than those with the proxy-based inventory.
The results for secondary inorganic aerosols are quite different.From Fig. 7 and Table 3 we can see that the SO 2− 4 concentrations are lower in most areas in the simulation with the unit-based inventory compared with concentrations with the proxy-based inventory, which is due to the fact that the sensitivity of SO 2− 4 concentrations to SO 2 concentrations is positive during all months (B.Zhao et al., 2017).The differences in the concentrations of SO 2− 4 are similar to those of SO 2 , which is shown in Fig. 5.The difference in the NH + 4 concentration is relatively small compared with other components.As for NO − 3 , the concentration of NO − 3 in the simulation with the unit-based inventory is much higher than that with the proxy-based inventory in winter, whereas the differences between the results with two inventories vary with location in summer.SO 2− 4 concentrations in the unit-based approach are much lower than in the proxy-based approach, whereas NH + 4 is almost constant, as shown in Fig. 7.In this case, more HNO 3 is converted to NO − 3 with excess NH + 4 , whereas these processes depend on abundance of HNO 3 or NH 3 .Taking all of the chemical components into account, the primary components account for most of the differences in the PM 2.5 concentrations between the simulations with the two inventories.However, the complex responses of various secondary components often counteract each other (especially in January), leading to an overall smaller contribution of secondary components to the PM 2.5 concentration differences.

Conclusion
In this study, we developed a high-resolution emission inventory of major pollutants for the BTH region for the year 2014 Atmos.Chem.Phys., 19,[3447][3448][3449][3450][3451][3452][3453][3454][3455][3456][3457][3458][3459][3460][3461][3462]2019 www.atmos-chem-phys.net/19/3447/2019/The emissions in the unit-based emission inventory are lower than those in the proxy-based emission inventory in most urban centers in the BTH region due to the concentrated emissions from point sources.The application of the unit-based emission inventory improves model-observation agreement for most pollutants.The accurate location of point sources leads to lower concentrations of primary pollutants in urban areas and higher concentrations in suburban areas.Plume rise accounts for the lower concentrations over the whole region.For SO 2 , NO 2 , and PM 2.5 , the concentrations in urban areas decrease significantly and become closer to the observations, mostly due to the decrease in urban emissions.For O 3 , the concentrations in urban areas increase slightly and also show better agreement with observations, mainly due to the more reasonable allocation of NO x emissions.The improvement is particularly significant for the urbansuburban concentration gradients.For PM 2.5 , the concentration gradients for the simulations with the unit-based inventory, the proxy-based inventory, and the observations in Beijing are 1.6, 2.1, and 1.5 in January and 1.3, 1.5, and 1.3 in July, respectively.For O 3 , the corresponding values are 0.7, 0.5, and 0.9 in January and 0.9, 0.8, and 1.1 in July, implying that the unit-based emission inventory better reproduces the distributions of pollutant emissions between the urban and suburban areas.
The unit-based industrial emission inventory enables more accurate source apportionment and more reliable research on the mechanism of air pollution formation; therefore, it contributes to the development of more precisely targeted control policies.To further improve the emission inventory, it is necessary to improve the spatial allocation of emissions from non-industrial sectors, such as the residential and commercial sectors.Our previous study provides an example of the development of a village-based residential emission inventory in rural Beijing (Cai et al., 2018).Such studies on high-resolution emission inventories, for both industrial and non-industrial sources, are highly needed and should also be extended to other provinces and/or regions.In addition, the plume-in-grid approach might help to further improve model performance, which merits further in-depth study.
Data availability.All data needed to evaluate the conclusion of this paper are provided in the main text and the Supplement.Additional related data are available upon request.
Author contributions.SW and BZ designed the research.HZ, SC, and XC performed the research.HZ, SC, BZ, SW, and XC analyzed the results.HZ, SC, BZ, SW, XC, and JH wrote the paper.HZ and SC contributed equally to this study.
Competing interests.The authors declare that they have no conflict of interest.Special issue statement.This article is part of the special issue "Regional assessment of air pollution and climate change over East and Southeast Asia: results from MICS-Asia Phase III".It is not associated with a conference.

Figure 1 .
Figure 1.The three-domain nested CMAQ domains (a) and the observational sites in the BTH region (b).

Figure 2 .
Figure 2. Sectoral contributions to emissions in the BTH region in 2014.

Figure 3 .
Figure 3. Locations and emissions of industrial sources in the BTH region.The industrial plants are divided into four groups for the sake of clarity.

Figure 4 .
Figure 4. Emission rate of PM 2.5 , NO x , and SO 2 emissions of the proxy-based (a) and the unit-based (b) inventories and their differences (unit-based minus proxy-based)( c).Note that the emissions are the same in provinces other than Beijing, Tianjin, and Hebei.

Figure 5 .
Figure 5. Spatial distribution of the monthly (January and July) mean concentrations of SO 2 , NO 2 , O 3 , 1 h-peak O 3 , MDA8 O 3 , and PM 2.5 with the proxy-based inventory, and the differences between the other two simulations and the proxy-based inventory (Diff1 refers to the hypo-unit-based minus proxy-based, and Diff2 denotes the unit-based minus the proxy-based).The units are µg m −3 for all panels.

Figure 6 .
Figure 6.Observed and simulated concentration gradients of SO 2 , NO 2 , O 3 , and PM 2.5 with the proxy-based and the unit-based inventories in Beijing (a) and Tianjin (b).

Figure 7 .
Figure 7.The differences (unit: µg m −3 ) in the simulation results of the components of PM 2.5 between the results from the two inventories (unit-based minus proxy-based).

Table 1 .
The statistics for model performance of PM 2.5 , NO 2 , SO 2 , 1 O 3 and daily maximum 8 h averaged (MDA8) O 3 in January and July of 2014 with the proxy-based and the unit-based inventories.The following abbreviations are used in the table: simulated (SIM), observed (OBS), normalized mean bias (NMB), normalized mean error (NME), mean fractional bias (MFB), and mean fractional error (MFE).

Table 2 .
The statistics for model performance of 1 h-peak NO 2 and daily maximum 8 h averaged (MDA8) O 3 concentrations in July of 2014 with the proxy-based and the unit-based inventories.See the caption of Table1for an explanation of the abbreviations used.

Table 3 .
The mean concentrations (unit: µg m −3 ) of the components of PM 2.5 with the proxy-based and the unit-based inventories and their differences.The following species are listed in this table: sulfate (SO 2− 4 ), nitrate (NO − 3 ), ammonium (NH + 4 ), elemental carbon (EC), and organic carbon (OC).