Next Article in Journal
Data-Driven Air Quality and Environmental Evaluation for Cattle Farms
Previous Article in Journal
Deactivation Mechanism and Anti-Deactivation Measures of Metal Catalyst in the Dry Reforming of Methane: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Case Study on the Impact of East Asian Summer Monsoon on Surface O3 in China

1
China Academy of Meteorological Sciences, China Meteorological Administration, Beijing 100081, China
2
Department of Earth System Science, Tsinghua University, Beijing 100084, China
3
Key Laboratory of Radiometric Calibration and Validation for Environmental Satellites, National Satellite Meteorological Center (National Center for Space Weather) and Innovation Center for FengYun Meteorological Satellite, China Meteorological Administration, Beijing 100081, China
4
School of Environmental Science and Engineering, Qilu University of Technology (Shandong Academy of Sciences), Jinan 250300, China
*
Authors to whom correspondence should be addressed.
Atmosphere 2023, 14(5), 768; https://doi.org/10.3390/atmos14050768
Submission received: 21 March 2023 / Revised: 18 April 2023 / Accepted: 21 April 2023 / Published: 23 April 2023
(This article belongs to the Section Air Quality)

Abstract

:
The East Asian summer monsoon (EASM) was extremely strong in 2018, which substantially affected surface ozone (O3) in China. Taking 2018 and the average synthesis of 2003 and 2010 to represent the strong and weak EASM cases, respectively, GEOS-Chem with constant anthropogenic emission was employed to investigate the impact of the EASM on surface O3 in the east of China. Simulations show that surface O3 decreased in the northeast and the eastern coast of China and increased in most of the remaining regions during strong EASM. The difference in surface O3 between strong and weak EASM was around −15~7 ppbv. After analyzing relevant meteorological fields, it is found that the decrease in northeast China was mainly attributed to the large increase in vertical upward transport. The considerable decrease in the Huang-Huai-Hai region depended on the dilution and diffusion of eastward anomalous horizontal circulation. The increase in Hunan-Hubei-Guangdong Province was largely due to input from the north. In addition, the vast areas between the Yangtze River and the Yellow River were supported by higher temperatures and stronger shortwave solar radiation that promoted photochemical reactions. The reasons for changes in Shanxi-Sichuan-Yunnan Province were relatively more complex and thus require more in-depth exploration.

1. Introduction

China has a greater level of ozone (O3) exposure to humans and vegetation than other developed nations [1]. Since 2020, O3 has surpassed PM2.5 to take the lead as the most significant air pollutant in the highly urbanized and industrialized regions of China [2]. Surface O3 is a potent oxidant and phytotoxin [3] that poses direct risks to human health [4,5,6] and reduces crop yields [7,8]. O3 may even lessen terrestrial carbon sinks by inhibiting vegetation photosynthesis, causing more severe global warming [9].
Surface O3 is mainly produced by photochemical reactions of precursors, such as volatile organic compounds (VOCs) and nitrogen oxides (NOx), in the presence of solar radiation [10,11,12,13]. Thus, O3 pollution, significantly influenced by precursor emissions, typically occurs in densely populated and industrially developed areas [14,15]. However, high concentrations of O3 and its precursors could be transported by atmospheric circulation over long distances. Therefore, in addition to factors related to photochemical reactions, regional O3 pollution is also regulated by atmospheric circulation transport and meteorological conditions associated with it [16,17].
China, bordered by the Pacific Ocean, has one of the most significant monsoon climates in the world [18]. The East Asian summer monsoon (EASM) is one of the most important components of the Asian monsoon system, which controls summer weather and climate in the area to the east of 100° E in China [19]. In addition, the EASM significantly affects the level and frequency of summertime O3 pollution in the east of China by changing its external-natural factors [20,21,22]. Firstly, the EASM transports relatively clean marine air masses which are, in general, free from O3 and its precursors from the northwest Pacific, favoring the dilution of O3 in southeastern China [23,24]. Secondly, increased cloud cover and precipitation under the influence of humid air masses impedes photochemical reactions, causing O3 concentrations to further decline in coastal areas [25,26]. Finally, the EASM circulation directly changes the horizontal distribution of O3 by transporting the polluted air mass to downwind regions [23]. Pollutants in highly industrialized southeastern China are transported northward by the EASM, contributing to north-south differences in the background level of summertime O3 in the east of China [27,28]. Numerous studies have focused on the effect of the EASM on O3 distribution [22,24,29,30,31]; however, the relevant influencing pathways and dominant processes in different regions have not been thoroughly investigated yet, and are thus worth exploring further.
In this study, we utilize a global three-dimensional model of atmospheric chemistry driven by meteorological input from the Goddard Earth Observing System (GEOS-Chem) to simulate the impact of the EASM on surface O3 concentrations in China during strong and weak EASM years selected by the East Asian summer monsoon index (EASMI). The interference of anthropogenic emission changes was eliminated since anthropogenic emissions were fixed during the simulations.

2. Materials and Methods

2.1. GEOS-Chem Model

