Journal List > J Korean Med Sci > v.38(2) > 1516080856

Yang, Jeong, Kweon, and Shin: Association Between Regional Levels of Particulate Matter and Recurrent Falls in Korea

Abstract

Background

We investigated the extent of regional disparity of recurrent falls. In addition, we examined the association between particulate matter (PM) and recurrent falls and the association between regional disparity of recurrent falls and regional PM levels.

Method

We used data from Korea Community Health Survey 2019 that included 204,395 participants from 237 municipal districts. The independent variables were the annual average PM10 and PM2.5 concentrations measured at the air quality measuring stations in each municipal district. The outcome variable was the experience of falls more than twice in the previous year. Multilevel analyses were conducted to estimate the association between regional PM10 and PM2.5 levels and recurrent falls.

Results

The regional variation was greater in the young people than that in the older people. PM10 and PM2.5 levels were positively associated with recurrent falls after adjusting for individual and regional covariates. These associations were more evident in the older group than in the young. PM10 and PM2.5 explained 2.82% and 3.33% of the remaining regional variance in models with individual and regional confounders, respectively. These proportions were greater in the older group (PM10 and PM2.5; 4.73% and 5.27%) than those in the younger age group (PM10 and PM2.5, 0.80% and 1.39%).

Conclusion

PM concentration was associated with recurrent falls even after accounting for other regional variables and individual-level differences. Moreover, there were regional differences in the occurrence of falls, and the PM concentration explained a part of the gap, but the gap was explained more in the older group than in the young.

Graphical Abstract

jkms-38-e19-abf001.jpg

INTRODUCTION

Falls are one of the leading causes of hospitalization and premature death, accounting for most of the disability due to unintentional injuries in the older population.1 The substantial burden of falls appears to be a global phenomenon as per the Global Burden of Diseases, Injuries and Risk Factors Study of 2017, and falls are the second leading cause of mortality owing to unintentional injuries.2 Additionally, in Korea, disability-adjusted life years (DALY) and years lived with disability (YLD) due to falls are higher as compared to those for other injuries, and DALY and YLD due to falls show increasing trends.3 In the older population, demographic factors, gait disturbance, impaired vision, dizziness, impaired cognition, and cardiovascular disease have been suggested as individual risk factors for falls.4
A variety of adverse health outcomes can be attributed to ambient air pollutants. Particulate matter (PM), particles with diameters of 10 microns or less, is a common proxy indicator of air pollutants and is known to have the greatest impact on human health among the causes of air pollution.5 Previous studies have reported associations between exposure to air pollutants and respiratory disease, cardiovascular morbidity and mortality, and cognitive decline.67 The accumulation of small particles that are not filtered out by barriers present in the respiratory tract may cause a pulmonary or systemic inflammatory response.8 The exposure to PM can contribute to falls at the proximal level because it is associated with the individual-level risk factors of falls.
The fundamental goal of healthcare activity is equality across different socioeconomic groups and geographical areas. Therefore, investigating the inequality of health outcomes and related factors is important. However, a traditional approach to estimate the difference between group health outcomes, mostly represented as averages, may connote the risk of erroneous estimates of regional differences because the method does not consider individual-level heterogeneity.9 Especially in the case of PM, the exposure can be measured and the intervention can be applied at the regional level. Therefore, the size of the regional difference and the effect of the risk factors should be considered simultaneously in terms of primary prevention.
Previous studies have suggested that the PM levels may be associated with the risk of falls or related health outcomes. In a multilevel analysis conducted in six middle- or low-income countries, the three-year PM2.5 concentration was positively associated with fall-related injuries. 1 However, this study did not include comorbidities at the individual level and confounders at the regional level; therefore, the air pollutant levels may act as a proxy for other environmental factors. In Korea, in a study on the association between air quality and hip fracture, PM10 was positively associated, while others were not, with hip fracture.10 A study in seven cities of Korea found that the levels of air pollutants such as SO2, CO, and NO2 were positively associated with death by unintentional injury, most of which were caused by falls. Moreover, PM levels did not show a significant association with unintentional injury deaths.11 These results reflected neither a direct association nor regional differences.
In the present study, we analyzed the association between air pollutant levels and falls and investigated the variation of falls between regions and how much of the variation was explained by air pollutant levels.

METHODS

Study participants

The participants of Korean Community Health Survey (KCHS), conducted by the Korea Centers for Disease Control and Prevention in 2019, were included in this study. The KCHS was a nationwide survey to investigate health statistics for adults aged 19 and over by regional district.12 In 229,099 participants of 2019 KCHS, 11,519 subjects of 13 districts that did not have an air quality monitoring station or did not measure PM10 or PM2.5 were excluded. Furthermore, 13,185 subjects who had missing values in covariates were excluded. Finally, 204,395 participants (139,971 participants aged under 65 and 64,424 participants aged 65 or over) from 237 municipal districts were included in the analysis.

Variables

Data on age, sex, marital status, education status, household income, drinking, smoking, physical activity, categorized BMI, depressive symptoms, cognitive impairment, diabetes, hypertension, urban forest area ratio, urban forest of life zone per capita, population density, PM10, and PM2.5 were included in the analysis. All the covariates were obtained from the KCHS dataset. Marital status was classified into married, unmarried, and others. Education status was defined as public education lasting > 9 years. Household income was measured with questionnaires, and the household income quartiles were defined by dividing the household incomes of all the participants into four quartiles. Drinking status was categorized as never drinker, former drinker, occasional drinker, and current drinker. Subjects who had never consumed alcohol in their lifetime were categorized as never drinker, and those who had a history of drinking, but did not currently drink were categorized as former drinkers. Among the participants who were currently consuming alcohol, those who drank at a frequency of less than once a month were categorized as occasional drinkers, and those who drank at a higher frequency were categorized as current drinkers. Smoking status was classified into never smoker, former smoker, and current smoker according to whether the participants were currently smoking and whether they had smoked in their lifetime. Physical activity was defined as walking for at least 5 days a week and for 30 minutes or more at a time. BMI was calculated with the heights and weights, and subjects were divided into four groups: underweight, normal weight, overweight, and obesity.13 Depressive symptoms were measured using the Korean version of the Patient Health Questionnaire-9 (PHQ-9) which is a simple screening tool for depression.14 Presence of depressive symptoms was defined as a score of 10 or higher on the PHQ-9. Cognitive impairment was defined as experiencing an exacerbated memory loss or confusion. Diabetes and hypertension were determined by treatment or a previous diagnosis. Urban forest area ratio and urban forest of life zone per capita in 2019, provided by Korea Forest Service, were used as regional-level covariates, and the participants were divided into quartiles according to them. Population density data of each district was obtained from the Ministry of the Interior and Safety of Korea. For PM10 and PM2.5, data of the annual exposure of each area from the annual air quality report of 2019 provided by the ministry of environment was used. If two or more measuring stations were located in the area, the average of the values of the stations was used. More than two falls in the past 12 months was defined as the outcome using the following two questions about fall experiences as follows: “Have you fallen in the past year?” and “How many times have you fallen in the past year?”

