Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Public Health, 30 September 2025

Sec. Environmental Health and Exposome

Volume 13 - 2025 | https://doi.org/10.3389/fpubh.2025.1661937

Long-term exposure to fine particulate matter, ozone, and greenness and the risk of lung cancer: a retrospective cohort analysis within a national sample cohort

Nam Eun Kim

&#x;
Nam Eun Kim1 Ji-Young Lee

&#x;
Ji-Young Lee2 Ga Young Lee
Ga Young Lee1 Chai Young Lee
Chai Young Lee1 Sojung Park
Sojung Park1 Yon Ju Ryu
Yon Ju Ryu1 Jin Hwa Lee



Jin Hwa Lee1 *
  • 1Division of Pulmonary and Critical Care Medicine, Department of Internal Medicine, College of Medicine, Ewha Womans University, Seoul, Republic of Korea
  • 2Inflammation-Cancer Microenvironment Research Center, College of Medicine, Ewha Womans University, Seoul, Republic of Korea

Background and objective: Given the rising incidence of lung cancer among never smokers and growing concerns about environmental risk factors, this study investigated the association between long-term exposure to air pollution and greenness and the risk of lung cancer.

Methods: We conducted a retrospective cohort analysis using data from the Korean National Sample Cohort (2002–2019), including 7,155 lung cancer patients and 28,620 propensity score-matched controls (matched by age, sex, and enrollment year). Long-term exposure to air pollution (quantified by PM2.5 and O₃ concentrations) and greenness (quantified by the normalized difference vegetation index, NDVI) was estimated based on residential area. Cox proportional hazards models were used to assess associations between exposure and lung cancer risk.

Results: Among 35,775 participants, lung cancer patients had lower BMI, higher smoking exposure, lower household income, and higher comorbidity scores than controls. PM2.5 exposure showed a modest association with increased lung cancer risk in the highest tertile (aHR = 1.06; 95% CI = 1.01–1.13). O₃ exposure was consistently associated with elevated risk across all tertiles (aHR = 1.42; 95% CI = 1.34–1.50). Greenness exposure demonstrated a protective effect (aHR = 0.89; 95% CI = 0.86–0.91). Subgroup analyses indicated that PM2.5 effects were more pronounced among male never smokers, O₃ exposure was associated with higher risk in female never smokers and males overall, and NDVI showed protective associations across all subgroups.

Conclusion: Long-term exposure to air pollution, particularly O₃, was significantly associated with an increased risk of lung cancer, independent of other confounders. In contrast, PM2.5 showed only a modest and inconsistent association, while high greenness exposure demonstrated a protective effect. These findings emphasize the need for region-specific environmental policies aimed at improving air quality and enhancing access to green spaces to reduce lung cancer risk.

Introduction

Particulate matter (PM) is classified by aerodynamic diameter, which determines how deeply particles penetrate into the respiratory tract. Fine particles (PM2.5; ≤2.5 μm) can reach the alveoli and enter systemic circulation (1, 2). Long term exposure to PM2.5 adversely affects multiple organ systems (3), contributing to respiratory diseases such as chronic obstructive pulmonary disease (4), cardiovascular disease (5), neurodegenerative disorders (6). These health effects are mediated by biological mechanisms including oxidative stress, inflammation, DNA damage, and epigenetic alterations (5, 79).

Ozone (O₃), another major ambient air pollutant, is a highly reactive gas with well-documented adverse effects on both respiratory and cardiovascular health (10, 11). Large cohort studies have demonstrated that long-term exposure is associated with increased respiratory mortality, with a 10 ppb rise in ozone concentration conferring approximately a 4% higher risk of death from respiratory causes, independent of PM2.5 exposure (12). In addition, ozone exposure has been linked to new-onset asthma in children and exacerbation of asthma symptoms in affected individuals (13). The underlying biological mechanisms involve oxidative stress and the generation of reactive oxygen species, which damage DNA, impair antioxidant defenses, and induce chronic airway inflammation and epithelial injury (14, 15). Collectively, these processes may contribute to cellular proliferation, mutagenesis, and the initiation of lung carcinogenesis (16).

Multiple epidemiological studies have demonstrated a significant association between exposure to PM2.5 and elevated lung cancer risk, particularly in cases with a higher concentration and longer duration of exposure. In the European ESCAPE study, each 5 μg/m3 increase in PM2.5 concentrations was associated with an 18% increase in lung cancer risk (17). A large U. S. cohort study similarly reported increased rates of lung adenocarcinoma among never-smokers exposed to PM2.5 (18). In another long-term investigation with a median follow-up of 10 years, PM2.5 exposures was associated with a 12% increase in lung cancer risk (19). Meta-analyses further support these findings, consistently demonstrating that exposure to PM2.5 or PM₁₀ is linked to elevated lung cancer risk (20). Evidence from East Asia aligns with these observations; large-scale cohort studies have reported significant associations between PM2.5 exposures and both lung cancer incidence and mortality (21). A nationwide Chinese study also demonstrated significant associations between ambient PM2.5 concentrations and cause-specific mortality, including deaths from lung cancer (22). In Korea, an NHIS-based cohort study further reported elevated lung cancer mortality among individuals exposed to ozone alone or in combination with PM2.5 with odds ratios ranging from 1.15 to 1.27 (23).

Based on accumulating evidence, the International Agency for Research on Cancer (IARC), a specialized agency of the World Health Organization (WHO), classified outdoor air pollution as a Group 1 carcinogen in 2013, indicating sufficient evidence of carcinogenicity in humans (24). Notably particulate matter, a major component of outdoor air pollution, was evaluated separately and was also classified as carcinogenic to humans. Subsequent longitudinal cohort studies have strengthened this conclusion by providing temporal evidence of the health burden attributable to air pollution. For example, examined temporal trends in lung cancer mortality attributable to PM2.5 exposures in China over a 30-year span using age-period-cohort analysis, demonstrating increasing burdens in older cohorts (25). Another investigation provided longitudinal insights into the health benefits of greenness, showing that reduced mortality was partly mediated by decreases in PM2.5 and NO₂ exposures, thereby highlighting the complex and time-varying interactions among environmental factors (26).

