3.1 Results of chemical compositions
PM2.5 was collected in Suncheon during the summer (June 2, 2023, to June 11, 2023) and winter (January 15, 2024, to January 21, 2024). The analysis results of the carbonaceous components (OC, EC), secondary ionic components (NH4+, NO3−, SO42−), QDTT-OP, and VOCs in the collected PM2.5 are summarized in Table 2. The time series of temperature, relative humidity, PM2.5, OC, QDTT-OPv, and QDTT-OPm in summer and winter are shown in Fig. 3. Meteorological conditions were observed during the measurement period, with average summer temperatures and relative humidity of 23.6°C and 58.3%, respectively, and winter values of 6.5°C and 62.9%. The average PM2.5 concentrations in summer and winter were 19.54 (3.88–53.35) µg/m3 and 18.76 (1.32–90.76) µg/m3, respectively, with slightly higher concentrations in summer. Generally, summer PM2.5 is formed through active photochemical oxidation processes involving high concentrations of gaseous precursors, while winter PM2.5 is produced by increased air pollutants from heating and stable atmospheric conditions (Cheng et al., 2016; Elser et al., 2016). As shown in Figs. 5a and 5g, summer PM2.5 concentrations tended to increase steadily in the afternoon until sunset due to photochemical reactions, whereas winter concentrations peaked during morning and evening rush hours. PM2.5 levels gradually increased from 3 AM on January 17 to 9 PM the following day and then decreased from the early morning of January 18. High relative humidity generally plays a crucial role in the formation of secondary pollutants through heterogeneous reaction processes in the atmosphere (Yu et al., 2018; Liu et al., 2015). Ma et al. (2017) observed that higher PM2.5 concentrations in winter were associated with increased relative humidity and decreased wind speed. On January 17, the average wind speed and relative humidity were 0.5 m/s and 71%, respectively, suggesting that low wind speeds led to the stagnation and increase of PM2.5.
Table 2
Overall average results of major chemical compounds in PM2.5, BTEX, and QDTT-OP
Compounds | unit | summer | winter |
average | S.D | average | S.D |
PM2.5 | µg/m3 | 19.54 | 10.63 | 18.76 | 19.04 |
OC | µg/m3 | 5.10 | 1.72 | 4.49 | 3.84 |
EC | µg/m3 | 0.20 | 0.17 | 0.82 | 0.71 |
NO3− | µg/m3 | 1.41 | 0.82 | 2.97 | 2.97 |
SO42− | µg/m3 | 3.80 | 1.48 | 1.91 | 1.93 |
NH4+ | µg/m3 | 1.49 | 0.50 | 1.43 | 1.21 |
QDTT-OPv | µM/m3 | 0.12 | 0.04 | 0.09 | 0.06 |
Benzene | ppb | 0.30 | 0.05 | 0.40 | 0.22 |
Toluene | ppb | 1.47 | 0.52 | 0.67 | 0.49 |
Ethylbenzene | ppb | 0.25 | 0.09 | 0.74 | 0.66 |
m&p-Xylene | ppb | 0.45 | 0.15 | 2.18 | 1.91 |
Styrene | ppb | 0.10 | 0.03 | 0.01 | 0.01 |
o-Xylene | ppb | 0.18 | 0.06 | 0.79 | 0.79 |
During the measurement period, OC and EC showed contrasting results between winter and summer (Figs. 5d and 5j). The average OC concentration was approximately 1.1 times higher in summer (5.10 µg/m3) than in winter (4.49 µg/m3), while the average EC concentration was about 4.0 times lower in summer (0.20 µg/m3) than in winter (0.82 µg/m3). Summer OC was highest during the evening hours (18:00–21:00), while winter OC peaked during the morning hours (06:00–09:00). EC increased during rush hours in summer and peaked in the morning and afternoon in winter. EC acts as a primary emission tracer from incomplete combustion processes based on fossil fuels used in heating activities, vehicles, and industrial processes, while OC can be emitted directly or formed secondarily (Wu et al., 2016; Karanasiou et al., 2015). The correlation between OC and EC was higher in winter (r2 = 0.77) than in summer (r2 = 0.45). The lower r2 value between OC and EC in summer suggests that other formation processes, such as secondary OC formation, contributed to the increase in OC concentration. In contrast, the higher r2 value in winter indicates that significant amounts of measured OC and EC were closely related to primary combustion sources. The OC/EC ratio is used to infer sources of carbonaceous aerosols and to determine the contributions of POA and SOA, with higher OC/EC ratios indicating greater SOA formation (Paraskevopoulou et al., 2014). The observed OC/EC ratio was higher in summer than in winter and increased gradually in the afternoon during summer, while it decreased in the afternoon during winter. The high OC/EC ratio in summer indicates that more secondary organic aerosols were generated due to active photochemical reactions, which is characteristic of the season.
In the atmosphere, NH4+, NO3−, and SO42− in PM2.5 are secondary products formed from the homogeneous and heterogeneous oxidation of large amounts of gaseous precursors (NH3, NOx, SO2) emitted from industrial complexes (Lin, 2002; Chow et al., 1996). The average concentrations of NH4+, NO3−, and SO42− in summer were 1.49 (0.63–2.39), 1.41 (0.39–3.99), and 3.80 (1.45–7.26) µg/m3, respectively, while in winter, they were 1.43 (0.04–5.32), 2.97 (0.07–12.24), and 1.91 (0.05–10.51) µg/m3, respectively. NH4+ showed no seasonal difference in this study, while NO3− and SO42− exhibited contrasting results. On January 17, 2024, the peak concentrations of secondary ionic components were observed, indicating their significant contribution to the high PM2.5 levels.
3.2 Stagnation and Long-Range Transport Period on January 17, 2024
Currently, the Korea Air Quality Forecasting Center operates a high-performance air quality forecasting modeling system at the level of a supercomputer, with numerical models run four times a day (www.airkorea.or.kr). The air quality forecasting modeling system is divided into the modeling process and the processing and refining of the driven model data. The modeling process is subdivided into the weather model (Weather Research and Forecasting, WRF), the emissions model (Sparse Matrix Operator Kernel Emissions, SMOKE), and the air quality model (Community Multiscale Air Quality, CMAQ). According to the results, on January 17, high concentrations of fine particulate matter in the measurement area were analyzed to be due to the stagnation of residual fine particulate matter from the previous day along with the influx of externally generated fine particulate matter (Fig. 4).
3.3 Results of benzene, toluene, ethylbenzene, and xylene
During the summer observation period, the concentrations were as follows: benzene 0.30 ppb, toluene 1.47 ppb, ethylbenzene 0.25 ppb, m&p-xylene 0.45 ppb, and o-xylene 0.18 ppb. When comparing the total concentrations of BTEX between seasons, the summer period recorded 2.63 ppb, while the winter period recorded 4.79 ppb, indicating higher concentrations in winter. Generally, the concentrations of VOCs, including BTEX, are influenced by seasonal temperature and atmospheric boundary layer conditions. In detail, the atmospheric boundary layer expands during the unstable summer atmosphere due to higher temperatures, leading to the diffusion of BTEX, which results in lower BTEX concentrations despite the same emission amount (Latif et al., 2019; Ghaffari et al., 2021). Conversely, higher temperatures in summer enhance the evaporation of BTEX, increasing its atmospheric concentration (Yurdakul et al., 2018; Ghaffari et al., 2021). Therefore, to understand the seasonal differences in BTEX concentrations, it is essential to examine the emission characteristics of the major BTEX sources.
To investigate the emission characteristics of BTEX, the ratios of m,p-xylene/ethylbenzene (X/E) and toluene/benzene (T/B) as well as the correlation between toluene and ethylbenzene (T/E) were analyzed. The X/E ratio provides insight into the photochemical aging of BTEX. Specifically, m,p-xylene and ethylbenzene, when emitted from the same source, undergo different degradation rates due to varying oxidation rates by OH radicals in the atmosphere, leading to a decreasing X/E ratio over time (Monod et al., 2001). An X/E ratio of 2 or more indicates local emissions (Nelson and Quigley, 1983). The T/B ratio is used to infer emission sources: T/B ≤ 1 suggests biomass burning, 0.1–5 indicates coal combustion, 0.5–10 suggests vehicle emissions from gasoline and diesel, and 1–10,000 implies industrial processes and solvent use (Zhang et al., 2016). The T/E correlation helps trace emission sources, with high correlations (> 0.9) suggesting combustion sources such as vehicles and heating, and low correlations (< 0.5) indicating evaporative sources like solvents (Monod et al., 2001; Sigsby et al., 1987). The X/E ratio in the study area averaged 2.53 in summer and 4.30 in winter, indicating that BTEX emissions were influenced by local sources rather than external transport. The T/B and T/E correlations are presented in Fig. 6. The T/B ratio was used to estimate BTEX sources under the assumption that toluene and benzene share the same sources. This was validated by their correlation, which varied over time. During the summer observation period, the correlation coefficient (r²) between toluene and benzene was significant (> 0.5) from 09:00 to 24:00, enabling source estimation. The average T/B ratio from 09:00 to 21:00 was 6.16, suggesting BTEX emissions from coal, diesel, gasoline combustion, and industrial processes during these hours. In winter, significant correlations were observed from 03:00 to 24:00, with a T/B ratio of 1.77, indicating identical local sources.
As shown in Fig. 6 (c & d), the T/E ratio analysis during the summer observation period revealed two distinct groups: from 21:00 to 09:00, the average r² was 0.64, and from 09:00 to 21:00, the average r² was 0.02. This pattern was also observed in winter, with r² values of 0.92 from 15:00 to 09:00 and 0.32 from 09:00 to 15:00. A high T/E correlation suggests emissions from combustion activities, while a low correlation indicates evaporative emissions from solvents. These findings suggest that the study area is influenced by a mix of combustion and evaporative sources. Specifically, the summer daytime (09:00–21:00) and winter daytime (09:00–15:00) periods were more influenced by solvent evaporation, while summer nighttime (21:00–09:00) and winter nighttime (15:00–09:00) were more influenced by combustion activities. The lower T/B ratio in winter (1.67) compared to summer (4.04) suggests that winter heating activities primarily involved low T/B fuels such as coal or LNG (Alsbou and Omari, 2020; Mokammel et al., 2022; Zhang et al., 2016).
The health risks associated with BTEX exposure were evaluated using the LTCR for benzene and ethylbenzene and the HQ for non-carcinogenic risk. The ECi was calculated based on a 24-hour exposure period to align with the 3-hour BTEX measurement intervals. The results are presented in Fig. 6 (e & f). The cumulative HQ for BTEX HI was 0.017 in summer and 0.049 in winter, indicating a 2.8 times higher non-carcinogenic risk in winter compared to summer. However, all HQ values were below 1, suggesting acceptable risk levels (Dehghani et al., 2018; Mokammel et al., 2022; Rostami et al., 2021). The LTCR for benzene and ethylbenzene was 2.07×10− 6 in summer and 2.81×10− 6 in winter, showing a 1.3 times higher carcinogenic risk in winter. According to WHO, LTCR values between 1×10− 5 and 1×10− 6 are considered acceptable for humans, while USEPA recommends values below 1×10− 6 (Dehghani et al., 2018; Delikhoon et al., 2018; Nabizadeh et al., 2020). The study area's LTCR met WHO standards but exceeded USEPA recommendations. Benzene was the primary contributor to LTCR, with r² values of 0.99 for both seasons, indicating the need to manage benzene to reduce carcinogenic risk. Despite benzene's high HQ, the overall non-carcinogenic risk (HI) was more closely correlated with toluene, ethylbenzene, and xylene, especially xylene (r² > 0.94).
Human health risks were also assessed based on work hours, distinguishing between daytime (09:00–18:00) and nighttime (21:00–06:00). The study area, consisting of small-scale businesses, showed higher Hi and LTCR values during the daytime for both summer and winter observation periods. Specifically, daytime HI and LTCR were 1.17 and 1.10 times higher, respectively, in summer, and 1.35 and 1.44 times higher, respectively, in winter compared to nighttime. This indicates higher carcinogenic and non-carcinogenic risks during the day, likely due to solvent evaporation. The results suggest the need to manage BTEX emissions from solvent evaporation and to prioritize benzene for carcinogenic risk and xylene for non-carcinogenic risk.
3.4 Association with Oxidative Potential
The QDTT-OP was quantified to evaluate the health risks of PM2.5 observed during the measurement period. The average concentration of oxidative potential per unit volume (QDTT-OPv) was higher in summer (0.12 µM/m³) than in winter (0.09 µM/m³). The highest QDTT-OPv concentration in summer was observed on June 5, when OC levels were high, while in winter, it was observed on January 17, during a high PM2.5 event. The diurnal variation of QDTT-OPv showed a decreasing trend in the afternoon during summer, contrary to the high PM2.5 concentrations due to active photochemical reactions (Fig. 5b). This result differs from previous studies showing a high correlation between DTT-OP and PM2.5 mass concentration (Kurihara et al., 2022; Janssen et al., 2014). Oxidative potential per unit mass of PM2.5 considers the intrinsic properties of PM2.5 that affect its oxidative capacity. Thus, the quantified QDTT-OPv was divided by the PM2.5 mass concentration to calculate the QDTT-OPm (µM/µg) and shown in Figs. 5b and 5h. As a result, the diurnal variations of QDTT-OPv and QDTT-OPm in summer were different. Additionally, the QDTT-OPv and OC/PM2.5 ratio in summer both peaked at 9 PM, showing a very similar diurnal trend and the highest correlation (r²=0.97) among major chemical components in PM2.5 (Fig. 7a). As the contribution of OC to PM2.5 increased, the oxidative potential also increased. Borlaza et al. (2018) observed that the oxidative potential increased with OC due to the influence of secondary organic aerosols in summer, consistent with the results of this study. These findings suggest that OC significantly impacts the oxidative potential of PM2.5 in summer and can be used to assess the health risks of PM2.5 from industrial complexes. Consequently, when evaluating the concentration of reactive oxygen species per unit mass, PM2.5 health risks were highest at night, rather than during rush hours or the active photochemical reaction periods in the afternoon. Saffari et al. (2013) also found a strong correlation between oxidative potential and carbonaceous components (OC, EC), especially with primary organic matter from vehicle emissions in winter. Furthermore, they confirmed the contributions of both primary and secondary sources to the oxidative potential of PM2.5. Similarly, this study found a good correlation between winter QDTT-OPm and unit mass of OC, EC, and NO3− (r²=0.80) (Fig. 7b). Recent studies have revealed that combustion-related sources (primarily oil, coal, and biomass burning) can directly emit NO3− (Chen et al., 2010; Cui et al., 2017). Unlike in summer, the diurnal variation in winter showed increased concentrations during rush hours, indicating the influence of vehicle emissions and biomass burning, consistent with the chemical component results. Overall, these results indicate that oxidative potential is more influenced by chemical composition than physical characteristics. Additionally, from a health risk perspective, regulating only PM2.5 concentration may be insufficient.
According to Song et al. (2023), VOCs, gaseous precursors of secondary organic aerosols, showed a positive correlation with QDTT-OPv, indicating direct and indirect impacts. Similarly, this study found a positive correlation between observed VOCs and QDTT-OPv, with winter showing relatively better correlations than summer, especially with benzene (r²=0.41). These results suggest that VOCs during the study period acted as major precursors of secondary organic aerosols, directly and indirectly influencing QDTT-OPv.