Statistical analyses

General characteristics were presented according to the presence or absence of recurrent fall experiences, and additionally stratified according to age group. Characteristics of participants included in the analysis were shown as mean ± standard deviation or number (%). Student’s t-test and chi-square tests were used to show the significance of difference according to recurrent fall experiences. Multilevel logistic regression with individuals nested within 237 municipal districts was conducted to estimate the association between air pollutants and recurrent falls using a four-step model. The null model included only the intercept. Model 1 included individual-level variables only. Model 2 included Model 1 variables plus urban forest area ratio and urban forest area of life zone per capita and population density. Finally, in Model 3, all the individual-level and regional-level variables were included. To estimate the regional effects in the total variance of recurrent fall experiences, intra-class correlation coefficient (ICC) was calculated using the latent variable approach.15 To assess the attribution of adjusted variables to the regional variance of recurrent fall experiences, proportional change in variance (PCV) was calculated in all the models except the null model. The PCV of the model was calculated as the proportional change in the area level variance between the null model and Model 1; Model 1 and Model 2; and Model 2 and Model 3. Because the characteristics and effect of falls on overall health may vary by age, age-stratified analyses with a cutoff of 65 years were conducted. All the analyses were performed using R (ver. 4.0.2, Vienna, Austria).

Ethics statement

This study was covered under the review list pursuant to Article 2.2 of the Enforcement Rule of Bioethics and Safety Act in Korea, and the data were exempted from Institutional Review Board review (IRB no. 1040198-221208-HR-151-01). All data were fully anonymized before the data were provided.

RESULTS

Table 1 presents the characteristics of the subjects by recurrent fall experience in each age group. In total population, subjects that experienced recurrent falls were younger, had higher proportions of women, participants living alone, less educated, and had lower incomes as compared to the group who did not experience falls. Furthermore, depressive symptom, hypertension, diabetes, and abnormal BMI were more prevalent in the group that experienced recurrent falls. In young group, subjects who experienced recurrent falls had a higher proportion of urban residents, current drinkers, current smokers, hypertension, and physically active people, while in older group, subjects who experienced recurrent had lower proportion of urban residents, current drinkers, current smokers, hypertension, and physically active people. In older group, cognitive impairment was more prevalent in fall group than in normal group.
Table 1

Descriptive statistics of study participants by recurrent fall experiences

