黃開(kāi)龍，林錦春，馬 盼，黃文靜，陸俊翔，唐小新，王式功

(1.成都信息工程大學(xué)大氣科學(xué)學(xué)院，高原大氣與環(huán)境四川省重點(diǎn)實(shí)驗(yàn)室，四川 成都 610225；2.深圳大學(xué)第三附屬醫(yī)院，廣東 深圳 518000；3.廣東省深圳市氣象服務(wù)中心，廣東 深圳 518000；4.廣東省汕頭市氣象局，廣東 汕頭 515041)

引 言

近年來(lái)全球氣候變化和災(zāi)害性天氣的頻發(fā)對(duì)人體健康造成了不可忽視的影響，因此關(guān)于天氣氣候變化對(duì)人群健康影響的研究日益增多[1-3]，隨著“健康中國(guó)2030”規(guī)劃綱要的提出，如何應(yīng)對(duì)全球氣候變化背景下的人群健康風(fēng)險(xiǎn)已成為當(dāng)下研究的重心之一[4]。呼吸系統(tǒng)作為人體的第一道防線，常受多種細(xì)菌病毒的侵蝕，較其他身體器官更易受到來(lái)自外界的刺激[5-6]。上呼吸道感染(簡(jiǎn)稱“上感”)屬于呼吸系統(tǒng)疾病中最易感、最常見(jiàn)的一類疾病，由于其發(fā)病狀況與天氣氣候條件密切相關(guān)，近年來(lái)備受流行病學(xué)、氣象學(xué)領(lǐng)域相關(guān)學(xué)者的關(guān)注[7-9]。研究表明，秋冬季為我國(guó)北方地區(qū)上感的主要高發(fā)期，氣溫變化是誘發(fā)上感的主要因素之一[10-13]；南方地區(qū)上感多于冬春季爆發(fā)，且相對(duì)濕度、氣壓等要素與上感的發(fā)病存在明顯的聯(lián)系[14-15]；伴隨冷空氣的活動(dòng)，上感發(fā)病人數(shù)增加明顯，老人、兒童在此期間的健康風(fēng)險(xiǎn)也有所上升[16-18]。深圳地處我國(guó)華南沿海地區(qū)，具有亞熱帶海洋性氣候特征，夏半年常受熱帶氣旋和季風(fēng)低壓的影響，冬半年受冷空氣活動(dòng)影響[19]。以往對(duì)深圳地區(qū)乃至華南地區(qū)上感疾病的研究較少[20-21]，而國(guó)內(nèi)外大多數(shù)相關(guān)研究側(cè)重于分析氣象條件對(duì)疾病死亡率或傳染病的影響[22-24]，資料多來(lái)源于疾病預(yù)防控制中心或綜合醫(yī)院，對(duì)社區(qū)人群患病情況的相關(guān)研究較少，而社區(qū)醫(yī)院多為上感等易感疾病的第一就診。本文基于深圳市獨(dú)特天氣氣候條件，選取基層醫(yī)院(社區(qū)健康服務(wù)中心)作為切入點(diǎn)，探究氣象要素與深圳地區(qū)上感發(fā)病的關(guān)系，分析疾病發(fā)病規(guī)律與我國(guó)其他地區(qū)的異同，更好地了解人群應(yīng)對(duì)天氣氣候變化的真實(shí)反饋情況，以期為當(dāng)?shù)芈毮懿块T和醫(yī)療機(jī)構(gòu)提供理論依據(jù)，從而有助于促進(jìn)深圳健康城市建設(shè)。

1 資料與方法

1.1 資 料

