- Research
- Open access
- Published:

# Lagging effects and prediction of pollutants and their interaction modifiers on influenza in northeastern China

*BMC Public Health*
**volume 23**, Article number: 1826 (2023)

## Abstract

### Background

Previous studies have typically explored the daily lagged relations between influenza and meteorology, but few have explored seasonally the monthly lagged relationship, interaction and multiple prediction between influenza and pollution. Our specific objectives are to evaluate the lagged and interaction effects of pollution factors and construct models for estimating influenza incidence in a hierarchical manner.

### Methods

Our researchers collect influenza case data from 2005 to 2018 with meteorological and contaminative factors in Northeast China. We develop a generalized additive model with up to 6 months of maximum lag to analyze the impact of pollution factors on influenza cases and their interaction effects. We employ LASSO regression to identify the most significant environmental factors and conduct multiple complex regression analysis. In addition, quantile regression is taken to model the relation between influenza morbidity and specific percentiles (or quantiles) of meteorological factors.

### Results

The influenza epidemic in Northeast China has shown an upward trend year by year. The excessive incidence of influenza in Northeast China may be attributed to the suspected primary air pollutant, NO_{2}, which has been observed to have overall low levels during January, March, and June. The Age 15–24 group shows an increase in the relative risk of influenza with an increase in PM_{2.5} concentration, with a lag of 0–6 months (ERR 1.08, 95% CI 0.10–2.07). In the quantitative analysis of the interaction model, PM_{10} at the level of 100–120 μg/m^{3}, PM_{2.5} at the level of 60–80 μg/m^{3}, and NO_{2} at the level of 60 μg/m^{3} or more have the greatest effect on the onset of influenza. The GPR model behaves better among prediction models.

### Conclusions

Exposure to the air pollutant NO_{2} is associated with an increased risk of influenza with a cumulative lag effect. Prioritizing winter and spring pollution monitoring and influenza prediction modeling should be our focus.

## Introduction

Influenza is an acute viral respiratory illness in humans, usually characterized by fever, headache, muscle pain, weakness, nasal congestion, sore throat and cough. Seasonal epidemics of influenza viruses can spread rapidly and cause significant morbidity and mortality worldwide [1,2,3]. Globally, influenza is estimated to cause 3 to 5 million severe cases and 290,000 to 650,000 deaths related to respiratory infections every year [4]. While certain non-pharmaceutical interventions that can effectively control influenza in the early stages exist, such as the use of masks, hand washing, and other hygiene measures, or even closing schools, influenza is still a big threat to human health [5]. As evidenced by the trends of influenza incidence in such countries as Europe, the United States and Japan [6, 7], the situation remains grim and global preventive measures have little impact on the trend of influenza outbreaks. Currently few models can effectively predict influenza outbreaks. By monitoring influenza morbidity indicators on a daily basis after environmental factors that tend to predict outbreaks can be identified and modelled, we could predict influenza virus outbreaks in advance to reduce or even prevent influenza with associated costs.

In early studies we can assume that air pollution may contribute to influenza-induced morbidity. An epidemiological investigation suggests that particulate matter ≤ 10 μm (PM_{10}) and ozone (O_{3}) should be considered when forecasting the incidence of influenza [8, 9]. Influenza viruses have been detected in polluted waters, possibly originating from bird excretion carrying the virus, as per some studies [10]. As shown in the 2002–2003 SARS pandemic and the 2009 H_{1}N_{1} influenza pandemic, influenza viruses are mainly transmitted through respiratory droplets. So air pollutants such as particulate matter (PM) and carbon monoxide (CO) may influence the transmission and prevalence of influenza viruses [11, 12]. In addition, secondary human-to-human transmission may occur and the outbreak may lead to the closure of schools and workplaces. Previous studies have investigated how meteorological factors facilitate the transmission of influenza among regions worldwide. A study in the UK showed that influenza viruses prefer low temperatures in temperate regions [13], while researchers in Canada have observed that the increase of influenza viruses is associated with low temperatures and high relative humidity [14]. The different subtypes of influenza viruses that have infected humans in recent decades include the H_{10}N_{8}, H_{5}N_{6} and H_{9}N_{2}, most of which were firstly reported in China [15].

