Skip to main content

NO2 and PM2.5 air pollution co-exposure and temperature effect modification on pre-mature mortality in advanced age: a longitudinal cohort study in China



There is a discourse on whether air pollution mixture or air pollutant components are causally linked to increased mortality. In particular, there is uncertainty on whether the association of NO2 with mortality is independent of fine particulate matter (PM2.5). Furthermore, effect modification by temperature on air pollution-related mortality also needs more evidence.


We used the Chinese Longitudinal Healthy Longevity Study (CLHLS), a prospective cohort with geographical and socio-economic diversity in China. The participants were enrolled in 2008 or 2009 and followed up in 2011-2012, 2014, and 2017-2018. We used remote sensing and ground monitors to measure nitrogen dioxide (NO2), fine particulate matter (PM2.5) , and temperature. We used the Cox-proportional hazards model to examine the association between component and composite air pollution and all-cause mortality, adjusted for demographic characteristics, lifestyle, geographical attributes, and temperature. We used the restricted cubic spline to visualize the concentration–response curve.


Our study included 11 835 individuals with an average age of 86.9 (SD: 11.4) at baseline. Over 55 606 person-years of follow-up, we observed 8 216 mortality events. The average NO2 exposure was 19.1 μg/m3 (SD: 14.1); the average PM2.5 exposure was 52.8 μg/m3 (SD: 15.9). In the single pollutant models, the mortality HRs (95% CI) for 10 μg/m3 increase in annual average NO2 or PM2.5 was 1.114 (1.085, 1.143) and 1.244 (1.221, 1.268), respectively. In the multi-pollutant model co-adjusting for NO2 and PM2.5, the HR for NO2 turned insignificant: 0.978 (0.950, 1.008), but HR for PM2.5 was not altered: 1.252 (1.227, 1.279). PM2.5 and higher mortality association was robust, regardless of NO2. When acccounting for particulate matter, NO2 exposure appeared to be harmful in places of colder climates and higher seasonal temperature variation.


We see a robust relationship of PM2.5 exposure and premature mortality in advance aged individuals, however, NO2 exposure and mortality was only harmful in places of colder climate such as northeast China, indicating evidence of effect modification by temperature. Analysis of NO2 without accounting for its collinearity with PM2.5, may lead to overestimation.

Peer Review reports


Nitrogen dioxide (NO2) has harmful health effects. Epidemiological studies indicate NO2 is associated with bronchitis in asthmatic children and reduced pulmonary function [1]. NO2 is part of the US EPA six Criteria Air Pollutants (along with carbon monoxide, ground-level ozone, particulate matter, sulfur dioxide, and lead) [2] and the WHO Air quality guidelines (along with particulate matter, ozone, sulfur dioxide) [3]. New 2021 WHO Global Air Quality Guidelines (AQGs) recommendations have lower AQG levels based on mortality or cardiovascular mortality studies [4]. The transportation sector mainly drives NO2 exposure. In North America and Europe, urban areas have higher NO2 despite low PM. Because PM2.5 and NO2 tend to be co-exposures, there is no agreement on the causal relationship between NO2 and health, particularly mortality. Particulate matter is identified as a causal agent for total, including cardiovascular and respiratory mortality, while nitrogen oxide is only suggestive as a causal agent for total mortality, but it is identified as an acute trigger of poor respiratory function and asthma in the WHO AQG. In China, the interplay of NO2 and PM2.5, along with temperature effect modification, needs further study.

NO2 is a highly reactive gas known as oxides of nitrogen or nitrogen oxides (NOx). NO2 forms from emissions from automobile exhaust, power plants, and machinery. Exposure to NO2 irritates the airways in the human respiratory system. The health effects of short-term and long-term exposure to NO2 are studied separately. Acute NO2 exposure has been associated with aggravated respiratory diseases, particularly asthma and pulmonary symptoms [5]. Exposure to NO2 may contribute to asthma incidences and other respiratory infections. Several epidemiological studies have linked NO2 to mortality. A meta-analysis of 23 studies found a pooled effect on mortality was 1.04 (95% CI 1.02–1.06) with an increase of 10 μg/m3 in the annual NO2 concentration, independent of the effect of PM2.5 [6]. A more recent review found associations between NO2 and mortality were attenuated upon adjustment for co-pollutants in some studies while not in others [7]. Furthermore, the WHO AQG specified future research needs on air pollution interaction with high and low temperature, or climatic conditions, on health.