利用深圳市羅湖區(qū)49家社區(qū)健康服務(wù)中心2014—2018年逐日上感就診人數(shù)資料，根據(jù)國(guó)際疾病分類標(biāo)準(zhǔn)第10版(ICD-10)[25]，從中篩選出上感(ICD-10編碼：J00.0-J06.9)就診病例資料。由于可獲取資料的局限性，利用國(guó)家氣象科學(xué)數(shù)據(jù)中心(http://data.cma.cn/)提供的深圳市國(guó)家基準(zhǔn)站2014—2018年逐日氣溫、氣壓、相對(duì)濕度、風(fēng)速等氣象要素，并使用五點(diǎn)平滑方法補(bǔ)齊缺測(cè)值，該測(cè)站位于深圳市福田區(qū)，與羅湖區(qū)相鄰，中心距離約10 km。利用觀測(cè)數(shù)據(jù)計(jì)算日較差、水汽壓，具體公式[26]如下：

(1)

式中：DTR(℃)為氣溫日較差；Tmax、Tmin(℃)分別為日最高、最低氣溫；T、t(℃)分別為絕對(duì)溫度、攝氏溫度；e0(T)(hPa)為飽和水汽壓；ea(hPa)為水汽壓，采用改進(jìn)馬格納斯公式[27]計(jì)算，RH(%)為相對(duì)濕度。

1.2 方 法

為保證數(shù)據(jù)的可信性，對(duì)采集的就診數(shù)據(jù)去除錄入錯(cuò)誤、數(shù)據(jù)不完整條目，在此基礎(chǔ)上，合并同一季度出現(xiàn)的重復(fù)就診序列以規(guī)避基層醫(yī)院由于復(fù)診增加的病患數(shù)，統(tǒng)計(jì)得到上感就診人數(shù)的逐日序列資料。

分布滯后非線性模型(the distributed lag non-linear model，DLNM)著重關(guān)注不同氣象、環(huán)境要素與疾病的暴露-反應(yīng)關(guān)系，并考慮其中的滯后關(guān)系和累積效應(yīng)，對(duì)健康效應(yīng)的評(píng)估符合實(shí)際變化規(guī)律[28]，引入廣義線性模型(generalized linear models，GLM)建立變量的期望值函數(shù)(連接函數(shù))與預(yù)測(cè)變量之間的線性關(guān)系。選取與上感就診人數(shù)相關(guān)的氣象因子(通過(guò)α=0.05的顯著性檢驗(yàn))，考慮到每日上感就診人數(shù)為小概率事件，為避免可能存在的過(guò)度離散化，采用廣義泊松回歸(generalized Poisson regression)進(jìn)行模型擬合，建立基于氣象要素基變量嵌套GLM的DLNM模型。在控制節(jié)假日效應(yīng)和周末效應(yīng)的同時(shí)，確定以中位數(shù)為基準(zhǔn)分析氣象要素與上感就診人數(shù)的暴露-反應(yīng)關(guān)系[29]，并計(jì)算相對(duì)風(fēng)險(xiǎn)(relative risk，RR)作為主要的表征量。相對(duì)風(fēng)險(xiǎn)代表一個(gè)群體暴露在一定風(fēng)險(xiǎn)下與未暴露在該風(fēng)險(xiǎn)下某事件發(fā)生概率的比值，比值等于1，表明暴露與否與該事件發(fā)生無(wú)關(guān)；大于1表明暴露使該事件發(fā)生風(fēng)險(xiǎn)增加；小于1表明暴露使該事件發(fā)生風(fēng)險(xiǎn)降低。基本模型[28]可表示為：

LnE(Yt|X)=NS(Time, df)+

NS(M,df,lag,df)+holiday+DOW+β(2)

式中：Yt是第t日上感就診人數(shù)；X是引入模型的不同氣象因子；NS是自然立方樣條基函數(shù)；Time是時(shí)間變量(Time=1，2，…，1826)；df是自由度參數(shù)，利用赤池信息準(zhǔn)則(akaike’s information criterion，AIC)調(diào)整氣象要素的非線性自由度和滯后自由度；M是相應(yīng)氣象因子的交叉基矩陣；lag(d)是滯后天數(shù)；β為常數(shù)?？紤]每年節(jié)假日特別是春節(jié)期間深圳地區(qū)人口大量外流，在模型中添加啞變量holiday和DOW，包括研究期內(nèi)的所有法定節(jié)假日和周末，在建模分析中，避免了因人為主觀因素影響而造成暴露風(fēng)險(xiǎn)的誤差。在氣溫分析中，引入累積效應(yīng)作為疊加表征量，以溫度為例：氣溫累積效應(yīng)表征持續(xù)暴露在低溫(高溫)環(huán)境后冷(熱)效應(yīng)的疊加效應(yīng)，其中低溫(高溫)環(huán)境是指低于(高于)研究期10%(90%)的氣溫范疇[30]，其余氣象要素也進(jìn)行相同劃分。

