Global trends in the burden of tracheal, bronchus, and lung cancer attributable to household air pollution
Original Article

Global trends in the burden of tracheal, bronchus, and lung cancer attributable to household air pollution

Shen Lao, Jianfu Li, Wangzhong Li, Hengrui Liang, Wenhua Liang, Fuli Li, Wei Wang

Department of Thoracic Surgery and Oncology, The First Affiliated Hospital of Guangzhou Medical University, State Key Laboratory of Respiratory Disease and National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Healthy, Guangzhou, China

Contributions: (I) Conception and design: S Lao, W Wang, W Liang; (II) Administrative support: F Li, W Wang, W Liang; (III) Provision of study materials or patients: None; (IV) Collection and assembly of data: J Li, W Li, H Liang; (V) Data analysis and interpretation: S Lao, J Li, W Li; (VI) Manuscript writing: All authors; (VII) Final approval of manuscript: All authors.

Correspondence to: Wei Wang, MD; Fuli Li, MD. Department of Thoracic Surgery and Oncology, The First Affiliated Hospital of Guangzhou Medical University, State Key Laboratory of Respiratory Disease and National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Healthy, 151 Yanjiang Rd., Guangzhou 510030, China. Email: bbmcwei@126.com; 337399983@qq.com.

Background: Tracheal, bronchial, and lung (TBL) cancers, largely attributable to household air pollution (HAP), remain a critical public health challenge, particularly in low- and middle-income countries (LMICs). Using data from the Global Burden of Disease (GBD) 2021, this study evaluates global trends in HAP-related TBL cancer burden from 1990 to 2021 and projects future trends to 2051.

Methods: The study extracted data on deaths and disability-adjusted life years (DALYs) from the GBD 2021 and assessed the burden of TBL cancer attributable to HAP by sex, age, region, and socio-demographic index (SDI). The age-period-cohort (APC) model examined the independent effects of age, period, and cohort, while the Bayesian age-period-cohort (BAPC) model predicted future trends over the next 30 years.

Results: Between 1990 and 2021, deaths and DALYs from TBL cancer due to HAP decreased by 39.2% and 44.1%, respectively, with the largest declines in high SDI regions, while rates increased in Oceania and sub-Saharan Africa. The global age-standardized mortality rate (ASMR) peaked in the 85–89 years age group, and the age-standardized disability rate (ASDR) peaked in the 70–74 years age group. Males experienced a higher burden than females, with significant differences observed across SDI regions. Projections suggest that mortality and DALYs rates will rise over the next 30 years, with a higher burden forecast for females.

Conclusions: Despite significant reductions, TBL cancer from HAP remains a major issue in low SDI regions, particularly in China. Projections indicate increasing mortality and DALYs rates over the next 30 years, disproportionately affecting women, highlighting the urgent need for targeted policies to reduce HAP exposure.

Keywords: Household air pollution (HAP); tracheal, bronchial, and lung cancer (TBL cancer); Global Burden of Disease 2021 (GBD 2021); socio-demographic index (SDI)


Submitted Jun 28, 2025. Accepted for publication Sep 19, 2025. Published online Nov 26, 2025.

doi: 10.21037/jtd-2025-1309


Highlight box

Key findings

• Between 1990 and 2021, deaths and disability-adjusted life years (DALYs) from household air pollution (HAP)-related tracheal, bronchial, and lung (TBL) cancer decreased by 39.2% and 44.1%, respectively, with significant declines in high-socio-demographic index (SDI) regions but a worsening burden in Oceania and sub-Saharan Africa. Projections indicate a reversal of this trend, with rising global mortality and DALYs rates over the next 30 years, and a higher burden forecast for females.

What is known and what is new?

• It is known that HAP is a major risk factor for TBL cancer, with a geographically heterogeneous burden.

• This study provides a comprehensive 30-year trend analysis using Global Burden of Disease 2021 data and projects the burden to 2051. It newly identifies a predicted shift to increasing global rates, highlighting a disproportionately higher future risk for women.

What is the implication, and what should change now?

• The findings signal that current interventions are insufficient. Targeted policies to reduce HAP exposure are urgently needed, particularly in low-SDI regions (e.g., China), Oceania, and sub-Saharan Africa. Action should include strategies specifically designed to protect women, to counter the projected increasing trend.


Introduction

Household air pollution (HAP) is a leading environmental risk factor for chronic diseases, including tracheal, bronchial, and lung (TBL) cancer, particularly in low- and middle-income countries (LMICs), significantly contributing to global morbidity and mortality. Among the 3.2 million deaths attributed to HAP, 6% are due to lung cancer, with 11% of adult lung cancer deaths linked to exposure to carcinogenic pollutants from traditional fuels (1,2). Despite efforts to transition to cleaner energy sources, over 2.4 billion people worldwide continue to use polluting fuels, exacerbating health disparities, especially in regions with lower socio-economic development (3). The use of solid fuels, such as wood, coal, and agricultural waste, in poorly ventilated indoor environments exposes millions of individuals to high levels of harmful pollutants, including particulate matter (PM), carbon monoxide, and carcinogenic compounds such as benzene and formaldehyde.

However, the global, regional, and national epidemiological burden of TBL cancer due to HAP remains poorly understood. Most previous studies were limited to specific regions or periods and lacked comprehensive predictive models (4,5), which hinders the development of effective policy interventions. Given the rapid demographic and environmental changes expected in the coming decades, particularly in rapidly developing countries, there is an urgent need to provide a comprehensive assessment of the current and future burden of TBL cancer attributable to HAP.

