098 CARCINO GENESIS，TERATO GENESIS&MUTA GENESIS 收稿日期：2021-08-09；修订日期：2021-03-02
基金项目：福建省自然科学基金（高校联合资助，2021J01722）
作者信息：屈莹，E-mail：*******************。*通信作者，向建军，E-mail：**********************
新冠肺炎前后
福州市区大气
颗粒物污染与居民非意外死亡的关系
屈莹1，周权2，陈一峰1，杨焱1，康淑玲1，何陈周1，姜雨1，吴传城1，刘宝英1，向建军1，*
(1．福建医科大学公共卫生学院预防医学系，福建福州350122；2．福州市疾病预防控制中心，福建福州350004)
Relationships between particulate matter pollution and
non-accidental deaths among
residents in Fuzhou city before
and after COVID-19
QU Ying1,ZHOU Quan2,CHEN Yifeng1,YANG Yan1, KANG Shuling1,HE Chenzhou1,JIANG Yu1, WU Chuancheng1,LIU Baoying1,XIANG Jianjun1,* (1.Department of Preventive Medicine,School of Public Health,Fujian
Medical University,Fuzhou350122;2.Fuzhou Center for Disease Control
and Prevention,Fuzhou350004,Fujian,China)
【摘要】目的：研究新冠肺炎前后福州市区大气颗粒物污染与居民非意外死亡的关系。方法：收集2016—2020年福州市大气污染数据、气象数据和死因监测数据，采用基于Quasi-Poisson回归广义线性模型的时间序列法，并调整长期趋势、季节性等混杂因素，分析比较新冠肺炎前后大气颗粒物-非意外死亡关系的变化。结果：新冠肺炎发生后，福州市PM2.5浓度显著降低。新冠肺炎前，PM2.5致居民非意外总死亡数的超额危险度于滞后2d内有显著性，而且在第1天效应最大，分别为[ER=1.69%，95%CI (0.79%，2.59%)]和[ER=0.93%，95%CI(0.49%，1.38%)]。而新冠肺炎后PM2.5对居民非意外总死亡数无显著影响。新冠肺炎前，PM2.5对呼吸系统疾病死亡于滞后第1天有显著影响[ER=3.01%，95%CI(0.3
5%，5.74%)]，新冠后则无显著影响。新冠肺炎前后，PM2.5浓度的降低对循环系统和心血管系统死亡风险的改变无明显影响。结论：新冠肺炎后PM2.5浓度的降低对居民非意外死亡尤其是呼吸系统疾病死亡风险的降低具有显著影响，提示大气污染治理的必要性和有效性。
【关键词】新冠肺炎；大气颗粒物；死亡率；时间序列分析全球十大污染城市
中图分类号：R122.7文献标志码：A文章编号：1004-616X(2022)02-0098-07doi：10.3969/j.issn.1004-616x.2022.02.004
【ABSTRACT】OBJECTIVE：To investigate relationships between atmospheric particulate matter pollution and non-accidental deaths among residents in Fuzhou City before and after the COVID-19pandemic.METHODS：Daily air pollution data，meteorological data，and death registry data between2016—2020were collected for Fuzhou City.Time-series analysis with a quasi-Poisson regression based generalized linear model to the relationships before and after the COVID-19pandemic after controlling for confounding factors including long-term trend，seasonality，day of week，and meteorological factors.RESULTS：After the occurrence of COVID-19，concentrations of PM2.5decreased significantly.Before the COVID-19，significant impacts of PM2.5on non-accidental death were observed between lag1-2days，with the greatest effects being lag day1，namely[ER= 1.6
9%，95%CI(0.79%，2.59%)]and[ER=0.93%，95%CI(0.49%，1.38%)]，respectively.By contrast，PM2.5had no significant impacts on non-accidental death after the COVID-19pandemic.Before the COVID-19，significant impact of PM2.5on respiratory-related death was observed on lag day1[ER=3.01%，95%CI(0.35%，5.74%)]，however，there was no significant effect after COVID-19.No noticeable changes were found on the relationships between PM2.5and circulatory-and cardiovascular-related deaths.CONCLUSION：Decreased PM2.5 concentrations after COVID-19reduced the risk of non-accidental death，particularly for respiratory-related death，indicating the necessity and effectiveness of air pollution control and prevention.
【KEY WORDS】COVID-19；atmospheric particulate；mortality rate；time-series
analysis
099
CARCINO GENESIS ，TERATO GENESIS &MUTA GENESIS
在世界卫生组织(World Health Organization ，WHO)发布的2020年全球十大健康威胁榜单中，“空气污染和气候变化”位列榜首，全球每10人中有9人会受到空气污染的影响，空气污染每年约导致全球320万人过早死亡[1]，是全球共同面对的重要公共卫生问题。而在大气污染物的组分中，颗粒物对人健康的影响超过臭氧、二氧化硫等其他气态污染物。美国的Apte 等[2]对全球185个国家的颗粒物污染相关疾病负担进行评估，发现2016年大气细颗粒物(atmospheric particulat matter 2.5，PM 2.5)暴露使得全球预期平均寿命缩短近1年，而对于空气污染较为严重的亚洲和非洲，平均寿命缩短1.2~1.9年；若全球PM 2.5均低于WHO 规定的限制标准(10μg/m 3)，全球人均预期寿命可增加近0.59年，堪比肺癌和乳腺癌彻底根除。因此，减少PM 2.5污染和暴露的公共卫生意义重大。
国内外大量流行病学研究表明，PM 2.5长期暴露不但会增加呼吸系统、心脑血管系统(如冠心病、心肌梗死、心衰、心律失常、中风)等疾病的发病率和死亡率[3-6]，而且会影响母婴健康[7-8]，也与癌症、认知功能降低、认知迟缓、阿尔兹海默氏症、帕金森氏病等存在关联[9-12]。Ferreira 等[13]2016年对巴西
某城市的研究发现，大气PM 2.5和PM 10浓度每增加10μg/m 3，相应的住院量分别增加23.5%和12.8%；尤其是循环系统相关疾病，PM 2.5浓度每增加10μg/m 3可导致住院量增加19.6%。Jo 等[14]对韩国春川县的研究发现，随着PM 2.5浓度的增加，COPD 相关的医院就诊人数也增加。
国内复旦大学公共卫生学院的Liu 等[3]分析了22个国家或地区的652个城市的大气颗粒物污染与逐日死亡率之间的关系，结果显示短期颗粒物暴露即可增加居民死亡的风险。陈仁杰等[15]分析国内272个城市PM 2.5-10与非意外死亡和心肺等系统疾病死亡的关系，在调整其他气态污染物后，关联仍具有统计学意义；而且，不同城市的PM 2.5-10与死亡风险的关联因各城市污染物水平不同而存在差异。
新冠肺炎暴发流行期间，因疫情防控影响社会经济生产和生活方式，空气监测数据显示全国88％的站点PM 2.5浓度平均降低约12.6μg/m 3，PM 10与PM 2.5变化一致[16]。新冠肺炎前后大气颗粒物浓度的改变，为当前和今后大气污染治理的健康效应评估提供良好机会。因此，我们收集2016—2020年福州市污染数据、气象数据、死亡数据，采用基于Quasi-Poisson 回归广义线型模型的时间序列法，并调整长期趋势、季节性等混杂因素，研究新冠肺炎前后大气颗粒物污染与居民非意外死亡的关系是否有改变。
1
材料与方法
1.1
研究对象
收集2016年1月1日—2020年12月31日期间福州市气象数据、大气污染物数据以及居民死因监测数据，以上监测数据分别由福州市气象局和疾病预防控制中心提供。死因数据以国际疾病分类10(ICD-10)进行编码，本研究分析的死因类别包括非意外事件死亡(A00-R99)、呼吸系统疾病死亡(J00-J98)和循环系统疾病死亡(I00-I99)，循环系统疾病死亡中又观察对心血管系统疾病死亡的影响(I05-09，I11，I20-27，I30-52)。同期福州市大气污染物数据包括PM 2.5、PM 10、NO 2、SO 2、CO 、O 3,1h 及O 3,8h 的浓度等7项指标，单位为μg/m 3。
1.2分析方法
采用基于Quasi-Poisson 回归的广义线性模型(generalized linear model ，GLM)分析福州市新冠肺炎前后每日非意外总死亡数、呼吸系统、循环系统及心血管系统死亡数与大气颗粒污染物的关系。新冠肺炎前后的时间划分点为2020年1月22日，即福州市出现第1例新冠肺炎患者。本研究中，新冠肺炎前是指2016年1月1日—2020年1月22日，新冠肺炎后指2020年1月23日—12月31日。自变量为PM 2.5日均浓度，因变量为每日非意外总死亡、呼吸系统、循环系统及心血管系统的死亡数。采用自然样条平滑函数(ns)拟合时间序列的长期趋势和季节趋势，自由度选择7/年。控制时间长期趋势、气象因素、星期几效应(DOW )、节假日效应(Holiday )等混杂因素。自由度的选择基于赤池信息法则(AIC )确定，DOW 、Holiday 效应以分类变量的形式纳入模型。
GLM 模型公式为：log E (Yt )=βZt +ns(Time ，7×3)+ns(Temp ，
df )+ns(RH ，df )+ns(Pressure df )+as.factor(DOW )+as.factor(Holiday )+Intercept
根据模型分析结果的显著性，分析当日(lag0)的PM 2.5浓度与居民非意外总死亡、呼吸系统、循环系统及心血管系统死亡数，然后再进行滞后效应分析，并且对性别及年龄进行分层分析。根据历史文献，滞后效应区间选为滞后第1~7天(lag1~lag7)。计算PM 2.5浓度每升高10μg/m 3，死亡发生的超额危险度(excess risk ，ER)及其95%CI(confidence interval)，采用最大效应值作为颗粒污染物对居民非意外死亡效应的暴露风险估计值。
利用SPSS 20.0统计软件分析新冠肺炎前后的气象因素和大气污染物浓度变化，用x ˉ±s 、最小值(Min)、最大值(Max)、四分位距(IQR)和百分位数(P 25
、
100
CARCINO GENESIS ，TERATO GENESIS &MUTA GENESIS
P 50、P 75)进行统计描述。福州市区居民逐日非意外死亡数也采用上述指标进行统计描述。大气污染物与福州市区居民每日非意外总死亡、呼吸系统、循环系统及心血管系统死亡数的相关性分析采用Spearman 相关分析。新冠肺炎前后气象、污染物浓度的改变采用t 检验。应用R 4.0.3分析软件作新冠肺炎前后福州市区PM 2.5暴露对居民非意外总死亡、呼吸系统、循环系统及心血管系统死亡效应的暴露风险的时间序列分析，并进行年龄、性别的分层分析。
2
结果
2.1
新冠肺炎前后福州市区大气污染物浓度的改变新冠肺炎前，福州市区每日PM 2.5、PM 10、NO 2、SO 2、CO 、O 3,1h 、O 3,8h 浓度均值分别为26.17、60.10、28.26、5.81、0.70、103.26、86.82μg/m 3(见表1)，按照《环境空气质量标准》(GB 3095-2012)，新冠肺炎前福州市PM 2.5、PM 10年均浓度均符合国家二级标准(标准限值分别为35、70μg/m 3)，NO 2、SO 2、CO 、O 3,1h 及O 3,8h 年均浓度均符合国家一级标准(标准限值分别为40、20、4、160、100μg/m 3)。
新冠肺炎后，福州市区每日PM 2.5、PM 10、NO 2、SO 2、CO 、O 3,1h 、O 3,8h 浓度均值分别为19.35、47.60、
21.29、4.55、0.66、106.52、88.42μg/m 3(见表1)，按照《环境空气质量标准》(GB 3095-2012)，新冠肺炎后福州市PM 2.5、PM 10年均浓度均符合国家二级标准(标准限值分别为35、70μg/m 3)，NO 2、SO 2、CO 、O 3,1h 及O 3,8h 年均浓度均符合国家一级标准(标准限值分别为40、20、4、160、100μg/m 3)。
2.2新冠肺炎前后福州市区非意外总死亡及其他系统疾病死亡的描述性分析
研究期间福州市区每日报告非意外死亡人数基本比较稳定，呈现一定的季节波动，冬春季死亡人数相对较多，且每日因呼吸、循环、心血管系统疾病死亡人数波动模式与非意外总死亡人数基本一致(见图1)。
新冠肺炎前(2016年1月1日—2020年1月22日)，福州市区非意外总死亡人数为44787人，其中男性为25290人，占非意外总死亡人数的56.5%；≥65岁年龄组死亡人数为35770人，占非意外总死亡人数的79.8%。新冠肺炎后(2020年1月23日—12月31日)，福州市区非意外总死亡人数为10927人，其中男性为6205人，占非意外总死亡人数的56.8%；≥65岁年龄组死亡人数为9002人，占非意外总死亡人数的82.4%。二者在性别、年龄分层上具有类似的分布趋势。在呼吸、循环、心血管系统疾病死亡人数方面，男性和≥65岁年龄组也表现出明显的特征(见表2)。
3,1h 3,8h 新冠肺炎前，其差异均具有统计学意义(P <0.05或
0.01).
图1
新冠肺炎前后福州市逐日非意外总死亡人数散点图
101
CARCINO GENESIS ，TERATO GENESIS &MUTA GENESIS
2.3
相关性分析
新冠肺炎前，福州市区空气污染物与居民逐日非意外总死亡、呼吸系统、循环系统及心血管系统疾病死亡数的Spearman 相关性分析(见表3)表明，非意外总死亡数与PM 2.5、PM 10、NO 2、CO 日均浓度均呈正相关，其中与PM 2.5的相关性系数最大(r =0.220)；与O 3,1h 、O 3,8h 日均浓度均呈负相关，且O 3,1h 的相关性稍大(r =-0.132)。呼吸系统死亡数与PM 2.5、NO 2、SO 2、CO 日均浓度均呈正相关，其中与NO 2的相关性系数最大(r =0.105)。循环系统死亡数与PM 2.5、PM 10、NO 2、SO 2、CO 日均浓度均呈正相关，其中与NO 2的相关性系数最大(r=0.211)；与O 3,1h 和O 3,8h 日均浓度均呈负相关，其中与O 3,1h 的相关性稍大(r=-0.153)。心血管系统死亡数与PM 2.5、PM 10、NO 2、SO 2、CO 日均浓度均呈正相关，其中与NO 2的相关性系数最大(r =0.135)；与O 3,1h 和O 3,8h 日均浓度均呈负相关，其中与O 3,1h 的相关性稍大(r =-0.096)。
新冠肺炎后，福州市区空气污染物与居民逐日非意外总死亡、呼吸系统、循环系统及心血管系统疾病死亡数的Spearman 相关性分析(见表4)表明，非意外总
死亡数与PM 2.5、SO 2日均浓度均呈正相关，而与SO 2的相关性相对较强(r=0.180)。呼吸系统死亡数与O 3,1h 日均浓度均呈负相关(r =-0.225)。心血管系统死亡数与O 3,8h 日均浓度均呈负相关(r =-0.135)。
2.4时间序列分析
2.4.1新冠肺炎前后福州市PM 2.5暴露与居民非意外总死亡风险及分层分析新冠肺炎前，福州市PM 2.5对居民非意外总死亡的超额风险度(ER 值)存在两天的滞后效应，且于滞后第1天达到最大值[ER=1.69%，95%CI(0.79%，2.59%)]。分层分析结果(见表5)显示：PM 2.5致男性居民非意外总死亡的ER 值也存在两天的滞后效应，于第1天达到最大值[ER=2.18%，95%CI (1.00%，
3.38%)]；女性居民非意外总死亡的ER 值先升高后降低再升高，于滞后第1天达到最大值[ER=1.07%，95%CI(-0.17%，2.32%)]，但不具有统计学意义；<65岁年龄组非意外总死亡数的ER 值先升高后降低，于滞后第4天达到最大值[ER=0.98%，95%CI (-0.84%，2.83%)]，但不具有统计学意义；≥65岁年龄组非意外总死亡数的ER 值先升高后降低，于滞后第1天达到最大值[ER=2.07%，95%CI(1.09%，3.07%)]。
而新冠肺炎后，福州市PM 2.5浓度对居民总非意外死亡及性别和年龄分层分析结果(见表6)均无显著性影响。比较新冠肺炎前后福州市PM 2.5对居民非意外总死亡风险可知，PM 2.5对死亡影响的滞后效应趋势基本一致，但新冠肺炎后死亡风险低于新冠肺炎前。
2.4.2新冠肺炎前后福州市PM 2.5暴露与居民呼吸系统、循环系统及心血管系统死亡超额风险度及性别年龄分层分析结果新冠肺炎前，福州市PM 2.5对居民呼吸系统疾病死亡的滞后效应先升高后降低，并于滞后第1天达到最大值(ER=
3.01%，95%CI ：0.35%~5.74%)。性别年龄分层分析结果(见表7)显示：PM 2.5对
102
CARCINO GENESIS ，TERATO GENESIS &MUTA GENESIS
男性居民呼吸系统疾病死亡风险无显著影响；PM 2.5对女性居民呼吸系统疾病死亡风险的影响先升高后降低再升高，于滞后第1天达到最大值[ER=3.90%，95%CI (0.32%，8.30%)]，但不具有统计学意义；<65岁年龄组呼吸系统疾病死亡数的ER 值先升高后降低，于滞后第4天达到最大值[ER=5.30%，95%CI(-2.86%，14.14%)]，但均不具有统计学意义；≥65岁年龄组呼吸系统疾病死亡数的ER 值先升高后降低，于滞后第1
天达到最大值[ER=2.92%，95%CI(0.11%，5.80%)]。
新冠肺炎后，福州市PM 2.5对居民呼吸系统疾病死亡风险影响无显著影响，性别、年龄分层分析结果也均无统计学意义(见表8)。新冠肺炎前后福州市PM 2.5对居民循环系统和心血管系统死亡风险均有显著影响，也即是新冠肺炎后PM 2.5浓度的降低短期内未能降低循环系统和心血管系统疾病死亡风险(因篇幅原因，数据未呈现)。
3讨论
本研究旨在探究新冠肺炎前后福州市区大气颗粒物污染与居民非意外死亡关系的改变。此前已有大量研究评估了大气颗粒物污染与居民非意外死亡的关照系[8，11，17]，突如其来的新冠肺炎疫情，让我们有一个难得的机会来评估在同一地区、同一人、相同经济水平下空气污染物浓度改变对居民非意外总死亡及呼吸系统等其他各系统疾病死亡的影响，为大气污染的防治提供参考依据。福州市新冠肺炎前逐日PM 2.5、
PM 10、浓度均值分别为26.17、60.10μg/m 3，而新冠肺炎后逐日PM 2.5、PM 10浓度均值分别为19.35、47.60μg/m 3，新冠肺炎后明显低于新冠肺炎前，这是由于新冠期间交通和生产活动的减少，与新冠肺炎疫情期间中国其他地区大气污染物变化基本一致[18]。我们对数据分析的结果表明，新冠肺炎后PM 2.5浓度的降低对居民非意外总死亡尤其是呼吸系统系统疾病死亡风险的降低具有显著影响，但新
冠肺炎后PM 2.5浓度降低对循
环系统和心血管系统疾病死亡风险的降低无明显影
版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。
发表评论