jkms-38-e19-i001
Variables Total population Young group (aged under 65) Older group (aged 65 or over)
No Yes P value No Yes P value No Yes P value
No. (%) 193,692 (94.8) 10,703 (5.2) 133,128 (95.1) 6,843 (4.9) 60,564 (94.0) 3,860 (6.0)
Age, yr 54.63 (17.19) 52.91 (21.30) < 0.001 45.76 ± 12.63 39.82 ± 14.46 < 0.001 74.14 ± 6.40 76.11 ± 6.70 < 0.001
Sex (female) 105,427 (54.4) 6,316 (59.0) < 0.001 70,749 (53.1) 3,756 (54.9) 0.005 34,678 (57.3) 2,560 (66.3) < 0.001
Marital status < 0.001 < 0.001 < 0.001
Unmarried 29,474 (15.2) 2,848 (26.6) 29,169 (21.9) 2,828 (41.3) 305 (0.5) 20 (0.5)
Married 130,787 (67.5) 5,447 (50.9) 91,381 (68.6) 3,347 (48.9) 39,406 (65.1) 2,100 (54.4)
Other (divorced, etc.) 33,431 (17.3) 2,408 (22.5) 12,578 (9.4) 668 (9.8) 20,853 (34.4) 1,740 (45.1)
Education (< 9 yr) 46,066 (23.8) 3,433 (32.1) < 0.001 9,183 (6.9) 575 (8.4) < 0.001 36,883 (60.9) 2,858 (74.0) < 0.001
Household income < 0.001 < 0.001 < 0.001
First quartile 52,266 (27.0) 3,924 (36.7) 16,496 (12.4) 1,306 (19.1) 35,770 (59.1) 2,618 (67.8)
Second quartile 43,809 (22.6) 2,202 (20.6) 30,795 (23.1) 1,532 (22.4) 13,014 (21.5) 670 (17.4)
Third quartile 51,700 (26.7) 2,460 (23.0) 43,925 (33.0) 2,066 (30.2) 7,775 (12.8) 394 (10.2)
Fourth quartile 45,917 (23.7) 2,117 (19.8) 41,912 (31.5) 1,939 (28.3) 4,005 (6.6) 178 (4.6)
Depressive symptom (yes) 5,118 (2.6) 1,174 (11.0) < 0.001 2,850 (2.1) 695 (10.2) < 0.001 2,268 (3.7) 479 (12.4) < 0.001
Cognitive impairment (yes)a - - - - - - 16,699 (27.6) 1,800 (46.6) < 0.001
Diagnosis of diabetes 21,702 (11.2) 1,505 (14.1) < 0.001 8,845 (6.6) 481 (7.0) 0.213 12,857 (21.2) 1,024 (26.5) < 0.001
Diagnosis of hypertension 54,337 (28.1) 3,208 (30.0) < 0.001 21,196 (15.9) 979 (14.3) < 0.001 33,141 (54.7) 2,229 (57.7) < 0.001
Residence area (urban) 111,488 (57.6) 6,151 (57.5) 0.855 85,495 (64.2) 4,677 (68.3) < 0.001 25,993 (42.9) 1,474 (38.2) < 0.001
Alcohol consumption < 0.001 < 0.001 < 0.001
Never drinker 36,066 (18.6) 1,899 (17.7) 15,986 (12.0) 502 (7.3) 20,080 (33.2) 1,397 (36.2)
Former drinker 31,974 (16.5) 1,936 (18.1) 16,600 (12.5) 792 (11.6) 15,374 (25.4) 1,144 (29.6)
Occasional drinker 28,373 (14.6) 1,632 (15.2) 20,829 (15.6) 1,180 (17.2) 7,544 (12.5) 452 (11.7)
Current drinker 97,279 (50.2) 5,236 (48.9) 79,713 (59.9) 4,369 (63.8) 17,566 (29.0) 867 (22.5)
Smoking status < 0.001 < 0.001 < 0.001
Never smoker 124,983 (64.5) 6,922 (64.7) 85,533 (64.2) 4,207 (61.5) 39,450 (65.1) 2,715 (70.3)
Former smoker 37,096 (19.2) 1,785 (16.7) 21,222 (15.9) 945 (13.8) 15,874 (26.2) 840 (21.8)
Current smoker 31,613 (16.3) 1,996 (18.6) 26,373 (19.8) 1,691 (24.7) 5,240 (8.7) 305 (7.9)
Physical activity (yes) 79,169 (40.9) 4,295 (40.1) 0.127 55,903 (42.0) 3,184 (46.5) < 0.001 23,266 (38.4) 1,111 (28.8) < 0.001
Categories of BMI < 0.001 < 0.001 < 0.001
Underweight 7,998 (4.1) 619 (5.8) 5,513 (4.1) 417 (6.1) 2,485 (4.1) 202 (5.2)
Normal 118,061 (61.0) 6,060 (56.6) 82,155 (61.7) 3,892 (56.9) 35,906 (59.3) 2,168 (56.2)
Overweight 57,741 (29.8) 3,241 (30.3) 38,155 (28.7) 1,973 (28.8) 19,586 (32.3) 1,268 (32.8)
Obesity 9,892 (5.1) 783 (7.3) 7,305 (5.5) 561 (8.2) 2,587 (4.3) 222 (5.8)
All values are presented as the mean ± standard deviation or n (%). Student’s t test and χ2 tests were conducted for continuous and categorical variables, respectively.
BMI = body mass index.
aCognitive impairment was presented in participants aged 65 or over.
Table 2 presents the results of the multilevel analysis of recurrent fall experiences. In the null model, the proportion of variance of recurrent fall experiences explained by region was 5.00%. Individual-level variables decreased the ICC of the null model to 4.67%, and regional confounders decreased the ICC of model 1 to 4.25%. PM10 and PM2.5 levels reduced the ICC of model 2 by 4.14% and 4.11%, respectively. After adjusting for individual and urban forest variables, PM10 (odds ratio [OR], 1.131; 95% confidence interval [CI], 1.026–1.247) and PM2.5 (OR, 1.202; 95% CI, 1.051–1.374) levels increased the risk of recurrent falls significantly. PM10 and PM2.5 accounted for 2.82% and 3.33% of the unexplained regional variance in Model 2, respectively. Among the individual-level covariates included in the analysis, old age, female, unmarried, low education, low income, drinking, smoking, underweight or obesity, depressive symptom, and diabetes increased the risk of recurrent fall experiences.
Table 2

Results of the multilevel logistic regression analysis of recurrent fall experiences