Whether NO2 is directly responsible for the health effects or is only an indicator of other pollutants, including particulate matter, evidence from more geographic areas is needed to better understand the concentration–response curve and the generalizability of adverse health effects. First, our study aims to assess the relationship between NO2 and mortality in diverse climatic regions of China. Second, we aim to determine the dose–response relation under the new WHO guidelines from 10 to 40 μg/m3 annual average exposures. Our unique cohort allows us to assess both high and low NO2 exposure regions throughout urban and rural areas of the country. Third, heterogeneous exposures also will enable us to create a multi-pollutant model to assess the collinearity and interaction between NO2 and PM2.5 and how NO2 modifies PM2.5 effects and vice versa. Lastly, we aim to look for effect modification by demographic variables such as age, gender, socioeconomic factors, and temperature, to find the most vulnerable group to NO2 exposure.


Study population

We used the Chinese Longitudinal Healthy Longevity Survey (CLHLS) datasets. It is a longitudinal cohort designed to study healthy longevity. This cohort aims to gather information of the elderly aged 65 and older in 23 provinces of China. The cohort was initially conceived as a survey to study the senior population’s health status, quality of life, socioeconomic characteristics, family, lifestyle, and demographic profile. We overlaid environmental exposure data based on the residential area with remote sensing. Health endpoints include respondents’ health conditions, daily functioning, self-perceptions of health status and quality of life, life satisfaction, mental attitude, and feelings about aging. We used the 2008-2009 wave of Chinese Longitudinal Healthy Longevity Study (CLHLS) with urban and rural coverage in 23 provinces. The participants were enrolled in 2008 or 2009 and followed up to 2018 roughly every two years.

Among the 16 954 participants in the 2008/2009 cohort, we excluded 3109 participants who were lost in the first follow-up or did not have death time records, 267 participants without matched NO2 or PM2.5, 1611 participants with missing values in covariates, and 132 participants aged younger than 65 years. We finally included 11 835 participants.

The CLHLS study was approved by research ethics committees of Peking University (IRB00001052-13074) and Duke University. Written informed consent was obtained from each respondent.

Air pollutant exposure assessment

The concentrations of nitrogen dioxide (NO2) concentration levels (μg/m3) were obtained at an area-level with spatial resolution up to one-kilometer [8]. Land-use regression model corrected for satellite pass time and cloud coverage was directly used for urban areas. For rural areas, NO2 concentrations were adjusted by using surface NO2 concentrations derived from the Ozone Monitoring Instrument satellite NO2 columns. Model performance differed regionally and the  coefficient of determination (R2) was 0.52 in Asia, approximately matched the global average (0.54) [9, 10].

We calculated PM2.5 concentration at an area-level, with baselayer data at 0.01° × 0.01° resolution obtained from the Atmospheric Composition Analysis Group. Exposure assessment techniques utilized monitors at the ground-level for PM2.5 between 1998–2020 (V5.GL.02) by combining Aerosol Optical Depth (AOD) retrievals from the NASA MODIS, MISR, and SeaWIFS instruments with the GEOS-Chem chemical transport model, and subsequently calibrating to global ground-based observations using a Geographically Weighted Regression (GWR), as detailed in the reference [11]. We matched the annual exposure of NO2 and PM2.5 in the year closest to the mortality. We further used the WHO air quality guidelines [4] and the median as the cut-off point to classify NO2 and PM2.5 into different categories.

Mortality outcome assessment

The immediate family members of subjects reported the mortality information in the follow-up surveys. We measured the survival time in months from the first interview to the recorded death date or last interview date.

Covariates measurements