Although O3 has been classified by the International Agency for Research on Cancer (IARC) as Group 3, indicating that its carcinogenicity in humans is not classifiable, emerging evidence suggests a potential association between O₃ exposure and lung cancer risk. Long-term exposure to ambient O₃ has been linked to lung tissue injury and chronic airway inflammation, processes that may increase the susceptibility to various pulmonary diseases, including cancer (14, 27). In addition, several studies have suggested that O₃ may contribute to carcinogenesis when combined with PM2.5, through mechanisms involving accelerated lung function decline and enhanced oxidative stress (28).

In this study, we aimed to analyze the effects of long-term exposure to two air pollutants (PM2.5 and O3) on the development of lung cancer using nationwide data from the Korean National Health Insurance Service (NHIS) between 2002 to 2019. We further assessed the potential protective effects of residential greenness and examined whether the associations of air pollutants with lung cancer differed according to smoking status and sex. We hypothesized that long-term exposure to PM2.5 and O₃ would be associated with an increased risk of lung cancer, whereas residential greenness would be associated with a reduced risk.

Materials and methods

Study population

This study conducted a retrospective analysis utilizing the National Health Insurance Service-National Sample Cohort (NHIS-NSC) database, a comprehensive dataset managed by the NHIS, which encompasses a representative sample of the Korean population. The NHIS in Korea maintains records of all covered inpatient and outpatient visits, procedures, and prescriptions. The NHIS established the target population using the National Health Information Database (NHID) in 2002 and created the NHIS-National Sample Cohort (NHIS-NSC) by randomly selecting a representative sample of 1,137,896 individuals, corresponding to approximately 2.2% of the eligible Korean population at that time. The NHIS-NSC is a nationwide, retrospective cohort that spans from 2002 to 2019.

The study population comprised adults aged 20 years or older who had undergone at least one national health examination and for whom data on smoking status were available. The primary endpoint of this study was the occurrence of lung cancer identified during the follow-up period. Lung cancer cases were classified based on the International Classification of Diseases, 10th Revision (ICD-10) code C34.x. To reduce the possibility of reverse causation, individuals with a lung cancer diagnosis before 2003 were excluded from analysis. Incident cases were defined as those with a first recorded diagnosis of lung cancer between 2004 and 2019. Cases were defined as individuals receiving a first-time diagnosis of lung cancer within the study period. For the comparison group, control subjects without lung cancer were randomly selected and matched to cases in a 1:4 ratio using propensity scores. Matching variables included age, sex, and year of cohort entry. Participants were tracked until the earliest of lung cancer diagnosis, death, or the end of 2019. Participants were censored at the time of death or at the end of the follow-up period if they did not develop lung cancer.

This study was approved by the Institutional Review Board (IRB) of Ewha Womans University Medical Center, Seoul, Republic of Korea (IRB number: SEUMC2021-08-003). The IRB waived the need to obtain informed consent considering the retrospective nature of the study. All procedures were conducted in accordance with the relevant guidelines and regulations outlined in the latest revision of the Declaration of Helsinki.

Air pollution and green space exposure

The exposure variables used in this study were 5-year average concentrations of PM2.5 and O3 and the normalized difference vegetation index (NDVI). PM2.5 and O3 concentrations were estimated using a satellite-based spatiotemporal model based on aerosol optical depth data from the National Aeronautics and Space Administration (NASA). Estimates were calculated at a spatial resolution of 1 km × 1 km for each participant based on their residential address.

To determine the level of exposure to greenness, we used the NDVI, a satellite-derived metric that reflects vegetation density and plant health. NDVI data were collected from the Moderate Resolution Imaging Spectroradiometer (MODIS), as well as the Landsat 7 and 8 collections provided by the United States Geological Survey. The MODIS is the primary sensor for ground surface monitoring and is mounted on the Earth observation satellites Terra and Aqua. Its data are widely used to examine green space and other environmental factors. NDVI values were averaged for each participant based on the date of enrollment in the cohort. Long-term exposure to air pollution and green space was defined as the five-year average at participants’ residential addresses prior to cohort enrollment.

Statistical analysis

Descriptive statistics are expressed as the mean (standard deviation) for continuous variables and number (%) for categorical variables. Differences between groups were analyzed using the t-test for continuous variables and chi-square test for categorical variables. Kaplan–Meier analysis and the log-rank test were used to evaluate differences in lung cancer incidence across exposure categories. Cox proportional hazards models were applied to estimate the effects of air pollution and greenness on lung cancer risk by adjusting for confounding factors such as age, sex, smoking status, household income, residential area, body mass index (BMI), and Charlson Comorbidity Index (CCI). Subgroup analyses were performed to evaluate the associations of PM2.5, O3, and green space exposure with lung cancer risk according to gender and smoking status, and the impact of O3 was further examined by residential area. Adjusted covariates were selected based on established or suspected confounders identified in previous literature and known risk factors for the outcomes. Observations with missing values for any variable were excluded to preserve data integrity and ensure the validity of the results.

To assess the potential nonlinear association between exposure and lung cancer risk, we modeled the relationship using restricted cubic splines within the Cox proportional hazards framework. Hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated across the full range of exposure. To enhance the precision of individual exposure assessment to air pollution and green space, we conducted a sensitivity analysis. For workers, who may be exposed to environmental factors at both their workplace and residence, these exposures were explicitly considered. Using data from the NHIS, we compared lung cancer risk by classifying health insurance subscribers into workplace-based and community-based groups. Hazard ratios (HRs) with 95% confidence intervals (CIs) were calculated for lung cancer risk. SAS version 9.4 (SAS Institute, Cary, NC, United States) was used for large-scale data management, and R software version 4.0.3 (R Foundation for Statistical Computing, Vienna, Austria) was used for data analysis. The primary R packages utilized were survival (v3.2-7), survminer (v0.4.8), splines (included in base R 4.0.3), dplyr (v1.0.2), and ggplot2 (v3.3.2), which facilitated comprehensive survival analysis, data processing, and high-quality visualization. The significance level was set at p < 0.05.

