Exploring 2016–2017 surface ozone pollution over China: source contributions and meteorological influences

Severe surface ozone pollution over major Chinese cities has become an emerging air quality concern, raising a new challenge for emission control measures in China. In this study, we explore the source contributions to surface daily maximum 8 h average (MDA8) ozone over China in 2016 and 2017, the 2 years with the highest surface ozone averaged over Chinese cities in record. We estimate the contributions of anthropogenic, background, and individual natural sources to surface ozone over China using the GEOS-Chem chemical transport model at 0.25×0.3125 horizontal resolution with the most up-to-date Chinese anthropogenic emission inventory. Model results are evaluated with concurrent surface ozone measurements at 169 cities over China and show generally good agreement. We find that background ozone (defined as ozone that would be present in the absence of all Chinese anthropogenic emissions) accounts for 90 % (49.4 ppbv) of the national March–April mean surface MDA8 ozone over China and 80 % (44.5 ppbv) for May– August. It includes large contributions from natural sources (80 % in March–April and 72 % in May–August). Among them, biogenic volatile organic compound (BVOC) emissions enhance MDA8 ozone by more than 15 ppbv in eastern China during July–August, while lightning NOx emissions and ozone transport from the stratosphere both lead to ozone enhancements of over 20 ppbv in western China during March–April. Over major Chinese city clusters, domestic anthropogenic sources account for about 30 % of the May–August mean surface MDA8 ozone and reach 39– 73 ppbv (38 %–69 %) for days with simulated MDA8 ozone >100 ppbv in the North China Plain, Fenwei Plain, Yangtze River Delta, and Pearl River Delta city clusters. These high ozone episodes are usually associated with high temperatures, which induce large BVOC emissions and enhance ozone chemical production. Our results indicate that there would be no days with MDA8 ozone >80 ppbv in these major Chinese cities in the absence of domestic anthropogenic emissions. We find that the 2017 ozone increases relative to 2016 are largely due to higher background ozone driven by hotter and drier weather conditions, while changes in domestic anthropogenic emissions alone would have led to ozone decreases in 2017. Meteorological conditions in 2017 favor natural source contributions (particularly soil NOx and Published by Copernicus Publications on behalf of the European Geosciences Union. 8340 X. Lu et al.: Exploring 2016–2017 surface ozone pollution over China BVOC ozone enhancements) and ozone chemical production, increase the thermal decomposition of peroxyacetyl nitrate (PAN), and further decrease ozone dry deposition velocity. More stringent emission control measures are thus required to offset the adverse effects of unfavorable meteorology, such as high temperature, on surface ozone air quality.