We included baseline characteristics including age, gender, marital status, education, smoking status, drinking status, physical activity, household income, BMI,  residence, geographical region of residence, and temperature. We classified marital status into two categories: Currently married and living with spouse as “married” and widowed/separated/divorced/Never married/married but not living with spouse as “not married.” We used the schooling year to evaluate education level and further classified the schooling year into three groups: 0 years (without formal education), 1–6 years (with primary education), and > 6 years (with higher education). We divided the regular exercise, smoking, and alcohol drinking status into three categories: “Current”, “Former”, and “Never” (See supplementary methods). We also quantified the current alcohol drinker based on the kind of alcohol and how much they drank per day. We defined those who drank equal or less than 14 g pure alcohol per day for the female or 28 g per day for the male as light drinkers, otherwise heavy drinkers (See supplementary methods). There were four categories for the annual household income (yuan): < 4000, < 10,000, < 20,000, and ≥ 20,000. We calculated BMI as body weight divided by the square of the body height (unit: kg/m2). We used the WHO standard of BMI, which defined a BMI of < 18.5 kg/m2 as underweight, a BMI of ≥ 18.5 to < 25 kg/m2 as normal weight, a BMI of ≥ 25 to < 30 kg/m2 as overweight, and a BMI of ≥ 30 kg/m2 as obese. We followed the CLHLS residence categories: “Urban” (including “City” and “Town”) and “Rural.” We divided the geographical region on the basis of residential address to account for climate and dietary differences: central China (Henan, Anhui, Jiangxi, Hubei, Hunan provinces), eastern China (Shandong, Shanghai, Jiangsu, Zhejiang, Fujian provinces), northeastern China (Heilongjiang, Jilin, and Liaoning provinces), northern China (Beijing, Tianjin, Hebei, Shanxi, Shaanxi provinces), southern China (Guangdong, Guangxi, and Hainan provinces), and southwestern China (Chongqing and Sichuan province). Daily meteorological data of the weather monitoring stations across China between 2008 to 2018 was obtained from China Meteorological Administration. Each study participant was matched with meteorological data collected from a monitoring station closest to their area. We used the annual average and standard deviation of the daily mean temperature as the two variables in our analyses.

Statistical analysis

Given the open cohort nature of our cohort with various subjects contributing different person-times to analysis, we decided to use the Cox proportional hazards model to examine the association between long-term NO2 exposure and all-cause mortality. We also calculated the Cumulative Risk Index (CRI) in the two-pollutant model [12]. These models are adjusted for potential confounders or predictors of outcome: age, gender, marital status, education, smoking status, drinking status, physical activity, household income, BMI,  residence, and geographical region of residence. We tried to avoid adjusting for mediators so that we do not reduce the explanatory power of exposure variables, recognizing that some variables are time-varying. To assess for non-linearity, we used the restricted cubic spline to describe the concentration–response relationship. Possible effect modifiers such as age and gender were tested via interaction terms and stratified analyses where needed. We also added the temperature mean and temperature variability (SD) in the same year of the air pollution in the model as a sensitivity analysis. We used R 4.0.0 to run all the analyses.


Those excluded due to the missing of NO2 data had similar age, gender, marriage, education, smoking, and alcohol drinking characteristics. Our study included 11 835 individuals, totaling 55 606 person-years of follow-up. During this time, we counted 8 216 mortality events. This high mortality is expected given the average age of our study participants of 86.9 (SD: 11.4) years old at baseline. Representative of demographic distributions on gender, we had a slightly higher proportion of female participants (57.0%). A more significant proportion of our study participants lived in rural areas (63.6%). Many of the study participants received no formal education, which is typical for the historical period of their birth years. The majority of the study participants were currently not married or living with a spouse at baseline (including having a decreased partner), were never smokers, and never consumed alcohol regularly.