Numerous studies, both domestically and internationally, have investigated the correlation between the seasonal distribution patterns of influenza and meteorological and pollution factors. In China, Yuzhou Zhang et al. [16] explored the effect of different meteorological factors on influenza incidence in Shanghai by developing a distribution lagged nonlinear model (DLNM). Some researches indicated that in the multi-day lag model, there was a statistically significant correlation between SO_{2}, NO_{2} and O_{3} concentrations and influenza risk between lags 0 and 1 [17]. Similarly, in Nanjing, Lei Huang et al. [18] found that PM_{2.5} and NO_{2} were associated with an increase in influenza cases. Air pollutants significantly affect the susceptibility of human respiratory epithelial cells to influenza virus infection by increasing virus attachment and entry.

The overall objective of this study is to explore influenza epidemic characteristics, lagged effects of pollutants and develop models suitable for predicting influenza virus outbreaks. Our specific objectives are to: a) screen environmental predictors of influenza outbreaks; b) evaluate the lagged and interaction effects of pollution factors; c) construct models for estimating influenza incidence in a hierarchical manner, selecting appropriate models for different characteristics.

## Materials and methods

Figure 1 shows the geographical location of the study area—Heilongjiang, Jilin and Liaoning provinces, which lie between 120° and 135° E and 40° and 53° N latitude in China. These three provinces are located in Northeast China with a medium level of economic development and population density.

We collected influenza case surveillance data from the National Public Health Data Centre of China between 2005 and 2018. All patients are diagnosed according to the criteria of influenza management issued by the Ministry of Health of the People's Republic of China. We obtained the corresponding daily weather data including air temperature, dew point temperature etc. from the China Meteorological Data Sharing Service. Pollutant information is originally from the National Oceanic and Atmospheric Administration (NOAA) including CO, NO_{2}, O_{3} etc.

## Statistical methods

To address missing values in the influenza epidemic and meteorological pollution data, we use multiple interpolation to fill them. LASSO regression analysis is used for feature selection in response to the effects of meteorological and pollution factors. We develop quantile regression models and generalized additive models [13] with a maximum lag of 6 months to assess the extreme effects of pollution and meteorological factors on influenza cases, lags, and interactions between pollutants, respectively. Finally, we make prediction by complex regression models. All analysis in our study is performed in R software (version 4.1.3).

### The screening of lasso regression with environmental variables

There are N sets of observations, each consisting of a total response variable *y* and *p* associated characteristic variables *x*_{i} = (*x*_{i1}, …, *x*_{ip}) ^{T}. A linear regression model can be set as follows:

where* β*_{0} and *β* = (*β*_{1}, *β*_{2}, …, *β*_{p}) are unknown parameters and *e*_{i} is the error term. The introduction of some variables in practical problems not only complicates the calculation, but also risks increasing the data covariance, thus affecting the model fit. We can use lasso regression to estimate the parameters by solving the following problem:

where \(\parallel\beta\parallel_1=\sum_{j=1}^p\left|\beta_j\right|\) is the \({l}_{1}\) parametrization of \(\beta\) and *t* is the specified tuning parameter. Overall, the lasso method improves the overall prediction accuracy, and the inclusion of the constraint term compresses the coefficients of some of the eigenvariables in the model to zero, thus enabling the selection of important variables among the many eigenvariables.

###
*The lagging and interaction effect of generalized additive model *[13]

The models are listed as follows (Model 1):

Here, *Y*_{t} is the number of monthly counts of influenza cases in monthly t; *α1* is the intercept of the whole model; *S* () is a smoothing function, and the penalty spline method is often used to smooth the parameters; *M* represents the estimated environmental variable related to influenza; *β* is the regression coefficients. The optimal degrees of freedom (*df*) for the spline function are estimated by Akaike information criterion for Poisson (AIC) and Minimum partial regression coefficient (PACF_{min}) criteria.