利用SPSS19.0對(duì)氣象要素與上感就診人數(shù)進(jìn)行描述性統(tǒng)計(jì)及相關(guān)性分析，采用R軟件3.4.4版“dlnm”程序包進(jìn)行分布滯后非線性模型(DLNM)分析，檢驗(yàn)水平為α=0.05。

2 結(jié)果分析

2.1 基本特征

表1列出2014—2018年深圳市不同氣象要素統(tǒng)計(jì)特征。可以看出，2014—2018年深圳市日平均氣溫、相對(duì)濕度、氣壓、風(fēng)速分別為23.5 ℃、75.4%、1005.6 hPa、1.9 m·s-1；主要表征量日均氣溫在3.5～33 ℃之間，日均相對(duì)濕度在19%～100%之間。研究期內(nèi)上感就診人數(shù)共計(jì)1 395 553例，日均764.3人(表2)；按年齡段分組，未成年組(0～18歲)病案共423 433例，占比30.3%，成年組(19～59歲)病案共858 430例，占比61.5%，老年組(60歲及以上)病案共113 690例，占比8.2%。

表1 2014—2018年深圳市不同氣象要素統(tǒng)計(jì)特征Tab.1 The statistical characteristics of meteorological elements in Shenzhen from 2014 to 2018

表2 2014—2018年深圳市羅湖區(qū)上感就診人數(shù)統(tǒng)計(jì)Tab.2 Summary statistics of the number of upper respiratory tract infections (URI) visits in Luohu district of Shenzhen from 2014 to 2018 單位：人

圖1是2014—2018年深圳市羅湖區(qū)上感就診人數(shù)的日變化?？梢钥闯?，就診人數(shù)極大值出現(xiàn)在2018年1月2日(2792人)，極小值出現(xiàn)在2015年2月19日(6人)；逐日上感就診人數(shù)變化呈明顯的季節(jié)特征，春(3—5月)、夏(6—8月)、秋(9—11月)、冬(12月至次年1月)日均就診人數(shù)分別為827.6、770.5、665.3、786.3人，分別占總上感就診人數(shù)的27.24%、25.40%、21.71%、25.65%。

圖1 2014—2018年深圳市羅湖區(qū)上感就診人數(shù)日變化Fig.1 The daily variation of the number of URI visits in Luohu district of Shenzhen from 2014 to 2018

2.2 不同時(shí)間尺度上感就診人數(shù)變化

圖2為2014—2018年深圳市羅湖區(qū)逐日上感就診人數(shù)的月平均變化?？梢钥闯?，深圳市羅湖區(qū)逐日上感就診人數(shù)月平均值存在3個(gè)高峰，分別是1月、4月、7月，逐日就診人數(shù)平均值分別為1007.9、935.1、968.9人；逐日上感就診人數(shù)月平均的最低值出現(xiàn)在9月(547.3人)；總體上男性患者較女性患者多4.44%。

圖2 2014—2018年深圳市羅湖區(qū)逐日上感就診人數(shù)月平均變化Fig.2 The monthly mean variation of the number of daily URI visits in Luohu district of Shenzhen from 2014 to 2018

研究發(fā)現(xiàn)節(jié)氣變化會(huì)直接或間接影響人體，從而產(chǎn)生相應(yīng)的病理反應(yīng)[31]，因此引入24節(jié)氣分析2014—2018年不同節(jié)氣上感就診人數(shù)變化(圖3)?？梢钥闯觯细邪l(fā)病的第一高峰出現(xiàn)在冬至、小寒期間(小寒1103.6人、冬至981.0人)，對(duì)應(yīng)深圳地區(qū)初冬季節(jié)；次高峰(993.3人)出現(xiàn)在小暑期間，對(duì)應(yīng)深圳持續(xù)高溫時(shí)節(jié)；第三高峰(953.9人)出現(xiàn)在清明期間，對(duì)應(yīng)深圳地區(qū)冷熱交替的春季；最低值(501.0人)出現(xiàn)在處暑期間。男女就診人數(shù)與總就診人數(shù)的節(jié)氣變化趨勢(shì)相似，在清明、小暑、小寒期間出現(xiàn)就診高峰，各節(jié)氣男性患者均比女性患者多，最大差值出現(xiàn)在夏至期間，男性就診人數(shù)較女性高8.16%。

結(jié)合不同時(shí)間尺度的上感分布圖發(fā)現(xiàn)，深圳地區(qū)初春、初夏和初秋時(shí)節(jié)上感發(fā)病人數(shù)較少，而冬春季大部分時(shí)間和盛夏時(shí)節(jié)的上感就診人數(shù)較多，與我國(guó)北方城市上感多為秋冬季爆發(fā)[32]不同，這與深圳地區(qū)海洋性氣候特點(diǎn)關(guān)系密切，溫暖濕潤(rùn)的環(huán)境有利于病菌滋生，且時(shí)斷時(shí)續(xù)的冷空氣活動(dòng)和夏季持續(xù)高溫天氣也是上感發(fā)病的重要誘因。

2.3 氣象要素對(duì)上感就診人數(shù)的影響

表3列出2014—2018年深圳市羅湖區(qū)不同人群上感就診人數(shù)與氣象要素的Spearman相關(guān)系數(shù)?？梢钥闯觯细芯驮\人數(shù)與氣溫、相對(duì)濕度、水汽壓、日累計(jì)降水呈顯著負(fù)相關(guān)，與日均氣壓、日較差成顯著正相關(guān)。深圳地區(qū)女性上感就診人數(shù)與氣象要素的相關(guān)系數(shù)較男性更高；成年組(18～60歲)和老年組(≥60歲)上感就診人數(shù)與氣象要素變化的關(guān)系顯著，與氣溫、相對(duì)濕度、水汽壓、日累計(jì)降水呈顯著負(fù)相關(guān)，與日均氣壓、日較差成顯著正相關(guān)。

2.3.1 氣溫與上感就診人數(shù)的暴露-反應(yīng)關(guān)系

選取日均氣溫的中位數(shù)(24.8 ℃)作為參考值，并以此為分界點(diǎn)分別討論溫度小于和大于該溫度時(shí)的冷效應(yīng)和熱效應(yīng)，研究深圳市羅湖區(qū)氣溫與上感就診人數(shù)的暴露-反應(yīng)關(guān)系(圖4)，同時(shí)計(jì)算持續(xù)暴露在低溫(高溫)環(huán)境中的累積效應(yīng)。由圖4可以看出，氣溫效應(yīng)對(duì)上感發(fā)病的影響以冷效應(yīng)為主，在滯后4～5 d達(dá)到最高，熱效應(yīng)以即時(shí)反饋為主。選取冷熱效應(yīng)關(guān)鍵表征值(日均氣溫的第1、5、10、90、95、99個(gè)百分位數(shù))的相對(duì)風(fēng)險(xiǎn)剖面(圖5)進(jìn)一步分析冷熱效應(yīng)對(duì)上感發(fā)病的影響。當(dāng)日均氣溫小于24.8 ℃時(shí)，上感的發(fā)病風(fēng)險(xiǎn)隨氣溫下降逐漸升高，RR在滯后4 d達(dá)到最高(RR為1.041，95%置信區(qū)間為1.022～1.060)，冷空氣活動(dòng)對(duì)深圳地區(qū)影響周期較長(zhǎng)；高溫的影響以即時(shí)效應(yīng)為主，滯后效應(yīng)不顯著(RR<1)。

圖4 2014—2018年深圳市羅湖區(qū)氣溫與上感就診人數(shù)的暴露-反應(yīng)關(guān)系三維圖Fig.4 Three-dimensional map of exposure-response relationship between temperature and number of URI visits in Luohu district of Shenzhen from 2014 to 2018

圖5 2014—2018年深圳市羅湖區(qū)氣溫關(guān)鍵表征值的相對(duì)風(fēng)險(xiǎn)剖面Fig.5 Relative risk (RR) profile of temperature key characteristic value in Luohu district of Shenzhen from 2014 to 2018

由圖5可見(jiàn)，冷熱效應(yīng)曲線均在滯后2～4 d達(dá)到拐點(diǎn)，故分別選取9、31 ℃(日均氣溫的第1個(gè)、99個(gè)百分位數(shù))進(jìn)行累積效應(yīng)分析(表4)。總體上，低溫的累積效應(yīng)隨時(shí)間遞增，累積期為10 d時(shí)相對(duì)風(fēng)險(xiǎn)達(dá)到最高(全人群RR為1.239，95%置信區(qū)間為1.088～1.411)；而高溫的累積效應(yīng)總體不顯著，累積期為4 d時(shí)全人群RR為0.947，95%置信區(qū)間為0.833～1.076，未通過(guò)顯著性檢驗(yàn)。性別分組可見(jiàn)，男性較女性更易受冷效應(yīng)的累積影響，累積期為10 d時(shí)相對(duì)風(fēng)險(xiǎn)為1.241，95%置信區(qū)間為1.088～1.415；年齡分組可見(jiàn)，成年組(19～59歲)與老年組(≥60歲)上感發(fā)病風(fēng)險(xiǎn)受冷效應(yīng)的累積影響較大，累積期為10 d時(shí)成年組、老年組的相對(duì)風(fēng)險(xiǎn)分別為1.409、1.597，且95%置信區(qū)間分別為1.243～1.598、1.376～1.854，未成年人對(duì)冷效應(yīng)的累積響應(yīng)不顯著。

表4 2014—2018年深圳市羅湖區(qū)溫度效應(yīng)對(duì)不同人群上感就診人數(shù)的影響Tab.4 The effects of temperature on the number of URI visits in Luohu district of Shenzhen from 2014 to 2018

2.3.2 相對(duì)濕度與上感就診人數(shù)的暴露-反應(yīng)關(guān)系

選取日均相對(duì)濕度(RH)的中位數(shù)(77%)作為參考值并以此為分界點(diǎn)分別討論小于和大于該相對(duì)濕度時(shí)的低濕效應(yīng)和高濕效應(yīng)，建立模型研究相對(duì)濕度與上感就診人數(shù)的暴露-反應(yīng)關(guān)系(圖6)?？梢钥闯?，濕度影響以低濕效應(yīng)為主，當(dāng)相對(duì)濕度小于77%時(shí)，隨濕度下降上感發(fā)病逐漸增加，且RR在當(dāng)天達(dá)到峰值(RR為1.058，95%置信區(qū)間為1.049～1.068)，且在低濕條件下，在滯后0～4 d都存在上感發(fā)病的風(fēng)險(xiǎn)；在濕潤(rùn)環(huán)境下(77%≤RH≤90%)上感發(fā)病相對(duì)較少。

圖6 2014—2018年深圳市羅湖區(qū)濕度與上感就診人數(shù)的暴露-反應(yīng)關(guān)系三維圖Fig.6 Three-dimensional map of exposure-response relationship between humidity and number of URI visits in Luohu district of Shenzhen from 2014 to 2018

選取相對(duì)濕度關(guān)鍵表征值(日均相對(duì)濕度的第1、5、10、90、95、99個(gè)百分位數(shù))的相對(duì)風(fēng)險(xiǎn)剖面(圖7)分析濕度效應(yīng)對(duì)上感發(fā)病的影響。低濕效應(yīng)曲線在當(dāng)天達(dá)到峰值，而后逐漸下降并在滯后3 d后趨于水平，高濕效應(yīng)曲線在滯后1～2 d變化明顯。當(dāng)RH<77%時(shí)，RH每下降1%，上感發(fā)病的相對(duì)風(fēng)險(xiǎn)上升12.00%(95%置信區(qū)間上升9.86%～14.19%)；高濕條件(RH>90%)對(duì)上感發(fā)病的貢獻(xiàn)以滯后效應(yīng)為主，RR在滯后2 d達(dá)到高峰，RR為1.010，95%置信區(qū)間為1.006～1.015；當(dāng)RH>90%時(shí)，RH每上升1%，上感發(fā)病的相對(duì)風(fēng)險(xiǎn)上升14.36%(95%置信區(qū)間上升11.82%～16.97%)。

圖7 2014—2018年深圳市羅湖區(qū)濕度關(guān)鍵表征值的相對(duì)風(fēng)險(xiǎn)剖面圖Fig.7 Relative risk profile of humidity key characteristic value in Luohu district of Shenzhen from 2014 to 2018

由相對(duì)濕度的關(guān)鍵表征值中分別選取36%、96%(相對(duì)濕度的第1個(gè)、99個(gè)百分位數(shù))進(jìn)行濕度分層效應(yīng)研究，得到濕度效應(yīng)分層表(表5)。從性別分組可以看出，女性較男性更易受低濕效應(yīng)的影響，女性、男性即時(shí)(滯后0 d)的相對(duì)風(fēng)險(xiǎn)(RR)分別為1.082、1.064，95%置信區(qū)間分別為1.065～1.100、1.047～1.081；在高濕環(huán)境下，男性對(duì)高濕條件的耐受程度較女性高，男性、女性滯后14 d的RR分別為1.027、1.031，95%置信區(qū)間分別為1.020～1.033、1.024～1.037，隨著時(shí)間累積發(fā)病風(fēng)險(xiǎn)都有不同程度的上升。從年齡分組可以看出，在低濕環(huán)境中，成年組(16～60歲)上感即時(shí)發(fā)病的RR較其他年齡組高(RR為1.077，95%置信區(qū)間為1.061～1.092)；未成年組(0～18歲)相對(duì)較低，滯后3 d時(shí)RR為0.998，95%置信區(qū)間為0.991～1.005，說(shuō)明未成年組對(duì)持續(xù)低濕天氣的耐受程度較好；老年組(>60歲)對(duì)高濕天氣的響應(yīng)明顯，RR在滯后2 d左右達(dá)到峰值(RR為1.035，95%置信區(qū)間為1.018～1.051)，未成年組則是隨著濕度累積，發(fā)病風(fēng)險(xiǎn)愈高，當(dāng)RH>90%時(shí)，RH每上升1%，上感發(fā)病的相對(duì)風(fēng)險(xiǎn)上升36.13%(95%置信區(qū)間上升30.74%～41.77%)，青少年兒童有必要提防持續(xù)高濕天氣帶來(lái)的不利影響。

表5 2014—2018年深圳市羅湖區(qū)濕度效應(yīng)對(duì)上感就診人數(shù)的影響Tab.5 The effect of humidity on the number of URI visits in Luohu district of Shenzhen from 2014 to 2018

2.3.3 氣壓、風(fēng)速與上感就診人數(shù)的暴露-反應(yīng)關(guān)系

氣壓、風(fēng)速對(duì)上感發(fā)病的影響雖沒(méi)有溫濕變化帶來(lái)的影響直接，但在天氣過(guò)程影響時(shí)，這兩個(gè)要素的變化較為劇烈，其反饋也十分重要。分別選取日均氣壓、日均風(fēng)速的中位數(shù)(1005.6 hPa、1.8 m·s-1)為模型參考值并分析氣壓、風(fēng)速與上感就診人數(shù)的暴露-反應(yīng)關(guān)系(圖8)?？梢钥闯?，氣壓影響主要表現(xiàn)為高壓效應(yīng)，冷空氣活動(dòng)常常伴隨高氣壓過(guò)境，與上感的發(fā)病高峰有著強(qiáng)對(duì)應(yīng)關(guān)系；風(fēng)速影響以大風(fēng)效應(yīng)為主，滯后1 d的相對(duì)風(fēng)險(xiǎn)隨風(fēng)速增大呈準(zhǔn)線性增加。