The average exposure of NO2 in the mortality year was 19.1 μg/m3 (SD: 14.1), higher in urban, northern, and eastern regions of China. Participants with higher education or income were also more likely to live in places with higher NO2. There were no large variations of NO2 exposure by different age groups, gender, marriage, exercise, smoking, and alcohol drinking. The average exposure of PM2.5 was 52.8 μg/m3 (SD: 15.9), which was similar between urban and rural areas (52.9 vs. 52.8), higher in northern, central, and southwestern regions of China. Participants with older age, no formal education, not married, not exercising currently, not heavy smokers, with higher BMI tended to live in higher PM2.5 places compared to their counterparts. There was no noticeable difference in PM2.5 exposure for different gender, alcohol drinking status, or household income. Among all the participants, 23.7% (n = 2807) were exposed to NO2 below the WHO recommended AQG level (< 10 μg/m3), and 92.7% (n = 10 967) lived in places that reached the interim target 1 (< 40 μg/m3). However, only 11.9% (n = 1410) of the participants had a PM2.5 exposure lower than the interim target 1 level (< 35 μg/m3) (Table 1).

Table 1 Study population exposure levels relative to 2021 WHO air quality guideline levels

The Pearson correlation coefficient between NO2 (μg/m3) and PM2.5 (μg/m3) was 0.37 (95% CI: 0.35, 0.38). To look for a contrast between NO2 and PM2.5, we looked for concordance and discordance statistics using the 16 μg/m3 for NO2 and 51 μg/m3 for PM2.5 as cut-off points, indicative of median concentrations. There were 17.0% (n = 2007) living in places with high NO2 and low PM2.5, 17.7% (n = 2094) living under low NO2 and high PM2.5, 32.3% (n = 3820) living under low NO2 and low PM2.5, and 33.1% (n = 3914) living under high NO2 and high PM2.5. Additionally, the annual average NO2 and PM2.5 were both significantly negatively associated with the annual average temperature [Pearson coefficient (95%CI): -0.44 (-0.45, -0.42) and -0.30 (-0.32, -0.28) respectively], and significantly positively associated with annual temperature SD [Pearson coefficient (95%CI): 0.44 (0.42, 0.45) and 0.31 (0.30, 0.33) respectively].