jkms-38-e19-i002
Recurrent fall experience Model 1 Model 2 Model 3 (PM10) Model 3 (PM2.5)
Individual level
Age (65 yr or over) 1.131 (1.062–1.204) 1.130 (1.061–1.203) 1.130 (1.061–1.203) 1.130 (1.061–1.203)
Sex (female) 1.394 (1.315–1.478) 1.393 (1.314–1.477) 1.393 (1.314–1.477) 1.393 (1.314–1.477)
Marital status (ref. unmarried)
Married 0.376 (0.357–0.397) 0.377 (0.357–0.397) 0.377 (0.357–0.397) 0.377 (0.357–0.397)
Other (divorced, etc.) 0.461 (0.429–0.495) 0.462 (0.430–0.496) 0.462 (0.430–0.496) 0.462 (0.430–0.496)
Education (9 yr or over) 0.709 (0.666–0.754) 0.708 (0.665–0.753) 0.707 (0.665–0.752) 0.707 (0.665–0.753)
Household income (ref. first quartile)
Second quartile 0.778 (0.733–0.826) 0.778 (0.733–0.825) 0.777 (0.732–0.825) 0.777 (0.732–0.825)
Third quartile 0.756 (0.711–0.805) 0.756 (0.710–0.804) 0.754 (0.709–0.803) 0.755 (0.709–0.803)
Fourth quartile 0.720 (0.674–0.770) 0.719 (0.673–0.769) 0.718 (0.672–0.767) 0.718 (0.672–0.768)
Alcohol consumption (ref. never drinker)
Former drinker 1.186 (1.108–1.270) 1.186 (1.108–1.269) 1.185 (1.108–1.269) 1.185 (1.108–1.269)
Occasional drinker 1.182 (1.099–1.270) 1.181 (1.099–1.270) 1.181 (1.099–1.269) 1.181 (1.099–1.269)
Current drinker 1.237 (1.162–1.317) 1.237 (1.162–1.316) 1.236 (1.162–1.316) 1.236 (1.162–1.316)
Smoking status (ref. never smoker)
Former smoker 1.199 (1.118–1.285) 1.199 (1.119–1.285) 1.199 (1.118–1.285) 1.199 (1.118–1.285)
Current smoker 1.381 (1.293–1.475) 1.381 (1.293–1.475) 1.381 (1.293–1.475) 1.381 (1.293–1.475)
Physical activity (yes) 1.005 (0.963–1.048) 1.003 (0.962–1.046) 1.003 (0.962–1.046) 1.003 (0.962–1.046)
BMI categories (ref. normal)
Underweight 1.158 (1.060–1.265) 1.157 (1.059–1.264) 1.157 (1.059–1.264) 1.157 (1.059–1.264)
Overweight 1.125 (1.075–1.177) 1.125 (1.075–1.177) 1.125 (1.075–1.177) 1.125 (1.075–1.177)
Obesity 1.361 (1.257–1.475) 1.361 (1.256–1.475) 1.361 (1.256–1.474) 1.361 (1.256–1.474)
Depressive symptom (yes) 3.483 (3.247–3.736) 3.479 (3.243–3.731) 3.476 (3.241–3.728) 3.475 (3.240–3.727)
Diagnosis of diabetes (yes) 1.221 (1.149–1.299) 1.222 (1.149–1.299) 1.221 (1.149–1.299) 1.221 (1.149–1.299)
Diagnosis of hypertension (yes) 0.973 (0.925–1.024) 0.973 (0.925–1.024) 0.973 (0.924–1.024) 0.973 (0.924–1.024)
Residence area (urban) 0.997 (0.931–1.069) 1.027 (0.955–1.105) 1.035 (0.962–1.113) 1.030 (0.957–1.108)
Regional level
PM10 (per 10 µg/m3) 1.131 (1.026–1.247)
PM2.5 (per 10 µg/m3) 1.202 (1.051–1.374)
Urban forest area ratio (ref. first quartile)a
Second quartile 1.278 (1.096–1.491) 1.277 (1.097–1.487) 1.278 (1.098–1.488)
Third quartile 1.094 (0.937–1.278) 1.094 (0.938–1.275) 1.097 (0.942–1.279)
Fourth quartile 0.950 (0.804–1.123) 0.981 (0.830–1.160) 0.985 (0.833–1.163)
Urban forest of life zone per capita (ref. first quartile)a
Second quartile 1.084 (0.924–1.273) 1.064 (0.907–1.247) 1.063 (0.907–1.246)
Third quartile 1.052 (0.885–1.252) 1.059 (0.892–1.257) 1.043 (0.879–1.237)
Fourth quartile 1.022 (0.854–1.223) 1.060 (0.885–1.268) 1.034 (0.866–1.234)
Population density (per person/km2) 1.000 (1.000–1.000) 1.000 (1.000–1.000) 1.000 (1.000–1.000)
ICC 4.67% 4.25% 4.14% 4.11%
PCVb 6.97% 9.40% 2.82% 3.33%
Values are presented as beta coefficients (95% confidence interval). In the null model, the ICC was 5.00%
PM = particulate matter, BMI = body mass index, ICC = intraclass correlation coefficient, PCV = proportional change in variance.
aUrban forest area ratio and urban forest of life zone per capita were acquired for 237 regions and categorized into quartiles; bThe proportional change in the area level variance between the empty model and model 1; model 1 and model 2; and, model 2 and model 3.
Tables 3 and 4 present the results of multilevel analyses in the older and younger people, respectively. In all the models, the ICC in the young group was larger than that in the older group. PM10 (OR, 1.164; 95% CI, 1.050–1.292) and PM2.5 (OR, 1.247; 95% CI, 1.081–1.438) levels were associated with recurrent falls in the older population. However, these associations were attenuated in the young (PM10 [OR, 1.081; 95% CI, 0.967–1.208]; PM2.5 [OR, 1.147; 95% CI, 0.984–1.338]). Furthermore, the PCV of model 3, the proportion of remaining regional variance of model 2 that was explained by PM, were larger in the older people (PM10: 4.73%; PM2.5: 5.27%) than that in the younger people (PM10: 0.80%; PM2.5: 1.39%). The direction of the effect of age and physical activity on recurrent falls was different depending on the age group. The risk of recurrent falls increased as age decreased in the young group and as age increased in the older group. Physical activity increased the risk of recurrent falls in the young group and decreased the risk of recurrent falls in the older group. Sex, education level, income level, smoking, diabetes, obesity and depressive symptoms were associated with an increased risk of recurrent falls in both age groups. Marital status, alcohol consumption and hypertension were associated with an increased risk of recurrent fall in the younger age group, but not in the older group.
Table 3

Results of the multilevel logistic regression analysis of recurrent fall experiences in older people (≥ 65 years)

