Short-term exposure to ozone and nonaccidental mortality in Northeast China＊

ZHANG Yi Fan, MA Yu Xia,#, LI He Ping, CHENG Bo Wen, FENG Feng Liu, and JIAO Hao Ran,2

Ambient pollutants are a major environmental risk factor for human health[1]. Ozone (O3) is a highly reactive, oxidative gas that is formed through the interaction of volatile organic compound precursors and nitrogen oxides (NOx). Several studies have demonstrated that exposure to ambient O3in the troposphere is associated with increased mortality. A study in the United States found that an increase in nonaccidental mortality among all age groups was associated with increased O3concentrations[2].Research in some European countries has also shown that short-term O3exposure contributes to nonaccidental death[3]. Related studies have been conducted in Asian countries. A multicenter time series study conducted in 21 major cities in East Asia provided strong evidence of a significant association between short-term O3exposure and mortality[4]. In China, due to the increase in motor vehicle ownership, O3pollution has become a major public health concern. However, few epidemiological studies on the health effects of O3in Chinese cities have been conducted due to a lack of long-term observational data. Given that the aforementioned studies were only conducted in large cities, the generalizability of the results is unclear. Therefore, in this study, a time series analysis was carried out in Harbin, a city in northeast China, to assess the effect of O3on mortality using data from 2014 to 2016. The findings of this study may be useful as reference for policies that mitigate the negative health effects of O3in northeast China.

This study included 12 districts in Harbin(Supplementary Figure S1 available in www.besjournal.com). Daily nonaccidental mortality data were acquired from Heilongjiang Provincial Center for Disease Control and Prevention. The causes of death were classified on the basis of theInternational Classification of Diseases, 10th edition(ICD-10). Only nonaccidental diseases (ICD-10: A00-R99) were included in this study. Air pollutant data were obtained from the Harbin Environmental Monitoring Center. The data collected by the Harbin Environmental Monitoring Center were retrieved from air quality monitoring stations located in Harbin. Less than 3.5% of the data from the monitoring stations were missing. Daily pollutant data were the average of the observations from all stations in Harbin. Missing data were replaced with the average of the observations from the remaining monitoring stations on the same day. The O3concentration used in this study was the daily maximum of the 8-hour running average concentration, and for other pollutants, 24-hour average concentrations were used. Daily meteorological data during the study period were collected from Harbin Meteorological Bureau.

The association between daily nonaccidental mortality and O3concentration during 2014-2016 was evaluated using the following model (Equation 1):

whereE(Yi) refers to the predicted number of deaths on dayi;αis the intercept term;HolidayandDOWdenote adjustments for public holidays and the day of the week, respectively;βis the regression coefficient;Xirepresents the O3concentration on dayi;s(time,df) ands(Zi,df) refer to the spline function of calendar time and meteorological factors, respectively; anddfin the spline function represents the degree of freedom. We used the Akaike information criterion to determine the df corresponding to the optimal model fitting.Considering the potential lagged effects of pollutants on mortality, both single-day lag (lag 0 to lag 6: lag 0 is the exposure of the current day and lag 6 is the exposure from 6 days prior) and multiday lag (lag 01 to lag 06: the moving average of the current and prior 1-6 days) structures were analyzed. We also fitted multipollutant models to assess the collinearity between pollutants. After identifying the optimal number of days of lag, the effect of O3on mortality was calculated by adding other pollutants to the model using natural cubic splines.Furthermore, the data were stratified by sex (female and male) and season (warm season: April-September, and cold season: October-March) for subgroup analyses. Two-sampleZtests were used to determine whether subgroup differences were statistically significant. All analyses were conducted using R version 3.6.1 (R Foundation for Statistical Computing, Vienna, Austria) with the mgcv package.Results are presented as relative risks (RRs) per interquartile range (IQR) increase in O3concentration. All statistical tests were two-sided,andPvalues of < 0.05 were considered statistically significant.

From 2014 to 2016 (1,096 days), 142,251 nonaccidental deaths (58,711 for females and 83,540 for males) were recorded in the study area.Supplementary Figure S2A (available in www.besjournal.com) shows the spatial distribution of mortality. During the study period, the mean daily number of nonaccidental deaths was 129.79 with a standard deviation of 0.53, and the mean daily maximum 8-hour concentration of O3was 64.99 μg/m3with a standard deviation of 0.93(Supplementary Table S1 available in www.besjournal.com), which is less than the national grade II standard (160 μg/m3); however, the daily maximum was higher than the air quality standard of the World Health Organization (100 μg/m3) on 13%of the days in the study period. Time series analysis revealed a clear seasonal pattern in daily O3concentrations and mortality. The concentration of O3was higher in warm seasons than in cold seasons,whereas the daily death count was slightly lower in warm seasons than in cold seasons (Supplementary Figure S2B).

