Effects of meteorology and emissions on urban air quality : a quantitative statistical approach to long-term records ( 1999 – 2016 ) in Seoul , South Korea

Together with emissions of air pollutants and precursors, meteorological conditions play important roles in local air quality through accumulation or ventilation, regional transport, and atmospheric chemistry. In this study, we extensively investigated multi-timescale meteorological effects on the urban air pollution using the long-term measurements data of PM10, SO2, NO2, CO, and O3 and meteorological variables over the period of 1999–2016 in Seoul, South Korea. The long-term air quality data were decomposed into trend-free short-term components and long-term trends by the Kolmogorov–Zurbenko filter, and the effects of meteorology and emissions were quantitatively isolated using a multiple linear regression with meteorological variables. In terms of short-term variability, intercorrelations among the pollutants and meteorological variables and composite analysis of synoptic meteorological fields exhibited that the warm and stagnant conditions in the migratory high-pressure system are related to the high PM10 and primary pollutant, while the strong irradiance and low NO2 by high winds at the rear of a cyclone are related to the high O3. In terms of long-term trends, decrease in PM10 (−1.75 μg m−3 yr−1) and increase in O3 (+0.88 ppb yr−1) in Seoul were largely contributed by the meteorology-related trends (−0.94 μg m−3 yr−1 for PM10 and +0.47 ppb yr−1 for O3), which were attributable to the subregional-scale wind speed increase. Comparisons with estimated local emissions and socioeconomic indices like gross domestic product (GDP) growth and fuel consumptions indicate probable influences of the 2008 global economic recession as well as the enforced regulations from the mid-2000s on the emission-related trends of PM10 and other primary pollutants. Change rates of local emissions and the transport term of long-term components calculated by the tracer continuity equation revealed a decrease in contributions of local emissions to the primary pollutants including PM10 and an increase in contributions of local secondary productions to O3. The present results not only reveal an important role of synoptic meteorological conditions on the episodic air pollution events but also give insights into the practical effects of environmental policies and regulations on the long-term air pollution trends. As a complementary approach to the chemical transport modeling, this study will provide a scientific background for developing and improving effective air quality management strategy in Seoul and its metropolitan area.


Introduction
Over the past few decades, rapid urbanization, population and economic growth, and increase in energy consumption have exacerbated air pollution in the developing countries of South and East Asia (Sun et al., 2016;Liu et al., 2017;Shi et al., 2018).In response to growing concerns about air quality deterioration, several East Asian countries have implemented strict regulations on emissions in recent decades Published by Copernicus Publications on behalf of the European Geosciences Union.and achieved some degree of improvement in the primary air pollutants.For example, ambient concentrations of carbon monoxide (CO), sulfur dioxide (SO 2 ), and nitrogen oxides (NO x ) have decreased since the 1990s in South Korea and since the mid-2000s in China owing to various efforts to reduce their emissions (e.g., desulfurization in coal-fired power plants and industries, installation of selective catalytic reduction equipment, use of low-sulfur fuels and natural gas, and enhancement of vehicle emission standards; Shon and Kim, 2011;Klimont et al., 2013;Ray and Kim, 2014;van der A et al., 2017;C. Li et al., 2017;Kim and Lee, 2018).Despite the efforts to reduce the primary pollutants and secondary precursors, however, East Asian countries still suffer from frequent severe haze pollution (Huang et al., 2014) and experience continuous increasing of ozone (O 3 ) levels (Seo et al., 2014;Sun et al., 2016).
Such a discrepancy between the past reduction of emissions and the current continuous air pollution can be minimized by considering effects of meteorological conditions on the air quality.The meteorological conditions often play important roles in local air quality through accumulation or ventilation of pollutants, regional transport of polluted or clean air, and atmospheric chemistry for the formation of secondary species (Zheng et al., 2015;Seo et al., 2017).For example, during the Chinese severe-haze episode in January 2013, active secondary aerosol formation in the stagnant surface conditions and shallow boundary layer over the North China Plain induced extremely high particulate matter (PM) concentrations without abrupt changes in emissions (Huang et al., 2014;Zheng et al., 2015;Cai et al., 2017;Zou et al., 2017).The meteorologically adjusted long-term PM trend in Beijing shows that unfavorable meteorological conditions have reduced efficiency of recent emission control (Zhang et al., 2018).Although the meteorological effects on the current increase of O 3 in China are unclear (Ma et al., 2016), the O 3 pollution is expected to become a more important issue in the future warmer climate (Jacob and Winner, 2009;Lin et al., 2008;Schnell et al., 2016) due to its high dependence on temperature (Sillman and Samson, 1995;Lin et al., 2001).
Considering the important role of meteorological conditions for local air pollution levels, long-term urban air quality data must be interpreted carefully with examination of both meteorological effects and changes in local and regional emissions.The meteorological conditions contribute to the various timescale changes in urban air quality by synopticscale weather in the short term (Seo et al., 2017) and climatic variabilities in the long term (Cai et al., 2017;Zou et al., 2017;Oh et al., 2018), as well as its inherent seasonality (Kim et al., 2018).On the other hand, the changes in emissions mostly occur on a long-term timescale due to the implementation of policy and regulations (van der A et al., 2017) and economic booms or recessions (Vrekoussis et al., 2013;Tong et al., 2016).Thus, temporal decomposition of air pollution measurement data onto different timescales provides useful information on the meteorological effects on day-to-day fluctuations, seasonal characteristics, and the long-term trend of the air pollution levels.In addition, the effectiveness of the past regulations on emissions and the influences of socioeconomic changes can be more reasonably assessed by considering and isolating the meteorological effects on air pollution.
As one of the highly populated megacities in East Asia, the Seoul metropolitan area (SMA), which has a population of 25 as well as 9 million vehicles, is obviously a large emission source of various air pollutant species (NIER, 2018).To mitigate air pollution in the SMA by regulation on primary emissions, the South Korean government enacted the Special Act on the Improvement of Air Quality in Seoul Metropolitan Area in 2003 and enforced the act in 2005 (Kim and Lee, 2018).In fact, long-term measurements in Seoul show that concentrations of PM 10 (PM with aerodynamic diameter ≤ 10 µm) and other primary pollutants have decreased over the decades, and these long-term decreasing trends have been regarded as resulting from the enhanced environmental regulations (Ghim et al., 2015;Kim and Lee, 2018).However, considering that the air quality of Seoul is largely affected by the transport of regional air pollutants and the synoptic meteorological conditions (Seo et al., 2017), the efficiency of the current emission control policy should be evaluated after understanding and quantitative isolating the meteorological effects on the long-term measurement data.Using regional air quality simulation, Kim et al. (2017) recently reported that interannual variability in winds could cause the PM 10 levels in Seoul to fluctuate despite the emission control efforts.However, quantification of the meteorological effects based on a direct comparison between measured air pollutants and meteorological variables is further needed to demonstrate the effects of the emission controls.
In this study, we aimed to extensively investigate the role of meteorological conditions on the air pollution in Seoul, South Korea, and quantitatively examine the contributions of meteorology and emissions using the statistical approach to 18-year-long (1999-2016) records of ground PM 10 , NO 2 , SO 2 , CO, and O 3 together, as well as various local meteorological variables.To decompose the air pollution time series onto different timescales and isolate the meteorological effects on the air pollution, we employed a simple statistical concept using the Kolmogorov-Zurbenko (KZ) filter and multiple-linear-regression model with meteorological factors, which has been used in previous studies (Wise and Comrie, 2005;Seo et al., 2014;Henneman et al., 2015;Ma et al., 2016;P. Li et al., 2017;Zhang et al., 2018).Synoptic meteorological conditions related to episodic air pollution events were investigated using composite analyses of the trend-free short-term components, and isolation of the meteorology-and emission-related components from the long-term air pollution trends were conducted.We further explored the long-term changes in contributions of local emissions and transport to the air pollution in Seoul using a simplified tracer continuity equation.Finally, the identified emission-related trends of each pollutant species were compared with the estimated emission inventories, as well as socioeconomic indices that could affect local emissions.