jkms-38-e19-i003
Recurrent fall experience Model 1 Model 2 Model 3 (PM10) Model 3 (PM2.5)
Individual level
Age, yr 1.025 (1.019–1.031) 1.025 (1.019–1.031) 1.025 (1.019–1.030) 1.025 (1.019–1.030)
Sex (female) 1.204 (1.068–1.358) 1.205 (1.069–1.358) 1.203 (1.067–1.357) 1.204 (1.068–1.358)
Marital status (ref. unmarried)
Married 0.810 (0.508–1.290) 0.813 (0.510–1.294) 0.816 (0.512–1.300) 0.815 (0.512–1.299)
Other (divorced, etc.) 0.901 (0.564–1.437) 0.905 (0.567–1.444) 0.909 (0.569–1.451) 0.909 (0.569–1.450)
Education (9 yr or over) 0.758 (0.696–0.827) 0.758 (0.695–0.827) 0.756 (0.694–0.825) 0.757 (0.694–0.825)
Household income (ref. first quartile)
Second quartile 0.899 (0.820–0.986) 0.900 (0.821–0.986) 0.896 (0.817–0.982) 0.897 (0.818–0.983)
Third quartile 0.941 (0.838–1.056) 0.940 (0.838–1.055) 0.935 (0.833–1.050) 0.936 (0.834–1.051)
Fourth quartile 0.856 (0.728–1.007) 0.856 (0.728–1.007) 0.851 (0.723–1.001) 0.852 (0.724–1.002)
Alcohol consumption (ref. never drinker)
Former drinker 1.123 (1.030–1.225) 1.125 (1.031–1.227) 1.123 (1.029–1.225) 1.123 (1.029–1.225)
Occasional drinker 1.003 (0.896–1.124) 1.005 (0.897–1.126) 1.003 (0.895–1.123) 1.003 (0.895–1.123)
Current drinker 0.966 (0.875–1.067) 0.966 (0.875–1.067) 0.966 (0.874–1.066) 0.965 (0.874–1.066)
Smoking status (ref. never smoker)
Former smoker 1.066 (0.945–1.203) 1.067 (0.946–1.204) 1.066 (0.945–1.203) 1.067 (0.946–1.204)
Current smoker 1.200 (1.034–1.393) 1.200 (1.034–1.393) 1.198 (1.033–1.391) 1.199 (1.034–1.392)
Physical activity (yes) 0.789 (0.731–0.851) 0.789 (0.731–0.851) 0.789 (0.731–0.851) 0.789 (0.732–0.852)
BMI categories (ref. normal)
Underweight 1.045 (0.895–1.221) 1.045 (0.894–1.220) 1.045 (0.894–1.220) 1.045 (0.894–1.220)
Overweight 1.105 (1.026–1.190) 1.104 (1.025–1.189) 1.104 (1.025–1.189) 1.104 (1.025–1.189)
Obesity 1.246 (1.073–1.447) 1.245 (1.072–1.446) 1.242 (1.069–1.442) 1.243 (1.070–1.443)
Depressive symptom (yes) 2.415 (2.160–2.699) 2.410 (2.156–2.694) 2.402 (2.149–2.685) 2.400 (2.147–2.682)
Cognitive impairment (yes) 1.958 (1.825–2.101) 1.953 (1.820–2.096) 1.954 (1.821–2.097) 1.953 (1.820–2.096)
Diagnosis of diabetes (yes) 1.298 (1.201–1.402) 1.298 (1.201–1.402) 1.297 (1.200–1.401) 1.297 (1.200–1.401)
Diagnosis of hypertension (yes) 0.968 (0.903–1.038) 0.968 (0.903–1.038) 0.967 (0.902–1.037) 0.967 (0.902–1.037)
Residence area (urban) 1.069 (0.968–1.181) 1.113 (0.993–1.248) 1.131 (1.009–1.268) 1.119 (0.998–1.253)
Regional level
PM10 (per 10 µg/m3) 1.164 (1.050–1.292)
PM2.5 (per 10 µg/m3) 1.247 (1.081–1.438)
Urban forest area ratio (ref. first quartile)a
Second quartile 1.214 (1.024–1.439) 1.214 (1.026–1.435) 1.217 (1.029–1.439)
Third quartile 1.046 (0.882–1.239) 1.045 (0.884–1.235) 1.047 (0.886–1.238)
Fourth quartile 0.911 (0.760–1.092) 0.952 (0.794–1.141) 0.951 (0.794–1.138)
Urban forest of life zone per capita (ref. first quartile)a
Second quartile 0.946 (0.790–1.134) 0.923 (0.772–1.104) 0.924 (0.773–1.105)
Third quartile 1.031 (0.853–1.246) 1.033 (0.857–1.245) 1.013 (0.840–1.220)
Fourth quartile 0.999 (0.823–1.212) 1.036 (0.855–1.255) 1.007 (0.832–1.218)
Population density (per person/km2) 1.000 (1.000–1.000) 1.000 (1.000–1.000) 1.000 (1.000–1.000)
ICC 3.94% 3.65% 3.49% 3.47%
PCVb 4.98% 7.53% 4.73% 5.27%
Values are presented as beta coefficients (95% confidence interval).
In the null model, the ICC was 4.14%.
PM = particulate matter, BMI = body mass index, ICC = intraclass correlation coefficient, PCV = proportional change in variance.
aUrban forest area ratio and urban forest of life zone per capita were acquired for 237 regions and categorized into quartiles; bThe proportional change in the area level variance between the empty model and model 1; model 1 and model 2; and, model 2 and model 3.
Table 4

Results of the multilevel logistic regression analysis of recurrent fall experiences in the young (< 65 years)