Subsequently, we explore the interaction of pollutants on the prevalence of influenza. The model can be written as follows (Model 2):

*α*2 is the intercept; *X*_{1} indicates two of the interaction factors whereas *X*_{3} denote the other one; *S*(*X*_{1}, *X*_{2}) is a spline function of the interaction between the parameters *X*_{1} and *X*_{2}*.* M is the meteorological factors.

### The establishment of Gradient boosting regression tree (GBRT) and Gaussian process regression (GPR) model

If the input training set is: T = {(× 1, y1), (× 2, y2), …, (xN, yN)}, the training samples i = 1, 2, …, N, the number of iteration rounds t = 1, 2, …, T, and the loss function is L, then the GBRT algorithm is divided into the following three steps:

First, initialize the weak learner: \(f_0\left(x\right)=\mathrm{argmin}\;\sum\limits_{i=1}^N\;L\left(y_i,\;c\right)\)

Next, calculate the negative gradient r_{ti} and the output of each leaf node region of the regression tree, R_{tm} output value c_{tm}, and update the strong learner:

Finally, get strong learners:

From a function space perspective, a Gaussian process [19] is defined to describe the distribution of the function (f(x)). The GP is the set of any finite number of random variables that have a joint Gaussian distribution, and its properties are determined entirely by the mean and covariance functions, that is:

where x,x' ∈ R are arbitrary random variables. Thus GP can be defined as f (x) ~ GP(m(x), k(x, x')) and the mean function is generally taken to be 0 (m(x) = 0). For the regression model as follows:

where *x* is the input value, *f* is the function value, and *y* is the observation plus the observation affected by noise, if noise \(\varepsilon \sim \mathrm{N}\left(0,{\upsigma }_{\mathrm{n}}^{2}\right)\) yields a priori distribution of the observation *y* as follows:

## Results

### Influenza surveillance in Northeast China

From 2005 to 2018, a total of 32,989 influenza cases were reported in the three eastern provinces of China, showing an increasing trend every year (Table 1 and Fig. 1). Heilongjiang province exhibited a significantly high level of epidemic in the first few years, followed by Liaoning province, which has consistently been the main epidemic area since then and had 14,921 reported cases of influenza by 2018. Young children aged 5–14 years and young adults aged 25–59 years had the highest incidence of influenza, accounting for 55.35% of all reported cases (Table 1). Significant differences in the incidence of influenza were observed in terms of seasonality, age, and region (*P* < 0.05).

### The screening and extreme effect for meteorological and pollutant factors among influenza prevalence

We apply a fivefold cross-validation to select a model with small and stable error fluctuations and a parameter λ of 28.3327. The results of the runs are shown in Fig. 2 and Table S1. After the initial selection of the lasso method, six of the variable variables including air temperature (AT), dew point temperature (DPT), sea level pressure (SLP), NO_{2}, PM_{10} and PM_{2.5} are selected and the coefficients of the other independent variables are contracted to zero.

### Exposure–response relationships for pollutants with different lag times

From the line graph of Figure S1 regression coefficients, we can see that SLP does not fluctuate at different levels, while PM_{2.5} fluctuates at different concentrations, but the overall effect increases with the quantile. In addition, PM_{10}, NO_{2} and AT has a negative effect on influenza incidence at different levels, and the overall effect also increases with the quantile. Significant differences were observed between the age groups of 25–59 years and 60 + year. In Table 2 and Fig. 3, in the single-pollutant model we find a negative association between short-term exposure to NO_{2} (within 1 month) and monthly influenza incidence (ERR-2.68% (-4.72%, -0.60%)), and this implies that low levels of NO_{2} may be the most responsible air pollutant for excess influenza incidence. The Age 25–59 years group is the most susceptible to NO_{2}, followed by the Age 0–4 years group, and the ERR increases with lag in both groups, with essentially no lag in the 60 + years age group. NO_{2} showed a positive correlation with influenza incidence in patients aged 15–24 and 25–59 years at a 3-month lag. At a 5-month lag, there was a positive correlation between influenza incidence and NO_{2} in patients aged 0–4, 15–24, and 25–59 years. At a maximum lag time of 6 months, the influenza incidence in patients aged 15–24 shows a positive correlation with PM_{2.5}.