Data
The Korean Ministry of Environment has measured PM 10 and four gaseous air pollutants of SO 2 , nitrogen dioxide (NO 2 ), CO, and O 3 at 264 urban air quality monitoring sites in South Korea and provides their 1 h average concentrations (NIER, 2017;Korea Environment Corporation, 2018).Daily average concentrations for PM 10 , SO 2 , NO 2 , and CO and daily maximum 8-h average concentration for O 3 (O 3 8 h ) at each site were derived from the hourly data.We choose 18 sites located in Seoul to obtain daily air pollutant concentrations representative for the air quality of Seoul (Figs. 1a  and S1a), based on data availability of more than 90 % for 1999-2016.Although the selected 18 sites are spread over Seoul, concentrations from individual sites show fair spatial homogeneity (e.g., coefficient of variations for PM 10 at the selected sites are ∼ 0.17; Fig. S1b and c in the Supplement).Since Seoul is located in a basin surrounded by mountainous terrain (higher than ∼ 500 m above sea level) except its western exit of the Han River toward the Yellow Sea, air pollutants from both emissions and transport can stagnate and mix in the basin in the prevailing westerlies.We averaged daily concentration data from the 18 sites in Seoul that were selected and utilized in this study.
The most noticeable features in the air pollutants data in Seoul are decreasing PM 10 and increasing O 3 levels on a long-term timescale.Annual average PM 10 concentrations and exceedance days of the South Korean air quality standard (AQS; 100 µg m −3 ) were decreased from 69 to 41 µg m −3 and from 64 to 2 days between 1999 and 2015 (Fig. 1b).Note that we excluded episodic Asian dust (AD) days (183 days for the analysis period; KMA, 2018a) from the PM 10 analysis to focus on the anthropogenic sources, although AD did not much affect the long-term PM 10 trend.On the other hand, annual average O 3 8 h levels and exceedance days of the AQS (60 ppbv) have been increased from 25 ppbv to 39 ppb and from 5 to 58 days between 2002 and 2016, respectively (Fig. 1c).
In terms of meteorological characteristics in Seoul, we used daily averages of temperature ( • C), sea level pressure (hPa), relative humidity (%), wind speed (m s −1 ), and solar irradiance (W m −2 ) at the Seoul weather station (Fig. 1a) managed by the Korea Meteorological Administration (KMA, 2018b).Note that the air quality monitoring sites in Seoul are spread over the area within a radius of ∼ 15 km from the weather station, and such spatial size is enough to resolve the influence of synoptic conditions on the local meteorological factors.To investigate the synoptic meteorological conditions, the geopotential height and wind fields at 850 hPa, 10 m wind fields, total cloud cover, and surface solar radiation were derived from the European Centre for Medium-Range Weather Forecasts Reanalysis Interim (ERA-Interim) data (Dee et al., 2011; http:// apps.ecmwf.int/datasets/data/interim-full-daily/,last access: 31 October 2018).
We additionally employed estimated local emission data and several socioeconomic indices that could affect the emissions to compare with the separated emission-related air pollution trends.The annual estimated emissions of sulfur oxides (SO x ), NO x , CO, volatile organic compounds (VOCs), and PM 10 in Seoul from the Clean Air Policy Support System (CAPSS) inventory (Lee et al., 2011;NIER, 2018); the national gross domestic product (GDP) growth (IMF, 2017); and the annual consumptions of the final energy, petroleum, and anthracite in Seoul (KEEI, 2016) were used in this study.
3 Decomposition of air pollutant time series

Temporal decomposition of air pollutant time series
The time series of daily air pollutant concentration at a given place comprises mainly three components: trend, seasonality, and white noise (Rao and Zurbenko, 1994).The trend (long-term component) is attributable to the long-term variations in local and regional emissions related to socioeconomic status and policies (Mijling et al., 2013) and the longterm changes in meteorological conditions, which can affect atmospheric chemistry and regional transport patterns (Cai et al., 2017;Zou et al., 2017).The seasonality (seasonal component) arises from the seasonal variations of meteorological conditions (Kim et al., 2018) and energy consumption pattern (Zhu et al., 2013).The white noise (short-term component) is a day-to-day variation mainly related to synoptic-scale weather changes, which control accumulation/ventilation of local air pollutants and transport of regional air pollutants (Seo et al., 2017), and is partly associated with short-term fluctuations in local emissions (Russell et al., 2010).
To decompose the air pollutant time series into these three components on different timescales, we used the KZ filter, which has been employed in many air pollution time series studies, particularly on O 3 and PM (Wise and Comrie, 2005;Seo et al., 2014;Ma et al., 2016;P. Li et al., 2017).The KZ filter, here denoted as KZ (m, p) , represents p times iteration of the moving average of time width m (days) and removes the high-frequency component of periods that are smaller than the effective filter width, N (≥ m × p 1/2 ).The KZ filter method is applicable to the time series with missing data owing to the iterative moving average process and provides a high accuracy level comparable to that of the wavelet transform method, despite its simplicity (Eskridge et al., 1997).Here we applied KZ (15,5) and KZ (365, 3) filters, which can remove variabilities of periods shorter than 33 days and 1.7 years, to filter out the short-term component and to leave the long-term component, respectively.The long-term component can be further separated into the emission-related trend and the meteorology-related trend by isolation of the emission-related component using a multiplelinear-regression model with representative meteorological variables.
Note that the original concentration (χ) was transformed into its natural logarithm values (X = lnχ) prior to the decomposition.Because the number distribution of daily air pollutant concentrations is usually lognormal (e.g., Fig. S2), the natural log transformation of the original concentration data is required for proper temporal decomposition with the KZ filter (Eskridge et al., 1997).The detailed decomposition procedure is described in the following and schematically summarized with PM 10 time series in Seoul in Fig. 2.
The time series of log-scaled pollutant concentration can be expressed as the sum of short-term (X ST ), seasonal (X SN ), and long-term (X LT ) components.
The sum of seasonal and long-term components is a baseline (X BL = X SN + X LT ).X BL and X ST can be easily decomposed by applying the KZ (15, 5) filter to X, which filters out the white-noise-like X ST , as follows: Now the baseline can be assumed to consist of its repeated climatological seasonal cycle X clm BL and residuals (ε).
Although X clm BL obviously occupies most of the seasonality in X BL , ε also contains some minor seasonal variability unconsidered in X clm BL together with the long-term trend.To obtain the long-term component (X LT ) by filtering out the minor seasonality, the KZ (365, 3) filter is applied to ε as follows: Then the seasonal component (X SN ), which represents the sum of the pure seasonal climatology X clm BL and the minor seasonality (ε − KZ (365, 3) [ε]), can be obtained as the difference between X BL and X LT .
Note that, if we define χ BL = exp (X BL ) and χ ST = exp (X ST ) and employ the same concept to the relationship between the original concentration and its log transformation (χ = exp (X)), χ BL represents the baseline concentration of the air pollutant and χ ST becomes the ratio of original concentration to baseline concentration (χ /χ BL ).Similarly, exp (X SN ) and exp (X LT ) can be defined as χ SN and χ LT , respectively.Then χ SN represents the seasonal change in concentration without trend, and χ LT becomes the ratio of baseline concentration to detrended seasonal concentration (χ BL /χ SN ).
As shown in the example with PM 10 in Seoul, the KZ (15, 5) filter effectively removes PM 10 ST of periods shorter than 33 days and leaves both the seasonality of high PM 10 concentrations in winter and spring and the long-term decreasing trend in PM 10 BL (Figs. 2 and S4b).PM 10 SN and PM 10 LT represent well the seasonal variation, i.e., periods that are between 33 days and 1.7 years with representative periodicities of 0.5 and 1 years, and the long-term variations, i.e., periods that are longer than 1.7 years, respectively (Figs. 2 and S4cd).The high levels in winter and spring of PM 10 SN in Seoul are attributable to the shallow boundary layer that traps local pollutants near the ground and frequent regional transport from China during the cold season (Kim et al., 2018).

Separation of emission-and meteorology-related trends
Since the long-term variability in air pollutant concentrations can be affected not only by changes in local and regional emissions but also by changes in meteorological conditions, the long-term trend is assumed to consist of meteorologically adjusted (emission-related) long-term components X emis LT and meteorology-related long-term components X met LT .Therefore, the baseline can be represented as follows: To isolate the term X emis LT in Eq. ( 5), we built a multiplelinear-regression model employing the baseline time series of the five representative meteorological variables (MET BL ), such as temperature (T BL ), sea level pressure (P BL ), relative humidity (RH BL ), wind speed (WS BL ), and solar irradiance (SI BL ), which are obtained by the KZ (15, 5) filter as follows: where ε is the sum of the non-meteorological long-term variability X emis LT and the minor seasonal variability unexplained by the multiple-linear-regression model ε − X emis LT .Note that daily maximum temperature (T max BL ) was used for the model of O 3 8 h BL instead of daily average temperature (T BL ) considering daytime O 3 as in a previous study (Seo et al., 2014).By removing the minor seasonality from ε using the KZ (365, 3) filter, X emis LT can be isolated as follows: Then X met LT can be simply obtained by differencing between X LT and X emis LT .In the example displayed in Fig. 2, PM 10 LT in Seoul shows a continuous decrease between 2003 and 2012.Such a decreasing trend in PM 10 has been recognized as a result of the reduction in diesel vehicle emissions and fugitive dust in Seoul (Ghim et al., 2015;Kim and Lee, 2018).However, a recent modeling study suggested that the long-term increase in wind speed might additionally contribute to the past improvement of PM 10 air quality in Seoul (Kim et al., 2017).In fact, both PM emis 10 LT and PM met 10 LT in Fig. 2 show decreasing patterns, and this supports probable influences of both emission controls and meteorology on the PM 10 trend in Seoul.

Contributions of local emissions and transport to
the long-term trends and changes in transport of regional emissions X emis(T) LT , and thus X LT can be represented as Eq. ( 8).
Then the rate of change of X LT should satisfy a simple continuity equation as follows: where the advection term, −V LT •∇X LT , represents the transport of regional emissions mainly by horizontal winds, and the sources and sinks term, S LT , can be regarded as the longterm changes by both local emissions and chemical production, accumulation, and dissipation by meteorological factors.Because X LT and X met LT were already known, the rates of change of X emis(T) LT and X emis(L) LT can be derived as follows (see Sect.S1 in the Supplement): where ) is long-term zonal and meridional winds components representative for the target area, and ∇X LT is the horizontal gradient of long-term components of air pollution for the larger area, which are identical to the difference between ∇X BL and ∇X SN .If ∇X LT > 0 at a specific time, then the horizontal gradient of baseline (∇X BL ) at that time is steeper than that of seasonal climatology (∇X SN ).Note that the above method for evaluating the changes in X Year Low-pass filter (KZ (15,5) filter) Climatological seasonal cycle of !+$ (! +$ 34' ) and Low-pass filter (KZ (365,3) filter)

Multiple linear regression model and
Low-pass filter (KZ (365,3) filter) . Schematic flowchart of temporal decomposition of air pollution time series (Seoul average daily PM 10 concentration) into shortterm, seasonal, and emission-related and meteorology-related long-term components.
X BL and X SN at those sites were utilized to determine ∇X BL and ∇X SN by linear regressions of X BL and X SN on the zonal and meridional distances from the weather station and finally to obtain ∇X LT (= ∇X BL − ∇X SN ; Figs.S5 and S6c, e, g, i, and k).Also, u LT and v LT in Seoul are calculated by the KZ (365, 3) filter with wind direction and speed at the Seoul weather station (Fig. S6a and b).From V LT and ∇X LT of each species, we calculated the long-term transport terms for each air pollutant species in Seoul (Fig. S6d, f, h, j, and l).
In Eq. ( 11 ∼ −2.5 × 10 −9 s −1 .If we further assume that the local PM 10 emission in Seoul (area of 605 km 2 ) has the same order of magnitude as the transport term for the concentration u ∂χ ∂x ∼ 1 × 10 −4 µg m −3 s −1 and the boundary layer height over the area is about 1 km, the total amount of yearly PM 10 emission in Seoul will be approximately 1.9 kt, which is consistent with the estimation by the CAPSS emission inventory (Fig. 5a; NIER, 2018).In fact, because air pollutant concentrations are close to the steady state for the long-term period, its tendency is a small residual between two large terms: influx of clean or polluted air (transport) and local emissions/production or dissipation/deposition (sources and sinks).It should be noted that the transport term and the source/sink term, however, are not balanced against each other on a short-term timescale because the day-to-day rate of change of the air pollution concentration is comparable to both transport and source/sink terms.

Application to air pollutants in Seoul
Variances of each timescale component of PM 10 , SO 2 , NO 2 , and CO in Seoul were generally largest in the short-term fluctuation (∼ 46-68 %), while that of O 3 8 h was largest in its seasonal cycle (Table 1).This indicates an important role of synoptic-scale weather on the day-to-day variations of the primary air pollutants concentrations in Seoul.The changes in synoptic flow pattern control local air quality by accumulation/ventilation and regional transport of air pollutants (Zheng et al., 2015;Seo et al., 2017).On the other hand, O 3 (a) -3 lag days (b) -2 lag days is a photochemical product and thus is controlled largely by the annual cycle of solar irradiance in such a NO x -and VOCrich urban area, although its short-term variability, which is caused by the changes in synoptic weather and related precursor concentrations and irradiance, is comparable to the seasonal variability.
In terms of the long-term trends, variabilities related to the long-term local and regional emission changes X emis LT are comparable to those induced by the long-term meteorological effects X met LT for PM 10 , NO 2 , CO, and O 3 8 h , while X emis LT is dominant for total long-term variabilities of SO 2 (Table 1).Since the local-scale meteorology could affect more the local emission-related long-term trend rather than the regional background-related long-term trend, the smaller variability of X met LT implies a smaller influence of the local emission changes compared to the regional background level changes on X emis LT .

Meteorological effects on seasonality in air pollution
In the process of isolating X emis LT from X LT , the multiplelinear-regression model predicting X BL from meteorological predictors (MET BL ) was used.The model using the meteorological variables explains well O 3 8 h BL (adjusted R 2 = 0.86), while the baseline of PM 10 (PM 10 BL ) is much less explained by the model (adjusted R 2 = 0.51; Table 2).This indicates that the long-term non-meteorological (or emission-related) variability in PM 10 BL is much larger than that in O 3 8 h BL .In the linear regression model, PM 10 BL and the baselines of primary gaseous pollutants of SO 2 , NO 2 , and CO commonly show strong positive correlation with P BL but strong negative correlations with T BL and RH BL , while O 3 8 h BL shows strong positive correlations (p < 0.05) with SI BL and T max BL but strong negative correlation (p < 0.05) with P BL (Table 2).Meteorological conditions over the Korean Peninsula in winter are affected by a cold and dry Siberian high that causes both the shallow boundary layer to trap the primary air pollutants near the surface and the northwesterly winds to transport the regional pollutants from the continent (Kim et al., 2018) and thus elevates PM 10 BL in winter to early spring (Fig. S3b).On the other hand, enhanced photochemistry by strong irradiance together with temperature effects on O 3 formation (increasing of biogenic hydrocarbons and hydroxyl radicals and enhanced thermal decomposition of peroxyacetyl nitrate in warm conditions; e.g., Sillman and Samson, 1995) elevates O 3 8 h levels in late spring and summer in Seoul (Fig. S3b).

Synoptic influences on short-term air pollution events
The short-term air pollution variability is closely related to the day-to-day variations of local meteorological factors.Intercorrelations among X ST of air pollutant species and meteorological variables in Seoul are summarized in Table 3.Note that using X ST instead of the original daily concentrations has the advantage of providing short-term features unbiased to seasonal characteristics and background levels because exp (X ST ) is equivalent to the ratio of measured concentration to filtered baseline concentration as aforementioned in Sect.3.1.Significant positive intercorrelations (p < 0.01) among PM 10 ST and primary gaseous pollutants (SO 2 ST , NO 2 ST , and CO ST ) together with their strong negative relationships to T ST and WS ST (p < 0.01) indicate that the high-PM 10 episodes in Seoul occurred in warm and stagnant weather conditions.Such warm and stagnant conditions can be induced by the migratory high-pressure system (e.g., negative relationships of WS ST to T ST and P ST in Table 3) or by the warm front of the extratropical cyclone (e.g., negative correlations of T ST to WS ST and P ST but its positive correlation to RH ST in Table 3).The lag composites of the geopotential height and wind anomalies at 850 hPa (about 1.5 km altitude) relative to each day of exp PM 10 ST , which is equivalent to the ratios of measured concentration to the filtered baseline concentration of PM 10 , exceeding the sum of its mean and 1 standard deviation (> 1.50) indicate that the high-PM 10 episodes in Seoul are associated with the former conditions (Fig. 3).As shown in Fig. 3a-c, a high-pressure anomaly develops over southern China and moves eastward in 3 days before the high-PM 10 days.This kind of synoptic pattern can induce both slow regional transport of secondary precursors from China along the northern boundary of the highpressure system (Fig. 3b and c) and gradual accumulation of primary and secondary aerosols in a local area in the highpressure system (Fig. 3d; Seo et al., 2017).Strong correlations of PM 10 ST with SO 2 ST and NO 2 ST (p < 0.01; Table 3) imply possible large contributions of secondary species such as sulfate and nitrate aerosols to the haze episodes in Seoul.
In fact, previous studies have reported that the proportion of secondary inorganic aerosols to fine-mode PM (PM 2.5 ) was measured as high as ∼ 50 % on average and even reached ∼ 75 % during the severe-haze event (Seo et al., 2017;Kim et al., 2018).
In terms of short-term variability in O 3 8 h levels in Seoul, O 3 8 h ST is strongly correlated positively with SI ST and T max ST but negatively with RH ST (p < 0.01; Table 3).This likely reflects favorable meteorological environment for O 3 formation: strong irradiance with warm and dry conditions.Interestingly, O 3 8 h ST shows positive correlation with WS ST , although stagnant conditions are regarded as conducive to O 3 formation in general.This is related to the strong negative relationship between O 3 8 h ST and NO 2 ST and that between WS ST and NO 2 ST .Since the reaction with NO, which is produced by photolysis of NO 2 , is a major sink of O 3 , high NO 2 concentration reduces O 3 on high-irradiance days, while the low NO 2 concentration can increase O 3 .In addition, NO 2 ST and CO ST are strongly correlated (r = +0.86)because both CO and NO x are mainly emitted by transportation exhaust in Seoul (82 % of CO and 67 % of NO x emissions; NIER, 2018), and thus O 3 8 h ST has a negative relationship with CO ST , similar to that with NO 2 ST .The composites of the 850 hPa geopotential height and wind anomalies for high-O 3 8 h days that exp O 3 8 h ST exceeds the sum of its mean and 1 standard deviation (> 1.41) show that the high-O 3 8 h event frequently occurs at the transition of weather from the cyclonic anomaly to the anticyclonic anomaly (Fig. 4a).As the anomalous low-pressure system moves eastward out of the Korean Peninsula, the decrease in cloud cover (Fig. 4b) results in the increase in surface irradiance (Fig. 4c) and thus the O 3 level.Since this kind of anomalous synoptic pattern enhances the mean westerly flow by the anomalous northwesterly (Figs.4a and S7a), the wind speed in this region is increased.As shown by the strong positive correlations between WS ST and NO 2 ST (Table 3), NO 2 can be reduced by the high winds and thus could contribute to the high O 3 concentration on short-term timescales.

Long-term trends of air pollution in Seoul
Times series of X LT and its two decomposed components, X emis LT and X met LT , for each air pollutant species are shown in Fig. 5, and their linear trends are summarized in Table 4.Note that these time series range from July 2000 to June 2015 because 546 days at the beginning and ending of data were lost due to the truncation effect of the KZ (365, 3) filter.The linear trend of X LT represents a fractional change rate (% yr −1 ) of the baseline concentration (χ BL ) because ∂X LT ∂t is equivalent to 1 χ BL ∂χ BL ∂t .The fractional change rate can be converted into an equivalent concentration change rate by multiplying with temporal-mean χ BL for the analysis period.
In Seoul, there are significant decreasing trends (p < 0.05) in PM 10 LT (−3.6 % yr −1 ) and CO LT (−2.9 % yr −1 ) and an increasing trend (p < 0.1) in O 3 8 h LT (+3.1 % yr −1 ) for the last 15 years, while NO 2 LT (−1.4 % yr −1 ) and SO 2 LT (+0.8 % yr −1 ) do not show statistically significant trends (Table 4).The decrease of PM 10 concentration and increase of O 3 8 h level in their long-term linear trends can be clearly identified in temporal variation of X LT (Fig. 3a and e) and are consistent with the temporal characteristics in their annual average concentrations (Fig. 1b and c).In terms of CO, abrupt decrease in the early 2000s, which was reported to be associated with introduction of natural gas vehicle supply and improvement of fuel quality (Kim and Shon, 2011), affects such a strong linear trend (Fig. 5b).On the other hand, SO 2 concentration in Seoul already stabilized at ∼ 5 ppb during the last decade (Fig. 5c), although it decreased significantly in the 1980s and 1990s owing to the government's efforts to control the SO x emissions by use of low-sulfur fuel and natural gas for industry and transportation (Ray and Kim, 2014;NIER, 2017).NO 2 concentration varied between 30 and 40 ppb in the last decade and also shows a relatively weak trend (Fig. 5d; NIER, 2017).Interestingly, the NO 2 level has stabilized despite an increase of the number of vehicles in Seoul from 2.3 million in 1999 to 3.1 million in 2016 probably owing to implementation of natural gas vehicles and low-emission diesel engines, and the NO x (= NO+NO 2 ) level has even decreased from ∼ 70 to ∼ 50 ppb for the same period (Shon and Kim, 2011;Kim and Lee, 2018).Such an increase of NO 2 to NO x ratio implies that additional conversion of the NO-to-NO 2 occurs somewhere before the emission (e.g., exhaust line of the vehicle) or in the atmosphere.Although further evidence is required, this could be attributable to the expansion of diesel particulate filter (DPF) and diesel oxidation catalyst (DOC) usage for diesel vehicles or the increase of atmospheric oxidative potential in the SMA (Alvarez et al., 2008;Kim and Lee, 2018).Since Seoul is a NO x -saturated regime area (Jin et al., 2012), such a decreasing trend in ambient NO x concentration in recent decades may be one of the causes of the long-term increase of O 3 8 h in Seoul.
Although the decreasing trends in PM 10 LT , CO LT , and NO 2 LT in the last decade have been regarded as the result of efforts to reduce the local emissions (Kim and Lee, 2018), their X emis LT linear trends are less than half of X LT linear trends (Table 4).This indicates the more important role of the longterm changes in local meteorology in the long-term air pollution trends in Seoul.
In contrast to the comparable linear trends of X emis LT and X met LT of PM 10 , CO, and NO 2 , the linear trend of SO 2 LT is dominantly contributed by that of SO emis 2 LT (Table 4).This suggests a larger influence of the regional background SO 2 than that of the local SO x emission on SO emis 2 LT because the long-term effects of local meteorological conditions on the local emission-related air pollution trend must be limited if the long-term change in local emission is negligible.In fact, the estimated SO x emission intensity of Seoul (7.4 t km −2 ) in the reference year 2010 was much lower than that of an industrial port city in the west adjacent to Seoul, Incheon (18.1 t km −2 ), or the Chinese eastern coastal region of Jing-Jin-Ji (Beijing-Tianjin-Hebei), Shandong, Jiangsu, and Shanghai (14.3 t km −2 ; M. Li et al., 2017;NIER, 2018).In consideration of the prevailing westerlies in this region (Fig. S7a), the long-term SO 2 concentration trend in Seoul could be more affected by the long-term change in regional background level than by that in local emission.
In the reference year 2010, emission intensities of CO and NO x in Seoul in 2010 were estimated as 215.3 t km −2 and 117.4,respectively, and these are obviously much higher than those in Incheon (44.1 t km −2 for CO and 51.3 t km −2 for NO x ) or the Chinese eastern coastal region (98.1 t km −2 for CO and 17.4 t km −2 for NO x ; M. Li et al., 2017;NIER, 2018).Therefore, CO emis LT and NO emis 2 LT in Seoul are probably much more affected by the long-term changes in local emissions than by the changes in regional background levels.

Meteorology-related long-term trends
For most of the species in Seoul except SO 2 , larger decreasing trends of X met LT than of X emis LT have been found (Table 4).For example, the equivalent linear trend of PM 10 LT was −17.5 µg m −3 decade −1 , while that of PM emis 10 LT was only −8.1 µg m −3 decade −1 .Thus, the recent achievement of ambient PM 10 reduction in Seoul would be less successful without the contribution by the PM met 10 LT trend (−9.4 µg m −3 decade −1 ; Table 4 and Fig. 5a).In terms of O 3 , the equivalent linear trend of O 3 8 h LT (+8.8 ppb decade −1 ) was also contributed to slightly more by O met 3 8 h LT (+4.7 ppb decade −1 ) than by O emis 3 8 h LT (+4.1 ppb decade −1 ; Table 4 and Fig. 5e).Therefore, the long-term changes in meteorology that help improve PM air quality may have made the O 3 air quality worse.Since the only meteorological parameter that shows a significant (p < 0.1) linear trend in Seoul is wind speed (+0.57m s −1 decade −1 ; Table 4 and Fig. 6a), the meteoro-logical effects on the long-term decreases of PM 10 , CO, and NO 2 are probably related to the long-term increase of wind speed, which could result in enhancing ventilation of the primary pollutants.The long-term component of wind speed (WS LT ) in Seoul was increased by ∼ 0.9 m s −1 from 2002 to 2012 (Fig. 6a).Previous statistical research on the dispersion effects of winds reported that the fraction of concentrations reduced by wind speed rising 2 to 3 m s −1 was ∼ 24 % for PM 10 and ∼ 9 % for NO 2 (Lim et al., 2012).These ratios are comparable to the decrease of meteorology-related long-term concentrations χ met LT in PM 10 and NO 2 by ∼ 30 % and ∼ 10 % for the period, respectively ( X met LT (= χ met LT /χ met LT ) is −0.3 for PM 10 and −0.1 for NO 2 ; Fig. 5a and d).Such a role of wind speed in recent changes in PM concentrations in Seoul was also indicated by regional air quality model simulations (Kim et al., 2017).Interestingly, the recent increasing of wind speed in Seoul is a subregional-scale change rather than synoptic-or larger-scale climate variability.Linear trend distribution of observed wind speed over South Korea for 1999-2016 shows a rough pattern of increasing in the inland area and decreasing in the coastal area (Fig. 6b).However, the significant long-term increases of wind speed near Seoul seem to be somewhat limited to the Han River basin.In addition, although the linear trend pattern of atmospheric circulation at 850 hPa (∼ 1.5 km of altitude) over the western North Pacific/East Asia region shows weakening of the Aleutian low during the period (Fig. S7a and b), no significant trend can be observed either in 850 hPa wind fields or at 10 m wind speeds over and near the Korean Peninsula (Fig. S7b and c).
Since the major increase in wind speed was observed during two periods (2001-2004and 2010-2011;Fig. 6a), changes in PM met 10 LT , CO met LT , and NO met 2 LT mostly occurred in these periods.Such an increase in wind speed could result in more ventilation of NO x and elevation of O 3 levels in the NO x -saturated regime.However, O met 3 8 h LT shows additional variability related to the long-term variation of solar irradiance (e.g., decrease for 2005-2006 and increase for 2007-2009 in both O met 3 8 h LT and SI LT , Figs. 5e and 6a).

Emission-related long-term trends
Although PM 10 LT continuously decreased between 2003 and 2012, a major decrease of PM emis 10 LT occurred between 2008 and 2009 (Fig. 5a).Since the decrease of X emis LT during the period can be found in other primary gaseous pollutants (Fig. 5b-d), it could be related to the reduction of local and/or regional primary emissions.For example, the CAPSS emission inventory shows abrupt decreases in estimated SO x , NO x , and PM 10 emissions in Seoul from 2007 to 2009 (Fig. 7a).One plausible reason for such a reduction of primary emissions during the period could be enforcement of the Special Act on the Improvement of Air Quality in Seoul Metropolitan Area, which took effect in January show that the horizontal gradients of X LT (∇X LT ) obtained by linear regression are statistically significant at the 95 % level or higher (p < 0.05).
trations to meet 40 µg m −3 and 22 ppb by the year 2014, by regulations of NO x , SO x , PM 10 , and VOC emissions mainly from industry and transportation (Ghim et al., 2015;Kim and Lee, 2018).Although the act includes the introduction of the emission cap-and-trade system and strengthening VOC management, most of the budget has been allocated for the expansion of DPF and DOC usage for old diesel vehicles, which was considered to be effective for reduction of NO x and primary PM emissions (Kim and Lee, 2018).The impact of this policy on the decrease in PM 10 and NO x emissions (Fig. 7a) is not easily distinguishable from the influence of the decrease in diesel consumption (Fig. 7c).However, the generalized use of DPF and DOC in diesel vehicles may be one reason for relatively stable trends of PM 10 and NO x emissions despite the rapid increase in diesel consumption mainly by vehicles since 2012.
Socioeconomic status can also be an important factor for emissions and air pollution (Vrekoussis et al., 2013;Tong et al., 2016;Stern and van Dijk, 2017).In terms of local emissions, the rapid reduction of primary emissions in Seoul between 2008 and 2009 coincided with the rapid drop of South Korean GDP growth (Fig. 7b), which was associated with the 2008 global financial crisis triggered by the US subprime-mortgage meltdown (Kim et al., 2017).Final energy consumption in Seoul shrank during the recession, primarily due to a decrease in consumption of petroleum products (especially diesel and liquefied petroleum gases; Fig. 5b; KEEI, 2016).The decrease in PM 10 , NO x , and SO x emissions in Seoul in 2007-2009 was attributable to the reductions of diesel consumption and anthracite consumption (Fig. 7a and c).In terms of regional emissions, temporal variation of SO emis 2 LT , which shows a large decrease in 2008-2009 and slight recovery in the following brief period (Fig. 5c), corresponds to interannual variation of satellite SO 2 observation over China (van der A et al., 2017) omy, as well as the desulfurization program of the Chinese government (van der A et al., 2017;C. Li et al., 2017).Temporal variation of O emis 3 8 h LT shows a general increasing trend throughout the analysis period (Fig. 5e).Considering NO x -rich conditions in Seoul, this may result from the long-term decrease of NO emis 2 LT (Fig. 5d) and relatively greater reduction of NO x emission compared to that of VOC emission in the last decade (Fig. 7a).However, the factors that caused its steep increase in 2005, slight decrease in 2010, and following recovery between 2011 and 2012 are somewhat www.atmos-chem-phys.net/18/16121/2018/vague, without further information on concentration and temporal variation of VOCs in Seoul.

Long-term effects of local emissions vs. transport
The change rates of X emis(L) LT and X emis(T) LT have nearly the same magnitude but opposite signs.Since magnitudes of the change rates of X LT , X met LT , and X emis LT are much smaller compared to those of X emis(L) LT and X emis(T) LT on a long-term timescale, characteristics of their temporal variation is not easy to explain using X > 0 shows gradual reduction of local emissions and compensation by an increase in transport of regional background emissions.
In Seoul, ∂ ∂t X

emis(L) LT
> 0 for PM 10 and primary gaseous pollutants (CO, SO 2 , and NO 2 ) in the early 2000s, but the change rates got close to zero owing to the efforts to reduce the primary emissions (Fig. 5).Unlike CO and NO 2 , of which emission intensities are large in Seoul, SO 2 shows obvious positive ∂ ∂t SO emis(T) 2 LT in recent years (Fig. 5c), and this indicates the large effect of regional transport on the variability of SO 2 concentrations in Seoul.Similar to SO 2 , PM 10 also shows a recent increase in PM  (Shin et al., 2012).However, because O emis(L) 3 8 h LT has gradually changed from the decreasing phase to the increasing phase over the analysis period (Fig. 5e), the recent changes in O emis 3 8 h LT in the 2010s are probably more related to the local secondary production than to the background O 3 transport.
Interestingly, the change rates of X emis(L) LT reflect some changes in socioeconomic factors.For example, PM emis(L) 10 LT and SO emis(L) 2 LT switched to the decreasing phase in 2008 (Fig. 5a and c), corresponding to the aforementioned descriptions in the previous subsection.In addition, the variation of  5d), and these changes correspond to the diesel consumption in Seoul, which decreased from 2007 to 2012 but has increased since then due to policy support for the diesel vehicle market to reduce greenhouse gas emissions (Fig. 7c).

Conclusions
Based on the statistical approach to the long-term air pollutant measurement data from the urban air quality monitoring sites in Seoul, the present study revealed the important role of synoptic weather conditions on the episodic air pollution events and the meteorological effects on the long-term air pollution trends.Temporal decomposition using the KZ filter technique and multiple linear regression with meteorological factors allowed separation of the trend-free short-term variability and isolation of the meteorology-and emissionrelated long-term trends in the daily air pollution data.The simplified continuity equation using the surface wind data in Seoul and air pollutant measurement data from the wider area around Seoul further provided the approximate estimation of the long-term changes in contributions of local emissions and transport to the air pollution in Seoul.
In terms of the short-term variabilities, which occupy the largest portion of the total air pollution variabilities, the high PM 10 and primary gaseous pollutants concentrations are related to the influence of migratory high-pressure systems, which induce both regional transport and local accumulation of the pollutants in warm and stagnant conditions.On the other hand, the high-O 3 episodes are related to the weather transitions at the rear of cyclones, which accompany the NO 2 reduction by high winds and the strong solar irradiance.
In terms of the long-term trends, Seoul experienced a decrease in PM 10 concentrations (−17.5 µg m −3 decade −1 ) and an increase in O 3 8 h levels (+8.8 ppb decade −1 ) over the period of 1999-2016.Such long-term changes are largely contributed by the meteorology-related trends (e.g., −9.4 µg m −3 decade −1 for PM 10 and +4.7 ppb decade −1 for O 3 8 h ) related to the decadal increase in surface wind speeds (+0.57m s −1 decade −1 ).Therefore, the recent improvement of particulate air quality in Seoul was not achieved solely by the emission control policies but was also induced by changes in local meteorological conditions, especially wind speeds.Although no evidence of the influence of synopticor larger-scale climate variability can be found, the long-term wind speed increase was also observed in the Han River basin around Seoul and probably enhanced the ventilation of local primary pollutants and secondary precursors.Since Seoul is a NO x -saturated regime, the long-term decrease in NO 2 by the enhanced winds could be a cause for the long-term increase in O 3 8 h levels.Unlike other species, SO 2 does not show a significant meteorology-related trend because of the small influence of local meteorology on its small local emissions.
The isolated emission-related air pollution trends of PM 10 and other primary gaseous pollutants commonly showed a major decrease between 2008 and 2009.Although the enforcement of the Special Act on the Improvement of Air Quality in Seoul Metropolitan Area in 2005 might affect the reduction of local emissions, socioeconomic indices like GDP growth and energy consumptions also indicate probable influence of the 2008 global economic recession on the changes in the emission-related air pollution trends.Contributions of the local emissions and the transport of the regional background air pollutants to the emission-related long-term trends almost balanced each other out because the change rate of the air pollution trend is negligibly small compared to that of the local emissions and the transport in the long term.The contributions of local emissions to PM 10 and primary gaseous pollutants changed to the decreasing phase in the mid-2000s in Seoul, and those species, especially PM 10 and SO 2 , have become more affected by the regional background levels in recent years.The contribution of local emissions to NO 2 changed to the increasing phase in the 2010s due to the recent policy support for diesel vehicles.Also, the recent increasing phase of local contributions to O 3 implies an important role of local secondary production in the increasing trend of O 3 in Seoul.
Although the results from the emission-related long-term air pollution trends showed general agreement with estimated emissions and socioeconomic factors, the statistical technique employed in this study overlooked the role of longterm changes in the chemical conditions and reactions on the urban aerosols and O 3 (e.g., atmospheric oxidants, aerosol acidity and hygroscopicity, and secondary organic aerosol formation) and their relationship with the long-term changes in relevant meteorological factors (e.g., atmospheric water, irradiance, and temperature).Nevertheless, this simple concept of considering the emission-related air pollution trends separately from the meteorology-related trends could give insights into the assessment of past regulations on emissions and help the improvement of current environmental policies on urban air quality.As a complementary approach to chemical transport modeling, our analysis will provide a scientific background for effective air quality management strategy in the SMA.
Data availability.The hourly data of PM 10 , SO 2 , NO 2 , CO, and O 3 concentrations over South Korea for 1999-2016 are available on the website managed by the Korea Environment Corporation (2018; https://www.airkorea.or.kr/last_amb_hour_data).The hourly data of temperature, sea level pressure, relative humidity, wind speed, and solar irradiance at the Seoul weather station for the same period can be found on the website of the KMA (2018b; https: //data.kma.go.kr/data/grnd/selectAsosRltmList.do?pgmNo=36).The ERA-Interim data (Dee et al., 2011) can be accessed via the European Centre for Medium-Range Weather Forecasts (ECMWF) data server (http://apps.ecmwf.int/datasets/data/interim-full-daily/). The daily records of the Asian dust events at Seoul can be obtained from the KMA website (2018a; http://www.weather.go.kr/weather/asiandust/observday.jsp?type= 2&stnId=108&year=2016&x=20&y=11).
Author contributions.JS and DY designed and initiated the study.JS and DSRP performed the data analysis and interpreted the results.YBL, JYK, and YK provided additional feedback on the manuscript.JS prepared the manuscript with contributions from all co-authors.
Competing interests.The authors declare that they have no conflict of interest.

Figure 1 .
Figure 1.(a) Geographical locations of air quality monitoring sites in the SMA (circles) and the Seoul weather station (a blue solid diamond).Red solid circles denote 18 air quality monitoring sites utilized to obtain average daily air quality data in Seoul.Yellow solid circles show 52 air quality monitoring sites, which together with 18 sites (red solid circles) were used to calculate spatial gradients of the long-term components of air pollutants.Boundary of Seoul is marked with green line, and the satellite image is a courtesy of Google Earth Pro.Annual average concentrations and exceedances of the World Health Organization (WHO) guidelines and the South Korean AQS for (b) PM 10 (excluding AD days) and (c) O 3 8 h in Seoul.Annual average PM 10 concentrations including AD days are additionally represented as an orange line in (b).The vertical bars on the annual average concentrations indicate annual standard deviations.

Figure 3 .
Figure3.Lag composites of seasonal anomalies of geopotential height (contours with interval of 2.5 gpm) and wind (arrows with reference scale of 1 m s −1 ) at 850 hPa relative to each first day that exp PM 10 ST exceeds the sum of its mean and 1 standard deviation.Total number of events, mean, and standard deviation is 283, 1, and 0.50, respectively.Geopotential height and wind statistically significant at the 95% confidence level are represented as light gray shading and arrows.Location of Seoul is marked by solid red circles.

Figure 4 .
Figure 4. Composites of seasonal anomalies of (a) geopotential height (contours with interval of 2.5 gpm) and wind (arrows with reference scale of 1 m s −1 ) at 850 hPa, (b) total cloud cover fraction (contours with interval of 2 %), (c) downward solar radiation anomaly at the surface (contours with interval of 5 W m −2 ) for each first day that exp O 3 8 h ST exceeds the sum of the mean and 1 standard deviation.Total number of events, mean, and standard deviation is 346, 1, and 0.41, respectively.Geopotential height and wind, total cloud cover, and downward solar radiation statistically significant at the 95 % confidence level are represented by light gray shading and arrows.Location of Seoul is marked by a solid red circle.

Figure 5 .;
Figure5.Long-term components that are unadjusted for meteorological variables (X LT ; violet lines), meteorology-related X met LT ; blue lines and emission-related X emis LT ; redlines long-term components, and contributions of local emissions ∂ ∂t Xemis(L) LT; orange lines and trans-

Figure 6 .
Figure 6.(a) Long-term components of daily average temperature, sea level pressure, relative humidity, wind speed, and solar irradiance.Vertical scales represented by gray horizontal lines are equivalent to 0.3 standard deviation of each original time series of the meteorological variables.(b) Linear trends of wind speed observed at 75 weather stations over South Korea for the period of 1999-2016.Upward and downward triangles with colors indicate the increasing and decreasing trends which are statistically significant at or above the 95 % confidence level (p < 0.05), respectively.Location of Seoul is marked by a red open circle.
these two terms could provide some physical and quantitative insights into temporal changes in each contribution of local emissions and transport of regional background emissions to the local air quality, albeit without considering atmospheric chemistry.For example, implications in the long term, increasing contributions of local emissions and, at the same time, flushing out air pollutants.On the other hand, during the 2000s (Fig.5e) implies that the transport of regional background O 3 played an important role in the meteorologically adjusted long-term O 3 8 h trend O emis 3 8 h LT during the period.This result is consistent with the previous study on O 3 in Seoul for 2002-2006 using the OZone Isopleth Plotting Package for Research (OZIPR) and the KZ filter well the changes in diesel consumption, which were not clearly shown in the CAPSS inventory.The change rate of NO emis(L) 2 LT was negative between 2006 and 2011 but turned positive in 2012 (Fig.

Table 1 .
Total variances (Var (X)) of log-scale times series for Seoul average concentrations of PM 10 , SO 2 , NO 2 , CO, and O 3 8 h , and relative contributions (%) of variances of and covariances (Cov) among each component to Var (X).Daily data for the period of July 2000 to June 2015 during which X LT data are available were used.
* Values of variances of each pollutant time series Var

Table 2 .
Correlation coefficients (r) between baseline time series of pollutants (PM 10 , SO 2 , NO 2 , CO, and O 3 8 h ) and meteorological variables (T , T max , P , RH, WS, and SI), and adjusted R 2 between the baseline time series of pollutants (X BL ) and multiple-linear-regression model (a 0 + i a i MET BL i ).
*The correlation is statistically significant at the 95 % level or higher (p < 0.05).

Table 3 .
Correlation coefficients (r) among short-term components of each pollutant and meteorological variable in Seoul for the period between February 1999 and November 2016.

Table 4 .
Linear trends of long-term components of air pollutants and meteorological variables in Seoul for the period between July 2000 and June 2015.
*The slope of linear trend is statistically significant at the 90 % level or higher (p < 0.1).