Results

We identified 672,951 individuals aged ≥20 years who underwent health examinations between 2002 and 2019 in the NHIS-NSC database. After excluding participants younger than 20 years, those diagnosed with lung cancer or who died during the 2002–2003 washout period, and those with missing smoking or BMI data, 7,155 patients with lung cancer remained. Using 1:4 propensity score matching, we selected 28,620 controls without lung cancer, yielding a final analytic cohort of 35,775 participants (Figure 1). The standardized mean differences for the propensity score matching variables, including age, sex, and enrollment year, were all below 0.1, indicating adequate covariate balance between the groups. A Love plot illustrating the covariate balance before and after matching is presented in the Supplementary Figure S1. Time-to-event analyses were conducted using a stratified Cox proportional hazards regression model based on the propensity score matching. The matched set ID was specified as strata, allowing the baseline hazard to vary across matched sets.

Figure 1
Flowchart detailing participant selection in a study from NHIS-NSC database (2002-2019). Starts with 1,137,896 individuals. Excludes 464,945 for being under 20, having lung cancer diagnosis or death in 2002-2003, or lacking smoking data. Leaves 672,951 eligible subjects. 7,155 diagnosed with lung cancer as cases. 665,796 not diagnosed; 28,620 selected as control after excluding 637,176 via propensity score matching by age, sex, and enrollment year.

Figure 1. Flow diagram for selection of the study population. This figure illustrates the selection of study subjects from the NHIS-NSC database (2002–2019). After applying exclusion criteria and propensity score matching, the final cohort comprised 7,155 lung cancer cases and 28,620 matched controls. NHIS-NSC, National Health Insurance Service–National Sample Cohort.

Baseline characteristics are summarized in Table 1. Compared with controls, patients with lung cancer had lower mean BMI, were more often current smokers with higher pack-years, and had lower household income, higher comorbidity scores, and were more likely to live in rural areas, and have higher CCI scores. The mean (standard deviation, SD) exposure levels for PM2.5, O₃, and greenness (NDVI) were 29.01 (2.42) μg/m3, 34.79 (3.41) ppb, and 0.16 (0.11), respectively (Supplementary Table S1). Correlation analysis of PM2.5, O3, and NDVI showed that PM2.5 and O3 showed a weak positive correlation, and PM2.5 and NDVI showed a weak negative correlation. There was almost no correlation between O3 and NDVI (Supplementary Figure S2).

Table 1
www.frontiersin.org

Table 1. Baseline characteristics.

Significant associations were identified between lung cancer risk and exposures to air pollution and greenness (Figure 2). After adjusting for major covariates, O3 was associated with increased lung cancer risk (aHR = 1.612; 95% CI = 1.297–1.978), and the NDVI demonstrated a protective effect (aHR = 0.885; 95% CI = 0.858–0.912). When exposures were categorized into tertiles, PM2.5 was significantly associated with increased lung cancer risk in the 3rd tertile (aHR = 1.064; 95% CI = 1.006–1.127). O3 showed a positive association across all tertiles (aHR = 1.421; 95% CI = 1.342–1.503). A high NDVI indicated significant protective effects in the 3rd tertile (aHR = 0.709; 95% CI = 0.658–0.764). Hazard ratios for exposure by tertiles of air pollution and NDVI are provided in the Supplementary Table S2.

Figure 2
Three graphs labeled A, B, and C display hazard ratios with 95% confidence intervals. Graph A shows PM2.5 tertiles, Graph B shows O3 tertiles, and Graph C shows NDVI tertiles. The x-axes indicate increases or tertiles, while the y-axes represent hazard ratios. Horizontal red lines mark a ratio of one. Each graph compares different tertiles, indicating changes in hazard ratios across environmental factors.

Figure 2. Adjusted hazard ratios for lung cancer risk according to tertiles of PM2.5, O₃, and NDVI. This figure shows adjusted hazard ratios (aHRs) and 95% confidence intervals (CIs) for lung cancer risk according to tertiles of (A) PM2.5, (B) O₃, and (C) NDVI. “All” represents the overall association in the total study population. “T1 (Ref.)” indicates the reference group (lowest tertile), while T2 and T3 indicate the middle and highest tertiles, respectively. The red horizontal line denotes aHR = 1.0 (the null value). Adjusted hazard ratios were adjusted for age, sex, smoking status, body mass index, household income level, residential area, and Charlson comorbidity index. aHR, adjusted hazard ratio; BMI, body mass index; CCI, Charlson Comorbidity Index; NDVI, normalized difference vegetation index.

Kaplan–Meier analysis demonstrated higher cumulative incidence of lung cancer with high O₃ exposure and with low NDVI (both p < 0.001) (Figure 3). In subgroup analyses, PM2.5 exposure was associated with increased lung cancer risk particularly among male never smokers and ex-smokers. O3 exposure showed a significant association with lung cancer risk among all male subjects and female never smokers, and the NDVI demonstrated the protective effects of greenness exposure regardless of sex and smoking status (Table 2). Results stratified by IQR increases in air pollution and NDVI, respectively, are presented in the Supplementary Table S3. There is a trend toward an increased risk of lung cancer associated with ozone exposure in rural areas, and the difference in effect between regions is of borderline statistical significance (Supplementary Figure S3).

Figure 3
Three Kaplan-Meier curves show cumulative events over time based on exposure levels. Graph A compares low and high levels of PM2.5, with a p-value of 0.097. Graph B compares low and high levels of O3, with a significant p-value of less than 0.0001. Graph C compares low and high levels of NDVI, also with a significant p-value of less than 0.0001. Each graph includes a number at risk table over days.