### Interaction and comparison of multiple-pollutant model

We develop a multi-pollutant model with a single-day lagged ERR maximum and significance test, yielding a lagged day of 5 days for PM_{10}, PM_{2.5} and 4 days for NO_{2}. Figure S2 revealed that PM_{10} has interactive effects with PM_{2.5} and NO_{2} on influenza incidence. Additionally, there is a weak positive correlation between pollutants and the risk of incidence. Moreover, ambient temperature (AT) exhibits positive correlation with the risk of incidence at low temperatures and inverse correlation at high temperatures, while dew point temperature (DPT) shows an inverse trend compared to AT. Figure 4 indicated a non-linear effect of pollutants on influenza onset, with PM_{10} levels of 100–120 μg/m^{3} and PM_{2.5} levels of 60–80 μg/m^{3}, and NO_{2} levels above 60 μg/m^{3} exhibiting the greatest impact. Results from the statistical test presented in Table 3 suggest that the PM_{10} and PM_{2.5} interaction model is better (R^{2} = 99.1%). These findings demonstrate the significant impact of air pollutants on influenza onset.

From the comparison of the parameters of the two modelling approaches in Table 4, the model fit is the best in Liaoning Province among different regions (R^{2} > 70%), and the model fit is the best among Age 25–59 groups, while the GPR model shows the same fit as the GBRT model.

## Discussion

Based on this study, it is determined that influenza epidemics in the northeastern region of China exhibit pronounced seasonality of winter-spring and demonstrate an upward trend annually. Prior to 2018, the province of Liaoning continued to be the primary epicenter of influenza outbreaks within the northeastern region, with a higher prevalence of cases observed among young children and adolescents, which could be attributed to their relatively weaker immune systems rendering them more susceptible to influenza viral infections [20].

Our single-pollutant lagged modeling yields intriguing results, indicating that lower concentrations of NO_{2} in January, March, and June may be the primary contributing factor to excessive influenza incidence. Exposure to NO_{2} can lead to reduced virus-specific immunity and increased cellular inflammation, potentially causing the onset of influenza virus, regardless of whether it occurs before or after respiratory virus infection. The relative risk of influenza associated with NO_{2} exposure increased with higher NO_{2} concentrations in the age groups of 0–4, 15–24, and 25–59. This suggests that young adults may have a higher susceptibility to influenza under NO_{2} exposure due to the rapid release of immune cells stimulated by the virus, resulting in an increased relative release of immune cells. The disruption of immune cell homeostasis caused by the rapid release of immune cells stimulated by influenza viruses may explain the higher susceptibility of young adults to influenza under NO_{2} exposure. This could be the result of a relative increase in immune cell release in this age group [21]. In the age group of 15–24, the relative risk of influenza associated with PM_{2.5} exposure increased with higher PM_{2.5} concentrations at lag 0–6 months. This suggests that long-term exposure to high levels of PM_{2.5} beyond 6 months may increase the risk of influenza. This was similar to the results of a monitoring study that found: Age 0–4 were significantly susceptible to PM_{10} and NO_{2}; Age 5–14 were significantly susceptible to PM_{2.5} and PM_{10}; and Age 15–24 were significantly susceptible to all air pollutants analyzed [18]. During the remaining lag months with low concentrations of PM_{2.5}, PM_{10} and NO_{2}, the onset of influenza is mainly attributed to factors other than air pollution, as indicated by ERR < 0. This may be due to low levels of external pollution, which increase population activity during high influenza season (spring, winter) and result in aggregated activity, thereby enhancing the risk of influenza.