圖8 2014—2018年深圳市羅湖區(qū)氣壓(a)、風(fēng)速(b)與上感就診人數(shù)的暴露-反應(yīng)關(guān)系三維圖Fig.8 Three-dimensional map of exposure-response relationship between air pressure (a), wind speed (b) and the number of URI visits in Luohu district of Shenzhen from 2014 to 2018

選取高壓和大風(fēng)關(guān)鍵表征值(日均氣壓和風(fēng)速的第90、95、99個(gè)百分位數(shù))的相對(duì)風(fēng)險(xiǎn)剖面做進(jìn)一步分析(圖9)?？梢钥闯觯?019 hPa(日均氣壓的第99個(gè)百分位數(shù))在滯后1 d時(shí)上感發(fā)病的相對(duì)風(fēng)險(xiǎn)(RR)為1.075，95%置信區(qū)間為1.065～1.085，這與冷效應(yīng)有強(qiáng)對(duì)應(yīng)關(guān)系；低壓效應(yīng)整體對(duì)上感發(fā)病的影響相對(duì)不大，但在極端低壓天氣時(shí)(日均氣壓小于991.83 hPa，日均氣壓的第1個(gè)百分位數(shù))，上感的相對(duì)風(fēng)險(xiǎn)也有所提高(圖略)。日均風(fēng)速為4.4 m·s-1(日均風(fēng)速的第99個(gè)百分位數(shù))滯后1 d時(shí)RR為1.030，95%置信區(qū)間為1.022～1.038，這與冷空氣大風(fēng)天氣有強(qiáng)對(duì)應(yīng)關(guān)系；另外，在靜穩(wěn)天氣下，即當(dāng)日均風(fēng)速小于1.1 m·s-1(日均風(fēng)速的第10個(gè)百分位數(shù))時(shí)，上感發(fā)病的風(fēng)險(xiǎn)隨時(shí)間升高，在滯后14 d達(dá)到最高，日均風(fēng)速為0.8 m·s-1(日均風(fēng)速的第1個(gè)百分位數(shù))滯后14 d時(shí)RR為1.017，95%置信區(qū)間為1.012～1.023(圖略)。

圖9 2014—2018年深圳市羅湖區(qū)高壓(a)、大風(fēng)(b)效應(yīng)對(duì)上感就診人數(shù)的相對(duì)風(fēng)險(xiǎn)剖面圖Fig.9 Relative risk profile of high-pressure (a) and strong wind (b) conditions to the number of URI visits in Luohu district of Shenzhen from 2014 to 2018

3 討論與結(jié)論

3.1 討 論

對(duì)比不同時(shí)間尺度分析結(jié)果，深圳市冬季上感的高發(fā)期與國(guó)內(nèi)相關(guān)研究趨于一致[14-15]，伴隨冷空氣南下，氣溫驟降，人群對(duì)低溫天氣有明顯響應(yīng)。與北方地區(qū)有所不同的是，春季特別是清明期間出現(xiàn)的上感高發(fā)期，可能與在此期間天氣異常變化過(guò)程有關(guān)，人群對(duì)冷熱交替變化的不適應(yīng)導(dǎo)致上感爆發(fā)[17]。夏季發(fā)病高峰也同其他地區(qū)有所差異，這可能與該地區(qū)夏季流動(dòng)人口激增有關(guān)[33]，人群接觸時(shí)間的增加、密閉的室內(nèi)空氣流通等因素[8]，都會(huì)在一定程度上加劇呼吸道疾病傳染的發(fā)生。考慮深圳地處華南沿海地區(qū)，緯度的差異造成了該地區(qū)夏季持續(xù)的高溫天氣；沿海地區(qū)汛期期間，特別是熱帶氣旋與季風(fēng)低壓影響時(shí)[19,34]，風(fēng)速激增，氣壓、氣溫驟降，從病理學(xué)角度來(lái)講，暴露在惡劣天氣和潮濕環(huán)境中，會(huì)使人體呼吸道局部溫度降低，毛細(xì)血管急劇收縮，黏膜上皮的纖毛活動(dòng)逐漸減慢，氣管排出細(xì)菌的功能越來(lái)越弱，極易誘發(fā)呼吸系統(tǒng)疾病[35-36]。