Abstract. Severe surface ozone pollution over major Chinese cities has become an emerging air quality concern, raising a new challenge for emission control measures in China. In this study, we explore the source contributions to surface daily maximum 8 h average (MDA8) ozone over China in 2016 and 2017, the 2 years with the highest surface ozone averaged over Chinese cities in record. We estimate the contributions of anthropogenic, background, and individual natural sources to surface ozone over China using the GEOS-Chem chemical transport model at 0.25 • ×0.3125 • horizontal resolution with the most up-to-date Chinese anthropogenic emission inventory. Model results are evaluated with concurrent surface ozone measurements at 169 cities over China and show generally good agreement. We find that background ozone (defined as ozone that would be present in the absence of all Chinese anthropogenic emissions) accounts for 90 % (49.4 ppbv) of the national March-April mean surface MDA8 ozone over China and 80 % (44.5 ppbv) for May-August. It includes large contributions from natural sources (80 % in March-April and 72 % in May-August). Among them, biogenic volatile organic compound (BVOC) emis-sions enhance MDA8 ozone by more than 15 ppbv in eastern China during July-August, while lightning NO x emissions and ozone transport from the stratosphere both lead to ozone enhancements of over 20 ppbv in western China during March-April. Over major Chinese city clusters, domestic anthropogenic sources account for about 30 % of the May-August mean surface MDA8 ozone and reach 39-73 ppbv (38 %-69 %) for days with simulated MDA8 ozone > 100 ppbv in the North China Plain, Fenwei Plain, Yangtze River Delta, and Pearl River Delta city clusters. These high ozone episodes are usually associated with high temperatures, which induce large BVOC emissions and enhance ozone chemical production. Our results indicate that there would be no days with MDA8 ozone > 80 ppbv in these major Chinese cities in the absence of domestic anthropogenic emissions. We find that the 2017 ozone increases relative to 2016 are largely due to higher background ozone driven by hotter and drier weather conditions, while changes in domestic anthropogenic emissions alone would have led to ozone decreases in 2017. Meteorological conditions in 2017 favor natural source contributions (particularly soil NO x and

Introduction
Ozone near the surface is a major air pollutant that harms human health and vegetation growth (Monks et al., 2015). China has become a global hot spot of urban ozone pollution Lu et al., 2018a). The present-day (2013The present-day ( -2017 summertime fourth-highest daily maximum 8 h average (MDA8) ozone levels over eastern China are typically higher than 100 ppbv, inducing significantly larger human health and vegetation damage than those in other industrialized regions (i.e., Japan, Korea, Europe, and the United States) as estimated by different ozone exposure metrics (Fleming et al., 2018;Lu et al., 2018a). In particular, 2016 and 2017 witnessed the most severe surface ozone pollution in record over most Chinese cities. The summertime surface MDA8 ozone levels in 2016-2017 averaged over the major Chinese cities increased by about 20 % compared to 2013-2014 , despite the fact that the Chinese Action Plan on Air Pollution Prevention and Control implemented in 2013 has significantly reduced anthropogenic emissions of nitrogen oxides (NO x = NO + NO 2 ) and alleviated winter haze pollution (Cheng et al., 2019;. Exploring source attributions to surface ozone in China is therefore critical to understanding the recent ozone increase and to reexamining the air pollution control strategy. Surface ozone is mainly produced by the sunlight-driven photochemical oxidation of volatile organic compounds (VOCs) and carbon monoxide (CO) in the presence of NO x . These ozone precursors are emitted intensively in China from anthropogenic sources, including emissions from industry, plant powers, residential areas, and transportation (Li et al., 2017a;Zheng et al., 2018;Wang et al., 2019). Significant ozone production driven by intense anthropogenic emissions has been observed and modeled in many urban regions in China (Wang et al., 2006;Ran et al., 2009;Xue et al., 2014;Li et al., 2017b;Tan et al., 2018a). While anthropogenic contribution can be reduced by emission control measures, background ozone, defined as ozone concentration in the absence of domestic anthropogenic sources, presents a base below which ozone air quality cannot be improved by domestic emission controls. Background ozone includes ozone produced from natural emissions or originated from the stratosphere (together referred to as natural background ozone) and ozone produced from foreign anthropogenic emissions via long-range transport, as described in previous studies in the context of North American background ozone (McDonald-Buller et al., 2011;Zhang et al., 2011;Jaffe et al., 2018). Both natural emissions (e.g., from the biosphere or lightning) and transport processes are sensitive to meteorological conditions, such as temperature and wind, and can therefore significantly alter the spatiotemporal variability of ozone (Lin et al., 2008;Jacob and Winner, 2009;Fiore et al., 2012).
Source contributions are not observable and are therefore preferably quantified through modeling studies. Li et al. (2017a), using regional chemical models, estimated that anthropogenic emissions from industry, residential areas, and transportation together contributed up to 80 ppbv ozone (more than half of total ozone) during a severe ozone pollution episode over eastern China in May 2015. For the seasonal and nationwide average, however, background sources typically have a larger contribution. Wang et al. (2011) estimated that annual background sources contributed 89 % (44.1 ppbv) of mean surface ozone over China in 2006. Ni et al. (2018 also showed a large contribution from background sources to springtime surface ozone in 2008, with more than 40 ppbv over the North China Plain and 70 ppbv over western China. Long-range transport of ozone from Southeast Asia, Europe, and North America can enhance regional surface ozone by more than 2 ppbv in China (Ni et al., 2018). Several recent studies also examined the contribution from biogenic VOC emissions to regional surface ozone (Li et al., 2018a;Mo et al., 2018). No previous studies so far, however, provide a comprehensive view of anthropogenic vs. individual natural source contributions to present-day surface ozone levels in China. Their contributions to recent surface ozone increases over China remain unexplored.
In this study, we explore the sources contributing to surface ozone over China in 2016-2017 using the highresolution (0.25 • latitude × 0.3125 • longitude) GEOS-Chem chemical transport model with the most up-to-date yearspecific Chinese anthropogenic emission inventory to interpret nationwide hourly ozone observations over Chinese cities. We estimate ozone contributions from background, natural, and anthropogenic sources and further decompose the individual contributions from major natural sources. These include biogenic emissions, soil emissions, lightning emissions, biomass burning emissions, and stratospheric ozone transport. The source attributions and their interactions with meteorology are analyzed both seasonally from a nationwide perspective and at a daily scale in individual cities for high ozone episodes. We further compare model results for the 2 years (2016 vs. 2017) to investigate changes in source attributions contributing to recent ozone changes over China. We obtain nationwide hourly ozone observations from the China National Environmental Monitoring Center (CNEMC) (http://www.cnemc.cn/en/, last access: 22 June 2019). The observational network became operational in 2013, and by 2017 it was monitoring six surface air pollutants (including ozone and particulate matter with an aerodynamic diameter less than or equal to 2.5 µm, PM 2.5 ) in 454 cities (in total 1597 non-rural sites) in mainland China (Fig. 1). The datasets have been widely used to examine urban air quality issues over China in recent studies Li et al., 2017b;Gao et al., 2018;Shen et al., 2019;Li et al., 2019). Quality controls to remove unreliable hourly observations are applied following our previous work .
The Chinese Ministry of Environmental Protection (now the Chinese Ministry of Ecology and Environment) categorized a total of 74 major cities (provincial capital cities and major cities in the North China Plain (NCP), Yangtze River Delta (YRD), and Pearl River Delta (PRD); these are three well-developed city clusters in China with severe air pollution) as State I key cities for air pollution monitoring in 2013. The number of key cities for air pollution monitoring has expanded to 169 since July 2018 (CNEMC, 2018). The list not only includes more cities in the NCP (now 55 cities), YRD (41 cities), and PRD (9 cities), but also covers regions where air pollution is emerging. As shown in Table 1 and Fig. 1, they are the Fenwei Plain city cluster (FWP, 11 cities), Chengdu-Chongqing city cluster (or Sichuan Basin, SCB, 16 cities), central Yangtze River city cluster (CYR, 22 cities), and other provincial capital cities in northeastern China (NECH, 4 cities), western China (WCH, 7 cities), and southern China (SCH, 4 cities). We focus the analyses on these 169 cities grouped into 9 city clusters (Table 1) in this study. Each city contains several monitoring sites, and we average them hourly to represent air quality at the city level.

Model description
We use the GEOS-Chem chemical transport model (CTM) (v11-02rc; http://geos-chem.org, last access: 22 June 2019) to interpret surface ozone measurements in China. The model is driven by assimilated meteorological data obtained from the Goddard Earth Observing System (GEOS) of the NASA Global Modeling and Assimilation Office (GMAO). The GEOS-FP datasets are available at a native horizontal resolution of 0.25 • latitude by 0.3125 • longitude, a temporal resolution of 1 h for surface variables, and boundary layer height and 3 h for others. We use the nested-grid version of GEOS-Chem that has 0.25 • × 0.3125 • horizontal resolution over East Asia (70-140 • E, 15-55 • N) (Chen et al., 2009;Zhang et al., 2016;Shao et al., 2019), with boundary conditions archived from the global simulation at 2 • latitude × 2.5 • longitude resolution.
The model includes a detailed mechanism of tropospheric HO x -NO x -VOC-ozone-halogen-aerosol chemistry (Wang et al., 1998;Bey et al., 2001;Park et al., 2004;Mao et al., 2013). The chemical kinetics are obtained from the Jet Propulsion Laboratory (JPL) and International Union of Pure and Applied Chemistry (IUPAC) (Sander et al., 2011;IU-PAC, 2013). Photolysis rates are calculated using the Fast-JX scheme (Bian and Prather, 2002). The linearized ozone parameterization (LINOZ) is used to simulate stratospheric ozone as described by McLinden et al. (2000). Other species in the stratosphere are calculated in the model based on archived monthly mean production and loss rates provided by the Global Modeling Initiative (GMI; Murray et al., 2013). Dry deposition of both gas and aerosols are calculated online based on the resistance-in-series algorithm (Wesely, 1989;Zhang et al., 2001). The wet deposition for water-soluble aerosols and gas in GEOS-Chem is described by  and Amos et al. (2012). The non-local scheme for the boundary layer mixing process and the relaxed Arakawa-Schubert scheme for cloud convection are described in Lin and McElroy (2010) and Moorthi and Suarez (1992), respectively.
Emissions in GEOS-Chem are processed through the Harvard-NASA Emission Component (HEMCO) (Keller et al., 2014). Global anthropogenic emissions in this study are from the Community Emissions Data System (CEDS; the latest 2014 condition is used for the model simulation) (Hosely et al., 2018) overwritten by regional emission inventories over the US (National Emission Inventory, NEI), Canada (Canadian Criteria Air Contaminant), Mexico (Kuhns et al., 2005), Europe (European Monitoring and Evaluation Program, EMEP), Africa (DICE-Africa inventory) (Marais and Wiedinmyer, 2016), and East Asia and South Asia (MIX inventory) (Li et al., 2017a). In particular, we apply the latest Chinese anthropogenic emissions for 2016 and 2017 from the Multi-resolution Emission Inventory for China (MEIC; http: //www.meicmodel.org, last access: 22 June 2019) . The MEIC is a bottom-up emission inventory with particular improvements in the accuracy of unit-based power plant emission estimates , vehicle emission modeling (Zheng et al., 2014), and the NMVOC speciation method . The bimonthly spatiotemporal distributions of anthropogenic NO x , CO, and non-methane VOC (NMVOC) emissions averaged for 2016-2017 are shown in Fig. S1 in the Supplement. The highest emissions are in the populated city clusters in central eastern China (typically includes NCP, FWP, YRD, CYR, and adjacent regions) with little seasonal variation. Total annual Chinese anthropogenic emissions of NO x , NMVOC, and CO are, respectively, 22.5, 28.4, and 141.9 Tg in 2016and 22.0, 28.6, and 136.2 Tg in 2017. The annual Chinese NO x and CO emissions both decrease by approximately 20 % in 2017 compared to 2013 when the Chinese State Council initiated the Action Plan on Air Pollution Prevention and Control, but NMVOC emissions show a slight increase of 2 % . All anthropogenic emissions are emitted at the lowest layer in the model with a thickness of 120 m. A recent study found that uncertainties in industrial emission injection heights could affect surface SO 2 concentrations by about 80 % ; this may further affect surface ozone concentrations and requires further study. GEOS-Chem has implemented a number of natural emissions, some of which are calculated online in the model. Lightning NO x emissions are parameterized as a function of cloud-top height (Price and Rind, 1992), vertically distributed following Ott et al. (2010), and spatially constrained by climatological observations of lightning flash rates from the Lightning Imaging Sensor (LIS) and the Optical Transient Detector (OTD) satellite instruments (Murray et al., 2012). Following previous studies (Hudman et al., 2007;Zhang et al. 2014), the amount of NO released per flash is 500 moles for the lightning north of 35 • N in Eurasia and 23 • N in North America; it is 260 moles for the rest of the world. Soil NO x emissions are calculated based on nitrogen (N) availability in soil and edaphic conditions such as soil temperature and moisture as described in Hudman et al. (2010Hudman et al. ( , 2012. Biogenic VOC (BVOC) emissions are calculated following the Model of Emissions of Gases and Aerosols from Nature (MEGAN version v2.1) algorithm (Guenther et al., 2012). Monthly mean biomass burning emissions for the year 2014 from the Global Fire Emissions Database version 4 (GFED4) are used in the simulation (van der Werf et al., 2017). Bimonthly mean total emissions of these natural processes averaged over China are summarized in Table 2. Mixing ratios of methane are prescribed in the model based on spatially interpolated monthly mean surface methane observations from the NOAA Global Monitoring Division for 1983-2016 and are extended to 2020 using the linear extrapolation of local 2011-2016 trends (Murray, 2016). We find that it leads to an approximate 1 % increase in methane concentration over eastern China in 2017 relative to 2016.

Model configurations
Model configurations are summarized in Table 3. The standard simulation (BASE) includes all anthropogenic and natural emissions as described above. We then conduct sensi-  (1) with all anthropogenic emissions (all emitted pollutants except methane, which is prescribed in all simulations here) turned off over China (noCH) and (2) with all anthropogenic emissions turned off globally (noGLOBE). Ozone concentrations in noCH therefore represent total background ozone in China (hereafter called background ozone), and those in noGLOBE represent natural background ozone (hereafter called natural ozone). The differences in ozone concentration between BASE and noCH denote the domestic anthropogenic ozone enhancements (CH anthropogenic ozone). We also conduct several additional sensitivity simulations by turning off individual sources to estimate their ozone enhancements in the presence of all other sources. The ozone contribution from each source can therefore be estimated as the ozone difference between the BASE simulation and each sensitivity simulation. This difference quantifies the source contributions to present-day surface ozone, which includes the interaction of each specific source with all other sources rather than the pure impact from the specific source alone (Li et al., 2018a). Four such sensitivity simulations are conducted by turning off (1) BVOC emissions (noBVOC), (2) lightning NO x emissions (noLIGHT), (3) soil NO x emissions (noSOIL), and (4) biomass burning emissions (noBB).
Due to the long lifetime of ozone in the lower stratosphere (approximately years; Wang et al., 1998), the stratospheric contribution to tropospheric ozone is typically quantified using the tagged ozone simulation (TagO 3 ) (Wang et al., 1998;Zhang et al., 2014) instead of perturbing the stratospheric ozone chemistry, which requires a long spinup time. The tagged ozone simulation labels ozone produced above the tropopause from the photolysis of molecular oxygen as stratospheric ozone and then simulates its transport and chemical loss in the troposphere . This tagged stratospheric ozone is calculated based on ozone production and loss frequency archived from the BASE simulation.
For all simulations (except for the TagO 3 simulation, which was spun up for 3 years), the global simulation (2 • latitude × 2.5 • longitude) was first conducted from April 2015 to November 2017. The results on 1 February 2016 were then interpolated to high resolution (0.25 • latitude × 0.3125 • longitude) over the nested domain and were used to initialize the nested model simulation. Results from the nested model for March-October in both 2016 and 2017 are analyzed. The simulation labels stratospheric ozone (ozone produced in the stratosphere from the photolysis of molecular oxygen) as a tagged tracer. The simulation is driven by ozone production rates and loss frequencies archived from the BASE simulation.
3 Sources contributing to surface ozone pollution in We use MDA8 throughout the analysis as it is the ozone air quality standard in China and also an important metric of human health exposure (Turner et al., 2016). The model reproduces the spatial distribution of observed surface MDA8 ozone in the warm season (May-August) with a high spatial correlation coefficient (r) of 0.81-0.82 and a relatively small positive mean bias of 4.5-7.4 ppbv (7 %-13 %) in the 169 cities. In particular, the model captures the ozone hot spots over eastern China (including the NCP, YRD, FWP, and CYR city clusters) in May-June when both observations and model results show MDA8 ozone > 60 ppbv in a number of cities (Fig. 2). In spring (March-April) and autumn (September-October), model results have larger positive biases (5.3-9.2 ppbv, 12 %-20 %) while still capturing the higher surface ozone in eastern China (r = 0.32-0.46 for all 169 cities). The simulated daily MDA8 values are in good agreement with the observed values, with r ranging from 0.53 to 0.75 in the nine city clusters (Fig. S2). The model captures 57 % of the days with observed MDA8 > 75 ppbv (approximately the grade II Chinese national air quality standard for the protection of residential areas) for all nine city clusters (and more than 60 % for NCP, YRD, and FWP cities). The model also reproduces the seasonal variation of surface ozone over mainland China, i.e., the peak in summertime and lower values in other months (Fig. 2). The seasonality of surface ozone typically reflects the dominant role of temperature and solar radiation in enhancing natural emissions of ozone precursors and accelerating photochemical ozone production under high precursor concentrations, with modulations by the arrival of the summer monsoon (Ding et al., 2013;Li et al., 2018b;Lu et al., 2018a).
Surface ozone levels at Chinese urban sites typically show large day-to-day variations driven by intense regional anthropogenic emissions, rapid chemical conversion, and shifts of meteorological patterns (Wang et al., 2006;Lu et al., 2010;Tan et al., 2018b). We find in Fig. S3 that observed MDA8 ozone concentrations in many cities can shift by more than 60 ppbv in 2 d. Such variations, however, are challenging for Eulerian chemical models to capture due to numerical diffusion, even at 10 km grid resolution, and representative issues when comparing gridded simulated results to observations at urban sites (Young et al., 2018). Artificial mixing of ozone precursors in model grids may cause higher ozone production efficiency and therefore positive ozone biases (Wild and Prather, 2006;Yu et al., 2016;Young et al, 2018). In addition, models have limited skills in presenting local emissions and meteorological conditions, particularly over complex terrain (e.g., in WCH and SCB cities) Ni et al., 2018). These limitations largely explain the ozone overestimates at Chinese urban sites and the deficiencies in capturing extremely high ozone values in our model, as well as many other modeling studies using GEOS-Chem at coarser resolutions (e.g., Wang et al., 2011;Ni et al., 2018) and regional models such as the Community Multiscale Air Quality (CMAQ) model (Chen et al., 2018;Liu et al., 2018). We find that by increasing the resolution from 2 • × 2.5 • to 0.25 • × 0.3125 • , model results show significant improvements in capturing the ozone daily variability and summertime high ozone episodes (model vs. observation cor-  (Fig. 3b). Natural sources contribute on average 80 % (44.1 ppbv) of the surface MDA8 ozone in March-April and 71 % (39.6 ppbv) in May-August (Fig. 3c). Higher background ozone concentrations are distributed in western and northern China where both natural ozone and ozone imported from foreign anthropogenic sources have larger influences than in eastern and southern China. The peak background contribution in spring reflects a mixed effect of larger transboundary transport from foreign anthropogenic sources driven by strong westerly winds HTAP, 2010;Ni et al., 2018) and larger transport from the stratosphere, as will be discussed in Sect. 3.3. In contrast, ozone contributions from domestic anthropogenic sources are much smaller when averaged over China, with a peak contribution of about 24 % (12.1 ppbv) in July-August and less than 15 % (5-9 ppbv) in other months (Fig. 3d). The largest anthropogenic ozone enhancements are spatially centered in central eastern China following the hot spots of anthropogenic emissions (Fig. S1) and population density; they seasonally peak in summer due to active photochemistry at higher temperature and solar radiation conditions (Fig. S4). There is over 30 ppbv of ozone enhancement by domestic anthropogenic sources in central eastern China during July-August, which is comparable to background ozone there. Figure 3e shows the spatial distribution of foreign anthropogenic ozone enhancements estimated as the difference between the background ozone (Fig. 3b) and natural ozone (Fig. 3c). Larger contributions are found near the national borders and in the springtime. The foreign anthropogenic contributions we estimate here (∼ 5 ppbv for all seasons averaged over China), however, may underestimate their true contribution, since they are derived in the absence of Chinese domestic anthropogenic emissions and thus do They used daily mean ozone as a metric and found that background ozone accounted for about 93 % of total ozone in March-May and 80 % in June-August. In good agreement with , our model results when converting to daily mean ozone as a metric estimate background ozone contributions of 95 % in March-April and 80 % in July-August. These are slightly higher than values derived from the MDA8 metric, reflecting higher domestic anthropogenic contributions to surface ozone at higher ozone levels as will be discussed later.

Identifying ozone enhancements from individual natural sources
We now illustrate in Fig. 4   The source attributions are diagnosed as the difference between the BASE simulation and a simulation with the specific source turned off, expect for ozone from the stratosphere, which is diagnosed by the TagO 3 simulation as described in the text. averaged over China (Fig. 4a). The enhancements are particularly large over central eastern China where BVOCs enhance surface MDA8 ozone by more than 10 ppbv in July-August. These values are considered to be the actual contributions of BVOCs to ozone, which include the pure contributions, i.e., contributions of BVOCs to ozone in the absence of all other sources and the interactional effects of BVOCs with other sources. Li et al. (2018a) decomposed interactional and pure contributions and found that the interactional contributions of BVOCs and anthropogenic emissions enhanced summertime (August 2011) surface ozone in urban Xi'an by 14.3 ppbv, which is comparable with our results in this city, while the pure BVOC contributions alone are only about 2.6 ppbv. The seasonal and spatial variations of BVOC ozone enhancements (Fig. 4a) are mainly driven by the exponential dependency of BVOC emissions on temperature. As shown in Figs. S4 and S5, high biogenic isoprene emissions occur over central eastern China in July-August when daily maximum 2 m air temperature (T MAX ) generally exceeds 30 • C. Previous studies have shown that estimates of BVOC ozone enhancements can be influenced by the model simulation of yields and fates of organic nitrates (RONO 2 ), a by-product from the oxidation of isoprene and its carbonyl (Ito et al., 2009;Fu et al., 2015). Here the GEOS-Chem model assumes a 9 % yield of RONO 2 from the reaction of isoprene peroxy radicals with NO, following recent laboratory observations (Paulot et al., 2012;Lee et al., 2014).
Lightning NO x emissions increase bimonthly mean surface MDA8 ozone by 6.5-9.9 ppbv averaged over China, with the largest contributions (typically more than 12 ppbv) found over the Tibetan Plateau (Fig. 4b). The large lightning ozone enhancements over western China (annual mean > 7 ppbv) were also simulated by Murray (2016), and these values were higher than those over the western US (annual mean of 3-5 ppbv). However, the model may overestimate springtime lightning ozone enhancements over the Tibetan Plateau as the model shows high surface ozone biases of more than 15 ppbv in spring over this region (Fig. 2). In central eastern China, lightning activities also enhance surface MDA8 ozone by about 2-8 ppbv for all seasons. The spatiotemporal patterns of lightning ozone enhancements, however, do not follow the lightning NO x emission patterns that typically peak in summer over central eastern China in the middle troposphere (600-400 hPa) driven by strong convection (Figs. S5-S6). The large lightning ozone enhancements for all seasons over the Tibetan Plateau and the Altai Mountains (in western Mongolia) are mainly due to their high elevations. Although lightning NO x emissions are larger in July-August, the shorter ozone lifetime and stronger upward transport over central eastern China in these months can suppress downward mixing of lightning ozone enhancements to the surface (Fig. S4), resulting in their minimum influences there in the period. This also partly explains the decreases in surface ozone from spring to summer in western China as shown in Fig. 2.
Stratosphere-troposphere exchange (STE) is typically active in the boreal spring at northern midlatitudes associated with synoptic-scale and mesoscale processes, such as tropopause folds, gravity wave breaking, and deep convection (Stohl et al., 2003). We show in Fig. 4c that stratospheric ozone contributions at the surface are 8.7 ppbv in March-April and 6.7 ppbv in May-June averaged for China. Western China, particularly the Tibetan Plateau, is strongly influenced by deep stratospheric intrusions with more than 20 ppbv of surface MDA8 ozone originating from the stratosphere, consistent with the estimate of Xu et al. (2018) using a different CTM. We also find important stratospheric influences (6-10 ppbv) on surface ozone over the NCP and FWP regions in March-June. Stratospheric contributions to surface ozone are at a minimum in summer when STE is weakest and upward transport is strongest, limiting downward mixing of stratospheric ozone.
Soil NO x emissions increase national mean surface MDA8 ozone by 1.6 to 3.1 ppbv (Fig. 4d). Figure S5 shows that the model simulates the highest soil NO x emissions in May-August over central eastern China and northern Asia, consistent with those derived from satellite observations (Vinken et al., 2014). Large emissions in these regions can be explained by the combined effect of high temperature and frequent pulsing emissions of soil NO x after rainfall . The largest soil ozone enhancements are not colocated with the emission hot spots but in regions such as northern China (2-6 ppbv in May-August) where surface ozone production is more sensitive to NO x as indicated by the modeled H 2 O 2 /HNO 3 values (Fig. S7). The lower bimonthly mean soil NO x ozone contributions in central eastern China compared to other regions are due to averaging some negative daily values, as will be discussed in Sect. 3.4.
Compared to other natural sources, biomass burning emissions lead to lower ozone enhancements averaged over China (less than 2 ppbv). This is because biomass burning emissions are generally small over China (less than 2 % of the total NO x emissions over China) (Wang et al., 2007;Lin, 2012). Larger influences are found during boreal spring (March-April) in southwestern China than other seasons and regions in China. The enhancements are particularly large in the southern part of Yunnan Province with over 15 ppbv of ozone enhancement driven by intense fires in Southeast Asia (Fig. S5d). Previous studies have shown that Eulerian chemical models such as GEOS-Chem tend to overestimate ozone enhancements from biomass burning emissions in regions adjacent to emission sources (i.e., in southwestern China) but underestimate them in downwind regions Lu et al., 2016). This can reflect a number of limitations, including inadequate representation of wildfire emission and chemistry (e.g., missed short-lived VOCs, fewer PAN emissions, and improper emission height), as well as the stretched-flow numerical diffusion of narrow plumes in the model even at fine resolutions Eastham et al., 2017). It should be noted that the total contributions of the above individual natural sources (sum of Fig. 4) together account for about half (from 47.4 % in September-October to 53.4 % in July-August) of natural ozone levels (Fig. 3c). The missing rest can be largely attributed to contributions from global methane. Fiore et al. (2008) found that 1 Tg of methane emissions would lead to an increase in global mean surface ozone by about 15 pptv (∼ 9 ppbv with present-day annual emissions of 600 Tg), and such sensitivity doubled at high NO x levels. Further analyses are therefore needed to quantify the contributions of global methane levels to background ozone in China. The discrepancies may also be due to nonlinear interactional effects between different sources that are not captured by sensitivity simulations with individual emissions turned off.

Sources contributing to surface ozone in city clusters
It is of particular importance to assess sources contributing to surface MDA8 ozone in major Chinese cities. Figure 5 shows the ozone contributions from background, anthropogenic, and natural sources in nine city clusters (Fig. 1, Table 1) averaged for March-October (Fig. 5a) and for May-August (Fig. 5b) (Fig. 5b), the domestic anthropogenic contribution fractions increase by about 6 %-9 % in all city clusters except PRD and reach 31 % in the NCP, 43 % in the YRD, and 32 % in the FWP, reflecting the increasing importance of domestic anthropogenic contributions to surface ozone in the peak ozone season. The domestic anthropogenic contribution fractions in PRD remain the same for the March-October and May-August averages as surface ozone there typically peaks in September-October (Fig. 2). In Fig. 6, we present the source contributions to daily MDA8 surface ozone in 2017 categorized by different surface ozone levels, which further illustrates the higher contributions of domestic anthropogenic sources to higher daily MDA8 ozone values in Chinese cities. For all city clusters, domestic anthropogenic contributions are much smaller (typically less than 15 % of the total simulated ozone) when simulated MDA8 ozone concentrations are lower than 50 ppbv. These lower contributions are more prevalent during nonsummer seasons when ozone is dominantly controlled by background ozone, except for PRD and SCH. Figure 6 also shows that domestic anthropogenic emissions can cause ozone decreases over central eastern China and northern China due to the NO x -saturated chemical regime in spring and fall (Fig. S7), leading to conditions with background ozone higher than simulated total ozone.
When simulated MDA8 ozone levels are higher than 75 ppbv, domestic anthropogenic contributions significantly increase to 28.8 ± 6.8 ppbv (35 %) in NCP, 27.8 ± 7.1 ppbv (33 %) in FWP, 38.8 ± 10.7 ppbv (46 %) in YRD, 37.7 ± 6.7 ppbv (46 %) in SCB, and 54.7 ± 11.8 ppbv (61 %) in PRD. For the extremely high ozone episodes (simulated MDA8 ozone > 100 ppbv), their contributions are as high as 42.0 ± 5.3 ppbv (41 %) in NCP, 39.1 ± 4.8 ppbv (38 %) in FWP, 61.4 ± 7.7 ppbv (58 %) in YRD, 52.4 ± 10.7 ppbv (48 %) in SCB, and 73.2 ± 7.0 ppbv (69 %) in PRD. We also find in Fig. 6 that there would be no days with MDA8 > 80 ppbv if there were no domestic anthropogenic influences in all 169 major Chinese cities. Our results show that background ozone has a dominant role in the monthly mean surface ozone, on top of which domestic anthropogenic contributions are of particular significance in ozone pollution episodes. Reducing anthropogenic contributions can therefore be effective in alleviating surface ozone pollution. One exception is the WCH region, where background ozone contributions are dominant for all seasons and even lead to high MDA8 ozone concentrations approaching 80 ppbv.
Figures 7-9 further show the time series of observed vs. simulated daily MDA8 ozone and the source contributions in representative cities from the three most polluted city clusters (Beijing in NCP, Shanghai in YRD, and Xi'an in FWP) in March-October 2017. The model is able to capture the daily variability of MDA8 ozone as well as the ozone peaks in Shanghai and Xi'an (r = 0.78 and 0.87, respectively). However, it fails to reproduce the high ozone episodes in Beijing during May-June 2017, although still with a high temporal correlation coefficient (r = 0.78), reflecting some model limitations as discussed in Sect. 3.1.
Analyses of the daily MDA8 ozone concentrations in the three cities indicate that both background and domestic anthropogenic ozone enhancements are significantly higher during high ozone episodes than clean days, consistent with Fig. 6. This may reflect their common correlations with meteorological parameters such as temperature. As shown in Figs. 7-9, the deseasonalized summertime (May-August) daily T MAX values are strongly correlated with both observed and simulated ozone, with r (T MAX vs. observed ozone) of 0.72 in Beijing, 0.54 in Shanghai, and 0.75 in Xi'an. We find that such positive temperature-ozone correlations (r > 0.60) are widespread in the cities over central eastern China in both observations and model results, as shown in Fig. 10. High summertime temperatures in China are typically associated with anticyclone systems (e.g., the West Pacific Subtropical High) that bring a clear sky with strong solar radiation and increased air stagnancy (Pu et al., 2017;. Such weather conditions tend to enhance ozone chemical production, weaken pollution ventilation, and favor natural emissions of ozone precursors. We find large temperaturedriven daily variations of BVOC ozone enhancements during May-August in all three cities (Figs. 7-9). Taking Shanghai (Fig. 8) as an example, observed and simulated high ozone episodes are usually associated with large BVOC ozone enhancements of over 20 ppbv. Nevertheless, in the absence of anthropogenic emissions, the total MDA8 ozone would be below 70 ppbv despite high temperatures and background ozone in those cities. Implementing emission control measures under high temperature conditions can be effective to avoid ozone exceedance in Chinese cities.
Other natural processes such as lightning, soil emissions, and stratospheric transport also contribute to MDA8 ozone daily variations and can episodically lead to large ozone enhancements. For example, in Xi'an, a stratospheric intrusion enhanced surface ozone by over 25 ppbv in early May 2017 and pushed the observed MDA8 ozone to approach 80 ppbv (Fig. 9). Surface ozone enhancements by lightning are episodically higher than 10 ppbv and show a similar temporal evolution as stratospheric ozone enhancements, as downward air motion delivers both ozone enhancements originating from high altitudes. Ozone enhancements from soil NO x emissions become notable when there are pulsing soil emissions after rainfall (e.g., ∼ 10 ppbv in mid-August in Xi'an). They can be either positive or negative, reflecting the monthly or even daily shifts of the ozone chemical production regime, as also reported from satellite observations (Jin and Holloway, 2015) and several field campaigns . Panel (b) shows the ozone enhancement from biogenic VOCs (BVOCs, black), lightning NO x (LITNO, red), soil NO x (SOILNO, blue), and the stratosphere (STRAT, purple). The estimation of ozone source attributions is described in the text. Also shown is daily maximum 2 m air temperature (TMAX, grey, right y axis) from the GEOS-FP data. The black pluses denote days with precipitation recorded from GEOS-FP data. in urban China (e.g., Lu et al., 2010). Biomass burning contributions to surface MDA8 ozone are generally small (less than 2 ppbv, not shown) in the three cities.

Factors contributing to the ozone increase from 2016 to 2017
The 2-year model results allow us to compare the surface ozone concentrations and their source contributions between 2016 and 2017. Here we focus on the warm season (May-August) when most cities in eastern China show high surface ozone levels (Fig. 2). We examine the differences in anthro-pogenic and natural emissions, meteorological conditions, and their roles in the surface ozone changes. As shown in Table 2, domestic anthropogenic emissions of NO and CO in May-August decline by 0.12 Tg (2.6 %) and 1.44 Tg (3.7 %) in 2017 compared to 2016, respectively. The decreases in NO and CO emissions are widespread over eastern China, with the largest reductions in Hebei, Shandong, and Jiangsu provinces, except in Beijing, Shaanxi, and Ningxia provinces where both NO and CO emissions show slight increases ( Fig. 11a and b). On the other hand, anthropogenic NMVOC emissions increase by 0.1 Tg (1.1 %) averaged over China, Figure 9. Same as Fig. 7 but for Xi'an. which is up to 4 % in the NCP (except in Hebei Province) and FWP regions (Fig. 11c). Figure 12 presents differences in May-August mean meteorological conditions between the 2 years. T MAX in central eastern China averaged for May-August is distinctly higher (> 2 K) in 2017 than 2016, together with drier air and more frequent clear sky as indicated by lower specific humidity and cloud cover fraction in 2017. Such warmer and drier weather conditions in 2017 are favorable for higher natural emissions of BVOCs and soil NO x . We can see that the model calculates larger biogenic isoprene (1.36 Tg larger) and soil NO (0.1 Tg larger) emissions in May-August 2017 than those in 2016 over China (Table 2), mainly over eastern China (Fig. 11d and f) following the spatial pattern of temperature differences. We find that regionally in eastern China the magnitude of soil NO emission enhancement exceeds the decrease in anthropogenic NO emissions, and their total changes over China are also comparable ( Table 2). Lightning NO x emissions in 2017 are lower than those in 2016 (Fig. 11e).
We show in Fig. 13 the observed and simulated May-August mean surface MDA8 ozone differences between 2016 and 2017. Observed May-August mean surface ozone concentrations show significant increases in 2017 in the cities over central eastern China (Fig. 13a). The regional increases are 7.5 ± 6.7 ppbv (mean ± standard deviation) averaged for NCP cities, 12.3 ± 6.7 ppbv for YRD cities, and 12.2 ± 12.2 ppbv for FWP cities. Several cities in PRD show ozone decreases during this period. Although the positive changes are widely distributed in eastern China, there are large variabilities among cities, indicating some local characteristics of ozone changes. Model results with year-specific meteorology and Chinese anthropogenic emissions successfully reproduce the pattern of ozone increases in central eastern China, yet the magnitudes (about 2-6 ppbv) are significantly lower than observations (Fig. 13a). This is likely due to the model limitations, as discussed in Sect. 3.1, and uncertainties in meteorological inputs and anthropogenic emissions. In western China and PRD cities, model results show that ozone decreases by about 2 ppbv in 2017 relative to 2016. Figure 13 also shows the surface ozone changes contributed from different sources as estimated by the model. We find that changes in background ozone alone lead to ozone increases in eastern China, particularly in the NCP, FWP, and YRD regions, and decreases in western China (Fig. 13b). The changes in background ozone concentrations are in good agreement with the changes in total ozone for both the spatial pattern and magnitude and are therefore identified as the main driver of surface ozone changes in 2017 from the 2016 levels. The changes in background ozone are entirely due  to changes in meteorological conditions between the 2 years (except for slight increases in prescribed methane concentrations that are not addressed in the present study). By contrast, changes in domestic anthropogenic emissions generally lead to slight ozone decreases in central eastern China and in the PRD, except for some regions in the Shandong, Hubei, Fujian, and Sichuan provinces (Fig. 13c).
The higher background ozone concentrations in China in May-August 2017 than 2016 can be partly explained by higher contributions from natural sources. Figure 13d and e show that the larger emissions of biogenic VOCs and soil NO x driven by the warmer weather (Figs. 11-12) each could have enhanced surface ozone by 1-2 ppbv over central eastern China in 2017. We find that the soil ozone enhancements are more important in driving the surface ozone changes in northern Asia (e.g., Mongolia and Inner Mongolia), where ozone chemical regimes are highly NO x sensitive (Fig. S7), while both soil and biogenic emissions enhance surface ozone in central eastern China where ozone production is sensitive to both NO x and VOCs during May-August. In addition, the abnormal northwestern winds over Mongolia in 2017 (Fig. 12a) promote the transport of large meteorology-driven soil ozone enhancement to the NCP regions. Although lightning NO x emissions are smaller in 2017 than 2016, stronger subsidence in the lower troposphere (Fig. 12d) in 2017 allows for stronger transport of lightningproduced ozone from the middle-upper troposphere to the surface, contributing slightly higher surface ozone in the NCP and FWP regions by about 1 ppbv (Fig. 13f). This also leads to slightly higher stratospheric contributions to surface ozone by 0.6 ppbv in 2017.
Furthermore, changes in meteorology can alter tropospheric ozone by modulating the chemical kinetics as discussed in Sect. 3.4. The warmer and drier weather conditions in May-August 2017 are expected to accelerate ozone production rates and also favor the thermal decomposition of peroxyacetyl nitrate (PAN, a NO x reservoir species whose decomposition is strongly temperature dependent), providing additional NO x to produce ozone (Steiner et al., 2010;Doherty et al., 2013). This is evident from the lower background PAN concentrations (model results in noCH only reflect changes in meteorology) over eastern China simulated in 2017 (Fig. 14a), which follow the spatial pattern of higher temperature (Fig. 12a). Additional influences can come from changes in ozone dry deposition. The model parametrization of dry deposition includes suppression of the stomatal uptake of ozone and thus the suppression of ozone dry deposition due to the closure of stomata to protect plants from desiccation in warm air and soil temperature (typically > 293 K) (Wesely, 1989). We show in Fig. 14b that the ozone dry deposition velocities decrease by about 2-10 % in May-August 2017 compared to 2016, so more ozone can remain in the surface layer.
The changes in Chinese anthropogenic ozone as shown in Fig. 13c result from the combined effects of changes in domestic anthropogenic emissions and meteorological conditions. While the overall decreases in CH anthropogenic ozone mainly reflect decreasing anthropogenic NO emissions in 2017, the expected ozone reductions due to domestic anthropogenic emission controls may have been largely offset by the higher production rates and lower deposition velocities induced by meteorology as discussed above. Figure 13c also accounts for the impact of decreasing anthro- pogenic PM 2.5 on ozone by affecting heterogeneous chemistry on the aerosol surface (i.e., reactive uptake of HO 2 , N 2 O 5 , NO 2 , and NO 3 ) and photolysis rates (Lou et al., 2014;Li et al., 2019). A recent study by Li et al. (2019) shows that reducing PM 2.5 levels by 30-50 µg m −3 (observed changes in 2017 compared to 2013), as well as the associated aerosol optical depths and surface areas in GEOS-Chem, leads to about 3-5 ppbv of ozone enhancement in eastern China. We find in Fig. S8 that modeled May-August mean anthropogenic PM 2.5 levels decrease by about 4-10 µg m −3 in central eastern China in 2017 from 2016 levels, which may then lead to about 0.4-1 ppbv of ozone enhancement. The above results reflect strong interannual variability of ozone induced by changes in meteorological conditions and also indicate that more stringent emissions controls on ozone precursors are required in the summertime to attain the air quality standard under unfavorable weather conditions. In this study, we have estimated the sources contributing to surface MDA8 ozone over China in 2016-2017 using the state-of-art GEOS-Chem CTM at the 0.25 • × 0.3125 • horizontal resolution with the latest Chinese anthropogenic emission inventory. These 2 years show the highest surface ozone concentrations in record over many major Chinese cities, and ozone increases are observed in most Chinese cities in 2017 compared to 2016. We show that with the latest anthropogenic emission inventory and high horizontal resolution, the model captures the spatial variability of surface ozone in Chinese cities well, particularly in the peak ozone season (May-August, r = 0.81-0.82), with a small positive mean bias (6.6 ppbv), although it has some limitations in reproducing the large daily variability in individual cities.
We quantify the contributions from background, natural, and anthropogenic sources and examine their relative roles in driving MDA8 ozone changes. We find that at the national scale the model identifies large contributions from background sources to surface MDA8 ozone concentrations, i.e., 90 % (49.4 ppbv) in March-April and 80 % (44.5 ppbv) in May-August averaged over China, including large proportions from natural sources (80 % in March-April and 72 % in May-August). Further diagnosing the major natural sources affecting background ozone, we find that biogenic VOC emissions alone enhance surface MDA8 ozone by more than 15 ppbv over central eastern China in July-August, mainly driven by high temperature. Lightning and stratospheric contributions are large (together over 20 ppbv of ozone on average during springtime) in western China due to the high elevation and can also enhance ozone in northern China by over 8 ppbv. Their contributions to ground-level ozone are linked to large-scale vertical transport and are less important in summer than other seasons. Soil emissions can lead to surface ozone increases by more than 5 ppbv in northern China during May-August. Biomass burning emissions are an important ozone contributor in southwestern China (up to 10 ppbv in Yunnan Province) during March-April but have low ozone enhancements (less than 2 ppbv) in other regions. In general, we find higher background ozone concentrations in western and northern China during springtime (March-April) due to the combined effect of larger natural contributions from the stratosphere, lightning emissions, and higher ozone imported from foreign anthropogenic sources.
MDA8 ozone contributions from domestic anthropogenic sources are about 19.5 % (10.8 ppbv) in May-August and less than 15 % in other months averaged over China. For the NCP, YRD, and FWP city clusters, domestic anthropogenic contributions are, respectively, 23.0 %, 34.1 %, and 26.5 % averaged for March-October. The contributions increase to 30.7 % in the NCP, 43.1 % in the YRD, and 32.2 % in the FWP in May-August, reflecting the increasing importance of domestic anthropogenic contributions to urban surface ozone in peak ozone seasons. For the extremely high ozone episodes with simulated MDA8 ozone > 100 ppbv, domestic anthropogenic contributions reach 42.0 ± 5.3 ppbv (41 %) in NCP, 39.1 ± 4.8 ppbv (38 %) in FWP, 61.4 ± 7.7 ppbv (58 %) in YRD, and 73.2 ± 7.0 ppbv (69 %) in PRD. Our model would predict no days with MDA8 ozone > 80 ppbv in the absence of domestic anthropogenic emissions in all 169 major Chinese cities. We also find that day-to-day ozone evolution in Beijing, Shanghai, and Xi'an shows significant positive correlations with the daily maximum temperature, and such significant correlations are widespread over central eastern China. The high ozone episodes in these cities are typically associated with high temperature conditions, which cause intensive biogenic VOC emissions and accelerate ozone production from both background and anthropogenic sources.
We further examine the relative contributions from different sources to surface MDA8 ozone changes in May-August 2017 compared to 2016. The model successfully reproduces the MDA8 ozone increases over central eastern China by 2-5 ppbv, although it underestimates the absolute change compared to the observation. The model attributes the ozone increases to background ozone changes. Higher background ozone concentrations are found in 2017 driven by hotter and drier weather conditions, which accelerate ozone chemical production and increase ozone enhancements from both natural soil NO x and BVOC emissions by 1-2 ppbv. Additional impacts of the 2017 meteorology include increasing air subsidence (therefore enhancing lightning and stratospheric ozone contributions), increasing thermal decomposition of PAN, and decreasing ozone dry deposition velocity. The 2017-2016 changes in anthropogenic emissions would have decreased surface MDA8 ozone by about 1 ppbv over central eastern China and 2 ppbv in the PRD regions. Our results indicate that more stringent emission control measures on both NO x and NMVOCs are needed to reduce present-day surface ozone pollution over China, particularly on days and in years with unfavorable meteorological conditions. The focus of this study is to understand the sources contributing to surface ozone over China at present-day (2016 and 2017) emission levels, and it does not address the nonlinear changes in source attribution and the effect on surface ozone when anthropogenic emissions levels change in the future. To further estimate the nonlinear response of ozone to changes in sources, we conduct an additional sensitivity simulation with Chinese domestic anthropogenic emissions reduced by 20 % in the model for July 2017. As shown in Fig. S9, compared to the zero-out method, the 20 % perturbation method estimates much lower domestic anthropogenic contributions (6.8 ppbv vs. 11.8 ppbv averaged over China in July 2017) with similar spatial distributions, which is consistent with Ni et al. (2018). Such strong nonlinear responses of ozone to precursor emissions in China reflect the fact that more stringent anthropogenic emission control measures are required to mitigate ozone pollution. Although the 2016-2017 period is warmer than 2013-2015 over most regions in China (Fig. S10), we do not see continuous increases in temperature over central eastern China from 2013 to 2017. This implies that temperature changes alone are not enough to explain the continuously increasing ozone trends in Chinese cities over 2013-2017 . Other meteorological factors and modes, such as the summer monsoon Lu et al., 2018b) and the West Pacific Subtropical High , may also modulate the interannual variability of ozone over China. More modeling studies are still needed to better understand ozone changes over China in a longer time period and to evaluate emission control strategies.
Data availability. The datasets, including measurements and model simulations, used in this study can be accessed by contacting the corresponding authors (Lin Zhang, zhanglg@pku.edu.cn; Xiao Lu, luxiao_atchem@pku.edu.cn).
Author contributions. LZ and XL designed the study. XL performed model simulations and conducted data analysis with the assistance of YFC, MZ, BZ, KL, and YML. JTL and TMF assisted in the interpretation of the results. BZ and QZ provided the Chinese anthropogenic emissions inventories. XL and LZ wrote the paper. All authors contributed to the discussion and improvement of the paper.
Competing interests. The authors declare that they have no conflict of interest.