A numerical study of the contribution to the air pollutant in Beijing during CAREBeijing-2006

Q. Z. Wu, Z. F. Wang, A. Gbaguidi, C. Gao, L. N. Li, and W. Wang State Key Laboratory of Atmospheric Boundary Layer Physics and Atmospheric Chemistry (LAPC), Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing, China College of Global Change and Earth System Science, Beijing Normal University, Beijing, China Beijing Municipal Environmental Protection Bureau, Beijing, China Graduate University of Chinese Academy of Sciences, Beijing, China


Introduction
In order to improve the air quality in Beijing, the Campaign of Air Quality Research in Beijing (CAREBeijing), which was initiated by the Beijing Municipal Environmental Protection Bureau (Beijing EPB) and Peking University (PKU) in conjunction with other eighteen institutes and research centers, was conducted in August from 2006 to 2008.The present numerical study is to provide an in-deep quantitative analysis on the contributions from local and regional sources to the air pollution in Beijing.Introduction

Conclusions References
Tables Figures

Back Close
Full Air pollution in Beijing is regional in nature and not attributable only to local emissions (Streets et al., 2007), that is due to its semi-basin geophysical feature and unique weather pattern (Ren et al., 2003).Making insights into the main pathways of the surrounding air pollutants contribution to Beijing is therefore vital for the establishment of effective control strategies on major pollutants reduction in the mega-city.Xu et al. (2005) investigated the transport pathways and found a convergence zone of air pollutants along the Yanshan-Taihang Mountains.However, accurate identification of the major pollutants sources and transport processes to thoroughly understand the mechanisms behind air pollution over Beijing remains a challenge task so far.
Recently, An et al. (2007) applied CMAQ model to study a heavy particulate and sulfur episode observed in spring in Beijing by using emission switch on/off method, in his results, the contribution to the fine particles in Beijing was about 39%∼53% from the northwest (notably from Zhangjiakou), and 15% from the southwest (notably from Baoding).Similarly, Streets et al. (2007) reported that PM 2.5 concentrations from Shandong Province might reach a maximum contribution of 63% to Beijing under certain meteorological conditions, Tianjin also exerted a significant influence on Beijing PM 2.5 concentration with the maximum contribution of 42% and the average contribution of 13%.Comparatively, the influences from regional sources on ozone were less pronounced, with maximum contributions of 28% from Hebei, 27% from Tianjin, and 24% from Shandong province, but not at the same time.According to Chen et al. (2007) of the MM5-CMAQ modeling with five emission scenarios in the months of January, April, July and October in 2002, much efforts should be placed on demanding more pollution reduction and better emission control from surrounding provinces, while Beijing needed to take positive steps to reduce itself own emissions.These methods commonly used to estimate the contribution of targeted emissions, highlight significant bias in non-linear source-receptor relationships analysis.And the potential contributions at the different layers were not described.This study, through simulations carried on the NAQPMS model, adopts air pollutant tagged module to estimate the contributions of local and regional emissions to the Introduction

Conclusions References
Tables Figures

Back Close
Full air pollutants in Beijing at the different layers.More detailed scales (by dissociating the whole Beijing and counties) analysis on air pollutants transport from regional emissions to Beijing, is also presented.That remains an extremely important issue for air quality management and emissions control efficiency.Section 2 describes the model and evaluation, Sect. 3 presents results analysis and Sect. 4 brings out conclusions and environmental implications.

