Effect modification of consecutive high concentration days on the association between fine particulate matter and mortality: a multi-city study in Korea

Hyungryul Lim; Sanghyuk Bae; Jonghyuk Choi; Kyung-Hwa Choi; Hyun-Joo Bae; Soontae Kim; Mina Ha; Ho-Jang Kwon

doi:10.4178/epih.e2022052

Articles

Page Path: HOME > Epidemiol Health > Volume 44; 2022 > Article

Original Article Effect modification of consecutive high concentration days on the association between fine particulate matter and mortality: a multi-city study in Korea: Hyungryul Lim¹, Sanghyuk Bae², Jonghyuk Choi¹, Kyung-Hwa Choi¹, Hyun-Joo Bae³, Soontae Kim⁴, Mina Ha¹, Ho-Jang Kwon¹; Epidemiol Health 2022;44:e2022052.
DOI: https://doi.org/10.4178/epih.e2022052
Published online: June 9, 2022

3,590 Views
272 Download

¹Department of Preventive Medicine, Dankook University College of Medicine, Cheonan, Korea

²Department of Preventive Medicine, College of Medicine, The Catholic University of Korea, Seoul, Korea

³Korea Environment Institute, Sejong, Korea

⁴Department of Environmental and Safety Engineering, Ajou University, Suwon, Korea

Correspondence: Ho-Jang Kwon Department of Preventive Medicine, Dankook University College of Medicine, 119 Dandae-ro, Dongnam-gu, Cheonan 31116, Korea E-mail: hojangkwon@gmail.com

• Received: February 18, 2022 • Accepted: June 9, 2022

This is an open-access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Full Article

Download PDF

Abstract
GRAPHICAL ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
SUPPLEMENTARY MATERIALS
NOTES
ACKNOWLEDGEMENTS
REFERENCES

Abstract

OBJECTIVES
Although there is substantial evidence for the short-term effect of fine particulate matter (PM_2.5) on daily mortality, few epidemiological studies have explored the effect of prolonged continuous exposure to high concentrations of PM_2.5. This study investigated how the magnitude of the mortality effect of PM_2.5 exposure is modified by persistent exposure to high PM_2.5 concentrations.
METHODS
We analyzed data on the daily mortality count, simulated daily PM_2.5 level, mean daily temperature, and relative humidity level from 7 metropolitan cities from 2006 to 2019. Generalized additive models (GAMs) with quasi-Poisson distribution and random-effects meta-analyses were used to pool city-specific effects. To investigate the effect modification of continuous exposure to prolonged high concentrations, we applied categorical consecutive-day variables to the GAMs as effect modification terms for PM_2.5.
RESULTS
The mortality risk increased by 0.33% (95% confidence interval [CI], 0.16 to 0.50), 0.47% (95% CI, -0.09 to 1.04), and 0.26% (95% CI, -0.08 to 0.60) for all-cause, respiratory, and cardiovascular diseases, respectively, with a 10 μg/m³ increase in PM_2.5 concentration. The risk of all-cause mortality per 10 μg/m³ increase in PM_2.5 on the first and fourth consecutive days significantly increased by 0.63% (95% CI, 0.20 to 1.06) and 0.36% (95% CI, 0.01 to 0.70), respectively.
CONCLUSIONS
We found increased risks of all-cause, respiratory, and cardiovascular mortality related to daily PM_2.5 exposure on the day when exposure to high PM_2.5 concentrations began and when exposure persisted for more than 4 days with concentrations of ≥35 μg/m³. Persistently high PM_2.5 exposure had a stronger effect on seniors.
Keywords: Air pollution; Particulate matter; Epidemiology; Mortality; Time-series analysis

GRAPHICAL ABSTRACT

INTRODUCTION

The background level of fine particulate matter (PM_2.5) concentration is much higher in Korea than in North America and Western Europe [1]. According to the Korean Ministry of Environment, the average annual concentration of PM_2.5 in all regions of Korea in 2019 was 23 μg/m³. In 2011, the Korean government established new air quality guidelines for PM_2.5, which were strengthened in March 2018 to advise a daily average of 35 μg/m³ instead of 50 μg/m³ and an annual average of 15 μg/m³ instead of 25 μg/m³ [2]. Consequently, the public has been notified on days when this standard is exceeded. Moreover, when the hourly average concentration exceeds 75 μg/m³ or 150 μg/m³ for more than 2 hours, an advisory or alert is issued regionally, which has increased public concern about the health effects of PM_2.5.

In epidemiological studies on the health effects of air pollution, exposure is generally classified as both short-term and long-term. The former refers to exposure over a period ranging from a few hours to a month (mostly within the same day to 1 week), and the latter refers to exposure over a period of more than 1 month to several years [3].

Despite its heterogeneity across different cities worldwide, several multi-city time-series studies of the short-term effects of ambient PM_2.5 exposure on daily mortality have consistently observed statistically significant effects [4-7], and evidence for causality has been reported in recent studies [8,9]. Although the background PM_2.5 concentration in Korea has shown a decreasing trend since these levels started to be measured in 2015 [10], prolonged exposure to a high concentration of PM_2.5 has frequently been observed. The average duration of exposure to a PM_2.5 concentration of 35 μg/m³ or above in the first quarter of the year in Korea reportedly increased from 16.2 hours in 2015 to 26.5 hours in 2018, and the average 1-hour concentration also increased slightly from 95.4 μg/m³ in 2015 to 102 μg/m³ in 2018 [11].

