Abstract

Compared with other regions in China, air pollution on the North China Plain (NCP) is serious. Fine particle pollution has been studied in-depth, but there is less research on long-term troposphere ozone (O3) variation. This study focuses on the summer interannual tropospheric O3 variation on the NCP and its influential factors. Our analysis relies on satellite observations (O3, nitrogen dioxide (NO2), sulfur dioxide (SO2), carbon monoxide (CO), and formaldehyde (HCHO), determined as vertical column density of the troposphere) and dynamical processes (El Niño-Southern Oscillation (ENSO), potential vorticity (PV), the quasibiennial oscillation (QBO), and East Asian summer monsoon index (EASMI)). Our results show the vertical column density of tropospheric O3 has a transition from the increasing trend to decreasing trend during the summer of 2005–2016. The summer series of tropospheric O3 show two distinct phases: the first phase (2005–2011), with an average growth rate of 0.55 ± 0.20 DU/yr, and a second phase (2012–2016), with an average reduction rate of 0.16 ± 0.23 DU/yr. The tropospheric NO2 column in the NCP also has a transition from the increasing trend to decreasing trend during the summer of 2005–2016. Tropospheric NO2 and CO column concentrations obtained from satellite observations indicate that emission reductions might be the main cause of the tropospheric O3 decrease. Particularly, the reduction of nitrogen oxides (NOx) is more significant, and NO2 decreased by (0.45 ± 0.11) × 1015 molec·cm−2 per year in summer since 2012. However, tropospheric column HCHO shows an increase of 0.05 × 1015 molec·cm−2 per year during the whole period of 2005 to 2016. An O3-NOx-VOC sensitivity experiment in the NCP showed that the O3 is still in a NOx-saturated state in some heavily polluted cities, although the NOx emissions are decreasing overall. In addition to the chemical reactions, atmospheric dynamic processes also have an effect on tropospheric O3. Finally, we built a model to analyze the contributions of chemical processes and dynamic processes to the tropospheric O3 column in the NCP. For the chemical process variables, 69.73% of the observed trend of tropospheric O3 could be explained by the NO2 tropospheric column. Therefore, the reduction of tropospheric O3 since 2012 is associated with the reduction of NOx. For the dynamical process variables, ENSO, PV, and EASMI can explain 60.64% of the observed trend of tropospheric O3. This result indicates that the atmospheric circulation of the western Pacific Ocean in summer has a significant impact on the interannual trends of tropospheric O3 in the NCP. It is also found that chemical processes had a more important impact on interannual tropospheric O3 than dynamic processes, although the dynamic processes cannot be neglected.

1. Introduction

Atmospheric ozone (O3) is mostly found within the stratosphere, and tropospheric O3 is approximately one-tenth of the atmospheric column O3 [1]. However, tropospheric O3 has direct and detrimental impacts on ecosystems [2, 3]. The production of tropospheric O3 is chemically controlled by the nonlinear relationship of its precursors, volatile organic compounds (VOCs), and nitrogen oxides (NOx), implying the complexity of the O3 pollution. The temporal and spatial changes of tropospheric O3 are also controlled by meteorological conditions [1].

Urban agglomerations in China have been experiencing O3 pollution in recent years. The Yangtze River Delta (YRD) is one of the regions experiencing serious O3 pollution, with the highest frequency occurring in late spring and early summer [4]. The Pearl River Delta (PRD) is another region with serious O3 pollution [5]. The NCP has not only been suffering from severe hazy weather but is also one of the regions with serious O3 pollution in summer [68]. It was reported that O3 pollution episodes reached 286 ppbv (1-h O3 = 286 ppbv) observed on June 30, 2005, in mountainous areas of north Beijing [9]. In the PRD, surface O3 increased by 0.86 ppbv/year from 2006 to 2011 [10]. In the NCP, aircraft data indicated a boundary-layer O3 increase of 2%/year in the summer time during 1995–2005; the surface daily 1-hour maximum O3 in urban Beijing increased by 1.3%/year during 2001–2006 [11], and the daily 8-hour maximum O3 in rural Beijing (Shangdianzi) increased at a rate of 1.1 ppbv/year during 2003–2015 [12]. However, due to the environmental protection regulations in China, emissions of precursor NO have decreased since 2011 and 2012. For 2010 and 2014, NO emissions were 1.6 and 1.5 Gg/d in the PRD, 3.9 and 3.0 Gg/d in the YRD, and 15.6 and 14.3 Gg/d in the NCP, respectively. Model and OMI HCHO data show the uptrends in East Asia resulting from anthropogenic activities in China during the last decade and future [13, 14]; however, the trends are negative in the PRD. Much work has been done via case studies of the O3–VOC–NOx system sensitivity [15]. Areas around the Bohai Sea became more NO saturated [16]. However, the long-term trend of O3 is less noticed and studied. Most of the research on O3 in the NCP was based on model simulations and short-term site observations [1719] and lacks long-term sequence presentation. Research on long-term changes in O3 pollution is very limited due to scarce ground observations.