jkms-38-e19-i004
Recurrent fall experience Model 1 Model 2 Model 3 (PM10) Model 3 (PM2.5)
Individual level
Age, yr 0.967 (0.964–0.970) 0.967 (0.964–0.970) 0.967 (0.964–0.970) 0.967 (0.964–0.970)
Sex (female) 1.447 (1.353–1.547) 1.446 (1.352–1.546) 1.446 (1.352–1.546) 1.446 (1.352–1.546)
Marital status (ref. unmarried)
Married 0.665 (0.618–0.716) 0.666 (0.619–0.717) 0.666 (0.619–0.717) 0.666 (0.619–0.717)
Other (divorced, etc.) 0.833 (0.746–0.931) 0.835 (0.747–0.932) 0.834 (0.747–0.932) 0.834 (0.747–0.932)
Education (9 yr or over) 0.550 (0.496–0.610) 0.549 (0.495–0.609) 0.549 (0.495–0.608) 0.549 (0.495–0.609)
Household income (ref. first quartile)
Second quartile 0.661 (0.609–0.716) 0.660 (0.609–0.716) 0.660 (0.609–0.716) 0.660 (0.609–0.716)
Third quartile 0.625 (0.578–0.677) 0.625 (0.577–0.676) 0.624 (0.577–0.675) 0.624 (0.577–0.676)
Fourth quartile 0.611 (0.563–0.662) 0.609 (0.562–0.661) 0.609 (0.561–0.660) 0.609 (0.561–0.660)
Alcohol consumption (ref. never drinker)
Former drinker 1.367 (1.216–1.537) 1.365 (1.214–1.535) 1.365 (1.214–1.535) 1.365 (1.214–1.535)
Occasional drinker 1.450 (1.299–1.620) 1.449 (1.297–1.618) 1.448 (1.297–1.617) 1.448 (1.297–1.617)
Current drinker 1.460 (1.321–1.613) 1.458 (1.320–1.611) 1.458 (1.320–1.610) 1.458 (1.320–1.610)
Smoking status (ref. never smoker)
Former smoker 1.368 (1.254–1.493) 1.368 (1.254–1.493) 1.368 (1.254–1.492) 1.368 (1.254–1.492)
Current smoker 1.525 (1.416–1.643) 1.525 (1.416–1.643) 1.525 (1.416–1.642) 1.525 (1.416–1.642)
Physical activity (yes) 1.158 (1.100–1.219) 1.155 (1.097–1.217) 1.155 (1.097–1.216) 1.155 (1.097–1.216)
BMI categories (ref. normal)
Underweight 1.066 (0.956–1.189) 1.065 (0.955–1.188) 1.065 (0.955–1.188) 1.065 (0.954–1.188)
Overweight 1.176 (1.109–1.247) 1.176 (1.109–1.247) 1.176 (1.109–1.247) 1.176 (1.109–1.247)
Obesity 1.323 (1.201–1.457) 1.323 (1.201–1.457) 1.323 (1.201–1.457) 1.323 (1.201–1.457)
Depressive symptom (yes) 3.657 (3.334–4.011) 3.652 (3.329–4.005) 3.650 (3.328–4.003) 3.649 (3.327–4.002)
Diagnosed diabetes (yes) 1.286 (1.158–1.429) 1.287 (1.159–1.429) 1.287 (1.159–1.429) 1.287 (1.159–1.429)
Diagnosed hypertension (yes) 1.185 (1.094–1.285) 1.185 (1.094–1.285) 1.185 (1.093–1.285) 1.185 (1.093–1.284)
Residence area (urban) 0.973 (0.894–1.059) 1.005 (0.918–1.100) 1.011 (0.923–1.107) 1.008 (0.920–1.103)
Regional level
PM10 (per 10 µg/m3) 1.081 (0.967–1.208)
PM2.5 (per 10 µg/m3) 1.147 (0.984–1.338)
Urban forest area ratio (ref. first quartile)a
Second quartile 1.265 (1.065–1.502) 1.264 (1.066–1.500) 1.265 (1.067–1.501)
Third quartile 1.103 (0.926–1.313) 1.103 (0.927–1.313) 1.107 (0.930–1.316)
Fourth quartile 0.991 (0.820–1.197) 1.011 (0.835–1.223) 1.018 (0.841–1.232)
Urban forest of life zone per capita (ref. first quartile)a
Second quartile 1.126 (0.942–1.346) 1.114 (0.932–1.332) 1.111 (0.930–1.328)
Third quartile 1.062 (0.874–1.290) 1.067 (0.878–1.296) 1.056 (0.870–1.282)
Fourth quartile 1.077 (0.879–1.320) 1.104 (0.899–1.356) 1.089 (0.889–1.333)
Population density (per person/km2) 1.000 (1.000–1.000) 1.000 (1.000–1.000) 1.000 (1.000–1.000)
ICC 5.34% 4.97% 4.94% 4.91%
PCVb 15.53% 7.27% 0.80% 1.39%
Values are presented as beta coefficients (95% confidence interval). In the null model, the ICC was 6.26%.
PM = particulate matter, BMI = body mass index, ICC = intraclass correlation coefficient, PCV = proportional change in variance.
aUrban forest area ratio and urban forest of life zone per capita were acquired for 237 regions and categorized into quartiles; bThe proportional change in the area level variance between the empty model and model 1; model 1 and model 2; and model 2 and model 3.

DISCUSSION

In the present study, PM10 or PM2.5 levels were positively associated with fall experiences in the total population after adjusting for regional-level and individual-level covariates. However, the association between PM and falls was attenuated in the young. The regional difference in recurrent falls was larger in the younger than that in the older people; however, the PM levels accounted for a larger proportion of regional variance in the older people than that in the young.
A previous multilevel analysis of participants from six countries in the SAGE study showed a positive association between long-term exposure of PM2.5 and fall-related injuries; this association was consistent with the results of our study.1 However, the results of the age-stratified analysis of this study differed from those of ours. In the SAGE study, the strength of the association between PM2.5 exposure and fall-related injury did not differ by age. This discrepancy may be due to differences in the adjusted variables in the analysis. Furthermore, there are no studies investigating the association between PM levels and falls in Korea. However, there is one study on the association between PM levels and hip fracture and one study on the association between PM levels and unintentional injury death. These studies may suggest a relationship between PM and falls because most hip fractures and unintentional injury deaths are caused by falls. A retrospective cohort study from 2011 to 2015 on air pollutants and hip fracture was conducted among participants aged ≥ 30 years using the National Health Insurance service cohort database in Korea. This study revealed positive associations between PM and hip fracture, similar to our study results.10 In contrast, a previous study on unintentional injury deaths from 2002 to 2008 in 7 metropolitan cities in Korea did not find a significant association between PM and unintentional injury deaths. This study was conducted as a time-stratified case-crossover analyses over the previous 3 days and investigated the short-term effects of PM. The differences in the study design and outcome variables may have caused the difference between the results of this study and the present one.11
In this study, the regional difference in recurrent falls was larger in the young than that in the older population. This result may be attributed to the characteristics of the falls. According to a previous study, as age increased, the proportion of indoor falls increased and the proportion of falling during vigorous physical activity decreased.16 Indoor falls are associated with worse physical function and poor health, whereas outdoor falls are associated with increased and better physical activity and above-average health.17 The subgroup analysis of young participants in the present study revealed that age and physical activity were inversely associated with falls. Falls in young participants may include a larger proportion of outdoor falls than that in the older participants; therefore, the regional disparities in recurrent falls in the young may be greater than that in the older population. On the other hand, PM levels explained the larger proportion of regional difference in the older population than that in the young, and the strength of the association between PM and falls was greater in the older population than that in the young. Previous studies have reported that the risk of PM exposure for adverse health outcomes increased with age.1819
The biological pathway explaining the induction of falls due to PM exposure has not been elucidated. However, there are some potential mechanisms. First, PM exposure can lead to impaired visual acuity and ophthalmic conditions such as eye irritation or conjunctivitis through the decreased protective function of tear film, cornea, and conjunctiva and the deterioration in the microcirculation that nourishes the retina.2021 Second, fall-related medical conditions may be attributed to systemic oxidative stress and inflammation due to pulmonary accumulation and translocation of PM to the systemic circulation. Particularly, systemic inflammation inducing impaired hemostatic and vascular function results in cardiovascular and neurologic diseases including atherosclerosis, stroke, and neurodegenerative disease. In a human study, mild cognitive impairment was associated with living along crowded roads, suggesting that chronic environmental air pollutant exposure may affect cognitive function.22 In a previous study using the SAGE dataset that showed the positive association between PM2.5 levels and fall-related injury in low- and middle-income countries, a stratified analyses was conducted by dietary antioxidant consumption, including the consumption of fruits and vegetables. In the high-consumption group, the association between PM2.5 levels and fall-related injury was attenuated as compared to that in the low-consumption group. The result suggested an oxidative stress-related mechanism between PM exposure and falls despite the limitation in terms of green leafy vegetables, and fruits intake in low- and middle-income countries as a result of socioeconomic status (SES).1
In our study, low SES, smoking, abnormal BMI, cognitive decline, depressive symptom, and diabetes were associated with increased risk of recurrent falls. These results are consistent with previous studies. SES is known as a risk factor for falls, and the association was explained by an indirect pathway through lifestyle, psychosocial effect, and resilience to sudden health changes.23 Smoking was associated with unintentional injury, because smoking decreases physical function.24 A report suggested U-shaped association between BMI and falls, because underweight may be a marker of frailty and obesity can cause poor gait and postural balance.25 Cognitive function have bidirectional relationship with falls. Falls may be the consequences of cognitive decline and also may deteriorate cognitive function. In the same way, depression may be caused by fall experience and may cause fall through exacerbating cognitive and physical function.26 Diabetes affects sensory, cognitive function and balance secondary to diabetic neuropathy.27 In our study, hypertension and alcohol consumption were associated with recurrent falls in younger age groups, but not in older age groups. Hypertension may be linked with falls as it is associated with several cardiovascular diseases. Although in older group, reduced cerebral blood flow from hypotension is also the risk factor of falls, therefore direction of association between hypertension and falls might be attenuated.28 Previous studies of association between alcohol consumption and falls reported conflicting results and usually did not find significant association in older group. This inconsistency may be accounted for protective effect of alcohol consumption to cognitive function.2930
There are some limitations worth noting in our study. First, air pollutant levels might be estimated crudely. Ambient air pollutant levels could be acquired only in the areas where the measuring stations were located in 2019, and one value was used as the representative value of the administrative region. Second, the effect of air pollutant levels on fall experiences may surrogate the effects of other unknown regional variables. Therefore, in our study, we adjusted for the urban forest variables and other regional-level variables that may be associated with the general condition and lifestyle including physical activity and effect of air pollutant on subjects. Urban forestry may mitigate PM levels in surrounding areas. At the plant level, the foliage surface can capture and accumulate PM, thus acting as a biological filter. At the canopy level, clustered plants have been reported to block PM.31 Third, our study did not measure an individual exposure level and used the annual mean levels of PM as the proxy of actual exposure levels of PM. However, previous study reported high correlations between ambient PM levels and both indoor and personal PM exposure.32 In addition, World Health Organization designates annual mean concentrations of PM as one of the indicators of the global health observatory.33 Finally, it is difficult to investigate temporal effects such as the lag time effect of PM due to the limitations of the cross-sectional study design.
In conclusion, PM concentration was associated with repeated falls even after considering other regional variables and individual-level differences. In addition, there was a regional difference in the occurrence of falls, and the PM concentration explained a part of the gap; however, the gap was explained more in the older population than in the young.

Notes

Funding: This research was supported by a grant (ISSN:2733-5488) from the Korea Disease Control and Prevention Agency.

Disclosure: The authors have no potential conflicts of interest to disclose.

Author Contributions:

  • Conceptualization: Shin MH.

  • Formal analysis: Yang JH, Shin MH.

  • Investigation: Kweon SS, Jeong JA, Shin MH.

  • Methodology: Shin MH.

  • Supervision: Shin MH.

  • Writing - original draft: Yang JH, Shin MH.

  • Writing - review & editing: Yang JH, Jeong JA, Kweon SS, Shin MH.

References

1. Guo Y, Lin H, Shi Y, Zheng Y, Li X, Xiao J, et al. Long-term exposure to ambient PM2.5 associated with fall-related injury in six low- and middle-income countries. Environ Pollut. 2018; 237:961–967. PMID: 29128246.
2. James SL, Lucchesi LR, Bisignano C, Castle CD, Dingels ZV, Fox JT, et al. The global burden of falls: global, regional and national estimates of morbidity and mortality from the Global Burden of Disease Study 2017. Inj Prev. 2020; 26(Supp 1):i3–11. PMID: 31941758.
3. Kim YE, Park H, Jo MW, Oh IH, Go DS, Jung J, et al. Trends and patterns of burden of disease and injuries in Korea using disability-adjusted life years. J Korean Med Sci. 2019; 34(Suppl 1):e75. PMID: 30923488.
4. Ambrose AF, Paul G, Hausdorff JM. Risk factors for falls among older adults: a review of the literature. Maturitas. 2013; 75(1):51–61. PMID: 23523272.
5. World Health Organization. Ambient (Outdoor) Air Pollution. Geneva, Switzerland: World Health Organization;2021.
6. Kim KH, Kabir E, Kabir S. A review on the human health impact of airborne particulate matter. Environ Int. 2015; 74:136–143. PMID: 25454230.
7. Shin J, Han SH, Choi J. Exposure to ambient air pollution and cognitive impairment in community-dwelling older adults: the Korean Frailty and Aging Cohort Study. Int J Environ Res Public Health. 2019; 16(19):3767. PMID: 31591354.
8. Scapellato ML, Lotti M. Short-term effects of particulate matter: an inflammatory mechanism? Crit Rev Toxicol. 2007; 37(6):461–487. PMID: 17661213.
9. Merlo J, Wagner P, Leckie G. A simple multilevel approach for analysing geographical inequalities in public health reports: the case of municipality differences in obesity. Health Place. 2019; 58:102145. PMID: 31195211.
10. Oh TK, Song IA. Exposure to air pollution and risk of hip fracture: a population-based cohort study with a 6-year follow-up in South Korea. J Occup Environ Med. 2020; 62(12):1034–1039. PMID: 33065728.
11. Ha KH, Cho J, Cho SK, Kim C, Shin DC. Air pollution and unintentional injury deaths in South Korea. Environ Sci Pollut Res Int. 2015; 22(10):7873–7881. PMID: 25598159.
12. Kang YW, Ko YS, Kim YJ, Sung KM, Kim HJ, Choi HY, et al. Korea community health survey data profiles. Osong Public Health Res Perspect. 2015; 6(3):211–217. PMID: 26430619.
13. Obesity: preventing and managing the global epidemic. Report of a WHO consultation. World Health Organ Tech Rep Ser. 2000; 894:i–ixii. PMID: 11234459.
14. Kim E, Shin MH, Yang JH, Ahn SK, Na BJ, Nam HS. Associations of regional-level perceived stress and depression with health-related quality of life in Korean adults: a multilevel analysis of 2017 Korea Community Health Survey data. Epidemiol Health. 2021; 43:e2021062. PMID: 34525500.
15. Austin PC, Merlo J. Intermediate and advanced topics in multilevel logistic regression analysis. Stat Med. 2017; 36(20):3257–3277. PMID: 28543517.
16. Talbot LA, Musiol RJ, Witham EK, Metter EJ. Falls in young, middle-aged and older community dwelling adults: perceived cause, environmental factors and injury. BMC Public Health. 2005; 5(1):86. PMID: 16109159.
17. Kelsey JL, Berry SD, Procter-Gray E, Quach L, Nguyen US, Li W, et al. Indoor and outdoor falls in older adults are different: the maintenance of balance, independent living, intellect, and Zest in the Elderly of Boston Study. J Am Geriatr Soc. 2010; 58(11):2135–2141. PMID: 20831726.
18. Fischer P, Hoek G, Brunekreef B, Verhoeff A, van Wijnen J. Air pollution and mortality in The Netherlands: are the elderly more at risk? Eur Respir J Suppl. 2003; 40:34s–38s. PMID: 12762572.
19. Kaufman JD, Adar SD, Barr RG, Budoff M, Burke GL, Curl CL, et al. Association between air pollution and coronary artery calcification within six metropolitan areas in the USA (the Multi-Ethnic Study of Atherosclerosis and Air Pollution): a longitudinal cohort study. Lancet. 2016; 388(10045):696–704. PMID: 27233746.
20. Adar SD, Klein R, Klein BE, Szpiro AA, Cotch MF, Wong TY, et al. Air Pollution and the microvasculature: a cross-sectional assessment of in vivo retinal images in the population-based multi-ethnic study of atherosclerosis (MESA). PLoS Med. 2010; 7(11):e1000372. PMID: 21152417.
21. Łatka P, Nowakowska D, Nowomiejska K, Rejdak R. How air pollution affects the eyes — a review. Ophthalmology Journal. 2018; 3(2):58–62.
22. Loane C, Pilinis C, Lekkas TD, Politis M. Ambient particulate matter and its potential neurological consequences. Rev Neurosci. 2013; 24(3):323–335. PMID: 23612538.
23. Liu H, Hu T. Impact of socioeconomic status and health risk on fall inequality among older adults. Health Soc Care Community. 2022.
24. Choi KC, Kim SA, Kim NR, Shin MH. Association between smoking and unintentional injuries among Korean adults. Chonnam Med J. 2018; 54(3):184–189. PMID: 30288375.
25. Ogliari G, Ryg J, Andersen-Ranberg K, Scheel-Hincke LL, Masud T. Association between body mass index and falls in community-dwelling men and women: a prospective, multinational study in the Survey of Health, Ageing and Retirement in Europe (SHARE). Eur Geriatr Med. 2021; 12(4):837–849. PMID: 33817765.
26. Trevisan C, Ripamonti E, Grande G, Triolo F, Ek S, Maggi S, et al. The association between injurious falls and older adults’ cognitive function: the role of depressive mood and physical performance. J Gerontol A Biol Sci Med Sci. 2021; 76(9):1699–1706. PMID: 33649763.
27. Vinik AI, Camacho P, Reddy S, Valencia WM, Trence D, Matsumoto AM, et al. Aging, diabetes, and falls. Endocr Pract. 2017; 23(9):1117–1139. PMID: 28704101.
28. Klein D, Nagel G, Kleiner A, Ulmer H, Rehberger B, Concin H, et al. Blood pressure and falls in community-dwelling people aged 60 years and older in the VHM&PP cohort. BMC Geriatr. 2013; 13(1):50. PMID: 23692779.
29. Ortolá R, García-Esquinas E, Galán I, Guallar-Castillón P, López-García E, Banegas JR, et al. Patterns of alcohol consumption and risk of falls in older adults: a prospective cohort study. Osteoporos Int. 2017; 28(11):3143–3152. PMID: 28725986.
30. Chen CM, Yoon YH. Usual alcohol consumption and risks for nonfatal fall injuries in the United States: results from the 2004–2013 National Health Interview Survey. Subst Use Misuse. 2017; 52(9):1120–1132. PMID: 28524713.
31. Zhang Z, Liu J, Wu Y, Yan G, Zhu L, Yu X. Multi-scale comparison of the fine particle removal capacity of urban forests and wetlands. Sci Rep. 2017; 7(1):46214. PMID: 28393900.
32. Monn C. Exposure assessment of air pollutants: a review on spatial heterogeneity and indoor/outdoor/personal exposure to suspended particulate matter, nitrogen dioxide and ozone. Atmos Environ. 2001; 35(1):1–32.
33. Vardell E. Global health observatory data repository. Med Ref Serv Q. 2020; 39(1):67–74. PMID: 32069199.
TOOLS
Similar articles