Thus, we reviewed the existing evidence on whether the mortality effect of short-term exposure to high concentrations of PM_2.5 differs from that of short-term exposure to non-high concentrations of PM_2.5, and, if so, whether the effect differed according to the length of exposure. However, relatively few observational studies have addressed the durational effects of continuous exposure to prolonged high concentrations of particulate matter (PM). Two existing previous studies found additional positive durational effects on daily mortality over several days [12,13].

We aimed to investigate the short-term effect of ambient PM_2.5 on daily mortality in Korea and assess the effect of prolonged exposure to a high concentration of PM_2.5 by examining whether the effect sizes were modified based on the duration of exposure.

Study area and population: The study areas were Seoul, the capital city of Korea, and 6 other metropolitan areas: Busan, Daegu, Incheon, Gwangju, Daejeon, and Ulsan (Supplementary Material 1). The study population from these 7 areas was approximately 22.6 million in 2019, covering approximately 43.7% of the nation’s total population. The study was conducted from January 1, 2006, to December 31, 2019, for a total of 5,113 days.
Mortality data: We obtained cause of death statistics data from the MicroData Integrated Service (https://mdis.kostat.go.kr/) of Statistics Korea, which publishes open-source mortality data available to the public for research purposes. We included daily non-accidental all-cause deaths (International Classification of Diseases, 10th edition [ICD-10], codes A00 to R99), deaths from respiratory disease (ICD-10 codes J00 to J98), and deaths from cardiovascular disease (I00 to I99) in 7 major cities from January 1, 2006, to December 31, 2019.
Fine particulate matter mass concentration data: In Korea, ambient PM_2.5 levels have been measured through the national air pollution monitoring network since 2015. Therefore, we used simulated data that considered weather conditions, anthropogenic and biogenic emissions, and chemical transport to estimate PM_2.5 concentrations extending back to 2006 [14,15], as has been done in several recent epidemiological studies of air pollution in Korea [16,17].; We used the Community Multiscale Air Quality (CMAQ) system, version 4.7.1, with the AERO5 aerosol module and Statewide Air Pollution Research Center 99. Weather simulations were performed using the Weather Research and Forecasting model, version 3.3.1, with the National Center for Environmental Protection final data as the initial field. The Meteorology Interface Processor version 3.6 was used to prepare the CMAQ-ready meteorological inputs. The Clean Air Policy Support System 2010, which is the Korean national emissions inventory, was processed through the Sparse Matrix Operator Kernel Emission, version 3.1, to estimate anthropogenic emissions, and biogenic emissions were estimated using the Model of Emissions of Gases and Aerosols from Nature. We applied a 27-km (covering Northeast Asia, including China, Japan, and Korea) and a 9-km modeling domain. The 27-km modeling domain simulation values were applied to the boundary condition of the 9-km modeling domain. The CMAQ system simulated gridded hourly PM_2.5 concentrations from the 9-km modeling domain for each of the 7 cities, which were resampled and averaged to allocate city-specific concentrations of PM_2.5. Subsequently, we calculated the daily average PM_2.5 concentration for each of the cities from January 1, 2006, to December 31, 2019.
Meteorological data for covariates: We collected daily average temperature and relative humidity data (1 measurement point per city) for 7 major cities from January 1, 2006 to December 31, 2019, from an open-source dataset published by the National Climate Data Center of the Korea Meteorological Administration (https://data.kma.go.kr/).
Statistical analysis: We first evaluated the association between daily PM_2.5 concentration and daily non-accidental all-cause, respiratory disease, and cardiovascular disease mortality using time-series analyses as basic procedures. A generalized additive model (GAM) with a quasi-Poisson distribution was used as a statistical model, and the commonly applied model equation was as follows:; Where E(Y_t,c) is the expected death count on day t in city c; PM_{2.5, t, l, c} is the average 24-hour mean PM_2.5 concentration (μg/m³) on day t applied lag l (0 to 6 days before, and the 2-day to 7-day moving average, which means the average over the current and previous day and the current day to before the sixth day, respectively) in city c; β_{PM 2.5, l, c} is the corresponding coefficient; temp_{t, c} is the average 24-hour mean temperature (°C) on day t and day t-1 in city c; rh_{t, c} is the mean relative humidity (%) on day t in city c, time_t is a continuous variable-processed of day t to adjust the long-term trend, and dow_t is a categorical day-of-week and holiday variable of day t, which has 4 levels (weekday, Sunday and holiday, the day after Sunday and holiday, and Saturday) [14], where s represents the smooth spline function, df represents the degrees of freedom applied to each function, and η represents the coefficients of the dummy variables of dow_t. For temperature, the 2-day moving average (lag 0-1) was determined since the generalized cross-validation (GCV) score was the lowest among the various lag-applied models; dfs for temperature, relative humidity, and time trend were determined based on the consensus model methodology of a recent large-scale multi-city time series study [4]. For the df of the long-term trend, we also constructed GAMs applying 1-14 dfs per year. The GCV score was the lowest when 7 df was applied in the all-cause mortality models (Supplementary Material 2).; The coefficients (β_{PM 2.5}) obtained from the 7 cities were pooled with the same exposure lag structure using random-effects meta-analyses, and the mortality risk percentage change per 10 μg/m³ increase in PM_2.5 concentration was estimated using the following equation:; For the second step, when the average 24-hour mean PM_2.5 concentration was ≥ 35 μg/m³, it was classified as a high-concentration day. To investigate whether the effect size of PM_2.5 changed during prolonged exposure to a high PM_2.5 concentration, we added a 6-stratum consecutive day variable (no high-concentration days, the first to fourth day, and the fifth day and beyond; a coding example is shown in Supplementary Material 3) and its effect modification term to PM_2.5 concentration in the quasi-Poisson GAM was as follows:; Where con_day_{t, l, c, k} are 6 strata dummy variables corresponding to whether day t is not a high-concentration day (k=1), the first (k=2) to fourth (k=5) high consecutive days, and the fifth day or more (k=6) applied lagged exposure (l) in city c, respectively. β_{1, l, c} to β_{6, l, c} are the effects of PM_2.5 in the corresponding high consecutive day strata. The rest of the variables are the same as those in the basic GAM. These models allow a different mortality effect within the predefined consecutive day strata [18,19]. We excluded the durational effect itself from the models since we did not observe a difference in model fitness based on whether an independent effect was applied. Thereafter, we pooled the PM_2.5 coefficients paired for each lag structure and level of consecutive day variables using random-effects meta-analyses and estimated the mortality risk percentage change using the process described above.; To investigate the possible effect modification by age group, the above analyses were performed for 2 stratified age groups: 20-64 years old and ≥ 65 years old. For sensitivity analyses, we designed 2-pollutant models with daily ozone and nitrogen dioxide concentrations to measure the robustness of the effects of PM_2.5 and their patterns within a prolonged duration of exposure to a high concentration of PM_2.5. As a final step, in order to examine whether the effect modification of consecutive high-concentration days on mortality was caused by different background concentrations on each consecutive day, we also constructed an effect modification model in the same manner as described above based on a 10 μg/m³ unit daily background concentration interval (< 10 to ≥ 60 μg/m³) and compared the effect size of PM_2.5 in the 2 models.; All dataset setups and statistical analyses were performed using R version 3.6.3 (Foundation for Statistical Computing, Vienna, Austria). We used the “mgcv” package for GAM modeling and the “metafor” package for random-effects meta-analyses in R software. The statistical significance level for the 2-tailed tests was set at 0.05.
Ethics statement: The Institutional Review Board of Dankook University, Korea, exempted this study from review since it exclusively used anonymous secondary data (IRB No. DKU 2021-03-042).

RESULTS

Table 1 summarizes the non-traumatic all-cause, respiratory, and cardiovascular mortality rates in 7 major cities from 2006 to 2019. In the study area, the total numbers of non-traumatic all-cause, respiratory, and cardiovascular deaths were 1,314,829 (daily mean, 257.2), 120,094 (daily mean, 23.5), and 326,034 (daily mean, 63.8), respectively.

The mean concentration of PM_2.5 throughout the study period ranged from 23.0 μg/m³ to 27.4 μg/m³, and the proportion of high PM_2.5 concentration (an average 24-hour mean PM_2.5 concentration of ≥ 35 μg/m³) in all cities during the study period ranged from 13.6% to 22.9%. Five or more consecutive high-concentration days accounted for 0.8-2.7% (Table 2). A declining trend was observed over the past 14 years. However, the annual average concentration in all 7 cities exceeded 20 μg/m³ in 2019, surpassing the annual governmental air quality guideline of 15 μg/m³ (Supplementary Material 4). Although the number of high-concentration days decreased every year up to 2019 (Supplementary Material 5), prolonged episodes of high concentrations did not disappear, and there were still consecutive days lasting more than a week in Seoul, Incheon, Gwangju, and Daejeon in 2019 (Figure 1).

The percentages of risk changes per 10 μg/m³ increase in PM_2.5 concentration for all-cause, respiratory, and cardiovascular mortality in the lag 0-1 models were 0.33% (95% confidence interval [CI], 0.16 to 0.50), 0.47% (95% CI, -0.09 to 1.04), and 0.26% (95% CI, -0.08 to 0.60), respectively. The largest effects per 10 μg/m³ increase in all-cause, respiratory, and cardiovascular mortality were 0.33% (95% CI, 0.18 to 0.47, lag 0 model; more precise than lag 0-1), 0.62% (95% CI, -0.01 to 1.26, lag 0-2 model), and 0.39% (95% CI, 0.10 to 0.68, lag 0 model), respectively (Table 3).

In the models that applied the effect modification of consecutive high-concentration days, applying lag 0-1 exposure, we found increases in non-accidental all-cause mortality increases per 10 μg/m³ of 0.63% (95% CI, 0.20 to 1.06) on the first day and 0.36% (95% CI, 0.01 to 0.70) on the fourth day, which were higher than the effect size we found in the basic model. This pattern was similar for respiratory and cardiovascular mortality (Table 4). Considering the lag structures, we found higher effect sizes and clearer modified patterns in the lag 0-1 model (Supplementary Material 6).

In sensitivity analyses, the PM_2.5 mortality effects were attenuated when we analyzed the 2-pollutant models that included daily ozone and nitrogen dioxide concentration. Nevertheless, the effect modification patterns of consecutive high-concentration days persisted (Supplementary Material 7).

Although not statistically significant, we observed patterns of an increased effect on the fourth consecutive day in the 20-64 years age group for all-cause and respiratory mortality (Figure 2, Supplementary Material 8). The effect modification pattern was more pronounced in the ≥ 65 years age group (Figure 2, Supplementary Material 9).

In the sensitivity analysis with the effect modification according to the background concentration of PM_2.5, we found that the effect sizes in all background concentration intervals were smaller than the effects of PM_2.5 on the first and fourth consecutive days for all-cause mortality in the lag 0-1 model (Supplementary Material 10).

DISCUSSION

We found a short-term positive effect of PM_2.5 on mortality in this multi-city time-series study of 7 major Korean cities from 2006 to 2019. With continuous exposure to prolonged high concentrations of PM_2.5 exceeding the daily mean of 35 μg/m³, the effects of PM_2.5 on daily all-cause, respiratory, and cardiovascular mortality were higher on the first and fourth consecutive high-concentration days. In addition, this effect was mainly observed among the elderly (aged 65 years or older).

There have been several studies from Korea on the short-term effects of PM_2.5 on mortality compared to the long-term effects; however, the effect sizes vary greatly depending on the spatiotemporal background of each study, the applied statistical model, and the exposure assessment [20]. Moreover, as of this study, only 2 multi-city designs have been published on the effect of particulate matter with a diameter of 10 microns or less (PM₁₀) on mortality [4,21].

Two previous studies reported the short-term effects of PM_2.5 on daily mortality for each city. Jung et al. [22] found that the total daily mortality risk of individuals aged ≥ 60 years increased by 0.36% per 10 μg/m³ increase of PM_2.5 from 2000 to 2012 in Seoul. Kim et al. [23] found that the total daily mortality risk increased with a 10 μg/m³ increase of PM_2.5 by 0.34%, 1.18%, and 0.43% in Seoul, Busan, and Incheon, respectively, from 2006 to 2012. We observed somewhat different effect sizes when compared to previous estimates. There were differences in the study period, which we expanded to 2019. The number of media reports in Korea of issues related to PM_2.5 exposure has increased rapidly since 2012 [24], which was reflected in the increased rate of the everyday use of face masks due to improved risk awareness [25].

Two previous studies focused on the short-term effects of prolonged continuous exposure to a high concentration of PM. The first study observed an increased risk of cardiovascular and respiratory mortality after continuous exposure to high concentrations of PM in Beijing, China. Using GAM, the PM_2.5 concentration variable was not added, and the categorical variable instead specified the duration of exposure to the high concentration of PM_2.5 based on daily averages of 75 μg/m³, 85 μg/m³, 105 μg/m³, and 115 μg/m³, which were applied to the models. On the ninth consecutive day of exposure to a high concentration of PM_2.5 of 105 μg/m³ or more, a 53% increase in the risk of cardiovascular mortality was reported among outdoor workers [12].

The second study observed a short-term durational effect of prolonged exposure to high concentrations of PM₁₀ in 28 cities in China, Japan, and Korea. Using quasi-Poisson GAM, the effects of PM₁₀ concentration and duration (number of consecutive days of exposure to 75 μg/m³ or more PM₁₀) were separated. In Korea, the estimated increases in risk for each additional consecutive day of exposure to a high concentration of PM₁₀ were 0.48% for non-traumatic all-cause mortality, 0.48% for cardiovascular mortality, and 1.13% for respiratory mortality [13].

Based on these 2 previous studies, we thought that it would be necessary to investigate whether the short-term mortality effect of PM_2.5 was modified to a greater degree when exposure to a high concentration of PM persisted. To the best of our knowledge, this is the first study to report the effect modification of the short-term effect of PM_2.5 after continued exposure to high concentrations.

A novel finding of this study was that the short-term effects of PM_2.5 on daily mortality on the first and fourth consecutive days of exposure to high concentrations of PM_2.5 were greater than the estimated effect for the entire period. From a short-term perspective, it was notable that the risk did not increase linearly as continuous exposure to high concentrations of PM_2.5 increased; instead, the risk increased starting at the first transition to a day of high concentrations and at a later point, such as, for example, the fourth consecutive day of exposure.

We did not observe any reported biological mechanism related to effect modification. However, mortality displacement has been observed in previous studies on short-term air pollution epidemiology. This means that air pollution can cause mortality in frail individuals several days or weeks sooner than if they had not been exposed to air pollution. Thus, the mortality effect would be lower than expected after an initial risk increase [26,27]. In fact, we found that the effect on the first consecutive day of exposure to a high concentration of PM_2.5 was remarkable among those aged ≥ 65 years. Among those aged 20-64 years, we did not observe a significant change in the effect due to the relatively low number of deaths. Nevertheless, even in the less vulnerable group, the point estimates were the greatest on the fourth day for all-cause and respiratory mortality. We speculate that another mortality displacement due to exposure to high concentrations of PM_2.5 was revealed in this study. Furthermore, when high concentrations persisted for several days, relatively less frail individuals could also be affected by air pollution.

There is still insufficient understanding of intermittent episodes of high concentrations of PM_2.5 in Korea, as well as spatial differences in the contribution sources and components of PM_2.5 [11]. Increases in sulfate, nitrate, and ammonium ion concentrations have recently been more pronounced than organic carbon, inorganic carbon, and heavy metals in most consecutive high-concentration episodes without yellow dust storms in the western, central, and southeastern parts of the Korean Peninsula [11,28].

Despite the findings reported in Korea, evidence of the effects of sulfate, nitrate, and ammonia ions is still insufficient and inconsistent [29-32]. However, a recent meta-analysis found a statistically significant short-term increase in the risk of cardiovascular mortality due to nitrate and sulfate exposure and respiratory hospitalization due to nitrate exposure [33]. The effect of PM_2.5 on daily mortality may increase with consecutive days of exposure to high concentrations of PM_2.5 due to variations in the PM composition, particularly increases in nitrates and sulfates.

In the sensitivity analysis, we estimated the relative risks for all-cause mortality across 10 μg/m³ intervals of background PM_2.5 concentration. A recent study in China showed that the relative risks of total and cardiovascular mortality increased as the background concentration increased [34]. We observed higher effect estimates on the first and fourth consecutive days than those derived from high PM_2.5 concentration intervals of > 40 μg/m³ (Supplementary Material 10); therefore, the effect modification we observed was likely not due to the background concentration effect on the duration of exposure to high concentrations of PM_2.5.

This study has several limitations. First, there may have been misclassifications in exposure assessment. In other words, local variations in the PM_2.5 concentration within a single metropolitan area could not be reflected. In addition, the specific causes of regional and temporal variations in the effect of PM_2.5 exposure could not be explained in our study. Future studies are needed to investigate the effect of applying the components of PM_2.5 to the models as proposed in our study.

Despite these shortcomings, this study has several strengths. Among time series studies of the effect of PM_2.5 exposure on short-term mortality in Korea, this study examined the longest time period at 14 years and was the first multi-city study of this type, including more than half of the entire Korean population. Moreover, by focusing on changes in the effect size of PM_2.5 after persistent exposure to high concentrations, which is a challenge in East Asia, we found that the mortality effect of PM_2.5 exposure may increase on a short-term basis even more during periods of exposure to high concentrations. Further epidemiological studies in other East Asian regions are needed.

We found a greater effect on daily mortality for the day when the duration of exposure to a high concentration of PM_2.5 began and when exposure lasted for approximately 4 days. The elderly may be more affected by persistent exposure to high concentrations of PM_2.5. Health authorities should encourage seniors or frail individuals to refrain from outdoor activity and wear a mask on days when there is a high concentration of PM_2.5 that is expected to persist for several days. When the causes of episodes of high concentrations of PM_2.5 and their components are identified, our results can be used as scientific evidence to support public risk communication and policy-making.

SUPPLEMENTARY MATERIALS

Supplementary materials are available at https://www.e-epih.org/.

Supplementary Material 1.

The map of seven major cities (Seoul, Busan, Daegu, Incheon, Gwangju, Daejeon, and Ulsan) in Korea.

epih-44-e2022052-suppl1.docx

Supplementary Material 2.

Sensitivity analysis for the degree of freedom of time trend in generalized additive model for all-cause, respiratory, and cardiovascular mortality

epih-44-e2022052-suppl2.docx

Supplementary Material 3.

An example of setting consecutive day variables in the effect modification models according to the daily mean concentration of PM_2.5 (reference level=35 μg/m³).

epih-44-e2022052-suppl3.docx

Supplementary Material 4.

Daily mean PM_2.5 concentration in seven major cities in Korea from 2006 to 2019.

epih-44-e2022052-suppl4.docx

Supplementary Material 5.

Days and proportion of the daily mean PM_2.5 concentration of 35 μg/m³ or more in seven major cities in Korea from 2006 to 2019 (5,113 days).

epih-44-e2022052-suppl5.docx

Supplementary Material 6.

Mortality effects of short-term ambient PM_2.5 exposure and their modification by consecutive days of exposure to high concentrations of PM_2.5 in seven major cities in Korea from 2006 to 2019.

epih-44-e2022052-suppl6.docx

Supplementary Material 7.

Mortality effects of daily ambient PM_2.5 exposure and their modification by consecutive days of exposure to high concentrations of PM_2.5 in seven major cities in Korea from 2006 to 2019, with and without adjustment for daily ozone and nitrogen dioxide concentration.

epih-44-e2022052-suppl7.docx

Supplementary Material 8.

Mortality effects of short-term ambient PM_2.5 exposure and their modification by consecutive days of exposure to high concentrations of PM_2.5 in seven major cities in Korea from 2006 to 2019 among those aged 20 and 64 years.

epih-44-e2022052-suppl8.docx

Supplementary Material 9.

epih-44-e2022052-suppl9.docx

Supplementary Material 10.

Comparison of the effect modifications of pooled effect estimates (all-cause mortality risk percent changes per 10 μg/m³ increase of daily lag 0-1) of the average of same-day exposure to PM_2.5 and PM_2.5 exposure 1 day prior. Left is according to consecutive days, and right is to background mean concentration interval (μg/m³).

epih-44-e2022052-suppl10.docx

NOTES

CONFLICT OF INTEREST

The authors have no conflicts of interest to declare for this study.

FUNDING

This research was funded by the Ministry of Environment, Republic of Korea, for the project “Mid/long-term health effect assessment research on particulate matter exposure.”

AUTHOR CONTRIBUTIONS

Conceptualization: Lim H, Bae S, Bae HJ, Kwon HJ. Data curation: Lim H. Formal analysis: Lim H. Funding acquisition: Kwon HJ. Methodology: Lim H, Bae S, Choi J, Choi KH, Bae HJ. Project administration: Kwon HJ. Visualization: Lim H, Choi J. Writing – original draft: Lim H. Writing – review & editing: Choi J, Bae S, Choi KH, Bae HJ, Kim S, Ha M, Kwon HJ.

ACKNOWLEDGEMENTS

None.

Figure 1.

Distribution of consecutive days with high concentrations (35 μg/m³ or more) of fine particulate matter (PM_2.5) in seven cities (A: Seoul, B: Busan, C: Daegu, D: Incheon; E: Gwangju, F: Daejeon, and G: Ulsan) in Korea by year from 2006 to 2019. The x-axis represents the year, and the y-axis represents the number of days with a high PM_2.5 concentration (reference: 35 μg/m³) by year. A darker shade of blue indicates a greater number of days of exposure to high concentrations of PM_2.5.

Figure 2.

Non-traumatic all-cause (A), respiratory (B), and cardiovascular (C) mortality risk percent changes per 10 μg/m³ increase in the daily average fine particulate matter concentration and their modification by consecutive days of high concentration (35 μg/m³ or more) in a model with applied lag 0-1 exposure. This graph shows the percent of risk among the total population (left column), those aged 20-65 years (middle column), and those aged 65 years or older (right column). The red square boxes are estimates of the effects over the whole study period in basic models without applying consecutive day variables. Others represent estimates of six coefficients in the strata of consecutive days (no high-concentration days, the first to fourth days, and the fifth day and beyond) in the effect modification models. Vertical lines with point estimates indicate 95% confidence intervals. ¹According to International Classification of Diseases, 10th edition code.

Table 1.

All-cause (non-traumatic), respiratory, and cardiovascular mortality in 7 major cities in Korea from 2006 to 2019

City	Population in 2019	All-cause mortality (non-traumatic)		Respiratory mortality		Cardiovascular mortality
City	Population in 2019	Total	Daily	Total	Daily	Total	Daily
Seoul	9,729,107	517,265	101.2±13.7	43,043	8.4±4.1	117,750	23.0±5.5
Busan	3,413,841	256,534	50.2±8.8	23,763	4.6±2.7	72,014	14.1±4.2
Daegu	2,438,031	157,211	30.7±6.6	15,373	3.0±2.1	42,042	8.2±3.0
Incheon	2,957,026	160,818	31.5±7.0	15,425	3.0±2.1	40,721	8.0±2.9
Gwangju	1,456,468	84,898	16.6±4.6	9,342	1.8±1.6	19,419	3.8±2.0
Daejeon	1,474,870	81,039	15.8±4.4	7,798	1.5±1.3	18,759	3.7±2.0
Ulsan	1,148,019	57,064	11.2±3.6	5,350	1.0±1.1	15,329	3.0±1.8
Total	22,617,362	1,314,829	257.2±31.5	120,094	23.5±9.8	326,034	63.8±10.4

Values are presented as number or mean±standard deviation.

Table 2.

Daily fine particulate matter (PM_2.5) concentrations and the proportion of exposure to high concentrations in seven major cities in Korea for 5,113 days from 2006 to 2019

City	Daily PM_2.5 (μg/m³)	High-concentration days¹	High-concentration duration¹
City	Daily PM_2.5 (μg/m³)	High-concentration days¹	1st day	2nd day	3rd day	4th day	5th day or more
Seoul	27.4±15.8	1,171 (22.9)	492 (9.6)	281 (5.5)	169 (3.3)	92 (1.8)	137 (2.7)
Busan	26.6±12.9	1,039 (20.3)	484 (9.5)	266 (5.2)	140 (2.7)	70 (1.4)	79 (1.5)
Daegu	25.3±13.2	919 (18.0)	426 (8.3)	238 (4.7)	132 (2.6)	65 (1.3)	57 (1.1)
Incheon	26.0±14.4	1,008 (19.7)	457 (8.9)	257 (5.0)	130 (2.5)	69 (1.3)	93 (1.8)
Gwangju	23.0±12.9	693 (13.6)	321 (6.3)	175 (3.4)	89 (1.7)	50 (1.0)	57 (1.1)
Daejeon	23.7±13.4	778 (15.2)	365 (7.1)	199 (3.9)	103 (2.0)	53 (1.0)	57 (1.1)
Ulsan	23.9±12.8	825 (16.1)	397 (7.8)	222 (4.3)	110 (2.2)	53 (1.0)	41 (0.8)

Values are presented as mean±standard deviation or number (%).

¹ When the average 24-hour mean PM_2.5 concentration was 35 μg/m³ or more, it was defined as a high-concentration day; We then specified the number of days of high concentrations for the duration.

Table 3.

Mortality effects of short-term ambient fine particulate matter (PM_2.5) exposure in seven major cities in Korea from 2006 to 2019

Lag structure²	All-cause (A00-R99)¹			Respiratory (J00-J98)¹			Cardiovascular (I00-I99)¹
	% Change³	95% CI		% Change³	95% CI		% Change³	95% CI
	% Change³	LL	UL	% Change³	LL	UL	% Change³	LL	UL
Single-day
0	0.33	0.18	0.47	0.27	-0.21	0.77	0.39	0.10	0.68
1	0.16	0.02	0.31	0.42	-0.06	0.91	0.00	-0.29	0.29
2	-0.05	-0.20	0.09	0.37	-0.09	0.84	-0.33	-0.60	-0.05
3	-0.05	-0.19	0.09	-0.02	-0.64	0.61	-0.31	-0.58	-0.04
4	-0.20	-0.41	0.00	-0.02	-0.47	0.43	-0.30	-0.87	0.28
5	-0.33	-0.55	-0.11	-0.53	-0.97	-0.08	-0.25	-0.62	0.12
6	-0.29	-0.42	-0.15	-0.01	-0.45	0.44	-0.22	-0.49	0.05
Moving average
0-1	0.33	0.16	0.50	0.47	-0.09	1.04	0.26	-0.08	0.60
0-2	0.25	0.05	0.44	0.62	-0.01	1.26	0.01	-0.37	0.40
0-3	0.19	-0.02	0.40	0.62	-0.07	1.32	-0.17	-0.59	0.25
0-4	0.06	-0.17	0.29	0.52	-0.26	1.31	-0.34	-0.96	0.28
0-5	-0.14	-0.43	0.16	0.24	-0.55	1.04	-0.47	-1.19	0.26
0-6	-0.28	-0.61	0.04	0.22	-0.62	1.07	-0.55	-1.24	0.15

City-specific estimated PM_2.5 effects were from quasi-Poisson generalized additive models and pooled with the same exposure lag structure using random-effects meta-analyses.

CI, confidence interval; LL, lower limit; UL, upper limit.

¹ According to International Classification of Diseases, 10th edition code.

² Single-day: the average 24-hour mean PM_2.5 concentration from 0 (same day) to 6 days before; moving average: the average PM_2.5 concentration over the current and previous day and the current to before the sixth day.

³ The % increase in mortality risk per 10 μg/m³ increase in PM_2.5 concentration.

Table 4.

Mortality effects of short-term ambient fine particulate matter (PM_2.5) exposure and their modification by consecutive days of high PM_2.5 concentration in seven major cities in Korea from 2006 to 2019

Variables	All-cause (A00-R99)¹			Respiratory (J00-J98)¹			Cardiovascular (I00-I99)¹
	% Change²	95% CI		% Change²	95% CI		% Change²	95% CI
	% Change²	LL	UL	% Change²	LL	UL	% Change²	LL	UL
Basic model	0.33	0.16	0.50	0.47	-0.09	1.04	0.26	-0.08	0.60
Effect modification model
Consecutive days³
No	0.35	-0.13	0.84	0.17	-0.84	1.18	0.61	0.00	1.23
1st day	0.63	0.20	1.06	0.82	-0.01	1.65	0.59	0.09	1.10
2nd day	0.24	0.02	0.46	0.18	-0.56	0.93	0.36	-0.09	0.80
3rd day	0.16	-0.10	0.41	0.23	-0.60	1.06	0.09	-0.40	0.58
4th day	0.36	0.01	0.70	1.33	0.36	2.31	0.47	-0.12	1.06
5th day or more	0.28	-0.10	0.65	-0.002	-0.85	0.85	0.35	-0.32	1.02

City-specific estimated PM_2.5 effects were from quasi-Poisson generalized additive models with the 2-day moving average (the average over the current and previous day, lag 0-1) of PM_2.5 concentration and pooled with the same consecutive-day strata using random-effects meta-analyses.

CI, confidence interval; UL, upper limit; LL, lower limit.

¹ According to International Classification of Diseases, 10th edition code.

² The % increase in mortality risk per 10 μg/m³ increase in PM_2.5 concentration.

³ The 6-strata categorical variable designating the number of consecutive days with daily mean PM_2.5 concentrations of 35 μg/m³ or more.

REFERENCES

1. Health Effects Institute. State of global air/2020: a special report on global exposure to air pollution and its health impacts; 2020 [cited 2022 May 20]. Available from: https://www.stateofglobalair.org/sites/default/files/documents/2020-10/soga-2020-report-10-26_0.pdf.
2. Korean Ministry of Environment. Enhancement of fine dust (PM2.5) environmental standards to that of the US and Japan; 2018 [cited 2022 Feb 7]. Available from: http://www.me.go.kr/home/web/board/read.do?menuId=286&boardMasterId=1&boardCategoryId=39&boardId=849570 (Korean).
3. United States Environmental Protection Agency. Integrated Science Assessment (ISA) for particulate matter; 2019 [cited 2022 Feb 7]. Available from: https://www.epa.gov/isa/integrated-science-assessment-isa-particulate-matter.
4. Liu C, Chen R, Sera F, Vicedo-Cabrera AM, Guo Y, Tong S, et al. Ambient particulate air pollution and daily mortality in 652 cities. N Engl J Med 2019;381:705-715.Article PubMed PMC
5. Lee H, Honda Y, Hashizume M, Guo YL, Wu CF, Kan H, et al. Short-term exposure to fine and coarse particles and mortality: a multicity time-series study in East Asia. Environ Pollut 2015;207:43-51.Article PubMed
6. Baxter LK, Crooks JL, Sacks JD. Influence of exposure differences on city-to-city heterogeneity in PM2.5-mortality associations in US cities. Environ Health 2017;16:1.Article PubMed PMC PDF
7. Zanobetti A, Schwartz J. The effect of fine and coarse particulate air pollution on mortality: a national analysis. Environ Health Perspect 2009;117:898-903.Article PubMed PMC
8. Schwartz J, Bind MA, Koutrakis P. Estimating causal effects of local air pollution on daily deaths: effect of low levels. Environ Health Perspect 2017;125:23-29.Article PubMed
9. Wei Y, Wang Y, Wu X, Di Q, Shi L, Koutrakis P, et al. Causal effects of air pollution on mortality rate in Massachusetts. Am J Epidemiol 2020;189:1316-1323.Article PubMed PMC PDF
10. National Institute of Environmental Research. Annual report of air quality in Korea 2020; 2021 [cited 2022 May 20]. Available from: https://www.airkorea.or.kr/web/detailViewDown?pMENU_NO=125 (Korean).
11. National Council on Climate and Air Quality. Report on the current state and solutions for fine particulate air pollution; 2020 [cited 2022 Feb 7]. Available from: https://www.ctis.re.kr/en/downloadBbsFile.do?atchmnflNo=4962 (Korean).
12. Wang J, Yin Q, Tong S, Ren Z, Hu M, Zhang H. Prolonged continuous exposure to high fine particulate matter associated with cardiovascular and respiratory disease mortality in Beijing, China. Atmos Environ 2017;168:1-7.Article
13. Kim SE, Bell ML, Hashizume M, Honda Y, Kan H, Kim H. Associations between mortality and prolonged exposure to elevated particulate matter concentrations in East Asia. Environ Int 2018;110:88-94.Article PubMed
14. Kim BU, Bae C, Kim HC, Kim E, Kim S. Spatially and chemically resolved source apportionment analysis: case study of high particulate matter event. Atmos Environ 2017;162:55-70.Article
15. Kim HC, Kim E, Bae C, Cho JH, Kim BU, Kim S. Regional contributions to particulate matter concentration in the Seoul metropolitan area, South Korea: seasonal variation and sensitivity to meteorology and emissions inventory. Atmos Chem Phys 2017;17:10315-10332.Article
16. Lee KS, Lim YH, Choi YJ, Kim S, Bae HJ, Han C, et al. Prenatal exposure to traffic-related air pollution and risk of congenital diseases in South Korea. Environ Res 2020;191:110060.Article PubMed
17. Lim YH, Kim S, Han C, Bae HJ, Seo SC, Hong YC. Source country-specific burden on health due to high concentrations of PM2.5. Environ Res 2020;182:109085.Article PubMed
18. Roberts S. Interactions between particulate air pollution and temperature in air pollution mortality time series studies. Environ Res 2004;96:328-337.Article PubMed
19. Kim SE, Lim YH, Kim H. Temperature modifies the association between particulate air pollution and mortality: a multi-city study in South Korea. Sci Total Environ 2015;524-525:376-383.Article PubMed
20. Bae S, Kwon HJ. Current state of research on the risk of morbidity and mortality associated with air pollution in Korea. Yonsei Med J 2019;60:243-256.Article PubMed PMC PDF
21. Kim SE, Honda Y, Hashizume M, Kan H, Lim YH, Lee H, et al. Seasonal analysis of the short-term effects of air pollution on daily mortality in Northeast Asia. Sci Total Environ 2017;576:850-857.Article PubMed
22. Jung EJ, Na W, Lee KE, Jang JY. Elderly mortality and exposure to fine particulate matter and ozone. J Korean Med Sci 2019;34:e311.Article PubMed PMC PDF
23. Kim TY, Kim H, Yi SM, Cheong JP, Heo J. Short-term effects of ambient PM2.5 and PM2.5-10 on mortality in major cities of Korea. Aerosol Air Qual Res 2018;18:1853-1862.Article
24. Kim MG, Lee SJ, Park D, Kim CH, Lee KH, Hwang JM. Relationship between the actual fine dust concentration and media exposure that influenced the changes in outdoor activity behavior in South Korea. Sci Rep 2020;10:12006.Article PubMed PMC PDF
25. Korea Society Opinion Institute. Public opinion poll on wearing a fine dust mask; 2019 [cited 2022 Feb 7]. Available from: https://opengov.seoul.go.kr/council/18004939 (Korean).
26. Zeger SL, Dominici F, Samet J. Harvesting-resistant estimates of air pollution effects on mortality. Epidemiology 1999;10:171-175.Article PubMed
27. Schwartz J. Is there harvesting in the association of airborne particles with daily deaths and hospital admissions? Epidemiology 2001;12:55-61.Article PubMed
28. National Institute of Environmental Research. Case study of high PM episodes observed in intensive monitoring station (2017); 2018 [cited 2022 Feb 7]. Available from: http://dl.nanet.go.kr/law/SearchDetailView.do?cn= MONO1201834601#none (Korean).
29. Heo J, Schauer JJ, Yi O, Paek D, Kim H, Yi SM. Fine particle air pollution and mortality: importance of specific sources and chemical species. Epidemiology 2014;25:379-388.Article PubMed
30. Son JY, Lee JT, Kim KH, Jung K, Bell ML. Characterization of fine particulate matter and associations between particulate chemical constituents and mortality in Seoul, Korea. Environ Health Perspect 2012;120:872-878.Article PubMed PMC
31. Hwang SH, Lee JY, Yi SM, Kim H. Associations of particulate matter and its components with emergency room visits for cardiovascular and respiratory diseases. PLoS One 2017;12:e0183224.Article PubMed PMC
32. Yoo SE, Park JS, Lee SH, Park CH, Lee CW, Lee SB, et al. Comparison of short-term associations between PM2.5 components and mortality across six major cities in South Korea. Int J Environ Res Public Health 2019;16:2872.Article PubMed PMC
33. Yang Y, Ruan Z, Wang X, Yang Y, Mason TG, Lin H, et al. Short-term and long-term exposures to fine particulate matter constituents and health: a systematic review and meta-analysis. Environ Pollut 2019;247:874-882.Article PubMed
34. Yan M, Wilson A, Bell ML, Peng RD, Sun Q, Pu W, et al. The shape of the concentration-response association between fine particulate matter pollution and human mortality in Beijing, China, and its implications for health impact assessment. Environ Health Perspect 2019;127:67007.Article PubMed