The longitudinal study shows that the highest overall risk (ERR) of influenza onset is observed for pollutants with a lag of 5–6 months, indicating that long-term exposure to pollutants may primarily promote influenza onset. However, in the age groups of 15–24 and 60 + , NO_{2} exposure is associated with a higher risk of influenza onset at an early stage. The lagged effects of PM_{2.5}, PM_{10} and NO_{2} are characterized by a bimodal distribution, with a significant decrease in the risk of influenza onset during the first 2 months of exposure. However, the risk of influenza onset is found to be higher during the 2nd-3rd and 4th-5th months of exposure to these pollutants. Several studies have suggested that influenza viruses exhibit higher sensitivity and pathogenicity during the winter season, based on experimental findings [22]. This might also be consistent with the findings of some surveys: The correlation between air pollutants and influenza varies by season and region, with higher effects estimated for the cold season, eastern and central regions, and provinces with wetter conditions and larger populations [23]. In the multi-pollutant lagged interaction model, ambient temperature (AT) is found to be positively correlated with the risk of influenza onset at low temperatures, while pollutants are only weakly positively correlated with the risk of morbidity. This suggests that lower temperatures may facilitate the spread of pollutants, thereby exacerbating the spread of influenza viruses and causing the onset of influenza. These findings are consistent with the seasonal characteristics of influenza onset, which predominantly occur during the winter and spring seasons in this study. The quantitative analysis of interaction models reveals that the interaction between PM_{10} and PM_{2.5} has a more significant effect on influenza onset. Several studies have suggested that particulate matter (PM) may stimulate macrophage apoptosis in lung tissue, which could exacerbate the damage caused by influenza viruses to the respiratory tract. This may explain why the combined effect of PM_{10} and PM_{2.5} is more detrimental to influenza onset [24].

During heavy pollution, reduced outdoor activities and increased indoor activities can heighten the risk of influenza [25]. Hence, the significance of the current study lies in the investigation of the link between pollutants and the development of influenza, which is in line with various domestic and international studies [17, 26, 27].

Temperature and sea level pressure are the most relevant meteorological factors in our study, as they can affect the transmission of pathogens and impact human immune function, leading to respiratory disease development [28,29,30]. In our final GAM model, we did not find a significant association between meteorological factors and influenza incidence. This may be due to the potential confounding effect of regional environment and lifestyle habits on the transmission of influenza.

Our analysis involves multiple models exploring the relationship between environmental factors and influenza incidence, and subsequent subgroup analyses demonstrate significant differences in the predictions made by region and age group. Interestingly, the Gaussian process regression (GPR) model outperforms the Gradient Boosting Regression Tree (GBRT) model in terms of predictive accuracy. In conclusion, the study suggests that Liaoning Province is proficient in predicting influenza outbreaks and accounting for environmental factors. Moreover, it serves as a suitable representative region for the northeastern part of China. Our study also finds that the model's fit and validation are satisfactory for the age group of 25–59 years, who are susceptible to influenza outbreaks. However, the predictive stability is suboptimal for the age group of 5–14 years, possibly due to the clustering of young individuals during the cold season.

It is notable that our study quantifies the impact of various factors on influenza incidence. The use of the GAM model allows us to control important confounding factors and examine the long-term monthly lagged effects of co-exposure to PM_{10} and PM_{2.5}. However, it should be noted that our study is conducted at an aggregate level and does not involve individual-level analysis. And the data is only collected from the Northeast region, so caution is needed when extending our findings to other regions. However, our study provides a starting point for future population epidemiology studies with larger samples and broader geographic coverage.

## Availability of data and materials