In this study, O3 long-term variations are investigated based on atmospheric compositions obtained from satellite observations. In recent years, satellite data have been used to study air pollutants [20, 21]. Atmospheric environmental satellite loads have nadir and limb scan modes. Limb mode instruments provide vertical column density and vertical profile data. The Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY) has three different viewing geometries: nadir, limb, and sun/moon occultations. Tropospheric Emission Spectrometer (TES) operates in a combination of limb and nadir mode. The Ozone Monitoring Instrument (OMI), Measurement of Pollution in the Troposphere (MOPITT), and Total Ozone Monitoring Spectrometer (TOMS) are nadir instruments and provide total vertical column (O3, SO2, NO2, HCHO, CO, and CH4). These data have been used to study air pollution [22], and a little research on greenhouse gases has also been involved [23]. Satellite data of column density for SO2, NO2, and CO are often used to study air pollution directly. Many studies indicated that OMI observations are credible for assessments of regional and global characterizations of NO2 and SO2 spatiotemporal variability [2426]. However, due to the particular characteristics of the vertical distribution of O3 (the peak in the stratosphere), it is not appropriate to use the total amount of the nadir column data alone. The nadir observation instrument OMI can obtain the O3 vertical column concentration of the whole atmosphere, and the limb observation instrument Microwave Limb Sounder (MLS) can obtain the vertical distribution information of O3. The tropospheric vertical column O3 can be obtained by combining the two instruments, which is the OMI/MLS satellite tropospheric O3 [27, 28]. The temporal and spatial distributions of trace gases are also discussed based on satellite.

Ozone-NOx-VOC sensitivity and explicative variables of chemical and dynamic processes are discussed in this study. The atmospheric circulation system parameters affecting the weather in China are used to measure the impact of meteorological conditions, such as ENSO, PV, QBO, and EASMI.

Descriptions of satellite observational data, meteorological data, and methods related are presented in Section 2. Section 3 presents the evaluation of the OMI/MLS tropospheric O3 interannual summer variation and trends of tropospheric O3 and trends of relevant trace gases and ozone-NOx-VOC sensitivity. Especially, Section 3.5 shows the summer O3 fraction of explained variation with respect to the two issues listed above (emissions and meteorological conditions). Finally, the conclusions are summarized in Section 4.

2. Materials and Methodology

2.1. Satellite Data and Methods