As expected in the single pollutant model, higher NO2 was associated with a greater risk for mortality, with the hazard ratio (HR, 95% CI) of 1.114 (1.085, 1.143) for per 10 μg/m3 increase after adjusting for demographics, lifestyles, living regions, BMI, annual average temperature, and annual temperature SD. However, after adjusting for PM2.5, the association between NO2 and mortality was reversed but not significant [HR (95% CI) for per 10 μg/m3 increase: 0.978 (0.950, 1.008). Higher PM2.5 was consistently associated with higher mortality risk [HR (95% CI) for per 10 μg/m3 increase in the single pollutant and two-pollutant model: 1.244 (1.221, 1.268) vs. 1.252 (1.227, 1.279)] (Tables 2 and 3). As we can see, after adjusting for annual average temperature and annual temperature SD in the single pollutant model, the association between NO2 and mortality became stronger while there was no significant change for PM2.5. We calculated the HR (95% CI) for cumulative risk estimates from the two-pollutant model as 1.23824 [1.23823, 1.23825].

Table 2 Association between NO2, PM2.5 concentrations and all-cause mortality by model saturation
Table 3 HR (95% CI) for all-cause mortality considering both NO2 and PM2.5

There was a significant negative interaction between NO2 and PM2.5 (βinteraction = -0.06, Pinteraction < 0.001), and higher NO2 was associated with higher mortality risk only when PM2.5 was lower than 53.3 μg/m3.

The restricted cubic spline for NO2 was supralinear, which means there were larger changes in risk for low concentrations compared with higher concentrations. Meanwhile, the spline also showed a reverse before and after adjusting for PM2.5 (Fig. 1).

Fig. 1
figure 1

Restricted cubic spline describing effect of the NO2 of the year closest to death on mortality before and after adjusting for PM2.5. Note: The left figure adjusted for age, gender, education, household income, marital status, smoking status, drinking status, physical activity, residence, geographical region of residence, and BMI. The right figure additionally adjusted for PM2.5

The association between NO2 and mortality remained significantly positive after adjusting for PM2.5 in some subgroups. It was even stronger than PM2.5 for those living in northeastern China [HR (95% CI: 1.125 (1.035, 1.224) vs. 0.943 (0.858, 1.035)]. NO2 was also still positively associated with mortality after adjusting for PM2.5 in areas with relatively high NO2, low annual temperature mean, or high annual temperature SD [HR (95% CI: 1.047 (1.003, 1.093), 1.083 (1.031, 1.138), and 1.134 (1.08, 1.191) respectively], but negatively associated with mortality in areas with low NO2, high annual temperature mean, or low annual temperature SD [HR (95% CI: 0.749 (0.681, 0.823), 0.843 (0.793, 0.896), and 0.876 (0.833, 0.922) respectively] (Table 4). Perhaps, the NO2 harm effect is more pronounced in colder regions, with mechanisms to be explored in future studies. The spline stratified by exposure level also showed that NO2 tended to have a harmful effect in areas with high NO2 (Figure S2). We found those living under exposure to high NO2 and low PM2.5 were mostly in relatively prosperous regions near the Yangtze Delta (Shanghai, Zhejiang, Jiangsu), while those under the exposure of high NO2 and high PM2.5 were more likely to be in the northern Jing-jin-ji (Beijing, Tianjin, Hebei) areas with a higher concentration of heavy industry.

Table 4 The association between NO2, PM2.5 and mortality in subgroups

Using the four categorical combination terms of NO2 and PM2.5 as the independent variable in cox model: those under high NO2 and low PM2.5 exposure had higher mortality risk than those living in low NO2 and low PM2.5 areas, and those under high NO2 and high PM2.5 exposure also had higher risk than those living in low NO2 and high PM2.5 areas (Table 3s).


In our analysis, we saw that the effect of NO2 disappeared or reversed after adjusting for PM2.5. Several explanations are possible. First, there may be high collinearity of NO2 and PM2.5, and the explanatory power of one variable was overshadowed statistically by the other variable. However, we found this explanation unsatisfactory because the effect of PM2.5 did not change drastically, and only NO2 changed. Second, the phenomenon may manifest the reversal paradox, where the association between an outcome variable and an explanatory (predictor) variable is reversed when another explanatory variable is added to the analysis [13, 14]. This may be a form of Simpson’s Paradox, which refers to a phenomenon whereby the association between a pair of variables (X, Y) reverses sign upon conditioning of a third variable, Z, regardless of the value taken by Z [15]. In our case, the effect estimate direction of NO2 reverses when adjusted for PM2.5. However, in the two pollutant model, the PM2.5 effect was not severely modified by NO2, but the presence of a significant interaction suggests a complex relationship. Other descriptions of this phenomenon are termed Lord’s Paradox or suppression effects.

We attempted to examine this relationship visually in Figure S2. PM2.5 and NO2 may rise and fall together in urban areas, thus making it difficult to separate the effects. To explore possible Simpson’s Bias, we categorized PM2.5 and NO2 by high and low and looked at the HR within each category. We also quantified the number of participants who fall into the WHO recommended AQG levels and interim targets. We can see that the association between PM2.5 and mortality persists regardless of NO2 levels. However, NO2 appears to be protective against mortality at low levels when adjusted for PM2.5. We used a directed acyclic graph (DAG) to explore possible explanations [16] (Figure S3). If we assume there is no causal relationship between NO2 and mortality, and there exists a common pollution source. Conditioning, or adjusting for, PM2.5 should yield a null effect of NO2 on mortality. However, conditioning, or adjusting for, PM2.5 yielded an overall protective effect of NO2 and mortality; we hypothesize there may be a pollution source that produces NO2, but at the same time is an indicator of road traffic access or other indicators beneficial to health.

Although exposure to ambient nitrogen dioxide (NO2) has been linked to increased mortality in several epidemiological studies [6, 7], the question remains whether NO2 is directly responsible for the health effects or is only an indicator of other pollutants, including particulate matter. A systematic review using pooled data from Asia, North America, and Europe found evidence of long-term NO2 on mortality. They found greater similar risk estimates for total mortality of effects of PM2.5 than of NO2 for cardiovascular (20% versus 13%) and respiratory (5% versus 3%) mortality, per 10 μg/m3 of pollutants [6]. Another review got random-effects summary relative risks (RR) ranging from 1.02 to 1.06 for NO2 (per 10 μg/m3) and all-cause (24 cohorts), respiratory (15 cohorts), chronic obstructive pulmonary disease (COPD) (9 cohorts), and acute lower respiratory infection (5 cohorts) mortality. Meanwhile, it identified high levels of heterogeneity for all causes of death except COPD [7]. A time-series analysis using MCC data in 398 cities in 22 countries or regions found an association of NO2 and total cardiovascular and respiratory mortality [17]. The study used death records using ICD-9 or -10 codes. This study found the pooled concentration–response curves for all three causes were almost linear without discernible thresholds. A study in northern China also found higher NO2 was associated with lower all-cause mortality risk no matter adjusting or not adjusting for PM10 or SO2, and it only showed a harmful effect on lung cancer mortality when adjusting for PM10 [18]. Another study based on Dutch national databases found the positive association between NO2 and mortality remained for non-accidental and lung cancer mortality, but reversed for circulatory diseases mortality and disappeared for respiratory diseases mortality after adjusting for PM10 [19].

Our study has many strengths. First, we utilized a prospective cohort originally designed to ascertain determinants of healthy longevity, and thus our study benefited from having access to a wide range of confounders for adjustment. Our study’s sample size and various regions allowed us to see a wide spectrum of air pollution exposure levels, or having heterogenous exposures. Our study contains some limitations typical of observational epidemiologic studies and also specific to our study design. First, our exposure ascertainment of air pollutants relied in part on remote sensing modeling techniques, and we do not have a personalized air pollution monitor for each individual. Some people may live in households with biomass for cooking and heating and could suffer from high indoor air pollutants. Second, we cannot estimate cause-specific mortality because death was ascertained from next-of-kin, who could not report clinically accurate mortality causes. Third, there is potential unaccounted for residual confoundings, such as underlying social-economic statuses that lead to differential air pollution exposure and are related to the health outcome. Lastly, a possible exposure misclassification may arise NO2 may be a regional pollutant and could vary substantially in space, a difficulty for accurate exposure assessment. As the elderly population may spend most of their lives indoors, NO2’s penetration coefficients from outdoor to indoor are lower than PM2.5’s [20], which may have impacted our findings.

In the 2021 World Health Organization Air Quality Guidelines, the annual guideline for nitrogen dioxide (NO2) is four times tighter than the 2005 limit value and is down to 10 μg / m3, from 40 μg / m3. The 24-h guideline of 25 μg / m3 has been additionally introduced. PM2.5 threshold has been limited to a mere five μg / m3, from 10 µg per cubic meter (μg / m3). The 24-h average is 15 μg / m3, from 25 μg / m3. Both annual and 24-h average guidelines for PM10 are lower by five μg / m3 each. The US EPA integrated Science Assessment advises that NO2 is a suggested but not causal factor with mortality. The WHO global air quality guideline also is inconclusive on the causality of NO2 and mortality. Our study does not indicate an association of NO2 on mortality independent of PM2.5. It is possible that NO2 does not have a strong relationship with mortality at low levels, as there is evidence that the health effect of NO2 is on respiratory health only. Perhaps prior observed associations between ill health and NO2 at low concentrations in the ambient air result from co-exposure by particulate matter [21].


Our findings indicate a consistent harmful effect of PM2.5 on all-cause mortality in a cohort of advanced aging population in China. We do not see harmful effects of NO2 when adjusted for PM2.5 on all-cause mortality. The results of our analysis suggest a complex interplay of these two air pollutants. We see consistent harmful effects of PM2.5 with all-cause pre-mature mortality among a cohort of elderly individuals. But NO2 was only harmful in some subgroups, namely colder regions. There are atmospheric explanations, such as the transfer of NO2 to PM2.5, exposure assessment situation where there is measurement error of the air pollution exposure, or that NO2 is more of a proxy for commercial activity, which in turn leads to better health, at least under the developing country context. Alternatively, it is possible that the NO2 does indeed have a null effect independent of PM2.5. Future studies of multiple pollutant models, along with temperature effect modification, are needed to determine the causal mechanisms for air pollution mixtures and health.

Availability of data and materials

The CLHLS datasets are available from the Peking University Open Research Data ( and Inter-university Consortium at University of Michigan (


NO2 :

Nitrogen dioxide

PM2.5 :

Fine particulate matter


Chinese Longitudinal Healthy Longevity Study


Cumulative Risk Index


Hazard ratio


Directed acyclic graph


Relative risks


Chronic obstructive pulmonary disease


  1. Mortimer K, Neugebauer R, Lurmann F, Alcorn S, Balmes J, Tager I. Air pollution and pulmonary function in asthmatic children: effects of prenatal and lifetime exposures. Epidemiology. 2008;19:550–2.

    Article  Google Scholar 

  2. United States Environmental Protection Agency. Criteria Air Pollutants. 2021. 03 Oct 2022.

  3. World Health Organization. WHO Air quality guidelines for particulate matter, ozone, nitrogen dioxide and sulfur dioxide. 2005.

  4. World Health Organization. WHO global air quality guidelines. Particulate matter (PM2.5 and PM10), ozone, nitrogen dioxide, sulfur dioxide and carbon monoxide. Geneva; 2021. 03 Oct 2022.

  5. Hesterberg TW, Bunn WB, McClellan RO, Hamade AK, Long CM, Valberg PA. Critical review of the human data on short-term nitrogen dioxide (NO2) exposures: evidence for NO2 no-effect levels. Crit Rev Toxicol. 2009;39:743–81.

    Article  CAS  Google Scholar 

  6. Faustini A, Rapp R, Forastiere F. Nitrogen dioxide and mortality: Review and meta-analysis of long-term studies. Eur Respir J. 2014;44:744–53.

    Article  CAS  Google Scholar 

  7. Huangfu P, Atkinson R. Long-term exposure to NO2 and O3 and all-cause and respiratory mortality: A systematic review and meta-analysis. Environ Int. 2020;144:105998.

  8. Mohegh A, Anenberg S. Global surface NO2 concentrations 1990-2020. figshare. 2020. Accessed 03 Oct 2022.

  9. Larkin A, Geddes JA, Martin RV, Xiao Q, Liu Y, Marshall JD, et al. Global Land Use Regression Model for Nitrogen Dioxide Air Pollution. Environ Sci Technol. 2017;51:6957–64.

    Article  CAS  Google Scholar 

  10. Anenberg SC, Mohegh A, Goldberg DL, Kerr GH, Brauer M, Burkart K, et al. Long-term trends in urban NO2 concentrations and associated paediatric asthma incidence: estimates from global datasets. Lancet Planet Heal. 2022;6:e49-58.

    Article  Google Scholar 

  11. van Donkelaar A, Hammer MS, Bindle L, Brauer M, Brook JR, Garay MJ, et al. Monthly Global Estimates of Fine Particulate Matter and Their Uncertainty. Environ Sci Technol. 2021;55:15287–300.

    Article  CAS  Google Scholar 

  12. Crouse DL, Peters PA, Hystad P, Brook JR, van Donkelaar A, Martin RV, et al. Ambient PM2.5, O3, and NO2 exposures and associations with mortality over 16 years of follow-up in the canadian census health and environment cohort (CanCHEC). Environ Health Perspect. 2015;123:1180–6.

  13. Nickerson CA, Brown NJL. Simpson’s Paradox is suppression, but Lord’s Paradox is neither: Clarification of and correction to Tu, Gunnell, and Gilthorpe (2008). Emerg Themes Epidemiol. 2019;16:1–11.

    Article  Google Scholar 

  14. Tu Y-K, Gunnell D, Gilthorpe MS. Simpson’s Paradox, Lord’s Paradox, and Suppression Effects are the same phenomenon – the reversal paradox. Emerg Themes Epidemiol. 2008;5:2.

    Article  Google Scholar 

  15. Comment PJ. Understanding Simpson’s Paradox. Am Stat. 2014;68:8–13.

    Article  Google Scholar 

  16. Weisskopf MG, Seals RM, Webster TF. Bias amplification in epidemiologic analysis of exposure to mixtures. Environ Health Perspect. 2018;126:047003.

  17. Meng X, Liu C, Chen R, Sera F, Vicedo-Cabrera AM, Milojevic A, et al. Short term associations of ambient nitrogen dioxide with daily total, cardiovascular, and respiratory mortality: Multilocation analysis in 398 cities. BMJ. 2021;372:n534.

  18. Chen X, Zhang L wen, Huang J ju, Song F ju, Zhang L ping, Qian Z min, et al. Long-term exposure to urban air pollution and lung cancer mortality: A 12-year cohort study in Northern China. Sci Total Environ. 2016;571:855–61.

  19. Fischer PH, Marra M, Ameling CB, Hoek G, Beelen R, De Hoogh K, et al. Air pollution and mortality in seven million adults: The dutch environmental longitudinal study (DUELS). Environ Health Perspect. 2015;123:697–704.

    Article  CAS  Google Scholar 

  20. Hu Y, Yao M, Liu Y, Zhao B. Personal exposure to ambient PM2.5, PM10, O3, NO2, and SO2 for different populations in 31 Chinese provinces. Environ Int. 2020;144:106018. doi:

  21. Seaton A, Dennekamp M. Hypothesis: III health associated with low concentrations of nitrogen dioxide - An effect of ultrafine particles? Thorax. 2003;58:1012–5.

    Article  CAS  Google Scholar 

Download references


The authors thank all the participants and data collection staff of the CLHLS study.


The Chinese Longitudinal Healthy Longevity Study (CLHLS) datasets analyzed in this paper are jointly supported by the National Key R&D Program of China (2018YFC2000400), National Natural Sciences Foundation of China (72061137004,71490732), and the U.S. National Institute of Aging of National Institute of Health (P01AG031719). The funders had no role in this study analysis, interpretation of data, or writing the manuscript.

Author information

Authors and Affiliations



J.S.J. and LX.L. conceptualized the study, conducted statistical analysis, drafted and edited the article; Y.Z. supervised the CLHLS data collections; JF.Z., HD.K., B.Z., Y.Z., and K.G.B. interpreted the results and revised the article. All authors provided critical insights and reviewed the article. The author(s) read and approved the final manuscript. Y.Z. is the senior authors of the study. 

Corresponding authors

Correspondence to John S. Ji or Yi Zeng.

Ethics declarations

Ethics approval and consent to participate

The research ethics committees of Peking University (IRB00001052-13074) and Duke University approved the study. All participants in the study have given informed consents.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have 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

Additional file 1: Supplementary methods.

Table S1. Population Characteristics by High and Low PM2.5 and NO2 Exposure Levels (unit: μg/m3). Table S2. Association between NO2, PM2.5 levels and all-cause mortality adjusting for different variables. Table S3. The association between NO2, PM2.5 and mortality stratified by different exposure level groups (unit: μg/m3). Figure S1. The scatter plot of NO2 and PM2.5 of the year closest to outcome assessment. Figure S2. The restricted cubic spline of NO2 and PM2.5 on mortality. Figure S3. DAGs of NO2 and PM2.5 Relationship on Mortality.

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 The Creative Commons Public Domain Dedication waiver ( 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

Ji, J.S., Liu, L., Zhang, J.(. et al. NO2 and PM2.5 air pollution co-exposure and temperature effect modification on pre-mature mortality in advanced age: a longitudinal cohort study in China. Environ Health 21, 97 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: