Skip to main content

Short-term effects and economic burden of air pollutants on acute lower respiratory tract infections in children in Southwest China: a time-series study

Abstract

Background

There are few studies on the effects of air pollutants on acute lower respiratory tract infections (ALRI) in children. Here, we investigated the relationship of fine particulate matter (PM2.5), inhalable particulate matter (PM10), sulfur dioxide (SO2), and nitrogen dioxide (NO2) with the daily number of hospitalizations for ALRI in children in Sichuan Province, China, and to estimate the economic burden of disease due to exposure to air pollutants.

Methods

We collected records of 192,079 cases of childhood ALRI hospitalization between January 1, 2017 and December 31, 2018 from nine municipal/prefecture medical institutions as well as the simultaneous meteorological and air pollution data from 183 monitoring sites in Sichuan Province. A time series-generalized additive model was used to analyze exposure responses and lagged effects while assessing the economic burden caused by air pollutant exposure after controlling for long-term trends, seasonality, day of the week, and meteorological factors.

Results

Our single-pollutant model shows that for each 10 μg/m3 increase in air pollutant concentration (1 μg/m3 for SO2), the effect estimates of PM2.5, PM10, SO2, and NO2 for pneumonia reached their maximum at lag4, lag010, lag010, and lag07, respectively, with relative risk (RR) values of 1.0064 (95% CI, 1.0004–1.0124), 1.0168(95% CI 1.0089–1.0248), 1.0278 (95% CI 1.0157–1.0400), and 1.0378 (95% CI, 1.0072–1.0692). By contrast, the effect estimates of PM2.5, PM10, SO2, and NO2 for bronchitis all reached their maximum at lag010, with RRs of 1.0133 (95% CI 1.0025–1.0242), 1.0161(95% CI 1.0085–1.0238), 1.0135 (95% CI 1.0025–1.0247), and 1.1133(95% CI 1.0739–1.1541). In addition, children aged 5–14 years were more vulnerable to air pollutants than those aged 0–4 years (p < 0.05). According to the World Health Organization’s air quality guidelines, the number of ALRI hospitalizations attributed to PM2.5, PM10, and NO2 pollution during the study period was 7551, 10,151, and 7575, respectively, while the incurring economic burden was CNY 2847.06, 3827.27, and 2855.91 million.

Conclusion

This study shows that in Sichuan Province, elevated daily average concentrations of four air pollutants lead to increases in numbers of childhood ALRI hospitalizations and cause a serious economic burden.

Peer Review reports

Background

Air pollution is a major environmental hazard as well as a leading cause of morbidity and mortality worldwide. The number of deaths caused by air pollution has substantially increased globally over the past 20 years, with ambient air pollution causing 2.9 million deaths in 2000, 4.2 million deaths in 2015, and 4.5 million deaths in 2019 [1]. Children are an especially vulnerable group, and a large body of research has shown that air pollution is associated with poor birth outcomes and respiratory diseases in children and poor mental health in adulthood [2,3,4,5,6]. Acute lower respiratory tract infections (ALRI), including lung and alveolar (pneumonia) as well as airway infections (bronchitis and bronchiolitis), are the leading cause of death in children under the age of 5 [7]. In the Global Burden of Diseases, Injuries, and Risk Factors Study 2019, lower respiratory tract infections ranked second among the top 10 causes of disability-adjusted life years in children under the age of 10 [8]. The World Health Organization (WHO) reports that approximately 93% of children under the age of 15 worldwide suffer from air pollution on a daily basis and that in 2016, 600,000 children died from ALRI caused by air pollution [9].

Children are particularly vulnerable to air pollution because their lung immune system is still developing. A number of studies conducted in different countries have demonstrated an association between long-term or short-term exposure to air pollutants and childhood ALRI hospitalization [10,11,12,13]. There are also several studies from China that have focused on the effects of air pollutants on the hospitalization of children with ALRI. For examples, a study conducted in 25 Chinese cities found that each 10 μg/m3 increase in the concentrations of the gaseous pollutants SO2 and NO2 resulted in an increase of 0.54% (95% CI, 0.30–0.79) and 0.60% (95% CI, 0.22–0.99) in respiratory hospitalizations among children aged 0–14 years, respectively [14]. Yang et al. investigated the association between short-term exposure to air pollutants and hospitalizations for ALRI in children aged 0–14 years in four Chinese cities, Guangzhou, Shanghai, Wuhan, and Xining [15]. Another study conducted in Shanghai found that with each interquartile range (IQR) increase in PM2.5, emergency department visits for pediatric bronchitis and pneumonia increased by 1.53% (95% CI, 0.01–3.08) and 1.90% (95% CI, 0.30–3.52) on lag3, respectively [16]. Moreover, a study conducted in Guangzhou by Liang et al. made a similar conclusion that an elevation in particulate matter concentrations is significantly associated with an increase in outpatient visits for pneumonia and bronchitis in children [17]. However, the above studies were conducted mostly in central and eastern parts of China, while there are relatively few studies on the effects of air pollution on childhood ALRI hospitalization in the western regions; moreover, most of the existing ones are single-city studies [18,19,20]. It has also been shown that children in the western regions of China are most vulnerable to air pollution [21]. Sichuan Province, which constitutes the largest economy in western China, is located in the southwest and has approximately 90 million inhabitants [22]. In terms of topography, it is roughly divided into the Sichuan Basin and the Western Sichuan Plateau. As the degradation of pollutants is subject to geographical factors, the alpine terrain of the basin and surrounding plateaus causes relatively slow diffusion of atmospheric pollutants, making the Sichuan Basin the fourth most polluted area, after the Beijing-Tianjin-Hebei region, the Yangtze River Delta, and the Pearl River Delta [23]. In a study conducted in Sichuan by Pu et al., only the effects of particulate matter were analyzed and the health effects of common air pollutants on lower respiratory tract diseases in children were not fully reflected [24]. Here, we therefore performed a time-series analysis to more comprehensively assess the effects of short-term exposure to PM2.5, PM10, SO2, and NO2 on childhood hospitalization for ALRI (including pneumonia and bronchitis) based on hospitalization data of children aged 0–14 years from urban medical institutions in nine cities and prefectures of Sichuan Province between 2017 and 2018. This study provides more scientific basis for a comprehensive evaluation of the impact of air pollution on children’s health across different regions. Moreover, we analyzed the effects of sex, age, and season, while estimating the economic burden of diseases due to exposure to air pollutants.

Methods

Data on ALRI in children

Sichuan Province comprises a total of 21 cities and prefectures, including 18 in the Sichuan Basin and three in the western Sichuan Plateau. For this study, we collected data on a total of 192,079 cases of pediatric ALRI inpatients aged 14 years and younger from January 1, 2017 to December 31, 2018, including each child’s age, sex, home address, admission date, discharge date, disease diagnosis, disease code, and total hospitalization cost, from nine city or prefecture hospitals in Chengdu, Guang’an, Luzhou, Liangshan Yi Autonomous Prefecture, Mianyang, Meishan, Nanchong, Yibin, and Zigong of Sichuan Province. The causes of hospital visits were coded according to the International Classification of Diseases, 10th Revision (ICD-10), as follows: pneumonia (J12–J18), bronchitis (J20–J21), and other acute lower respiratory tract infections (J22). The latter were not analyzed separately because we counted only nine children with such infections in the entire sample. Of the above nine cities and prefectures, eight are located in the Sichuan Basin and one in the Western Sichuan Plateau. Our samples therefore represent, to a certain extent, the overall situation in Sichuan Province.

Meteorological and pollutant data

We collected daily meteorological data from the Sichuan Meteorological Bureau (http://www.scdata.net.cn/) between January 1, 2017 and December 31, 2018, including daily average temperature and relative humidity. Particulate matter with an aerodynamic diameter ≤ 10 μm (PM10), particulate matter with an aerodynamic diameter ≤ 2.5 μm (PM2.5), sulfur dioxide (SO2), and nitrogen dioxide (NO2) were chosen as indicators of outdoor air pollution, as these four pollutants are closely associated with respiratory diseases. The daily average concentrations of PM2.5, PM10, SO2, and NO2 were collected from environmental monitoring stations in Sichuan Province (including a total of 183 air monitoring stations) during the same period.

Air pollutant exposure was assessed by using the inverse distance weighting (IDW) method. Specifically, the locations of all monitoring stations and home addresses of hospitalized cases of pediatric ALRI were geo-coded using the Gaudet Map API (https://lbs.amap.com/). For each hospitalized case and monitoring station, the inverse distance (1/distance2) weighted average of the concentrations from all monitoring stations was used to estimate air pollutant exposure on the hospitalization day (lag0), single-day lags from the current day (lag0) and each 1–10 days before the ALR events (lag1, lag2, lag3, lag4, lag5, lag6, lag7, lag8, lag9, and lag10), as well as multi-day moving average lag exposures (lag01, lag02, lag03, lag04, lag05, lag06, lag07, lag08, lag09, and lag010).

Health efects of air pollutants exposure in overall and subgroup

populationThe correlation between air pollutants and meteorological indicators was analyzed using the Spearman correlation test, with absolute values of the correlation coefficient r closer to 1 indicating a stronger correlation.

Time series analysis methods applied to the generalized additive model (GAM) have been repeatedly used to assess the association between air pollutants and hospitalizations for respiratory diseases. Since the daily number of hospitalizations for pediatric ALRI usually obeys an over-dispersed Poisson distribution, a quasi-Poisson GAM model was adopted in this study. As mentioned above, our analysis of the correlations of PM2.5, PM10, SO2, and NO2 with the daily number of hospitalizations was based on single-day lags exposures (lag1, lag2, lag3, lag4, lag5, lag6, lag7, lag8, lag9, and lag10) and multi-day moving average lags exposures (lag01, lag02, lag03, lag04, lag05, lag06, lag07, lag08, lag09, and lag010). Daily average temperature and relative humidity were introduced into the model as control variables. The GAM model was formulated as follows:

$${\displaystyle \begin{array}{l} Log\left[E(Yi)\right]=\alpha +\beta Zi+s\left( time,k= df+1\right)+\textrm{s}\ \left(\textrm{temperature},\textrm{k}=\textrm{df}+1\right)\ \\ {}\kern5.5em +\textrm{s}\ \left(\textrm{humidity},\textrm{k}=\textrm{df}+1\right)+\textrm{as}.\textrm{factor}\left(\textrm{dow}\right)\end{array}}$$
(1)

where E (Yi) is the expected number of children hospitalized for ALRI on day i; α stands for the model intercept; Zi represents the air pollutant concentration on day i (μg/m3); β represents the exposure - response coefcient, which is the increase in the number of daily hospitalizations caused by the per unit increase in pollutant concentration; s is a non-smooth parameter item; and df is the degree of freedom. The time stands for a date variable with a degree of freedom of 7/year; dow is an indicator variable of “day of the week”; temperature and humidity stand for the average daily temperature and relative humidity, whose degrees of freedom are both 3.

The exposure-response coefficient β was estimated by the generalized additive model (Eq. 1). Subsequently, the relative risk (RR) and 95% CI of children hospitalized with ALRI for every 10 μg/m3 (1μg/m3 for SO2) increase in the concentration of air pollutants were calculated.

To identify potentially susceptible populations, we performed stratified analyses according to sex (boys and girls), age (0–1 years, 2–4 years, and 5–14 years), and season (hot season, June to August; transitional season, April, May, September, and October; cold season, November to March). The 95% confidence interval (95% CI) for the difference in the effect estimates between different categories in each subgroup (e.g., boys and girls) was calculated by the following formula to test whether a difference was statistically significant:

$$\left(\hat{Q}1-\hat{Q}2\right)\pm \sqrt{{\left(S\hat{E}1\right)}^2+{\left(S\hat{E}2\right)}^2}$$
(2)

where \(\hat{Q}1\) and \(\hat{Q}2\) are the estimates of different categories in each subgroup, and \(S\hat{E}1\) and \(S\hat{E}2\) are the corresponding standard errors for each estimate [25].

Attributable health risks and economic costs due to air pollution

The attributable number of ALRI hospitalizations due to air pollutants exposure was calculated according to the attributable risk equation [26]. The air pollutant concentrations in the air quality guidelines issued by the WHO were used as a standard reference (24-hour average: 15 μg/m3 for PM2.5, 45 μg/m3 for PM10, 40 μg/m3 for SO2, and 25 μg/m3 for NO2) [27]. The attribution eqs. 3, 4, and 5 were formulated as follows:

$${\textrm{AN}}_{\textrm{i}}=\left(\exp \left(\upbeta \ast {\Delta \textrm{AP}}_{\textrm{i}}\right)\hbox{-} 1\right)/\exp \left(\upbeta \ast {\Delta \textrm{AP}}_{\textrm{i}}\right)\ast {\textrm{N}}_{\textrm{i}}\kern1em$$
(3)
$$meanC={c}_h+\textrm{dPCDI}\times {\textrm{meanT}}_{\textrm{h}}$$
(4)
$$\varDelta C= AN\times meanC$$
(5)

In these equations, ANi is the number of hospitalizations attributed to air pollutants exposure on day i, β is the exposure-response coefficient (Eq. 1) between pollutants and the number of ALRI hospitalizations, ΔAPi is the difference between the observed and reference concentrations of air pollutants on day i, Ni is the number of ALRI hospitalizations on day i, meanC is the average total economic cost of hospitalization per case, ch is the average total hospitalization cost per case during the study period, dPCDI is the per capita daily disposable income in Sichuan Province, meanTh is the average number of hospitalization days per case, and ΔC is the total economic cost attributed to air pollution.

All statistical analyses in this study were carried out using R3.4.3, and the quasi-Poisson regression model was constructed using the package “mgcv”. The test level was α = 0.05.

Sensitivity analysis

Three types of sensitivity analyses were performed to verify the stability of the model. First, if a significant association between a given air pollutant and childhood ALRI was observed, a two-pollutant models was further fitted to assess the robustness of our results. Second, the stability of the time trend was tested and the model was fitted by varying the degrees of freedom of the time series (df = 5, 6, 8, 9/year). Finally, a sensitivity analysis was performed by assessing the number of ALRI cases within a circular 50 km area around the air monitoring station. By doing this we were able to assess the potential impact of the distance between the air pollution monitoring site and the patient’s home address [28].

Results

Table 1 provides air pollution levels, meteorological variables, and daily ALRI hospitalizations. The daily average concentrations of PM2.5, PM10, SO2, and NO2 during the study period were 46.01, 71, 12.01, and 28.93 μg/m3, respectively. The daily average concentrations of PM2.5, PM10, and NO2 exceeded the standards of the air quality guidelines issued by the WHO (15, 45, and 25 μg/m3) [27], and the numbers of days on which the concentrations of PM2.5, PM10, and NO2 exceeded the standards were 698, 497, and 439, respectively. The daily average temperature and relative humidity were 17.81 °C and 76.82%. A total of 192,079 hospitalized cases of childhood ALRI were enrolled during the study period, including 129,870 cases of pneumonia and 62,200 of bronchitis. The average daily number of hospitalizations was higher in boys than in girls, and the age group from 0 to 1 years displayed the highest average daily hospitalizations (111 cases). Regarding the three seasons we defined for our analyses, the highest average daily number of hospitalizations (336 cases) was recorded during the cold season. Moreover, the daily number of hospitalizations for pneumonia was larger than that for bronchitis.

Table 1 Description of daily air pollutants, meteorological parameters, and hospitalized characteristics of children with ALRI in the Sichuan Province, China, during 2017–2018

We observed significant positive correlations between PM2.5 and PM10, SO2, and NO2, with the correlation coefficient r ranging from 0.3465 to 0.9442. In addition, PM10 showed a strong positive correlation with SO2 and NO2, with correlation coefficients of 0.3428 and 0.6392, respectively. The coefficient for the correlation between SO2 and NO2 was 0.2, and our analyses revealed a negative correlation between air pollutants and meteorological variables (p < 0.05) (Table S1).

Effects of pollutants on health of subgroups and overall population

Fig. 1 shows the relationship between each 10 μg/m3 increase in the concentrations of PM2.5, PM10, SO2, and NO2 and hospitalizations for overall ALRI, pneumonia or bronchitis, at different lag days. While the effects of PM2.5, PM10, SO2, and NO2 on hospitalizations all reached their maximum at lag010. The differences in these effects were statistically significant. For pneumonia, the effect estimates for PM2.5, PM10, SO2, and NO2 reached their maximum at lag4, lag010, lag010, and lag07, respectively, with RR values of 1.0064 (95% CI, 1.0004–1.0124), 1.0168(95% CI 1.0089–1.0248), 1.0278 (95% CI 1.0157–1.0400), and 1.0378 (95% CI, 1.0072–1.0692). For bronchitis, the differences in the effect estimates of the effect estimates of PM2.5, PM10, SO2, and NO2 all reached their maximum at lag010, with RR values of 1.0133 (95% CI 1.0025–1.0242), 1.0161(95% CI 1.0085–1.0238), 1.0135 (95% CI 1.0025–1.0247), and 1.1133(95% CI 1.0739–1.1541).

Fig. 1
figure 1

The associations between the concentrations of air pollutants and the number of children hospitalized for ALRI in the single-pollutant model in the Sichuan Province, China, during 2017–2018

Fig. 2 illustrates how each 10 μg/m3 increase in the concentrations of PM2.5, PM10, SO2, and NO2 affects incidences of pneumonia and bronchitis in children of different sexes and ages and in different seasons after adjusting for temperature, relative humidity. The effects of the four air pollutants were significant overall, but there were no statistically significant differences between girls and boys (p > 0.05). SO2 and NO2 caused significantly more hospitalizations for pneumonia or bronchitis in children aged 5–14 years than in those aged 0–1 years and 2–4 years (p < 0.05). Both PM2.5 and PM10 displayed a positive and significant effect on hospitalizations for pediatric pneumonia during the warm season, which statistically differed from the effects observed in the transitional and the cold season (p < 0.05). Similarly, PM2.5, PM10, and NO2 had a positive and significant effect on hospitalizations for bronchiolitis during the warm season, which statistically differed from the effects observed in the transitional and the cold season (p < 0.05) (Table S2).

Fig. 2
figure 2

RR (95% CI) of stratified analyses for each air pollutant in the single-pollutant model in the Sichuan Province, China, during 2017–2018

Based on the largest effect estimates in single-pollutant model, pneumonia PM2.5 is lag4, PM10 and SO2 are lag010,NO2 is lag07; In bronchitis, PM2.5, PM10, SO2 and NO2 are lag010.

Table 2 summarizes the attributable number (AN), attributable risk (AR), and economic cost of ALRI hospitalizations related to PM2.5, PM10, and NO2 pollution in Sichuan Province from 2017 to 2018 using the air quality guidelines of the WHO as a reference standard. Since the daily recorded average concentration of SO2 at lag010 was far lower than the reference concentration (40 μg/m3), we did not measure the AN, AR, or economic costs related to SO2. According to the WHO reference concentrations, the number of ALRI hospitalizations attributable to PM2.5, PM10, and NO2 pollution was 7551, 10,151, and 7575, respectively. The total economic cost of ALRI attributable to PM2.5, PM10, and NO2 pollution during the study period was CNY 2847.06, 3827.27, and 2855.91 million, respectively. Among the specific diseases, pneumonia caused a higher disease burden than bronchitis.

Table 2 The attributable number, attributable risk and economic costs of hospital admissions related to air pollution using WHO’s air quality guidelines in Sichuan Province 2017–2018

Based on the largest effect estimates in single-pollutant model, in ALRI and pneumonia, pneumonia PM2.5 is lag4, PM10 and SO2 are lag010,NO2 is lag07; In bronchitis, PM2.5, PM10, SO2 and NO2 are lag010.

Sensitivity analysis

First, the correlation between particulate matter and hospitalizations for pneumonia or bronchitis weakened and became insignificant when we controlled for gaseous pollutants in the two-pollutant models. When controlling for PM2.5, PM10, and NO2, the association between SO2 and pneumonia hospitalizations alone was strong and significant. Likewise, the association between NO2 and bronchitis hospitalizations was strong and significant when we controlled for PM2.5, PM10, and SO2 (Table 3). Second, the acute effect of air pollutants on hospitalizations did not change substantially after replacing the annual degrees of freedom of the time series with 5, 6, 8 or 9(Table S3). Finally, there were 190 patients whose home was more than 50 km away from the nearest monitoring station, including 146 with pneumonia and 44 with bronchitis. After removing the above data, the effect values for the different pollutants did not change significantly, and our observations from the sensitivity analysis did not significantly differ from the results obtained from the original model (Table S4, Table S5).

Table 3 RR (95% CI) of children hospital admissions per 10 μg/m3 increase in concentrations of pollutants in the single and two-pollutant models

Based on the largest effect estimates in single-pollutant model, pneumonia PM2.5 is lag4, PM10 and SO2 are lag010,NO2 is lag07; In bronchitis, PM2.5, PM10, SO2 and NO2 are lag010.

Discussion

In this study, we performed a time-series analysis to investigate the relationship between air pollutant exposure and the risk of hospitalization for ALRI in children in Sichuan Province, China, and found that short-term exposure to air pollutants increased the risk of hospitalization when data were controlled for confounders such as daily temperature and relative humidity. Stratified analyses revealed that older children were more vulnerable to outdoor air pollutants, as were all children during the hot season. Moreover, the present study identified a serious economic burden of disease due to exposure to air pollutants.

The findings of our single-pollutant model analysis are largely consistent with previously reported findings. These findings support hypothesis that short-term exposure to particulate matter is associated with ALRI hospitalizations. Specifically, it has been shown that exposure to PM2.5 and PM10 increases the risk of ALRI hospitalization in children. In this study, we observed that for each 10 μg/m3 increase in particulate matter concentration of PM2.5, an effect on pneumonia occurred only at lag4, with an RR of 1.0064(95% CI 1.0004–1.0124) and the effect on bronchitis reached its maximum at lag010 with an RR of 1.0133 (95% CI 1.0025–1.0242). By contrast, the effect of PM10 on both pneumonia and bronchitis reached its maximum at lag010, with RR of 1.0168 (95% CI 1.0089–1.0248)、1.0161 (95% CI 1.0085–1.0238), respectively. Multiple studies from several regions, including Taiwan and Jinan in China, found that elevated particulate matter concentrations were associated with increased outpatient visits and hospitalizations for pneumonia [29, 30]. A study conducted in South Korea revealed that the concentration of PM2.5 at lag06 caused the greatest risk of ALRI hospitalization in children [12]. In a case-crossover study of 112,467 children in the U.S., 1-week exposure to PM2.5 was found to be associated with hospitalization for ALRI in children aged 0–2 years and 3–17 years [31]. It has been demonstrated that the mechanisms underlying destructive effects of PM2.5 on the respiratory system involve free radical peroxidative damage, imbalance of intracellular calcium homeostasis, and inflammatory injury [32]. Besides, exposure to particulate matter has been shown to cause adverse effects on lung function in preschool children. A study from China found that for each 10 μg/m3 increase in the concentration of PM2.5, decreases in 25–75% forced expiratory flow (FEF25–75%), forced expiratory volume/forced vital capacity (FEV1/FVC) per second, 75% forced expiratory flow (FEF75%), and forced expiratory volume in 1 second (FEV1) reached their maximum at lag01, at 80.44 mL/s, 35.85%, 78.58 mL/s, and 61.34 mL, respectively [33]. Particulate pollution can also affect the spread and prevalence of respiratory viruses [34, 35]. Among these, respiratory syncytial virus (RSV) is considered to be the most important viral pathogen causing ALRI in young children [36]. A study conducted in Italy revealed that PM10 exposure is associated with increased hospitalizations for RSV bronchiolitis in infants, while the level of PM10 concentrations 2 weeks before hospital admission was closely correlated with an increased risk of hospitalization for RSV bronchiolitis [37].

We also observed that for each 1 μg/m3 increase in SO2 concentration and each 10 μg/m3 increase in NO2 concentration, the effect of SO2 on both pneumonia and bronchitis reached its maximum at lag010 with RRs of 1.0278 (95% CI 1.0157–1.0400), 1.0135 (95% CI 1.0025–1.0247), respectively, and the effect of NO2 on pneumonia reached its maximum at lag07 with an RR of 1.0378 (95% CI 1.0072–1.0692) and the effect on bronchitis reached its maximum at lag010 with an RR of 1.1133 (95% CI 1.0739–1.1541). Similarly, one study on the relationship between short-term air pollutant exposure and childhood ALRI in Nanjing, China found that for each interquartile range (IQR) increase in the concentrations of SO2 and NO2, the estimated cumulative effects of both pollutants reached their maximum at lag05 (5.6%, 2.6–8.6%; 4.1%, 1.2–7.0%) [38]. Moreover, a study in 25 cities in China revealed an association of short-term exposure to SO2 and NO2 with pneumonia and bronchitis [14]. SO2 is readily hydrated into sulfite (SO32−) and bisulfate (HSO3−) ions, and exposure to SO2 can lead to mitochondrial dysfunction in the lungs, which in turn causes cellular disorders and subsequently lung diseases [39, 40]. NO2, a major traffic-related pollutant, has low solubility and can penetrate into the lungs [41]. NO2 can be absorbed throughout the respiratory tract, and incurring damage may occur in the trachea, bronchi, bronchioles, alveolar ducts, and proximal airways, depending on its concentration and dose [42, 43].

In our sex-stratified analysis, we found that while both boys and girls were vulnerable to air pollution, no significant differences were observed between the two groups. This finding is consistent with the results of other studies [44]. Meanwhile, the age-specific analysis in this study showed that SO2 and NO2 had a greater effect on hospitalizations for pneumonia and bronchitis in children aged 2–4 years and 5–14 years than in those aged 0–1 years. A multi-city study from Colombia found that both PM2.5 and NO2 had the greatest impact on emergency room visits for respiratory diseases among children aged 5–9 years [45]. Another study from China also observed that PM2.5, SO2, and NO2 had a greater impact on respiratory diseases in children aged 4–14 years [15]. This observation may be attributed to the fact that children aged 2–4 years and 5–14 years spend more time in school than those aged aged 0–1 years, and data shows that the average concentration of NO2 in the school environment exceeds the WHO guidelines [46]. Notably, an average of 82 asthma attacks per school per year could be avoided by reducing outdoor NO2 concentrations [47], since school-age children begin to participate in more outdoor activities and are exposed to more traffic-related pollutants. At the same time, breastfeeding improves the infant’s resistance and, with relatively quiet infants, the physical activity of larger children leads to deeper and more frequent aspiration [48, 49]. In the present study, a specific season analysis showed that compared with the transitional and the cold season, PM2.5 and PM10 had greater adverse effects on childhood pneumonia during the warm season, as had PM2.5, PM10, and NO2 on childhood bronchitis. Similarly, Cheng et al. found that children are more vulnerable to NO2 during the warm season [29]. Lv et al. also found a greater impact on child hospital admissions during the warm season [30]. The above seasonal differences could be attributed to the fact that the correlation between personal air pollutant exposure and ambient air pollutant concentrations is higher in summer than in winter, and there are differences in personal exposure between different seasons [50]. Alternatively, the concentrations of different components of particulate matter may vary among seasons [51]. The discrepancy could also be explained by the fact that high temperatures lead to more NO2 emissions, especially in cities with larger traffic volumes, and that higher concentrations of NO2 can occur in economically developed areas [52]. In this case, the spatial and temporal distribution of air pollutants could potentially be affected.

The present study not only discusses the relationship between air pollution and ALRI hospitalization in children, but also assesses the economic burden caused by ALRI hospitalizations; it therefore has important implications for improving air quality and preventing respiratory diseases. During the study period, the greatest economic burden of hospitalizations for ALRI was caused by PM10, followed by NO2, while the smallest burden was caused by PM2.5. Likewise, a study conducted in Guiyang, China, found that PM10 caused higher hospitalization costs than PM2.5 [53]. This finding could be explained by the fact that PM2.5 is included in PM10 concentrations, and that the two types of particulate matter are highly correlated (r = 0.9442). An earlier study that employed machine learning to verify the accuracy of linking particulate matter concentrations to upper respiratory infections found higher accuracy for PM10 than PM2.5 [54]. Overall, ALRI pose a serious burden on children and their families.

We observed that the association of PM2.5 with ALRI hospitalizations weakened and became insignificant after adjusting for SO2 and NO2 in our two-pollutant models. Similarly, Zheng et al. found that after controlling for SO2 and NO2, the estimates associated with PM2.5 in the two-pollutant models were not significant (PM2.5, 1.50, 95% CI, 0.35–2.66; PM2.5 adjusted for SO2, 0.17, 95% CI, − 2.55–0.43; PM2.5 adjusted for NO2, 0.17, 95% CI, − 1.47–1.85) [55]. Moreover, the authors observed that the association of SO2 and NO2 with ALRI hospitalizations became partially statistically insignificant after mutual adjustments. This observation is similar to the results of a number of previous studies [56]. There may be an interaction between SO2 and NO2; moreover, both particulate matter and NO2 are traffic-related pollutants and are highly correlated [44]. Thus, it is difficult to determine the individual effects of each pollutant. The above results imply that the effects of pollutants are not simply superimposed. Instead, collinearity effects between different pollutants may occur and may have a synergistic effect on the acute lower respiratory tract, thereby affecting the authenticity of the model.

The present study has the following three advantages: First, this is the first study on disease burden due to air pollutant exposure in children aged 14 years and younger in nine cities and prefectures in Sichuan Province. Second, cities and prefectures from the Sichuan Basin as well as from the western Sichuan Plateau were selected, and our sample therefore represents the overall situation in Sichuan Province. Third, the IDW interpolation method was used to construct a high spatial resolution for the estimation of pollutant concentrations, which improved the spatial accessibility of pollutants. However, this study also has some limitations. First, this is an ecological study, and our findings may be slightly inconsistent with the real situation, due to the presence of confounding factors. Second, we collected only 2 years’ worth of data for our analysis of the associations between air pollutants and hospitalizations for ALRI in children, which may have caused some instability in the model we used [57].

Conclusion

Based on a time-series analysis, this study assessed the short-term effects of short-term air pollutant exposure on the daily number of hospitalizations for ALRI in children in Sichuan Province between 2017 and 2018, and found that elevated daily average concentrations of PM2.5, PM10, SO2, and NO2 increase hospitalizations, showing lagged effects. Moreover, we found that SO2 and NO2 have more significant effects on older children and that the economic burden due to ALRI can be partially attributed to excessive pollutant exposure. In order to protect children’s health, authorities in Sichuan Province should take effective measures to reduce the emission of harmful substances in the air.

Abbreviations

PM2.5:

Particulate matter < 2.5 μm in aerodynamic diameter

SO2:

Sulfur dioxide

NO2:

Nitrogen dioxide

ALRI:

Acute lower respiratory tract infections

IDW:

Inverse distance weighting

GAM:

Generalized additive model

SD:

Standard deviation

RR:

Relative risk

CI:

Confidence interval

FEF25–75% :

Forced expiratory flow

FEV1/FVC:

Forced expiratory volume/forced vital capacity

FEF75% :

75% forced expiratory flow

RSV:

Respiratory syncytial virus.

References

  1. Fuller R, Landrigan PJ, Balakrishnan K, Bathan G, Bose-O'Reilly S, Brauer M, et al. Pollution and health: a progress update. Lancet Planet Health. 2022;6(6):e535–47.

    Article  Google Scholar 

  2. Stieb DM, Chen L, Eshoul M, Judek S. Ambient air pollution, birth weight and preterm birth: a systematic review and meta-analysis. Environ Res. 2012;117:100–11.

    Article  CAS  Google Scholar 

  3. Pedersen M, Giorgis-Allemand L, Bernard C, Aguilera I, Andersen AM, Ballester F, et al. Ambient air pollution and low birthweight: a European cohort study (ESCAPE). Lancet Respir Med. 2013;1(9):695–704.

    Article  CAS  Google Scholar 

  4. Weinmayr G, Romeo E, De Sario M, Weiland SK, Forastiere F. Short-term effects of PM10 and NO2 on respiratory health among children with asthma or asthma-like symptoms: a systematic review and meta-analysis. Environ Health Perspect. 2010;118(4):449–57.

    Article  CAS  Google Scholar 

  5. Mok PLH, Antonsen S, Agerbo E, Brandt J, Geels C, Christensen JH, et al. Exposure to ambient air pollution during childhood and subsequent risk of self-harm: a national cohort study. Prev Med. 2021;152(Pt 1):106502.

    Article  Google Scholar 

  6. Milojevic A, Dutey-Magni P, Dearden L, Wilkinson P. Lifelong exposure to air pollution and cognitive development in young children: the UK millennium cohort study. Environ Res Lett. 2021;16(5):055023.

    Article  CAS  Google Scholar 

  7. GBD 2015 LRI Collaborators. Estimates of the global, regional, and national morbidity, mortality, and aetiologies of lower respiratory tract infections in 195 countries: a systematic analysis for the Global Burden of Disease Study 2015. Lancet Infect Dis. 2017;17(11):1133–61.

    Article  Google Scholar 

  8. GBD 2019 Diseases and Injuries Collaborators. 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(10258):1204–22.

    Article  Google Scholar 

  9. World Health Organization. More than 90% of the world’s children breathe toxic air every day. https://www.who.int/news/item/29-10-2018-more-than-90-of-the-worlds-children-breathe-toxic-air-every-day. Accessed 4 July 2022.

  10. Larson PS, Espira L, Glenn BE, Larson MC, Crowe CS, Jang S, et al. Long-term PM2.5 exposure is associated with symptoms of acute respiratory infections among children under five years of age in Kenya, 2014. Int J Environ Res Public Health. 2022;19(5):2525.

    Article  Google Scholar 

  11. Agudelo-Castañeda DM, Calesso Teixeira E, Alves L, Fernández-Niño JA, Rodríguez-Villamizar LA. Monthly-term associations between air pollutants and respiratory morbidity in South Brazil 2013-2016: a Multi-City, time-series analysis. Int J Environ Res Public Health. 2019;16(20):3787.

    Article  Google Scholar 

  12. Oh J, Han C, Lee DW, Jang Y, Choi YJ, Bae HJ, et al. Short-term exposure to fine particulate matter and hospitalizations for acute lower respiratory infection in Korean children: a time-series study in seven metropolitan cities. Int J Environ Res Public Health. 2020;18(1):144.

    Article  Google Scholar 

  13. Cheng J, Su H, Xu Z. Intraday effects of outdoor air pollution on acute upper and lower respiratory infections in Australian children. Environ Pollut. 2021;268(Pt A):115698.

    Article  CAS  Google Scholar 

  14. Xu H, Wang X, Tian Y, Tian J, Zeng Y, Guo Y, et al. Short-term exposure to gaseous air pollutants and daily hospitalizations for acute upper and lower respiratory infections among children from 25 cities in China. Environ Res. 2022;212(Pt D):113493.

    Article  CAS  Google Scholar 

  15. Yang H, Yan C, Li M, Zhao L, Long Z, Fan Y, et al. Short term effects of air pollutants on hospital admissions for respiratory diseases among children: a multi-city time-series study in China. Int J Hyg Environ Health. 2021;231:113638.

    Article  CAS  Google Scholar 

  16. Wu Y, Jin T, He W, Liu L, Li H, Liu C, et al. Associations of fine particulate matter and constituents with pediatric emergency room visits for respiratory diseases in Shanghai. China Int J Hyg Environ Health. 2021;236:113805.

    Article  CAS  Google Scholar 

  17. Liang Z, Meng Q, Yang Q, Chen N, You C. Size-specific particulate matter associated with acute lower respiratory infection outpatient visits in children: a counterfactual analysis in Guangzhou. China Front Public Health. 2021;9:789542.

    Article  Google Scholar 

  18. Wang ZB, Ren L, Lu QB, Zhang XA, Miao D, Hu YY, et al. The impact of weather and air pollution on viral infection and disease outcome among pediatric pneumonia patients in Chongqing, China, from 2009 to 2018: a prospective observational study. Clin Infect Dis. 2021;73(2):e513–22.

  19. Liu Y, Wang Y, Dong J, Wang J, Bao H, Zhai G. Association between air pollution and emergency department visits for upper respiratory tract infection in Lanzhou. China Environ Sci Pollut Res Int. 2022;29(19):28816–28.

    Article  CAS  Google Scholar 

  20. Ma Y, Yue L, Liu J, He X, Li L, Niu J, et al. Association of air pollution with outpatient visits for respiratory diseases of children in an ex-heavily polluted northwestern city, China. BMC Public Health. 2020;20(1):816.

    Article  Google Scholar 

  21. Chen Y, Yang Y, Yao Y, Wang X, Xu Z. Spatial and dynamic effects of air pollution on under-five children's lower respiratory infections: an evidence from China 2006 to 2017. Environ Sci Pollut Res Int. 2022;29(17):25391–407.

    Article  Google Scholar 

  22. Liu L, Chen YY, Wu T, Li HM. The drivers of air pollution in the development of western China: the case of Sichuan province. J Clean Prod. 2018;197(PTa1):1169–76.

    Article  CAS  Google Scholar 

  23. Zhang XY, Wang YQ, Niu T, Zhang XC, Gong SL, Zhang YM, et al. Atmospheric aerosol compositions in China: spatial/temporal variability, chemical signature, regional haze distribution and comparisons with global aerosols. Atmos Chem Phys. 2012;12(258):779–99 (in Chinese).

    Article  CAS  Google Scholar 

  24. Pu X, Wang L, Chen L, Pan J, Tang L, Wen J, et al. Differential effects of size-specific particulate matter on lower respiratory infections in children: a multi-city time-series analysis in Sichuan. China Environ Res. 2021;193:110581.

    Article  CAS  Google Scholar 

  25. Altman DG, Bland JM. Interaction revisited: the difference between two estimates. BMJ. 2003;326(7382):219.

    Article  Google Scholar 

  26. Zhang P, Zhou X. Health and economic impacts of particulate matter pollution on hospital admissions for mental disorders in Chengdu. Southwestern China Sci Total Environ. 2020;733:139114.

    Article  CAS  Google Scholar 

  27. World Health Organization. WHO global air quality guidelines: particulate matter (P-M2.5 and PM10), ozone, nitrogen dioxide, sulfur dioxide and carbon monoxide. https://www.who.int/publications-detail-redirect/9789240034228. Accessed 10 July 2022.

  28. Liu Y, Pan J, Zhang H, Shi C, Li G, Peng Z, et al. Short-term exposure to ambient air pollution and asthma mortality. Am J Respir Crit Care Med. 2019;200(1):24–32.

    Article  CAS  Google Scholar 

  29. Cheng CY, Cheng SY, Chen CC, Pan HY, Wu KH, Cheng FJ. Ambient air pollution is associated with pediatric pneumonia: a time-stratified case-crossover study in an urban area. Environ Health. 2019;18(1):77.

    Article  Google Scholar 

  30. Lv C, Wang X, Pang N, Wang L, Wang Y, Xu T, et al. The impact of airborne particulate matter on pediatric hospital admissions for pneumonia among children in Jinan, China: a case-crossover study. J Air Waste Manag Assoc. 2017;67(6):669–76.

    Article  CAS  Google Scholar 

  31. Horne BD, Joy EA, Hofmann MG, Gesteland PH, Cannon JB, Lefler JS, et al. Short-term elevation of fine particulate matter air pollution and acute lower respiratory infection. Am J Respir Crit Care Med. 2018;198:759–66.

    Article  Google Scholar 

  32. Xing YF, Xu YH, Shi MH, Lian YX. The impact of PM2.5 on the human respiratory system. J Thorac Dis. 2016;8(1):E69–74.

    Google Scholar 

  33. Yang XY, Wen B, Han F, Wang C, Zhang SP, Wang J, et al. Acute effects of individual exposure to fine particulate matter on pulmonary function in schoolchildren. Biomed Environ Sci. 2020;33(9):647–59.

    Google Scholar 

  34. Vandini S, Corvaglia L, Alessandroni R, Aquilano G, Marsico C, Spinelli M, et al. Respiratory syncytial virus infection in infants and correlation with meteorological factors and air pollutants. Ital J Pediatr. 2013;39(1):1.

    Article  Google Scholar 

  35. Wrotek A, Badyda A, Czechowski PO, Owczarek T, Dąbrowiecki P, Jackowska T. Air Pollutants' concentrations are associated with increased number of RSV hospitalizations in polish children. J Clin Med. 2021;10(15):3224.

    Article  CAS  Google Scholar 

  36. Nair H, Nokes DJ, Gessner BD, Dherani M, Madhi SA, Singleton RJ, et al. Global burden of acute lower respiratory infections due to respiratory syncytial virus in young children: a systematic review and meta-analysis. Lancet. 2010;375(9725):1545–55.

    Article  Google Scholar 

  37. Carugno M, Dentali F, Mathieu G, Fontanella A, Mariani J, Bordini L, et al. PM10 exposure is associated with increased hospitalizations for respiratory syncytial virus bronchiolitis among infants in Lombardy. Italy Environ Res. 2018;166:452–7.

    Article  CAS  Google Scholar 

  38. Zhu L, Ge X, Chen Y, Zeng X, Pan W, Zhang X, et al. Short-term effects of ambient air pollution and childhood lower respiratory diseases. Sci Rep. 2017;7(1):4414.

    Article  Google Scholar 

  39. Meng Z, Qin G, Zhang B, Bai J. DNA damaging effects of sulfur dioxide derivatives in cells from various organs of mice. Mutagenesis. 2004;19(6):465–8.

    Article  CAS  Google Scholar 

  40. Qin G, Wang J, Sang N. Sulfur dioxide inhibits expression of mitochondrial oxidative phosphorylation genes encoded by both nuclear DNA and mitochondrial DNA in rat lungs. Environ Sci Pollut Res Int. 2017;24(3):2527–34.

    Article  CAS  Google Scholar 

  41. Schraufnagel DE, Balmes JR, Cowl CT, De Matteis S, Jung SH, Mortimer K, et al. Air pollution and noncommunicable diseases: a review by the forum of international respiratory Societies' environmental committee, part 1: the damaging effects of air pollution. Chest. 2019;155(2):409–16.

    Article  Google Scholar 

  42. Goldstein E, Peek NF, Parks NJ, Hines HH, Steffey EP, Tark-ington B. Fate and distribution of inhaled nitrogen dioxide in rhesus monkeys. Am Rev Respir Dis. 1977;115(3):403–12.

    CAS  Google Scholar 

  43. Moldeus P. Toxicity induced by nitrogen dioxide in experimental animals and isolated cell systems. Scand J Work Environ Health. 1993;19:28–36.

    CAS  Google Scholar 

  44. Bai L, Su X, Zhao D, Zhang Y, Cheng Q, Zhang H, et al. Exposure to traffic-related air pollution and acute bronchitis in children: season and age as modifiers. J Epidemiol Community Health. 2018;72(5):426–33.

    Article  Google Scholar 

  45. Odríguez-Villamizar LA, Rojas-Roa NY, Blanco-Becerra LC, Herrera-Galindo VM, Fernández-Niño JA. Short-term effects of air pollution on respiratory and circulatory morbidity in Colombia 2011-2014: a Multi-City, time-series analysis. Int J Environ Res Public Health. 2018;15(8):1610.

    Article  Google Scholar 

  46. Salonen H, Salthammer T, Morawska L. Human exposure to NO2 in school and office indoor environments. Environ Int. 2019;130:104887.

    Article  CAS  Google Scholar 

  47. Guerriero C, Chatzidiakou L, Cairns J, Mumovic D. The economic benefits of reducing the levels of nitrogen dioxide (NO2) near primary schools: the case of London. J Environ Manag. 2016;181:615–22.

    Article  CAS  Google Scholar 

  48. Duijts L, Jaddoe VW, Hofman A, Moll HA. Prolonged and exclusive breastfeeding reduces the risk of infectious diseases in infancy. Pediatrics. 2010;126(1):e18–25.

    Article  Google Scholar 

  49. World Health Organization. World health statistics 2013. Geneva: World Health Organization Press; 2013.

    Google Scholar 

  50. Sarnat JA, Brown KW, Schwartz J, Coull BA, Koutrakis P. Ambient gas concentrations and personal particulate matter exposures: implications for studying the health effects of particles. Epidemiology. 2005;16(3):385–95.

    Article  Google Scholar 

  51. Ueda K, Yamagami M, Ikemori F, Hisatsune K, Nitta H. Associations between fine particulate matter components and daily mortality in Nagoya. Japan J Epidemiol. 2016;26(5):249–57.

    Article  Google Scholar 

  52. Xiao ZY, Jiang H. Characteristics of atmospheric NO2 over Sichuan basin. China Environ Sci. 2011;31(11):1782–8 (in Chinese).

    CAS  Google Scholar 

  53. Zhou H, Wang T, Zhou F, Liu Y, Zhao W, Wang X, et al. Ambient air pollution and daily hospital admissions for respiratory disease in children in Guiyang, China. Front Pediatr. 2019;7:400.

    Article  Google Scholar 

  54. Chen MJ, Yang PH, Hsieh MT, Yeh CH, Huang CH, Yang CM, et al. Machine learning to relate PM2.5 and PM10 concentrations to outpatient visits for upper respiratory tract infections in Taiwan: a nationwide analysis. World J Clin Cases. 2018;6(8):200–6.

    Article  Google Scholar 

  55. Zheng PW, Wang JB, Zhang ZY, Shen P, Chai PF, Li D, et al. Air pollution and hospital visits for acute upper and lower respiratory infections among children in Ningbo, China: a time-series analysis. Environ Sci Pollut Res Int. 2017;24(23):18860–9.

    Article  CAS  Google Scholar 

  56. Chen J, Shi C, Li Y, Ni H, Zeng J, Lu R, et al. Effects of short-term exposure to ambient airborne pollutants on COPD-related mortality among the elderly residents of Chengdu city in Southwest China. Environ Health Prev Med. 2021;26(1):7.

    Article  CAS  Google Scholar 

  57. Strickland MJ, Gass KM, Goldman GT, Mulholland JA. Effects of ambient air pollution measurement error on health effect estimates in time-series studies: a simulation-based analysis. J Expo Sci Environ Epidemiol. 2015;25(2):160–6.

    Article  CAS  Google Scholar 

Download references

Acknowledgments

Not applicable.

Availability of data and material

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Funding

This study was funded by the Calculation and Research of Total Health Expenditure in Sichuan Province in 2021 (Grant NO. 301021002).

Author information

Authors and Affiliations

Authors

Contributions

HY initiated the idea for the study and was involved in writing of the original draft. JW and GX developed the formal analysis and software. LCW and LJ was involved in data curation. YL was involved in reviewing and editing. YL was the PI for the fund and designed ideas of research. All the authors participated in revising manuscript according to the comments of reviewers and approved the final report. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lian Yang.

Ethics declarations

Ethics approval and consent to participate

All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki. This study was approved by the Ethics Committee of Hospital of Chengdu University of Traditional Chinese Medicine (approval no. 2020KL-001).

Consent for publication

Not applicable.

Competing interests

No financial or non-financial benefits have been received or will be received from any party related directly or indirectly to the subject of this article.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1 Table S1

. Pearson correlation coefficients between daily meteorological factors and air pollutants for Sichuan (2017–2018). Table S2. RR (95% CI) of stratified analyses for each air pollutant in the single-pollutant model. Based on the largest effect estimates in single-pollutant model, pneumonia PM2.5 is lag4, PM10 and SO2 are lag010,NO2 is lag07; In bronchitis, PM2.5, PM10, SO2 and NO2 are lag010. PM2.5 fine particulate matter, PM10 inhalable particulate matter, SO2 sulfur dioxide, NO2 nitrogen dioxide, CI confidence interval. Table S3. Association between air pollutants (10 μg/m3 increase) and the daily hospitalization in children with ALRI by degrees of freedom per year. Table S4. Associations between air pollutants (every 10 μg/m3 increase in the later period of retention) and hospitalization in children with pneumonia (Eliminate the data from the home address to the monitoring station greater than 50 km). Table S5. Associations between air pollutants (every 10 μg/m3 increase in the later period of retention) and hospitalization in children with bronchiolitis (Eliminate the data from the home address to the monitoring station greater than 50 km.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

He, Y., Jiang, W., Gao, X. et al. Short-term effects and economic burden of air pollutants on acute lower respiratory tract infections in children in Southwest China: a time-series study. Environ Health 22, 6 (2023). https://0-doi-org.brum.beds.ac.uk/10.1186/s12940-023-00962-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12940-023-00962-3

Keywords