Tropospheric O3 data are obtained from combined observations of two satellite instruments, the Ozone Monitoring Instrument (OMI) and Microwave Limb Sounder (MLS). Daily OMI/MLS tropospheric O3 data were obtained by subtracting MLS stratospheric column O3 from OMI total column O3. Stratospheric column O3 from MLS was spatially interpolated each day to infill the actual along-track measurements. The monthly means were then determined by averaging all available daily data within each month. The spatial resolution of the data files is 1° × 1.25° (latitude × longitude) (https://acdext.gsfc.nasa.gov/Data_services/cloud_slice/new_data.html) [29, 30].

Daily SO2 and NO2 and monthly HCHO data are from the OMI aboard the EOS Aura spacecraft, which was launched on July 15, 2004. The NO2 tropospheric vertical column densities (cloud-screened at 30%) are from the OMI Level 3 daily global products with a spatial resolution of 0.25° × 0.25°. The NO2 slant column is retrieved in the spectral range of 405–465 nm with a fitting error of 0.3–1 × 1015 cm−2 before the row anomaly. The noise plumes from strong anthropogenic sources of SO2 (such as smelters and coal burning power plants) and from strong regional pollution can be detected in scene data. The noise standard deviation (sigma) of vertical column SO2 data products is about 0.4 DU (1 DU = 2.69 × 1016 molec·cm−2) and decreases with an increasing average number of scenes. The monthly values of NO2 and SO2 are calculated from the daily average. Detailed data descriptions are provided in the references [3134].

The daily total column CO is retrieved from the Measurement of Pollution in the Troposphere (MOPITT) available at https://terra.nasa.gov/about/terra-instruments/mopitt [35]. Its resolution is 1° × 1°. In our analysis, the missing values are eliminated. MOPITT’s spatial resolution is 22 km at nadir, and its swaths are 640 km wide. Moreover, it can measure the concentrations of CO in 5 km layers down a vertical column of atmosphere. The satellite is deployed in a polar Sun‐synchronous orbit with a 10 : 30 local equator crossover time. Retrieved CO total column day data were applied in the study. The MOPITT observations reported here are version 3 data, which have a 10% precision on both column and mixing ratio.

To find the inflection point of the time series, we perform quadratic (second-order polynomial) nonlinear regression. The normal calculation formula is as follows:where a1 and a2 are fitting coefficients, a3 is the intercept, is the time index, is the error term, i is the time index, and are the annual pollutant variables (NO2, HCHO, SO2, CO) in summer.

Based on long-term observation data, we first examined long-term interannual variations of O3 concentrations in the NCP during 2005–2016. We use the linear regression equation to calculate the growth trends of O3 and pollution gas emissions and divide time into two phases (according to O3 trend variation). The results are shown in Table 1. The linear regression equation iswhere a is the slope from time-series calculation, b is the corresponding intercept, and is the error term.

2.2. Dynamical Variables and Regression Model

Figure 1 shows the study area, which is the NCP. Summer (June, July, and August, i.e., JJA) was selected as the study period.

The considered NCP domain ranges between 34.5°N and 42.5°N in latitude and between 110.5°E and 120.5°E in longitude (Figure 1). The BTT domain ranges between 38.5°N and 40.5°N in latitude and from 115.5°E to 119.0°E in longitude (Figure 1). The other representative cities in the NCP include Shijiazhuang, Jinan, Taiyuan, and Chengde, the city center of which is located within the domain. Annual summer concentrations are obtained based on the regional averages of monthly concentrations calculated from the daily averages. Then, we have standardized (i.e., normalized) NO2 and HCHO concentrations and meteorological parameters.

A multiple linear regression model was established to assess the contribution of each parameter, as described in the literature [36]:where O3 is the summer monthly mean troposphere O3, is the slope calculated based on regional time series, i is the time index, is the regularized time series parameter, is the corresponding fitting coefficient, and b is the intercept. is the error term. is R squared goodness of fit. is adjusted R square. n is the number of samples. p is the number of parameters. Fraction of explained variation is the adjusted R squared of the multivariate regression model (explained variance/total variance). When a meaningful variable is added, adjusted R square will increase. We consider a 95% confidence range ( is considered to be significant). The standardization method used in this study involved normalizing to (0, 100) [37]. The standardized parameters are calculated by the following formula:where is the normalized time series parameter and and are the minimum and maximum values in the series, respectively [38].

Some large-scale dynamics influence the tropospheric O3 such as the quasibiennial oscillation (QBO) and the El Niño-Southern Oscillation (ENSO) related to the dynamical processes leading to a modulation of the stratospheric circulation and of the stratospheric-tropospheric exchanges (STE). EASM affects the horizontal transmission of tropospheric O3. To study the impacts of atmospheric dynamics on tropospheric O3, several atmospheric circulation factors have also been taken into account, as described below.

The QBO data at 10 hPa (QBO10) and 30 hPa (QBO30) were taken from http://www.geo.fu-berlin.de/met/ag/strat/produkte/qbo/singapore.dat. The QBO dataset adopted has been produced since 1987 from the Singapore data (monthly mean zonal wind components at Singapore, 1N/104E) by using the daily vertical wind profiles to obtain a higher vertical resolution. The time period used in this article is from the summer of 2005 to 2016 (June to August), and the other parameters are in the same period. The multivariate ENSO index (MEI) was taken from https://www.esrl.noaa.gov/psd/enso/mei/table.html. The bimonthly multivariate El Niño/Southern Oscillation (ENSO) index is the time series of the leading combined Empirical Orthogonal Function (EOF) of five different variables (sea level pressure (SLP), sea surface temperature (SST), zonal and meridional components of the surface wind, and outgoing longwave radiation (OLR)) over the tropical Pacific basin (30°S-30°N and 100°E–70°W). The Nino index was from https://http://www.esrl.noaa.gov/psd/gcos_wgsp/Timeseries/Nino34/. The East Asian summer monsoon (EASM) index is defined as an area-averaged seasonally (JJA) dynamical normalized seasonality at 850 hPa within the East Asian monsoon domain (10°−40°Ν, 110°−140°Ε) [3941]. The data were taken from http://ljp.gcess.cn/dct/page/65577. The QBO, ENSO, Nino indices, and EASM are all monthly data. The tropopause height is characterized by the potential vorticity at 300 hPa. These parameters were taken from the ERA-Interim reanalysis at http://apps.ecmwf.int/datasets/data/interim-full-daily/.

2.3. Ozone-NOx-VOC Sensitivity

The formaldehyde nitrogen ratio (FNR) has been used to indicate spatially and temporally heterogeneous ozone-NOx-VOC sensitivity [42]. FNR is defined as the ratio of tropospheric HCHO column concentration to the NO2 column concentration and is used here to study the ozone-NOx-VOC sensitivity in the NCP:

The HCHO and NO2 column concentrations retrieved from the Ozone Monitoring Instrument (OMI) aboard the NASA Aura satellite were used to calculate FNR during the O3 season. Based on previous experience, the O3 production regime is thought as VOC-limited if FNR < 1.0, NOx-limited if FNR > 2.0, and transitional (either NOx changes or VOC changes would be expected to change O3 in this interval) if 1.0 < FNR < 2.0 [43].

2.4. Evaluation Data

The World Ozone and Ultraviolet Data Center (WOUDC) ozonesonde provides an independent approach of ground observations of tropospheric O3 and is used here to evaluate the OMI/MLS tropospheric O3. Another independent tropospheric O3 observation is balloon radiosonde via an electrochemical concentration cell (ECC-Ozone Sensor) in Beijing (39.80°N, 116.47°E) during April to May 2005, and this approach is also used in this paper for evaluation of the OMI/MLS tropospheric O3.

3. Results and Discussion

3.1. Evaluation of the OMI/MLS Tropospheric O3

Column O3 (CO, in Dobson Units, DU) was determined by standard log-pressure integration of O3 partial pressure: , where is O3 volume mixing ratio in unit parts per million by volume (ppmv) and P is pressure in units hPa. The satellite and sondes are for a tropospheric average. Tropospheric O3 is the entire column in the troposphere [29, 30]. The validation is from OMI/MLS and ozonesonde monthly data. The correlation coefficient R between OMI/MLS tropospheric O3 and WOUDC ozonesonde tropospheric O3 is 0.936 for ozonesonde station locations lying between 25°S and 50°N, and the Root Mean Square Error (RMSE) is 3.18 ppbv. The deviation is smaller in the lower latitudes. This result also indicates closely similar signatures for seasonal cycles and spatial variability from the comparisons of OMI/MLS tropospheric O3 between the climatology and other data products [44]. Comparisons with Beijing (39.80°N, 116.47°E) ozonesonde data during April to May 2005 show a typical error of ±10%. If there are enough samples, it will be more convincing.

3.2. Interannual Summer Variation and Trends of Tropospheric O3

Figure 2 shows the monthly tropospheric O3 distributions for the study area from 2005 to 2016. The tropospheric O3 levels are lowest in winter (December/January) and highest in summer (June to August). Spring and autumn are transitional periods. The peak monthly O3 occurred in June/July. This pattern is consistent with the seasonal variations of temperature and solar radiation. The tropospheric O3 levels in July increased from 2005 to 2016 at an average rate of 0.2 DU per year, while no such apparent trend is found in winter. The lowest O3 levels generally occur in January. For instance, the lowest O3 level during this period occurred in January 2006, which accounted for only 5.65% of the accumulated monthly O3 of 2006. Similar low values are often recorded in the period of December to February; for example, the O3 levels in February 2012 and December 2015 accounted for 5.50% and 5.87% of the accumulated monthly O3, respectively. In contrast, the O3 level in July 2016 accounted for 11.91% of the O3 levels for the whole year. Because of the increase in O3 pollution incidents in North China that occurs during most summers, we focus on the situation of the summer period. The summer time series (Figure 2) have shown that the trend of tropospheric O3 is increasing with an average rate of 0.36 ± 0.09 DU/yr over 2005–2016. However, the summer tropospheric O3 column in the NCP has a transition from the increasing trend to decreasing trend over 2005–2016. The tropospheric O3 trend shows two distinct phases: the first phase (2005–2011), with an average increase rate of 0.55 ± 0.20 DU/yr, and the second phase (2012–2016), with an average decrease rate of 0.16 ± 0.23 DU/yr. The decrease in the tropospheric O3 in the latter phase might be caused by the emission reductions in recent years. To investigate the impacts of trace gases on tropospheric O3, tropospheric column concentrations of several trace gases and meteorological factors are discussed in Section 3.2.

3.3. Interannual Summer Variation and Trends of Relevant Trace Gases

The atmospheric nitrogen oxides (NOx, including NO and NO2) are important air pollutants and among the O3 precursors [45, 46]. Previous studies showed that there is generally a positive correlation between atmospheric NOx and O3 concentrations [47].

According to satellite observations, NO2 pollution in North China is more serious than in adjacent areas [48]. This is generally attributed to the northern industrial areas. Annual averages of tropospheric NO2 for summer months are shown in Figure 3. Over these 12 years, tropospheric NO2 peaked in 2011 and 2012, reaching 6.97 × 1015 molec·cm−2 in June, 6.15 × 1015 molec·cm−2 in July, and 6.14 × 1015 molec·cm−2 in August. NO2 concentrations showed a downward trend from 2011/2012 to 2016. Two phases are discerned: the first phase (2005–2011), with an average increase rate of (0.25 ± 0.05) × 1015 molec·cm−2·yr−1, and the second phase (2012–2016), with an average decrease rate of (−0.45 ± 0.11) × 1015 molec·cm−2·yr−1. For the whole period (2005–2016), the NO2 concentrations decreased at (−0.02 ± 0.05) × 1015 molec·cm−2·yr−1. In terms of spatial distribution, the NO2 concentration in urban areas is higher than that in suburbs. Especially in the Beijing-Tianjin-Tangshan area, the NO2 concentration was extremely high over this period. This result suggests that the BTT area might probably have been in the NO-saturated regime, which is consistent with Souri et al. [16]. Further ozone-NOx-VOC sensitivity will be discussed in Section 3.3.

HCHO is produced directly and indirectly from both biogenic and anthropogenic sources. According to the research by Zhu et al. [33], satellite observations of HCHO columns provide emissions of highly reactive volatile organic compounds (HRVOCs). This approach has been used previously in the US to estimate isoprene emissions from vegetation. According to Souri et al., in East Asia, OMI HCHO levels were substantially higher in urban regions than in forested areas, indicating predominant anthropogenic VOC emissions. Its spatial distribution map also proves the pattern (Figure 4). The trends of total HCHO columns from the OMI showed widespread upward trends in eastern China and the majority of enhancements in the northern regions from the HCHO time series in NCP from 2005 to 2016. From the HCHO time series from the OMI, the annual minimum is in winter and the annual maximum is in summer.

HCHO is an important intermediate for VOCs, and VOCs are among the O3 precursors. Thus, the variability of HCHO could be indicative of the spatiotemporal distributions of reactive VOCs. During the first phase (2005–2011), HCHO in the NCP increased at an average rate of (0.24 ± 0.11) × 1015 molec·cm−2·yr−1, and during the second phase (2012–2016), HCHO decreased at an average rate of (0.19 ± 0.32) × 1015 molec·cm−2·yr−1. For the whole period (2005–2016), HCHO increased by (0.05 ± 0.07) × 1015 molec·cm−2·yr−1.

SO2 does not typically have a significant impact on O3, but it can characterize pollution emissions level [49]. The annual trend of summer-time SO2 in North China is shown in Figure 5(c). Quantitative calculations indicate a decreasing overall trend from 2005 to 2016. In the three months (June–August), the monthly average content of SO2 decreased by 49.1%, 43.2%, and 43.6%, respectively. For comparison purposes, the analysis of SO2 is also divided into two phases: the first phase (2005–2011), with an average rate of increase of 0.01 ± 0.02 DU·yr−1, and the second phase (2012–2016), with an average rate of decrease of 0.05 DU·yr−1. For the whole period (2005–2016), SO2 decreased by 0.03 ± 0.01 DU·per·year. The change in SO2 in summer proves that the emission reduction is significant once again.

The satellite observations of CO indicate that the spatial distribution of CO is basically similar to those of NO2 and SO2, all of which have elevated levels in North China compared with surrounding areas. Considering larger scales, their distribution is similar, while for smaller scales, the distribution is different because the life cycle of NO2 and SO2 is indeed much shorter than CO. However, the vertical column densities of CO show a clear interannual downward trend (Figure 5(d)) [50, 51], which is different from the trends of NO2 and SO2. During the first phase (2005–2011), CO decreased at an average rate of (0.03 ± 0.03) × 1018 molec·cm−2·yr−1, and during the second phase (2012–2016), it decreased faster, with an average rate of (0.05 ± 0.05) × 1018 molec·cm−2·yr−1. For the whole period (2005–2016), CO decreased at an average rate of (0.02 ± 0.01) × 1018 molec·cm−2·yr−1.

The sink of CO includes diffusion, transport, deposition, and chemical transformation. The diffusion and transport of CO in the atmosphere are controlled by wind speed [52]. Local and short-term O3 concentrations are indeed most severely affected by VOCs. The study region of this paper is the entire North China region, and CO represents a level of primary emissions to a certain extent. Therefore, the CO variable, due to its long lifetime, is considered to be a proxy for large-scale emission changes that may regionally affect O3 [53, 54].

3.4. FNR and Ozone-NOx-VOC Sensitivity

From the perspective of photochemical reactions, tropospheric NOx and VOCs are precursors that need be considered to elucidate the sensitivity of O3 production. OMI observations of NO2 and HCHO column densities were used for the exploration of the sensitivities of surface O3 pollution to NOx and VOC emissions over the NCP. However, in general, the vertically integrated column may not represent the near-surface environment. This study attempts to use FNR to assess the sensitivities of tropospheric O3 to NOx and VOC emissions over the NCP. NOx-limited and NOx-saturated conditions were discerned by FNR thresholds at 90% confidence [55].

FNR may vary in space and time, suggesting that there is a transition between NOx-limited and NOx-saturated conditions in space and time. Figure 6 shows the average summer (JJA) FNR in the NCP during 2005–2016. A decrease in average summer FNR is observed for the first phase (2005–2011), and this decrease may be caused by increased NOx emissions (with the NO2 increase rate reaching 0.25 molec·m−2·yr−1 (Table 1)). In the urban areas of major cities in the NCP (e.g., BTT, Taiyuan, Shijiazhuang, and Jinan) the summer FNR was extremely low, and this result suggests that in urban areas, O3 was probably in the NOx-saturated regime. However, during the second phase (2012–2016), the ratio increased significantly as a result of NOx emission reduction (with a NO2 decrease rate of −0.45 molec·m−2·yr−1 (Table 1)). In representative cities of the NCP, the summer O3 experienced a transition from the VOC-limited regime to the NOx-limited regime. Away from these urban areas, the summer O3 had been NOx-limited.

The regional O3-NOx-VOC sensitivities in summer are illustrated by O3-NO2-HCHO sensitivities in Figure 7. Generally, higher HCHO would result in higher O3, while higher NO2 would or would not, depending on FNR (Figure 7(a)). When the value of FNR is equal to 2, the tropospheric O3 is close to the vicinity minimum. Therefore, emission reductions of VOCs and NOx are effective for controlling O3. For the situation with FNR greater than 2 (the area on the right side of line FNR = 2), emission reduction of VOCs is effective for O3 control, and reduction of NOx may even increase the O3 concentration. This finding suggests that O3 is in a NOx -limited regime. For the situation with FNR less than 1 (the area on the left side of line FNR = 1), emission increase of NOx would not effectively increase O3 concentration, while reduction of VOCs would, thus implying that NOx is oversaturated. For the situation with FNR between 1 and 2 (the area between line FNR = 1 and line FNR = 2), emission reductions of both VOCs and NOx are effective for O3 control, suggesting that O3 is in the transition regime. As seen from Figure 6, the FNR was greater than 2 over most areas of the NCP, but in the major metropolitan cities, the FNR was less than 2 and might even have been less than 1 for some years. In these cities, NOx was still saturated for tropospheric O3 production in the second phase, although their NOx emissions have been decreasing (Figures 3 and 7(b)). Because of the high urban NOx base level, the emission reduction of NOx did not result in a distinct O3 decrease. Overall, the NCP has been experiencing the transition from the NOx-saturated state to NOx-limited state by comparison of two time periods (2005–2011 and 2012–2016).

Due to the limitations of the time resolution of satellite observations, the sensitivity experiments we performed did not reflect the continuous variation in a day and daily changes. Because the daily variation of O3 is relatively large, OMI scanning transit time over NCP is approximately local time 13 : 45, only with once a day [30]. This sensitivity analysis can only reflect the situation of satellite transit time, which is the limitation of this experiment. Another point to note is that tropospheric O3 is different from the concerned surface O3 (e.g., 8-hour max) because it directly affects human health.

3.5. Summer O3 Fraction of Explained Variation

The fraction of explained variation is from adjusted R square in the multiple regression model formula (3). It is a change caused by other variables. The multivariate linear regression model described above was used to calculate the fraction of explained variation of the trend observed by OMI/MLS in the tropospheric O3 for the summer season (JJA) of the entire period from 2005 to 2016. We applied the model to the variables in sequence to identify the variables that were significant for tropospheric O3.

During the entire period of 2005–2016, the significant variables included NO2, HCHO, CO, the QBO, the potential vorticity (PV) at 300 hPa, the ENSO index, and EASMI (Figure 8). We divided these variables into two types: atmospheric chemistry factors (NO2, HCHO, and CO) and atmospheric dynamics factors (QBO, PV at 300 hPa, the ENSO index, and EASMI). The regression model was applied again to the atmospheric chemistry factors, and the significance of their regression fittings was in the order of NO2 > CO > HCHO. The first significant variable we found for this period was the monthly NO2. The NO2 variable, the stable state of NOx, is considered to be a proxy for large-scale emission changes that may affect O3 regionally, and NO2 could explain 69.73% of the initial trend (Table 2). CO represents a regional change in emissions due to its longer life cycle to some extent. CO could explain 60.14% of the initial trend. HCHO is an intermediate product of VOCs, and they are closely related. HCHO could explain 54.73% of the initial trend. Their overall contribution was 63.50%.

Contribution to O3 is defined by the adjusted R square of the multivariate regression model. SSE is the sum of squared residuals (SSE). Because as the number of variables increases, R square will inevitably increase, and the accuracy cannot be truly quantified. Looking at R square alone does not infer whether the added features make sense, so here we applied adjusted R square instead of R square. When a variable is added, if this variable is meaningful, adjusted R square will increase. If this feature is redundant, adjusted R square will decrease, so the combined contribution of NO2, CO, and HCHO to O3 is lower than the contribution of NO2 alone, which indicates the change of NO2 plays a leading role in the variation of O3, and the information of the other two substances is covered by NO2. This analysis is also applied to atmospheric circulation.

The effects of atmospheric circulation on tropospheric O3 are characterized by the atmospheric dynamic factors: the QBO, the potential vorticity (PV) at 300 hPa, the ENSO index, and EASMI. Similar to the atmospheric chemistry factors, the regression model was also applied to the dynamic factors in sequence for another residual trend. The significance of the regression fittings of the dynamics factors was in the order of ENSO > PV > EASMI > QBO. The ENSO index showed high significance () for this period, which could explain 57.69% of the initial trend. The normalized ENSO index experienced an increase over the second phase (2012–2016) with strong El Niño events in 2015–2016 (Figure 8). PV reflects the height of the troposphere, and it could explain 55.79% of the initial trend. EASMI represents the intensity of the airflow between the continent and ocean. The effect of horizontal airflow could explain 54.72% of the initial trend. QBO reflects the contribution of O3 exchange between the troposphere and the stratosphere; the vertical transport can explain 54.61% of the initial trend. Their overall contribution to the observed trend was 55.32%. The above results show that it is meaningful to add EASMI, and ENSO covers the PV information.

4. Conclusions

Satellite observations are applied to study the variability of O3 in the troposphere in the NCP. The data used in this study combine the advantages of satellite nadir observation and limb observation modes; the nadir mode facilitates observation of the total column, and the limb mode facilitates observation of the vertical profile. Evaluation of the MLS/OMI observations has shown that the tropospheric O3 trends seem reliable fairly and can be used to study the tropospheric O3 trend over the NCP.

Based on long-term statistics, the tropospheric O3 in the NCP is highest in summer (from June to August). In contrast, the O3 content is lowest during the winter period (from December to February). Summer tropospheric O3 column concentrations increased at an average rate of 0.55 DU/year in the NCP from 2005 to 2016.

Spatiotemporal variations of several key atmospheric components are analyzed to evaluate their impacts on O3 formation in summer. Major air pollutants showed a trend from increasing in the first phase (2005–2011) to decreasing in the second phase (2012–2016) because of the implementation of emission reduction measures. Due to the reduction of nitrogen oxides, the NCP has been experiencing the transition of O3-NOx-VOC sensitivity from the NOx-saturated state to the NOx-limited state. However, NOx is still saturated in megacities.

A multivariate linear regression model is applied to identify the processes driving the observed trend of summer tropospheric O3. The contribution of atmospheric chemistry is 63.50%. The first significant variable is NO2, and it can explain 69.73% of the initial trend. The ENSO also shows the high significance for this period, and it can explain 57.69% of the initial trend. The contribution of large-scale dynamical processes to the observed trend is 55.32%. The results suggest that both large-scale dynamical processes and regional emissions affect the interannual tropospheric O3 variation in summer.

Our results indicate that chemical processes had a more important impact on interannual tropospheric O3 than dynamic processes. This is also reflected by other research results. The summer tropospheric O3 presented a growth trend with the growth rate of 1.28 DU·per decade over the North China Plain (NCP) for the period of 1979–2013 mainly due to increased NOx and VOC emissions [56]. Significant increasing trends of tropospheric O3 are also found for all seasons except for winter, with a maximum rate of 1.10 DU·per decade for summer during 1979–2005. They do not consider atmospheric circulations to be a major factor in the tropospheric O3 increasing trend over NCP [57].

Remote sensing data could provide long-term characteristics. It is difficult for the atmospheric chemical transport model due to emission uncertainty and limitation of computing resource [58]. However, the statistical relationship is not completely deterministic, although the results selected have passed the significance test, so we regard our conclusion as a hypothesis, which needs to be validated by a 3D model with atmospheric chemistry and meteorology that can provide mechanism analysis. The conclusion based on remote sensing data in this study is very macroscopic, and some details require further separate research. In addition, specific vertical distribution, convection processes, long-range transport, and chemical nonlinear processes require an atmospheric chemical transport model for sensitivity experiments to quantify its effects [59].

Data Availability

Data sources and links are listed in the data introduction section.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors are grateful to all the organizations and individuals who provided the data. Lihua Zhou acknowledges the financial support from the China Scholarship Council. This research was funded by the National Key R&D Program of China (grant no. 2017YFA0603603), National Science Foundation of China (grant no. 41575144), the Fundamental Research Funds for the Central Universities (grant no. 312231103), and Natural Science Foundation of China and Hong Kong Federation of Trade Unions (NSFC-RGC) (grant no. 403331160396).