GEOS-Chem, developed by Harvard University [32], has been widely used in the study of atmospheric composition and its transports [33,34,35,36]. In this study, GEOS-Chem (v11-01, http://wiki.seas.harvard.edu/geos-chem/index.php/GEOS-Chem_v11-01#v11-01_public_release, accessed on 1 May 2020) was utilized. The nested-grid GEOS-Chem over Asia (60–150° E, 10° S–55° N) with a horizontal resolution of 0.5° × 0.625° (latitude × longitude) was driven by the MERRA-2 reanalysis meteorological dataset, which also served as boundary conditions with a resolution of 2° × 2.5°. The chemical mechanism used in this study is GEOS-Chem’s tropospheric chemistry (Tropchem) mechanism, which was formerly referred to as “NOx-Ox-HC-aerosol-Br”. The Tropchem mechanism includes NOx-Ox-hydrocarbon-aerosol chemistry, Methyl proxy nitrate chemistry, Halogen chemistry, simplified secondary organic aerosol species, and updated isoprene and monoterpene chemistry [37,38]. The linearized O3 (LINOZ) chemical mechanism provided the stratospheric O3 [39]. The non-local scheme was employed for vertical mixing within the boundary layer [40].
The Community Emissions Data System (CEDS) inventory, the default global anthropogenic inventory in GEOS-Chem, was used for global anthropogenic emissions. The MIX inventory for Asia includes the Asian Regional Emissions Inventory (REAS2) for the whole of Asia [41], the Multi-resolution Emission Inventory for China (MEIC) by Tsinghua University, a high-resolution NH3 emission inventory by Peking University [42], an Indian emission inventory developed by Argonne National Laboratory (ANL-India) [43,44], and the official Korean emission inventory from the Clean Air Policy Support System (CAPSS) [45]. The biomass-burning emission inventory was provided by the Global Emission Database v4 (GFED4). Natural emissions were calculated online using the Model of Emissions of Gases and Aerosol from Nature (MEGAN v2.1). Lightning NOx emissions were parameterized by Price and Rind [46] and then driven by Lightning Imaging Sensor/Optical Transient Detector (LIS/OTD) climatology data sets. Both the lightning and soil NOx emissions were calculated online. Here, the anthropogenic emissions were fixed during the different years.

2.2. East Asian Summer Monsoon Index

The EASM has been measured by various EASMI [47,48,49]. Spatial mean Dynamical Normalized Seasonality (DNS) can indicate the seasonal variation in the horizontal wind field, causing it to be used in defining the EASMI [49,50]. The DNS index (DNSI) is given by
D N S I = V 1 ¯ V i V ¯ 2
where V 1 and V i denote horizontal wind vectors in January and month i , respectively. V ¯ is the average horizontal wind vector of January and July. The expression of the norm A is
A = (   | A | 2 d S ) 1 / 2
where S represents the domain of integration. For A at a point ( i , j ), it is defined as
A i , j Δ s [ ( | A i 1 ,   j 2 | + 4 | A i ,   j 2 | + | A i + 1 ,   j 2 | ) c o s φ j + | A i ,   j 1 2 | c o s φ j 1 + | A i ,   j + 1 2 | c o s φ j + 1 ] 1 2  
where φ j is the latitude of the point ( i , j ), Δ s is
Δ s = a Δ φ Δ λ 4
where a represents the mean Earth radius, Δ φ and Δ λ in radian are resolutions at meridional and zonal directions, respectively. Here, the EASMI was defined as summer (June-July-August, JJA) mean DNSI at 850 hPa for East Asia (10–40° N, 110–140° N).
As shown in Figure 1, 2018 was the strongest year for the EASM in at least the last 20 years. Chen, et al. [51] concluded that an extremely strong EASM appeared in 2018 based on multiple definitions of the EASMI [47,52,53]. To capture the characteristics of this extremely strong EASM in 2018, we selected it as a separate case of strong events. To increase the representativeness of weak EASM, we synthesized weak events of different intensities in 2003 and 2010 by using their average to represent a weak EASM.

2.3. Other Surface O3 Data

Observed maximum daily 8-h average (MDA8) O3 from the China National Environmental Monitoring Centre (CNEMC) was used to evaluate the performance of GEOS-Chem. Stations in Beijing, Shanghai, and Guangzhou, where extreme pollution events occur more frequently [54,55,56], were chosen for comparison. There are 12 stations in Beijing, 9 stations in Shanghai, and 11 stations in Guangzhou. A city’s mean MDA8 O3 level was represented by averaging observations from all stations within that city. The simulation data was selected from grids of neighboring stations to match the observations.
Tracking Air Pollution in China (TAP, http://tapdata.org.cn/, accessed on 31 March 2023) provides near real-time MDA8 O3 data products with high a resolution of 10 km from 2013 to the present [57]. TAP integrates multi-source data including ground observations, satellite remote sensing information, high-resolution emission inventory, air quality model simulation, and other ancillary data such as meteorological fields, land use data, population, and elevation.

2.4. Meteorological Data

The monthly mean meteorological data in JJA of 2003, 2010, and 2013–2022 were obtained from the European Centre for Medium-Range Forecasting (ECMWF) Reanalysis v5.0 (ERA5). Data on pressure levels between 1000–500 hPa (16 layers, https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-pressure-levels-monthly-means?tab=form, accessed on 31 March 2023), including u-component (U), v-component (V), and the vertical velocity ( ω ) of wind, as well as geopotential height (GH), were used. Meteorological variables, including 2 m temperature (T-2m), surface solar radiation downwards (SSRD), total precipitation (TP), total cloud cover (TCC), low cloud cover (LCC), surface pressure (SP), 1000 hPa relative humidity (RH), U, V, the speed of wind at 10m (U-10m, V-10m, and WS-10m, respectively), and boundary layer height (BLH), were selected from ERA5 monthly data on single levels (https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels-monthly-means?tab=form, accessed on 31 March 2023). All meteorological variables have a spatial resolution of 0.5° × 0.5°.

2.5. Regression Model

An exhaustive search with linear regression (ES-LR) was utilized to measure the influence of meteorological processes on MDA8 O3. Meteorological processes are divided into three categories: processes related to chemical reactions and horizontal and vertical transport processes. Meteorological parameters associated with each process are listed in Table 1. The linear regression equation is
MDA 8   O 3 t = α 0 + 1 i α i x i ( t ) + ε ( t )
where MDA 8   O 3 t is the detrended time series of JJA monthly mean MDA8 O3 from TAP (2013–2022), x i ( t ) is detrended time series of meteorological parameter x i , ε ( t ) is time series of the residual, α 0 is the intercept, and α i is the fitting coefficient of x i . Explanatory variables were screened from Table 1 by maximizing the adjusted R-square ( R a d j u s t e d 2 ) in ES-LR. R a d j u s t e d 2 represents the fraction of variance in MDA8 O3 explained by the regression model, and its expression is
R a d j u s t e d 2 = 1 ( 1 R 2 ) ( n 1 ) n m 1
where n is the number of samples, m is the number of parameters, and R 2 is the coefficient of determination.

3. Results

3.1. GEOS-Chem Model Evaluation

Figure 2 shows the comparison of observed and simulated MDA8 O3 in July 2015 in Beijing, Shanghai, and Guangzhou. GEOS-Chem performed the worst in Beijing, failing to present the fluctuation in the observation. The simulation and observation had the greatest agreement in Shanghai, particularly during the first three weeks. While not performing as remarkably as in Shanghai, the simulation in Guangzhou captured the approximate trend and extremums of observed MDA8 O3. The correlation coefficients between the simulated MDA8 O3 and the observations were 0.61 (p < 0.01), 0.89 (p < 0.01), and 0.76 (p < 0.01) in Beijing, Shanghai, and Guangzhou, respectively. Overall, GEOS-Chem reproduced the variations of surface O3 in China to some extent, with higher accuracy in eastern and southern China, which is consistent with the results of the previous evaluation [58,59,60]. As a result, the GEOS-Chem simulation results were applied to the following discussion.

3.2. Comparison of Simulated MDA8 O3 during Strong and Weak EASM

Figure 3 illustrates the spatial distribution of simulated MDA8 O3 during strong and weak EASM along with their relative changes. Central and eastern China, as well as the Sichuan Basin, were the most heavily polluted regions. MDA8 O3 in these regions reached approximately 80 ppbv during both weak and strong EASM (Figure 3a,b). Figure 3c reveals obvious changes in the spatial distribution of summertime surface O3 induced by different intensities of EASM. Strong EASM greatly reduced surface O3 mainly in three regions, which we refer to as Regions A−, B−, and C−, respectively (Figure 3c). Region A−, at the junction of southern Gansu Province, eastern Sichuan Province, southwestern Shanxi Province, western Chongqing Municipality, and northwestern Guizhou Province, had the smallest drop (around 3 ppbv) of the three. Covering the southeast of Hebei Province, Shandong-Anhui-Jiangsu Province, and eastern Zhejiang Province, Region B− experienced the largest reduction, reaching −15 ppbv in its coastal areas. Heilongjiang-Jilin-Liaoning Province and the northeast of Inner Mongolia Autonomous Region were divided into Region C−. The decline in MDA8 O3 exceeded 5ppbv in the northeast of Region C−. Surface O3 was increased by strong EASM in most of the remaining areas, particularly in the junction of Sichuan-Yunnan Province (Region A+), the junction of Hubei-Hunan Province (Region B+), and the western part of Guangdong Province (Region C+). The rise in MDA8 O3 reached 5~7 ppbv in these regions. The relative change in MDA8 O3 between strong and weak EASM relative to its level during weak EASM ((strong-weak)/weak) was −20~16% in the east of China. This variation pattern of surface O3 caused by EASM has high consistency with the simulation results of RegCM4-Chem [20] and other studies based on GEOS-Chem [21,31].

3.3. Comparison of Meteorological Conditions during Strong and Weak EASM

The difference in surface O3 between strong and weak EASM cases was mainly attributed to meteorological conditions since the anthropogenic emissions were fixed. The impact of meteorological conditions will be discussed from three perspectives: chemical reactions, horizontal diffusions and transports, and vertical processes. Before quantitatively evaluating changes in meteorological parameters caused by the EASM and qualitatively discussing their impacts on surface O3 in selected regions, we first assessed fractions of variance in regional MDA8 O3 explained by each meteorological process. According to R a d j u s t e d 2 based on ES-LR (Table 2), meteorological parameters totally explained 76%, 72%, 71%, 68%, 79%, and 75% of the summertime MDA8 O3 variation in Regions A+, B+, C+, A−, B−, and C−, respectively. In Regions A+, B+, A−, and B−, chemical reactions were more influential than transports [61], particularly in Region A− [62]. Regions C+ and C− were more affected by horizontal transports [63,64]. Vertical processes have little effect on Regions A+, B+, and A− and considerable effect on Region C+. The following analysis will focus on relative changes in the explanatory variables screened by ES-LR during strong and weak EASM.

3.3.1. Comparison of Chemical Reactions

Figure 4 shows the spatial distribution of the relative difference in meteorological parameters related to chemical reactions. According to Figure 4, strong EASM caused a positive-negative-positive tripolar response of RH, TP, TCC, and LCC, roughly bounded by the Yangtze and Yellow Rivers. With a consistent boundary, SSRD has an opposite response of negative-positive-negative. In general, the meteorological conditions in the regions between the Yangtze River and the Yellow River were favorable for photochemical reactions to produce O3 during strong EASM.
Table 3 lists the spatial mean of relative change in MDA8 O3 and meteorological elements related to photochemical reactions in selected regions. The meteorological conditions in Region B+ were more favorable for photochemical reactions during strong EASM than during weak EASM [61]. Higher T-2m (0.22%) and lower RH (−0.91%) favored the NO2 accumulation [65] and natural emissions of biogenic volatile organic compounds (BVOCs) [66,67]. TCC was lowered by −15.19% and SSRD was increased by 5.61%. Increased temperature and solar radiation further promoted O3 production by accelerating photochemical reactions [10,63]. In addition, considerably decreased TP (−34.29%) weakened the wet scavenging effects of rainfall on O3 precursors and enhanced solar radiation indirectly [68,69,70]. All these changes in meteorological conditions contributed to the increase in local surface O3. In comparison, Region A+ gained fewer benefits from meteorologic-related chemical reactions than Region B+ due to the inhibitory effect of increased RH (0.31%) and TP (7.19%) on chemical reactions [71,72]. However, Region A+ got the largest increase in surface O3 (7.72%) under the premise of a low level of emission [73,74,75], which indicates that Region A+ may be supported by dynamic input more greatly than Region B+ during strong EASM. The contributions of dynamic transport were even more dominant in Regions A−, B−, and C+. Though located in meteorological fields favorable to photochemical reactions with a high level of formaldehyde [76], Regions A− and B− experienced a decline in surface O3, particularly in Region B− (−7.19%). Similarly, given that photochemical reactions were greatly suppressed by changes in meteorological conditions in Region C+, the local surface O3 increase (6.88%) was attributed to dynamic processes. In Region C−, since the impact of higher T-2m (0.19%) was inconsistent with the weaker SSRD (−4.55%) during strong EASM, it is difficult to explain local changes in surface O3 (−5.33%) from the perspective of photochemical reactions regulated by meteorological conditions.

3.3.2. Comparison of Horizontal Diffusions and Transports

Figure 5 shows the spatial distribution of the difference in meteorological elements closely related to diffusions in the boundary layer between strong and weak EASM. The WS-10m (Figure 5b) and BLH (Figure 5c) in Region B− increased by 12.86% and 17.54% (Table 4), respectively, during strong EASM. While the increase in BLH was favorable for the mixing of precursors to produce O3 [77,78], its promotion of outward diffusion contributed more to the decline of surface O3 in Region B− [79,80]. In terms of wind vector at 10m (Figure 5a), the eastward wind anomaly located in Region B−, with a relative increase up to 175.59% in U-10m, brought the relatively clean air mass from the ocean and diluted local surface O3. Similarly, the southern part of Region C− (i.e., Liaoning Province), where V-10m increased by 48.5%, was also affected by the dilution of the northward ocean air. The enhanced southwestward diffusion extended from Region B− toward the south China coast (Figure 5a), further causing surface O3 to decrease in Region B− and increase in Regions B+ and C+ [21,81]. The influence of the EASM on horizontal diffusion conditions within the boundary layer was relatively insignificant in Regions A+ and A−, particularly in the former (Table 4). Complex geography may have played an important role. Previous research has attributed O3 pollution in Sichuan Province to downward transport from the Qinghai-Tibet Plateau [82,83].
Figure 6 further illustrates the large-scale cyclonic anomalous circulation over the Northwest Pacific during strong EASM. The east wind anomalies had incredible impacts on southeastern China, including Regions B+, C+, B−, and A+, where the U-850hPa experienced relative changes by around 200% (Table 4). Figure 6 shows that strong EASM induced a north-south dipole response of GH-500hPa with a transition latitude near 31° N, making the Northwest Pacific Subtropical High (WPSH) weakened. Interestingly, surface O3 in northern and southern China in turn responded to weakened WPSH with a similar dipole-like pattern [67]. Bounded by 32° N, MDA8 O3 rises significantly in southern China, particularly in the southeast part, and drops in northern China during weaker WPSH [67]. Thus, increased O3 in Regions A+, B+, and C+ was also additionally supported by weaker WPSH, which led to a significant rise in natural emissions in southern China [67]. According to Table 4, the relative change in the U-850hPa and V-850hPa in Region C− was notably smaller (−9.41% and −42.44%, respectively) than in other regions, indicating that local O3 change was relatively less dependent on horizontal transport in Region C−.

3.3.3. Comparison of Vertical Processes

The drop of surface O3 in Region C−, particularly in its northern part, which failed to be explained by changes in chemical reactions and horizontal circulation transport (Figure 4, Figure 5 and Figure 6), may be attributed to vertical processes. Figure 7 shows the vertical cross-section of the differences in regional zonal mean V-W between strong and weak EASM in selected regions. The upward W anomaly in Region C− extended from the ground surface to the top of the boundary layer (Figure 7f), which may greatly account for the decline in local MDA8 O3. The relative change in W-1000hPa in Region C− was the largest of the selected 6 regions, up to −244.81%, which was nearly 5.9~21.5 times that of the other regions (Table 5). The upward transport within the boundary layer was also found in Region A−, particularly in areas outside the Basin (Figure 7b). The upward transport in the northern part of Region A was even stronger than in Region C−. However, the downward input between 30–32° N partially offset the loss of surface O3 in Region A− and reduced its decline (Table 5). Strong downward transport in the south part of Region A+ seemed to be an important source of the regional surface O3 (Figure 7a), which supported the hypothesis that the Tibetan Plateau greatly contributed to high O3 in the Basin [83]. In contrast, Regions B+, C+, and B− were dominated by strong horizontal transport during strong EASM (Figure 7c–e), which counteracted the effects of vertical transport on regional O3 (Table 5).

4. Discussion

Based on the similarity in the changes between the surface O3 integrated process rate caused by chemical reactions and its concentration simulated by RegCM4-Chem, previous research pointed out that the chemical process is the uppermost factor in reducing surface O3 levels in the eastern part of China (i.e., Region B− in our study) during summer monsoon seasons [20]. However, our study illustrates the important effects of abnormal horizontal and vertical circulation transport, caused by changes in EASM intensity, on the distribution pattern of surface O3 in central and southeastern China. A variety of reasons can lead to inconsistent outcomes, such as the fact that the objects and methods of these two studies are not identical, the variations in physical and chemical processes between the two models, etc. More in-depth quantitative analysis is required even though our findings are in line with the model analysis [21,31] and viewpoints from space monitoring [30,84]. In addition, the complex response of natural emissions of BVOCs and NOx to the climate system [85] and their impact on surface O3 [86] need to be further explored in future studies.

5. Conclusions

In this study, we first evaluated the performance of GEOS-Chem in simulating summertime surface O3 concentrations against the monitoring data from the CNEMN. The results demonstrate that GEOS-Chem performs better in eastern and southern China. Correlation coefficients between simulations and observations are 0.61 (p < 0.01), 0.89 (p < 0.01), and 0.76 (p < 0.01) in Beijing, Shanghai, and Guangzhou, respectively.
To compare the impact of the EASM on surface O3 in the east of China, GEOS-Chem with fixed anthropogenic emissions was used to simulate summer surface O3 during strong and weak EASM. We found that strong EASM mainly greatly reduced surface O3 in the Heilongjiang-Jilin-Liaoning Province in the northeast, the Huang-Huai-Hai coast in the east, and the junction of Shanxi-Sichuan-Chongqing in the west of China. Surface O3 increased during strong EASM in most of the remaining areas, particularly in the junction of Sichuan-Yunnan Province, the border of Hubei-Hunan Province, and the western part of Guangdong Province. The spatial mean of the relative difference of MDA8 O3 between the strong and weak EASM in these regions reached approximately −7~7% of their levels during weak EASM. Utilizing MDA8 O3 data with high accuracy from TAP, it was confirmed by linear regression that the combination of meteorological parameters could explain 68~79% of summertime surface O3 variability in these regions.
The meteorological process analysis revealed that the decrease in surface O3 of northeast China during the strong EASM was mainly attributed to the large increase in vertical upward transport. The considerable decrease in the Huang-Huai-Hai region depended on the dilution and diffusion of the eastward anomalous horizontal circulation. The increase in Hunan-Hubei-Guangdong Province was largely due to input from the north. In addition, the vast area between the Yangtze River and the Yellow River has also been supported by higher temperatures and stronger shortwave solar radiation that promoted photochemical reactions [61,65]. The reasons for the changes in Shanxi-Sichuan-Yunnan Province were relatively more complex and required more in-depth exploration.

Author Contributions

Conceptualization, L.Z.; methodology, X.Z. (Xin Zhang) and L.Z.; software, L.Z. and L.S.; validation, L.Z.; formal analysis, X.Z. (Xin Zhang); investigation, X.Z. (Xin Zhang); writing—original draft preparation, X.Z. (Xin Zhang); writing—review and editing, X.Z. (Xin Zhang) and L.Z.; visualization, X.Z. (Xin Zhang) and L.Z.; supervision, X.Z. (Xingying Zhang) and Y.L.; funding acquisition, X.Z. (Xingying Zhang) and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 41921005, 41805098 and 41775028.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

ERA5 reanalysis data can be found here: https://cds.climate.copernicus.eu/, accessed on 31 March 2023; TAP data can be found here: http://tapdata.org.cn/, accessed on 31 March 2023.

Acknowledgments

We thank the team of ERA5 for allowing us to use the meteorological data. We thank the team of TAP for allowing us to use the surface ozone data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lu, X.; Hong, J.; Zhang, L.; Cooper, O.R.; Schultz, M.G.; Xu, X.; Wang, T.; Gao, M.; Zhao, Y.; Zhang, Y. Severe Surface Ozone Pollution in China: A Global Perspective. Environ. Sci. Technol. Lett. 2018, 5, 487–494. [Google Scholar] [CrossRef]
  2. MEEPRC. Report on the State of the Ecology and Environment in China. Ministry of Ecology and Environment of the People’s Republic of China. 2022. Available online: https://www.mee.gov.cn/hjzl/sthjzk/zghjzkgb/ (accessed on 14 February 2023).
  3. Lombardozzi, D.; Bonan, G.B.; Sparks, J.P. Integrating O3 Influences on Terrestrial Processes: Photosynthetic and Stomatal Response Data Available for Large-Scale Modeling. Biogeosciences 2013, 10, 6815–6831. [Google Scholar] [CrossRef]
  4. Manisalidis, I.; Stavropoulou, E.; Stavropoulos, A.; Bezirtzoglou, E. Environmental and Health Impacts of Air Pollution: A Review. Front. Public Health 2020, 8, 14. [Google Scholar] [CrossRef]
  5. Brown, J.; Bowman, C. Integrated Science Assessment for Ozone and Related Photochemical Oxidants; NTIS/12280203; US Environmental Protection Agency: Washington, DC, USA, 2013.
  6. Monks, P.S.; Archibald, A.T.; Colette, A.; Cooper, O.; Coyle, M.; Derwent, R.; Fowler, D.; Granier, C.; Law, K.S.; Mills, G.E.; et al. Tropospheric Ozone and Its Precursors from the Urban to the Global Scale from Air Quality to Short-Lived Climate Forcer. Atmos. Chem. Phys. 2015, 15, 8889–8973. [Google Scholar] [CrossRef]
  7. Avnery, S.; Mauzerall, D.L.; Liu, J.; Horowitz, L.W. Global Crop Yield Reductions Due to Surface Ozone Exposure: 1. Year 2000 Crop Production Losses and Economic Damage. Atmos. Environ. 2011, 45, 2284–2296. [Google Scholar] [CrossRef]
  8. Feng, Z.; Xu, Y.; Kobayashi, K.; Dai, L.; Zhang, T.; Agathokleous, E.; Calatayud, V.; Paoletti, E.; Mukherjee, A.; Agrawal, M.; et al. Ozone Pollution Threatens the Production of Major Staple Crops in East Asia. Nat. Food 2022, 3, 47–56. [Google Scholar] [CrossRef]
  9. Sitch, S.; Cox, P.M.; Collins, W.J.; Huntingford, C. Indirect Radiative Forcing of Climate Change through Ozone Effects on the Land-Carbon Sink. Nature 2007, 448, 791–794. [Google Scholar] [CrossRef] [PubMed]
  10. Preiss, P. Photochemical Ozone Formation. Life Cycle Impact Assess. 2015, 115–138. [Google Scholar] [CrossRef]
  11. National Research, Council. Rethinking the Ozone Problem in Urban and Regional Air Pollution; National Academies Press: Washington, DC, USA, 1992. [Google Scholar]
  12. Chameides, W.L.; Walker, J.C. A Time-Dependent Photochemical Model for Ozone near the Ground. J. Geophys. Res. 1976, 81, 413–420. [Google Scholar] [CrossRef]
  13. Tie, X.; Brasseur, G.; Ying, Z. Impact of Model Resolution on Chemical Ozone Formation in Mexico City: Application of the WRF-Chem Model. Atmos. Chem. Phys. 2010, 10, 8983–8995. [Google Scholar] [CrossRef]
  14. Wang, T.; Xue, L.; Brimblecombe, P.; Lam, Y.F.; Li, L.; Zhang, L. Ozone Pollution in China: A Review of Concentrations, Meteorological Influences, Chemical Precursors, and Effects. Sci. Total Environ. 2017, 575, 1582–1596. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, Y.; Cooper, O.R.; Gaudel, A.; Thompson, A.M.; Nédélec, P.; Ogino, S.Y.; West, J.J. Tropospheric Ozone Change from 1980 to 2010 Dominated by Equatorward Redistribution of Emissions. Nat. Geosci. 2016, 9, 875–879. [Google Scholar] [CrossRef] [PubMed]
  16. Lu, X.; Zhang, L.; Shen, L. Meteorology and Climate Influences on Tropospheric Ozone: A Review of Natural Sources, Chemistry, and Transport Patterns. Curr. Pollut. Rep. 2019, 5, 238–260. [Google Scholar] [CrossRef]
  17. Derwent, R.G.; Stevenson, D.S.; Collins, W.J.; Johnson, C.E. Intercontinental Transport and the Origins of the Ozone Observed at Surface Sites in Europe. Atmos. Environ. 2004, 38, 1891–1901. [Google Scholar] [CrossRef]
  18. Wang, B. The Asian Monsoon; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  19. Zhang, R. Changes in East Asian Summer Monsoon and Summer Rainfall over Eastern China During Recent Decades. Sci. Bull. 2015, 60, 1222–1224. [Google Scholar] [CrossRef]
  20. Li, S.; Wang, T.; Huang, X.; Pu, X.; Li, M.; Chen, P.; Yang, X.Q.; Wang, M. Impact of East Asian Summer Monsoon on Surface Ozone Pattern in China. J. Geophys. Res. Atmos. 2018, 123, 1401–1411. [Google Scholar] [CrossRef]
  21. Zhou, Y.; Yang, Y.; Wang, H.; Wang, J.; Li, M.; Li, H.; Wang, P.; Zhu, J.; Li, K.; Liao, H. Summer Ozone Pollution in China Affected by the Intensity of Asian Monsoon Systems. Sci. Total Environ. 2022, 849, 157785. [Google Scholar] [CrossRef]
  22. He, Y.J.; Uno, I.; Wang, Z.F.; Pochanart, P.; Li, J.; Akimoto. Significant Impact of the East Asia Monsoon on Ozone Seasonal Behavior in the Boundary Layer of Eastern China and the West Pacific Region. Atmos. Chem. Phys. 2008, 8, 7543–7555. [Google Scholar] [CrossRef]
  23. Zhou, D.; Ding, A.; Mao, H.; Fu, C.; Wang, T.; Chan, L.Y.; Ding, K.; Zhang, Y.; Liu, J.; Lu, A.; et al. Impacts of the East Asian Monsoon on Lower Tropospheric Ozone over Coastal South China. Environ. Res. Lett. 2013, 8, 044011. [Google Scholar] [CrossRef]
  24. Tang, H.; Liu, G.; Zhu, J.; Han, Y.; Kobayashi, K. Seasonal Variations in Surface Ozone as Influenced by Asian Summer Monsoon and Biomass Burning in Agricultural Fields of the Northern Yangtze River Delta. Atmos. Res. 2013, 122, 67–76. [Google Scholar] [CrossRef]
  25. Zhao, Z.; Wang, Y. Influence of the West Pacific Subtropical High on Surface Ozone Day to Day Variability in Summertime over the Yangtze River Delta Region. AGU Fall Meet. Abstr. 2016, 2016, A31B-0033. [Google Scholar]
  26. Han, H.; Liu, J.; Shu, L.; Wang, T.; Yuan, H. Local and Synoptic Meteorological Influences on Daily Variability in Summertime Surface Ozone in Eastern China. Atmos. Chem. Phys. 2020, 20, 203–222. [Google Scholar] [CrossRef]
  27. Ding, A.J.; Wang, T.; Thouret, V.; Cammas, J.-P.; Nédélec, P. Tropospheric Ozone Climatology over Beijing: Analysis of Aircraft Data from the MOZAIC Program. Atmos. Chem. Phys. 2008, 8, 1–13. [Google Scholar] [CrossRef]
  28. Liu, N.; Lin, W.; Ma, J.; Xu, W.; Xu, X. Seasonal Variation in Surface Ozone and Its Regional Characteristics at Global Atmosphere Watch Stations in China. J. Environ. Sci. 2019, 77, 291–302. [Google Scholar] [CrossRef] [PubMed]
  29. Worden, J.; Jones, D.B.; Liu, J.; Parrington, M.; Bowman, K.; Stajner, I.; Beer, R.; Jiang, J.; Thouret, V.; Kulawik, S. Observed Vertical Distribution of Tropospheric Ozone During the Asian Summertime Monsoon. J. Geophys. Res. Atmos. 2019, 114. [Google Scholar] [CrossRef]
  30. Zhao, C.; Wang, Y.; Yang, Q.; Fu, R.; Cunnold, D.; Choi, Y. Impact of East Asian Summer Monsoon on the Air Quality over China: View from Space. J. Geophys. Res. Atmos. 2010, 115. [Google Scholar] [CrossRef]
  31. Yang, Y.; Liao, H.; Li, J. Impacts of the East Asian Summer Monsoon on Interannual Variations of Summertime Surface-Layer Ozone Concentrations over China. Atmos. Chem. Phys. 2014, 14, 6867–6879. [Google Scholar] [CrossRef]
  32. Bey, I.; Jacob, D.; Yantosca, R.; Logan, J.; Field, B.; Fiore, A.; Li, Q.; Liu, H.; Mickley, L.; Schultz, M. Global Modeling of Tropospheric Chemistry with Assimilated Meteorology: Model Description and Evaluation. J. Geophys. Res. Atmos. 2001, 106, 23073–92305. [Google Scholar] [CrossRef]
  33. Fusco, A.C.; Logan, J.A. Analysis of 1970–1995 Trends in Tropospheric Ozone at Northern Hemisphere Midlatitudes with the GEOS-Chem Model. J. Geophys. Res. Atmos. 2003, 108. [Google Scholar] [CrossRef]
  34. Nassar, R.; Logan, J.A.; Megretskaia, I.A.; Murray, L.T.; Zhang, L.; Jones, D.B. Analysis of Tropical Tropospheric Ozone, Carbon Monoxide, and Water Vapor During the 2006 El Niño Using TES Observations and the GEOS-Chem Model. J. Geophys. Res. Atmos. 2009, 114. [Google Scholar] [CrossRef]
  35. Hu, L.; Jacob, D.J.; Liu, X.; Zhang, Y.; Zhang, L.; Kim, P.S.; Sulprizio, M.P.; Yantosca, R.M. Global Budget of Tropospheric Ozone: Evaluating Recent Model Advances with Satellite (OMI), Aircraft (IAGOS), and Ozonesonde Observations. Atmos. Environ. 2017, 167, 323–334. [Google Scholar] [CrossRef]
  36. Potts, D.A.; Marais, E.A.; Boesch, H.; Pope, R.J.; Lee, J.; Drysdale, W.; Chipperfield, M.P.; Kerridge, B.; Siddans, R.; Moore, D.P.; et al. Diagnosing Air Quality Changes in the UK During the COVID-19 Lockdown Using TROPOMI and GEOS-Chem. Environ. Res. Lett. 2021, 16, 054031. [Google Scholar] [CrossRef]
  37. Mao, J.; Jacob, D.J.; Evans, M.J.; Olson, J.R.; Ren, X.; Brune, W.H.; St Clair, J.M.; Crounse, J.D.; Spencer, K.M.; Beaver, M.R.; et al. Chemistry of Hydrogen Oxide Radicals (HOx) in the Arctic Troposphere in Spring. Atmos. Chem. Phys. 2010, 10, 5823–5838. [Google Scholar] [CrossRef]
  38. Lin, J.T.; Liu, Z.; Zhang, Q.; Liu, H.; Mao, J.; Zhuang, G. Modeling Uncertainties for Tropospheric Nitrogen Dioxide Columns Affecting Satellite-Based Inverse Modeling of Nitrogen Oxides Emissions. Atmos. Chem. Phys. 2012, 12, 12255–12275. [Google Scholar] [CrossRef]
  39. McLinden, C.A.; Olsen, S.C.; Hannegan, B.; Wild, O.; Prather, M.J.; Sundet, J. Stratospheric Ozone in 3-D Models: A Simple Chemistry and the Cross-Tropopause Flux. J. Geophys. Res. Atmos. 2000, 105, 14653–14665. [Google Scholar] [CrossRef]
  40. Lin, J.T.; McElroy, M.B. Impacts of Boundary Layer Mixing on Pollutant Vertical Profiles in the Lower Troposphere: Implications to Satellite Remote Sensing. Atmos. Environ. 2010, 44, 1726–1739. [Google Scholar] [CrossRef]
  41. Kurokawa, J.; Ohara, T.; Morikawa, T.; Hanayama, S.; Janssens-Maenhout, G.; Fukui, T.; Kawashima, K.; Akimoto, H. Emissions of Air Pollutants and Greenhouse Gases over Asian Regions During 2000–2008: Regional Emission Inventory in Asia (Reas) Version 2. Atmos. Chem. Phys. 2013, 13, 11019–11058. [Google Scholar] [CrossRef]
  42. Huang, X.; Song, Y.; Li, M.; Li, J.; Huo, Q.; Cai, X.; Zhu, T.; Hu, M.; Zhang, H. A High-Resolution Ammonia Emission Inventory in China. Glob. Biogeochem. Cycles 2012, 26. [Google Scholar] [CrossRef]
  43. Lu, Z.; Streets, D.G. Increase in NOx Emissions from Indian Thermal Power Plants During 1996–2010: Unit-Based Inventories and Multisatellite Observations. Environ. Sci. Technol. 2012, 46, 7463–7470. [Google Scholar] [CrossRef]
  44. Lu, Z.; Zhang, Q.; Streets, D.G. Sulfur Dioxide and Primary Carbonaceous Aerosol Emissions in China and India, 1996–2010. Atmos. Chem. Phys. 2011, 11, 9839–9864. [Google Scholar] [CrossRef]
  45. Lee, D.G.; Lee, Y.M.; Jang, K.W.; Yoo, C.; Kang, K.H.; Lee, J.H.; Jung, S.W.; Park, J.M.; Lee, S.B.; Han, J.S.; et al. Korean National Emissions Inventory System and 2007 Air Pollutant Emissions. Asian. J. Atmos. Environ. 2011, 5, 278–291. [Google Scholar] [CrossRef]
  46. Price, C.; Rind, D. A Simple Lightning Parameterization for Calculating Global Lightning Distributions. J. Geophys. Res. Atmos. 1992, 97, 9919–9933. [Google Scholar] [CrossRef]
  47. Zhang, Q.Y.; Tao, S.Y.; Chen, L.T. The Inter-Annual Variability of East Asian Summer Monsoon Indices and Its Association with the Pattern of General Circulation over East Asia. Acta Meteorol. Sin. 2003, 61, 559–569. [Google Scholar]
  48. Shi, N.; Lu, J.-J.; Zhu, Q.-G. East Asian Winter/Summer Monsoon Intensity Indices with Their Climatic Change in 1873–1989. Nanjing Inst. Meteor. 1996, 19, 168–177. [Google Scholar]
  49. Li, J.; Zeng, Q. A Unified Monsoon Index. Geophys. Res. Lett. 2002, 29, 115-1–115-4. [Google Scholar] [CrossRef]
  50. Li, J.P.; Zeng, Q.C. A New Monsoon Index, Its Interannual Variability and Relation with Monsoon Precipitation. Clim. Environ. Res. 2005, 10, 351–365. [Google Scholar]
  51. Chen, L.; Gu, W.; Li, W. Why Is the East Asian Summer Monsoon Extremely Strong in 2018? —Collaborative Effects of SST and Snow Cover Anomalies. J. Meteorol. Res. 2019, 33, 593–608. [Google Scholar] [CrossRef]
  52. Shi, N.; Zhu, Q.; Wu, B. The East Asian Summer Monsoon in Relation to Summer Large Scale Weather-Climate Anomaly in China for Last 40 Years. Chin. J. Atmos. Sci. 1996, 20, 575–583. (In Chinese) [Google Scholar]
  53. Zhu, C.W.; He, J.H.; Wu, G.X. East Asian Monsoon Index and Its Inter-Annual Relationship with Largescale Thermal Dynamic Circulation. Acta Meteor. Sin. 2000, 58, 391–402. (In Chinese) [Google Scholar]
  54. Fu, Y.; Liao, H.; Yang, Y. Interannual and Decadal Changes in Tropospheric Ozone in China and the Associated Chemistry-Climate Interactions: A Review. Adv. Atmos. Sci. 2019, 36, 975–993. [Google Scholar] [CrossRef]
  55. Huang, J.; Zhou, C.; Lee, X.; Bao, Y.; Zhao, X.; Fung, J.; Richter, A.; Liu, X.; Zheng, Y. The Effects of Rapid Urbanization on the Levels in Tropospheric Nitrogen Dioxide and Ozone over East China. Atmos. Environ. 2013, 77, 558–567. [Google Scholar] [CrossRef]
  56. Dufour, G.; Eremenko, M.; Beekmann, M.; Cuesta, J.; Foret, G.; Fortems-Cheiney, A.; Lachâtre, M.; Lin, W.; Liu, Y.; Xu, X.; et al. Lower Tropospheric Ozone over the North China Plain: Variability and Trends Revealed by IASI Satellite Observations for 2008–2016. Atmos. Chem. Phys. 2018, 18, 16439–16459. [Google Scholar] [CrossRef]
  57. Xue, T.; Zheng, Y.; Geng, G.; Xiao, Q.; Meng, X.; Wang, M.; Li, X.; Wu, N.; Zhang, Q.; Zhu, T. Estimating Spatiotemporal Variation in Ambient Ozone Exposure During 2013–2017 Using a Data-Fusion Model. Environ. Sci. Technol. 2020, 54, 14877–14888. [Google Scholar] [CrossRef] [PubMed]
  58. Lou, S.; Liao, H.; Yang, Y.; Mu, Q. Simulation of the Interannual Variations of Tropospheric Ozone over China: Roles of Variations in Meteorological Parameters and Anthropogenic Emissions. Atmos. Environ. 2015, 122, 839–851. [Google Scholar] [CrossRef]
  59. Lu, X.; Zhang, L.; Wu, T.; Long, M.S.; Wang, J.; Jacob, D.J.; Zhang, F.; Zhang, J.; Eastham, S.D.; Hu, L.; et al. Development of the Global Atmospheric Chemistry General Circulation Model BCC-GEOS-Chem V1. 0: Model Description and Evaluation. Geosci. Model Dev. 2020, 13, 3817–3838. [Google Scholar] [CrossRef]
  60. David, L.M.; Ravishankara, A.R.; Brewer, J.F.; Sauvage, B.; Thouret, V.; Venkataramani, S.; Sinha, V. Tropospheric Ozone over the Indian Subcontinent from 2000 to 2015: Data Set and Simulation Using GEOS-Chem Chemical Transport Model. Atmos. Environ. 2019, 219, 117039. [Google Scholar] [CrossRef]
  61. Xia, Y.; Hu, Y.; Huang, Y.; Bian, J.; Zhao, C.; Wei, J.; Yan, Y.; Xie, F.; Lin, J. Concurrent Hot Extremes and High Ultraviolet Radiation in Summer over the Yangtze Plain and Their Possible Impact on Surface Ozone. Environ. Res. Lett. 2022, 17, 064001. [Google Scholar] [CrossRef]
  62. Zhao, S.; Yu, Y.; Yin, D.; Qin, D.; He, J.; Dong, L. Spatial Patterns and Temporal Variations of Six Criteria Air Pollutants During 2015 to 2017 in the City Clusters of Sichuan Basin, China. Sci. Total Environ. 2018, 624, 540–557. [Google Scholar] [CrossRef]
  63. Chen, X.; Zhong, B.; Huang, F.; Wang, X.; Sarkar, S.; Jia, S.; Deng, X.; Chen, D.; Shao, M. The Role of Natural Factors in Constraining Long-Term Tropospheric Ozone Trends over Southern China. Atmos. Environ. 2020, 220, 117060. [Google Scholar] [CrossRef]
  64. Liang, Z.; Ju, T.; Dong, H.; Geng, T.; Duan, J.; Huang, R. Study on the Variation Characteristics of Tropospheric Ozone in Northeast China. Environ. Monit. Assess. 2021, 193, 1–16. [Google Scholar] [CrossRef]
  65. Camalier, L.; Cox, W.; Dolwick, P. The Effects of Meteorology on Ozone in Urban Areas and Their Use in Assessing Ozone Trends. Atmos. Environ. 2007, 41, 7127–7137. [Google Scholar] [CrossRef]
  66. Bai, J.; Guenther, A.; Turnipseed, A.; Duhl, T. Seasonal and Interannual Variations in Whole-Ecosystem Isoprene and Monoterpene Emissions from a Temperate Mixed Forest in Northern China. Atmos. Pollut. Res. 2015, 6, 696–707. [Google Scholar] [CrossRef]
  67. Jiang, Z.; Li, J.; Lu, X.; Gong, C.; Zhang, L.; Liao, H. Impact of Western Pacific Subtropical High on Ozone Pollution over Eastern China. Atmos. Chem. Phys. 2021, 21, 2601–2613. [Google Scholar] [CrossRef]
  68. Yoo, J.M.; Lee, Y.R.; Kim, D.; Jeong, M.J.; Stockwell, W.R.; Kundu, P.K.; Oh, S.M.; Shin, D.B.; Lee, S.J. New Indices for Wet Scavenging of Air Pollutants (O3, CO, NO2, SO2, and PM10) by Summertime Rain. Atmos. Environ. 2014, 82, 226–237. [Google Scholar] [CrossRef]
  69. Cuchiara, G.C.; Fried, A.; Barth, M.C.; Bela, M.M.; Homeyer, C.R.; Walega, J.; Weibring, P.; Richter, D.; Woods, S.; Beyersdorf, A.; et al. Effect of Marine and Land Convection on Wet Scavenging of Ozone Precursors Observed During a Seac4rs Case Study. J. Geophys. Res. Atmos. 2023, 128, e2022JD037107. [Google Scholar] [CrossRef]
  70. Shan, W.; Yin, Y.; Zhang, J.; Ding, Y. Observational Study of Surface Ozone at an Urban Site in East China. Atmos. Res. 2008, 89, 252–261. [Google Scholar] [CrossRef]
  71. Guo, L.C.; Zhang, Y.; Lin, H.; Zeng, W.; Liu, T.; Xiao, J.; Rutherford, S.; You, J.; Ma, W. The Washout Effects of Rainfall on Atmospheric Particulate Pollution in Two Chinese Cities. Environ. Pollut. 2016, 215, 195–202. [Google Scholar] [CrossRef]
  72. Elperin, T.; Fominykh, A.; Krasovitov, B.; Vikhansky, A. Effect of Rain Scavenging on Altitudinal Distribution of Soluble Gaseous Pollutants in the Atmosphere. Atmos. Environ. 2011, 45, 2427–2433. [Google Scholar] [CrossRef]
  73. Wu, K.; Yang, X.; Chen, D.; Gu, S.; Lu, Y.; Jiang, Q.; Wang, K.; Ou, Y.; Qian, Y.; Shao, P.; et al. Estimation of Biogenic VOC Emissions and Their Corresponding Impact on Ozone and Secondary Organic Aerosol Formation in China. Atmos. Res. 2020, 231, 104656. [Google Scholar] [CrossRef]
  74. Ren, J.; Guo, F.; Xie, S. Diagnosing Ozone–NOx–VOC Sensitivity and Revealing Causes of Ozone Increases in China Based on 2013–2021 Satellite Retrievals. Atmos. Chem. Phys. 2022, 22, 15035–15047. [Google Scholar] [CrossRef]
  75. Dou, X.; Liao, C.; Wang, H.; Huang, Y.; Tu, Y.; Huang, X.; Peng, Y.; Zhu, B.; Tan, J.; Deng, Z.; et al. Estimates of Daily Ground-Level No2 Concentrations in China Based on Random Forest Model Integrated K-Means. Adv. Appl. Energy 2021, 2, 100017. [Google Scholar] [CrossRef]
  76. Wenjing, S.; Liu, C.; Chan, K.L.; Hu, Q.; Liu, H.; Ji, X.; Zhu, Y.; Liu, T.; Zhang, C.; Chen, Y. An Improved TROPOMI Tropospheric HCHO Retrieval over China. Atmos. Meas. Tech. 2020, 13, 6271–6292. [Google Scholar]
  77. Zhao, W.; Fan, S.; Guo, H.; Gao, B.; Sun, J.; Chen, L. Assessing the Impact of Local Meteorological Variables on Surface Ozone in Hong Kong During 2000–2015 Using Quantile and Multiple Line Regression Models. Atmos. Environ. 2016, 144, 182–193. [Google Scholar] [CrossRef]
  78. White, A.B.; Darby, L.; Senff, C.; King, C.; Banta, R.; Koermer, J.; Wilczak, J.; Neiman, P.; Angevine, W.; Talbot, R. Comparing the Impact of Meteorological Variability on Surface Ozone During the Neaqs (2002) and Icartt (2004) Field Campaigns. J. Geophys. Res. Atmos. 2007, 112. [Google Scholar] [CrossRef]
  79. Miao, Y.; Che, H.; Zhang, X.; Liu, S. Relationship between Summertime Concurring PM2. 5 and O3 Pollution and Boundary Layer Height Differs between Beijing and Shanghai, China. Environ. Pollut. 2021, 268, 115775. [Google Scholar] [CrossRef]
  80. Stull, R. An Introduction to Boundary Layer Meteorology; Springer Science & Business Media: Berlin/Heidelberg, Germany, 1988. [Google Scholar]
  81. Jinhai; Ju, J.; Wen, Z.; Lü, J.; Jin, Q. A Review of Recent Advances in Research on Asian Monsoon in China. Adv. Atmos. Sci. 2007, 24, 972–992. [Google Scholar]
  82. Yang, X.; Wu, K.; Wang, H.; Liu, Y.; Gu, S.; Lu, Y.; Zhang, X.; Hu, Y.; Ou, Y.; Wang, S.; et al. Summertime Ozone Pollution in Sichuan Basin, China: Meteorological Conditions, Sources and Process Analysis. Atmos. Environ. 2020, 226, 117392. [Google Scholar] [CrossRef]
  83. Zhao, S.; Yu, Y.; Qin, D.; Yin, D.; Dong, L.; He, J. Analyses of Regional Pollution and Transportation of PM2.5 and Ozone in the City Clusters of Sichuan Basin, China. Atmos. Pollut. Res. 2019, 10, 374–385. [Google Scholar] [CrossRef]
  84. Safieddine, S.; Boynard, A.; Hao, N.; Huang, F.; Wang, L.; Ji, D.; Barret, B.; Ghude, D.S.; Coheur, P.-F.; Hurtmans, D. Tropospheric Ozone Variability During the East Asian Summer Monsoon as Observed by Satellite (IASI), Aircraft (MOZAIC) and Ground Stations. Atmos. Chem. Phys. 2016, 16, 10489–10500. [Google Scholar] [CrossRef]
  85. Hantson, S.; Knorr, W.; Schurgers, G.; Thomas; Pugh, A.; Arneth, A. Global Isoprene and Monoterpene Emissions under Changing Climate, Vegetation, CO2 and Land Use. Atmos. Environ. 2017, 155, 35–45. [Google Scholar] [CrossRef]
  86. Xiao, L.; Zhang, L.; Chen, Y.; Zhou, M.; Zheng, B.; Li, K.; Liu, Y.; Lin, J.; Fu, T.-M.; Zhang, Q. Exploring 2016–2017 Surface Ozone Pollution over China: Source Contributions and Meteorological Influences. Chemistry. A. Physics 2019, 19, 8339–8361. [Google Scholar]
Figure 1. Normalized EASMI from 2002 to 2018.
Figure 1. Normalized EASMI from 2002 to 2018.
Atmosphere 14 00768 g001
Figure 2. The comparison of observed and simulated daily MDA8 O3 in July 2015 in Beijing, Shanghai, and Guangzhou. The marked error bar represents the standard deviation.
Figure 2. The comparison of observed and simulated daily MDA8 O3 in July 2015 in Beijing, Shanghai, and Guangzhou. The marked error bar represents the standard deviation.
Atmosphere 14 00768 g002
Figure 3. Spatial distribution of MDA8 O3 simulated during weak and strong EASM with their difference. JJA mean MDA8 O3 of (a) the average of 2003 and 2010 (hereinafter referred to as the weak year) and (b) 2018 (hereinafter referred to as the strong year). (c) The relative change in MDA8 O3 between strong and weak EASM relative to its level during weak EASM ((strong−weak)/weak). The red (blue) box represents areas where MDA8 O3 is significantly reduced (increased), and they are marked by A−, B−, and C−, (A+, B+, and C+), respectively, the same below.
Figure 3. Spatial distribution of MDA8 O3 simulated during weak and strong EASM with their difference. JJA mean MDA8 O3 of (a) the average of 2003 and 2010 (hereinafter referred to as the weak year) and (b) 2018 (hereinafter referred to as the strong year). (c) The relative change in MDA8 O3 between strong and weak EASM relative to its level during weak EASM ((strong−weak)/weak). The red (blue) box represents areas where MDA8 O3 is significantly reduced (increased), and they are marked by A−, B−, and C−, (A+, B+, and C+), respectively, the same below.
Atmosphere 14 00768 g003
Figure 4. Spatial distribution of relative difference ((strong−weak)/weak) in meteorological elements related to chemical reactions, including (a) 2 m temperature (T-2m), (b) 1000 hPa relative humidity (RH), (c) total precipitation (TP), (d) surface solar radiation downwards (SSRD), (e) total cloud cover (TCC), and (f) low cloud cover (LCC).
Figure 4. Spatial distribution of relative difference ((strong−weak)/weak) in meteorological elements related to chemical reactions, including (a) 2 m temperature (T-2m), (b) 1000 hPa relative humidity (RH), (c) total precipitation (TP), (d) surface solar radiation downwards (SSRD), (e) total cloud cover (TCC), and (f) low cloud cover (LCC).
Atmosphere 14 00768 g004
Figure 5. Spatial distribution of the difference in meteorological elements related to diffusions in the boundary layer between strong and weak EASM, including absolute changes in (a) surface pressure (SP) and wind vector at 10m and relative changes in (b) 10 m wind speed (WS-10m) and (c) boundary layer height.
Figure 5. Spatial distribution of the difference in meteorological elements related to diffusions in the boundary layer between strong and weak EASM, including absolute changes in (a) surface pressure (SP) and wind vector at 10m and relative changes in (b) 10 m wind speed (WS-10m) and (c) boundary layer height.
Atmosphere 14 00768 g005
Figure 6. Spatial distribution of the difference in GH-500hPa and wind vector at 850 hPa between strong and weak EASM.
Figure 6. Spatial distribution of the difference in GH-500hPa and wind vector at 850 hPa between strong and weak EASM.
Atmosphere 14 00768 g006
Figure 7. Vertical cross-section of difference (strong−weak) in regional zonal mean GH (shaded) and zonal mean V-W (arrows) within (a) Region A+, (b) Region A−, (c) Region B+, (d) Region B−, (e) Region C+, and (f) Region C−. The change in W was magnified by a factor of 500 for easy visualization. The red and blue solid lines represent relative change ((strong − weak)/weak, %) in zonal mean MDA8 O3 and BLH, respectively. The black solid line shows the average level of zonal mean BLH during strong and weak EASM ((strong + weak)/2, hPa).
Figure 7. Vertical cross-section of difference (strong−weak) in regional zonal mean GH (shaded) and zonal mean V-W (arrows) within (a) Region A+, (b) Region A−, (c) Region B+, (d) Region B−, (e) Region C+, and (f) Region C−. The change in W was magnified by a factor of 500 for easy visualization. The red and blue solid lines represent relative change ((strong − weak)/weak, %) in zonal mean MDA8 O3 and BLH, respectively. The black solid line shows the average level of zonal mean BLH during strong and weak EASM ((strong + weak)/2, hPa).
Atmosphere 14 00768 g007
Table 1. Categories of meteorological processes affecting MDA8 O3 and their associated meteorological variables.
Table 1. Categories of meteorological processes affecting MDA8 O3 and their associated meteorological variables.
Process CategoriesRelated Meteorological Variables
chemical reactionT-2M, RH, TP, SSRD, TCC, and LCC
horizontal transportSP, U-10m, V-10m, WS-10m, BLH, GH-500hPa 1, U-850hPa 1, and V-850hPa 1
vertical transportW-950hPa 1,2, W-975hPa 1,2, and W-1000hPa 1,2
1 The suffix specifies the pressure layer where the variable is located. 2 To facilitate vector synthesis, ω (Pa/s) is converted to W (m/s).
Table 2. Explanatory variables screened by ES-LR used in linear regression models for different meteorological processes and corresponding R a d j u s t e d 2 .
Table 2. Explanatory variables screened by ES-LR used in linear regression models for different meteorological processes and corresponding R a d j u s t e d 2 .
RegionProcess CategoriesExplanatory Variables
R a d j u s t e d 2
A+chemical RH + TP + TCC0.51
horizontal V-10m + BLH + GH-500hPa + V-850hPa0.37
vertical W-975hPa + W-950hPa0.00
totalTCC + LCC + SP + U-10m + BLH + GH-500hPa + U-850hPa + W-950hPa0.76
B+chemical T-2m + TP + SSRD + TCC0.42
horizontal V-10m + GH-500hPa + U-850hPa0.39
vertical W-1000hPa + W-975hPa + W-950hPa0.03
totalRH + LCC + U-10m + WS-10m + BLH + GH-500hPa + V-850hPa0.72
C+chemical RH + TP0.26
horizontal U-10m + GH-500hPa + U-850hPa + V-850hPa0.56
vertical W-975hPa0.46
totalT-2m + RH + TCC + SP + U-10m + WS-10m + BLH + GH-500hPa + U-850hPa + W-1000hPa0.71
A−chemical T-2m + TP + SSRD + TCC + LCC0.46
horizontal V-10m + WS-10m + BLH + V-850hPa0.05
vertical W-1000hPa + W-975hPa0.00
totalT-2m + RH + SSRD + TCC + LCC + V-10m + WS-10m + V-850hPa + W-975hPa0.68
B−chemical T-2m + RH + SSRD0.69
horizontal U-10m + V-10m + BLH + U-850hPa + V-850hPa0.67
vertical W-1000hPa + W-950hPa0.28
totalRH + TP + SP + U-10m + V-10m + BLH + U-850hPa + W-1000hPa + W-975hPa + W-950hPa0.79
C−chemical T-2m + RH + SSRD + LCC0.68
horizontal U-10m + WS-10m + BLH + GH-500hPa + U-850hPa0.71
vertical W-975hPa0.18
totalRH + TP + U-850hPa + W-950hPa0.75
Table 3. The spatial mean of relative changes ((strong − weak)/weak, %) in MDA8 O3 and meteorological elements related to chemical reactions in selected regions with significant changes in surface O3.
Table 3. The spatial mean of relative changes ((strong − weak)/weak, %) in MDA8 O3 and meteorological elements related to chemical reactions in selected regions with significant changes in surface O3.
A+B+C+A−B−C−
MDA8 O37.72 6.92 6.88 −1.02 −7.19 −5.33
T-2m0.06 0.22 −0.08 0.30 0.45 0.19
RH0.31 −0.91 2.34 −0.05 −1.18 3.65
TP7.19 −34.29 42.24 −6.66 −18.68 14.94
SSRD3.20 5.61 −8.59 8.35 4.70 −4.55
TCC−3.27 −15.19 10.27 −9.35 −16.79 0.62
LCC−2.04 −21.63 5.25 −13.27 −15.62 14.65
Table 4. The spatial mean of relative changes ((strong − weak)/weak, %) in MDA8 O3 and meteorological elements related to horizontal diffusions and transports in selected regions with significant changes in surface O3.
Table 4. The spatial mean of relative changes ((strong − weak)/weak, %) in MDA8 O3 and meteorological elements related to horizontal diffusions and transports in selected regions with significant changes in surface O3.
A+B+C+A−B−C−
MDA8 O37.726.926.88−1.02−7.19−5.33
SP−0.07−0.12−0.31−0.08−0.120.01
U-10m−19.62−10.0890673.13175.59−2.81
V-10m15.78−114.89−60.650.76−7.2548.5
WS-10m−8.42−7.57−7.25−4.6212.861.21
BLH−2.9−1.01−11.820.5717.54−2.34
GH-500hPa−0.07−0.03−0.320.060.270.51
U-850hPa33.24−196.81−202.17269.72229.43−9.41
V-850hPa−30.66−71.55−36.7993.8813.36−42.44
Table 5. The spatial mean of relative changes ((strong − weak)/weak, %) in MDA8 O3 and W within the boundary layer in selected regions with significant changes in surface O3.
Table 5. The spatial mean of relative changes ((strong − weak)/weak, %) in MDA8 O3 and W within the boundary layer in selected regions with significant changes in surface O3.
A+B+C+A−B−C−
MDA8 O37.72 6.92 6.88 −1.02 −7.19 −5.33
W-950hPa−25.68 53.40 −125.69 −3.52 −413.15 −47.08
W-975hPa−14.14 9.44 −81.76 −24.87 549.37 −9.38
W-1000hPa13.24 −11.39 −50.35 −13.72 11.83 −244.81
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, X.; Zhou, L.; Zhang, X.; Luo, Y.; Sun, L. A Case Study on the Impact of East Asian Summer Monsoon on Surface O3 in China. Atmosphere 2023, 14, 768. https://doi.org/10.3390/atmos14050768

AMA Style

Zhang X, Zhou L, Zhang X, Luo Y, Sun L. A Case Study on the Impact of East Asian Summer Monsoon on Surface O3 in China. Atmosphere. 2023; 14(5):768. https://doi.org/10.3390/atmos14050768

Chicago/Turabian Style

Zhang, Xin, Lihua Zhou, Xingying Zhang, Yong Luo, and Lei Sun. 2023. "A Case Study on the Impact of East Asian Summer Monsoon on Surface O3 in China" Atmosphere 14, no. 5: 768. https://doi.org/10.3390/atmos14050768

APA Style

Zhang, X., Zhou, L., Zhang, X., Luo, Y., & Sun, L. (2023). A Case Study on the Impact of East Asian Summer Monsoon on Surface O3 in China. Atmosphere, 14(5), 768. https://doi.org/10.3390/atmos14050768

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop