• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    環(huán)境條件對遵義市呼吸系統(tǒng)疾病的影響及預(yù)測研究

    2018-11-28 03:47:06王式功謝佳君尚可政
    中國環(huán)境科學(xué) 2018年11期
    關(guān)鍵詞:滯后效應(yīng)遵義市氣象

    樂 滿,王式功,*,謝佳君,馬 盼,尚可政

    ?

    環(huán)境條件對遵義市呼吸系統(tǒng)疾病的影響及預(yù)測研究

    樂 滿1,王式功2,1*,謝佳君3,馬 盼2,尚可政1

    (1.蘭州大學(xué)大氣科學(xué)學(xué)院,氣象環(huán)境與人體健康研究中心,甘肅 蘭州 730000;2.成都信息工程大學(xué)大氣科學(xué)學(xué)院,高原大氣與環(huán)境四川省重點(diǎn)實(shí)驗(yàn)室,四川 成都 610225;3.貴州省遵義市氣象局,貴州 遵義 563000)

    為探究氣象和空氣污染等環(huán)境條件對呼吸系統(tǒng)疾病發(fā)病的影響,為遵義地區(qū)相關(guān)疾病預(yù)防提供科學(xué)依據(jù),采用分布滯后非線性模型和廣義線性、相加模型,利用當(dāng)?shù)貧庀蠛臀廴举Y料,分析了2012~2016年遵義市氣象環(huán)境要素對呼吸系統(tǒng)疾病發(fā)病的影響.結(jié)果顯示,遵義市呼吸系統(tǒng)疾病發(fā)病狀況與當(dāng)?shù)亻L期氣候狀態(tài)基本保持一致,氣候效應(yīng)對其的影響占主導(dǎo)作用,其中,冬季為呼吸系統(tǒng)疾病高發(fā)期,立秋至處暑前后其發(fā)病人數(shù)最少,表明此時(shí)間段內(nèi)當(dāng)?shù)貧夂驐l件對呼吸系統(tǒng)疾病患者有氣候療養(yǎng)效應(yīng).氣溫對呼吸系統(tǒng)疾病發(fā)病的影響以低溫滯后效應(yīng)為主,在其敏感閾值附近氣溫每變化1℃,發(fā)病人數(shù)將累積增加31.6%(95%CI:4.4%~65.8%);氣壓以高壓滯后效應(yīng)為主,相對濕度則在低濕部分同時(shí)有即時(shí)和滯后效應(yīng).舒適度對呼吸系統(tǒng)疾病的影響,在冷、熱斜脅迫下其發(fā)病人數(shù)明顯多于舒適狀況時(shí). PM2.5、SO2和NO2三種污染物的影響都以即時(shí)效應(yīng)為主,而CO則在累積滯后lag04時(shí)相對危險(xiǎn)度最高,PM2.5與呼吸系統(tǒng)疾病發(fā)病人數(shù)的暴露-反應(yīng)曲線呈單調(diào)線性分布,SO2、NO2和CO均為“J”型分布.低溫與高濃度NO2或者低濕與高濃度SO2的協(xié)同作用對呼吸系統(tǒng)疾病的影響較大.建立的全年和季節(jié)多元逐步回歸方程的試預(yù)報(bào)準(zhǔn)確率在75%以上(夏季除外),其中分季節(jié)建模預(yù)測效果顯著優(yōu)于全年預(yù)測效果.

    環(huán)境條件;氣象因子;空氣污染;呼吸系統(tǒng)疾病

    研究表明,氣象條件的變化和空氣污染的日益加劇與人們的健康息息相關(guān)[1-12].全球氣候變化對人類健康產(chǎn)生的直接或間接的影響進(jìn)一步加劇[13-14],導(dǎo)致相關(guān)天氣敏感性疾病發(fā)病率、乃至超額死亡率增加,如高溫?zé)崂说冉陙碓斐纱罅咳藬?shù)死亡[15-18].同時(shí),空氣污染對人類健康的影響[19-22]也越來越受到國際的關(guān)注,研究發(fā)現(xiàn)[23-24],受大氣污染的影響,呼吸系統(tǒng)疾病、心腦血管疾病及癌癥的發(fā)病率及死亡率逐年升高.但是多數(shù)研究主要著眼于易出現(xiàn)較極端環(huán)境條件的區(qū)域如京津冀或南京等展開天氣敏感性疾病的研究,對于氣候環(huán)境條件較溫和區(qū)域的不良健康影響評估卻較為缺失.因此,文章以遵義市作為研究目標(biāo)區(qū)域,探索氣象因子和空氣污染等環(huán)境條件對呼吸系統(tǒng)疾病發(fā)病與流行的影響,并建立預(yù)報(bào)模型,不僅能夠?yàn)榇龠M(jìn)當(dāng)?shù)睾侠淼拈_展健康氣候療養(yǎng)提供科學(xué)依據(jù),也可為這些區(qū)域更好的開展應(yīng)對氣候變化和空氣污染提供科學(xué)參考,對目前醫(yī)學(xué)氣象學(xué)的研究進(jìn)行補(bǔ)充.

    1 資料與方法

    1.1 資料

    氣象資料:遵義市12個(gè)觀測臺站的2012~2016年的地面日均氣象資料,包括氣溫、氣壓、相對濕度和風(fēng)速等氣象要素.

    污染資料:遵義市五個(gè)污染觀測站點(diǎn)的2014年1月1日~2016年12月31日的逐時(shí)觀測資料以及由國家環(huán)境保護(hù)部網(wǎng)站(http://www.zhb. gov.cn)提供的日均資料,主要包括AQI指數(shù)、首要污染物、PM2.5、PM10、SO2、CO、NO2等污染因子變量.

    疾病資料:遵義市2012年1月1日~2016年12月31日播州區(qū)各醫(yī)院就診病案記錄,其中包含2家省直屬三級甲等綜合醫(yī)院和4家貴州省二級甲等綜合醫(yī)院,病案記錄中主要包括患者家庭住址、就診日期、患者年齡、性別、疾病診斷和就診醫(yī)院等基本信息.根據(jù)國際疾病分類標(biāo)準(zhǔn)第10版(ICD-10)[25],共挑選出呼吸系統(tǒng)疾病就診病例資料264347例,并對呼吸系統(tǒng)(ICD-10編碼:J00-J99)疾病資料進(jìn)行分類整理,分別挑選上呼吸道感染和下呼吸道感染兩類疾病資料進(jìn)行分析.

    氣象資料和大氣污染物濃度資料均來自正規(guī)網(wǎng)站系統(tǒng),并有專人負(fù)責(zé)檢查、判別及管理,最終的資料無缺測.疾病資料涵蓋范圍廣,且與京津冀和北上廣等區(qū)域多流動(dòng)人口不同,遵義市研究對象基本為常駐人口,因而疾病資料對當(dāng)?shù)囟愿叽硇?并剔除了年齡、性別等信息不全和非遵義市居住的人群疾病資料,從而保證研究資料的準(zhǔn)確性和可靠性.

    1.2 方法

    文章利用DLNM(分布滯后非線性模型)、GLM(廣義線性模型)和GAM(廣義相加模型)進(jìn)行分析.其中DLNM模型可靈活地描述氣溫健康效應(yīng)的空間和其發(fā)生的滯后維度同時(shí)變化的影響[26-28],其核心是二維空間交叉基的構(gòu)建,而GLM模型[29]主要假設(shè)連接函數(shù)與預(yù)測變量之間的關(guān)系為線性關(guān)系.GAM模型則是在沒有選擇特定響應(yīng)函數(shù)的情況下使用.在控制了時(shí)間趨勢、星期幾效應(yīng)和節(jié)假日效應(yīng)的前提下:選擇相對危險(xiǎn)度(RR)作為主要的表征量,利用DLNM和GLM模型進(jìn)行氣象因子與相關(guān)天氣敏感性急病發(fā)病的關(guān)系研究,GAM模型則進(jìn)行空氣污染對呼吸系統(tǒng)疾病的影響及環(huán)境因子之間的協(xié)同效應(yīng)分析.以上模型參數(shù)自由度的選擇采用赤池信息準(zhǔn)則(Akaike’s lnformation Criterion,AIC).

    文章在構(gòu)建呼吸系統(tǒng)疾病預(yù)測模型、進(jìn)行預(yù)測模型優(yōu)度檢驗(yàn)時(shí),分別選取平均絕對誤差(MAE)、均方根誤差(RMSE)、平均絕對百分誤差(MAPE)、預(yù)測準(zhǔn)確率()四種參數(shù)進(jìn)行驗(yàn)證,其中和表示實(shí)測值和預(yù)測值,表示樣本數(shù)量,計(jì)算方法如下:

    主要利用R語言編程軟件和SPSS軟件進(jìn)行分析研究.

    2 結(jié)果與分析

    2.1 呼吸系統(tǒng)疾病發(fā)病人數(shù)和環(huán)境要素的描述性統(tǒng)計(jì)及相關(guān)性分析

    遵義市屬亞熱帶季風(fēng)性氣候,地形起伏大,地貌類型復(fù)雜,其平均氣溫較為舒適為16.1℃,平均相對濕度接近80%,平均風(fēng)速較小.整體上空氣質(zhì)量較好,各污染物濃度較低,PM2.5、PM10、SO2和NO2的年均濃度分別為47.4, 78.0, 17.1和30.4μg/m3(表1).

    表1 呼吸系統(tǒng)疾病就診人數(shù)與主要環(huán)境要素的描述性分析

    圖1 遵義市呼吸系統(tǒng)疾病發(fā)病人數(shù)的不同時(shí)間尺度分布

    圖1為不同時(shí)間尺度遵義市總呼吸系統(tǒng)、上呼吸道感染和下呼吸道感染(簡稱“上感”和“下感”,下同)發(fā)病人數(shù)的分布特征.從季節(jié)分布情況看,冬春季為遵義市呼吸系統(tǒng)疾病高發(fā)期,與多數(shù)北方城市呼吸系統(tǒng)疾病高發(fā)期在秋冬季有所不同,且在12、1、2和3月較多,8月份最少為102.3人/d,下呼吸道感染發(fā)病人數(shù)變化特征與之一致.上呼吸道感染的就診人數(shù)在夏季和冬季最多,需引起注意.

    我國24節(jié)氣的劃分充分考慮了季節(jié)、氣候和物候等自然現(xiàn)象的變化,在注重氣象因素的同時(shí)也考慮了氣候因素,對于人們根據(jù)中醫(yī)理論進(jìn)行養(yǎng)生防病具有更好的指導(dǎo)意義.圖1(b)為24節(jié)氣呼吸系統(tǒng)疾病發(fā)病人數(shù)分布(圖中黑色曲線為不同節(jié)氣的氣溫分布,黑色豎實(shí)線為發(fā)病人數(shù)年均值):呼吸系統(tǒng)疾病發(fā)病第一高峰期出現(xiàn)在隆冬季節(jié)的大雪至小寒期間,次峰期出現(xiàn)在冬末春初的雨水到驚蟄期間;而夏末的立秋至處暑前后遵義市呼吸系統(tǒng)疾病發(fā)病人數(shù)最少(102.7人/d),表明此時(shí)間段內(nèi)當(dāng)?shù)貧夂驐l件對呼吸系統(tǒng)疾病患者有氣候療養(yǎng)效應(yīng).與24節(jié)氣的氣溫分布對比,發(fā)現(xiàn)呼吸系統(tǒng)疾病發(fā)病人數(shù)與氣溫呈現(xiàn)相同周期的反相位變化,可見遵義市呼吸系統(tǒng)疾病發(fā)病主要與當(dāng)?shù)亻L期氣候狀態(tài)保持一致,氣候效應(yīng)對其的影響占主導(dǎo)作用,注意到上呼吸道感染24節(jié)氣發(fā)病人數(shù)主要夏至?xí)r達(dá)到最多(56.3人/d),在秋冬之后發(fā)病人數(shù)也逐漸增加,可見上呼吸道感染同時(shí)受到外界高溫、低溫和粉塵的顯著影響.

    表2 呼吸系統(tǒng)疾病發(fā)病人數(shù)與主要環(huán)境因子的Spearman相關(guān)

    注:“**”和“*”分別表示通過了顯著性水平=0.01和=0.05的顯著性檢驗(yàn),下同.

    表2為環(huán)境因子與遵義市呼吸系統(tǒng)疾病發(fā)病人數(shù)的Spearman相關(guān)表,污染因子變量主要與呼吸系統(tǒng)疾病發(fā)病人數(shù)呈顯著正相關(guān),表明污染物濃度的增加會(huì)明顯導(dǎo)致呼吸系統(tǒng)發(fā)病人數(shù)的上升;總呼吸和下感疾病發(fā)病人數(shù)與氣壓呈顯著正相關(guān),與氣溫、風(fēng)速等均呈顯著負(fù)相關(guān),上感發(fā)病人數(shù)與氣溫和水汽壓呈正相關(guān).

    2.2 氣象因子對呼吸系統(tǒng)疾病發(fā)病的影響與滯后性

    2.2.1 氣溫對呼吸系統(tǒng)疾病發(fā)病的影響 考慮到疾病對氣象要素的敏感性存在滯后效應(yīng),且一般一個(gè)完整的天氣過程時(shí)長約為1~2周左右,故而考慮最長累積滯后時(shí)長為14d,即共考慮最長累積滯后時(shí)長為15d,以平均氣溫的中位數(shù)17.1℃作為模型參考溫度,并以此為分界點(diǎn)(圖2中的中間豎虛線)分別討論溫度小于和大于該溫度時(shí)的冷效應(yīng)和熱效應(yīng).

    圖2中灰色區(qū)域?yàn)?5%的置信區(qū)間,兩側(cè)豎虛線分別為各氣溫指標(biāo)的1%和99%分位數(shù).氣溫對遵義市總呼吸系統(tǒng)疾病發(fā)病的累積滯后效應(yīng)相似:整體上呈“反J”型分布,主要以冷效應(yīng)為主,當(dāng)氣溫持續(xù)下降低于7℃左右時(shí),RR幾乎以指數(shù)形式增長迅速增大,對呼吸系統(tǒng)疾病發(fā)病的危害顯著增加;熱效應(yīng)部分對遵義市呼吸系統(tǒng)疾病發(fā)病的影響不大,其RR值接近于1.

    圖3 不同氣溫對呼吸系統(tǒng)疾病發(fā)病的滯后效應(yīng)

    為探究高溫和低溫對呼吸系統(tǒng)疾病發(fā)病的滯后效應(yīng),分別選取日最高氣溫的27℃、30℃、33℃和最低氣溫-1℃、5℃和10℃分析其滯后影響(圖3).高溫影響主要表現(xiàn)出即時(shí)效應(yīng),在累積滯后1d時(shí)對呼吸系統(tǒng)疾病發(fā)病有顯著影響;溫度偏低且約在lag=2~3時(shí)對呼吸系統(tǒng)疾病的危害最大,尤其是氣溫在0℃以下時(shí),滯后效應(yīng)可以累積到14d左右,溫度越低其危害越大.同時(shí),從頻數(shù)分布(圖2)可以看到,遵義市氣溫主要集中于5~27℃之間,日最高氣溫高于32℃和最低氣溫低于0℃的占比均很小.因此,接下來的分析主要以平均氣溫為代表分析溫度對不同性別人群發(fā)病的影響.

    不同氣溫段的時(shí)間滯后差異明顯(圖4(a)):平均氣溫對總呼吸系統(tǒng)疾病發(fā)病的影響以低溫滯后效應(yīng)為主,且低溫對女性患者的影響主要在lag=7時(shí)達(dá)到最大,同時(shí)女性患者對低溫的影響敏感性略低于男性患者.在滯后0d時(shí)(圖4(b))平均氣溫對呼吸系統(tǒng)疾病發(fā)病的影響主要以高溫即時(shí)效應(yīng)為主,低溫效應(yīng)不顯著,當(dāng)lag=2時(shí),低溫效應(yīng)開始凸顯,且隨著滯后時(shí)間的延長其危險(xiǎn)度進(jìn)一步增加,對呼吸系統(tǒng)疾病發(fā)病有顯著影響的低溫溫度也逐漸上升,因此遵義市低溫滯后的持續(xù)效應(yīng)對呼吸系統(tǒng)疾病發(fā)病影響持續(xù)時(shí)間較長.

    表3 低溫和高溫對呼吸系統(tǒng)疾病發(fā)病人數(shù)的影響統(tǒng)計(jì)表

    注:表中加粗部分表明通過了=0.05的顯著性檢驗(yàn),下同.

    表3為平均氣溫的1%和99%分位數(shù)(1.1℃和28.8℃)分別滯后09d和0d RR統(tǒng)計(jì)表,以研究其低溫滯后效應(yīng)和高溫即時(shí)效應(yīng)對呼吸系統(tǒng)疾病發(fā)病的影響.遵義市總呼吸系統(tǒng)疾病受低溫滯后效應(yīng)顯著,低溫每變化1℃,相應(yīng)的發(fā)病人數(shù)將累積增加31.6%(95%CI: 4.4%~65.8%).上感發(fā)病情況則同時(shí)受到低溫滯后和高溫即時(shí)效應(yīng)的顯著影響,高溫每升高1℃,上感發(fā)病人數(shù)將增加10.6%(95%CI: 0.9%~21.2%).此外,女性患者對低溫累積滯后效應(yīng)的響應(yīng)比男性的更加敏感,需引起注意.

    2.2.2 氣壓、風(fēng)速和相對濕度對呼吸系統(tǒng)疾病發(fā)病的影響 選取平均氣壓、平均風(fēng)速和相對濕度的中位數(shù)925hPa、1.2m/s和80.0%為相應(yīng)的參照標(biāo)準(zhǔn)分析對呼吸系統(tǒng)疾病發(fā)病的滯后效應(yīng).圖5為平均氣壓、風(fēng)速和相對濕度的滯后效應(yīng)RR曲面圖:平均氣壓主要表現(xiàn)為高壓效應(yīng),低壓對呼吸系統(tǒng)疾病發(fā)病的影響不顯著;同時(shí),風(fēng)速極小和極大時(shí)均表現(xiàn)為一定程度的即時(shí)效應(yīng),說明適當(dāng)風(fēng)速對呼吸系統(tǒng)疾病發(fā)病影響不大;相對濕度主要表現(xiàn)低濕效應(yīng),且存在兩個(gè)滯后危險(xiǎn)度峰值,第一峰值在滯后2~3d左右出現(xiàn),第二滯后峰值則出現(xiàn)在7~14d左右,此時(shí)主要是相對濕度的變化會(huì)對病毒或細(xì)菌等的繁殖產(chǎn)生影響從而進(jìn)一步誘發(fā)呼吸系統(tǒng)疾病發(fā)病,需引起重視,在高濕條件下,對遵義市呼吸系統(tǒng)疾病危害較小.

    圖6 平均氣壓、風(fēng)速和相對濕度特定值的RR值隨滯后時(shí)長的變化

    根據(jù)以上分析,分別選擇氣壓939.3hPa(99%分位數(shù))、風(fēng)速2m/s(95%分位數(shù))和相對濕度56%(1%分位數(shù))研究其對呼吸系統(tǒng)疾病發(fā)病的影響(圖6):高壓主要在滯后3d時(shí)對呼吸系統(tǒng)疾病發(fā)病有明顯促進(jìn)作用,且在lag=6~7d后危險(xiǎn)度最高;風(fēng)速的滯后效應(yīng)仍然不顯著,相對濕度的滯后效應(yīng)隨時(shí)間呈現(xiàn)“V”型變化,低濕天氣對呼吸系統(tǒng)疾病的影響同時(shí)具有即時(shí)和滯后效應(yīng).

    根據(jù)圖6結(jié)果,計(jì)算各因子的相應(yīng)累積滯后效應(yīng)(表4):風(fēng)速對呼吸系統(tǒng)疾病發(fā)病的效應(yīng)整體不顯著.高壓低濕條件對女性群體有顯著影響,氣壓每變化1hPa女性發(fā)病人數(shù)將累積增加20.0%(95%CI: 5%~43.3%),相對濕度每減小一個(gè)單位,女性患者將累積增加20.9%(95%CI:6.9%~37.0%),這兩種因子對男性人群沒有顯著影響.上感和下感分別與高壓和低濕條件有顯著關(guān)系,此種氣象條件下遵義市上呼吸道感染發(fā)病人數(shù)將增加34.8%(95%CI:8.7%~ 66.9%),下感發(fā)病人數(shù)將增加13.6%(95%CI:0.1%~ 29.4%).

    表4 遵義市高壓、風(fēng)速和低濕天氣對呼吸系統(tǒng)疾病發(fā)病的影響

    表5 呼吸系統(tǒng)疾病發(fā)病人數(shù)與體感溫度(Tg)的Spearman相關(guān)

    圖7 不同分位數(shù)體感溫度滯后影響曲線(左)和不同舒適度等級[31]發(fā)病人數(shù)的分布(右)

    2.2.3 舒適度對呼吸系統(tǒng)疾病發(fā)病的影響 多項(xiàng)研究表明呼吸系統(tǒng)疾病往往受多種氣象因子的綜合作用,因此選擇綜合考慮了溫度、濕度等多種氣象要素、基于黃金分割率的舒適度指數(shù)[30-31]分析呼吸系統(tǒng)疾病發(fā)病對氣象要素變化綜合效應(yīng)的響應(yīng)情況.表5中體感溫度與呼吸系統(tǒng)疾病發(fā)病具有顯著相關(guān)性.

    選取體感溫度的1%,5%和10%百分位數(shù)研究其冷效應(yīng),選取90%,95%和99%百分位數(shù)研究其熱效應(yīng),結(jié)果表明(圖7左)低體感溫度對呼吸疾病發(fā)病的滯后效應(yīng)顯著,體感溫度越低,越容易引發(fā)呼吸系統(tǒng)疾病疾病;高體感溫度對呼吸系統(tǒng)疾病的影響主要為即時(shí)效應(yīng),隨滯后時(shí)間增加其影響逐漸減小.整體上冷、熱斜脅迫下呼吸系統(tǒng)疾病發(fā)病人數(shù)明顯多于舒適情況(圖7右).上感發(fā)病人數(shù)在熱不舒適等級(微熱、炎熱、熱和酷熱等級,下同)約占48.1%,遠(yuǎn)遠(yuǎn)大于舒適和冷不舒適等級各自所占比例,下感發(fā)病人數(shù)在冷不舒適等級所占比例遠(yuǎn)超舒適和熱不舒適等級,約為53.3%.

    2.3 污染物濃度對呼吸系統(tǒng)疾病發(fā)病的影響與滯后性

    研究表明,污染物對人體健康的影響存在一定的滯后和累積效應(yīng)[32-33].經(jīng)過模型調(diào)試,GAM模型具有較好的描述效果,因此本部分主要利用GAM模型探討不同污染物濃度對呼吸系統(tǒng)疾病發(fā)病的滯后效應(yīng)(圖8),圖中各污染物濃度分別增加相應(yīng)單位濃度時(shí),呼吸系統(tǒng)疾病發(fā)病總?cè)藬?shù)增加的ER%值及其95%CI置信區(qū)間,均通過=0.01的顯著性檢驗(yàn),有統(tǒng)計(jì)學(xué)意義.

    當(dāng)遵義市PM2.5、SO2和NO2增加相應(yīng)濃度10μg/m3時(shí),均在滯后0d (lag0)與呼吸系統(tǒng)疾病就診人數(shù)關(guān)聯(lián)的ER%值最大,然后逐漸減小,以即時(shí)效應(yīng)為主,這是因?yàn)樽窳x市整體上空氣質(zhì)量優(yōu)良率高,當(dāng)?shù)厝巳簩諝馕廴緪夯那闆r較為敏感,當(dāng)出現(xiàn)空氣污染時(shí)隨即就表現(xiàn)出不適應(yīng),從而表現(xiàn)為即時(shí)效應(yīng).CO濃度增加10μg/m3時(shí)在累積滯后0~4d (lag04)時(shí)ER%值達(dá)到最大,此時(shí)對遵義市人體健康危害較大,即CO濃度對呼吸系統(tǒng)的危害主要為累積滯后效應(yīng).以上4種污染物濃度均增加10μg/m3時(shí),總呼吸系統(tǒng)發(fā)病人數(shù)將分別增加0.09%,0.44%,0.37%和0.016%,可見在遵義市酸性污染物SO2和NO2對呼吸系統(tǒng)疾病的危害最大.

    圖8 2014~2016年遵義市污染物日均濃度增加10μg/m3時(shí),呼吸系統(tǒng)疾病發(fā)病人數(shù)在不同滯后條件下的增加百分比(ER%)及95%的置信區(qū)間(95%CI)分布

    從以上4種污染物在最敏感滯后條件下與呼吸系統(tǒng)疾病發(fā)病人數(shù)的暴露-反應(yīng)關(guān)系圖中(圖9)看到,PM2.5的logRR值曲線大致呈單調(diào)線性分布,隨著其濃度的增加,RR值逐漸增大,由于遵義市主要污染物為PM2.5,因此需多加注意污染天氣時(shí)呼吸系統(tǒng)疾病的預(yù)防.其余3種污染物logRR值均隨著濃度的增加呈現(xiàn)“J”型分布,在較低濃度時(shí),污染物對人體健康的危害較小,隨著濃度的持續(xù)增高,logRR值快速上升,對人體健康的風(fēng)險(xiǎn)進(jìn)一步加大,同時(shí)對比logRR值大小發(fā)現(xiàn)SO2等酸性污染物對遵義市呼吸系統(tǒng)疾病發(fā)病的影響最大.

    表6為4種污染物對不同人群和不同類別呼吸系統(tǒng)疾病發(fā)病情況的影響:PM2.5、SO2和CO對女性健康的危害要高于男性,增加10μg/m3時(shí)ER%分別為0.13%,0.63%和0.024%,NO2則相反.對比不同類別的呼吸系統(tǒng)疾病發(fā)現(xiàn),上呼吸道感染發(fā)病人數(shù)主要對SO2酸性氣體的濃度變化有顯著反應(yīng),因?yàn)槿梭w呼吸時(shí)最先受到刺激的為鼻、喉等部分,而下呼吸道感染則在一定程度上均受到以上4種污染物的顯著危害.

    表6 污染物濃度每增加10μg/m3時(shí)呼吸系統(tǒng)疾病發(fā)病人數(shù)增加百分比(ER%)及95%的置信區(qū)間(95%CI)的統(tǒng)計(jì)表

    2.4 氣象和空氣污染因子協(xié)同效應(yīng)對呼吸系統(tǒng)疾病發(fā)病的影響

    由于不同氣象條件對空氣污染的形成和分布有不同的影響,氣象條件的改變可以通過影響大氣中光化學(xué)反應(yīng)從而影響大氣中二次污染物的形成,間接影響人體健康.

    表7 2014~2016年遵義市氣象因子與污染物的Spearman相關(guān)性

    表7為遵義市氣象因子和空氣污染變量之間的Spearman相關(guān)系數(shù)統(tǒng)計(jì)表,氣象因子和空氣污染變量之間呈顯著相關(guān)性,表明二者之間存在固有的理化特征,且對人體健康的影響并不是單因素的,而是存在多因素的協(xié)同作用.

    經(jīng)模型計(jì)算分析,發(fā)現(xiàn)對遵義市呼吸系統(tǒng)疾病發(fā)病有顯著影響的主要是氣溫和NO2、相對濕度和SO2之間的協(xié)同效應(yīng)(表8).

    從β系數(shù)分析,氣溫和NO2協(xié)同效應(yīng)對女性人群的作用比男性強(qiáng);RH與SO2的協(xié)同作用則是對男性患者的影響要大些,且對下呼吸道感染也有顯著影響,但是對上呼吸道感染影響不顯著.

    表8 氣象因子和空氣污染因子的協(xié)同效應(yīng)對呼吸系統(tǒng)疾病發(fā)病的影響

    圖10中(a)-(c)和(d)-(f)分別為氣溫與NO2、RH和SO2對呼吸系統(tǒng)疾病有顯著協(xié)同作用的三維曲面圖.整體上在低溫、高濃度NO2的條件下,呼吸系統(tǒng)疾病發(fā)病人數(shù)最多,即當(dāng)冬半年出現(xiàn)污染天氣時(shí),二者協(xié)同效應(yīng)對呼吸系統(tǒng)疾病發(fā)病的影響增強(qiáng),對上感和下感發(fā)病的影響與總呼吸基本一致.

    圖10 氣象和空氣污染因子的協(xié)同作用對呼吸系統(tǒng)疾病發(fā)病人數(shù)影響的平滑曲面圖

    低濕、高濃度SO2的協(xié)同效應(yīng)對遵義市呼吸系統(tǒng)疾病發(fā)病的影響最大,發(fā)病人數(shù)隨著相對濕度減小(或SO2濃度增加)逐漸增加.此外,男性患者對相對濕度與SO2之間的協(xié)同作用更敏感,且對低濃度SO2也較敏感(圖10(e));同時(shí),低濕與高濃度SO2或低濃度SO2對下呼吸道感染發(fā)病均有顯著影響.

    2.5 呼吸系統(tǒng)疾病發(fā)病的預(yù)測模型研究

    最后,文章采用多元逐步回歸方法對遵義市2014~2016年呼吸系統(tǒng)疾病進(jìn)行預(yù)測模型構(gòu)建研究.

    在進(jìn)行逐步回歸預(yù)報(bào)模型的構(gòu)建中,根據(jù)環(huán)境因子與呼吸系統(tǒng)疾病發(fā)病的相關(guān)性,遴選出“氣溫、氣壓、相對濕度、風(fēng)速” 4個(gè)基本氣象因子和PM2.5、SO2、NO2和CO等污染因子變量.并充分考慮滯后效應(yīng),對挑選出的環(huán)境因子變量相對于患者發(fā)病時(shí)間進(jìn)行提前一周考慮,對所選的氣象因子分別考慮單滯后和累積滯后一周考慮,在累積滯后的基礎(chǔ)上再次進(jìn)行單滯后一周,即考慮累積滯后時(shí)所用到的因子可追溯到兩周前,這與大氣行星尺度的天氣變化周期相對應(yīng),即將大尺度天氣過程的變化周期也予以了考慮,更合理地研究其與疾病發(fā)病與流行的關(guān)系.污染因子則僅考慮其單滯后效應(yīng).最后共派生出416個(gè)可能的影響因子進(jìn)行逐步回歸,構(gòu)建預(yù)報(bào)方程,優(yōu)選出參與構(gòu)建預(yù)報(bào)方程的因子變量及物理意義統(tǒng)計(jì)如表9:

    表9 優(yōu)選出的因子變量及物理意義

    注: (d1)d2為氣象因子連續(xù)d1天內(nèi)進(jìn)行平均,同時(shí)相對發(fā)病人數(shù)提前考慮d2天;()d2為污染因子單滯后d2天,相對發(fā)病人數(shù)提前考慮d2天;(Δ24)d2、(Δ48)d2、(Δ72)d2分別為當(dāng)天和前天的要素值之差,當(dāng)天和前兩天的要素值之差,... ,以此類推,即不考慮日變化的要素變幅,同時(shí)相對發(fā)病人數(shù)提前考慮d2天;(Δ24)d2、(Δ48)d2、(Δ72)d2分別為當(dāng)天的變幅,連續(xù)兩天的變幅,連續(xù)三天的變幅,即是考慮了日變化之后的各因子最大變幅,相對發(fā)病人數(shù)提前考慮d2天.表中,各氣象因子單位分別為:(℃) ,(hPa) , RH (1) , W (m/s).

    選取遵義市2014~2015年共兩年的數(shù)據(jù)資料構(gòu)建逐步回歸預(yù)報(bào)模型,以2016年的數(shù)據(jù)作為試預(yù)報(bào)驗(yàn)證.此外,在建立不同季節(jié)的回歸方程時(shí),考慮到滯后效應(yīng)以及逐季節(jié)預(yù)報(bào)方程轉(zhuǎn)換時(shí)會(huì)產(chǎn)生不連續(xù)的問題,故在構(gòu)建逐季節(jié)預(yù)報(bào)方程的過程中將每個(gè)季節(jié)的疾病和氣象資料均向前和向后延伸半個(gè)月從而保證月際之間要素的連續(xù)性以及緩解月份之間的跳躍性,同時(shí)還增加了建模樣本量,從而使得預(yù)報(bào)方程的連續(xù)性和穩(wěn)定性更好.

    構(gòu)建的呼吸系統(tǒng)疾病逐步回歸預(yù)報(bào)方程見表10,所有預(yù)報(bào)方程均通過了=0.01的顯著性檢驗(yàn).從預(yù)報(bào)方程中的組成變量來看,各因子的作用存在一定的滯后效應(yīng),主導(dǎo)因素因季節(jié)等不同而存在差異.呼吸系統(tǒng)疾病發(fā)病人數(shù)主要與CO濃度、PM2.5、變溫與變壓成正比,與濕度和平均氣溫呈反比,意味著遵義市呼吸系統(tǒng)疾病發(fā)病更易受干冷空氣活動(dòng)的影響.

    表10 遵義市呼吸系統(tǒng)疾病發(fā)病人數(shù)全年和四季預(yù)報(bào)方程

    對所建立的預(yù)報(bào)方程分別進(jìn)行回代和試預(yù)報(bào)效果優(yōu)度檢驗(yàn)(表11和圖11),逐步回歸預(yù)報(bào)方程的回代預(yù)測準(zhǔn)確率均在82%以上,春、秋季節(jié)的準(zhǔn)確率達(dá)87%以上,同時(shí)各項(xiàng)誤差指標(biāo)結(jié)果與之保持一致;在試預(yù)報(bào)檢驗(yàn)中,夏季預(yù)測效果最差(P約為44.4%),全年和其他季節(jié)的準(zhǔn)確率均在75%以上.

    表11 遵義市呼吸系統(tǒng)疾病預(yù)報(bào)方程的回代和試預(yù)報(bào)優(yōu)度檢驗(yàn)

    此外,全年和四季預(yù)報(bào)方程的回代曲線與實(shí)測曲線擬合效果良好(圖11),對于2016年試預(yù)報(bào)檢驗(yàn)而言,年預(yù)報(bào)方程結(jié)果主要在9月份谷值處與實(shí)測值存在較大差異,夏季預(yù)報(bào)方程的曲線擬合效果最差,春、秋和冬季的效果則明顯優(yōu)于全年的.整體上除夏季外,季節(jié)效果要優(yōu)于全年預(yù)測效果,因此建議對遵義市呼吸系統(tǒng)疾病發(fā)病人數(shù)進(jìn)行季節(jié)預(yù)報(bào).

    3 討論

    隨著“健康中國2030”戰(zhàn)略的提出,展開不同區(qū)域氣象和空氣污染健康效應(yīng)評估、對戰(zhàn)略的實(shí)施和氣象敏感性疾病預(yù)防具有重要意義.同時(shí),我國地域廣闊,地貌地形復(fù)雜多樣,氣象條件和空氣污染區(qū)域性特征差異較大,以往醫(yī)學(xué)氣象學(xué)研究多集中于京津冀、南京等環(huán)境條件變化顯著區(qū)域,對遵義等地區(qū)的健康效應(yīng)評估較為缺失.本文選取遵義市2012~ 2016年的疾病資料及同時(shí)期的環(huán)境因子資料,分析了遵義市呼吸系統(tǒng)疾病發(fā)病高峰期主要在冬春季,并從24節(jié)氣的角度揭示了立秋到處暑期間的氣候環(huán)境條件對呼吸系統(tǒng)疾病有一定的氣候療養(yǎng)效應(yīng).利用DLNM、GLM和GAM模型全面的探討了環(huán)境要素對呼吸系統(tǒng)疾病發(fā)病的影響及彼此之間的協(xié)同作用,并利用逐步回歸方法進(jìn)行了發(fā)病人數(shù)預(yù)測模型的構(gòu)建與檢驗(yàn),

    從氣象要素對遵義市呼吸系統(tǒng)疾病的影響和滯后性來看,本文研究結(jié)果與其他不同區(qū)域的研究結(jié)果具有較好的一致性[34-37],氣溫對呼吸系統(tǒng)疾病的影響主要表現(xiàn)為高溫即時(shí)效應(yīng)和低溫滯后效應(yīng),由于遵義市常年高溫日數(shù)不多,因此對當(dāng)?shù)囟砸缘蜏販笮?yīng)為主.遵義市相對濕度相對較高,當(dāng)?shù)厝巳簩Ω邼癍h(huán)境具有較好的適應(yīng)性,因此當(dāng)相對濕度偏低時(shí),反而容易引起呼吸系統(tǒng)疾病發(fā)病,且以滯后效應(yīng)為主.綜合考慮各氣象要素之間的相互作用,從舒適度的角度揭示了冷、熱斜脅迫下呼吸系統(tǒng)疾病發(fā)病人數(shù)明顯多于舒適情況,且熱強(qiáng)迫下女性發(fā)病人數(shù)略多于男性,冷強(qiáng)迫下則相反.

    在空氣污染因子與呼吸系統(tǒng)疾病的暴露-反應(yīng)關(guān)系研究中,除CO外,其余污染變量主要以即時(shí)效應(yīng)為主,而研究表明京津冀以及蘭州等污染較重區(qū)域的污染因子對呼吸系統(tǒng)疾病發(fā)病有一定的滯后性[38-41].這是因?yàn)檎w上遵義市空氣質(zhì)量優(yōu)良率很高,人們對污染物濃度的變化更加敏感,當(dāng)污染物濃度升高時(shí),能在當(dāng)天就表現(xiàn)出不適,使呼吸系統(tǒng)疾病就診人數(shù)上升,污染物濃度越高危害越大.同時(shí),酸性污染物對呼吸道的刺激較大,其對呼吸系統(tǒng)疾病的影響最顯著.

    氣象因子和空氣污染因子之間存在緊密聯(lián)系,一方面附著有有毒有害物質(zhì)的大氣污染物可直接進(jìn)入血液,作用于呼吸系統(tǒng),另一方面可作為細(xì)菌病毒等傳播媒介,而氣象要素則可通過影響原微生物的生存、繁殖和傳播等從而間接誘發(fā)疾病,已有證據(jù)指出氣象要素和大氣污染在不良健康效應(yīng)上可能存在協(xié)同作用[42-43].利用GAM模型分析優(yōu)選,發(fā)現(xiàn)在低溫、高濃度NO2的條件下,呼吸系統(tǒng)疾病發(fā)病人數(shù)較多,女性人群為易感脆弱人群;在低濕、高濃度SO2時(shí)呼吸系統(tǒng)疾病發(fā)病人數(shù)最多,且發(fā)病人數(shù)隨著相對濕度減小(或SO2濃度增加)逐漸增加,同時(shí)低濃度SO2對下呼吸道感染發(fā)病也有顯著影響.

    最后,建立了遵義市全年和季節(jié)呼吸系統(tǒng)疾病發(fā)病人數(shù)的預(yù)測方程,預(yù)測值與實(shí)測值曲線趨勢高度一致,季節(jié)效果要優(yōu)于全年預(yù)測效果,故建議對遵義市呼吸系統(tǒng)疾病發(fā)病人數(shù)進(jìn)行不同季節(jié)預(yù)報(bào).

    4 結(jié)論

    4.1 冬春季為遵義市呼吸系統(tǒng)疾病發(fā)病高峰期,立秋至處暑前后發(fā)病人數(shù)最少.

    4.2 高溫對呼吸系統(tǒng)疾病發(fā)病以即時(shí)效應(yīng)為主,低溫主要表現(xiàn)為滯后效應(yīng),氣壓和相對濕度則分別以高壓和低濕的滯后效應(yīng)為主,平均風(fēng)速的影響不大;在最敏感滯后條件下,PM2.5與呼吸系統(tǒng)疾病發(fā)病的暴露關(guān)系反應(yīng)曲線呈單調(diào)線性分布,SO2、NO2和CO則呈“J”型分布;低溫與高濃度NO2或者低濕與高濃度SO2的協(xié)同作用對呼吸系統(tǒng)疾病發(fā)病均有顯著影響.

    4.3 呼吸系統(tǒng)疾病發(fā)病人數(shù)的季節(jié)預(yù)報(bào)方程的預(yù)報(bào)效果顯著優(yōu)于年預(yù)報(bào)方程的效果.

    [1] Liu L, Breitner S, Pan X, et al. Associations between air temperature and cardio-respiratory mortality in the urban area of Beijing, China: a time-series analysis. [J]. Environmental Health, 2011,10(1):1-11.

    [2] Tang J W, Lai F Y L, Nymadawa P, et al. Comparison of the incidence of influenza in relation to climate factors during 2000~2007 in five countries [J]. Journal of Medical Virology, 2010,82(11):1958-1965.

    [3] 趙立宇,余勇夫,席俊文,等.上海市氣溫變化對醫(yī)院兒科日門診人次的影響 [J]. 中國醫(yī)院統(tǒng)計(jì), 2010,17(4):313-316.

    [4] Nascimento-Carvalho C M, Cardoso M R, Barral A, et al. Seasonal patterns of viral and bacterial infections among children hospitalized with community-acquired pneumonia in a tropical region [J]. Scandinavian Journal of Infectious Diseases, 2010,42(11/12):839-844.

    [5] 錢孝琳,闞海東.大氣顆粒物污染對心血管系統(tǒng)影響的流行病學(xué)研究進(jìn)展 [J]. 中華流行病學(xué)雜志, 2005,26(12):999-1001.

    [6] 闞海東,陳秉衡,賈 健.大氣污染物暴露與人群健康效應(yīng)的暴露-反應(yīng)關(guān)系分析 [J]. 環(huán)境與健康雜志, 2004,21(7):253-254.

    [7] 張 瑩,王式功,賈旭偉,等.氣溫與PM2.5協(xié)同作用對疾病急診就診人數(shù)的影響 [J]. 中國環(huán)境科學(xué), 2017,37(8):3175-3182.

    [8] 王敏珍,鄭 山,王式功,等.氣溫與濕度的交互作用對呼吸系統(tǒng)疾病的影響 [J]. 中國環(huán)境科學(xué), 2016,36(2):581-588.

    [9] Bleakley H. Disease and development: evidence from hookworm eradication in the American South [J]. The Quarterly Journal of Economics, 2007,122(1):73-117.

    [10] 杜鵬瑞,杜 睿,任偉珊.城市大氣顆粒物毒性效應(yīng)及機(jī)制的研究進(jìn)展 [J]. 中國環(huán)境科學(xué), 2016,36(9):2815-2827.

    [11] 郝吉明,程 真,王書肖.我國大氣環(huán)境污染現(xiàn)狀及防治措施研究 [J]. 環(huán)境保護(hù), 2012,(9):16-20.

    [12] 代玉田,于鳳英,常 平,等.德州市高血壓發(fā)病氣象條件與預(yù)報(bào)模式 [J]. 氣象與環(huán)境學(xué)報(bào), 2012,28(5):79-82.

    [13] Climate change 2007-the physical science basis: Working group II contribution to the Fourth assessment report of the IPCC [M]. Cambridge University Press, 2007.

    [14] Climate change 2013-the physical science basis: Working group I contribution to the Fifth assessment report of the IPCC [M]. Cambridge University Press, 2013.

    [15] http://www.cma.gov.cn/2011xzt/2015zt/qmttt/20150813/2015041004/201508/t20150813_290458.html

    [16] Michelozzi P, Accetta G, De Sario M, et al. High temperature and hospitalizations for cardiovascular and respiratory causes in 12European cities [J]. American Journal of Respiratory and Critical Care Medicine, 2009,179(5):383-389.

    [17] 劉建軍,鄭有飛,吳榮軍.熱浪災(zāi)害對人體健康的影響及其方法研究 [J]. 自然災(zāi)害學(xué)報(bào), 2008,17(1):151-156.

    [18] Semenza J C, McCullough J E, Flanders W D, et al. Excess hospital admissions during the July 1995 heat wave in Chicago [J]. American journal of preventive medicine, 1999,16(4):269-277.

    [19] Morelli X, Rieux C, Cyrys J, et al. Air pollution, health and social deprivation: a fine-scale risk assessment [J]. Environmental Research, 2016,147:59-70.

    [20] Lin Y K, Chang S C, Lin C S, et al. Comparing ozone metrics on associations with outpatient visits for respiratory diseases in Taipei Metropolitan area [J]. Environmental Pollution, 2013,177:177-184.

    [21] Shen Y, Wu Y, Chen G, et al. Non-linear increase of respiratory diseases and their costs under severe air pollution [J]. Environmental Pollution, 2017,224:631-637.

    [22] 謝元博,陳 娟,李 巍.霧霾重污染期間北京居民對高濃度PM2.5持續(xù)暴露的健康風(fēng)險(xiǎn)及其損害價(jià)值評估 [J]. 環(huán)境科學(xué), 2014,35(1): 1-8.

    [23] 張 瑩.我國典型城市空氣污染特征及其健康影響和預(yù)報(bào)研究 [D]. 蘭州:蘭州大學(xué), 2016.

    [24] Smith K R, Corvalán C F, Kjellstrom T. How much global ill health is attributable to environmental factors? [J]. Epidemiology-Baltimore, 1999,10(5):573-584.

    [25] WHO. International classification of diseases (ICD)-10online [EB/ OL].http://www.whoint/classifications/apps/icd/icd10online//,2007.1. 10.

    [26] 楊 軍,歐春泉,丁 研,等.分布滯后非線性模型 [J]. 中國衛(wèi)生統(tǒng)計(jì), 2012,29(5):772-773.

    [27] Gasparrini A. Distributed lag linear and non-linear models in R: the package dlnm [J]. Journal of Statistical Software, 2011,43(8):1.

    [28] Gasparrini A, Armstrong B, Kenward M G. Distributed lag non-linear models [J]. Statistics in Medicine, 2010,29(21):2224-2234.

    [29] 汪建均.基于廣義線性模型的變量選擇與穩(wěn)健參數(shù)設(shè)計(jì) [D]. 南京:南京理工大學(xué), 2011.

    [30] 王式功,馬 盼,尚可政等.一種基于黃金分割法構(gòu)建人體舒適度計(jì)算方法 [P]. 中國: 2012105963660,2013-05-15.

    [31] 張志薇.基于舒適度研究的中國大陸健康氣候評估 [D]. 蘭州:蘭州大學(xué), 2014.

    [32] Kan H, Chen B, Hong C. Health impact of outdoor air pollution in China: current knowledge and future research needs [J]. Environmental Health Perspectives, 2009,117(5):A187.

    [33] Dominici F, Daniels M, Zeger S L, et al. Air pollution and mortality: estimating regional and national dose-response relationships [J]. Journal of the American Statistical Association, 2002,97(457): 100-111.

    [34] 李 娟,張志薇,于庚康,等.氣象要素對南京市呼吸系統(tǒng)疾病的影響研究 [J]. 氣象科學(xué), 2017,37(3):409-415.

    [35] González S, Díaz J, Pajares M S, et al. Relationship between atmospheric pressure and mortality in the Madrid Autonomous Region: a time-series study [J]. International Journal of Biometeorology, 2001,45(1):34-40.

    [36] Díaz J, García R, Lopez C, et al. Mortality impact of extreme winter temperatures [J]. International Journal of Biometeorology, 2005,49(3):179-183.

    [37] 馬 盼,王式功,尚可政,等.氣象舒適條件對呼吸系統(tǒng)疾病的影響 [J]. 中國環(huán)境科學(xué), 2018,38(1):374-382.

    [38] 馬 盼,李若麟,趙笑顏,等.北京市空氣質(zhì)量指數(shù)與氣象條件的關(guān)系及其對呼吸系統(tǒng)疾病的滯后影響 [J]. 蘭州大學(xué)學(xué)報(bào):自然科學(xué)版, 2017,53(3):388-393.

    [39] 張 瑩,邵 毅,王式功,等.北京市空氣污染物對呼吸系統(tǒng)疾病門診人數(shù)的影響 [J]. 中國環(huán)境科學(xué), 2014,34(9):2401-2407.

    [40] 吳 彧,張衍燊,劉 玲,等.鄭州市空氣污染對每日門診人次影響的滯后效應(yīng)和收獲效應(yīng) [J]. 環(huán)境與健康雜志, 2011,28(3):233-237.

    [41] 陶燕,羊德容,蘭嵐,等. 蘭州市空氣污染對呼吸系統(tǒng)疾病入院人數(shù)的影響[J]. 中國環(huán)境科學(xué), 2013, 33(1): 175-180.

    [42] Green R S, Basu R, Malig B, et al. The effect of temperature on hospital admissions in nine California counties [J]. International Journal of Public Health, 2010,55(2):113-121.

    [43] 伍燕珍,張金良.北京市大氣顆粒物與溫度對人群非意外死亡影響的交互作用 [J]. 環(huán)境科學(xué)研究, 2009,(12):1403-1410.

    Study about the impact of environmental conditions on respiratory diseases and prediction in Zunyi City.

    YUE Man1, WANG Shi-gong2,1*, XIE Jia-jun3, MA Pan2, SHANG Ke-zheng1

    (1.Center for Meteorological Environment and Human Health, College of Atmospheric Sciences, Lanzhou University, Lanzhou 730000, China;2.College of Atmospheric Sciences, Chengdu University of Information Technology, Plateau Atmosphere and Environment Key Laboratory of Sichuan Province, Chengdu 610225, China;3.Zunyi Meteorological Bureau, Guizhou, Zunyi 563000, China)., 2018,38(11):4334~4347

    To explore environmental impacts such as the meteorological and air pollution factors on respiratory diseases and to provide scientific basis for the prevention of these diseases in Zunyi City, distributed lag non-linear model together with generalized linear and additive models are applied to analyse the exposure-response relationship between environmental factors and respiratory diseases from 2012 to 2016 in Zunyi City. Results show that the changes in respiratory diseases are mainly consistent with the local long-term climatic conditions, and the impact of climate change is dominant. Winter and spring are the peak periods with high respiratory diseases number, and during the Start of Autumn and Stopping the heat periods, the respiratory diseases number is the lowest, indicating that local climatic conditions have positive climatic effects on patients with respiratory diseases during this time period. The impact of temperature on respiratory disease is mainly low-temperature lagged effect with the patients increasing by 31.6% (95%CI: 4.4%~65.8%) if the temperature changes by 1°C.The pressure mainly has the high-pressure lagged effect on respiratory diseases, while the relative humidity has both immediate and lagged effects in lower humidity. The number of respiratory diseases is significantly higher under cold and hot uncomfortable levels than that of comfort levels. PM2.5, SO2, and NO2mainly show immediate effects on respiratory diseases, while the CO had the highest risk if lagged for four days. The exposure-response relationship between respiratory diseases and PM2.5shows an monotonously linear distribution, while that of SO2, NO2, and CO are J-type distribution, and the synergistic effects between low temperature and high concentration NO2or low humidity and high concentration SO2both have significant impact on respiratory diseases. The accuracy of annual and seasonal regression equations is over 75% (except for summer equation), and the seasonal equations’ prediction effect is better than that of the annual equation.

    environmental conditions;meteorological factors;air pollution;respiratory diseases

    X503.1,P49

    A

    1000-6923(2018)11-4334-14

    樂 滿(1992-),女,湖北咸寧人,碩士,主要從事氣象環(huán)境與健康研究.發(fā)表論文2篇.

    2018-04-19

    國家重點(diǎn)研發(fā)計(jì)劃“全球變化及應(yīng)對”重點(diǎn)專項(xiàng)(2016YFA0602004);國家自然科學(xué)基金重大研究計(jì)劃重點(diǎn)支持項(xiàng)目(91644226);國家基礎(chǔ)科技條件平臺建設(shè)項(xiàng)目(NCMI-SBS17-201707、NCMI-SJS15-201707)

    * 責(zé)任作者, 教授, wangsg@lzu.edu.cn

    猜你喜歡
    滯后效應(yīng)遵義市氣象
    我為旗幟增輝
    ——貴州省遵義市綏陽縣老年大學(xué)校歌
    氣象樹
    《內(nèi)蒙古氣象》征稿簡則
    氣象戰(zhàn)士之歌
    心聲歌刊(2019年5期)2020-01-19 01:52:52
    遵義市關(guān)工委:圍繞示范抓鞏固“五好”組織顯活力
    晚晴(2019年4期)2019-07-08 03:49:42
    大國氣象
    基于環(huán)境保護(hù)的企業(yè)社會(huì)責(zé)任與企業(yè)財(cái)務(wù)績效的關(guān)系研究
    商情(2016年32期)2017-03-04 00:54:25
    城鎮(zhèn)化中人口結(jié)構(gòu)變化與經(jīng)濟(jì)增長的關(guān)系
    喜迎新年
    ——遵義市春節(jié)聯(lián)歡晚會(huì)
    遵義(2016年4期)2016-08-23 11:47:15
    習(xí)近平總書記在遵義市楓香鎮(zhèn)花茂村考察
    日韩一区二区视频免费看| eeuss影院久久| 欧美不卡视频在线免费观看| 日韩一区二区视频免费看| 国产aⅴ精品一区二区三区波| 男人狂女人下面高潮的视频| 欧美xxxx性猛交bbbb| 一区二区三区免费毛片| 俄罗斯特黄特色一大片| 九九热线精品视视频播放| 久久精品国产亚洲av香蕉五月| 国产 一区精品| 寂寞人妻少妇视频99o| 久久久久国产网址| 亚洲av电影不卡..在线观看| 一个人观看的视频www高清免费观看| 国产一区二区亚洲精品在线观看| 色吧在线观看| 最新中文字幕久久久久| 人人妻,人人澡人人爽秒播| 人妻夜夜爽99麻豆av| 国产伦一二天堂av在线观看| 国产在线精品亚洲第一网站| 国产91av在线免费观看| 韩国av在线不卡| 能在线免费观看的黄片| 日本欧美国产在线视频| 美女内射精品一级片tv| 精品免费久久久久久久清纯| 淫秽高清视频在线观看| 色综合站精品国产| 久久人人爽人人爽人人片va| 午夜老司机福利剧场| 三级毛片av免费| 超碰av人人做人人爽久久| 亚洲av免费高清在线观看| 国产亚洲欧美98| 国产一区亚洲一区在线观看| 激情 狠狠 欧美| 97热精品久久久久久| 人妻久久中文字幕网| 成人一区二区视频在线观看| 天天一区二区日本电影三级| 日韩精品青青久久久久久| 国产国拍精品亚洲av在线观看| 一区二区三区四区激情视频 | 22中文网久久字幕| 麻豆国产97在线/欧美| av视频在线观看入口| 99在线视频只有这里精品首页| 又爽又黄无遮挡网站| 久久久久精品国产欧美久久久| 小蜜桃在线观看免费完整版高清| 九九在线视频观看精品| av中文乱码字幕在线| 91午夜精品亚洲一区二区三区| АⅤ资源中文在线天堂| 性色avwww在线观看| 久久草成人影院| 日韩大尺度精品在线看网址| 国产女主播在线喷水免费视频网站 | 黄色日韩在线| 国产精品福利在线免费观看| 哪里可以看免费的av片| 亚洲国产精品成人久久小说 | av黄色大香蕉| 日本黄大片高清| 人人妻,人人澡人人爽秒播| 舔av片在线| 麻豆精品久久久久久蜜桃| 老女人水多毛片| 国产精品亚洲一级av第二区| 亚洲真实伦在线观看| 舔av片在线| 97超级碰碰碰精品色视频在线观看| 日本一本二区三区精品| 三级经典国产精品| 日日啪夜夜撸| 亚洲精品粉嫩美女一区| 久久久久九九精品影院| 又爽又黄a免费视频| 亚洲国产高清在线一区二区三| 美女大奶头视频| 国产av不卡久久| 亚洲国产欧美人成| av在线播放精品| 可以在线观看的亚洲视频| aaaaa片日本免费| 久久久欧美国产精品| 人妻夜夜爽99麻豆av| 精品一区二区免费观看| 少妇猛男粗大的猛烈进出视频 | 99热这里只有是精品在线观看| 中文字幕人妻熟人妻熟丝袜美| 中国美白少妇内射xxxbb| 亚洲欧美精品自产自拍| 亚洲成人中文字幕在线播放| 男女边吃奶边做爰视频| 久久午夜福利片| 免费看日本二区| 美女被艹到高潮喷水动态| eeuss影院久久| 亚洲中文字幕日韩| 久久精品国产鲁丝片午夜精品| 亚洲欧美日韩卡通动漫| 亚洲高清免费不卡视频| 精品欧美国产一区二区三| 精品无人区乱码1区二区| 一级av片app| 亚洲人成网站在线播| 99久国产av精品| 久久精品国产亚洲av天美| 久久久精品大字幕| 国产精品亚洲美女久久久| 九九在线视频观看精品| 亚洲人成网站在线播放欧美日韩| 中文字幕av在线有码专区| 久久久欧美国产精品| h日本视频在线播放| 国产伦精品一区二区三区四那| 成年女人毛片免费观看观看9| 老司机影院成人| 久久久久性生活片| 精品日产1卡2卡| 真实男女啪啪啪动态图| 啦啦啦韩国在线观看视频| 中文字幕av在线有码专区| 国产伦一二天堂av在线观看| 18+在线观看网站| 国产精品亚洲美女久久久| 亚洲专区国产一区二区| 1000部很黄的大片| ponron亚洲| 婷婷精品国产亚洲av在线| 国产精品综合久久久久久久免费| 麻豆久久精品国产亚洲av| 精品一区二区三区视频在线观看免费| 一级黄色大片毛片| 亚洲人成网站在线播放欧美日韩| 日本-黄色视频高清免费观看| 精品一区二区三区视频在线观看免费| 日韩制服骚丝袜av| 99久久精品热视频| 欧美性感艳星| 日日摸夜夜添夜夜爱| 99riav亚洲国产免费| 国模一区二区三区四区视频| 久久久久国内视频| 亚洲人成网站在线播| 国产免费男女视频| 看非洲黑人一级黄片| 欧美精品国产亚洲| 国产亚洲欧美98| 国产精品人妻久久久久久| 99久久久亚洲精品蜜臀av| 色吧在线观看| 一夜夜www| 日本三级黄在线观看| 久久鲁丝午夜福利片| 特大巨黑吊av在线直播| 亚洲av二区三区四区| 国产真实伦视频高清在线观看| 少妇裸体淫交视频免费看高清| 精品99又大又爽又粗少妇毛片| 在线a可以看的网站| 亚洲av成人av| 老司机福利观看| 国产精品久久久久久久电影| 有码 亚洲区| 最后的刺客免费高清国语| 人妻丰满熟妇av一区二区三区| 国产视频一区二区在线看| 菩萨蛮人人尽说江南好唐韦庄 | 成人漫画全彩无遮挡| av卡一久久| 国产成人a区在线观看| 三级男女做爰猛烈吃奶摸视频| eeuss影院久久| 搡老妇女老女人老熟妇| 亚洲国产精品久久男人天堂| 1000部很黄的大片| 国产乱人视频| 国产精品1区2区在线观看.| 亚洲国产精品久久男人天堂| 老女人水多毛片| 亚洲欧美日韩高清在线视频| 在线播放无遮挡| 天天一区二区日本电影三级| 精品一区二区三区视频在线观看免费| 日韩一本色道免费dvd| 一级a爱片免费观看的视频| 99热只有精品国产| 九九热线精品视视频播放| 亚洲av熟女| 毛片女人毛片| 99在线视频只有这里精品首页| 国产 一区 欧美 日韩| 在线观看66精品国产| 精品国产三级普通话版| 久久久久国内视频| 99热这里只有是精品50| 亚洲av二区三区四区| 搡女人真爽免费视频火全软件 | av视频在线观看入口| 夜夜夜夜夜久久久久| 亚洲经典国产精华液单| 最新中文字幕久久久久| 卡戴珊不雅视频在线播放| 人妻久久中文字幕网| 国产精品国产高清国产av| av中文乱码字幕在线| 91在线精品国自产拍蜜月| 韩国av在线不卡| 欧美不卡视频在线免费观看| 亚洲性久久影院| 久久精品夜色国产| 免费在线观看影片大全网站| 伊人久久精品亚洲午夜| 欧美潮喷喷水| 中文在线观看免费www的网站| 热99re8久久精品国产| 久久综合国产亚洲精品| 精品不卡国产一区二区三区| 婷婷精品国产亚洲av| 成年女人毛片免费观看观看9| 啦啦啦观看免费观看视频高清| 国产精品乱码一区二三区的特点| 又粗又爽又猛毛片免费看| 亚洲美女黄片视频| 亚洲美女视频黄频| 淫秽高清视频在线观看| 国产精品久久久久久av不卡| 观看免费一级毛片| 三级国产精品欧美在线观看| 天堂√8在线中文| 国产高清有码在线观看视频| 中文字幕人妻熟人妻熟丝袜美| av.在线天堂| 国产精品,欧美在线| 国产亚洲av嫩草精品影院| 国产精品一区二区免费欧美| 国产三级中文精品| 久久精品影院6| 女的被弄到高潮叫床怎么办| 在线天堂最新版资源| 国产亚洲91精品色在线| 亚洲性夜色夜夜综合| 乱码一卡2卡4卡精品| 人妻夜夜爽99麻豆av| 久久久久久大精品| 精品午夜福利视频在线观看一区| 日本成人三级电影网站| 国产一区二区在线av高清观看| 色视频www国产| 国产大屁股一区二区在线视频| h日本视频在线播放| 国产私拍福利视频在线观看| 能在线免费观看的黄片| 国产伦在线观看视频一区| 国产爱豆传媒在线观看| 99riav亚洲国产免费| 青春草视频在线免费观看| 大型黄色视频在线免费观看| 综合色丁香网| 一级毛片aaaaaa免费看小| av中文乱码字幕在线| av视频在线观看入口| 天堂影院成人在线观看| 亚洲中文字幕一区二区三区有码在线看| 99热网站在线观看| 亚洲精品久久国产高清桃花| 麻豆av噜噜一区二区三区| 国产伦精品一区二区三区四那| 亚洲中文字幕日韩| 婷婷六月久久综合丁香| 97在线视频观看| 免费搜索国产男女视频| 亚洲欧美中文字幕日韩二区| 高清毛片免费看| 国产精品久久电影中文字幕| 午夜久久久久精精品| 一个人看视频在线观看www免费| 日韩欧美 国产精品| 男女做爰动态图高潮gif福利片| 国产精品精品国产色婷婷| 好男人在线观看高清免费视频| 欧美成人免费av一区二区三区| 啦啦啦观看免费观看视频高清| 免费av不卡在线播放| 波多野结衣高清作品| 春色校园在线视频观看| 国产精品一区二区三区四区免费观看 | 亚洲av中文av极速乱| 校园人妻丝袜中文字幕| 十八禁国产超污无遮挡网站| 国产女主播在线喷水免费视频网站 | 精品午夜福利在线看| 亚洲精品影视一区二区三区av| 日本黄色片子视频| 小蜜桃在线观看免费完整版高清| 两个人的视频大全免费| 色综合亚洲欧美另类图片| 免费大片18禁| 精品无人区乱码1区二区| av黄色大香蕉| 国产精品乱码一区二三区的特点| 在线观看av片永久免费下载| 国产精品久久久久久精品电影| 18禁在线播放成人免费| 51国产日韩欧美| 真实男女啪啪啪动态图| 亚洲色图av天堂| 久久欧美精品欧美久久欧美| eeuss影院久久| 亚洲最大成人av| 熟女人妻精品中文字幕| 亚洲国产日韩欧美精品在线观看| 婷婷六月久久综合丁香| 国产成人aa在线观看| 伊人久久精品亚洲午夜| 国产成人福利小说| 在线免费观看不下载黄p国产| 日韩成人伦理影院| 女生性感内裤真人,穿戴方法视频| 成人鲁丝片一二三区免费| 亚洲专区国产一区二区| 18禁在线播放成人免费| 国产综合懂色| 成人国产麻豆网| 亚洲av免费高清在线观看| 亚洲性久久影院| 美女黄网站色视频| 在线看三级毛片| 观看免费一级毛片| 成人三级黄色视频| 1024手机看黄色片| 亚洲三级黄色毛片| 亚洲精品在线观看二区| 一个人观看的视频www高清免费观看| 精品久久久久久久久久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 少妇的逼水好多| 美女 人体艺术 gogo| 99久国产av精品| 午夜视频国产福利| 精品熟女少妇av免费看| 91麻豆精品激情在线观看国产| 国产精品嫩草影院av在线观看| 免费观看精品视频网站| 亚洲激情五月婷婷啪啪| 欧美一区二区国产精品久久精品| 天天躁日日操中文字幕| 欧美中文日本在线观看视频| 亚洲久久久久久中文字幕| 中文资源天堂在线| 波多野结衣高清无吗| 岛国在线免费视频观看| 亚洲在线观看片| 午夜福利18| 国产精华一区二区三区| av在线亚洲专区| 日韩欧美三级三区| 麻豆精品久久久久久蜜桃| 欧美最新免费一区二区三区| 国产精品三级大全| 一级av片app| 三级国产精品欧美在线观看| а√天堂www在线а√下载| 国内精品久久久久精免费| 久久精品久久久久久噜噜老黄 | 色播亚洲综合网| 成年版毛片免费区| 亚洲婷婷狠狠爱综合网| 亚洲电影在线观看av| 高清毛片免费观看视频网站| 久久精品91蜜桃| 高清毛片免费观看视频网站| 人妻夜夜爽99麻豆av| 亚洲国产精品成人久久小说 | 久久久久久久午夜电影| 国产精品一区二区免费欧美| 久久久成人免费电影| 国产淫片久久久久久久久| 色哟哟哟哟哟哟| 亚洲婷婷狠狠爱综合网| 国产精品乱码一区二三区的特点| 国产男人的电影天堂91| 免费黄网站久久成人精品| 97碰自拍视频| 给我免费播放毛片高清在线观看| 国产伦一二天堂av在线观看| 日本爱情动作片www.在线观看 | 99久久中文字幕三级久久日本| 麻豆国产av国片精品| 精品久久久久久久人妻蜜臀av| 91久久精品国产一区二区三区| 久久人妻av系列| 婷婷色综合大香蕉| 在线免费观看不下载黄p国产| 欧美+亚洲+日韩+国产| 亚洲va在线va天堂va国产| 国产午夜精品论理片| 国产精品,欧美在线| 国产成人a区在线观看| 少妇丰满av| 在线看三级毛片| 亚洲成人精品中文字幕电影| 能在线免费观看的黄片| 亚洲欧美精品自产自拍| 毛片一级片免费看久久久久| 免费看日本二区| 自拍偷自拍亚洲精品老妇| 亚洲国产色片| 国产高清不卡午夜福利| 在线观看美女被高潮喷水网站| 成人高潮视频无遮挡免费网站| 丰满乱子伦码专区| 国内精品久久久久精免费| 九九久久精品国产亚洲av麻豆| 少妇人妻精品综合一区二区 | 欧美绝顶高潮抽搐喷水| 久久久精品94久久精品| 久久久久久久久中文| 国产国拍精品亚洲av在线观看| 亚洲国产精品久久男人天堂| 亚洲中文字幕日韩| 亚洲七黄色美女视频| 久久亚洲国产成人精品v| 五月玫瑰六月丁香| 人妻久久中文字幕网| 91午夜精品亚洲一区二区三区| 亚洲国产精品成人综合色| av在线亚洲专区| 精品午夜福利在线看| 久久精品国产99精品国产亚洲性色| 国产真实伦视频高清在线观看| 免费观看的影片在线观看| 青春草视频在线免费观看| 久久精品久久久久久噜噜老黄 | 欧美极品一区二区三区四区| 美女免费视频网站| 成人av在线播放网站| 成人高潮视频无遮挡免费网站| 女的被弄到高潮叫床怎么办| 欧美区成人在线视频| 乱人视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲成av人片在线播放无| 美女大奶头视频| 麻豆乱淫一区二区| 国产免费一级a男人的天堂| 丝袜美腿在线中文| 成人精品一区二区免费| 婷婷精品国产亚洲av在线| 在线看三级毛片| 中文亚洲av片在线观看爽| 观看美女的网站| 国内精品久久久久精免费| 三级经典国产精品| 欧美性感艳星| 日韩一区二区视频免费看| 热99在线观看视频| 国产中年淑女户外野战色| 久久久久免费精品人妻一区二区| 日韩欧美免费精品| 免费一级毛片在线播放高清视频| 亚洲一级一片aⅴ在线观看| 中出人妻视频一区二区| 女同久久另类99精品国产91| 亚洲四区av| 亚洲熟妇熟女久久| 精品一区二区三区视频在线| 淫妇啪啪啪对白视频| 久久99热这里只有精品18| 国产精品久久久久久久久免| 在线观看66精品国产| 中国国产av一级| 亚洲中文日韩欧美视频| a级毛片免费高清观看在线播放| 亚洲国产精品久久男人天堂| 在线免费观看的www视频| 国产综合懂色| av免费在线看不卡| 成年av动漫网址| 人人妻人人澡欧美一区二区| 亚洲成人久久爱视频| 国内精品美女久久久久久| 午夜福利视频1000在线观看| 欧美激情久久久久久爽电影| 天美传媒精品一区二区| av在线亚洲专区| 免费高清视频大片| 成人美女网站在线观看视频| 国产成年人精品一区二区| 桃色一区二区三区在线观看| 亚洲精品一卡2卡三卡4卡5卡| av免费在线看不卡| 国产黄色视频一区二区在线观看 | 六月丁香七月| 99精品在免费线老司机午夜| 亚洲成人精品中文字幕电影| 亚洲av熟女| 成熟少妇高潮喷水视频| 国产精品一区二区三区四区免费观看 | 精华霜和精华液先用哪个| 插阴视频在线观看视频| 日韩一区二区视频免费看| 久久热精品热| 成人午夜高清在线视频| 国产在线精品亚洲第一网站| 伊人久久精品亚洲午夜| 日本与韩国留学比较| 色播亚洲综合网| 高清毛片免费看| 国产三级中文精品| 狂野欧美白嫩少妇大欣赏| 91午夜精品亚洲一区二区三区| 国产成人aa在线观看| 看十八女毛片水多多多| 久久午夜亚洲精品久久| 成人毛片a级毛片在线播放| 成人性生交大片免费视频hd| 亚洲一区二区三区色噜噜| 嫩草影院入口| 亚洲色图av天堂| 亚洲aⅴ乱码一区二区在线播放| 国产成人freesex在线 | 免费不卡的大黄色大毛片视频在线观看 | 亚洲欧美日韩无卡精品| 亚洲天堂国产精品一区在线| 在线观看一区二区三区| 国产伦精品一区二区三区四那| 欧美在线一区亚洲| 身体一侧抽搐| 亚洲av电影不卡..在线观看| 亚洲熟妇熟女久久| 亚洲欧美日韩高清在线视频| 综合色丁香网| 网址你懂的国产日韩在线| 免费人成视频x8x8入口观看| 亚洲成人久久性| 可以在线观看的亚洲视频| 日韩精品青青久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 国国产精品蜜臀av免费| 伦理电影大哥的女人| 99在线人妻在线中文字幕| 长腿黑丝高跟| 99久久精品一区二区三区| 免费人成视频x8x8入口观看| 日韩av不卡免费在线播放| 亚洲人成网站在线播| 久久精品综合一区二区三区| 成人永久免费在线观看视频| 丝袜美腿在线中文| 中文亚洲av片在线观看爽| 亚洲欧美中文字幕日韩二区| 免费看a级黄色片| 国产视频内射| 2021天堂中文幕一二区在线观| 国产亚洲欧美98| 精品久久久久久久久亚洲| 国产91av在线免费观看| 日韩欧美一区二区三区在线观看| 色尼玛亚洲综合影院| 在线免费观看的www视频| 可以在线观看毛片的网站| 黄色视频,在线免费观看| 欧美+日韩+精品| 日日摸夜夜添夜夜爱| 亚洲人成网站在线播放欧美日韩| 99久久精品国产国产毛片| 99久久久亚洲精品蜜臀av| 熟女人妻精品中文字幕| 欧美日本亚洲视频在线播放| 午夜福利成人在线免费观看| 国产成人91sexporn| 大又大粗又爽又黄少妇毛片口| 久久久久久大精品| 两性午夜刺激爽爽歪歪视频在线观看| 欧美激情在线99| 午夜福利高清视频| 欧美另类亚洲清纯唯美| 国产一区二区三区在线臀色熟女| 久久精品91蜜桃| 亚洲欧美精品综合久久99| 变态另类丝袜制服| 在线a可以看的网站| 干丝袜人妻中文字幕| 亚洲内射少妇av| 亚洲成人av在线免费| 夜夜夜夜夜久久久久| 亚洲av成人av| .国产精品久久| 97超级碰碰碰精品色视频在线观看| 国产一区亚洲一区在线观看| 色综合色国产| 在线看三级毛片| 久久人人爽人人爽人人片va| 黄色一级大片看看| 麻豆乱淫一区二区| 亚洲av成人av| 一级黄色大片毛片| 久久草成人影院| 综合色av麻豆| 亚洲电影在线观看av| 国产精品久久电影中文字幕| 欧美性猛交╳xxx乱大交人| 亚洲成人久久性| av黄色大香蕉| 国产高清有码在线观看视频| 男人舔奶头视频| 国产色爽女视频免费观看| 色播亚洲综合网| 亚洲自拍偷在线|