在各氣象要素建模分析中，上感發(fā)病有明顯的氣溫效應(yīng)，但暴露風(fēng)險(xiǎn)在不同地區(qū)存在差異[37-39]。在北京地區(qū)，低溫效應(yīng)的RR峰值為1.811，95%置信區(qū)間為1.229～2.667[10]，而在本研究區(qū)低溫效應(yīng)的RR峰值為1.041，95%置信區(qū)間為1.022～1.060，說(shuō)明冷空氣在南下過(guò)程中強(qiáng)度逐漸減弱，氣溫降幅減小，人群呼吸道疾病發(fā)病風(fēng)險(xiǎn)也有所降低。深圳地區(qū)濕度的發(fā)病風(fēng)險(xiǎn)與已有研究結(jié)論[40-41]相似，其曲線均呈“U”型分布，并在低濕干燥的環(huán)境下表現(xiàn)為較高的發(fā)病風(fēng)險(xiǎn)；不同的是，南方地區(qū)上感的濕度閾值較北方高[42]，北方地區(qū)上感發(fā)病的相對(duì)濕度低值區(qū)為45%～55%，而研究區(qū)發(fā)病的相對(duì)濕度低值區(qū)為70%～80%，說(shuō)明南北方人群呼吸道黏膜對(duì)空氣濕度的耐受程度不同。

在實(shí)際情況中各氣象要素的影響機(jī)制比較復(fù)雜，后續(xù)也將綜合考慮不同大氣環(huán)流背景下各氣象要素的疊加效應(yīng)，特別是在極端天氣過(guò)程期間，氣溫變幅、氣壓變幅及它們之間的相互作用也會(huì)對(duì)上感的發(fā)病產(chǎn)生影響[40-41]。上感的發(fā)病與空氣質(zhì)量狀況存在一定關(guān)聯(lián)[43-44]，本研究曾嘗試在模型中加入污染物(例如SO2、NO2)作為混雜因素，再量化氣象要素的影響，發(fā)現(xiàn)其暴露風(fēng)險(xiǎn)變化甚微，加之為避免模型過(guò)擬合，故在最終模型去除了污染變量，僅探討氣象要素誘發(fā)的相關(guān)風(fēng)險(xiǎn)。

3.2 結(jié) 論

(1)深圳市羅湖區(qū)上感發(fā)病的分布特征存在時(shí)間尺度的差異，春季上感就診人數(shù)最多，冬夏次之；24節(jié)氣分組發(fā)現(xiàn)，處暑期間上感就診人數(shù)最低，同時(shí)年內(nèi)有3個(gè)發(fā)病高峰期，分別是冬至小寒、清明與小暑期間。

(2)模型研究顯示，氣溫對(duì)于上感發(fā)病的影響以滯后效應(yīng)為主，在滯后4 d時(shí)相對(duì)風(fēng)險(xiǎn)達(dá)到最高(RR為1.041，95%置信區(qū)間為1.022～1.060)；冬春季的冷空氣活動(dòng)造成的大風(fēng)降溫天氣易引發(fā)上呼吸道疾病，同時(shí)，女性較男性、老年人較中青年人有更顯著的冷效應(yīng)反饋；濕度對(duì)上感發(fā)病的影響以低濕效應(yīng)為主，相對(duì)風(fēng)險(xiǎn)在當(dāng)天達(dá)到峰值(RR為1.058，95%置信區(qū)間為1.049～1.068)；氣壓對(duì)上感發(fā)病的影響以高壓效應(yīng)為主且RR在滯后1 d達(dá)到峰值；風(fēng)速變化對(duì)上感發(fā)病的影響以大風(fēng)效應(yīng)為主。