1 Introduction
O is a strong oxidant detrimental to human health (Lu et al., 2020; Liu et al., 2018a) and vegetation growth (Monks et al., 2015; Wang et al., 2017a). Furthermore, it is a crucial active compound influencing the earth's radiative balance, with an effective radiative forcing of up to 0.47 W m in 2019 (IPCC, 2021). Tropospheric O has garnered significant attention over the past few decades due to its crucial role in air quality and climate change (Duan et al., 2017; Li et al., 2019; Ashmore and Bell, 1991; Lu et al., 2018).
With the rapid development in China, emissions of O precursors have been on the rise, leading to an annual increase in O concentrations since the beginning of the 20th century (Liu and Wang, 2020a; Ma et al., 2016). Surface O pollution has become a severe air quality concern in China (Verstraeten et al., 2016; Xu et al., 2018), particularly in major urban areas such as the North China Plain (NCP), Fenwei Plain (FWP), Yangtze River Delta (YRD), Pearl River Delta (PRD), and the Sichuan Basin (SCB) (Y. Wang et al., 2020; T. Wang et al., 2017a; Yin and Ma, 2020; Shen et al., 2019; Zhao et al., 2018; X. Wang et al., 2009). Although the Chinese government initialized the Clean Air Action Plan in 2013 to control air pollution, the concentration of O precursors and PM has significantly decreased (Zhai et al., 2019). However, surface O concentrations continue to increase in major urban areas.
Recent studies have suggested that regional meteorological conditions influence surface O levels through various pathways (Jacob and Winner, 2009; Shen et al., 2016; Lin et al., 2008). Modeling studies have shown that O levels are sensitive to temperature, humidity, wind speed, mixing height, and other meteorological conditions (Pfister et al., 2014; Sanchez-Ccoyllo et al., 2006). For instance, temperature affects the chemical formation rate of O (Lee et al., 2014), while precipitation reduces surface O concentrations through wet removal (Fang et al., 2011). Additionally, the elevated planetary boundary layer (PBL) height enhances upward movement, resulting in lower surface O concentrations (Haman et al., 2014). Therefore, long-term modeling of surface O levels must consider changes in meteorological conditions.
CO is the primary anthropogenic radiative force of the climate system (Gauss et al., 2003; Schimel et al., 2015). CO can impact regional air temperature and precipitation, leading to changes in surface O concentrations (Lu et al., 2013; Yang et al., 2014).
On the other hand, biogenic volatile organic compounds (BVOCs) are significant O precursors, and isoprene is the primary species among BVOCs that vegetation emits (Zheng et al., 2009; Fiore et al., 2011). In most of China, O is volatile organic compound (VOC) limited in the summer, especially in industrial cities (Li et al., 2018; Wu et al., 2018). Thus, it plays a significant role in modulating O levels and positively correlates with O concentrations in major urban areas of China. It is known that CO can enhance vegetation's photosynthesis (Sun et al., 2013; Heald et al., 2009; Tai et al., 2013; Monson and Fall, 1989), which may directly increase isoprene emission (Rapparini et al., 2004). Based on the observation, Rosenstiel et al. (2003) found that the isoprene emissions of plants grew by about 21 % and 41 % when CO reached 800 and 1200 ppm, respectively. However, Wilkinson et al. (2009) indicated that different vegetation types show varying responses in isoprene emission when CO increases. Isoprene emission was decreased by 30 %–40 % in Populus tremuloides Michx but increased by about 100 % in Quercus rubra when CO concentrations were grown (Sharkey et al., 1991). High concentrations of CO may inhibit the emission of isoprene by reducing the activity of BVOCs synthetase or decreasing the synthesis of adenosine triphosphate (Possell et al., 2005). Guenther et al. (1991) also indicated that isoprene emissions were significantly reduced when CO was increased from 100 to 600 mol mol. In summary, the impact of elevated CO on isoprene emission may be positive or negative, mainly related to the relative size of the inhibition caused by elevated CO and promotion by enhanced photosynthesis.
Numerous studies have concluded that anthropogenic emissions are the primary drivers of surface O increases in different regions or years in China. Meanwhile, the effects of meteorological parameters can be substantial (Wang et al., 2019c; Lu et al., 2019; Dang et al., 2021; Liu and Wang, 2020a). For instance, Li et al. (2020) indicated that anthropogenic emissions were the primary cause of surface O increase in China from 2013 to 2019. Liu and Wang (2020a) suggested that anthropogenic emissions play a dominant role in the O variety in China, but the effects of meteorological conditions could be more significant in some regions. Han et al. (2020) analyzed the O changes in summer and suggested that meteorology can explain about 43 % of that in eastern China.
Previous studies have mainly focused on the impact of anthropogenic emissions and meteorological factors on the rise of O levels, with limited attention given to the role of CO variations. However, due to the rapid socioeconomic growth in China and the subsequent surge in energy consumption, CO emissions and concentrations have also increased significantly, particularly in the eastern coastal region (Lv et al., 2020; Ren et al., 2014). Furthermore, given the significant impact of CO on O, it is crucial to evaluate the influence of changes in CO concentration on the maximum daily 8 h average (MDA8) O concentrations at the surface. Thus, a comprehensive assessment of the impact of anthropogenic emissions, meteorological factors, and CO on surface O is imperative.
Here, we employed an advanced regional climate–chemistry–ecology model to assess the impact of anthropogenic emissions, meteorological factors, and carbon dioxide variations during the summer monsoon period (May, June, July, and August) on surface O levels. Our findings can facilitate the development of a comprehensive O improvement strategy. Section 2 describes the methods and data, and the results and discussion are given in Sect. 3; finally, the conclusions are shown in Sect. 4.
2 Methods and data
2.1 Measurement data
We compared the simulated regional meteorological factors with the European Centre for Medium-Range Weather Forecasts Interim reanalysis data (ERA-Interim) at 37 vertical levels, which included temperature, relative humidity, and wind speed (Balsamo et al., 2015; Hoffmann et al., 2019). The observed surface O was taken from the China National Environmental Monitoring Center (CNEMC), which had more than 1400 environmental monitoring stations in 2018 (Wang et al., 2018; Kong et al., 2021; Zheng et al., 2014). The World Data Centre for Greenhouse Gases (WDCGG) data (Liu et al., 2009; Li et al., 2017) were applied to evaluate the simulated surface CO concentrations. The monitoring stations of O and CO are shown in Fig. 1.
Figure 1
Model domains for the RegCM-Chem-YIBs model. The regions with black boundaries are the North China Plain (34–41 N, 113–119 E), the Yangtze River Delta (30–33 N, 119–122 E), the Pearl River Delta (21.5–24 N, 112–115.5 E), the Sichuan Basin (28.5–31.5 N, 103.5–107 E), and the Fenwei Plain (33.5–39 N, 106–113 E) regions.
[Figure omitted. See PDF]
2.2 Model descriptionThe RegCM-Chem-YIBs is a regional climate–chemistry–ecology model developed from the RegCM model. RegCM is a regional climate model initially developed by the International Centre for Theoretical Physics (ICTP) (Giorgi et al., 2012). Shalaby et al. (2012) integrated the Chem chemistry model into the RegCM model and incorporated the condensed version of the Carbon Bond Mechanism Z (CBM-Z) to enhance the model's capabilities. To further improve the model's performance, Yin et al. (2015) added a volatility basis set (VBS) scheme to simulate secondary organic aerosols (SOAs). Xie et al. (2020) further modified the model by incorporating CO as a tracer, which is subject to regulation by sources, sinks, and atmospheric transport processes. The model represents the four sources and sinks of CO as surface fluxes, including emissions from fossil fuels and biomass burning, air–sea CO exchange, and terrestrial biosphere CO fluxes. Additionally, the model incorporated the Yale Interactive Terrestrial Biosphere (YIBs), a land carbon cycle model that enables the simulation of ecological processes, including carbon assimilation, allocation, and autotrophic and heterotrophic respiration (Yue and Unger, 2015).
The ecological model (YIBs) was fully coupled into the regional climate–chemical model (RegCM-Chem) to reproduce the interactions between atmospheric composition and the ecosystem in the actual atmosphere (Xie et al., 2019). The meteorological factors and air components simulated by RegCM-Chem were input into the YIBs model to simulate vegetation physiological processes and calculate land surface parameters such as carbon dioxide flux, BVOCs emissions, and stomatal conductance. Conversely, the simulations of the YIBs model were fed back into the RegCM-Chem model to adjust the air qualities, temperature, humidity, circulation, and other meteorological fields. The RegCM-Chem-YIBs has been extensively applied to study surface O, PM, CO, the summer monsoon, and the interactions between air quality and the ecosystem (Zhuang et al., 2018; Pu et al., 2017; Xu et al., 2022; Xie et al., 2018; Ma et al., 2023).
The RegCM model offers a variety of physical and chemical parameterization options. Here, the climatological chemical boundary conditions were driven by the Model for Ozone and Related Chemical Tracers (MOZART). The gas-phase chemistry employed the CBM-Z scheme (Zaveri and Peters, 1999). For the boundary layer scheme, the Holtslag PBL approach was utilized (Khayatianyazdi et al., 2021). The Grell cumulus convection scheme was employed to simulate convective processes (Grell, 1993). The CCM3 radiation scheme and CLM3.5 land surface module were used in the model (Collins et al., 2006; Giorgi and Mearns, 1999; Decker and Zeng, 2009).
Figure 2
Changes in the surface MDA8 O concentrations (units: ppb) during the summer monsoon period from 2009 (a), 2010 (b), 2011 (c), 2012 (d), 2013 (e), 2014 (f), 2015 (g), 2016 (h), 2017 (i), and 2018 (j) relative to 2008.
[Figure omitted. See PDF]
Table 1The numerical experiment in this study.
Experiment | Time | Meteorological fields | CO emissions | Anthropogenic emissions |
---|---|---|---|---|
Base | 2008–2018 | Varying | Varying | Varying |
SIM | 2009–2018 | 2008 | Varying | Varying |
SIM | Varying | 2008 | Varying |
2.3.1 Radiation
In the previous version of the RegCM-Chem-YIBs model, radiative calculations only accounted for changes in the spatiotemporal distribution of particulate matter. To simplify the radiation calculations, the atmospheric CO and O concentrations were assumed to be constant throughout the year. However, atmospheric CO and O are subject to modulation by various sources, sinks, physical processes, and chemical processes (Ballantyne et al., 2012; Wang et al., 2019a). Additionally, rapid urbanization in China has led to an annual increase in CO and O concentrations (Guan et al., 2021; Wei et al., 2022), with elevated concentrations and growth rates primarily distributed in the eastern regions where urbanization is most prominent (Shi et al., 2016; Wang et al., 2017b). To more accurately simulate the atmospheric radiation balance and East Asian monsoon climate, it is necessary to incorporate spatiotemporal variations of CO and O concentrations into the radiation module. Therefore, we included the varying CO and O concentrations simulated by the model in the radiation module to calculate the corresponding radiative forcing.
2.3.2 Photolysis
The photolysis process was simulated using the Tropospheric Ultraviolet and Visible (TUV) model, which is commonly used to compute photolysis rates in various models (Tie et al., 2003; Shetter et al., 2002; Borg et al., 2011). The TUV model employs input parameters such as zenith angle, altitude, ozone column, SO column, NO column, aerosol optical depth (AOD), single scattering albedo (SSA), and albedo, among others, to calculate photolysis rates (Singh and Singh, 2004). However, in the TUV module of the RegCM-Chem-YIBs model, AOD and SSA were held constant. This is problematic as accurate aerosol optical parameters, such as AOD and SSA, play a crucial role in the photolysis of O (Lefer et al., 2003). To address this issue, we incorporated temporally and spatially varying AOD and SSA simulated by the RegCM-Chem-YIBs model into the photolysis rate calculations in the TUV module. This enabled us to accurately incorporate the extinction effect of the varying particles into the photolysis reaction, leading to more realistic simulations of air components and regional meteorology.
2.4 Emissions and Experiment settings
Anthropogenic emissions from 2008 to 2018 were obtained from the Multi-resolution Emission Inventory for China (MEIC), which has been compiled and maintained by Tsinghua University since 2010 (Zheng et al., 2018; Wang et al., 2014). CO emissions and boundary conditions were derived from the NOAA CarbonTracker CT2019 dataset (Jacobs et al., 2021). The initial meteorological boundary data, such as temperature, relative humidity, and wind, are derived from the ERA-Interim reanalysis dataset with a horizontal resolution of 0.125, a temporal resolution of 6 h, and 37 vertical levels (Liu et al., 2018b). The weekly mean sea surface temperature dataset was obtained from the National Ocean and Atmosphere Administration (NOAA) (Reynolds et al., 2002).
The simulation domain was illustrated in Fig. 1, with the target region centered at 36 N and 107 E and with a grid resolution of 60 km by 60 km. The model used 18 vertical levels, ranging from the surface to 50 hPa.
In the Base experiment, we incorporated interannual variations in anthropogenic emissions, meteorological fields, and CO emissions. Meteorological conditions (CO emissions) were kept constant at 2008 levels over 10 years, referred to as the SIM (SIM) experiment.
The changes in O concentrations relative to 2008 between 2009 and 2018 were determined by comparing simulations of different years with 2008 in the Base experiment (Eq. 1). The impact of changed meteorological conditions on O concentrations relative to 2008 was assessed by comparing results between SIM and the Base experiment in the same year (Eq. 2). The contribution of changed CO emissions was similarly estimated (Eq. 3). Finally, the influence of anthropogenic emissions was calculated by excluding the impact of meteorological factors and CO from the changes in O concentrations (Eq. 4). Table 1 shows the results of the numerical experiments. O represents the changes in O concentrations in the year relative to 2008. Base represents the O concentrations in the Base experiment in the year . represents the changes in O concentrations in the year due to meteorological factors variations. SIM represents the O concentrations in the SIM experiment in the year . C represents the changes in O concentrations in the year due to CO variations. SIM represents the O concentrations in the SIM experiment in the year . represents the changes in O concentrations in the year due to anthropogenic emissions variations.
In this work, both meteorological and CO boundary conditions were kept consistent in base and sensitivity studies. We did not consider the impact of boundary conditions on O due to the following reasons. First, in general, the regional model was coupled with the global model to get a more realistic influence from the boundary. However, for long-term climate–chemistry modeling, such a coupling means a large computing resource. Second, the boundary conditions were derived from global models (Liu et al., 2017; Ban et al., 2014) and have to be prescribed in numerical experiments. Finally, fixed boundary conditions were widely used in some O studies in China (Liu and Wang, 2020a, b; Wang et al., 2019b). Moreover, regional emissions are the primary source of surface O in China, with contributions accounting for 80 % from May to August (Lu et al., 2019). Therefore, the impact of fixed boundary conditions can be ignored in the current stage.
Figure 3
The response of the surface MDA8 O mixing ratios (units: ppb) to variations in meteorological conditions during the summer monsoon period in 2009 (a), 2010 (b), 2011 (c), 2012 (d), 2013 (e), 2014 (f), 2015 (g), 2016 (h), 2017 (i), and 2018 (j) relative to 2008.
[Figure omitted. See PDF]
Figure 4
The MDA8 O (a–c, units: ppb), precipitation (d–f, units: mm d), clouds (g–i, units: %), shortwave flux (j–l, units: W m), wind speed (m–o, units: m s), temperature (p–r, units: K), and planetary boundary layer height (s–u, units: m) during the summer monsoon period in 2008 from the base simulations (the left column) and their responses due to variations in meteorological conditions in PreG (2009–2013, the central column) and PostG (2014–2018, the right column) relative to 2008.
[Figure omitted. See PDF]
3 Results and discussion3.1 Model evaluation
The ability of RegCM to reproduce East Asian climate and air quality has been widely evaluated in recent years. Previous studies have demonstrated that RegCM was capable of the essential characteristics and interannual variations of air components and meteorological fields in East Asia (Xu et al., 2023; Ma et al., 2023; Zhuang et al., 2018). Given that the monitoring of near-surface O levels by CNEMC was initiated only in late 2013, the monitoring sites in 2013 and 2014 were limited, and the monitoring period was disjointed. As a result, in this study, we compared the simulated meteorological fields, O, and CO levels with observations only from 2015 to 2018.
Figures S1–S4 demonstrated that the RegCM-Chem-YIBs model effectively captured the spatial distribution and magnitude of temperature, humidity, and wind over East Asia at 500, 850, and 1000 hPa between 2015 and 2018. However, due to the complex terrain's influence on the lower atmosphere, most models show better results at higher levels (Zhuang et al., 2018; Anwar et al., 2019; Xie et al., 2019). Thus, the simulations at 500 hPa were more consistent with the reanalysis data. At 1000 hPa, the simulated wind speed was slightly higher than the reanalysis data in eastern China. This difference may be due to common deficiencies in meteorological models, such as insufficient horizontal resolution, initial and boundary conditions, and physical parameterizations (Cassola and Burlando, 2012; Accadia et al., 2007), particularly in areas with low wind speeds (Carvalho et al., 2012).
Figures S5 and S6 demonstrated that the model accurately reproduced the observed increase in surface CO and O from 2015 to 2018, with high correlation coefficients ranging from 0.39 to 0.74 (Table 2). The model effectively captured the high concentrations of O in major urban areas such as the NCP, the YRD, the PRD, the SCB, and the FWP, while also successfully reproducing the gradient in CO concentrations between eastern and western China. However, the model slightly underpredicted MDA8 O concentrations (4.02 to 3.21 ppb) and overestimated CO levels (3.32–7.07 ppm). These discrepancies are mainly attributed to uncertainties in the emissions inventory (Hong et al., 2017). Overall, the simulated meteorological factors and surface CO and O concentrations were deemed acceptable.
Table 2
Evaluations of the surface CO (units: ppm) and MDA8 O (units: ppb) during the summer monsoon period in East Asia.
Species | Year | OBS | SIM | MB | RMSE | |
---|---|---|---|---|---|---|
CO (ppm) | 2015 | 402.82 | 406.98 | 4.16 | 9.37 | 0.44 |
2016 | 407.12 | 410.44 | 3.32 | 8.22 | 0.69 | |
2017 | 408.35 | 413.62 | 5.27 | 11 | 0.39 | |
2018 | 409.61 | 416.68 | 7.07 | 11.32 | 0.41 | |
MDA8 O (ppb) | 2015 | 48.77 | 44.75 | 4.02 | 29.39 | 0.57 |
2016 | 50.16 | 46.95 | 3.21 | 27.56 | 0.60 | |
2017 | 55.43 | 51.87 | 3.56 | 21.55 | 0.74 | |
2018 | 55.53 | 52.08 | 3.42 | 24.78 | 0.73 |
OBS: observation; SIM: simulation; MB: bias; RMSE: root mean square error; : correlation coefficient. MDA8 O: the maximum daily 8 h average O.
Table 3The changes of MDA8 O over the North China Plain, Fenwei Plain, Yangtze River Delta, Pearl River Delta, and Sichuan Basin during the summer monsoon period from 2009 to 2018 relative to 2008.
Regions | 2009 | 2010 | 2011 | 2012 | 2013 | 2014 | 2015 | 2016 | 2017 | 2018 | PreG | PostG |
---|---|---|---|---|---|---|---|---|---|---|---|---|
NCP | 0.14 | 2.85 | 4.53 | 6.13 | 2.7 | 4.78 | 10.1 | 14.25 | 30.92 | 32.03 | 3.27 | 18.42 |
FWP | 3.23 | 1.78 | 5.01 | 6.78 | 1.37 | 7.9 | 10.5 | 6.24 | 13.71 | 16.17 | 3.63 | 10.90 |
YRD | 8.33 | 1.47 | 1.46 | 0.5 | 3.12 | 6.04 | 3.46 | 7.09 | 17.64 | 16.12 | 2.98 | 10.07 |
PRD | 5.76 | 0.26 | 2.56 | 5.13 | 0.4 | 3.82 | 1.46 | 3.16 | 9.45 | 6.79 | 2.56 | 4.94 |
SCB | 4.92 | 1.03 | 3.46 | 5 | 3.94 | 8.54 | 9.27 | 9.78 | 13.67 | 14.8 | 3.67 | 11.21 |
Figure S7 illustrates the mean seasonal MDA8 O concentrations in East Asia during the summer monsoon period from 2008 to 2018. High O concentrations appeared in eastern China, which can be attributed to increased emissions, high temperatures, high humidities, and intense radiation in the region (Gao et al., 2020; Mousavinezhad et al., 2021; Wei et al., 2022). Surface O increased annually in most of China between 2008 and 2018, with megacity clusters experiencing a more significant increase.
We conducted a regional analysis of surface O increases in five target regions: the NCP, the YRD, the PRD, the SCB, and the FWP. In 2018, the surface MDA8 O concentrations averaged 74 ppb in the NCP region, while the other areas had lower concentrations (ranging from 42 to 67 ppb in the FWP, YRD, PRD, and SCB). The lower surface O levels in the SCB and FWP regions were attributed to lower anthropogenic emissions. The YRD and PRD regions were more affected by meteorological factors, with the East Asian summer monsoon bringing in cleaner air and precipitation from the sea, leading to lower air pollution concentrations (He et al., 2012). The spatial distribution and increasing trend of surface MDA8 O concentrations in China were consistent with previous studies (Li et al., 2020; Shen et al., 2022).
Figure 2 and Table 3 illustrate the changes in surface MDA8 O concentrations from 2009 to 2018 relative to 2008. The surface MDA8 O concentrations in China increased drastically over the past decade, particularly in 2017 and 2018 (6.79–32.03 ppb). We divided the period from 2009 to 2018 into two phases based on the Clean Air Action Plan implemented in 2013: the pre-governance period (PreG, 2009–2013) and the post-governance period (PostG, 2014–2018). The surface MDA8 O concentration increased significantly in NCP (18.42 ppb), followed by SCB (11.21 ppb), FWP (10.9 ppb), and the YRD (10.07 ppb), while it increased slightly in PRD (4.94 ppb) in PostG relative to 2008. Our results were consistent with previous studies by Lu et al. (2020), Ma et al. (2016), and Mousavinezhad et al. (2021).
Table 4
Response of the MDA8 O mixing ratios (units: ppb), precipitation (units: mm d), cloud fraction (units: %), shortwave flux (units: W m), wind speed (units: m s), temperature (units: K), and planetary boundary layer height (units: m) to the changes in meteorological conditions over the North China Plain, Fenwei Plain, Yangtze River Delta, Pearl River Delta, and Sichuan Basin during the summer monsoon period in PreG (2009–2013) and PostG (2014–2018) relative to 2008.
Regions | Period | MDA8 O | Precipitation | Clouds | SWF | Wind speed | Temperature | PBL |
---|---|---|---|---|---|---|---|---|
(ppb) | (mm d) | (%) | (W m) | (m s) | (K) | (m) | ||
NCP | PreG | 0.88 | 0.58 | 1.33 | 3.04 | 0.17 | 0.32 | 46.8 |
PostG | 0.04 | 0.6 | 0.93 | 3.06 | 0.26 | 0.6 | 14.5 | |
FWP | PreG | 1.41 | 1.68 | 2.86 | 10.63 | 0.06 | 0.1 | 108.5 |
PostG | 0.09 | 0.81 | 0.94 | 0.81 | 0.05 | 0.46 | 15.3 | |
YRD | PreG | 1.03 | 1.02 | 1.07 | 1.6 | 0.18 | 0.29 | 33.9 |
PostG | 0.96 | 0.48 | 1.18 | 4.85 | 0.08 | 0.45 | 21.9 | |
PRD | PreG | 0.23 | 2.39 | 1.93 | 2.24 | 0.02 | 0.36 | 29.6 |
PostG | 1.08 | 3.24 | 3.98 | 5.37 | 0.18 | 1.00 | 52.2 | |
SCB | PreG | 0.41 | 1.81 | 0.59 | 8.8 | 0.13 | 0.58 | 136.5 |
PostG | 0.71 | 0.37 | 2.23 | 3.2 | 0.03 | 0.14 | 76 |
Overall, the meteorological variations from 2008 to 2018 were unfavorable for the O increase during the EASM (East Asia summer monsoon) period, as illustrated in Fig. 3.
Based on Fig. 3 and Table 4, it is evident that meteorological conditions had a significant impact on surface MDA8 O in the NCP and FWP regions during the PostG period (0.09 to 0.04 ppb) compared to the PreG period (1.41 to 0.88 ppb). In the SCB region, the impact of meteorological fields was relatively weak (0.41–0.71 ppb), attributed to the basin topography and stable atmospheric conditions. However, in the eastern and southeastern coastal areas of China, due to the significant influence of the EASM, the impact of meteorological conditions may be more critical than that of anthropogenic emissions. For instance, in the YRD and PRD regions, meteorological conditions significantly changed O levels (1.29–1.3 ppb) compared to anthropogenic emissions (0.81–0.87 ppb) in 2013, indicating the significant influence of meteorological conditions on surface O.
Figure 5
Simulated responses of surface MDA8 O mixing ratios (units: ppb) to the variations in CO emissions during the summer monsoon period in 2009 (a), 2010 (b), 2011 (c), 2012 (d), 2013 (e), 2014 (f), 2015 (g), 2016 (h), 2017 (i), and 2018 (j) relative to 2008.
[Figure omitted. See PDF]
Figure 6
The simulated averaged MDA8 O (a–c, units: ppb), CO (d–f, units: ppm), precipitation (g–i, units: mm d), and isoprene mixing ratios (j–l, units: g m) in 2008 from the base simulations (the left column) and their changes due to variations in CO emissions in PreG (2009–2013, the central column) and PostG (2014–2018, the right column) relative to 2008.
[Figure omitted. See PDF]
Our findings are consistent with previous studies. Liu and Wang (2020a) reported a decrease in O in Shanghai from 2013 to 2017 due to changes in meteorological conditions. Chen et al. (2019) and Liu and Wang (2020a) also suggested that changed meteorological conditions had a negative impact on O formation in the NCP and FWP regions and that the influence of meteorology on surface-level O decreased in PostG. In addition, Cheng et al. (2019) found that the effects of meteorological conditions on long-term O variations were less than 3 %, which is similar to our study.
As we know, the formation of surface O is promoted by rising temperatures (Steiner et al., 2010). However, increased surface temperatures can also intensify turbulence within the planetary boundary layer (PBL), increasing PBL height (Guo et al., 2016). This increase in PBL height, coupled with the enhanced upward motion, can transport near-surface pollutants to the upper atmosphere, reducing their concentration in the lower atmosphere (Gao et al., 2016). Additionally, the upward motion can also facilitate cloud formation and precipitation, resulting in a reduction of near-surface atmospheric pollutants via precipitation washout (Yoo et al., 2014).
We have improved the accuracy of O photodissociation rate calculations by including varying AOD and SSA in the TUV module, as described in Sect. 2.3.2. As a result, the increase in cloud cover reduced the shortwave radiation flux and photochemical formation rates of near-surface O, leading to decreased formation. Thus, the increase in near-surface temperature is often accompanied by an elevation in PBL height, enhanced cloud cover, precipitation, and reduced shortwave radiation. Moreover, higher wind speeds can enhance the dispersion of O (Gorai et al., 2015).
The variations of MDA8 O, precipitation, clouds, shortwave flux (SWF), wind speed, temperature, and PBL height are presented individually in Fig. 4. The increase in SWF can accelerate O formation through photochemistry (Jiang et al., 2012; Lelieveld and Crutzen, 1990). Therefore, the increased cloud fraction reduced surface O by decreasing shortwave radiation, especially in NCP, FWP, YRD, and SCB in the PreG period (10.63 to 1.6 W m). Furthermore, the enhanced precipitation in these regions (0.37–1.81 mm d) reduced surface O levels significantly. The significant increase in wind speed (0.17–0.26 m s) also contributed to the reduction of surface O in the NCP region (Table 4).
Another crucial factor is the elevated surface temperature (0–5 K), which intensified upward motion and raised the planetary boundary layer (PBL) height (0–500 m) across much of East Asia. Consequently, the increased temperature and PBL height could disperse surface-level O, thereby reducing its concentration.
3.4The effect of CO in the 2008–2018 ozone increase
The surface O in southern China, which includes the YRD, PRD, and SCB regions, was characterized by high precipitation, temperatures, and vegetation cover and was significantly impacted by CO (Fig. 5). Figure 6e and f demonstrate a marked rise in CO levels across East Asia, particularly in eastern China, which was attributable to extensive human activity.
CO is a significant driver of climate change and alterations in biogenic emissions. As shown in Fig. 6b and c, the impact of CO on O levels varies across locations, with a positive effect of 0.28–0.46 ppb along the southeastern coast of China but a negative influence of 0.51 to 0.11 ppb in southwest and central China. CO affects O concentration by influencing both precipitation and isoprene emissions. In western and central China, CO primarily affects O concentration through its impact on precipitation (Table 5). Elevated CO concentrations lead to increased precipitation (0.06–0.64 mm d) in the FWP and SCB regions, resulting in a decrease in surface O (up to 0.51 ppb). In eastern and southern coastal China, where vegetation is abundant, CO has a greater impact on isoprene emissions. In the YRD region, decreased isoprene (0.58 to 0.32 g m) and increased precipitation (0.09–0.13 mm d) reduced MDA8 O levels (0.09–0.14 ppb). In PRD, increased isoprene levels (0.31–0.92 g m) and decreased precipitation (1.02 to 0.33 mm d) led to the enhancement of MDA8 O (0.28–0.46 ppb).
In some years, the impact of changed CO can be as significant as or even surpass that of anthropogenic emissions and meteorology (Fig. 10). For example, in 2013, CO caused an increase of 0.95 ppb in MDA8 O in the YRD region, which exceeded that of anthropogenic emissions (0.87 ppb). Similarly, in the PRD region in 2012, the effect of CO, anthropogenic emissions, and meteorology was 1.41, 1.77, and 1.95 ppb, respectively. Even in the NCP in 2010, the impact of CO (0.75 ppb) was comparable to that of anthropogenic emissions (1.5 ppb). In summary, CO has a significant impact on surface O concentrations by influencing radiation and isoprene emissions, with more prominent effects in regions with abundant vegetation.
Table 5Simulated responses of MDA8 O mixing ratios (units: ppb), CO mixing ratios (units: ppm), precipitation (units: mm d), and isoprene mixing ratios to the changes in CO emissions over the North China Plain, Fenwei Plain, Yangtze River Delta, Pearl River Delta, and Sichuan Basin in PreG (2009–2013) and PostG (2014–2018) relative to 2008.
Regions | Period | MDA8 O | CO | Precipitation | Isoprene |
---|---|---|---|---|---|
(ppb) | (ppm) | (mm d) | (g m) | ||
NCP | PreG | 0.07 | 3.19 | 0.27 | 0.1 |
PostG | 0.05 | 4.24 | 0.13 | 0.26 | |
FWP | PreG | 0.11 | 1.70 | 0.21 | 0.16 |
PostG | 0.51 | 2.05 | 0.06 | 0.33 | |
YRD | PreG | 0.09 | 4.1 | 0.13 | 0.32 |
PostG | 0.14 | 6.2 | 0.09 | 0.58 | |
PRD | PreG | 0.46 | 1.97 | 1.02 | 0.31 |
PostG | 0.28 | 3.20 | 0.33 | 0.92 | |
SCB | PreG | 0.30 | 2.80 | 0.64 | 0.78 |
PostG | 0.30 | 2.78 | 0.21 | 0.69 |
Finally, we calculated the anthropogenic emissions' effect on the 2008–2018 O increase. Figure S8 and Table S1 illustrate that the levels of PM, PM, SO, CO, and OC emissions remained consistently high during the PreG period. However, a linear decrease in emissions was observed after the implementation of the Clean Air Action Plan in 2013. Prior to 2013, the emission of VOCs increased steadily but subsequently stabilized. Similarly, the emission of nitrogen oxides (NO) exhibited an upward trend before 2013, but since then, the emissions have shown a linear decrease, with each subsequent year exhibiting lower levels of NO emissions. In comparison to other species, the emissions of ammonia (NH) remained relatively stable from 2008 to 2018. Our analysis results of the emissions of different species align with those of Zheng et al. (2018), who computed the changes of each species in the MEIC inventory from 2010 to 2017.
Figure 7
Simulated responses of the surface MDA8 O mixing ratios (units: ppb) to variations in anthropogenic emissions in 2009 (a), 2010 (b), 2011 (c), 2012 (d), 2013 (e), 2014 (f), 2015 (g), 2016 (h), 2017 (i), and 2018 (j) relative to 2008.
[Figure omitted. See PDF]
Figure 8
Changes in the simulated surface MDA8 O mixing ratios (units: ppb) from the base simulation (All, a, b); those due to variations in anthropogenic emissions (Emis, c, d), meteorological conditions (Met, e, f), and CO emissions (CO, g, h) in PreG (2009–2013, the left column) and PostG (2014–2018, the right column) relative to 2008.
[Figure omitted. See PDF]
Figure 7 illustrates that anthropogenic emissions have caused a notable increase in surface O levels across most of China, particularly in megacity clusters. The impact of anthropogenic emissions on O concentration ranged from 2.33 to 18.51 ppb in the five target regions.
Figure 8 and Table 6 illustrate that the changes in surface O caused by anthropogenic emissions are similar in magnitude and spatial distribution to the changes in the Base experiment. This suggests that anthropogenic emissions were the dominant factor driving the increase in surface O in China from 2008 to 2018. Notably, a high-impact center of anthropogenic emissions was simulated in northern China, with the NCP region experiencing the most significant increase in surface O (4.08–18.51 ppb), followed by the FWP, YRD, and SCB regions (4.10–11.5 ppb). In the PRD region, anthropogenic emissions led to a slight enhancement of O by 2.33–5.74 ppb.
The role of anthropogenic emissions increased linearly from 2008 to 2018, despite the implementation of the Clean Air Action Plan in 2013 (Table 6). For example, anthropogenic emissions significantly increased surface MDA8 O in the NCP region (4.08 ppb in PreG and 18.51 ppb in PostG). Similarly, FWP experienced increases of 5.15 and 11.5 ppb in the PreG and PostG periods, respectively. In the SCB region, the surface MDA8 O was mainly affected by variations in anthropogenic emissions due to the high levels of anthropogenic emissions in the complex basin topography. In the YRD and PRD regions, anthropogenic emissions resulted in changes to O of 2.33–11.17 ppb.
The reasons for this characteristic are multiple. Before 2013, the continuous increase in VOCs and NO emissions (Fig. S8b, c) facilitated the rise of O levels. Following the implementation of the Clean Air Action Plan in 2013, the emissions of VOCs and NO were regulated. However, with the decrease in PM levels, direct radiation increased, and scattered radiation decreased (Fig. 9), thereby promoting the photochemical formation of O (Bian et al., 2007). In addition, the reduced NO emission weakened the titration effect (Fig. S8b), thus increasing surface O (Li et al., 2022).
Our results are consistent with previous studies by Wang et al. (2019b) and Liu and Wang (2020b), which also showed the dominant and almost linear role of anthropogenic emissions in the increase in O from 2013 to 2015 in four major Chinese cities (Beijing, Shanghai, Guangzhou, and Chengdu).
Figure 9
The variations of the surface direct radiation (a, b, units: W m) and diffuse radiation (c, d, units: W m) in the PreG (2009–2013, a, c) and PostG (2014–2018, b, d) period relative to 2008.
[Figure omitted. See PDF]
Figure 10
Interannual variations of the surface MDA8 O mixing ratios (units: ppb) in the summer monsoon period (ALL) and the responses of variations in anthropogenic emissions (Emis), meteorological conditions (Met), and CO emissions (CO) in the (a) North China Plain, (b) Fenwei Plain, (c) Yangtze River Delta, (d) Pearl River Delta, and (e) Sichuan Basin in 2008–2018 relative to 2008.
[Figure omitted. See PDF]
Table 6Simulated response of the MDA8 O mixing ratios (units: ppb) to the changes in anthropogenic emissions (Emis), meteorological conditions (Met), and CO emissions (CO) over the North China Plain, Fenwei Plain, Yangtze River Delta, Pearl River Delta, and Sichuan Basin in PreG (2009–2013) and PostG (2014–2018) relative to 2008.
Regions | Period | ALL | Emis | Met | CO |
---|---|---|---|---|---|
(ppb) | (ppb) | (ppb) | (ppb) | ||
NCP | PreG | 3.27 | 4.08 | 0.88 | 0.07 |
PostG | 18.42 | 18.51 | 0.04 | 0.05 | |
FWP | PreG | 3.63 | 5.15 | 1.41 | 0.11 |
PostG | 10.9 | 11.5 | 0.09 | 0.51 | |
YRD | PreG | 2.98 | 4.10 | 1.03 | 0.09 |
PostG | 10.07 | 11.17 | 0.96 | 0.14 | |
PRD | PreG | 2.56 | 2.33 | 0.23 | 0.46 |
PostG | 4.94 | 5.74 | 1.08 | 0.28 | |
SCB | PreG | 3.67 | 4.38 | 0.41 | 0.30 |
PostG | 11.21 | 10.80 | 0.71 | 0.30 |
Finally, we presented an attribution diagram depicting the changes in O concentration from 2008 to 2018. The total variation in O concentration can be attributed to the combined effects of meteorological changes, changes in CO concentration, and anthropogenic emissions (Fig. 10).
The primary driver of the O concentration variation from 2008 to 2018 was the changes in anthropogenic emissions, particularly in regions with high emissions, such as the NCP and FWP. Although the Clean Air Action Plan was implemented in 2013, it did not reduce the contribution of anthropogenic emissions to the O increase. Even in the PostG period, with the development of urbanization and industrialization, the impact of changed anthropogenic emissions on O has gradually become more prominent than changed meteorology and CO. The contribution of changed meteorology to O was generally negative in the five regions, with a more significant impact in the YRD and PRD regions. This may be attributed to their proximity to the ocean and susceptibility to the summer monsoon influence. Changes in CO concentration affected O concentration by altering radiation and isoprene emissions, with a more significant impact in the YRD and PRD regions where vegetation was abundant. In some years, it even surpassed the effects of anthropogenic emissions. Therefore, we suggest that the influence of CO concentration changes on O concentration should be considered in regions with high vegetation coverage.
4 Conclusions
In this study. First, we improved the RegCM-Chem-YIBs model regarding the photolysis of O and the radiation effect of CO and O. Second, we assessed the impacts of anthropogenic emissions, meteorological factors, and CO on the surface O increase in China during the EASM season from 2008 to 2018.
In the NCP and FWP regions. The increased surface O was primarily attributed to the changes in anthropogenic emissions (4.08–5.15 ppb in PreG and 11.5–18.51 ppb in PostG). Furthermore, the impact of anthropogenic emissions has increased linearly, despite the Clean Air Action Plan being implemented in 2013. In contrast, the effects of meteorological factors and CO on O were weak during the study period.
In the YRD and PRD regions. Ignoring the principal contributions of anthropogenic emissions, CO significantly impacted the O variations (0.14–0.46 ppb). The varied CO led to surface MDA8 O changes of 0.09 to 0.14 ppb in the YRD and 0.28–0.46 ppb in the PRD by modulating the isoprene emissions and precipitation. On the other hand, the meteorological conditions played a more significant role in surface O than in the NCP, FWP, and SCB regions, resulting in a decrease in MDA8 O from 2008 to 2018 (4.42–3.25 ppb).
In the SCB region. The increase in surface O from 2008 to 2018 was primarily driven by variations in anthropogenic emissions. The effect of meteorological conditions was weak due to the high level of emissions and basin topography. However, the changes in CO significantly impacted surface O levels and were unfavorable to O formation during the study period (3.0 ppb in PreG and PostG).
In conclusion, anthropogenic emissions dominated the O increase in China from 2008 to 2018, and the effects of meteorological conditions on surface O could be more significant in some regions. Furthermore, we emphasize the significance of CO emissions, particularly in southern China, as a critical contributor to O variations. Therefore, it is vital to consider CO variability in future predictions of O concentrations. Such consideration would be helpful for designing long-term O control policies.
Data availability
ERA-Interim data are available at
The supplement related to this article is available online at:
Author contributions
DM: performed experiments; TW: designed the overall research; HW, YQ, JiL, JaL, and SL: reviewed and edited the manuscript; BZ, ML, and MX: contributed to the development of the RegCM-Chem-YIBs model.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements
We would like to acknowledge the anthropogenic emissions inventory support from Tsinghua University and the observed data from the China National Environmental Monitoring Center. We also gratefully acknowledge a wide range of other institutional partners.
Financial support
This work was supported by the National Natural Science Foundation of China (grant nos. 42077192 and 41621005), the National Key Basic Research and Development Program of China (grant nos. 2020YFA0607802 and 2019YFC0214603), the Creative talent exchange program for foreign experts in the Belt and Road countries, and the Emory University–Nanjing University Collaborative Research Grant.
Review statement
This paper was edited by Frank Dentener and reviewed by two anonymous referees.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2023. This work is published under https://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Despite the implementation of the Clean Air Action Plan by the Chinese government in 2013, the issue of increasing surface ozone (O
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details



1 School of Atmospheric Sciences, Nanjing University, Nanjing 210023, China
2 Key Laboratory of Transportation Meteorology of China Meteorological Administration, Nanjing Joint Institute for Atmospheric Sciences, Nanjing 210000, China
3 College of Intelligent Science and Control Engineering, Jinling Institute of Technology, Nanjing 211169, China
4 Key Laboratory for Virtual Geographic Environment of Ministry of Education, Jiangsu Center for Collaborative Innovation in Geographical Information Resource Development and Application, School of Geography Science, Nanjing Normal University, Nanjing 210023, China
5 Department of Geography and Planning, University of Toronto, Toronto M5S 2E8, Canada