华中某市大气污染物O3、PM2.5暴露对居民死亡的时间序列研究

彭星宇, 王彦丁, 张新民, 谭海萍, 郭庶

彭星宇, 王彦丁, 张新民, 谭海萍, 郭庶. 华中某市大气污染物O3、PM2.5暴露对居民死亡的时间序列研究[J]. 环境与职业医学, 2023, 40(3): 331-341. DOI: 10.11836/JEOM22254
引用本文: 彭星宇, 王彦丁, 张新民, 谭海萍, 郭庶. 华中某市大气污染物O3、PM2.5暴露对居民死亡的时间序列研究[J]. 环境与职业医学, 2023, 40(3): 331-341. DOI: 10.11836/JEOM22254
PENG Xingyu, WANG Yanding, ZHANG Xinmin, TAN Haiping, GUO Shu. Associations of ambient PM2.5 and O3 with human mortality: A time-series study in a city of central China[J]. Journal of Environmental and Occupational Medicine, 2023, 40(3): 331-341. DOI: 10.11836/JEOM22254
Citation: PENG Xingyu, WANG Yanding, ZHANG Xinmin, TAN Haiping, GUO Shu. Associations of ambient PM2.5 and O3 with human mortality: A time-series study in a city of central China[J]. Journal of Environmental and Occupational Medicine, 2023, 40(3): 331-341. DOI: 10.11836/JEOM22254

华中某市大气污染物O3、PM2.5暴露对居民死亡的时间序列研究

基金项目: 国家重点研发计划课题(2020YFC1806303、2019YFC1803401)
详细信息
    作者简介:

    彭星宇(1997—),男,硕士生;E-mail:1552511610@qq.com

    通讯作者:

    郭庶,E-mail:guoshu@scies.org

  • 中图分类号: R122

Associations of ambient PM2.5 and O3 with human mortality: A time-series study in a city of central China