The datasets analyzed during the current study are available in the National Public Health Data Centre of China (https://www.phsciencedata.cn/), the China Meteorological Data Sharing Service (www.data.cma.cn) and the National Oceanic and Atmospheric Administration (NOAA**)**.

## References

Nair H, Brooks WA, Katz M, Roca A, Berkley JA, Madhi SA, Simmerman JM, Gordon A, Sato M, Howie S, et al. Global burden of respiratory infections due to seasonal influenza in young children: a systematic review and meta-analysis. Lancet. 2011;378(9807):1917–30.

Moyes J, Cohen C, Pretorius M, Groome M, von Gottberg A, Wolter N, Walaza S, Haffejee S, Chhagan M, Naby F, et al. Epidemiology of respiratory syncytial virus-associated acute lower respiratory tract infection hospitalizations among HIV-infected and HIV-uninfected South African children, 2010–2011. J Infect Dis. 2013;208(Suppl 3):S217-226.

Cohen C, Walaza S, Treurnicht FK, McMorrow M, Madhi SA, McAnerney JM, Tempia S. In- and out-of-hospital mortality associated with seasonal and Pandemic influenza and respiratory syncytial virus in South Africa, 2009–2013. Clin Infect Dis. 2018;66(1):95–103.

Wang Y, Zhang L, Wu SS, Duan W, Sun Y, Zhang M, Zhang XX, Zhang Y, Ma CN, Wang QY, et al. Application of the moving epidemic method in the development of epidemic thresholds and tiered warning alert approachs for influenza prevention in Beijing. Zhonghua Liu Xing Bing Xue Za Zhi. 2020;41(2):201–6.

MacIntyre CR, Cauchemez S, Dwyer DE, Seale H, Cheung P, Browne G, Fasher M, Wood J, Gao Z, Booy R, et al. Face mask use and control of respiratory virus transmission in households. Emerg Infect Dis. 2009;15(2):233–41.

Leuba SI, Yaesoubi R, Antillon M, Cohen T, Zimmer C. Tracking and predicting U.S. influenza activity with a real-time surveillance network. PLoS Comput Biol. 2020;16(11):e1008180.

Sakai T, Suzuki H, Sasaki A, Saito R, Tanabe N, Taniguchi K. Geographic and temporal trends in influenzalike illness, Japan, 1992–1999. Emerg Infect Dis. 2004;10(10):1822–6.

Xu Z, Hu W, Williams G, Clements AC, Kan H, Tong S. Air pollution, temperature and pediatric influenza in Brisbane. Australia Environ Int. 2013;59:384–8.

Wong CM, Yang L, Thach TQ, Chau PY, Chan KP, Thomas GN, Lam TH, Wong TW, Hedley AJ, Peiris JS. Modification by influenza on health effects of air pollution in Hong Kong. Environ Health Perspect. 2009;117(2):248–53.

Rönnqvist M, Ziegler T, von Bonsdorff CH, Maunula L. Detection method for avian influenza viruses in water. Food Environ Virol. 2012;4(1):26–33.

Boone SA, Gerba CP. Significance of fomites in the spread of respiratory and enteric viral disease. Appl Environ Microbiol. 2007;73(6):1687–96.

Otter JA, Donskey C, Yezli S, Douthwaite S, Goldenberg SD, Weber DJ. Transmission of SARS and MERS coronaviruses and influenza virus in healthcare settings: the possible role of dry surface contamination. J Hosp Infect. 2016;92(3):235–50.

Price RHM, Graham C, Ramalingam S. Association between viral seasonality and meteorological factors. Sci Rep. 2019;9(1):929.

Pica N, Bouvier NM. Environmental factors affecting the transmission of respiratory viruses. Curr Opin Virol. 2012;2(1):90–5.

Wang X, Wang Q, Cheng W, Yu Z, Ling F, Mao H, Chen E. Risk factors for avian influenza virus contamination of live poultry markets in Zhejiang, China during the 2015–2016 human influenza season. Sci Rep. 2017;7:42722.

Zhang Y, Ye C, Yu J, Zhu W, Wang Y, Li Z, Xu Z, Cheng J, Wang N, Hao L, et al. The complex associations of climate variability with seasonal influenza A and B virus transmission in subtropical Shanghai. China Sci Total Environ. 2020;701:134607.

Meng Y, Lu Y, Xiang H, Liu S. Short-term effects of ambient air pollution on the incidence of influenza in Wuhan, China: a time-series analysis. Environ Res. 2021;192:110327.

Huang L, Zhou L, Chen J, Chen K, Liu Y, Chen X, Tang F. Acute effects of air pollution on influenza-like illness in Nanjing, China: a population-based study. Chemosphere. 2016;147:180–7.

Suntronwong N, Vichaiwattana P, Klinfueng S, Korkong S, Thongmee T, Vongpunsawad S, Poovorawan Y. Climate factors influence seasonal influenza activity in Bangkok, Thailand. PLoS ONE. 2020;15(9):e0239729.

Kondrich J, Rosenthal M. Influenza in children. Curr Opin Pediatr. 2017;29(3):297–302.

Jaspers I, Ciencewicki JM, Zhang W, Brighton LE, Carson JL, Beck MA, Madden MC. Diesel exhaust enhances influenza virus infections in respiratory epithelial cells. Toxicol Sci. 2005;85(2):990–1002.

Atkinson KV, Bishop LA, Rhodes G, Salez N, McEwan NR, Hegarty MJ, Robey J, Harding N, Wetherell S, Lauder RM, et al. Influenza C in Lancaster, UK, in the winter of 2014–2015. Sci Rep. 2017;7:46578.

Yang J, Yang Z, Qi L, Li M, Liu D, Liu X, Tong S, Sun Q, Feng L, Ou CQ, et al. Influence of air pollution on influenza-like illness in China: a nationwide time-series analysis. EBioMedicine. 2023;87:104421.

Rivas-Santiago CE, Sarkar S, Cantarella Pt, Osornio-Vargas Á, Quintana-Belmares R, Meng Q, Kirn TJ, Ohman Strickland P, Chow JC, Watson JG, et al. Air pollution particulate matter alters antimycobacterial respiratory epithelium innate immunity. Infect Immun. 2015;83(6):2507–17.

Sumpter C, Chandramohan D. Systematic review and meta-analysis of the associations between indoor air pollution and tuberculosis. Trop Med Int Health. 2013;18(1):101–8.

Liu K, Li S, Qian ZM, Dharmage SC, Bloom MS, Heinrich J, Jalaludin B, Markevych I, Morawska L, Knibbs LD, et al. Benefits of influenza vaccination on the associations between ambient air pollution and allergic respiratory diseases in children and adolescents: New insights from the Seven Northeastern Cities study in China. Environ Pollut. 2020;256:113434.

Woodby B, Arnold MM, Valacchi G. SARS-CoV-2 infection, COVID-19 pathogenesis, and exposure to air pollution: What is the connection? Ann N Y Acad Sci. 2021;1486(1):15–38.

Wang W. Progress in the impact of polluted meteorological conditions on the incidence of asthma. J Thorac Dis. 2016;8(1):E57-61.

Ding PH, Wang GS, Guo YL, Chang SC, Wan GH. Urban air pollution and meteorological factors affect emergency department visits of elderly patients with chronic obstructive pulmonary disease in Taiwan. Environ Pollut. 2017;224:751–8.

Fernandes FMC, Martins ES, Pedrosa D, Evangelista M. Relationship between climatic factors and air quality with tuberculosis in the Federal District, Brazil, 2003–2012. Braz J Infect Dis. 2017;21(4):369–75.

## Acknowledgements

Not applicable.

## Funding

No funding.

## Author information

### Authors and Affiliations

### Contributions

Ye Chen: Software, Conceptual, Methodology, Formal analysis, Investigation, Resources, Writing-eview & editing. Weiming Hou: Software, Conceptual, Methodology, Formal analysis, Investigation, Resources, Writing-eview & editing. Weiyu Hou: Methodology, Software, Writing—original draft, Visualization. Jing Dong: Conceptualization, Methodology, Formal analysis, Writing-review & editing, Funding acquisition, Supervision.

### Corresponding author

## Ethics declarations

### Ethics approval and consent to participate

Not applicable.

### Consent for publication

Not applicable.

### Competing interests

The authors declared no competing interests.

## Additional information

### Publisher’s Note

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

## Supplementary Information

## 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.

## About this article

### Cite this article

Chen, Y., Hou, W., Hou, W. *et al.* Lagging effects and prediction of pollutants and their interaction modifiers on influenza in northeastern China.
*BMC Public Health* **23**, 1826 (2023). https://doi.org/10.1186/s12889-023-16712-6

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/s12889-023-16712-6