International efforts spearheaded by agencies such as the World Health Organization (WHO) and the United Nations Children’s Fund (UNICEF) have been implemented to reduce HAP exposure and improve household air quality (6). These initiatives include worldwide programs promoting the adoption of clean energy and advanced cookstove technologies, which are designed to curb the use of solid fuels and mitigate associated health risks. Such interventions are vital for attaining the health-oriented objectives of the 2030 Sustainable Development Goals (SDGs), especially SDG 3 (to ensure healthy lives and promote well-being for all at all ages) and SDG 7 (to ensure access to affordable, reliable, sustainable, and modern energy (7). Achieving progress in these areas is crucial for alleviating the global disease burden related to HAP and enhancing overall quality of life.

This study used data from the Global Burden of Disease (GBD) 2021 to assess the global burden of TBL cancer from HAP between 1990 and 2021. It projects future trends up to 2051 using Bayesian age-period-cohort (BAPC) models. The analysis examines demographic and geographic variations, persistent inequalities in low SDI regions like China, and trends related to gender and aging. The findings underscore the need for targeted interventions to address energy poverty, gender inequalities, and regional disparities, aligning with sustainable development goals. This work builds on recent advancements in pollution-cancer epidemiology and provides a framework to reduce HAP’s health risks amidst evolving environmental and demographic challenges. We present this article in accordance with the STROBE reporting checklist (available at https://jtd.amegroups.com/article/view/10.21037/jtd-2025-1309/rc).


Methods

Data sources

The research leveraged data from the GBD 2021, conducted by the Institute for Health Metrics and Evaluation (IHME) (https://vizhub.healthdata.org/gbd) (4). The GBD 2021 provides comprehensive data on global health issues, including diseases, injuries, and risk factors, covering 204 countries and territories from 1990 to 2021. The study included 371 diseases and injuries, including TBL cancer (ICD-10 codes C33–C34) (8,9). This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. The data utilized in the study are publicly accessible and did not require additional ethical approval.

Disability-adjusted life years (DALYs)

Data on deaths and DALYs for TBL cancer from the GBD 2021. DALYs, which combine the years of life lost due to premature death and the years lived with disability, serve as a key metric for assessing disease burden and the impact of healthcare interventions (9-11). This metric is essential for evaluating the impact of TBL cancer and the effectiveness of health interventions, particularly in disease prevention and control.

Socio-demographic index (SDI)

The SDI is a composite indicator of national development, incorporating factors such as educational attainment, per capita income, and fertility rates (12,13). Scaled between 0 and 1, higher SDI values indicate greater socio-economic development. In GBD 2021, countries were grouped into five SDI quintiles, allowing for analysis of the relationship between the burden of TBL cancer attributable to HAP and socio-economic progress (14). Adaptive Loess regression was used to assess the correlation between SDI and age-standardized rates (ASRs) of TBL cancer, including age-standardized mortality rates (ASMRs) and age-standardized disability rates (ASDRs).

Age-period-cohort (APC) model and BAPC model

The APC model was employed to assess the independent effects of age, period, and cohort on the burden of TBL cancer attributable to HAP (14). Based on the Poisson distribution, this model estimates mortality or DALYs rates using the equation: log(λij) = µ + αi + βj + γk, where λij was the value for TBL cancer attributable to HAP, and µ, αi, βj, and γk were the mean effect, age, period, and cohort effects, respectively (15). Risk ratios (RRs) were calculated for period and cohort effects relative to a reference. The National Cancer Institute’s APC network tool was used to estimate these parameters. Additionally, the BAPC model, combined with Integrated Nested Laplace Approximations (INLA), was applied to project future trends in ASMR and ASDR for TBL cancer attributable to HAP over the next 30 years (16).

Data analysis

Temporal patterns and trends in TBL cancer attributable to HAP were assessed using annual percentage change and estimated annual percentage change (EAPC) from 1990 to 2021 (5). ASRs per 100,000 people were calculated using the formula (17).

ASR=i=1Aaiwii=1Awi×100000

Where ai represents the age-specific rate for the ith age group, Wi denotes the number of individuals within the corresponding age group, and A indicates the total number of distinct age groups.

The EAPC, derived from a regression model, represents the annual percentage change in ASR and is calculated as 100 × [exp(β) − 1]. An increasing trend was indicated when the EAPC and its 95% uncertainty interval (UI) lower bound were positive, while a decreasing trend was identified when the EAPC and its 95% UI upper bound were negative. If neither condition was met, the ASR was considered stable (18).

Data analysis was conducted using R software (version 4.4.1), with visualizations created using the ggplot2 package. Final figure modifications were completed using Adobe Illustrator.


Results

Global burden of TBL cancer attributable to HAP

From 1990 to 2021, there was a notable decline in the global burden of TBL cancer attributable to HAP. Mortality due to TBL cancer decreased by 39.2%, from 125,844 [95% uncertainty interval (UI): 77,843 to 177,114] to 76,482 (95% UI: 28,604 to 187,340). Concurrently, the ASMR fell by 71.97%, from 3.14 (95% UI: 1.94 to 4.43) to 0.88 (95% UI: 0.33 to 2.16) per 100,000 people (Table 1). Similarly, the number of DALYs decreased by 44.1%, from 3,517,878 (95% UI: 2,180,370 to 4,954,257) to 1,966,426 (95% UI: 758,405 to 4,632,387), with the ASDR declining by 73.37%, from 84.14 (95% UI: 52.12 to 118.39) to 22.41 (95% UI: 8.64 to 52.88) per 100,000 people (Table 2).

Table 1

Numbers and ASRs per 100,000 cases of deaths of TBL cancer attributable to HAP in 1990 and 2021, along with the relative changes and EAPC in ASRs per 100,000 cases from 1990 to 2021, categorized by global, SDI, and GBD regions

Location 1990 2021 1990 to 2021 EPAC (95% CI)
Number (95% CI) Age-standardized rate (95% UI) Number (95% CI) Age-standardized rate (95% UI) Relative change of numbers (%) Relative change of age-standardized rate (%)
Global 125,844
(77,843, 177,114)
3.14
(1.94, 4.43)
76,482
(28,604, 187,340)
0.88
(0.33, 2.16)
−39.22475 −71.97452 −4.6
(−5.16, −4.03)
High SDI 3,659
(1,109, 9,173)
0.33
(0.1, 0.84)
215
(0, 2,011)
0.01
(0, 0.1)
−94.12408 −96.9697 −12.15
(−12.67, −11.63)
High-middle SDI 42,140
(25,090, 65,176)
4.21
(2.51, 6.49)
8,765
(476, 44,605)
0.44
(0.02, 2.23)
−79.20028 −89.54869 −8.11
(−9.17, −7.03)
Middle SDI 59,497
(38,828, 81,149)
5.86
(3.84, 7.99)
34,338
(7,201, 96,083)
1.3
(0.27, 3.65)
−42.28617 −77.8157 −5.13
(−5.83, −4.41)
Low-middle SDI 15,585
(9,978, 21,466)
2.59
(1.66, 3.55)
23,813
(13,505, 37,394)
1.67
(0.95, 2.62)
52.794354 −35.52124 −1.55
(−1.73, −1.36)
Low SDI 4,865
(3,175, 6,934)
2.19
(1.43, 3.12)
9,294
(5,854, 13,002)
1.9
(1.21, 2.64)
91.038027 −13.24201 −0.59
(−0.65, −0.53)
Andean Latin America 420
(177, 748)
2.13
(0.89, 3.8)
206
(34, 672)
0.36
(0.06, 1.16)
−50.95238 −83.09859 −5.82
(−6.29, −5.35)
Australasia 5 (0, 48) 0.02 (0, 0.2) 1 (0, 4) 0
(0, 0.01)
−80 −100 −10
(−10.59, −9.4)
Caribbean 377
(190, 632)
1.46
(0.73, 2.47)
360
(180, 662)
0.67
(0.33, 1.23)
−4.509284 −54.10959 −2.62
(−2.69, −2.55)
Central Asia 1,151
(421, 2,545)
2.39
(0.88, 5.27)
346
(118, 890)
0.42
(0.14, 1.08)
−69.93918 −82.42678 −6.57
(−7.39, −5.74)
Central Europe 3,362
(878, 9,207)
2.21
(0.58, 6.04)
640
(19, 4,480)
0.29
(0.01, 2.03)
−80.96371 −86.87783 −7.64
(−8.29, −7)
Central Latin America 770
(288, 1,653)
0.97
(0.36, 2.09)
625
(216, 1,571)
0.25
(0.09, 0.64)
−18.83117 −74.2268 −4.56
(−4.81, −4.32)
Central sub-Saharan Africa 758
(453, 1,242)
3.44
(2.07, 5.67)
1,394
(733, 2,481)
2.6
(1.35, 4.56)
83.905013 −24.4186 −1.05
(−1.2, −0.91)
East Asia 85,356
(54,658, 117,796)
10.24
(6.55, 14.07)
35,484
(6,624, 113,906)
1.63
(0.3, 5.25)
-58.42823 −84.08203 −6.39
(−7.22, −5.55)
Eastern Europe 1,226
(239, 4,967)
0.43
(0.08, 1.73)
229
(33, 1,064)
0.06
(0.01, 0.3)
-81.32137 −86.04651 −8.05
(−9.34, −6.74)
Eastern sub-Saharan Africa 2,087
(1,397, 2,851)
2.81
(1.88, 3.81)
3,666
(2,364, 4,993)
2.31
(1.49, 3.15)
75.65884 −17.79359 −0.84
(−0.93, −0.76)
High-income Asia Pacific 76
(5, 490)
0.04
(0, 0.24)
6
(0, 42)
0
(0, 0.01)
−92.10526 −100 −10.57
(−11.45, −9.67)
High-income North America 12
(0, 91)
0
(0, 0.03)
2
(0, 13)
0
(0, 0)
−83.33333 0 −8.1
(−8.48, −7.72)
North Africa and Middle East 1,692
(666, 4,151)
1.02
(0.4, 2.51)
901
(487, 1,635)
0.2
(0.11, 0.36)
−46.74941 −80.39216 −5.61
(−5.69, −5.54)
Oceania 141
(79, 242)
5.13
(2.91, 8.66)
340
(187, 555)
4.86
(2.66, 7.99)
141.13475 −5.263158 −0.12
(−0.16, −0.09)
South Asia 11,367
(7,300, 15,453)
1.99
(1.28, 2.7)
17,018
(9,001, 27,968)
1.15
(0.61, 1.9)
49.714085 −42.21106 −1.98
(−2.12, −1.84)
Southeast Asia 13,851
(8,566, 19,405)
5.59
(3.46, 7.84)
11,937
(4,413, 23,356)
1.84
(0.68, 3.62)
−13.8185 −67.08408 −3.77
(−4.14, −3.39)
Southern Latin America 462
(104, 1,314)
1
(0.22, 2.83)
46
(0, 415)
0.05
(0, 0.47)
−90.04329 −95 −9.57
(-−9.85, −9.3)
Southern sub-Saharan Africa 483
(238, 891)
1.8
(0.89, 3.3)
580
(280, 1,066)
0.99
(0.48, 1.83)
20.082816 −45 −2.4
(−2.97, −1.82)
Tropical Latin America 1,181
(503, 2,302)
1.34
(0.58, 2.59)
543
(102, 1,626)
0.21
(0.04, 0.64)
−54.02202 −84.32836 −6.07
(−6.43, −5.71)
Western Europe 127
(1, 1,035)
0.02
(0, 0.18)
10
(0, 89)
0
(0, 0.01)
−92.12598 −100 −9.41
(−9.93, −8.88)
Western sub-Saharan Africa 939
(595, 1,338)
1.11
(0.71, 1.58)
2,149
(1,274, 3,265)
1.18
(0.7, 1.76)
128.86049 6.3063063 0.22
(0.1, 0.34)

ASRs, age-standardized rates; CI, confidence interval; EAPC, estimated annual percentage change; GBD, Global Burden of Disease; HAP, household air pollution; SDI, socio-demographic index; TBL, tracheal, bronchial, and lung.

Table 2

Numbers and ASRs per 100,000 cases of DALYs of TBL cancer attributable to HAP in 1990 and 2021, along with the relative changes and EAPC in ASRs per 100,000 cases from 1990 to 2021, categorized by global, SDI, and GBD regions

Location 1990 2021 1990 to 2021 EPAC (95% CI)
Number (95% CI) Age-standardized rate (95% UI) Number (95% CI) Age-standardized rate (95% UI) Relative change of numbers (%) Relative change of age-standardized rate (%)
Global 3,517,878
(2,180,370, 4,954,257)
84.14
(52.12, 118.39)
1,966,426
(758,405, 4,632,387)
22.41
(8.64, 52.88)
−44.10192736 −73.36581887 −4.77
(−5.3, −4.23)
High SDI 96,655
(29,164, 241,156)
9.09
(2.74, 22.67)
4,745
(8, 44,239)
0.24
(0, 2.25)
−95.09078682 −97.35973597 −12.53
(−13.04, −12.02)
High-middle SDI 1,163,420
(689,996, 1,803,939)
112.57
(66.87, 174.27)
206,102 (11,337, 1,050,686) 10.33
(0.57, 52.7)
−82.28481546 −90.82348761 −8.49
(−9.55, −7.41)
Middle SDI 1,671,512
(1,086,315, 2,285,378)
149.32
(97.36, 203.98)
839,797 (180,687, 2,339,653) 30.21
(6.47, 84.17)
−49.75824284 −79.76828288 −5.42
(−6.11, −4.73)
Low-middle SDI 445,032
(284,840, 612,946)
67.11
(42.95, 92.43)
652,895
(372,150, 1,019,360)
42.62
(24.26, 66.79)
46.70742778 −36.49232603 −1.59
(−1.77, −1.41)
Low SDI 138,525
(90,380, 198,515)
56.27
(36.73, 80.38)
261,337
(163,721, 369,459)
47.37
(29.74, 66.56)
88.65692113 −15.81659854 −0.73
(−0.79, −0.66)
Andean Latin America 10,904
(4,609, 19,183)
51.23
(21.65, 90.56)
4,909
(826, 15,808)
8.19
(1.38, 26.41)
−54.97982392 −84.01327347 −5.98
(−6.42, −5.54)
Australasia 118 (0, 1,139) 0.5 (0, 4.87) 12 (0, 78) 0.02 (0, 0.15) −89.83050847 −96 −10.21
(−10.83, −9.59)
Caribbean 10,041
(5,197, 16,622)
37.93
(19.6, 62.91)
9,837
(4,975, 17,716)
18.35
(9.28, 32.97)
−2.031670152 −51.62140786 −2.47
(−2.57, −2.36)
Central Asia 34,051
(12,422, 75,501)
67.83
(24.81, 150)
9,852
(3,390, 24,984)
10.97
(3.75, 28)
−71.06692902 −83.8272151 −6.86
(−7.68, −6.04)
Central Europe 94,441
(24,976, 257,761)
61.81
(16.35, 168.76)
15,362
(451, 107,021)
7.4
(0.22, 51.4)
−83.7337597 −88.02782721 −7.96
(−8.61, −7.3)
Central Latin America 20,065
(7,585, 42,805)
23.08
(8.68, 49.41)
15,357
(5,337, 37,584)
6.02
(2.09, 14.79)
−23.46374284 −73.91681109 −4.57
(−4.79, −4.36)
Central sub-Saharan Africa 21,971
(13,074, 35,776)
88.14
(52.92, 144.06)
41,009
(21,456, 73,267)
66.04
(34.64, 117.76)
86.65058486 −25.07374631 −1.08
(-1.22, -0.94)
East Asia 2,375,963
(1,518,184, 3,284,865)
255.17
(163, 352.18)
843,081
(161,745, 2,668,935)
37.66
(7.23, 119.33)
−64.51624036 −85.24121174 -6.64
(-7.45, -5.83)
Eastern Europe 34,747
(6,817, 140,915)
12.14
(2.4, 49.06)
5,989
(855, 27,104)
1.75
(0.25, 7.85)
−82.76397962 −85.58484349 −8.26
(−9.55, −6.94)
Eastern sub-Saharan Africa 59,757
(40,079, 82,103)
73.11
(48.95, 100)
102,700
(65,805, 140,126)
56.02
(36.07, 76.2)
71.86271064 −23.37573519 −1.1
(−1.19, −1.01)
High-income Asia Pacific 1,957 (135, 12,191) 0.95 (0.07, 5.92) 99 (0, 692) 0.02 (0, 0.15) −94.94123659 −97.89473684 −11.27
(−12.19, −10.34)
High-income North America 298
(0, 2,199)
0.09
(0, 0.66)
49
(0, 272)
0.01
(0, 0.04)
−83.55704698 −88.88888889 −8.51
(−8.88, −8.14)
North Africa and Middle East 47,716
(18,750, 117,627)
26.1
(10.27, 64.29)
26,142
(14,295, 46,604)
5.11
(2.79, 9.22)
−45.21334563 −80.42145594 −5.64
(−5.73, −5.54)
Oceania 4,008
(2,249, 6,966)
125.89
(70.82, 216.19)
9,632
(5,330, 15,696)
118.27
(64.97, 193.03)
140.3193613 −6.052903328 −0.15
(−0.2, −0.1)
South Asia 332,336
(213,197, 451,186)
52.19
(33.52, 70.93)
474,158
(251,259, 776,461)
29.98
(15.86, 49.2)
42.67428145 −42.55604522 −1.98
(−2.12, −1.84)
Southeast Asia 384,077
(236,845, 540,598)
140.76
(86.9, 197.31)
319,631
(118,242, 621,900)
45.5
(16.81, 88.79)
−16.77944787 −67.67547599 −3.81
(−4.17, −3.45)
Southern Latin America 12,080
(2,707, 34,680)
25.72
(5.76, 73.85)
1,038
(4, 9,338)
1.21
(0, 10.87)
−91.40728477 −95.29548989 −9.9
(−10.19, −9.61)
Southern sub-Saharan Africa 13,768
(6,638, 25,683)
47.31
(22.99, 87.73)
16,627
(8,037, 30,579)
26.22
(12.69, 48.24)
20.76554329 −44.57831325 −2.38
(−2.96, −1.81)
Tropical Latin America 31,544
(13,223, 62,227)
32.8
(13.85, 64.36)
13,219
(2,462, 39,517)
5.05
(0.94, 15.1)
−58.09345676 −84.60365854 −6.2
(−6.57, −5.83)
Western Europe 3,078
(26, 24,977)
0.57
(0, 4.6)
220
(0, 1,856)
0.03
(0, 0.22)
−92.85250162 −94.73684211 −9.64
(−10.14, −9.13)
Western sub-Saharan Africa 24,960
(15,770, 35,554)
27.09
(17.13, 38.6)
57,501
(33,950, 88,798)
27.68
(16.39, 42.15)
130.3725962 2.177925434 0.09
(−0.02, 0.2)

ASRs, age-standardized rates; CI, confidence interval; DALYs, disability-adjusted life years; EAPC, estimated annual percentage change; GBD, Global Burden of Disease; HAP, household air pollution; SDI, socio-demographic index; TBL, tracheal, bronchial, and lung.

The high SDI regions reported the lowest ASMR and ASDR in 2021, with 0.01 (95% UI: 0 to 0.1) and 0.24 (95% UI: 0 to 2.25), respectively. These regions also experienced the most significant reductions in ASMR and ASDR from 1990 to 2021, with decreases of 94.12% and 97.36%, respectively.

In 21 GBD regions, the Oceania region saw a sharp increase in TBL cancer-related deaths and DALYs, with mortality rising by 141.1%, from 141 to 340, and the number of DALYs increasing by 140.3%, from 4,008 to 9,632. Western sub-Saharan Africa experienced a slight increase in both ASMR, from 1.11 to 1.18 per 100,000, and ASDR, from 27.09 to 27.68 per 100,000.

The burden of TBL cancer due to HAP in 204 countries and territories

The ASMR of EAPC from 1990 to 2021 revealed significant regional variations. The Northern Mariana Islands experienced the highest annual increase in ASMR, with an EAPC of 7.1 (95% CI: 5.57 to 8.65), while the Republic of the Niger showed the smallest increase, with an EAPC of 0.28 (95% CI: 0.05 to 0.5). The Republic of Equatorial Guinea saw the largest decrease in ASMR, with an EAPC of −17.38 (95% CI: −19.3 to −15.42) (Figure 1A). For China, the EAPC of ASMR decreased by −4.2 (95% CI: −5.03 to −3.37).

Figure 1 ASMR (A) and ASDR (B) of the burden of TBL cancer due to HAP in 204 countries and territories, 1990–2021. The study is based on GBD 2021, in which Taiwan is treated as a province of China and its data are separately calculated and statistically analyzed. Therefore, there are differences in the raw data for the disease burden between Taiwan and mainland China, which is reflected in the color variations on the map. ASDR, age-standardized disability rate; ASMR, age-standardized mortality rate; EAPC, estimated annual percentage change; HAP, household air pollution; TBL, tracheal, bronchial, and lung.

Similarly, the EAPC for ASDR showed a significant increase in the Northern Mariana Islands (6.64, 95% CI: 5.14 to 8.17), while the People’s Republic of Vanuatu had the smallest increase (0.3, 95% CI: 0.19 to 0.42). The Republic of Equatorial Guinea experienced the largest decrease in ASDR, with an EAPC of −17.47 (95% CI: −19.41 to −15.49) (Figure 1B). In China, the EAPC for ASDR also showed a decreasing trend, with an EAPC of −4.73 (95% CI: −5.56 to −3.9).

The burden of TBL cancer due to HAP by SDI

A significant negative correlation was observed between both ASMR and ASDR and the SDI. Regions with an SDI below 0.4 showed relatively stable ASMR and ASDR, whereas regions with SDI values above 0.4 exhibited more pronounced declines in both indicators. Notably, East Asia, Southeast Asia, and Oceania regions, though showing similar trends, consistently reported higher levels of ASMR and ASDR at comparable SDI levels compared to other regions (Figure 2A). Additionally, the burden of TBL cancer attributable to HAP was notably lower in countries with higher SDI, where both ASMR and ASDR remained low and stable. However, despite China’s SDI approximating 0.7, its ASMR and ASDR for TBL cancer attributable to HAP remained higher than those of other countries with similar SDI levels, suggesting that, despite socio-economic development, China continues to face a disproportionately high burden (Figure 2B).

Figure 2 The burden of TBL cancer due to HAP in the 21 GBD regions and 204 countries and territories by SDI, 1990–2021. (A) ASMR (B) ASDR. ASDR, age-standardized disability rate; ASMR, age-standardized mortality rate; GBD, Global Burden of Disease; HAP, household air pollution; SDI, socio-demographic index; TBL, tracheal, bronchial, and lung.

The burden of TBL cancer due to HAP by age and sex

Globally, the burden of TBL cancer due to HAP varied significantly by sex, with ASMR and ASDR consistently higher in males than in females. In the 85–89 years age group, the ASMR for males was 10.51 per 100,000 (95% UI: 3.10 to 30.97), which was twice the rate in females (5.03 per 100,000; 95% UI: 1.42 to 13.84) (Figure 3A). In the 70–74 years age group, the ASDR was substantially higher in males (166.85 per 100,000; 95% UI: 56.79 to 427.96) than in females (81.44 per 100,000; 95% UI: 28.77 to 205.78) (Figure 3B). The ASMR and ASDR for TBL cancer due to HAP generally increased with age before declining. In males, ASMR peaked in the 85–89 years age group, while in females, it peaked in the 80–84 years age group at 5.12 per 100,000 (95% UI: 1.61 to 13.53). The lowest ASMR and ASDR were observed in the 25–29 years age group, with males showing an ASMR of 0.04 per 100,000 (95% UI: 0.02 to 0.08) and an ASDR of 2.41 per 100,000 (95% UI: 1.12 to 4.82), and females having an ASMR of 0.03 per 100,000 (95% UI: 0.02 to 0.06) and an ASDR of 2.00 per 100,000 (95% UI: 0.99 to 3.67) (Figure 3).

Figure 3 Trends in ASMR (A) and ASDR (B) for TBL cancer due to HAP by age and sex in 2021. ASDR, age-standardized disability rate; ASMR, age-standardized mortality rate; DALYs, disability-adjusted life years; HAP, household air pollution; TBL, tracheal, bronchial, and lung.

APC model

The APC model revealed distinct patterns in age effects on mortality and DALYs rates for TBL cancer due to HAP. Globally, mortality peaked at around 70 years of age, while DALYs rates peaked at approximately 60 years (Figure 4A, Figure S1A). This pattern was consistent across all SDI regions, although with regional variations. In low and middle SDI regions, mortality rates peaked at approximately 90 years, while DALYs rates peaked around 70 years. In high SDI regions, mortality peaked between 40 and 55 years, while DALYs peaked at 50 years. Period effects showed a general decline in RRs globally, except in low SDI regions, where a decrease was observed until 2015, followed by a slight increase thereafter (Figure 4B, Figure S1B). Cohort effects revealed a general decline in TBL cancer mortality and DALYs rates in high SDI regions, while low SDI regions experienced an upward trend for individuals born between 1900 and 1920, followed by a decline for later birth cohorts (Figure 4C, Figure S1C).

Figure 4 Age, period, and birth cohort effects of the APC model on mortality in the TBL cancer due to HAP by SDI. (A) Age effects. (B) Period effects. (C) Birth cohort effects. APC, age-period-cohort; HAP, household air pollution; RR, risk ratio; SDI, socio-demographic index; TBL, tracheal, bronchial, and lung.

Future projections of the BAPC model

The BAPC model projections indicated that, despite significant declines in ASMR and ASDR from 1990 to 2021, both are expected to increase over the next 30 years (Figure 5). A gender disparity is projected, with females experiencing a greater burden of TBL cancer due to HAP than males. By 2051, the projected ASMR and ASDR for males are 3.97 (95% UI: −3.35 to 11.30) and 92.32 (95% UI: −75.67 to 260.31), respectively. For females, the projections are 4.08 (95% UI: −3.76 to 11.92) and 116.66 (95% UI: −112.38 to 345.69), respectively.

Figure 5 Global trends in the burden of TBL cancer due to HAP by sex over the next 30 years as predicted by the BAPC model. (A) ASMR. (B) ASDR. ASDR, age-standardized disability rate; ASMR, age-standardized mortality rate; BAPC, Bayesian age-period-cohort; HAP, household air pollution; TBL, tracheal, bronchial, and lung.

Discussion

Our analysis utilizing the GBD 2021 data demonstrated a notable decrease in the global burden of TBL cancer attributable to HAP between 1990 and 2021. A study by Lu et al. (19) accentuated the escalating cancer burden linked to ambient particulate matter pollution (APMP) and the declining trend in HAP-related cancer burden on a global scale. This research underscored the imperative for targeted interventions to mitigate HAP in low SDI regions, where the burden persists at elevated levels.

Specifically, there was a 39.2% reduction in global deaths and a 44.1% decrease in DALYs cases due to HAP-related TBL cancer. This trend mirrors the heightened global awareness of HAP’s health risks and the efficacy of interventions designed to curb indoor air pollution exposure. Nonetheless, this general decline masks considerable regional disparities. High SDI regions, including North America and Western Europe, witnessed substantial declines in ASMR and ASDR, whereas low SDI regions, especially sub-Saharan Africa and Oceania, experienced rises in mortality and DALYs rates. These increases highlight the persistent challenges faced by these regions in adopting cleaner energy sources, compounded by economic limitations, restricted access to clean cooking technologies, and a dependence on traditional biomass fuels for cooking, which sustain high levels of HAP exposure (20,21).

Importantly, indoor air quality in these regions is shaped not only by fuel combustion but also by additional sources of pollution, such as tobacco smoking, which further compound overall exposure and health risks (22). Moreover, inadequate household ventilation exacerbates the accumulation of multiple pollutants, intensifying harmful emissions and worsening respiratory health outcomes (23). Effective interventions must therefore adopt an integrated approach that simultaneously addresses cleaner fuel adoption, ventilation improvements, and broader household air quality management.

Our findings align with previous studies showing significant reductions in HAP-related TBL cancer in high SDI regions, where clean cooking technologies and policies to reduce PM have been effective (24). In contrast, low SDI regions continued to face considerable barriers in reducing HAP exposure, particularly in sub-Saharan Africa, where solid fuels for cooking remained widely used, contributing to higher levels of indoor air pollution and, consequently, higher cancer burden (25). A particularly striking finding in our study was the situation in China. Despite rapid socio-economic development and a high SDI, the burden of TBL cancer attributable to HAP remained disproportionately high compared to other countries with similar SDI levels (26-28). This aligned with previous research suggesting that China’s reliance on coal for domestic heating and cooking, particularly in rural areas, had continued to contribute significantly to HAP-related cancer risks. This observation highlighted the need for targeted, region-specific interventions to address the persistent challenge of HAP in even the most economically advanced nations.

Our study also highlighted significant sex disparities in the burden of TBL cancer attributable to HAP. Males consistently experienced higher ASMR and ASDR compared to females, with the highest rates occurring in the 85–89 years age group for men and the 80–84 years age group for women. This finding was in line with existing literature, which attributed higher exposure levels among men to a combination of HAP and other risk factors such as smoking and occupational exposures to carcinogens (29). However, our data also revealed a higher EAPC for women, suggesting a more pronounced upward trend in cancer burden among females, particularly in older age groups. This increasing burden among women could be attributed to several factors. First, women were often more likely to be responsible for cooking in households relying on solid fuels, leading to prolonged exposure to HAP (30,31). Second, sociocultural factors in many LMICs, such as gendered roles in household chores, exacerbated women’s exposure to indoor air pollution (32,33). Additionally, the increased life expectancy among women in many regions may have contributed to the higher burden observed in older age groups.

Our findings suggested that the burden of TBL cancer attributable to HAP increased with age, with mortality rates peaking in the 85–89 years age group for men and 80–84 years for women. This pattern reflected the natural progression of cancer risk as individuals aged, compounded by cumulative exposure to risk factors such as HAP. The APC reveals higher mortality peaks in low SDI regions compared to high SDI regions, suggesting that chronic HAP exposure in low SDI regions significantly exacerbates the cancer burden among older adults. This suggested the need for targeted interventions to reduce HAP exposure in low SDI regions, including transitioning to cleaner fuels and improving indoor air quality (34). Enhancing healthcare services and promoting cancer prevention programs are also essential to mitigate the high cancer burden in these vulnerable populations.

Projections from the BAPC model indicated that the global burden of TBL cancer due to HAP was expected to rise over the next 30 years, particularly among women. These projections underscored the need for continued efforts to mitigate HAP exposure, especially in low SDI regions, where the burden remained high and was projected to increase in the coming decades. Given the sex disparities in burden and projections, public health policies needed to focus on reducing HAP exposure in a gender-sensitive manner.

The research outcomes highlight a pressing demand for tailored public health strategies in low SDI regions, which continue to grapple with substantial TBL cancer rates associated with HAP. There is an urgent call for localized approaches to manage these issues, including the promotion of eco-friendly cooking solutions, broader access to cleaner energy sources, and public awareness initiatives to inform about HAP risks (35,36). Although high-income nations have successfully reduced HAP exposure, ongoing vigilance and action are required, particularly for rural and indigenous communities that are often overlooked. The study also signals the importance of gender-responsive policies in light of the rising TBL cancer rates due to HAP, predominantly affecting women. It is vital for ongoing research to explore the interplay among HAP, socio-economic conditions, gender, and age to craft more precise interventions for high-risk groups (37). There is a clear need for both national and global health entities to focus on interventions that tackle the socio-economic and cultural aspects contributing to HAP-linked cancer risks, especially in areas where the disease burden is growing.

This study provided important insights into the global burden of TBL cancer attributable to HAP, but limitations include reliance on model-based estimates, particularly in regions with limited data (38). Future research should focus on improving exposure assessments and integrating more granular data on fuel use and health disparities. Longitudinal studies on the long-term effects of HAP on cancer incidence and mortality are also needed to enhance understanding of disease progression.

In addition to policy-level strategies, household-level actions play a pivotal role in reducing exposure to HAP and its associated cancer burden. A central priority is the gradual transition from solid fuels, such as wood, coal, and agricultural residues, to cleaner energy sources, including liquefied petroleum gas, electricity, and biogas, which markedly lower the emission of carcinogenic pollutants. Equally important is improving ventilation in cooking and living spaces; the installation of chimneys, the use of exhaust fans, or relocating cooking to semi-open areas can substantially diminish the accumulation of harmful indoor pollutants. Reducing secondary sources of exposure, particularly indoor tobacco use and the burning of incense or candles, further lessens overall pollutant load. Regular inspection and maintenance of cooking stoves and ventilation devices are also essential to ensure optimal performance and pollutant removal. Finally, raising awareness of the health risks of HAP, particularly among women and children, who often bear the highest exposure burden, can foster sustainable behavioral changes and community-level adoption of safer practices. Collectively, these practical recommendations, while relatively low-cost and feasible, can complement structural clean-energy policies and accelerate progress toward reducing the global burden of tracheal, bronchial, and lung cancers attributable to HAP.


Conclusions

In conclusion, our study highlights the significant reductions in the global burden of TBL cancer attributable to HAP in high SDI regions, while emphasizing the continued challenges faced by low SDI regions. The projected rise in TBL cancer burden, particularly among women, underscores the urgent need for targeted, gender-sensitive interventions to reduce HAP and its associated health impacts. Addressing these disparities through region-specific policies and interventions is essential to achieving global reductions in the burden of HAP-related TBL cancer.


Acknowledgments

None.


Footnote

Reporting Checklist: The authors have completed the STROBE reporting checklist. Available at https://jtd.amegroups.com/article/view/10.21037/jtd-2025-1309/rc

Peer Review File: Available at https://jtd.amegroups.com/article/view/10.21037/jtd-2025-1309/prf

Funding: This study was supported by Natural Science Foundation of Guangdong Province, China (grant No. 2024A1515013161).

Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://jtd.amegroups.com/article/view/10.21037/jtd-2025-1309/coif). H.L. serves as an unpaid editorial board member of Journal of Thoracic Disease. The other authors have no conflicts of interest to declare.

Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments.

Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.


References

  1. Household air pollution and health. Geneva: World Health Organization; 2022. Available online: https://www.who.int/news-room/fact-sheets/detail/household-air-pollution-and-health
  2. Puthumana JS, Ngaage LM, Borrelli MR, et al. Risk factors for cooking-related burn injuries in children, WHO Global Burn Registry. Bull World Health Organ 2021;99:439-45. [Crossref] [PubMed]
  3. Dutta A, Olopade CO. Transitioning to gaseous and liquid fuels: a right step towards clean cooking in low-income and middle-income countries. Lancet Respir Med 2024;12:257-8. [Crossref] [PubMed]
  4. Zhou RX, Liao HJ, Hu JJ, et al. Global Burden of Lung Cancer Attributable to Household Fine Particulate Matter Pollution in 204 Countries and Territories, 1990 to 2019. J Thorac Oncol 2024;19:883-97. [Crossref] [PubMed]
  5. Yang X, Zhang T, Zhang X, et al. Global burden of lung cancer attributable to ambient fine particulate matter pollution in 204 countries and territories, 1990-2019. Environ Res 2022;204:112023. [Crossref] [PubMed]
  6. World Health Organization. Burning opportunity: clean household energy for health, sustainable development, and wellbeing of women and children. Geneva: World Health Organization;2016.
  7. Migała-Warchoł A, Ziółkowski B, Lew A, et al. Factors Enabling Access to Affordable, Reliable, Sustainable and Modern Energy in the European Union. Energies 2025;18:722.
  8. Qin K, Qing J, Wang Q, et al. Epidemiological shifts in chronic kidney disease: a 30-year global and regional assessment. BMC Public Health 2024;24:3519. [Crossref] [PubMed]
  9. Kuang Z, Wang J, Liu K, et al. Global, regional, and national burden of tracheal, bronchus, and lung cancer and its risk factors from 1990 to 2021: findings from the global burden of disease study 2021. EClinicalMedicine 2024;75:102804. [Crossref] [PubMed]
  10. Global burden of 369 diseases and injuries in 204 countries and territories, 1990-2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet 2020;396:1204-22. [Crossref] [PubMed]
  11. Global, regional, and national burden of diabetes from 1990 to 2021, with projections of prevalence to 2050: a systematic analysis for the Global Burden of Disease Study 2021. Lancet 2023;402:203-34. [Crossref] [PubMed]
  12. Sun P, Yu C, Yin L, et al. Global, regional, and national burden of female cancers in women of child-bearing age, 1990-2021: analysis of data from the global burden of disease study 2021. EClinicalMedicine 2024;74:102713. [Crossref] [PubMed]
  13. Global incidence, prevalence, years lived with disability (YLDs), disability-adjusted life-years (DALYs), and healthy life expectancy (HALE) for 371 diseases and injuries in 204 countries and territories and 811 subnational locations, 1990-2021: a systematic analysis for the Global Burden of Disease Study 2021. Lancet 2024;403:2133-61. [Crossref] [PubMed]
  14. Bai Z, Wang H, Shen C, et al. The global, regional, and national patterns of change in the burden of nonmalignant upper gastrointestinal diseases from 1990 to 2019 and the forecast for the next decade. Int J Surg 2025;111:80-92. [Crossref] [PubMed]
  15. Yu W, Gargett T, Du Z. A Poisson distribution-based general model of cancer rates and a cancer risk-dependent theory of aging. Aging (Albany NY) 2023;15:8537-51. [Crossref] [PubMed]
  16. Bai Z, Han J, An J, et al. The global, regional, and national patterns of change in the burden of congenital birth defects, 1990-2021: an analysis of the global burden of disease study 2021 and forecast to 2040. EClinicalMedicine 2024;77:102873. [Crossref] [PubMed]
  17. Global, regional, and national burden of stroke and its risk factors, 1990-2021: a systematic analysis for the Global Burden of Disease Study 2021. Lancet Neurol 2024;23:973-1003. [Crossref] [PubMed]
  18. Pu L, Wang L, Zhang R, et al. Projected Global Trends in Ischemic Stroke Incidence, Deaths and Disability-Adjusted Life Years From 2020 to 2030. Stroke 2023;54:1330-9. [Crossref] [PubMed]
  19. Lu J, Zhao X, Gan S. Global, regional and national burden of tracheal, bronchus, and lung cancer attributable to ambient particulate matter pollution from 1990 to 2021: an analysis of the global burden of disease study. BMC Public Health 2025;25:108. [Crossref] [PubMed]
  20. Granier C, Liousse C, McDonald B, et al. Anthropogenic Emissions Inventories of Air Pollutants. In: Akimoto H, Tanimoto H. editors. Handbook of Air Quality and Climate Change. Singapore: Springer Nature Singapore; 2023:3-52.
  21. Josyula S, Lin J, Xue X, et al. Household air pollution and cancers other than lung: a meta-analysis. Environ Health 2015;14:24. [Crossref] [PubMed]
  22. Ni Y, Shi G, Qu J. Indoor PM(2.5), tobacco smoking and chronic lung diseases: A narrative review. Environ Res 2020;181:108910. [Crossref] [PubMed]
  23. Raju S, Siddharthan T, McCormack MC. Indoor Air Pollution and Respiratory Health. Clin Chest Med 2020;41:825-43. [Crossref] [PubMed]
  24. Lai PS, Lam NL, Gallery B, et al. Household Air Pollution Interventions to Improve Health in Low- and Middle-Income Countries: An Official American Thoracic Society Research Statement. Am J Respir Crit Care Med 2024;209:909-27. [Crossref] [PubMed]
  25. Jafta N, Shezi B, Buthelezi M, et al. Household air pollution and respiratory health in Africa: persistent risk and unchanged health burdens. Curr Opin Pulm Med 2025;31:89-97. [Crossref] [PubMed]
  26. Liu H, Qi J, Ren W, et al. Disparities in trends and drivers of the burden of tracheal, bronchus, and lung cancer among Chinese population during 1990-2021: a systematic analysis for the global burden of disease study 2021. Sci Bull (Beijing) 2025;70:496-9. [Crossref] [PubMed]
  27. Nie P, Sousa-Poza A, Xue J. Fuel for Life: Domestic Cooking Fuels and Women's Health in Rural China. Int J Environ Res Public Health 2016;13:810. [Crossref] [PubMed]
  28. Barone-Adesi F, Chapman RS, Silverman DT, et al. Risk of lung cancer associated with domestic use of coal in Xuanwei, China: retrospective cohort study. BMJ 2012;345:e5414. [Crossref] [PubMed]
  29. Abbasi-Kangevari M, Malekpour MR, Masinaei M, et al. Effect of air pollution on disease burden, mortality, and life expectancy in North Africa and the Middle East: a systematic analysis for the Global Burden of Disease Study 2019. Lancet Planet Health 2023;7:e358-69. [Crossref] [PubMed]
  30. Gall ET, Carter EM, Earnest CM, et al. Indoor air pollution in developing countries: research and implementation needs for improvements in global public health. Am J Public Health 2013;103:e67-72. [Crossref] [PubMed]
  31. Smith KR, Samet JM, Romieu I, et al. Indoor air pollution in developing countries and acute lower respiratory infections in children. Thorax 2000;55:518-32. [Crossref] [PubMed]
  32. Duflo E. Women Empowerment and Economic Development. J Econ Lit 2012;50:1051-79.
  33. Pachauri S, Rao ND. Gender impacts and determinants of energy poverty: are we asking the right questions? Curr Opin Environ Sustain 2013;5:205-15.
  34. Local Burden of Disease Household Air Pollution Collaborators. Mapping development and health effects of cooking with solid fuels in low-income and middle-income countries, 2000-18: a geospatial modelling study. Lancet Glob Health 2022;10:e1395-411. [Crossref] [PubMed]
  35. Smith KR, Mehta S, Maeusezahl-Feuz M. Indoor Air Pollution from Household Solid Fuels Use. In: Ezzati M, Lopez AD, Rodgers A et al. editors. Comparative Quantification of Health Risks: Global and Regional Burden of Disease Attributable to Selected Major Risk Factors. Geneva: World Health Organization; 2004:1435-93.
  36. Elehinafe FB, Ibukun FD, Lasebikan OL, et al. A review on indoor air pollution from lightning sources in developing countries—The pollutants, the prevention and the latest technologies for control. Applied Chemical Engineering 2023;6.
  37. Véliz KD, Alcantara-Zapata DE, Chomalí L, et al. Gender-differentiated impact of PM2.5 exposure on respiratory and cardiovascular mortality: a review. Air Qual Atmos Health 2024;17:1565-86.
  38. Wang F, Ma B, Ma Q, et al. Global, regional, and national burden of inguinal, femoral, and abdominal hernias: a systematic analysis of prevalence, incidence, deaths, and DALYs with projections to 2030. Int J Surg 2024;110:1951-67. [Crossref] [PubMed]
Cite this article as: Lao S, Li J, Li W, Liang H, Liang W, Li F, Wang W. Global trends in the burden of tracheal, bronchus, and lung cancer attributable to household air pollution. J Thorac Dis 2025;17(11):9943-9957. doi: 10.21037/jtd-2025-1309

Download Citation