Model description
The Nested Air Quality Prediction Model System (NAQPMS) is a fully modularized, three-dimensional chemical transport model over regional and urban-scale air pollution (Wang et al., 2002(Wang et al., , 2006)).The chemical transport module reproduces the physical and chemical evolution of reactive pollutants by solving the mass balance equation in the terrain-following coordinates, and includes advection and diffusion processes, gas/aqueous chemistry and dry/wet deposition processes.More details are given elsewhere (Wang et al., 2002(Wang et al., , 2006;;Wu;Li et al., 2011), and its gas chemistry module is updated to CBM-Z (Zaveri and Peters, 1999) in (Li et al., 2007).The air pollutant tagged module was developed and coupled into the NAQPMS model (Li et al., 2008) and described in Sect.2.2.The model has been applied to study the issues related to the transport and chemical processes of air pollutants and the interaction between mineral aerosols and acid rain over east Asia in the regional scale (Huang et al., 1995;Wang et al., 2002;Li et al., 2007), and also the impact of mobile sources on air quality in Beijing in the urban scale (Wu et al., 2010b).
The NCAR/PSU Fifth-Generation Mesoscale Model (MM5) (Grell et al., 1995) Version 3.6 is applied as the meteorological driver for the chemical transport module.The National Centers for Environmental Prediction (NCEP) global final analysis data (FNL), with 1 • × 1 • and four times a day, are used for the meteorological initial and boundary Introduction

Conclusions References
Tables Figures

Back Close
Full conditions.The select schemes in this study were the simple ice for explicit moisture, Grell cumulus, MRF for PBL and cloud scheme for atmospheric radiation according to the sensitivity numerical experiments in (Gao et al., 2007).Three nested domains for horizontal resolution of 81 km, 27 km and 9 km with the center location at (35.0 • N, 110 • E) were used for both MM5 and NAQPMS model (Fig. 1).Vertically, the NAQPMS model had 20 p-vertical layers extending from 1000 hPa to 100 hPa, and its meteorological field was interpolated from 23 sigma layers of the MM5 model using INTERPB module.
The simulation was performed to the period of 20 July∼31 August 2006, and the first twelve days were used as spin-up time to reduce the influence of the initial concentration to the model.The influence of the initial conditions would be decreasing while the simulation goes further (Seinfeld and Pandis, 2006).In addition, since the simulations of the MM5 model were initialized with a cold start, and might take at least 12 hours on average for the model to spin up (Colle et al., 1999;Lo et al., 2007), we finished the meteorological simulation with predicting cycle, consisting of making 36-h simulation and taking the last 24-h simulation as the meteorological driver of the chemical transport module in every predicting cycle.

Air pollutant tagged method
As mentioned above, an online air pollutant tagged module was developed in the NAQPMS model and applied in this study.In this module, we assume that different tagged species (S T ) are properly integrated in each grid, each tagged species S T shares the same loss coefficients with total species during outflow, chemical transformations and dry/wet removal processes, and the same as the total species (S) in the removal processes (Kondo et al., 2004;Davis et al., 2003).So all of the removal processes have no alteration on the fraction (F S T ) of the tagged species in the total species.Introduction

Conclusions References
Tables Figures

Back Close
Full Thus, the tendency of F S T depends on the produce processes and can be described as following: where j represents the grid.P T,j is the production of the tagged species in the j th grid, including the gross chemical production and emissions.When the j th grid outside the tagged emission region, the emission is equal to zero.(M T,j ) dif+adv+conv is the S T inflow from adjacent grid due to advection, diffusion and convection processes, and calculated by the flow fluxes and the S T fraction in those adjacent grids.Compared with the sensitive analysis with switch on/off emission, the online tagged method provides a different and more efficient measurement of relative importance of various tagged regions.Similar analysis can also be found in previous works (Davis et al., 2003;Kondo et al., 2004;Sudo and Akimoto, 2007;Li et al., 2008).The Supplement provides more details on the tagged method, for example how the PM 10 from Baoding and ozone from Tianjin has been tagged.This study tagged 22 emissions regions in the model domain, according to the regional administrative system and the position relative to Beijing (shown in Table  The lateral, top boundary and initial pollutants are also tagged.The GIS software (Arc-GIS 8.3) was applied to properly provide readable regional labeled database for the model.

Emission inventory
The emissions from the Beijing and its surrounding areas are especially large (Matsui et al., 2009) and very important for modeling studies.Several anthropogenic emission inventories for this study were developed for the year 2006 with a horizontal resolution of 0.5 × 0.5 • by revising those reported by Streets et al. (2003Streets et al. ( , 2007) ) using the method described by Zhang et al. (2009).Seasonal dependencies of emissions are not taken into account.Four emission inventories are considered in this study: the regional emissions with 10km resolution update from TRACE P emission inventory according to personal communications (Streets et al., 2003); the power plant emission database of Beijing and its surrounding provinces including Tianjin, Hebei, Shanxi, Inner Mongolia and Shandong provinces from the Chinese Research Academy of Environmental Sciences (CRAES); the main industrial emissions in Tianjin and Hebei provinces from Beijing University of Technology (BUT) and CRAES; the detailed local emissions database of Beijing city from Beijing EPB.The amount of emissions used is presented in Table 2.For all species listed in Table 2, the emissions used are greater than those reported by Streets et al. (2003) by a factor of 1.6 to 2.7 in the inner model domain (D3).These increases are mostly due to actual emissions growth since 2002, in addition to improvements in estimation methodology provided by Zhang et al. (2009).

Model evaluation
In order to estimate the model performances, the statistical parameters were introduced and defined in appendix.Sim(i ) and Obs(i ) represent the simulated and the corresponding observed variables respectively.Introduction

Conclusions References
Tables Figures

Back Close
Full The observed air quality data were provided from the campaign observation sites and the National Standard Air Quality Observation Stations (NSAQ stations) in urban and outer Beijing during CAREBeijing-2006(Garland et al., 2009).Details of the observation sites and measurements can be found in JGR special section paper about CAREBeijing-2006, such as (Garland et al., 2009;Takegawa et al., 2009).The model evaluation about SO 2 , NO 2 , PM 10 and O 3 with the observed hourly data at IAP tower stations (at 120 m and 280 m height), and SO 2 , NO 2 and O 3 at Yufa and Xinglong stations are presented in this study.The observation sites are shown in Fig. 2, of which the IAP tower station is located in the urban Beijing, and the Yufa station to the south of Beijing.Xinglong stations are located at the east of Beijing, where the air pollutants transport from the south and the east respectively were monitored.
Statistical parameters were used to evaluate the model results.Detailed comparisons between the hourly observed and simulated SO 2 , PM 10 , NO 2 and Ozone at the IAP tower station with the scatter plots are presented in (Fig. 3), illustrating properly the model performance in the urban Beijing.The scatter plot of the SO 2 concentration indicated that the NAQPMS model results overall agreed with the observation in urban Beijing with the correlation coefficients of 0.78 at Layer-2 and 0.73 at Layer-3, respectively.And the observed mean concentrations were 17.97 ppbv and 12.21 ppbv at the Layer-2 and Layer-3 respectively, which are close to the simulated values of 22.45 ppbv and 15.17 ppbv respectively.The vertical characteristics that the observed SO 2 concentration at the lower layer (Layer-2) was higher than the higher layer (Layer-3) was also simulated out.The model showed a high model performance on SO 2 .At the Layer-2, the simulated NO 2 concentration (mean: 35.16 ppbv) was close to the observed (mean: 35.53 ppbv) according to Fig. 3, and the simulated PM 10 and ozone concentration (shown in Fig. 3) were also highly correlated to the observations (r = 0.86, 0.88 and 0.85 respectively).However, the simulated PM 10 averaged concentration was 97.61 µg m −3 , which was much higher than the observed of 68.49 µg m −3 , implying a higher RMSE of about 35.89 µg m −3 .This controversy might be due to some environmental measures implemented in Beijing for reduction of PM 10 emissions in the urban Introduction

Conclusions References
Tables Figures

Back Close
Full Beijing (Hao et al., 2008) were not considered in the emissions inventory.Compared with the PM 10 , the average concentration of simulated ozone seems to be closer to the observation at the Layer-3.The mean bias and normal mean bias of ozone concentration at the Layer-2 are 6.42 ppbv and 0.27 respectively, in comparison to 1.72 ppbv and 0.06 at the Layer-3, suggested that the model had a better model performance at the Layer-2.
As shown in the Table 3 and Table 4, the simulation in this study appears to be more accurate than many previous numerical studies in Beijing, such as (An et al., 2007) and (Matsui et al., 2009).The updated emission data significantly reduced the model biases and errors, and the model results in urban Beijing were reasonable for the evaluation of the air pollutant contribution from different sources.
However, the model results showed worse performance at other sites (Yufa and Xinglong).Despite the high correlation coefficients of NO 2 (0.58) and O 3 (0.62) at Yufa station, most statistical parameters described in Table 5 for the model performance at Yufa and Xinglong stations were weaker than those at the IAP tower station, e.g. the normal mean bias at Xinglong site was greater than that at the IAP tower station.This suggested a modeling discrepancy between the urban and rural sites (Yufa and Xinglong were two rural sites).
Even so, from the time series plots of the observed/simulated SO 2 , NO 2 and ozone concentration (Fig. 4), we could find the NAQPMS model actually simulated out the trends of the pollutants, especially ozone.The model results not only showed good agreements in trends, but also the level of ozone concentration at both Yufa and Xinglong station.The ozone trend depends more on the daily radiation, but the level of concentration is limited by the amount of its precursor concentrations.Meanwhile, the model performance of ozone concentration was better than that of SO 2 and NO 2 .In linkage with trends in ozone concentrations, such characteristics indicated the complex photochemical relationships, particularly between NO 2 and ozone, which is supported by many previous studies about photochemical pollution (Streets et al., 2007;Su, 2005).Moreover, the simulated SO 2 at Xinglong station was much lower than the observed in the whole concentration distribution, indicating that the SO 2 emissions to the east of Beijing might be underestimated.

Results and discussion
To better understand the impact of regional transport to the air pollutants in Beijing, we identified the potential regional contribution at the different vertical layers as well as the principal direct entry-ways of non-Beijing pollutants in this study.The monthly mean contribution percentage was estimated at each layer, and the comparisons between different tendencies were presented.In addition, the counties scale analysis was performed to better characterize the regional contribution to the air pollution in Beijing.

The contribution of the SO 2 concentration to Beijing
In August 2006, four obvious pollution accumulation/extinction processes (3-10, 14-19, 22-26 and 27-29) can be separated through the daily variations of SO 2 contribution (Fig. 5).On 7, 15, 16 and 25 August, there were clear evidences showing significant SO 2 contribution from Tianjin (TJ) with the daily range of 12%∼25%.For example, the contribution from TJ under the southeasterly winds in the relative high pollution event on August 25 was about 25% (2/3 of the total surrounding areas contribution).The increase of SO 2 concentration on 22-23 August was associated with the contribution (30%∼40%) from Hebei Province (essentially from BD and LF); the non-Beijing contribution in the last sequence (on 27-29 August) seemed to be also dominated by the contribution from TJ, but with less amount (about 15%).Such important direct contribution from closer neighboring areas illustrated the important roles surrounding regions played in affecting the air pollution in Beijing.Particularly, the emission controls on industrial and vehicles sources in these surrounding regions are not as stringent as those in Beijing, and the emissions of pollutants are often high.Introduction

Conclusions References
Tables Figures

Back Close
Full Analysis of the monthly mean SO 2 concentrations and its contribution over Beijing is presented in Figs. 6 and 7. Results suggested a significant regional contribution of SO 2 from surrounding areas to Beijing at the upper layer (1.1 km), mostly from Shanxi, Shandong and Inner Mongolia.This indicated that there might have large SO 2 emissions in these regions.The local contribution to Beijing was more significant at the surface layer.The monthly mean concentration from local sources was about 5.53 ppbv at the surface layer, accounting for more than 65% of the total.In comparison, the contribution from local sources only accounted for 14% of the total at the upper layer.However, the surrounding areas which are closer to Beijing (especially TJ, BD, LF, ZJK and CD) also showed considerable contributions to the surface layer, with the total monthly contribution about 28%.
At the upper layer, the most significant contribution to the monthly mean SO 2 concentration in Beijing was from the southern regions (about 50% of total contribution), of which TJ accounted for 19%, SD for 11%, and XH for 10% respectively.The total contribution from surrounding areas at this layer was estimated to be about 86%.Consistent with the previous studies (An et al., 2007;Streets et al., 2007), the results in this study further confirm the long range transport of sulfur dioxide to Beijing, and the significant SO 2 emissions in those regions relatively far from Beijing.Our results also provided an estimation on the contributions from the surrounding areas to the different vertical layers, with the contribution of 86% at the upper layer and 30% at the surface layer.
These findings suggested the important roles that the SO 2 emissions in surrounding areas played in affecting the air pollution in Beijing.These regions are all industrialized, urbanized and coal burning areas within a few hundred kilometers of Beijing.Therefore, SO 2 directly emitted from power plant, industrial facilities, breweries etc., can be easily transported to Beijing under the favorable meteorological conditions.
Although Beijing local emissions were found to be dominant contributions at the surface layer, it is also important to know how the surrounding areas affect the air quality at the surface layer in Beijing.Sectoral analysis at the surface layer was then performed Introduction

Conclusions References
Tables Figures

Back Close
Full to gain insights into the direct impact of the surrounding regions on SO 2 pollution over certain counties in Beijing.Beijing administrative subdivisions were used (Fig. 7) to evaluate this impact.Consistent with the results above, the Beijing local contribution was dominant except Pinggu County.However, TJ and ZJK were distinct from the surrounding areas with a significantly direct contribution of SO 2 to its closer counties in Beijing, in particular, 86% and 68% of the surrounding contribution on average at Pinggu County and Yanqing County respectively (Table 6).The contribution from TJ SO 2 emissions to Pinggu County in Beijing was much higher than that from BJ local emissions.Overall, without the Beijing local contribution, the impact of surrounding regions (e.g., TJ, CD, ZJK, BD and LF) on the SO 2 pollution in Beijing are mainly via Pinggu, Miyun, Yanqing, Fangshan and Tongzhou counties respectively (Fig. 7, Table 6).
Basically, such characteristic of SO 2 contribution is closely associated with Beijing geographical position, highlighting the performance of the air pollutant tagged method in estimating the contributions to air pollution over Beijing.The results specify the principal enter-ways of the typical influence of the key surrounding contributors on SO 2 pollution over Beijing.It was undoubtedly helpful for further SO 2 pollution control strategy especially for an efficient long-term County-by-County SO 2 emissions control program over Beijing.

The contribution of the PM 10 concentration to Beijing
The PM 10 from primary emissions were tagged and its contribution to the PM 10 the PM 10 pollution in Beijing were also calculated in this study.The contributions of the secondary particle matters (often refers to secondary aerosol) to the PM 10 contribution, however, were not included in this study.
Similar to SO 2 , the Beijing local contribution to the PM 10 concentration at the surface layer (23.5 µg m −3 ) was much higher than that at the upper layer (1.0 µg m −3 ) (Fig. 8).
The Beijing local sources contributed more than 75% PM 10 concentration in Beijing at the surface layer, while only 13% at the upper layer (about 1.1km).That indicated that 5282 Introduction

Conclusions References
Tables Figures

Back Close
Full the surrounding areas contributed 87% PM 10 concentration in Beijing at the upper layer, with the most contribution from Hebei province (43%).Thus, the regional transport played an important role in controlling the PM 10 pollution and visibility in Beijing at the upper layer.In addition, the regional PM 10 contribution from the southern surrounding areas appeared to be higher than those from the northern regions, which is similar to that of SO 2 , but with relatively lower contribution.
The daily variations of the contribution percentage from different regions to PM 10 in Beijing are presented in Fig. 9. Similar to SO 2 in Fig. 5, the PM 10 pollution in Beijing were dominated by the contributions from BD, LF,TJ and ZJK on 7,12,17,20,23 and 25 August respectively.This suggested a typical flux in regional air pollutant transport, including SO 2 and PM 10 .The finding was also consistent with the typical heavy pollution sequences described in Sect.2.5.2 (model evaluation).Such convergence reflected the air pollutant tagged method capability to accurately reproduce the air pollutants transport.On the other hand, the lower contribution of the long-range transport from Shandong, Inner-Mongolia and Shanxi Provinces to the PM 10 in comparison to that of SO 2 was likely due to the fact that tracer gases are more transportable than particulate matter which has larger dry deposition velocity (An et al., 2007).
Similar to SO 2 pollution in Beijing, county scale analysis of the PM 10 contribution without the local contribution is presented in Fig. 10 and the results are summarized in Table 6.Again, our results showed the significant PM 10 contributions from TJ to Pinggu County (63% of 12.8 µg m −3 ), from LF to Daxing (76% of 8.7 µg m −3 ), Tongzhou (71% of 7.3 µg m −3 ) and Shunyi (68% of 5.0 µg m −3 ) counties, and from BD to Fangshan County (80% of 4.5 µg m contributor to urban Beijing at the surface layers was the pair LF and BD with 72% of the surrounding contributed PM 10 concentration. Our results are also consistent with Streets et al. (2007) analysis on PM 2.5 , but more detailed, despite the different methods, simulations periods and the meteorological conditions between these two studies.Given that PM 10 is often the primary pollutant in Beijing (Hao et al., 2008;An et al., 2007), the sectoral analysis in this study brought out the substantial information for relevant environmental implications by specifying areas of direct influence of each potential contributor.Our results support the strategies that controlling the PM 10 pollution in Beijing should be implemented at the county scale.

The contribution of the NO 2 and ozone (O 3 ) concentration to Beijing
Nitrogen oxides (NO x = NO + NO 2 ) play a key role in tropospheric chemistry, e.g., the production of ozone in the troposphere (Seinfeld and Pandis, 2006;Zhang et al., 2009).
The contribution from the surrounding areas to NO 2 pollution in Beijing showed overall similar patterns to those of PM 10 and SO 2 , that the percentage of the surrounding contribution at the upper layer were much higher than those at the surface layer, essentially from the southern regions (Fig. 11), while the contributions at the upper layer was lower than those at the surface layer.The contribution from the regional transport accounted for 10% of NO 2 at the surface layer, while up to 40% at the upper layer, which are all lower than those of SO 2 and PM 10 .Yet, TJ, BD and ZJK were still the most significant surrounding contributors, accounting for 11%, 9% and 8% respectively to the NO 2 in Beijing at the upper layer.
The county scale analysis of the NO 2 contribution is presented in Fig. 12.The local contribution was dominant in all counties of Beijing, including Pinggu County, which is slightly different from the SO 2 .The urban Beijing showed the highest NO 2 concentration of 36.13 ppbv (36.53 ppbv-0.40ppbv) due to the contribution from local NO x emissions.Without considering the contribution from Beijing local sources, TJ was the most important regional contributor to Pinggu County with an average contribution of 84% (monthly mean 3.99 ppbv); comparatively, LF was the primary contributor to Daxing County accounting for 68% of 2.66 ppbv.TJ and LF together contributed 97% of 5284 Introduction

Conclusions References
Tables Figures

Back Close
Full the surrounding NO 2 concentration (2.63 ppbv) to Tongzhou County (Fig. 12 and Table 6).Similar to PM 10 , the surrounding areas showed more impact on NO 2 pollution in the eastern counties, especially in Pinggu, Daxing and Tongzhou counties with the contribution range of 2.63 ppbv∼3.99ppbv.Therefore, the NO 2 pollution over Beijing was primarily from local NO x emission, consistent with the results reported previously (Streets et al., 2007;Zhang et al., 2009).However the environmental strategy in the eastern counties should also take efforts to control the NO x emission in the regions close to Beijing.
The ozone was also tagged after being produced by the photochemical reaction, and the results are shown in Fig. 13.It appeared that the contributions of ozone from both Beijing local sources and regional transport were strongly associated with the meteorological conditions and compatible with the present analysis.The mean monthly and maximum daily contribution to Beijing is presented in Table 7.The local production of ozone concentration to urban Beijing was about 53%, with maximum daily percentage of 75% at the surface layer.Similar to other species, the ozone pollution at surface layer is dominated by the local production.However, the contribution from surrounding areas at the surface layer was estimated to be about 47%, indicating a considerable contribution from the regional transport, especially from the southern regions.Without considering the local production, BD was the principal direct contributor to ozone pollution in urban Beijing, with a monthly average of 19% and maximum daily percentage of 44%.In total, 22 days with the contribution from BD exceeding 5% were observed.
Overall, a significant amount of ozone produced in surrounding regions can be transported to Beijing, which has large influence on the pollution level of NO x and VOCs.Our results also highlight the crucial role of transport and other meteorological conditions in photochemical pollution of ozone as previously reported by Streets et al. (2007), Hao et al. (2008) and Tang et al. (2010).The contribution of NO 2 from surrounding areas to Beijing discussed above were relatively smaller compared to that of ozone pollution.Further study on the transport of ozone precursors (NO x and VOCs) was therefore needed to properly characterize the ozone contribution.Introduction

Conclusions References
Tables Figures

Back Close
Full

Conclusions
The NAQPMS model coupled with an on-line air pollutant tagged module has been developed to investigate the contributions of local sources and regional transport to the air pollution in Beijing during CAREBeijing-2006.An in-depth analysis through this tagged method, with tagged 22 emissions regions, was performed to quantify the influence of various emissions regions on SO 2 , PM 10 , NO 2 and ozone pollution over Beijing.In addition to improved estimation of the surrounding areas contribution by separating pollutants transport layers, another important result from this study was the detailed information on the contribution from different surrounding regions to the air pollution over Beijing at county scale.Results were summarized as following: -The NAQPMS model showed highly improved model performance in simulating SO 2 , PM 10 , NO 2 and O 3 as suggested by the overall good agreements between the simulated and the observed data provided by CAREBeijing-2006.The model had better performance at the IAP tower station in the urban Beijing than Yufa and Xinglong station, which are located at the rural areas.And the simulated SO 2 at Xinglong station in the east of Beijing was much lower than the observed, likely due to the underestimated source emissions in this region.
-Analysis indicated that the local sources were dominant contributions at the surface layer, accounting for 65% of SO 2 concentration.On the episode day of 25 August, there was a significant contribution (about 25%) from TJ.The local contribution at the upper layer was lower than that at the surface layer.However, the surrounding areas, especially the southern regions, played a dominant contribution to the pollution at the upper layer, for example, 50% of the total SO 2 .
County scale analysis indicated that the regional transport from TJ contributed the most to SO 2 over Pinggu County, accounting for 86% of the SO 2 from the total surrounding areas.Introduction

Conclusions References
Tables Figures

Back Close
Full -Similar to SO 2 , the Beijing local contribution to the PM 10 at the surface layer (75%) was much higher than that at the upper layer.At the upper layer, the contributions from surrounding areas were more significant.The PM 10 contribution from the southern surrounding areas also seemed to be higher than those from the northern regions.In addition, the contributions from surrounding areas showed more impacts on PM 10 in eastern counties of Beijing than the western counties.
The most important surrounding contributor to Beijing at the surface layer was LF, which contributed more than 50% of the surrounding PM 10 contribution over Daxing, Tongzhou and Shunyi counties.
-The local contribution to the NO 2 concentration at the surface layer was dominant in all counties of Beijing, including Pinggu County, while the surrounding areas contributed 10% NO 2 concentration versus 40% at the upper layer, much lower than those of SO 2 and PM 10 .The surrounding areas also contributed more to the NO 2 pollution in the eastern counties, especially in Pinggu, Daxing and Tongzhou counties.
-The local production of ozone contributed about 53% of total ozone at urban Beijing at the surface layer.The contribution from surrounding areas was also significant and estimated to be about 47%, of which BD was the primary contributor.
In summary, the local contribution played the actual key role in controlling the air pollution over Beijing at the surface layer, while the surrounding areas contribution to air pollution in Beijing were more significant at the upper layer.Results in this study suggest that long-term efficient air quality control strategies in Beijing should be based on regional-scale collaborations, by not only abating the local emission, but also controlling emissions from the surrounding areas.However, achievement of sustainable ozone control over Beijing needs further long-term study on ozone sensitivity to precursor emissions distribution over Beijing and concerned surrounding areas.Introduction

Conclusions References
Tables Figures

Back Close
Full Fractional Bias: Supplementary material related to this article is available online at:  Full  Full Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

S
T ,F S T = S T /S 1 ): defining Beijing and Tianjin Municipalities, Shandong and Liaoning Provinces as four tagged individual regions, Hebei Province divided into eight regions, Inner Mongolia into five regions and Shanxi Province in four regions, and the resting neighboring areas are defined as one tagged emisions region and labeled into the NAQPMS model.Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | −3 ).The surrounding areas contributed more to the PM 10 pollution over the eastern counties (7.3 µg m −3 ∼12.8 µg m −3 ) of Beijing than the western (3.1 µg m −3 ∼3.9 µg m −3 ) counties, which might be associated with Beijing's semi-basin topography and the weather patterns of the summer monsoon.And the most regional contribution to PM 10 in Beijing at the surface layer is from LF (1.4 µg m −3 , Fig. 8), accounting for more than 50% of the surrounding contributed PM 10 concentration over Daxing, Tongzhou and Shunyi counties as mentioned above.The most important Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Fig. 4 .
Fig. 4. Temporal variations of observed and simulated SO 2 , NO 2 and Ozone (O 3 ) concentration at Xinglong and Yufa station in Aug 2006 25 Fig. 4. Temporal variations of observed and simulated SO 2 , NO 2 and Ozone (O 3 ) concentration at Xinglong and Yufa station in August 2006.

3 at Yufa Station
O