The exposure-response relationship between O3concentration on the current day and mortality is illustrated in Supplementary Figure S3 (available in www.besjournal.com). The curve indicates that mortality risk was increased with increasing O3concentration. Given the likely continuous effects of ambient pollutants on human health, the association between O3exposure and nonaccidental mortality on different lag days is demonstrated in Figure 1.Results indicate that nonaccidental mortality significantly increased by 2.0% [95% confidence interval (CI): 0.6%-3.4%, at lag 1] for an IQR increment in O3concentration. Previous studies have investigated the health effects of O3. For example, in the early years, the APHENA (Air Pollution and Health: A Combined European and North American Approach) project in Europe and North America showed that for every 10 μg/m3increment in O3concentration, the increase in nonaccidental mortality among all age groups was 0.17% (95%CI:0.09%-0.25%) in Europe, 0.56% (95%CI:0.28%-0.83%) in Canada, and 0.18% (95%CI:0-0.35%) in the United States[2]. In Shanghai, China,Kan et al. found a significant correlation between O3and daily nonaccidental mortality and reported a 0.31% (95%CI: 0.04%-0.58%) increase in mortality due to O3exposure[5]. Similar findings have been reported in other multicity studies and metaanalyses[6]. However, unlike the aforementioned findings, a few studies have indicated that no significant association exists between O3exposure and mortality, and some studies have even found a significant negative association. This discrepancy may be explained by the fact that O3is a highly reactive molecule that is easily consumed when inhaled; thus, its effects on the cardiovascular system remain unclear. Furthermore, studies have shown that oxidative stress may be beneficial at low doses and harmful at moderate to high doses[7]; this leads to uncertainty about the health effects of O3.Additionally, differences between study populations can influence the health burden associated with O3.Nonbiological factors, including climate, geographic location, availability and quality of health care,lifestyle, and differences in study methodology also contribute to discrepancies in the findings of different studies.

Figure 1. RRs and 95% CIs for nonaccidental mortality with each IQR increase in O3 concentration at single-day lags. IQR: per interquartile range.

The possible biological mechanisms through which O3contributes to mortality include oxidative stress, airway inflammation, and the triggering of autonomic reflexes that influence the regulation of heart rate and vascular tone. These changes in the human body may damage respiratory tissues and lung function, alter the autonomic control of cardiac frequency, increase blood pressure, and ultimately lead to fatal outcomes[8].

When stratified by sex, the strongest and most statistically significant association between O3exposure and nonaccidental mortality was found at lag 0 for females and lag 1 for males; mortality increased by 2.6% (95%CI: 0.3%-5.0%) for females and 1.9% (95%CI: 0.2%-3.7%) for males. O3had a slightly greater effect on females than on males,although the differences between subgroups were insignificant (Ztest:P＞ 0.05). This may be due to differences in physiological characteristics between females and males. Females have fewer red blood cells and a narrower respiratory tract than males;thus, females are more susceptible to pollutants and have higher airway reactivity than males[9]. When the study period was stratified into cold and warm seasons, the effect of O3on nonaccidental mortality was generally stronger and more urgent during warm seasons compared with the effect during cold seasons, although the difference was only statistically significant for males (Ztest: total:P＞0.05, female:P＞ 0.05, male:P< 0.05). This may partly be explained by climatic features and human activities. During warm seasons, in northern China,O3concentrations are typically higher, people spend more time outdoors, and windows tend to be open for longer, which may increase O3exposure[10].During cold seasons, however, people go out less frequently and windows are usually closed, which may decrease O3exposure.

Table 1 presents the results obtained frommultiday lag structures. Compared with the estimates calculated from single-day lags(Supplementary Table S2 available in www.besjournal.com),RRvalues were generally higher under the multiday lag structures. These results indicate that the risk of death from O3exposure accumulates over time. Because other pollutants in the environment might influence the effects of O3,we added NO2, SO2, CO, PM2.5(particulate matter< 2.5 μm diameter), and PM10(PM < 10 μm diameter) into the models to perform multipollutant analysis.

Table 1. RRs and 95% CIs for nonaccidental mortality with each IQR increase in O3 concentration at multiday lags

As shown in Table 2, almost allRRvalues of mortality changed marginally after adjusting for other pollutants, and the estimated effects remained statistically significant, indicating that the results were robust.

Table 2. Estimated effects of O3 exposure on nonaccidental mortality after adjusting for NO2, SO2, CO, PM2.5,and PM10

The present study has certain limitations. First,as with most time series studies, we used pollutant concentrations from environmental monitoring stations that do not reflect true exposure at the individual level. Second, due to limited data availability, confounding factors such as socioeconomic status and individual behaviors were not controlled for; such factors may mediate the health effects of O3. We were also unable to conduct a stratified analysis according to age due to lack of data. Finally, our study was limited to Harbin; the results may not be generalizable to other regions.Thus, caution should be exercised when extrapolating the results of this study.

The authors have no conflicts of interest to declare.

&These authors contributed equally to this work.

#Correspondence should be addressed to MA Yu Xia,Tel: 81-931-8915728, E-mail: mayuxia07@lzu.edu.cn

Biographical notes of the first authors: ZHANG Yi Fan,female, born in 1996, PHD Candidate, majoring in the health effects of ambient pollutants; MA Yu Xia, female,born in 1974, Professor, majoring in weather and climate effect on health.

Received: July 12, 2022;

Accepted: November 4, 2022