Abstract

Sulfur dioxide (SO2) in the planetary boundary layer (PBL) as a kind of gaseous pollutant has a strong effect regarding atmospheric environment, air quality, and climate change. As one of the most polluted regions in China, air quality in Beijing-Tianjin-Hebei (BTH) region has attracted more attention. This paper aims to study the characteristics of SO2 distribution and variation over BTH. Spatial and temporal variations for a long term (2006–2017) over BTH derived from OMI PBL SO2 products were discussed. The temporal trends confirm that the SO2 loading falls from average 0.88 DU to 0.16 DU in the past 12 years. Two ascending fluctuations in 2007 and 2011 appeared to be closely related to the economic stimulus of each five-year plan (FYP). The spatial analysis indicates an imbalanced spatial distribution pattern, with higher SO2 level in the southern BTH and lower in the northern. This is a result of both natural and human factors. Meanwhile, the SO2 concentration demonstrates a decreasing trend with 14.92%, 28.57%, and 27.43% compared with 2006, during the events of 2008 Olympic Games, 2014 Asia-Pacific Economic Cooperation (APEC) summit, and 2015 Military Parade, respectively. The improvement indicates that the direct effect is attributed to a series of long-term and short-term control measures, which have been implemented by the government. The findings of this study are desirable to assist local policy makers in the BTH for drawing up control strategies regarding the mitigation of environmental pollution in the future.

1. Introduction

Sulfur dioxide (SO2) is a short-lived gas primarily produced by volcanoes, power plants, refineries, metal smelting, and burning of fossil fuels. When SO2 remains near the Earth’s surface, it is toxic, causes acid rain, and degrades air quality. It forms sulfate aerosols that can alter cloud reflectivity and precipitation in the free troposphere [1, 2]. As a kind of important atmospheric pollutants, SO2 critically affects the global environment,climate change, and public health. SO2 has become one of the popular research topics in the past decades, to examine its changes over some of the world’s most polluted regions [310].

China, with its incredible economic growth has been the focus of many studies during the previous decade because of its increasing sulfur dioxide emissions’ contribution to the Earth’s atmosphere [1116]. The sources of SO2 are both natural (volcanic) and anthropogenic emissions. Natural emissions include intentional biomass burnings and volcanic eruptions. Anthropogenic emissions are mainly due to fossil fuel burning (e.g., coal and oil), which accounts for more than 75% of global emissions [3]. Anthropogenic SO2 emissions are predominantly in or slightly above the planetary boundary layer (PBL), impacting on regional variations of aerosol types [17].

Satellite measurements of trace gases have been widely used and been an essential way to provide global, consistent observations for detecting, monitoring, and quantifying the SO2. The first space-based quantitative data on SO2 mass of the El Chichon volcanic eruption in 1983 were obtained from Total Ozone Mapping Spectrometer (TOMS) on board Nimbus 7 [18]. Subsequently, anthropogenic SO2 sources from power plants in eastern Europe [19, 20] and smelters in Peru and Russia [21] were demonstrated through detection of SO2 emissions using Global Ozone Monitoring Experiment (GOME) measurements on the Earth Research Satellite 2 (ERS-2). The tropospheric SO2 were detected by the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) on board the ENVISAT [22] and the Global Ozone Monitoring Experiment-2 (GOME-2) instrument on MetOp-A [23]. The Ozone Monitoring Instrument (OMI) on NASA’s Aura spacecraft enables to provide daily, nearly global measurements of ozone columns and aerosols, and the trace gases with the highest spatial resolution and the longest data record currently available [24, 25].

OMI data have been applied to assess the effect of pollutant transmission, analyze pollutant source contribution, evaluate pollutant emission inventory, observe regional pollution changes, and quantify the reduction of power plant emissions [14, 16, 17, 2628], due to the higher spatial and temporal resolution. Three major air pollutants (NO2, SO2, and CO) in China before, during, and after the Olympic Games from Aura’s Ozone Monitoring Instrument (OMI) and Terra’s Measurements of Pollutants in the Troposphere (MOPITT) instrument have been measured [29]. Annual emissions by sector and fuel types calculated from satellite data show an increasing trend of SO2 during 1996–2008 and decreasing thereafter in China [16]. Substantial changes in SO2 emissions in the northern China for the period 2005–2008 were analyzed [14]. The spatiotemporal variation of SO2 concentration during 2005–2008 over China from the planetary boundary layer (PBL) SO2 column concentration retrieved from OMI has been analyzed [28]. Long-term SO2 pollution changes over China or region have been observed through OMI observations [29]. A long-term trend of NO2 and SO2 levels (2005–2014) of the Henan province in China has been retrieved from the OMI [30]. In the past decades, China has adopted different policies for air quality control consistently, such as carbon reduction, energy saving, and other measures to reduce SO2 emissions [31, 32]. The latest findings represent that large reductions in SO2 are benefiting from the effective control policies in China [33, 34].

In this study, we analyze trend variation and distribution in SO2 concentrations over the Beijing-Tianjin-Hebei (BTH) region observed by the OMI between January 2006 and December 2017. Figure 1 shows the multiyear average spatial distribution map of SO2 based on OMI data (2006–2017) over China, which clearly shows the hotspots of SO2 in the North China Plain. Meanwhile, it is noted that the BTH region is one of the most polluted regions in China. Accordingly, tight emission control arrangement of this area always adopted, the SO2 emission of Hebei rank 3rd in 2013 went down to 5th with 17.4% rate of decline among all provinces in China, and Beijing and Tianjin declined 24.1% and 17.2%, respectively. Regional SO2 time evolution and spatial distribution are discussed in the following. The datasets and monitoring area are presented in Section 2. The analysis and associated findings are described in Section 3. Finally, the main conclusions are summarized in Section 4.

2. Method Description

2.1. Data Sources

Data used in this study over the BTH are based on OMI SO2 products with Dobson Units (DU, 1 DU = 2.69 × 1016 molecules/cm2). The Ozone Monitoring Instrument (OMI) is a sun-synchronous polar orbiting Dutch/Finnish sensor on the AURA satellite launched on 15 July 2004. The science goals of OMI are directly related to these questions and focus on (1) measuring the ozone layer and its destroying trace gases BrO and OClO, (2) tropospheric pollution by ozone, nitrogen dioxide, tropospheric aerosols, SO2, and formaldehyde, and (3) detection of species important for climate change such as aerosols, clouds, and ozone. The OMI measures the radiation backscattered by the Earth’s atmosphere and surface over the entire wavelength range from 270 to 500 nm, with a spectral resolution of about 0.5 nm, and high spatial resolution (13 × 24 km2), and daily global coverage [24]. OMI data have 4 processing grade products: Level-0, Level-1, Level-2, and Level-3. For this study, we used the PBL SO2 vertical column density from the Level-3 0.25 × 0.25 degree gridded OMI/AURA SO2 data product. The data used here with the time span of January 1, 2006 to December 31, 2017, were obtained from Giovanni interface (http://giovanni.gsfc.nasa.gov/giovanni/), derived from the NASA Goddard Earth Sciences Data Active Archive Center (GES DISC; http://disc.sci.gsfc.nasa.gov) [35, 36]. These Level-3 products have been widely used [15, 34, 37, 38]. The original OMI PBL SO2 product employed the band residual difference (BRD) algorithm [39]. But this product has a high noise level and systematic artifacts that required empirical corrections [5, 6]. A new operational OMI PBL SO2 product produced with the principal component analysis (PCA) algorithm was released [40, 41]. Validation of these two algorithmic products has been analyzed [42].

Daily satellite observations were retrieved with the given longitude and latitude for China and BTH region, as shown in Figure 2, to gain insight into the distribution of SO2 columns in the BTH region. The data were gridded onto monthly 0.25 × 0.25 fields and then onto seasonal and yearly maps. The meteorological conditions used for analysis are from NCEP reanalysis data (https://www.esrl.noaa.gov/psd/).

2.2. Study Area

Our study area focuses on the Beijing-Tianjin-Hebei (BTH) region (Beijing, Tianjin, and Hebei integration), which is the most polluted industrialized regions in China (Figure 1). The BTH region is located in the northwest part of the North China Plain as shown in Figure 2 (36°05′–42°37′N, 113°11′–119°45′E), with a total area of 216,000 km2 and more than 110 million residential population. The BTH region includes two municipalities (Beijing and Tianjin) and one province (Hebei) (Table 1), which contain thirteen cities: Beijing, Tianjin, Baoding, Langfang, Tangshan, Shijiazhuang, Xingtai, Handan, Cangzhou, Hengshui, Qinhuangdao, Chengde, and Zhangjiakou.

As one of the most economically vibrant regions in China, the BTH region covers only 2.3% of the Chinese territory but generates over 10% of the total national gross domestic product (GDP) in 2016 (National Bureau of Statistics of China (NBSC)) [43]. As the main high-tech and heavy industry base of China, there are mainly the automotive industry, electronic industry, machinery industry, iron industry, and steel industry. The SO2 map (Figure 2) shows hotspots associated with the major coal-fired power plants and industrial activities. Figure 3 reveals that high sulfur coal-fired power plants are the major contributors to the SO2 concentrations over the BTH region [44, 45]. OMI-derived spatial distribution shows generally good agreement with these main anthropogenic emission sources from burning sulfur-contaminated fossil fuels as Figure 3. The last decade has seen frequent occurrences of severe air pollution episodes (haze), and the high SO2 loading observed certainly contributed to PM 2.5 problems, especially in winter. Atmospheric environment quality has attracted more and more attention related to air pollution prevention and control policy-making.

3. Results and Discussion

The temporal and spatial variations of SO2 concentration in the region of Beijing-Tianjin-Hebei during the period of 2006–2017 were analyzed based on the satellite OMI data to characterize the variation of SO2 columns.

3.1. Overall Temporal Trend of SO2
3.1.1. SO2 Decadal Change

Figure 4(a) clearly shows the annual average SO2 concentration time series change trend over BTH from 2006 to 2017. The SO2 loading has decreased over the recent years without clear regularity, which is in line with the study results by others [41, 46]. The plot describes the irregular upward and downward rule. It is necessary to identify the specific period for the pollution attenuation, to a certain extent, related to governmental actions. It was found that SO2 peak in 2007 with an upward of 14.96% compared to 2006, reflecting the total SO2 emission in China had substantially increased from 2000 to 2007. The phenomenon has been largely driven by expansion in manufacturing industries and fueled by coal for the Chinese economic growth [14]. And then, an obvious downward trend appeared from 2007 to 2010, and the decline was 12.92% in 2007-2008 and 13.23% in 2008-2009, respectively. Afterwards from 2009 to 2010, the downward trend slowed down with the rate of decline 0.68%. The decrease was mainly due to China’s 11th FYP requiring power plants to install FGD devices.

Nonetheless, there was a temporary rebound from 2010 to 2011 with 19.19% growth rate. The brief period of emission growth can probably be attributed to the government stimulus for resurgence of economy in response to the global financial crisis of 2007-2008. Subsequently, there was a sharp decrease with 37.48% reduction over the 7-year period during 2011–2017. Through the seasonal variation tendency (Figures 4(b)4(e)), the increase in 2011 mostly comes from summer, accompanied by the industrial production slowdown in the latter half of 2011.

As an industrialized and populated region, the levels of BTH air pollution are determined by population density, economic activity, type of power generation and fuel used, and regulatory policies [47]. Here, we examine the SO2 changes with the regional SO2 emission data derived from national statistics (http://www.stats.gov.cn/tjsj/). It is interesting to find out the yearly emission trend that agrees well with the trend based on OMI retrieved SO2 as mentioned earlier (Figure 5), along with the estimated SO2 emissions data from power plants. Some mismatch is expected considering the different observation means and uncertainties with satellite and ground measurements.

Two fluctuations occur in the long trends; therefore, we broke down the data into two periods, 2006–2010 and 2011–2017, to verify the coherence with precursors’ investigation [30, 48] as to the changes in SO2 loadings. Figure 6 shows the 5-year mean SO2 concentration distribution maps over the BTH region of China’s 11th five-year plan (2006–2010) and second 7 years mean SO2 meets China’s 12th five-year plan (2011–2015) and China’s new 13th five-year plan (2016–2017), respectively. The dramatic decrease in SO2 loading (Figure 6) well illustrates the achievements and improvements due to a series of air pollution control policies. The average SO2 was 0.87 DU for 2006–2010 (Figure 6(a)). The average concentration decreased to 0.72 DU for 2011–2015 and 0.23 DU for 2016-2017.

According to the above statistics (Figure 5), it is more clearly shown that the SO2 loading over the BTH peaked in 2007, and then presented overall decreasing trend. This can be attributed to the emission control measures taken by the government. The sharp decline from 2007 to 2010 is closely related to the installation of flue gas desulfurization (FGD) and follow-up effects of strict pollution reduction measures implemented before the 2008 Beijing Olympic Games. The results also show that atmospheric SO2 loadings in BTH have drastically decreased by 17.24% (2011–2015) and 75.29% (2016-2017) related to the average quantity in 11th five-year plan (2006–2010), due to more stricter emission reduction targets, such as flue gas desulfurization control of enterprises, new energy to replace polluting energy, 50% of privately owned vehicles were banned through an odd and even number system, and switching from coal to natural gas for heating. And it is noted that the SO2 loadings have been with a short-lived upswing in the early of 11th and 12th five-year plan that may be caused by the government’s economy stimulus.

3.1.2. Seasonal Variations

A more detailed temporal variation of SO2 concentration during 2006–2017 in the BTH region is analyzed. As shown in Figure 7, the SO2 seasonal mean distribution has an obvious seasonal change sequence as winter (December, January, and February)>autumn (September, October, and November)>spring (March, April, and May)>summer (June, July, and August). The average SO2 loading of four seasons (spring, summer, autumn, and winter) was 0.59 DU, 0.52 DU, 0.73 DU, and 0.89 DU, respectively. Figure 8 also shows seasonal variation by calculating the monthly SO2 columns from 2006 to 2017. SO2 levels are peak in winter and minimum in summer and they are mainly attributed to the difference in pollution diffusion, which is caused by anthropogenic activities and meteorological conditions. Coal heating is a major source accounting for high SO2 levels in winter. The SO2 rapid drawdown in 2016 and 2017 represents the fact that coal heating had been replaced by natural gas is successful in reducing sulfur dioxide emissions. In addition, the SO2 in the PBL has short lifetimes during the warm season, on the short time scale of days to months, and meteorology plays an important role in regional air pollution [6]. The pollutants in the atmosphere have a dilution effect. The meteorological data of 2007, the highest loading over the past twelve years, were chosen as an example to interpret the seasonal differences (Figure 9). In summer, near-surface temperature is higher and air convection is stronger, besides relative humidity and amount of precipitation reaches the annual maximum. All these aforementioned basics can speed up the diffusion of atmospheric pollutants, which have given birth to short lifetimes of SO2 loading. On the contrary, the temperature, relative humidity, and precipitation in winter go against the transformation from SO2 to sulfate, making the SO2 with the longest lifetime in winter [48].

3.2. Spatial Distribution Pattern

In Figure 10, the spatial distribution of the average SO2 column in the Beijing-Tianjin-Hebei region from 2006 to 2017 is presented. As shown in Figure 10, the high SO2 concentration in BTH is significantly distributed in the southwest and eastern regions.

As for the thirteen cities, multiyear average SO2 clearly presents the spatial distribution discrepancy. We set up two borderlines with SO2 column amount (1.0 DU and 0.5 DU) to identify hotspots over BTH (Figure 11). The SO2 concentrations of Handan, Xingtai, and Shijiazhuang in the southwest of BTH are more than 1.0 DU. Cities with SO2 concentrations exceeding 0.5 DU include Tangshan, Tianjin, and Qinhuangdao located in the coastal beach areas of Bohai Bay, Hengshui, and Cangzhou adjacent to the Shandong province, as well as Baoding and Langfang around Beijing. The SO2 loading of Zhangjiakou and Chengde in the north BTH region cities are lower than the other eleven cities less than 0.5 DU.

The characteristic of spatial distribution over BTH can be interpreted from both natural factors and human activities. As shown in Figure 12, the terrain declines semicircularly from the northwest to the southeast over the BTH region due to the mountains and plains landforms. Figure 13 also shows the landscapes in this region, plateaus, mountains, and hills account for 54% of the entire BTH region centering on the area of northwest. Surface pressures gradually decrease from eastern to the northwest over the mountains (Figure 14). The central and southeast plains account for 46% of the region. Forests and grasslands are representative natural ecological landscapes in the region and account for 38% of the study area, and they are mainly distributed in the Yanshan Mountains, the Taihang Mountains, and the northwest edge of the Inner Mongolia Plateau [49]. We also found that relative humidity is different in this region due to the terrain and surface coverage (Figure 9). The whole area is mostly above 30%, especially the northwest is with higher values above 50% in winter, and wind also plays a major part in contributing to the SO2 spatial distribution. Wind direction and wind speeds determine the pollutant transmission path and diffusion velocity. The southern part of BTH underlies the leeward area which makes the SO2 difficult to diffusion and dilution, leading to the SO2 accumulation. However, the area of northern BTH located in the upwind accompanied with high wind speeds in favor of SO2 diffusion as shown in Figure 14. From the point of human factors, coal, petrification, motor vehicles, and iron and steel industrial emissions are the major source of the BTH region, which has a close relationship with the spatial distribution of population. The spatial pattern of population over BTH exhibits that the southern plain area is more likely to cause human activities variance than the northern region. The southern region is so flat and densely populated that human activities and industrial emissions greatly influence the environment.

Besides natural factors discussed above, anthropogenic factors also take up significant role, contributing to heavy pollution. The BTH region is adjacent to Inner Mongolia, Shanxi, Shandong, and Henan, which is the top five major pollution emission provinces in China (Figure 12), where gathering the most heavily emission sources in China is shown in Figure 3. As shown in Figure 15, we calculate the annual SO2 columns for each city of the BTH region. Handan, Xingtai, and Shijiazhuang as the emission hotspots will be served for further analysis. At first, the high SO2 loadings in the southwest of BTH around the province of Shanxi has nearly a hundred of coal-fired power plants that help contribute to the higher value. Secondly, the BTH region has the largest iron and steel industrial scale in China, among that Beijing, Tianjin, Shijiazhuang, Handan, and Tangshan occupies the majority. Meanwhile, these cities have been the most polluted in BTH region due to unfavourable factors, such as meteorological and geographic conditions and population density. Tangshan, the iron and steel producer of Hebei province, is the second most polluted city because of its high smoke SO2 emission. Although population and enterprises of Beijing are high, the highest green coverage rate and strict pollution control measures make the pollution level relatively low in the BTH region. Zhangjiakou and Chengde had a low SO2 column with the least amount of population density, industrial enterprises as well as the terrain advantage. The amplitude of each urban curves variation is able to reveal the degree of being influenced by natural and anthropogenic sources.

3.3. Typical Pollution Control Cases

As analyzed in Sections 3.1 and 3.2, there is an overall large decreasing trend of SO2 emission since 2011. The increasing/decreasing trend of SO2 in this period has been discussed with decadal, seasonal, and district variability. The results suggest that long-term variations are attributed to stricter emission reduction. In particularly, three typical episodes (the Olympic Games in 2008, the Asia-Pacific Economic Cooperation (APEC) in 2014, and Military Parade in 2015) have been taken place in Beijing during 2006–2017. A series of strengthened emission reduction measures have played an important role in preventing and controlling air pollution, which significantly improved the air quality in Beijing and neighboring regions [5054]. Therefore, we evaluate the variation (Figure 8) by tracking the changes of SO2 to reveal the effective operations driven by environmental policy and measures. The 29th Olympic Games were held in Beijing, started on 8 August 2008, the 2014 Asia-Pacific Economic Cooperation (APEC) Economic Leaders' Meeting was held on November 10 and 11 in Beijing, and China Military Parade was held on 3 September 2015. Prior to the events, emission reduction measures were began to conduct. Monthly, SO2 column concentration of other years from 2006 to 2017 corresponding to these periods have been compared with the basis of 2006.

As shown in Figure 16, the SO2 average concentration for the period of Olympic Games decreased significantly to 14.92% and 7.76% compared to the neighboring years; the average SO2 concentration in the APEC conference was significantly lower than that of the same periods in past ten years, and it declined 28.57% compared to the average of November in other years between 2006 and 2015; and the SO2 concentration in September of the 2015 Military Parade, with 27.43% decline, reached the lowest value compared to ten years before. Standard deviations can reflect that more clearly. Besides, the standard deviations verify the previous analysis of seasonal characteristic and brief ascent in 2011.

During the three events, air pollution control policies have been reinforced. To achieve the goal of “green Olympic Games” [55], China has released a series of air pollution control policies to improve air quality. The government has implemented a series of long-term pollution reduction measures, such as coal-fired power plant in Beijing to install the desulfurization equipment, part of closure of small power plants in this area near Beijing during the Olympic Games, and about 94% of the small coal-fired boilers to use clean energy transformation. The government also implemented some short-term strategies; for example, from July 1st to September 20th of 2008, the vehicles with exhaust emissions that failed to meet the European No. 1 standard were all-day forbidden on the roads; from July 20th to September 20th, the odd/even license plate number rule was applied on personal vehicles in Beijing [12]; power generation facilities were run only 30% of the equipment to stop all construction activities; some heavy-polluting factories were closed during the Olympic Games; and some heavily polluting companies around Beijing city were closed [29, 31]. During the APEC meeting and the Military Parade, the government also implemented a series of measures to ensure the safety and environmental protection measures, and these measures are similar even more and stricter than the 2008 Olympic Games, as with Beijing, the surrounding six provinces also have taken similar measures [27, 56]. During the period from August 20th to September 3rd of 2015, Beijing adopted to strengthen urban transportation management to strictly limit the motor vehicle population [57].

Relevant conclusion comes out that any air pollution activities by human beings in the entire BTH region should be controlled with the most rigorous management. Beijing as China’s political and administrative center has strict pollution control measures and high execution efficiency. Meanwhile, the environment policies to a certain extent affect the surrounding cities of Beijing. We focused on several affairs, such as the 2008 Olympic Games, 2014 APEC (Asia-Pacific Economic Cooperation), and 2015 Victory Day Military Parade, by means of analyzing the concentration of SO2 during the events to confirm the effectiveness of environmental management measures. The relative changes clearly demonstrate that these measures were effective in reducing SO2 concentration during the periods of those affairs taken place.

4. Conclusions

In this study, the past 12 years (2006–2017), OMI observations have yielded profound insights into the spatial distribution and temporal trends in SO2 emission over the Beijing-Tianjin-Hebei region. The SO2 loading distribution has close correlation with their emission sources. Spatiotemporal variation characteristics over BTH can more clearly reflect the natural and anthropogenic emission sources, to provide references to air pollution prevention and control. The main conclusions are as follows:(1)The temporal changes (2006–2017) over the BTH region exhibit the upward and downward trend consistent with the national trend in China. According to the dipping and heaving, we find that the gridded SO2 data can be divided into two phases: the first (2006–2010) and second (2011–2017) period, to observe changes in SO2 loadings presenting the Chinese government actions and policies and accomplishments in addressing air pollution. The SO2 loadings have drastically decreased by more than 30% from the 2006–2010 period to the 2011–2017 period. SO2 peaked in 2007 and the secondary peak was in 2011, to a certain extent, referring to the economic policy stimulus in the early of each FYP. Meanwhile, OMI observations show generally good agreement with independent SO2 emission.(2)The annual cycles of SO2 show a pronounced seasonal pattern, with the highest values occurring in winter and the lowest values in summer. This seasonal variation can be explained mainly by the seasonality of emission strengths, lifetimes of these pollutants, and meteorological factors.(3)Spatial distribution of BTH is also characterized, which is interesting to find that progressive changes from the northern to the southern regions are associated with the terrain and surface coverage, besides industrial pattern and population distribution. In the cities of BTH, we found that rapid SO2 reductions generally correlate well with sharp reductions in industrial activity.(4)China made great progress in pollution control with implementation of a series of major policies. The concentration of SO2 during 2008 Olympic Games, 2014 APEC, and 2015 Military Parade have been compared to identical months in other years, and the result shows that concentration of SO2 column was relatively low in these episodes because of reinforced emission controls.

These findings demonstrate that SO2 concentration over BTH do not follow simple linear trends, but instead reflect a repercussion of environmental measures and political economic activities. In terms of temporal changes, the decreasing trend has been observed since 2011, mainly due to the government efforts to restrain emissions from the power and industrial sectors. On the other hand, spatial characteristics are regarding with natural and anthropogenic factors. One of the important goals in this work is to evaluate the effectiveness of a series of policy to control the pollution problems, and three typical time quantum status reflects that the national environmental pollution control measures have taken an evident effect. These findings can provide a basis for the development of environmental management measures during the Winter Olympics in 2022, also for the national environmental pollution prevention and control as a reference. The results clearly illustrate effectiveness of central governmental policies regarding emission mitigation of SO2 and aid in important policy implications for the future reduction action plans to provide better air quality of China.

Data Availability

The OMI SO2 data are downloaded from http://giovanni.gsfc.nasa.gov/giovanni/. The meteorological data are from https://www.esrl.noaa.gov/psd/. The SO2 emission data have been publicly released from http://www.stats.gov.cn/tjsj/. The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (Grant no. 41501404) and Natural Science Foundation of Hebei Province under Project no. F2017203220.