Funds: This study was funded.
More Information
  • 摘要:
    背景

    大气污染物O3、PM2.5暴露与人群死亡存在密切关联。研究往往集中在沿海地区,而中西部地区研究较少。

    目的

    探讨华中某市大气O3、PM2.5对居民死亡的急性影响。

    方法

    收集华中某市2015年1月1日—2021年6月30日的大气污染物数据、气象数据以及人群死亡资料。气象数据包括日均气温、气压、相对湿度;大气污染物数据包括PM2.5、PM10、SO2、NO2、CO的日均质量浓度(后简称为浓度)以及O3的日最大8 h平均浓度。采用基于泊松分布的广义相加模型定量探讨大气污染物(O3、PM2.5)暴露与居民死亡之间的关系,并进行年龄、性别和文化程度的分层分析。

    结果

    研究期间该市O3的日最大8 h平均浓度为94.38 μɡ·m−3,PM2.5的日均浓度为55.56 μɡ·m−3。在单污染物模型中,O3与总死亡以及呼吸系统、循环系统、高血压、冠心病、中风死亡的关联分别在lag02、lag2、lag02、lag0、lag02、lag0时最高,O3每升高10 μɡ·m−3,其ER(95%CI)分别为0.09%(−0.08%~0.25%)、0.35%(0~0.71%)、0.43%(0.18%~0.68%)、0.45%(0.02%~0.91%)、0.59%(0.16%~1.02%)、0.33%(0.01%~0.65%),其中O3对总死亡的影响没有统计学意义(P>0.05)。PM2.5与总死亡以及呼吸系统、循环系统、高血压、冠心病、中风死亡的关联分别在lag1、lag5、lag01、lag05、lag04、lag01时最高,PM2.5每升高10 μɡ·m−3,其ER(95%CI)分别为0.02%(−0.09~0.13%)、0.25%(0.01%~0.50%)、0.35%(95% CI:0.16%~0.54%)、1.18%(95% CI:0.59%~1.77%)、0.17%(−0.13%~0.40%)、0.65%(0.38%~0.92%),PM2.5对总死亡和冠心病死亡的影响没有统计学意义(P>0.05)。在暖季(5—10月)时,O3每升高10 μɡ·m−3,男性、6~65岁人群、>65岁人群、高中以下文化人群的总死亡的ER(95%CI)分别为0.46%(0.16%~0.75%)、0.38%(0.08%~0.68%)、0.41%(0.14%~0.66%)、0.38%(0.14%~0.61%),而在冷季(11月—次年4月)时,O3效应没有统计学意义(P>0.05)。

    结论

    大气污染物(O3、PM2.5)对该市居民死亡存在急性影响,老年人、高中文化以下的人群以及患有循环系统疾病的人群对O3、PM2.5暴露更为敏感。

     

    Abstract:
    Background

    Exposure to air pollutants O3 and PM2.5 is closely related to population mortality. Most of the domestic research findings are for residents in coastal areas, and less for those in the central and western regions.

    Objective

    To investigate the acute effects of O3 and PM2.5 on the mortality of residents in a city of central China.

    Methods

    Data were collected on atmospheric pollutants, meteorological data, and population mortality in a city of central China from January 1, 2015 to June 30, 2021. Meteorological data included daily average temperature, air pressure, and relative humidity. Atmospheric pollution data included daily mean concentrations of PM2.5, PM10, SO2, NO2, and CO and maximum 8 h O3. Generalized additive model with Poisson distribution was used for estimating the relationships between air pollutants (O3 and PM2.5) and population mortality, and further stratified by age, gender, and education.

    Results

    The daily maximum 8 h average concentration of O3 in the city during the study period was 94.38 μɡ·m−3 and the daily average concentration of PM2.5 was 55.56 μɡ·m−3. In the single-pollutant model, the correlations between O3 concentration and total deaths as well as deaths due to respiratory, circulatory, hypertension, coronary heart disease, and stroke were strongest at lag02, lag2, lag02, lag0, lag02, and lag0, and for every 10 μɡ·m−3 increase in concentration of O3, the associated ER (95%CI) values of daily mortality were increased by 0.09% (−0.08%–0.25%), 0.35% (0–0.71%), 0.43% (0.18%–0.68%), 0.45% (0.02%–0.91%), 0.59% (0.16%–1.02%), and 0.33% (0.01%–0.65%), respectively. The effect of O3 on total mortality was not statistically significant (P>0.05). The correlations between PM2.5 concentration and total deaths, as well as deaths due to respiratory, circulatory, hypertension, coronary heart disease, and stroke were strongest at lag1, lag5, lag01, lag05, lag04, and lag01, and for every 10 μɡ·m−3 increase in concentration of PM2.5, the associated ER (95%CI) values of daily mortality increased by 0.02% (−0.09–0.13%), 0.25% (0.01%–0.50%), 0.35% (0.16%–0.54%), 1.18% (0.59%–1.77%), 0.17% (−0.13%–0.40%), and 0.65% (0.38%–0.92%), respectively, with no statistically significant effects of PM2.5 on total mortality and mortality due to coronary heart disease (P>0.05). During warm season (from May to October), the ER (95%CI) values of total deaths per 10 μɡ·m−3 increase in O3 in male, people aged 6~65 years, people aged >65 years, and people below high school education were 0.46% (0.16%–0.75%), 0.38% (0.08%–0.68%), 0.41% (0.14%–0.66%), and 0.38% (0.14%–0.61%), respectively, while the O3 effect was not statistically significant (P>0.05) during cool season (from November to April).

    Conclusions

    Atmospheric pollutants (O3 and PM2.5) have acute effects on mortality in the city, with the elderly, people with less than a high school education, and those with circulatory disease being more sensitive to O3 and PM2.5 exposures.

     

  • 快速的城市化、工业化是目前全球气温剧变的主要促成因素,并且导致了严重的空气污染问题,威胁着人类的生命健康。世界卫生组织2022年4月发布的《WHO空气质量数据库(2022年更新)》显示[1],全球约99%人口呼吸的空气超过了世卫组织空气质量指南中的浓度限值。臭氧(ozone, O3)和细颗粒物(particles with diameter<2.5µm, PM2.5)是最常见的两种大气污染物,可导致人群过早死亡[2-3]。流行病学研究发现O3短期暴露对人体健康会产生不良影响,包括人群死亡率升高、呼吸系统疾病住院率和急诊科就诊率增加、以及肺功能下降[4-6]。PM2.5因为其对人体的危害大,世界卫生组织将PM2.5污染视为全球卫生面临的十项威胁之一[7]。据报道,PM2.5易沉积在肺泡区和肺的终末细支气管区[8],还可进入血液中,导致血栓形成的风险升高,甚至影响肺功能[9],除此,PM2.5暴露还与急性下呼吸道疾病[10-11]、肺癌[12]以及慢性阻塞性肺疾病[13-14]患病有关。

    由于污染物数据的缺乏,目前国内的研究大多聚焦于长三角、珠三角、京津冀等经济发达地区,而关于中西部地区研究的较少[15]。Yin等[5]、Liu等[16]的多城市研究结果表明,各个地区的污染物暴露与死亡率之间的暴露-反应系数都不同。研究发现,与其他臭氧浓度度量方式相比,日最大8 h平均臭氧浓度与人群死亡率之间有着更密切的联系[17]。因此,本研究将使用臭氧日最大8 h平均值,基于广义相加模型(generalized additive model, GAM),探讨华中某市大气污染物O3、PM2.5浓度与居民死亡之间的关系,并分析年龄和性别的影响。

    从该市疾病预防控制中心收集2015年1月1日至2021年6月30日居民死因数据资料,数据指标主要包括年龄、性别、文化程度、死因,对其具体死因按照第十版《国际疾病分类》(International Classification of Disease, ICD)分为总死亡(I00~Y99)、非意外死亡(I00~R99)、循环系统疾病(I00~I99)、呼吸系统疾病(J00~J99),以及特定疾病高血压(I00~I15)、冠心病(I20~I25)、中风(I60~I69)。本研究使用的数据已通过中国疾病预防控制中心环境与健康相关产品安全所伦理审查委员会审查(编号:202124)。

    大气污染资料来源于该市环保局5个环境监测站点日报数据,取均值作为全市居民的暴露水平。污染数据主要有O3日最大8 h平均质量浓度(后简称为浓度)和PM2.5、可吸入颗粒物(PM10)、二氧化硫(SO2)、二氧化氮(NO2)、一氧化碳(CO)的日均浓度。同期气象数据来源于该市气象局,数据资料主要包括日平均气温、湿度以及气压。

    对所收集到的人群死亡数据、大气污染数据以及气象数据采用$ \bar x \pm s $、最大值(Max)、最小值(Min)、百分位数(P25P50P75)进行统计描述。使用Spearman非参数相关分析大气污染物和气象指标之间的相关性。

    本研究使用基于泊松分布的GAM分析O3、PM2.5浓度与居民日死亡人数的相关性。同时在模型中控制时间趋势、气象因素(气温、湿度、气压)、星期几效应等混杂因素。采用平滑样条函数拟合时间序列的长期趋势,气象指标如日均气温、日均相对湿度、日均气压也均使用平滑样条函数拟合。星期几效应以分类变量纳入模型中。GAM模型公式如下:

    $$ \begin{split} &\mathrm{log}[E\left({Y}_{t}\right)]={\text{α}} +\beta {X}_{t}+\\ &{\rm{ns}}\left({\rm{time}},{V}_{df}\right)+ns\left(temperature,{V}_{df}\right)+\\ &{\rm{ns}}\left({\rm{RH}},{V}_{df}\right)+s\left(pressure,{\rm{df}}\right)+{V}_{DOW} \end{split}$$ (1)

    式(1)中:$ {Y}_{t} $为观察日$ t $日当天的实际死亡人数;$ E({Y}_{t} $)为$ t $日死亡人数的期望值;ns为自然样条平滑函数;$ {X}_{t} $为观察日$ t $日的污染物浓度;表示暴露-反应系数;$ temperature、RH、pressure $分别表示气温、相对湿度、气压;$ {V}_{df} $表示自由度。通过查阅文献,最终时间趋势的自由度选择7,温度的自由度选择6[16],湿度和气压的自由度选择3[18]

    首先建立单污染模型,在控制气象因素、星期几效应等混杂因素后,将大气污染物O3、PM2.5纳入模型中,分析探讨PM2.5、O3浓度与居民当日(lag0)死亡人数的相关性。通过查阅文献,发现臭氧暴露具有滞后效应,因此选取lag1~lag5(lag1是指前1天污染物平均浓度,以此类推)进行分析[19];考虑到单独的滞后效应可能会降低O3与居民每日死亡人数之间的相关性,故又选取lag01~lag05(lag01表示污染物在当天至前1天的平均浓度,以此类推)进行累积滞后效应分析。通过单独滞后效应和累计滞后效应分析,可以确定O3(或PM2.5)暴露对居民死亡率影响最大的滞后时间。最后按年龄、性别、文化程度进行O3与居民死亡人数的分层分析。敏感性分析:在确定最佳滞后时间后,建立双污染模型,选取单污染模型中的最佳滞后时间,纳入其他污染物(PM10、SO2、NO2、CO),观察臭氧与居民死亡人数之间关系的稳定性。

    以O3(或PM2.5)每升高10 μɡ·m−3,居民每日死亡人数的超额危险度(excess risk , ER)及其95%置信区间(95%CI)来表示影响程度。ER的计算公式如下所示:

    $$ \mathrm{E}\mathrm{R}=\left[\mathrm{E}\mathrm{X}\mathrm{P}\left(\beta \times 10\right)-1\right]\times 100\% $$ (2)

    运用R 4.1.2软件进行时间序列分析,研究污染物浓度每升高10 μɡ·m−3对人群死亡率的影响,并进行年龄(6~65岁、>65岁)、性别(男性、女性)、文化程度(高中以下文化、高中及以上文化)分层分析。检验水准为α=0.05。

    2015—2021年研究地区居民死亡数据的基本情况信息见表1。在研究期间内,该市总死亡117001人,平均每日总死亡人数为49.39人,男性死亡占总死亡人数的56.17%,>65岁以上老人占总死亡人数的74.46%,高中及以上人群占总死亡人数的6.50%。按死因分组,每天非意外、呼吸系统疾病、循环系统疾病死亡人数分别为42.90、5.29、20.66人。呼吸系统疾病、循环系统疾病死亡人数分别占总死亡人数的10.71%、41.83%,高血压、冠心病、中风分别占6.17%、14.34%、19.58%。

    表  1  2015—2021年华中某市人群每日死亡描述性分析
    Table  1.  Descriptive summary for daily deathsal data in a city of central China from 2015 to 2021
    因素n$\bar x \pm s $MinP25P50P75Max
    总体 117001 49.39±11.47 22.00 41.00 48.00 56.00 116.00
    性别
     男性 65716 27.74±7.02 11.00 23.00 27.00 32.00 76.00
     女性 51285 21.65±6.48 7.00 17.00 21.00 25.00 51.00
    年龄/岁
     6~65 29268 12.36±3.87 2.00 10.00 12.00 15.00 34.00
     >65 87124 36.78±9.98 13.00 30.00 36.00 42.00 98.00
    教育
     高中及以上 7609 3.21±1.91 0.00 2.00 3.00 4.00 14.00
     高中以下 109392 46.18±11.01 21.00 38.00 45.00 52.00 107.00
    死因
     非意外 101635 42.90±10.63 17.00 36.00 41.00 49.00 108.00
     呼吸系统疾病 12536 5.29±3.04 0.00 3.00 5.00 7.00 21.00
     循环系统疾病 48947 20.66±6.62 4.00 16.00 20.00 24.00 68.00
     高血压 7221 3.05±2.01 0.00 2.00 3.00 4.00 12.00
     冠心病 16782 7.08±3.25 0.00 5.00 7.00 9.00 22.00
     中风 22906 9.67±3.84 0.00 7.00 9.00 12.00 41.00
    下载: 导出CSV 
    | 显示表格

    2015—2021年研究地区大气污染和气象指标的基本情况可见表2。该市的O3 8 h平均浓度为94.38 μɡ·m−3,PM2.5、PM10、SO2、NO2、CO的日均浓度分别为55.56 μɡ·m−3、86.60 μɡ·m−3、15.38 μɡ·m−3、31.20 μɡ·m−3、0.86 mɡ·m−3。研究期间日均气温、湿度、气压分别为17.05 ℃、77.20%、1004.61 hPa。

    表  2  2015—2021年华中某市大气污染、气象因素描述性分析(n=547)
    Table  2.  Descriptive summary for air pollution and meteorological data in a city of central China from 2015 to 2021(n=547)
    因素$ \bar x \pm s $MinP25P50P75Max
    污染物
     O3/(μɡ·m−3) 94.38±36.95 11.50 66.50 91.75 119.50 223.25
     PM2.5/(μɡ·m−3) 55.56±39.59 2.50 28.25 44.00 71.46 298.00
     PM10/(μɡ·m−3) 86.60±50.76 5.00 52.00 76.00 108.56 413.50
     SO2/(μɡ·m−3) 15.38±11.22 0.75 7.81 12.50 20.25 115.25
     NO2/(μɡ·m−3) 31.20±14.98 6.50 20.00 28.25 39.25 117.33
     CO/(mɡ·m−3) 0.86±0.29 0.28 0.67 0.80 1.00 2.60
     气温/℃ 17.05±8.88 −5.50 9.20 17.70 24.80 33.80
     湿度/% 77.20±14.62 27.00 67.00 79.00 89.00 100.00
     气压/hPa 1004.61±9.54 983.70 996.00 1004.10 1011.40 1032.10
    暖季
     O3/(μɡ·m−3) 109.17±35.04 12.50 83.75 107.50 133.25 223.25
     PM2.5/(μɡ·m−3) 35.26±20.22 2.50 21.50 31.25 44.00 148.00
    冷季
     O3/(μɡ·m−3) 80.20±32.97 11.50 53.75 77.50 103.00 179.50
     PM2.5/(μɡ·m−3) 75.10±43.40 8.00 43.75 63.75 96.50 298.00
    下载: 导出CSV 
    | 显示表格

    经过分析可知,全年中5—10月的平均温度最高。因此,本研究将5—10月定义为暖季,11月—次年4月定义为冷季。该市在寒冷季节死亡人数较多,而暖季较少。研究期间该市的臭氧浓度呈周期性变化趋势(见图1),一年中呈倒“V”型分布,暖季O3、PM2.5平均浓度分别为109.17 μɡ·m−3、35.26 μɡ·m−3,冷季平均浓度为80.20 μɡ·m−3、75.10 μɡ·m−3(见表2)。暖季的臭氧浓度明显高于冷季,而PM2.5浓度的季节分布情况恰好跟臭氧相反。研究期间,O3的超标率为4.43%,PM2.5的超标率为22.67%。

    图  1  华中某市2015—2021年O3浓度(A)、PM2.5浓度(B)、人群死亡(C)的周期性变化
    Figure  1.  Time-series of O3 concentration (A), PM2.5 concentration (B), and population death (C) in a city in central China from 2015 to 2021

    O3与气温存在正相关,而与湿度、气压之间存在负相关(P<0.01);O3与PM2.5、NO2、CO存在负相关(P<0.01),与SO2不具有统计学意义的关联。PM2.5与除O3外的污染物呈正相关关系,与气温、湿度呈负相关,与气压之间存在正相关(P<0.01)(见表3)。通过对污染物进行共线性诊断,方差膨胀因子小于4,共线性可接受。

    表  3  2015—2021年华中某市大气污染物与气象指标的相关性(r
    Table  3.  Correlation between air pollutants and meteorological factors in a city of central China from 2015 to 2021 (r)
    指标O3PM2.5SO2NO2CO气温湿度气压
    O31.000
    PM2.5−0.198*1.000
    SO20.0340.430*1.000
    NO2−0.154*0.614*0.608*1.000
    CO−0.145*0.646*0.470*0.576*1.000
    气温0.579*−0.556*−0.210*−0.450*−0.291*1.000
    湿度−0.438*−0.167*−0.461*−0.392*−0.037−0.0331.000
    气压−0.434*0.547*0.322*0.489*0.248*−0.868*−0.195*1.000
    [注]*:P<0.01(双侧)。
    下载: 导出CSV 
    | 显示表格

    大气污染物O3每升高10 μɡ·m−3,总死亡(lag02)、呼吸系统疾病死亡(lag2)、循环系统疾病死亡(lag02)的ER(95%CI)分别为0.09%(−0.08%~0.25%)、0.35%(0~0.71%)、0.43%(0.18%~0.68%)。大气污染物PM2.5每升高10 μɡ·m−3,总死亡(lag1)、呼吸系统疾病死亡(lag5)、循环系统疾病死亡(lag01)的ER(95%CI)分别为0.02%(−0.09%~0.13%)、0.25%(0.01%~0.50%)、0.35%(0.16%~0.54%)。O3和PM2.5均对人群总死亡影响无统计学意义(P>0.05)。见表4。随着滞后时间的延长,在一定程度上O3和PM2.5均对人群总死亡具有保护作用。

    表  4  2015—2021年华中某市大气污染物(O3、PM2.5)每升高10 μɡ·m−3居民死亡滞后效应(ER及其95%CI)
    Table  4.  Lag effects on daily deaths associated with per 10 μɡ·m−3 increase in pollutant (O3 and PM2.5) concentrations in a city of central China from 2015 to 2021 (ER and its 95%CI) 单位:%
    滞后期O3PM2.5
    非意外总死亡呼吸系统疾病循环系统疾病非意外总死亡呼吸系统疾病循环系统疾病
    Lag00.04(−0.10~0.18)−0.56(−0.98~−0.13)0.36(0.14~0.58)*−0.03(−0.14~0.08)−0.62(−0.94~−0.30)0.29(0.13~0.46)*
    Lag10.07(−0.06~0.19)−0.04(−0.44~0.36)0.28(0.08~0.48)*0.02(−0.09~0.13)−0.08(−0.39~0.24)0.25(0.09~0.42)*
    Lag20.06(−0.06~0.17)0.35(0.00~0.71)0.10(−0.08~0.28)−0.06(−0.15~0.05)−0.06(−0.36~0.24)0.06(−0.10~0.21)
    Lag3−0.19(−0.30~−0.07)−0.71(−1.06~−0.36)*−0.12(−0.30~0.06)−0.15(−0.25~−0.04)−0.14(−0.44~0.16)−0.04(−0.20~0.12)
    Lag4−0.18(−0.30~−0.07)−0.66(−1.00~−0.31)−0.10(−0.25~0.08)−0.18(−0.30~−0.07)−0.66(−1.00~−0.31)−0.10(−0.28~0.08)
    Lag5−0.33(−0.45~−0.22)**−0.94(−1.29~−0.59)**−0.13(−0.31~0.04)−0.19(−0.30~−0.09)0.25(0.01~0.50)*−0.21(−0.36~−0.05)
    Lag010.07(−0.09~0.23)−0.38(−0.85~0.10)0.43(0.18~0.68)*−0.02(−0.14~0.11)−0.46(−0.83~−0.10)0.35(0.16~0.54)*
    Lag020.09(−0.08~0.25)−0.06(−0.56~0.10)0.37(0.12~0.63)*−0.05(−0.18~0.09)−0.41(−0.81~−0.01)0.31(0.10~0.51)
    Lag03−0.04(−0.21~0.13)−0.43(−0.94~0.09)0.23(−0.03~0.50)−0.11(−0.26~0.04)−0.43(−0.86~−0.01)0.25(0.03~0.48)
    Lag04−0.12(−0.30~0.05)−0.68(−1.21~−0.15)0.15(−0.12~0.42)−0.25(−0.41~−0.09)−0.33(−0.79~0.13)0.12(−0.12~0.36)
    Lag05−0.26(−0.44~−0.08)−1.01(−1.56~−0.46)0.07(−0.21~0.35)−0.32(−0.49~−0.15)−0.21(−0.70~0.28)0.03(−0.23~0.29)
    [注]*:P<0.05;**:P<0.01。
    下载: 导出CSV 
    | 显示表格

    在对居民特定疾病分析发现,大气O3每升高10 μɡ·m−3,高血压(lag0)、冠心病(lag02)、中风(lag0)死亡的ER(95%CI)分别为0.45%(0.02%~0.91%)、0.59%(0.16%~1.02%)、0.33%(0.01%~0.65%);大气中PM2.5每升高10 μɡ·m−3,高血压(lag05)、冠心病(lag4)、中风(lag01)死亡的ER(95%CI)分别为1.18%(0.59%~1.77%)、0.17%(−0.13%~0.40%)、0.65%(0.38%~0.92%),除冠心病外均具有统计学意义(P<0.05)。见表5

    表  5  2015—2021年华中某市大气污染物(O3、PM2.5)每升高10 μɡ·m−3居民特定疾病死亡滞后效应(ER及95%CI)
    Table  5.  Lag effects on specific disease deaths associated with per 10 μɡ·m−3 increase in pollutant (O3 and PM2.5) concentrations in a city of central China from 2015 to 2021 (ER and its 95%CI) 单位:%
    滞后期O3PM2.5
    高血压冠心病中风高血压冠心病中风
    Lag00.45(0.02~0.91)*0.30(−0.07~0.67)0.33(0.01~0.65)*0.93(0.53~1.32)*0.11(−0.17~0.39)0.43(0.19~0.67)*
    Lag1−0.07(−0.56~0.42)0.53(0.19~0.86)*0.09(−0.20~0.37)0.75(0.36~1.14)*0.01(−0.26~0.29)0.56(0.32~0.79)*
    Lag2−0.65(−1.10~−0.20)0.29(−0.02~0.59)0.05(−0.21~0.32)0.47(0.08~0.86)*−0.14(−0.40~0.13)0.24(0.01~0.47)*
    Lag3−0.29(−0.73~0.15)−0.03(−0.33~0.27)−0.30(−0.56~−0.05)0.09(−0.30~0.48)0.02(−0.25~0.28)0.05(−0.18~0.28)
    Lag4−0.05(−0.48~0.39)0.17(−0.13~0.46)−0.44(−0.69~−0.18)−0.05(−0.48~0.39)0.17(−0.13~0.46)−0.44(−0.69~−0.18)
    Lag5−0.06(−0.49~0.37)−0.41(−0.70~−0.11)−0.06(−0.31~0.19)0.51(0.12~0.90)*−0.51(−0.78~−0.24)−0.02(−0.25~0.21)
    Lag010.21(−0.38~0.80)0.57(0.16~0.97)*0.26(−0.10~0.62)1.06(0.62~1.50)*0.06(−0.26~0.38)0.65(0.38~0.92)*
    Lag02−0.24(−0.86~0.37)0.59(0.16~1.02)*0.21(−0.16~0.57)1.09(0.62~1.57)*−0.03(−0.37~0.31)0.63(0.34~0.93)*
    Lag03−0.35(−0.98~0.29)0.45(0.01~0.89)*0.06(−0.32~0.44)0.94(0.43~1.46)*−0.03(−0.40~0.34)0.57(0.25~0.89)*
    Lag04−0.31(−0.98~0.34)0.46(0.01~0.92)*−0.21(−0.60~0.18)1.03(0.48~1.59)*−0.13(−0.53~0.27)0.39(0.04~0.73)*
    Lag05−0.30(−0.96~0.37)0.23(−0.23~0.70)−0.21(−0.61~0.19)1.18(0.59~1.77)*−0.33(−0.75~0.10)0.38(0.01~0.74)*
    [注]*:P<0.05。
    下载: 导出CSV 
    | 显示表格

    冷季O3每升高10 μɡ·m−3,总死亡(lag1)以及呼吸系统疾病(lag1)、循环系统疾病(lag02)、高血压(lag01)、冠心病(lag02)、中风(lag5)死亡的ER(95%CI)分别为0.03%(−0.17%~0.23%)、0.15%(−0.45%~0.74%)、0.21%(−0.21%~0.62%)、0.55(−0.35%~1.45%)、0.50%(−0.18%~1.19%)、0.43%(0.05%~0.82%),除中风(P<0.05)外,臭氧对其他疾病死亡的影响均无统计学意义(P>0.05)。暖季O3每升高10 μɡ·m−3,总死亡(lag02)以及呼吸系统疾病(lag2)、循环系统疾病(lag01)、高血压(lag0)、冠心病(lag01)、中风(lag0)死亡的ER(95%CI)分别为0.33%(0.10%~0.55%)、0.35%(−0.13%~0.83%)、0.74%(0.40%~1.08%)、0.72%(0.02%~1.48%)、0.81%(0.25%~1.37%)、0.70%(0.25%~1.15%),除呼吸系统以外,其他结果均存在统计学意义(P<0.05)。见图2

    图  2  不同季节大气污染物O3每升高10 μɡ·m−3居民各死因死亡的滞后效应(ER及其95%CI)
    A:表示冷季(11月—次年4月);B:表示暖季(5—10月)。
    Figure  2.  Season-specific lag effects on daily deaths associated with per 10 μɡ·m−3 increase in O3 concentrations (ER and its 95%CI)

    冷季大气污染物PM2.5每升高10 μɡ·m−3,总死亡(lag02)以及呼吸系统疾病(lag5)、循环系统疾病(lag02)、高血压(lag05)、冠心病(lag3)、中风(lag02)死亡的ER(95%CI)分别为−0.04%(−0.16%~0.08%)、0.36%(0.04%~0.69%)、0.14%(−0.10%~0.37%)、1.76%(1.20%~2.32%)、0.01%(−0.28%~0.30%)、0.54%(0.20%~0.88%),臭氧对高血压和中风死亡的影响有统计学意义(P<0.05)。暖季PM2.5每升高10 μɡ·m−3,总死亡(lag01)以及呼吸系统疾病(lag2)、循环系统疾病(lag01)、高血压(lag05)、冠心病(lag01)、中风(lag02)死亡的ER(95%CI)分别为0.60%(0.26%~0.93%)、1.71%(0.88%~2.53%)、1.66%(1.16%~2.17%)、2.51%(0.95%~4.16%)、1.92%(1.10%~2.74%)、1.63%(0.85%~2.42%),均具有统计学意义(P<0.05)。暖季时,PM2.5对各个系统疾病或特定疾病死亡的影响最大值多出现在lag01或lag02,见图3

    图  3  不同季节大气污染物PM2.5每升高10 μɡ·m−3各死因死亡的滞后效应(ER及其95%CI)
    A:冷季(11月—次年4月);B:暖季(5—10月)。
    Figure  3.  Season-specific lag effects on daily deaths associated with per 10 μɡ·m−3 increase in PM2.5 concentrations (ER and its 95%CI)

    2015—2021年大气污染物O3、PM2.5暴露与居民死亡存在较为明显的季节性差异,污染物(O3、PM2.5)效应在暖季较为明显。此外,O3、PM2.5对特定疾病死亡的影响高于系统疾病死亡与总死亡。

    通过季节效应分析发现,污染物在暖季时效应更高,因此,此研究选择在暖季时进行分层分析。大气O3每升高10 μɡ·m−3,在lag02时,男性总死亡增加0.46%(95%CI:0.16%~0.75%),在lag01时女性总死亡增加0.16%(95%CI:-0.17%~0.49%),在lag4时6~65岁人群总死亡增加0.38%(95%CI:0.08%~0.68%),在lag01时>65岁人群总死亡增加0.41%(95%CI:0.14%~0.66%),在lag4时高中及以上文化人群总死亡增加-0.13%(95%CI:−0.85%~0.59%),在lag02时高中以下文化人群总死亡增加0.38%(95%CI:0.14%~0.61%)。 O3对女性总死亡的影响没有统计学意义(P>0.05)。男性死亡效应高于女性,>65岁人群死亡效应高于6~65岁人群,受教育程度低的人群死亡效应高于受教育程度高的人群。见图4

    图  4  暖季大气污染物O3每升高10 μɡ·m−3居民总死亡的滞后效应(ER及其95%CI)
    Figure  4.  Lag effects on total deaths associated with per 10 μɡ·m−3 increase in O3 concentrations in warm season (ER and its 95%CI)

    PM2.5每升高10 μɡ·m−3,在lag01时男性总死亡增加0.54%(95%CI:0.09%~0.98%),在lag04时女性总死亡增加1.01%(95%CI:0.35%~1.67%),在lag04时6~65岁人群总死亡增加0.41%(95%CI:0.13%~0.70%),在lag03时>65岁人群总死亡增加0.89%(95%CI:0.41%~1.37%),在lag4时高中以上文化人群总死亡增加−0.06%(95%CI:−0.60%~0.48%),在lag03时高中以下文化人群的总死亡增加0.73%(95%CI:0.31%~1.16%)。 PM2.5对高中以上文化人群总死亡的影响没有统计学意义(P>0.05)。在暖季时,PM2.5暴露导致女性的死亡效应高于男性,>65岁人群的死亡效应高于6~65岁人群,受教育程度低的人群死亡效应高于受教育程度高的人群。见图5

    图  5  暖季大气污染物PM2.5每升高10 μɡ·m−3居民总死亡滞后效应(ER及其95%CI)
    Figure  5.  Lag effects on total deaths associated with per 10 μɡ·m−3 increase in PM2.5 concentrations in warm season (ER and its 95%CI)

    O3与PM2.5的双污染物及多污染物模型结果与单污染物模型结果基本一致,且95%CI变化也不大,说明此模型稳定,结果可靠。见表6

    表  6  空气污染物(O3、PM2.5)每升高10 μɡ·m−3对人群死亡的双污染物及多污染物模型(ER及其95%CI)
    Table  6.  Daily deaths associated with per 10 μɡ·m−3 increase in pollutant (O3 and PM2.5) concentrations in double pollutant and multi-pollutant models (ER and its 95%CI) 单位:%
    调整因素总死亡呼吸系统疾病死亡循环系统疾病死亡
    单O3模型0.09(−0.08~0.25)0.35(0~0.71)0.43(0.18~0.68)*
    O3+PM2.50.10(−0.07~0.27)0.45(0.08,0.81)*0.35(0.10~0.60)*
    O3+SO20.09(−0.08~0.25)0.34(−0.02~0.70)0.42(0.18~0.67)*
    O3+NO20.05(−0.11~0.21)0.35(−0.01~0.71)0.41(0.16~0.66)*
    O3+CO0.08(−0.09~0.23)0.35(−0.01~0.71)0.39(0.14~0.64)*
    O3+PM2.5+SO2+NO2+CO0.10(−0.07~0.27)0.44(0.07~0.80)*0.35(0.10~0.61)*
    单PM2.5模型0.02(−0.09~0.13)0.25(0.01~0.50)*0.35(0.16~0.54)*
    PM2.5+O30.01(−0.10~0.12)0.25(−0.05~0.55)0.30(0.11~0.50)*
    PM2.5+SO20.02(−0.09~0.12)0.24(−0.06~0.53)0.34(0.15~0.53)*
    PM2.5+NO2−0.01(−0.11~0.10)0.25(−0.05~0.55)0.31(0.12~0.51)*
    PM2.5+CO0.01(−0.12~0.12)0.25(−0.05~0.55)0.32(0.10~0.54)*
    PM2.5+O3+SO2+NO2+CO0.01(−0.10~0.13)0.25(−0.04~0.55)0.28(0.05~0.51)*
    [注]在滞后效应最强时加入其他污染物;*:P<0.05。
    下载: 导出CSV 
    | 显示表格

    O3单污染模型单独引入PM2.5、SO2、NO2、CO或同时引入PM2.5、SO2、NO2、CO时,对居民总死亡率影响不大;与O3单污染模型相比在单独引入PM2.5或同时引入PM2.5、SO2、NO2、CO时,会增加呼吸系统疾病死亡;与O3单污染模型相比在单独引入PM2.5或同时引入PM2.5、SO2、NO2、CO时,会降低循环系统疾病死亡。PM2.5单污染模型在单独引入O3、SO2、NO2、CO或同时引入O3、SO2、NO2、CO时,对居民总死亡和呼吸系统疾病死亡影响不大;PM2.5单污染模型单独引入O3或者同时引入O3、SO2、NO2、CO时,与PM2.5单污染模型相比会降低循环系统疾病死亡。

    该市处于华中地区,既是长江中游城市群的重要成员,也是鄂中区域性城市。该市的工业污染主要来源于建材行业、化工行业和电力行业[20],因此会导致一定程度的的空气污染。呼吸系统作为空气污染物进入人体的最主要途径,可直接与污染物接触并受到侵害[21]。研究表明,臭氧暴露会导致很多循环生物标志发生明显变化,包括炎症、氧化应激、血管反应性等[22]。多项研究证明,O3和PM2.5会增加呼吸系统疾病和循环系统疾病的死亡率[4623-24],而心脑血管疾病已经是危害我国公众健康的重要公共卫生问题之一,但地区差异明显。也有研究报道,短期的臭氧暴露不会损害血管功能,升高血压或者影响心率[25-26]。因此,明确呼吸系统疾病和循环系统疾病的相关危险因素十分重要,故本次研究将O3、PM2.5和居民死亡效应作为观察对象。

    本研究通过收集2015—2021年华中某市大气污染、气象指标、居民死亡个案的日报数据,基于广义相加模型分析O3、PM2.5对居民死亡的影响。结果显示,2015—2021年O3的平均浓度为94.38 μɡ·m−3,这跟包头市(93.94 μɡ·m−3[27]、北京市顺义区(94.60 μɡ·m−3[28]的结果相似,但这一结果略低于广州市(98.10 μɡ·m−3[6]的研究;O3和PM2.5的分布具有明显的季节差异,O3浓度呈倒“V”型分布,冬季低夏季高;而PM2.5的分布跟O3相反,呈冬季高夏季低。

    本研究结果显示,在全年分析中,O3、PM2.5对居民总死亡的影响没有统计学意义,而随着滞后时间越长,O3、PM2.5对居民总死亡具有保护作用。这可能是由于大气污染物的短期滞后效应[5],或是本次样本量较小,结果易波动,还有其他不确定因素,需进行深入研究。但在暖季时,两种污染物对居民总死亡的影响均有统计学意义,可认为O3、PM2.5对居民总死亡的影响受到温度的影响。气温越高,颗粒物与温度对人群死亡率的联合作用就越大[29]

    本研究显示,O3浓度每升高10 μɡ·m−3,在lag2时,呼吸系统疾病死亡的ER为0.35%(0~0.71%),低于广州市2013—2018年研究结果0.79%(0.32%~1.26%)[23]、广州市2006—2016研究结果0.68%(0.39%~0.97%)[6],与上海市2001—2004年、香港1999—2003年研究结果相似[30-32]。对于循环系统,O3每升高10 μɡ·m−3,在lag01时,死亡ER(95%CI)为0.43%(0.18%~0.68%),低于包头市2015—2017年的研究结果1.62%(0.38%~2.88%)[27],也低于南京市2013—2017年的研究结果0.52%(0.31%~0.73%)[33]。产生这种差异的原因有:①每个地区大气污染物浓度状况存在差异[34];②由于每个地区的气候变化,各地的高温持续时长不一样,温度是影响疾病死亡的重要因素之一[35-37];③GAM模型中的参数选择差异,以及人群生活行为方式的不同[38]

    高血压、冠心病、中风是最常见的几种心脑血管疾病。本研究结果显示,O3每升高10 μɡ·m−3,三种疾病的死亡ER(95%CI)都会明显升高,最大效应出现在lag0或lag02时,说明臭氧对特定疾病死亡存在急性影响,这与芦静等[27]、陈琦等[39]研究结果类似。在季节分层分析时,暖季对居民死亡的影响明显高于冷季,表明臭氧效应存在明显的季节效应。研究发现高温与O3暴露会通过加重炎性反应、氧化应激以及诱导血管收缩等途径使高血压疾病加重[40]。另一项研究发现,高温会增加O3暴露所引起的心脏毒性,使心血管疾病加重[41]。相比于O3,PM2.5对心脑血管疾病死亡的影响更大,这与合肥市的结果相似[24]。大量研究证明,PM2.5可通过氧化应激反应、免疫调节、炎症反应等途径对心肌细胞产生影响[42],从而增加心血管系统疾病的死亡。本研究表明,O3、PM2.5暴露对循环系统疾病以及高血压、冠心病、中风死亡的影响明显高于总死亡,因此,患有这些基础疾病的人群为O3、PM2.5暴露的敏感人群。

    在对总死亡进行年龄、性别、教育程度分层分析时发现,O3每升高10 μɡ·m−3,男性死亡效应高于女性,与Niu等[43]的全国研究结果相似;此外,>65人群死亡效应高于6~65岁人群,高中以下文化人群死亡效应高于高中及以上文化人群。除了性别分层结果以外,其他结果均与Zhang等[4]、班婕等[38]、Lei等[44]的研究结果类似。PM2.5每升高10 μɡ·m−3,女性死亡效应高于男性,>65岁人群死亡效应高于6~65岁人群,高中以下文化人群死亡效应高于高中及以上文化人群,这与Cheng等[24]在合肥的研究结果类似。O3、PM2.5对高中以上文化人群的死亡率影响较低,且不具有统计学意义(P>0.05),可能是由于:①人群的生活方式跟工作地点差异,教育程度低的人群工作的地点污染源多一点,或者户外暴露时间较长;②本次研究人群中高中以上文化占比很低,结果可能不准确。女性死亡比男性高可能是由于生理特征、生活方式的差异[45],女性对大气污染物更敏感[46]。大气污染物对65岁以上人群死亡影响大,可能由于老年人器官衰弱以及功能减退,抵抗力以及机体修复能力降低[4]

    在敏感性分析中,双污染和多污染模型的结果与单污染模型的结果基本一致,说明模型的稳定性较好。当O3单污染物模型单独引入PM2.5时,会增加呼吸系统疾病的死亡,降低循环系统疾病的死亡。出现这种结果的原因有:PM2.5可直接进入呼吸道中,直接对呼吸系统产生影响,PM2.5对呼吸系统疾病的影响可能大于O3;O3可直接进入血液中影响循环系统,O3对循环系统疾病的影响可能大于PM2.5;黏附在人体皮肤表明的PM2.5(金属组分,氧化还原反应)可能会消耗一部分臭氧,从而导致臭氧对循环系统的影响降低。O3与PM2.5的浓度分布相反,两者之间为负相关,当O3浓度升高时,PM2.5浓度降低,且O3浓度在夏季时达到峰值,而PM2.5浓度在冬季时达到峰值,会降低两种污染物对人群健康影响的联合作用[25],这可能是ER下降的部分原因。也有研究发现,O3与PM2.5联合作用时,会加重小鼠的心血管疾病和肺损伤[47]

    此次研究存在一定的局限性:首先,使用的大气污染物浓度为大气监测站浓度,不能反映居民的真实暴露浓度,会对结果产生误差;其次,研究使用的人群死亡资料来源于疾控中心,可能存在漏报情况;最后,本次研究是基于群体水平而不是个体水平探讨大气污染物与居民死亡之间的关系,可能存在生态学谬论。

    综上所述,O3、PM2.5暴露对该市居民死亡存在急性影响,其中对循环系统疾病(包括高血压、冠心病、中风)的影响较大,大于65岁的老年人群以及受教育程度低的人群为敏感人群。因此,此次研究进一步可为当地环保部门管控大气污染提供科学依据;并提示政府部门可重点对敏感人群进行宣传教育,避免敏感人群暴露在高浓度O3或PM2.5中,从而降低大气污染物对居民死亡效应影响。

  • 图  1   华中某市2015—2021年O3浓度(A)、PM2.5浓度(B)、人群死亡(C)的周期性变化

    Figure  1.   Time-series of O3 concentration (A), PM2.5 concentration (B), and population death (C) in a city in central China from 2015 to 2021

    图  2   不同季节大气污染物O3每升高10 μɡ·m−3居民各死因死亡的滞后效应(ER及其95%CI)

    A:表示冷季(11月—次年4月);B:表示暖季(5—10月)。

    Figure  2.   Season-specific lag effects on daily deaths associated with per 10 μɡ·m−3 increase in O3 concentrations (ER and its 95%CI)

    图  3   不同季节大气污染物PM2.5每升高10 μɡ·m−3各死因死亡的滞后效应(ER及其95%CI)

    A:冷季(11月—次年4月);B:暖季(5—10月)。

    Figure  3.   Season-specific lag effects on daily deaths associated with per 10 μɡ·m−3 increase in PM2.5 concentrations (ER and its 95%CI)

    图  4   暖季大气污染物O3每升高10 μɡ·m−3居民总死亡的滞后效应(ER及其95%CI)

    Figure  4.   Lag effects on total deaths associated with per 10 μɡ·m−3 increase in O3 concentrations in warm season (ER and its 95%CI)

    图  5   暖季大气污染物PM2.5每升高10 μɡ·m−3居民总死亡滞后效应(ER及其95%CI)

    Figure  5.   Lag effects on total deaths associated with per 10 μɡ·m−3 increase in PM2.5 concentrations in warm season (ER and its 95%CI)

    表  1   2015—2021年华中某市人群每日死亡描述性分析

    Table  1   Descriptive summary for daily deathsal data in a city of central China from 2015 to 2021

    因素n$\bar x \pm s $MinP25P50P75Max
    总体 117001 49.39±11.47 22.00 41.00 48.00 56.00 116.00
    性别
     男性 65716 27.74±7.02 11.00 23.00 27.00 32.00 76.00
     女性 51285 21.65±6.48 7.00 17.00 21.00 25.00 51.00
    年龄/岁
     6~65 29268 12.36±3.87 2.00 10.00 12.00 15.00 34.00
     >65 87124 36.78±9.98 13.00 30.00 36.00 42.00 98.00
    教育
     高中及以上 7609 3.21±1.91 0.00 2.00 3.00 4.00 14.00
     高中以下 109392 46.18±11.01 21.00 38.00 45.00 52.00 107.00
    死因
     非意外 101635 42.90±10.63 17.00 36.00 41.00 49.00 108.00
     呼吸系统疾病 12536 5.29±3.04 0.00 3.00 5.00 7.00 21.00
     循环系统疾病 48947 20.66±6.62 4.00 16.00 20.00 24.00 68.00
     高血压 7221 3.05±2.01 0.00 2.00 3.00 4.00 12.00
     冠心病 16782 7.08±3.25 0.00 5.00 7.00 9.00 22.00
     中风 22906 9.67±3.84 0.00 7.00 9.00 12.00 41.00
    下载: 导出CSV

    表  2   2015—2021年华中某市大气污染、气象因素描述性分析(n=547)

    Table  2   Descriptive summary for air pollution and meteorological data in a city of central China from 2015 to 2021(n=547)

    因素$ \bar x \pm s $MinP25P50P75Max
    污染物
     O3/(μɡ·m−3) 94.38±36.95 11.50 66.50 91.75 119.50 223.25
     PM2.5/(μɡ·m−3) 55.56±39.59 2.50 28.25 44.00 71.46 298.00
     PM10/(μɡ·m−3) 86.60±50.76 5.00 52.00 76.00 108.56 413.50
     SO2/(μɡ·m−3) 15.38±11.22 0.75 7.81 12.50 20.25 115.25
     NO2/(μɡ·m−3) 31.20±14.98 6.50 20.00 28.25 39.25 117.33
     CO/(mɡ·m−3) 0.86±0.29 0.28 0.67 0.80 1.00 2.60
     气温/℃ 17.05±8.88 −5.50 9.20 17.70 24.80 33.80
     湿度/% 77.20±14.62 27.00 67.00 79.00 89.00 100.00
     气压/hPa 1004.61±9.54 983.70 996.00 1004.10 1011.40 1032.10
    暖季
     O3/(μɡ·m−3) 109.17±35.04 12.50 83.75 107.50 133.25 223.25
     PM2.5/(μɡ·m−3) 35.26±20.22 2.50 21.50 31.25 44.00 148.00
    冷季
     O3/(μɡ·m−3) 80.20±32.97 11.50 53.75 77.50 103.00 179.50
     PM2.5/(μɡ·m−3) 75.10±43.40 8.00 43.75 63.75 96.50 298.00
    下载: 导出CSV

    表  3   2015—2021年华中某市大气污染物与气象指标的相关性(r

    Table  3   Correlation between air pollutants and meteorological factors in a city of central China from 2015 to 2021 (r)

    指标O3PM2.5SO2NO2CO气温湿度气压
    O31.000
    PM2.5−0.198*1.000
    SO20.0340.430*1.000
    NO2−0.154*0.614*0.608*1.000
    CO−0.145*0.646*0.470*0.576*1.000
    气温0.579*−0.556*−0.210*−0.450*−0.291*1.000
    湿度−0.438*−0.167*−0.461*−0.392*−0.037−0.0331.000
    气压−0.434*0.547*0.322*0.489*0.248*−0.868*−0.195*1.000
    [注]*:P<0.01(双侧)。
    下载: 导出CSV

    表  4   2015—2021年华中某市大气污染物(O3、PM2.5)每升高10 μɡ·m−3居民死亡滞后效应(ER及其95%CI)

    Table  4   Lag effects on daily deaths associated with per 10 μɡ·m−3 increase in pollutant (O3 and PM2.5) concentrations in a city of central China from 2015 to 2021 (ER and its 95%CI) 单位:%

    滞后期O3PM2.5
    非意外总死亡呼吸系统疾病循环系统疾病非意外总死亡呼吸系统疾病循环系统疾病
    Lag00.04(−0.10~0.18)−0.56(−0.98~−0.13)0.36(0.14~0.58)*−0.03(−0.14~0.08)−0.62(−0.94~−0.30)0.29(0.13~0.46)*
    Lag10.07(−0.06~0.19)−0.04(−0.44~0.36)0.28(0.08~0.48)*0.02(−0.09~0.13)−0.08(−0.39~0.24)0.25(0.09~0.42)*
    Lag20.06(−0.06~0.17)0.35(0.00~0.71)0.10(−0.08~0.28)−0.06(−0.15~0.05)−0.06(−0.36~0.24)0.06(−0.10~0.21)
    Lag3−0.19(−0.30~−0.07)−0.71(−1.06~−0.36)*−0.12(−0.30~0.06)−0.15(−0.25~−0.04)−0.14(−0.44~0.16)−0.04(−0.20~0.12)
    Lag4−0.18(−0.30~−0.07)−0.66(−1.00~−0.31)−0.10(−0.25~0.08)−0.18(−0.30~−0.07)−0.66(−1.00~−0.31)−0.10(−0.28~0.08)
    Lag5−0.33(−0.45~−0.22)**−0.94(−1.29~−0.59)**−0.13(−0.31~0.04)−0.19(−0.30~−0.09)0.25(0.01~0.50)*−0.21(−0.36~−0.05)
    Lag010.07(−0.09~0.23)−0.38(−0.85~0.10)0.43(0.18~0.68)*−0.02(−0.14~0.11)−0.46(−0.83~−0.10)0.35(0.16~0.54)*
    Lag020.09(−0.08~0.25)−0.06(−0.56~0.10)0.37(0.12~0.63)*−0.05(−0.18~0.09)−0.41(−0.81~−0.01)0.31(0.10~0.51)
    Lag03−0.04(−0.21~0.13)−0.43(−0.94~0.09)0.23(−0.03~0.50)−0.11(−0.26~0.04)−0.43(−0.86~−0.01)0.25(0.03~0.48)
    Lag04−0.12(−0.30~0.05)−0.68(−1.21~−0.15)0.15(−0.12~0.42)−0.25(−0.41~−0.09)−0.33(−0.79~0.13)0.12(−0.12~0.36)
    Lag05−0.26(−0.44~−0.08)−1.01(−1.56~−0.46)0.07(−0.21~0.35)−0.32(−0.49~−0.15)−0.21(−0.70~0.28)0.03(−0.23~0.29)
    [注]*:P<0.05;**:P<0.01。
    下载: 导出CSV

    表  5   2015—2021年华中某市大气污染物(O3、PM2.5)每升高10 μɡ·m−3居民特定疾病死亡滞后效应(ER及95%CI)

    Table  5   Lag effects on specific disease deaths associated with per 10 μɡ·m−3 increase in pollutant (O3 and PM2.5) concentrations in a city of central China from 2015 to 2021 (ER and its 95%CI) 单位:%

    滞后期O3PM2.5
    高血压冠心病中风高血压冠心病中风
    Lag00.45(0.02~0.91)*0.30(−0.07~0.67)0.33(0.01~0.65)*0.93(0.53~1.32)*0.11(−0.17~0.39)0.43(0.19~0.67)*
    Lag1−0.07(−0.56~0.42)0.53(0.19~0.86)*0.09(−0.20~0.37)0.75(0.36~1.14)*0.01(−0.26~0.29)0.56(0.32~0.79)*
    Lag2−0.65(−1.10~−0.20)0.29(−0.02~0.59)0.05(−0.21~0.32)0.47(0.08~0.86)*−0.14(−0.40~0.13)0.24(0.01~0.47)*
    Lag3−0.29(−0.73~0.15)−0.03(−0.33~0.27)−0.30(−0.56~−0.05)0.09(−0.30~0.48)0.02(−0.25~0.28)0.05(−0.18~0.28)
    Lag4−0.05(−0.48~0.39)0.17(−0.13~0.46)−0.44(−0.69~−0.18)−0.05(−0.48~0.39)0.17(−0.13~0.46)−0.44(−0.69~−0.18)
    Lag5−0.06(−0.49~0.37)−0.41(−0.70~−0.11)−0.06(−0.31~0.19)0.51(0.12~0.90)*−0.51(−0.78~−0.24)−0.02(−0.25~0.21)
    Lag010.21(−0.38~0.80)0.57(0.16~0.97)*0.26(−0.10~0.62)1.06(0.62~1.50)*0.06(−0.26~0.38)0.65(0.38~0.92)*
    Lag02−0.24(−0.86~0.37)0.59(0.16~1.02)*0.21(−0.16~0.57)1.09(0.62~1.57)*−0.03(−0.37~0.31)0.63(0.34~0.93)*
    Lag03−0.35(−0.98~0.29)0.45(0.01~0.89)*0.06(−0.32~0.44)0.94(0.43~1.46)*−0.03(−0.40~0.34)0.57(0.25~0.89)*
    Lag04−0.31(−0.98~0.34)0.46(0.01~0.92)*−0.21(−0.60~0.18)1.03(0.48~1.59)*−0.13(−0.53~0.27)0.39(0.04~0.73)*
    Lag05−0.30(−0.96~0.37)0.23(−0.23~0.70)−0.21(−0.61~0.19)1.18(0.59~1.77)*−0.33(−0.75~0.10)0.38(0.01~0.74)*
    [注]*:P<0.05。
    下载: 导出CSV

    表  6   空气污染物(O3、PM2.5)每升高10 μɡ·m−3对人群死亡的双污染物及多污染物模型(ER及其95%CI)

    Table  6   Daily deaths associated with per 10 μɡ·m−3 increase in pollutant (O3 and PM2.5) concentrations in double pollutant and multi-pollutant models (ER and its 95%CI) 单位:%

    调整因素总死亡呼吸系统疾病死亡循环系统疾病死亡
    单O3模型0.09(−0.08~0.25)0.35(0~0.71)0.43(0.18~0.68)*
    O3+PM2.50.10(−0.07~0.27)0.45(0.08,0.81)*0.35(0.10~0.60)*
    O3+SO20.09(−0.08~0.25)0.34(−0.02~0.70)0.42(0.18~0.67)*
    O3+NO20.05(−0.11~0.21)0.35(−0.01~0.71)0.41(0.16~0.66)*
    O3+CO0.08(−0.09~0.23)0.35(−0.01~0.71)0.39(0.14~0.64)*
    O3+PM2.5+SO2+NO2+CO0.10(−0.07~0.27)0.44(0.07~0.80)*0.35(0.10~0.61)*
    单PM2.5模型0.02(−0.09~0.13)0.25(0.01~0.50)*0.35(0.16~0.54)*
    PM2.5+O30.01(−0.10~0.12)0.25(−0.05~0.55)0.30(0.11~0.50)*
    PM2.5+SO20.02(−0.09~0.12)0.24(−0.06~0.53)0.34(0.15~0.53)*
    PM2.5+NO2−0.01(−0.11~0.10)0.25(−0.05~0.55)0.31(0.12~0.51)*
    PM2.5+CO0.01(−0.12~0.12)0.25(−0.05~0.55)0.32(0.10~0.54)*
    PM2.5+O3+SO2+NO2+CO0.01(−0.10~0.13)0.25(−0.04~0.55)0.28(0.05~0.51)*
    [注]在滞后效应最强时加入其他污染物;*:P<0.05。
    下载: 导出CSV
  • [1]

    WHO Air Quality Database 2022. https://cdn.who.int/media/docs/default-source/air-pollution-documents/air-quality-and-health/who-air-quality-database-2022---v7.pdf?sfvrsn=c6d52e7b_7&download=true.

    [2]

    LANDRIGAN P J, FULLER R, ACOSTA N J R, et al. The Lancet Commission on pollution and health[J]. Lancet, 2018, 391(10119): 462-512. doi: 10.1016/S0140-6736(17)32345-0

    [3]

    LIU H, TIAN Y, XU Y, et al. Ambient particulate matter concentrations and hospitalization for stroke in 26 Chinese cities: a case-crossover study[J]. Stroke, 2017, 48(8): 2052-2059. doi: 10.1161/STROKEAHA.116.016482

    [4]

    ZHANG J, CHEN Q, WANG Q, et al. The acute health effects of ozone and PM2.5 on daily cardiovascular disease mortality: a multi-center time series study in China[J]. Ecotoxicol Environ Saf, 2019, 174: 218-223. doi: 10.1016/j.ecoenv.2019.02.085

    [5]

    YIN P, CHEN R, WANG L, et al. Ambient ozone pollution and daily mortality: a nationwide study in 272 Chinese cities[J]. Environ Health Perspect, 2017, 125(11): 117006. doi: 10.1289/EHP1849

    [6]

    WU R, SONG X, CHEN D, et al. Health benefit of air quality improvement in Guangzhou, China: results from a long time-series analysis (2006-2016)[J]. Environ Int, 2019, 126: 552-559. doi: 10.1016/j.envint.2019.02.064

    [7]

    ZHANG R, LIU G, JIANG Y, et al. Acute effects of particulate air pollution on ischemic stroke and hemorrhagic stroke mortality[J]. Front Neurol, 2018, 9: 827. doi: 10.3389/fneur.2018.00827

    [8]

    XIAO Q, YANG L, MULHOLLAND J A, et al. Pediatric emergency department visits and ambient Air pollution in the U. S. State of Georgia: a case-crossover study[J]. Environ Health, 2016, 15(1): 115. doi: 10.1186/s12940-016-0196-y

    [9]

    SHI L, ZANOBETTI A, KLOOG I, et al. Low-concentration PM2.5 and mortality: estimating acute and chronic effects in a population-based study[J]. Environ Health Perspect, 2016, 124(1): 46-52. doi: 10.1289/ehp.1409111

    [10]

    BURNETT R T, POPE CA III, EZZATI M, et al. An integrated risk function for estimating the global burden of disease attributable to ambient fine particulate matter exposure[J]. Environ Health Perspect, 2014, 122: 397-403. doi: 10.1289/ehp.1307049

    [11]

    RIVA D R, MAGALHÃES C B, LOPES A A, et al. Low dose of fine particulate matter (PM2.5) can induce acute oxidative stress, inflammation and pulmonary impairment in healthy mice[J]. Inhal Toxicol, 2011, 23(5): 257-267. doi: 10.3109/08958378.2011.566290

    [12]

    WANG N, MENGERSEN K, TONG S, et al. Short-term association between ambient air pollution and lung cancer mortality[J]. Environ Res, 2019, 179: 108748. doi: 10.1016/j.envres.2019.108748

    [13]

    QIU H, TAN K, LONG F, et al. The burden of COPD morbidity attributable to the interaction between ambient air pollution and temperature in Chengdu, China[J]. Int J Environ Res Public Health, 2018, 15(3): 492. doi: 10.3390/ijerph15030492

    [14]

    SONG C, HE J, WU L, et al. Health burden attributable to ambient PM2.5 in China[J]. Environ Pollut, 2017, 223: 575-586. doi: 10.1016/j.envpol.2017.01.060

    [15] 董继元, 刘兴荣, 张本忠, 等. 我国臭氧短期暴露与人群死亡风险的Meta分析[J]. 环境科学学报, 2016, 36(4): 1477-1485. doi: 10.13671/j.hjkxxb.2015.0555

    DONG J Y, LIU X R, ZHANG B Z, et al. Meta-analysis of association between short-term ozone exposure and population mortality in China[J]. Acta Sci Circumstantiae, 2016, 36(4): 1477-1485. doi: 10.13671/j.hjkxxb.2015.0555

    [16]

    LIU C, CHEN R, SERA F, et al. Ambient particulate air pollution and daily mortality in 652 cities[J]. N Engl J Med, 2019, 381(8): 705-715. doi: 10.1056/NEJMoa1817364

    [17]

    DARROW L A, KLEIN M, SARNAT J A, et al. The use of alternative pollutant metrics in time-series studies of ambient air pollution and respiratory emergency department visits[J]. J Expo Sci Environ Epidemiol, 2011, 21(1): 10-19. doi: 10.1038/jes.2009.49

    [18]

    ZANOBETTI A, SCHWARTZ J. The effect of fine and coarse particulate air pollution on mortality: a national analysis[J]. Environ Health Perspect, 2009, 117(6): 898-903. doi: 10.1289/ehp.0800108

    [19]

    YAN Y, CHEN X, GUO Y, et al. Ambient air pollution and cerebrovascular disease mortality: an ecological time-series study based on 7-year death records in central China[J]. Environ Sci Pollut Res Int, 2021, 28(21): 27299-27307. doi: 10.1007/s11356-021-12474-z

    [20] 雷武琴, 袁红明, 赵鸣, 等. 荆门市大气环境质量变化及防治策略[J]. 中国资源综合利用, 2021, 39(4): 156-159. doi: 10.3969/j.issn.1008-9500.2021.04.047

    LEI W Q, YUAN H M, ZHAO M, et al. Changes of atmospheric environmental quality and control strategies in Jingmen City[J]. China Resour Compr Util, 2021, 39(4): 156-159. doi: 10.3969/j.issn.1008-9500.2021.04.047

    [21]

    ARBEX M A, DE PAULA SANTOS U, MARTINS L C, et al. Air pollution and the respiratory system[J]. J Bras Pneumol, 2012, 38(5): 643-655. doi: 10.1590/S1806-37132012000500015

    [22]

    GOODMAN J E, PRUEITT R L, SAX S N, et al. Ozone exposure and systemic biomarkers: evaluation of evidence for adverse cardiovascular health impacts[J]. Crit Rev Toxicol, 2015, 45(5): 412-452. doi: 10.3109/10408444.2015.1031371

    [23]

    LI M, DONG H, WANG B, et al. Association between ambient ozone pollution and mortality from a spectrum of causes in Guangzhou, China[J]. Sci Total Environ, 2021, 754: 142110. doi: 10.1016/j.scitotenv.2020.142110

    [24]

    CHENG H, ZHU F, LEI R, et al. Associations of ambient PM2.5 and O3 with cardiovascular mortality: a time-series study in Hefei, China[J]. Int J Biometeorol, 2019, 63(10): 1437-1447. doi: 10.1007/s00484-019-01766-2

    [25]

    BARATH S, LANGRISH J P, LUNDBÄCK M, et al. Short-term exposure to ozone does not impair vascular function or affect heart rate variability in healthy young men[J]. Toxicol Sci, 2013, 135(2): 292-299. doi: 10.1093/toxsci/kft157

    [26]

    HOFFMANN B, LUTTMANN-GIBSON H, COHEN A, et al. Opposing effects of particle pollution, ozone, and ambient temperature on arterial blood pressure[J]. Environ Health Perspect, 2012, 120(2): 241-246. doi: 10.1289/ehp.1103647

    [27] 芦静, 张晓梅, 冯晓冬. 包头市大气臭氧污染对人群循环系统疾病死亡的急性效应[J]. 环境与职业医学, 2019, 36(4): 381-387. doi: 10.13213/j.cnki.jeom.2019.18472

    LU J, ZHANG X M, FENG X D. Acute effect of ambient ozone pollution on death from circulatory diseases in Baotou City[J]. J Environ Occup Med, 2019, 36(4): 381-387. doi: 10.13213/j.cnki.jeom.2019.18472

    [28] 谷亚亚, 甄国新, 谈敦芳, 等. 北京市顺义区大气臭氧对居民每日死亡的影响[J]. 环境与健康杂志, 2019, 36(4): 329-334. doi: 10.16241/j.cnki.1001-5914.2019.04.011

    GU Y Y, ZHEN G X, TAN D F, et al. Effects of ambient ozone on daily mortality in Shunyi district, Beijing[J]. J Environ Health, 2019, 36(4): 329-334. doi: 10.16241/j.cnki.1001-5914.2019.04.011

    [29] 伍燕珍, 张金良. 北京市大气颗粒物与温度对人群非意外死亡影响的交互作用[J]. 环境科学研究, 2009, 22(12): 1403-1410. doi: 10.13198/j.res.2009.12.53.wuyzh.009

    WU Y Z, ZHANG J L. Interactive effects of particulate matter and temperature on population non-accidental deaths in Beijing, China[J]. Res Environ Sci, 2009, 22(12): 1403-1410. doi: 10.13198/j.res.2009.12.53.wuyzh.009

    [30]

    ZHANG Y, HUANG W, LONDON S J, et al. Ozone and daily mortality in Shanghai, China[J]. Environ Health Perspect, 2006, 114(8): 1227-1232. doi: 10.1289/ehp.9014

    [31]

    WONG C M, OU C Q, CHAN K P, et al. The effects of air pollution on mortality in socially deprived urban areas in Hong Kong, China[J]. Environ Health Perspect, 2008, 116(9): 1189-1194. doi: 10.1289/ehp.10850

    [32]

    KAN H, LONDON S J, CHEN G, et al. Season, sex, age, and education as modifiers of the effects of outdoor air pollution on daily mortality in Shanghai, China: the public health and air pollution in Asia (PAPA) study[J]. Environ Health Perspect, 2008, 116(9): 1183-1188. doi: 10.1289/ehp.10851

    [33] 贾云飞, 陈春静, 唐彦钊, 等. 南京市大气臭氧与人群循环系统疾病死亡的时间序列分析[J]. 环境与健康杂志, 2020, 37(1): 42-46. doi: 10.16241/j.cnki.1001-5914.2020.01.010

    JIA Y F, CHEN C J, TANG Y Z, et al. Association between ambient ozone and circulatory diseases mortality in Nanjing: a time-series analysis[J]. J Environ Health, 2020, 37(1): 42-46. doi: 10.16241/j.cnki.1001-5914.2020.01.010

    [34] 王欣, 张星光, 高成花, 等. 中国大气PM2.5短期暴露对心血管疾病死亡率影响的meta分析[J]. 环境与职业医学, 2021, 38(1): 17-22.

    WANG X, ZHANG X G, GAO C H, et al. Effect of short-term atmospheric PM2.5 exposure on cardiovascular mortality in China: a meta analysis[J]. J Environ Occup Med, 2021, 38(1): 17-22.

    [35]

    CHEN R, YIN P, WANG L, et al. Association between ambient temperature and mortality risk and burden: time series study in 272 main Chinese cities[J]. BMJ, 2018, 363: k4306.

    [36]

    DIMITROVA A, INGOLE V, BASAGAÑA X, et al. Association between ambient temperature and heat waves with mortality in South Asia: systematic review and meta-analysis[J]. Environ Int, 2021, 146: 106170. doi: 10.1016/j.envint.2020.106170

    [37] 马守存, 张书余, 王宝鉴, 等. 气象条件对心脑血管疾病的影响研究进展[J]. 干旱气象, 2011, 29(3): 350-354,361.

    MA S C, ZHANG S Y, WANG B J, et al. Review about the impact of meteorology condition on cardiovascular and cerebrovascular diseases[J]. J Arid Meteorol, 2011, 29(3): 350-354,361.

    [38] 班婕, 李湉湉. 北京市不同度量方式下臭氧短期暴露人群急性健康效应研究[J]. 环境与健康杂志, 2016, 33(4): 287-291. doi: 10.16241/j.cnki.1001-5914.2016.04.002

    BAN J, LI T T. Short-term effects of different ozone metrics on daily mortality in Beijing[J]. J Environ Health, 2016, 33(4): 287-291. doi: 10.16241/j.cnki.1001-5914.2016.04.002

    [39] 陈琦, 孙宏, 陈晓东, 等. 南京市臭氧短期暴露人群急性健康效应研究[J]. 江苏预防医学, 2017, 28(4): 366-368,386. doi: 10.13668/j.issn.1006-9070.2017.04.02

    CHEN Q, SUN H, CHEN X D, et al. Acute health impacts of ozone exposure on daily mortality in Nanjing[J]. Jiangsu J Prev Med, 2017, 28(4): 366-368,386. doi: 10.13668/j.issn.1006-9070.2017.04.02

    [40] 梁婷婷. 高温与臭氧对高血压病影响效应的研究[D]. 兰州: 兰州大学, 2019.

    GLIANG T T. The study of effect of high temperature and ozone on hypertension[D]. Lanzhou: Lanzhou University, 2019.

    [41] 宋全全. 高温热浪与臭氧对冠心病的交互作用影响及机制研究[D]. 兰州: 兰州大学, 2019.

    SONG Q Q. Study on the interactive effects and the related mechanism of high temperature heat wave and ozone on coronary heart disease[D]. Lanzhou: Lanzhou University, 2019.

    [42] 高豆豆. 基于代谢组学的大气PM2.5染毒大鼠肺及心血管系统损伤的机制研究[D]. 太原: 山西医科大学, 2017.

    GAO D D. Metabolomic-based mechanism research of lung and cardiovascular system injury in rats with atmospheric PM2.5 exposure[D]. Taiyuan: Shanxi Medical University, 2017.

    [43]

    NIU Y, ZHOU Y, CHEN R, et al. Long-term exposure to ozone and cardiovascular mortality in China: a nationwide cohort study[J]. Lancet Planet Health, 2022, 6(6): e496-e503. doi: 10.1016/S2542-5196(22)00093-6

    [44]

    LEI R, ZHU F, CHENG H, et al. Short-term effect of PM2.5/O3 on non-accidental and respiratory deaths in highly polluted area of China[J]. Atmos Pollut Res, 2019, 10(5): 1412-1419. doi: 10.1016/j.apr.2019.03.013

    [45]

    BELL M L, ZANOBETTI A, DOMINICI F. Evidence on vulnerability and susceptibility to health risks associated with short-term exposure to particulate matter: a systematic review and meta-analysis[J]. Am J Epidemiol, 2013, 178(6): 865-876. doi: 10.1093/aje/kwt090

    [46]

    KÜNZLI N, JERRETT M, MACK W J, et al. Ambient air pollution and atherosclerosis in Los Angeles[J]. Environ Health Perspect, 2005, 113(2): 201-206. doi: 10.1289/ehp.7523

    [47]

    WONG E M, WALBY W F, WILSON D W, et al. Ultrafine particulate matter combined with ozone exacerbates lung injury in mature adult rats with cardiovascular disease[J]. Toxicol Sci, 2018, 163(1): 140-151. doi: 10.1093/toxsci/kfy018

  • 期刊类型引用(5)

    1. 刘洁,黄俊,徐晓文,甄玲燕,陈林利,应圣洁,杜喜浩. 臭氧污染对居民心脑血管疾病死亡的影响:基于上海市闵行区. 环境与职业医学. 2024(05): 467-473 . 本站查看
    2. 宋露露,于琦,刘楠楠,高宇卉,牛泽宇,张妍,郑会秋,田家瑜,刘军霞,赵利芳,张志红. 大气细颗粒物和臭氧短期暴露与肺炎患者外周血炎性指标的相关性. 上海预防医学. 2024(06): 551-558 . 百度学术
    3. 李小强,陈小岳,陈志永,武常倩,谷令,王珂. 空气污染对高血压患者血压控制的影响及其交互作用——基于中国健康营养调查的研究. 现代预防医学. 2024(23): 4279-4285 . 百度学术
    4. 余林玲,朱昱,肖长春. 合肥市2017—2020年大气臭氧短期暴露与人群非意外死亡风险的时间序列研究. 环境与职业医学. 2024(12): 1407-1413 . 本站查看
    5. 易仁玲,陈志永,王珂,陈小岳. 苏南某市大气污染物O_3、PM_(2.5)对呼吸系统门诊量的时间序列分析. 职业与健康. 2024(24): 3416-3421 . 百度学术

    其他类型引用(0)

图(5)  /  表(6)
计量
  • 文章访问数:  161
  • HTML全文浏览量:  23
  • PDF下载量:  42
  • 被引次数: 5
出版历程
  • 收稿日期:  2022-06-19
  • 录用日期:  2023-02-05
  • 网络出版日期:  2023-04-05
  • 刊出日期:  2023-04-11

目录

/

返回文章
返回