Figure 3. Cumulative incidence of lung cancer according to levels of PM2.5, O₃, and NDVI. Cumulative incidence curves for lung cancer according to low (blue) and high (red) levels of (A) PM2.5, (B) O₃, and (C) NDVI. The y-axis represents the cumulative event probability, and the x-axis shows follow-up time (days). Numbers at risk at each time point are indicated below each panel. p-values compare incidence between exposure groups. NDVI, Normalized Difference Vegetation Index.

Table 2
www.frontiersin.org

Table 2. Subgroup analysis of lung cancer risk associated with air pollution and greenness exposure stratified by sex and smoking status.

The dose–response analysis demonstrated a nonlinear association between exposure and lung cancer risk. To further elucidate this relationship, we applied a restricted cubic spline model, and the resulting spline-based hazard ratio curve is presented in the Supplementary Figure S4. The curve indicates that the HR increases above a certain exposure level. The analysis of air pollution and green space exposure in relation to lung cancer risk between workplace-based and community-based groups yielded comparable results, with no statistically significant interactions observed upon inclusion of interaction terms in the model (Supplementary Table S4).

Discussion

Overall, long-term PM2.5 exposure was not significantly associated with lung cancer risk; however, participants in the highest tertile of exposure showed a 6.4% higher risk compared with those in the lowest tertile. In contrast, O₃ exposure was associated with a 61.2% increased risk of lung cancer, even after adjustment for potential confounding factors. Exposure to residential greenness, measured by the NDVI, was associated with a 12% lower risk of lung cancer, with the strongest protective effects observed in the highest tertile. A clear dose–response relationship was observed for PM2.5, O₃, and greenness, with more pronounced effects at higher exposure levels.

Subgroup analyses indicated that never-smoking men were particularly susceptible to PM2.5 exposure. This may reflect biological factors, such as sex-specific differences in pulmonary responses or the lack of adaptive mechanisms that could be present in smokers, although further research is needed to clarify this vulnerability (29). Previous studies have also reported that lung cancer risk among never smokers is influenced by ambient air pollution (30). For instance, the AHSMOG-2 cohort, predominantly composed of never smokers, demonstrated an elevated risk of lung cancer associated with PM2.5 exposure among individuals with long-term residence or greater time spent outdoors, indicating a dose–response relationship (31). In contrast, this association was not observed among current smokers, likely because the strong carcinogenic effect of smoking may mask the association with PM2.5 exposure. Consistent with this, a meta-analysis reported that the effect of PM2.5 was attenuated after adjusting for smoking status (32).

Although recent studies have reported significant associations between PM2.5 exposure and lung cancer development among female never smokers, this association was not observed in our study. This discrepancy may be explained by sex-specific differences in susceptibility and exposure levels. For example, one cohort study found that males were more susceptible to lung cancer at lower PM2.5 concentrations (0–35 μg/m3), whereas females demonstrated greater susceptibility at higher concentrations (35–75, 75–115, and 115–150 μg/m3) (33). Similarly, a study conducted in Taiwan reported that residential PM2.5 exposure above 30 μg/m3 was associated with an increased risk of lung adenocarcinoma among females (34). In our cohort, the mean PM2.5 concentration was 29.01 μg/m3, which falls within the range associated with increased susceptibility among males but not females.

PM2.5 has been classified as a Group 1 carcinogen by the International Agency for Research on Cancer (IARC), a specialized agency of the World Health Organization (WHO), with proposed mechanisms involving oxidative stress, chronic inflammation, DNA damage, and epigenetic alterations. A dose–response relationship between PM2.5 exposure and lung cancer risk has been consistently demonstrated in epidemiological studies (35, 36). Although evidence regarding O₃ remains inconclusive, our analyses showed that O₃ was more strongly associated with lung cancer risk than PM2.5, warranting consideration of several possible explanations.

First, subgroup analyses revealed that O₃ exposure was significantly associated with increased lung cancer risk in all male participants regardless of smoking status, and in female never and current smokers, but not in former smokers. Mechanistically, unlike PM2.5, which can reach the peripheral alveoli, O₃ primarily affects the airway epithelium (16), leading to localized oxidative injury and chronic inflammation. Smokers, who often have pre-existing airway inflammation, may therefore exhibit heightened susceptibility to O₃ exposure. In addition, this increased risk may also reflect behavioral factors, as males generally spend more time outdoors, thereby increasing cumulative O₃ exposure.

Second, recent trends suggest that O₃ exposure may pose a greater health risk compared with PM2.5 exposure. In Korea, as well as in North America and Europe, nationwide air pollution control policies have reduced the concentrations of major pollutants, including PM2.5 (37). According to data from Statistics Korea, mean PM2.5 levels decreased from 26.1 μg/m3 in 2015 to 23.6 μg/m3 in 2019 (38). In contrast, annual O₃ concentrations increased by approximately 13 ppb, corresponding to a 42% rise over a similar period, with higher levels in rural compared with urban areas (39). A Chinese study quantifying PM2.5-O₃ interactions similarly reported a 25.9% reduction in the health burden attributable to PM2.5 but an 11.8% increase in the burden attributable to O₃, primarily affecting cardiovascular, cerebrovascular, and respiratory diseases (40). Consistent with these findings, lung cancer incidence in our study was higher in rural areas. Although rural regions generally contain more green space, ambient O₃ levels are often elevated in such areas. In urban environments, nitric oxide (NO) emitted from traffic and industrial sources reacts with O₃, leading to reduced ambient concentrations. In contrast, in rural areas, O₃ precursors such as NO₂ can be transported by wind and undergo photochemical reactions more readily, resulting in higher O₃ concentrations (39).

Lastly, global warming has been suggested to contribute to a climate penalty effect, characterized by elevated O₃ formation and adverse health outcomes. Mechanistically, higher temperatures facilitate photochemical reactions that increase ambient O₃ concentrations. Biologically, chronic O₃ exposure is associated with sustained airway inflammation and oxidative stress, which may promote carcinogenesis through DNA damage (41). From a public health perspective, vulnerable populations—including the older adult(s), socioeconomically disadvantaged individuals, and those with limited access to cooling resources—may be at heightened risk of lung cancer during periods of elevated temperature and O₃ levels (42).

In contrast to the adverse effects of air pollution, green environments may act as a protective factor against lung cancer. In our study, higher NDVI values were associated with reduced lung cancer risk, independent of sex and smoking status. Consistent with these findings, a meta-analysis reported significant reductions in lung cancer incidence and mortality with greater exposure to greenness (43). Moreover, another large-scale meta-analysis demonstrated that green space exerts a protective effect on respiratory health through multiple pathways, including improved air quality, reduced heat exposure, alleviated stress and inflammation, increased physical activity, and enhanced immune function (44).

The protective role of greenness has also been observed in relation to air pollution–related outcomes. For example, a previous study reported that green space exposure was associated with reduced PM2.5-related mortality, with urban residents experiencing greater benefits than rural residents (45). In our study, however, the incidence of lung cancer was higher in rural populations. This finding underscores the need for region-specific greening strategies, particularly because rural areas often experience elevated O₃ levels due to long-range transport of precursors and enhanced photochemical reactions (38). Conversely, in urban settings, reductions in NOₓ emissions can paradoxically increase O₃ concentrations by reducing the scavenging of O₃ by freshly emitted NO, as illustrated in a case study from Zaragoza, Spain (46).

This study has several limitations. First, because the study population was restricted to Korea, regional differences in air pollutant composition, climate, and urban planning may limit the generalizability of our findings. For example, the toxicity of PM2.5 can vary depending on dominant emission sources such as coal combustion or traffic emissions (47). O₃ levels are influenced by meteorological and climatic conditions (48), while access to green spaces differs substantially across countries (49). Therefore, multi-region cohort studies and meta-analyses are needed to confirm the applicability of our results in diverse environmental contexts.

Second, exposure estimates for PM2.5, O₃, and greenness were derived from satellite-based models linked to residential addresses, which may not fully capture individual exposure variability, including workplace environments and indoor air quality.

Third, other air pollutants such as NO₂ were not considered, which could confound the observed associations. Future studies incorporating multi-pollutant models are needed to better reflect real-world exposure conditions and to disentangle the independent and interactive effects of multiple pollutants on lung cancer risk. Such approaches are critical for advancing environmental epidemiology and informing effective public health policies (50, 51).

Fourth, residual confounding cannot be excluded. Unmeasured factors such as occupational exposures, lifestyle characteristics, dietary factors, and comorbidities may have influenced the results.

Fifth, NDVI reflects vegetation density but does not account for actual accessibility, quality, or individual utilization of green spaces (52). Future research should consider alternative metrics, such as proximity-based or quality-adjusted measures, to more accurately assess individuals’ true exposure to green environments.

Conclusion

In conclusion, while PM2.5 has been classified as a Group 1 carcinogen by the International Agency for Research on Cancer (IARC), our findings indicate that long-term ozone exposure may represent an independent and underappreciated risk factor for lung cancer. Given anticipated changes in climate and atmospheric composition leading to rising and fluctuating ozone levels, ozone should be considered not only a short-term respiratory irritant but also a potential long-term contributor to lung carcinogenesis. Moreover, increasing access to urban green spaces may help mitigate lung cancer risk. Targeted early detection and prevention programs, particularly for high-risk groups such as never smokers and residents of high-ozone regions, may be warranted based on environmental exposure data. Finally, further research is needed to elucidate the biological pathways underlying ozone-related carcinogenesis, to examine interactions between PM2.5 and O₃, and to clarify the protective role of greenness across diverse populations and urban settings. Collectively, these findings provide evidence to inform global environmental interventions aimed at reducing the burden of lung cancer attributable to air pollution.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving humans were approved by IRB of Ewha Womans University Medical Center. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

NK: Conceptualization, Writing – original draft, Writing – review & editing. J-YL: Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Software, Writing – original draft. GL: Writing – review & editing. CL: Writing – review & editing. SP: Investigation, Writing – review & editing. YR: Conceptualization, Investigation, Writing – review & editing. JL: Conceptualization, Investigation, Project administration, Resources, Supervision, Validation, Visualization, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the National Research Foundation of Korea (NRF-2020R1A5A2019210 and NRF-2022R1I1A1A01068568).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2025.1661937/full#supplementary-material

References

1. Kim, K-H , Kabir, E , and Kabir, S . A review on the human health impact of airborne particulate matter. Environ Int. (2015) 74:136–43. doi: 10.1016/j.envint.2014.10.005

PubMed Abstract | Crossref Full Text | Google Scholar

2. Oberdörster, G , Oberdörster, E , and Oberdörster, J . Nanotoxicology: an emerging discipline evolving from studies of ultrafine particles. Environ Health Perspect. (2005) 113:823–39. doi: 10.1289/ehp.7339

PubMed Abstract | Crossref Full Text | Google Scholar

3. Sangkham, S , Phairuang, W , Sherchan, SP , Pansakun, N , Munkong, N , Sarndhong, K, et al. An update on adverse health effects from exposure to PM2.5. Environ Adv. (2024) 18:100603. doi: 10.1016/j.envadv.2024.100603

Crossref Full Text | Google Scholar

4. Yan, M , Ge, H , Zhang, L , Chen, X , Yang, X , Liu, F, et al. Long-term PM2.5 exposure in association with chronic respiratory diseases morbidity: a cohort study in northern China. Ecotoxicol Environ Saf. (2022) 244:114025. doi: 10.1016/j.ecoenv.2022.114025

PubMed Abstract | Crossref Full Text | Google Scholar

5. Rinaldi, R , Russo, M , Bonanni, A , Camilli, M , Caffè, A , Basile, M, et al. Short-term air pollution exposure and mechanisms of plaque instability in acute coronary syndromes: an optical coherence tomography study. Atherosclerosis. (2024) 390:117393. doi: 10.1016/j.atherosclerosis.2023.117393

PubMed Abstract | Crossref Full Text | Google Scholar

6. Lee, J , Weerasinghe-Mudiyanselage, PDE , Kim, B , Kang, S , Kim, J-S , and Moon, C . Particulate matter exposure and neurodegenerative diseases: a comprehensive update on toxicity and mechanisms. Ecotoxicol Environ Saf. (2023) 266:115565. doi: 10.1016/j.ecoenv.2023.115565

PubMed Abstract | Crossref Full Text | Google Scholar

7. Thangavel, P , Park, D , and Lee, YC . Recent insights into particulate matter (PM(2.5))-mediated toxicity in humans: an overview. Int J Environ Res Public Health. (2022) 19:7511. doi: 10.3390/ijerph19127511

PubMed Abstract | Crossref Full Text | Google Scholar

8. Morales-Rubio, R , Amador-Muñoz, O , Rosas-Pérez, I , Sánchez-Pérez, Y , García-Cuéllar, C , Segura-Medina, P, et al. PM(2.5) induces airway hyperresponsiveness and inflammation via the AhR pathway in a sensitized Guinea pig asthma-like model. Toxicology. (2022) 465:153026. doi: 10.1016/j.tox.2021.153026

PubMed Abstract | Crossref Full Text | Google Scholar

9. Zhao, J , Gao, Z , Tian, Z , Xie, Y , Xin, F , Jiang, R, et al. The biological effects of individual-level PM2.5 exposure on systemic immunity and inflammatory response in traffic policemen. Occup Environ Med. (2013) 70:426–31. doi: 10.1136/oemed-2012-100864

Crossref Full Text | Google Scholar

10. Turner, MC , Jerrett, M , Pope, CA 3rd, Krewski, D , Gapstur, SM , Diver, WR, et al. Long-term ozone exposure and mortality in a large prospective study. Am J Respir Crit Care Med. (2016) 193:1134–42. doi: 10.1164/rccm.201508-1633OC

Crossref Full Text | Google Scholar

11. Malley, CS , Henze, DK , Kuylenstierna, JCI , Vallack, HW , Davila, Y , Anenberg, SC, et al. Updated global estimates of respiratory mortality in adults ≥30Years of age attributable to long-term ozone exposure. Environ Health Perspect. (2017) 125:087021. doi: 10.1289/EHP1390

PubMed Abstract | Crossref Full Text | Google Scholar

12. Jerrett, M , Burnett, RT , Pope, CA 3rd, Ito, K , Thurston, G , Krewski, D, et al. Long-term ozone exposure and mortality. N Engl J Med. (2009) 360:1085–95. doi: 10.1056/NEJMoa0803894

Crossref Full Text | Google Scholar

13. Nuvolone, D , Petri, D , and Voller, F . The effects of ozone on human health. Environ Sci Pollut Res Int. (2018) 25:8074–88. doi: 10.1007/s11356-017-9239-3

PubMed Abstract | Crossref Full Text | Google Scholar

14. Wiegman, CH , Li, F , Ryffel, B , Togbe, D , and Chung, KF . Oxidative stress in ozone-induced chronic lung inflammation and emphysema: a facet of chronic obstructive pulmonary disease. Front Immunol. (2020) 11:1957. doi: 10.3389/fimmu.2020.01957

PubMed Abstract | Crossref Full Text | Google Scholar

15. Pryor, WA , and Squadrito, GL . The chemistry of peroxynitrite: a product from the reaction of nitric oxide with superoxide. Am J Phys. (1995) 268:L699–722. doi: 10.1152/ajplung.1995.268.5.L699

PubMed Abstract | Crossref Full Text | Google Scholar

16. Russo, RC , Togbe, D , Couillin, I , Segueni, N , Han, L , Quesniaux, VFJ, et al. Ozone-induced lung injury and inflammation: pathways and therapeutic targets for pulmonary diseases caused by air pollutants. Environ Int. (2025) 198:109391. doi: 10.1016/j.envint.2025.109391

PubMed Abstract | Crossref Full Text | Google Scholar

17. Raaschou-Nielsen, O , Andersen, ZJ , Beelen, R , Samoli, E , Stafoggia, M , Weinmayr, G, et al. Air pollution and lung cancer incidence in 17 European cohorts: prospective analyses from the European study of cohorts for air pollution effects (ESCAPE). Lancet Oncol. (2013) 14:813–22. doi: 10.1016/S1470-2045(13)70279-1

PubMed Abstract | Crossref Full Text | Google Scholar

18. Puett, RC , Hart, JE , Yanosky, JD , Spiegelman, D , Wang, M , Fisher, JA, et al. Particulate matter air pollution exposure, distance to road, and incident lung cancer in the nurses' health study cohort. Environ Health Perspect. (2014) 122:926–32. doi: 10.1289/ehp.1307490

PubMed Abstract | Crossref Full Text | Google Scholar

19. Huang, YJ , Lee, PH , Chen, LC , Lin, BC , Lin, C , and Chan, TC . Relationships among green space, ambient fine particulate matter, and cancer incidence in Taiwan: a 16-year retrospective cohort study. Environ Res. (2022) 212:113416. doi: 10.1016/j.envres.2022.113416

PubMed Abstract | Crossref Full Text | Google Scholar

20. Ciabattini, M , Rizzello, E , Lucaroni, F , Palombi, L , and Boffetta, P . Systematic review and meta-analysis of recent high-quality studies on exposure to particulate matter and risk of lung cancer. Environ Res. (2021) 196:110440. doi: 10.1016/j.envres.2020.110440

PubMed Abstract | Crossref Full Text | Google Scholar

21. Liu, X , Mubarik, S , Wang, F , Yu, Y , Wang, Y , Shi, F, et al. Lung cancer death attributable to long-term ambient particulate matter (PM2.5) exposure in east Asian countries during 1990–2019. Front Med Lausanne. (2021) 8:742076. doi: 10.3389/fmed.2021.742076

Crossref Full Text | Google Scholar

22. Yang, J , Zhou, M , Li, M , Yin, P , Hu, J , Zhang, C, et al. Fine particulate matter constituents and cause-specific mortality in China: a nationwide modelling study. Environ Int. (2020) 143:105927. doi: 10.1016/j.envint.2020.105927

PubMed Abstract | Crossref Full Text | Google Scholar

23. Jin, T , Lee, S , Seo, J , Ye, S , Kim, S , Oh, JK, et al. Long-term ambient ozone exposure and lung cancer mortality: a nested case-control study in Korea. Environ Pollut. (2025) 375:126299. doi: 10.1016/j.envpol.2025.126299

PubMed Abstract | Crossref Full Text | Google Scholar

24. International Agency for Research on Cancer . IARC. Lyon: Outdoor air pollution a leading environmental cause of cancer deaths International Agency for Research on Cancer (2011).

Google Scholar

25. Wu, X , Zhu, B , Zhou, J , Bi, Y , Xu, S , and Zhou, B . The epidemiological trends in the burden of lung cancer attributable to PM2.5 exposure in China. BMC Public Health. (2021) 21:737. doi: 10.1186/s12889-021-10765-1

PubMed Abstract | Crossref Full Text | Google Scholar

26. Wu, C , Liu, J , Li, Y , Qin, L , Gu, R , Feng, J, et al. Association of residential air pollution and green space with all-cause and cause-specific mortality in individuals with diabetes: an 11-year prospective cohort study. EBioMedicine. (2024) 108:105376. doi: 10.1016/j.ebiom.2024.105376

PubMed Abstract | Crossref Full Text | Google Scholar

27. Guo, H , Liu, J , and Wei, J . Ambient ozone, PM(1) and female lung cancer incidence in 436 Chinese counties. Int J Environ Res Public Health. (2021) 18:10386. doi: 10.3390/ijerph181910386

PubMed Abstract | Crossref Full Text | Google Scholar

28. Valavanidis, A , Vlachogianni, T , Fiotakis, K , and Loridas, S . Pulmonary oxidative stress, inflammation and cancer: respirable particulate matter, fibrous dusts and ozone as major causes of lung carcinogenesis through reactive oxygen species mechanisms. Int J Environ Res Public Health. (2013) 10:3886–907. doi: 10.3390/ijerph10093886

PubMed Abstract | Crossref Full Text | Google Scholar

29. May, L , Shows, K , Nana-Sinkam, P , Li, H , and Landry, JW . Sex differences in lung cancer. Cancers (Basel). (2023) 15:3111. doi: 10.3390/cancers15123111

PubMed Abstract | Crossref Full Text | Google Scholar

30. Lim, JU , and Yoon, HK . Narrative review: association between lung cancer development and ambient particulate matter in never-smokers. J Thorac Dis. (2022) 14:553–63. doi: 10.21037/jtd-21-655

PubMed Abstract | Crossref Full Text | Google Scholar

31. Gharibvand, L , Shavlik, D , Ghamsary, M , Beeson, WL , Soret, S , Knutsen, R, et al. The association between ambient fine particulate air pollution and lung cancer incidence: results from the AHSMOG-2 study. Environ Health Perspect. (2017) 125:378–84. doi: 10.1289/EHP124

PubMed Abstract | Crossref Full Text | Google Scholar

32. Hamra, GB , Guha, N , Cohen, A , Laden, F , Raaschou-Nielsen, O , Samet, JM, et al. Outdoor particulate matter exposure and lung cancer: a systematic review and meta-analysis. Environ Health Perspect. (2014) 122:906–11. doi: 10.1289/ehp/1408092

PubMed Abstract | Crossref Full Text | Google Scholar

33. Ma, Z , Meng, X , Chen, C , Chao, B , Zhang, C , and Li, W . Short-term effects of different PM2.5 ranges on daily all-cause mortality in Jinan, China. Sci Rep. (2022) 12:5665. doi: 10.1038/s41598-022-09057-4

PubMed Abstract | Crossref Full Text | Google Scholar

34. Yang, SC , Lin, FY , Wu, TI , Wu, CD , and Wang, JD . PM(2.5) exposure and risk of lung adenocarcinoma in women of Taiwan: a case-control study with density sampling. Respirology. (2022) 27:951–8. doi: 10.1111/resp.14316

PubMed Abstract | Crossref Full Text | Google Scholar

35. Huang, F , Pan, B , Wu, J , Chen, E , and Chen, L . Relationship between exposure to PM2.5 and lung cancer incidence and mortality: a meta-analysis. Oncotarget. (2017) 8:43322–31. doi: 10.18632/oncotarget.17313

PubMed Abstract | Crossref Full Text | Google Scholar

36. Yoo, J , Lee, Y , Park, Y , Lee, J , Choi, JY , Lee, H, et al. Update in association between lung cancer and air pollution. Tuberc Respir Dis (Seoul). (2025) 88:228–36. doi: 10.4046/trd.2024.0092

PubMed Abstract | Crossref Full Text | Google Scholar

37. Amann, M , Kiesewetter, G , Schöpp, W , Klimont, Z , Winiwarter, W , Cofala, J, et al. Reducing global air pollution: the scope for further policy interventions. Philos Trans A Math Phys Eng Sci. (2020) 378:20190331. doi: 10.1098/rsta.2019.0331

PubMed Abstract | Crossref Full Text | Google Scholar

38. Kim, SW , Kim, KM , Jeong, Y , Seo, S , Park, Y , and Kim, J . Changes in surface ozone in South Korea on diurnal to decadal timescales for the period of 2001–2021. Atmos Chem Phys. (2023) 23:12867–86. doi: 10.5194/acp-23-12867-2023

Crossref Full Text | Google Scholar

39. Mongju, J , Woo Jin, K , Jongbae, H , and Sun-Young, K . Exploration of the spatial and temporal patterns of surface ozone concentrations for development of ozone prediction model in South Korea. J Korean Soc Atmos Environ. (2022) 38:100–25. doi: 10.5572/KOSAE.2022.38.1.100

Crossref Full Text | Google Scholar

40. Lu, Z , Guan, Y , Shao, C , and Niu, R . Assessing the health impacts of PM2.5 and ozone pollution and their comprehensive correlation in Chinese cities based on extended correlation coefficient. Ecotoxicol Environ Saf. (2023) 262:115125. doi: 10.1016/j.ecoenv.2023.115125

PubMed Abstract | Crossref Full Text | Google Scholar

41. Ware, LB , Zhao, Z , Koyama, T , May, AK , Matthay, MA , Lurmann, FW, et al. Long-term ozone exposure increases the risk of developing the acute respiratory distress syndrome. Am J Respir Crit Care Med. (2016) 193:1143–50. doi: 10.1164/rccm.201507-1418OC

PubMed Abstract | Crossref Full Text | Google Scholar

42. Bell, ML , Zanobetti, A , and Dominici, F . Who is more affected by ozone pollution? A systematic review and meta-analysis. Am J Epidemiol. (2014) 180:15–28. doi: 10.1093/aje/kwu115

PubMed Abstract | Crossref Full Text | Google Scholar

43. Li, J , Xie, Y , Xu, J , Zhang, C , Wang, H , Huang, D, et al. Association between greenspace and cancer: evidence from a systematic review and meta-analysis of multiple large cohort studies. Environ Sci Pollut Res. (2023) 30:91140–57. doi: 10.1007/s11356-023-28461-5

PubMed Abstract | Crossref Full Text | Google Scholar

44. Tang, M , Liu, W , Li, H , and Li, F . Greenness and chronic respiratory health issues: a systematic review and meta-analysis. Front Public Health. (2023) 11:1279322. doi: 10.3389/fpubh.2023.1279322

PubMed Abstract | Crossref Full Text | Google Scholar

45. Ji, JS , Zhu, A , Lv, Y , and Shi, X . Interaction between residential greenness and air pollution mortality: analysis of the Chinese longitudinal healthy longevity survey. Lancet Planet Health. (2020) 4:e107–15. doi: 10.1016/S2542-5196(20)30027-9

PubMed Abstract | Crossref Full Text | Google Scholar

46. Escudero, M , Lozano-Rojo, Á , Hierro, J , del Valle, J , and Mantilla, E . Urban influence on increasing ozone concentrations in a characteristic Mediterranean agglomeration. Atmos Environ. (2014) 99:322–32. doi: 10.1016/j.atmosenv.2014.09.061

Crossref Full Text | Google Scholar

47. Lelieveld, J , Evans, JS , Fnais, M , Giannadaki, D , and Pozzer, A . The contribution of outdoor air pollution sources to premature mortality on a global scale. Nature. (2015) 525:367–71. doi: 10.1038/nature15371

PubMed Abstract | Crossref Full Text | Google Scholar

48. Staten, PW , Lu, J , Grise, KM , Davis, SM , and Birner, T . Re-examining tropical expansion. Nat Clim Chang. (2018) 8:768–75. doi: 10.1038/s41558-018-0246-2

Crossref Full Text | Google Scholar

49. Nieuwenhuijsen, MJ , Khreis, H , Triguero-Mas, M , Gascon, M , and Dadvand, P . Fifty shades of green: pathway to healthy urban living. Epidemiology. (2017) 28:63–71. doi: 10.1097/EDE.0000000000000549

PubMed Abstract | Crossref Full Text | Google Scholar

50. Dominici, F , Peng, RD , Barr, CD , and Bell, ML . Protecting human health from air pollution: shifting from a single-pollutant to a multipollutant approach. Epidemiology. (2010) 21:187–94. doi: 10.1097/EDE.0b013e3181cc86e8

PubMed Abstract | Crossref Full Text | Google Scholar

51. Mauderly, JL , and Samet, JM . Is there evidence for synergy among air pollutants in causing health effects? Environ Health Perspect. (2009) 117:1–6. doi: 10.1289/ehp.11654

PubMed Abstract | Crossref Full Text | Google Scholar

52. Markevych, I , Schoierer, J , Hartig, T , Chudnovsky, A , Hystad, P , Dzhambov, AM, et al. Exploring pathways linking greenspace to health: theoretical and methodological guidance. Environ Res. (2017) 158:301–17. doi: 10.1016/j.envres.2017.06.028

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: particulate matter, ozone, greenness, lung cancer risk, air pollution

Citation: Kim NE, Lee J-Y, Lee GY, Lee CY, Park S, Ryu YJ and Lee JH (2025) Long-term exposure to fine particulate matter, ozone, and greenness and the risk of lung cancer: a retrospective cohort analysis within a national sample cohort. Front. Public Health. 13:1661937. doi: 10.3389/fpubh.2025.1661937

Received: 08 July 2025; Accepted: 09 September 2025;
Published: 30 September 2025.

Edited by:

Tong Wang, Duke University, United States

Reviewed by:

Worradorn Phairuang, Chiang Mai University, Thailand
Panpan Gu, Fudan University, China
Zixu Wang, Bicycle Therapeutics, United Kingdom

Copyright © 2025 Kim, Lee, Lee, Lee, Park, Ryu and Lee. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jin Hwa Lee, amluaHdhbGVlQGV3aGEuYWMua3I=

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.