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

    氣溫與臭氧協(xié)同作用對(duì)COPD住院人數(shù)的影響

    2019-12-26 04:32:34付桂琴田亞芹張成偉
    中國(guó)環(huán)境科學(xué) 2019年12期
    關(guān)鍵詞:平均氣溫石家莊市氣溫

    付桂琴,田亞芹,張成偉

    氣溫與臭氧協(xié)同作用對(duì)COPD住院人數(shù)的影響

    付桂琴*,田亞芹,張成偉

    (河北省氣象服務(wù)中心,河北省氣象與生態(tài)環(huán)境重點(diǎn)實(shí)驗(yàn)室,河北 石家莊 050021)

    為評(píng)價(jià)氣溫與臭氧(O3)協(xié)同作用對(duì)慢性阻塞性肺病(COPD)住院人數(shù)的影響,采用廣義相加模型(GAM)、平滑曲線閾值效應(yīng)和飽和效應(yīng)方法,在控制了時(shí)間序列長(zhǎng)期趨勢(shì)、季節(jié)效應(yīng)、節(jié)假日效應(yīng)和相對(duì)濕度、污染物PM2.5、SO2、NO2濃度混雜因素的影響后,分析了2013~2016年石家莊市氣溫與O3協(xié)同作用對(duì)COPD住院人數(shù)影響的暴露-反應(yīng)關(guān)系,并考慮了不同年齡、性別的分人群效應(yīng).結(jié)果表明,日平均氣溫的閾值為0.3,19℃,O3的分層臨界值100μg/m3;當(dāng)O33100mg/m3時(shí),隨著氣溫每上升1℃,COPD日均住院人數(shù)增加的相對(duì)危險(xiǎn)度1.042 (95%:1.020, 1.065),并且在319℃時(shí)其效應(yīng)值最大,為1.084 (95%:1.048, 1.121).較高氣溫、高O3污染環(huán)境對(duì)COPD疾病人群的影響存在協(xié)同加強(qiáng)效應(yīng),且對(duì)老年人影響更敏感,而男女性別上沒(méi)有顯著差異.

    慢性阻塞性肺病(COPD);氣溫和臭氧;協(xié)同作用;廣義相加模型(GAM)

    國(guó)內(nèi)外流行病學(xué)研究發(fā)現(xiàn),除細(xì)顆粒物PM2.5污染對(duì)人體健康不利外,臭氧(O3)污染與呼吸系統(tǒng)疾病發(fā)病率和死亡率增加有關(guān)[1-4].O3污染會(huì)造成呼吸系統(tǒng)的感染及損害,全球每年因O3污染引起的過(guò)早死高達(dá)70000人[5].Tao等[6]研究珠江三角洲O3污染的急性影響發(fā)現(xiàn),O3濃度每增加10mg/m3,總死亡率將增加0.81%,對(duì)呼吸系統(tǒng)死亡率的影響更強(qiáng).胡悅等[7]研究大氣O3濃度與居民急救關(guān)系發(fā)現(xiàn),O3濃度每升高10mg/m3,居民因呼吸系統(tǒng)疾病急救人次增加1.21%.

    與此同時(shí),全球氣候變暖,極端天氣事件增多,天氣氣候變化對(duì)人體健康產(chǎn)生的影響受到人們高度關(guān)注[8-11].大量研究表明,氣溫與呼吸疾病發(fā)病率或死亡率呈“J”型或“U”型關(guān)系[12-14],當(dāng)溫度超過(guò)了人體的承受極限(閾值溫度)時(shí),容易導(dǎo)致疾病的發(fā)生或加重,甚至死亡[15].

    氣溫和臭氧都與呼吸系統(tǒng)疾病有關(guān),但以往研究氣溫對(duì)呼吸系統(tǒng)疾病的影響時(shí),往往把大氣污染物作為混雜因素加以排除,而在研究污染物對(duì)呼吸系統(tǒng)疾病的影響時(shí),又把氣溫作為混雜因素來(lái)考慮[16].事實(shí)上,氣象條件與大氣污染存在著一定的關(guān)聯(lián),尤其是O3濃度和溫度存在顯著的正相關(guān),且在不同的溫度范圍,O3表現(xiàn)出不同的增長(zhǎng)速率[17].氣溫與O3對(duì)呼吸系統(tǒng)疾病影響的協(xié)同效應(yīng)值得關(guān)注.本文選擇O3污染和夏季高溫相對(duì)集中的石家莊,利用2013~2016年石家莊市慢性阻塞性肺病(COPD)住院醫(yī)保報(bào)銷(xiāo)病例資料,分析氣溫與O3協(xié)同作用對(duì)慢性阻塞性肺病住院人數(shù)的影響關(guān)系,為確立科學(xué)的預(yù)防和醫(yī)療技術(shù)提供參考,為提高敏感人群的主動(dòng)預(yù)防,推動(dòng)由疾病治療向健康管理的轉(zhuǎn)變提供科技支撐.

    1 資料與方法

    1.1 資料來(lái)源

    疾病資料:2013年1月1日~2016年12月31日石家莊市COPD住院病例資料來(lái)源于石家莊市醫(yī)療保險(xiǎn)管理中心,包括:患者年齡、性別、住院日期、出院日期、病種代碼、醫(yī)療費(fèi)用等信息.

    氣象資料由河北省氣象信息中心提供,選取2013年1月1日~2016年12月31日石家莊站地面觀測(cè)資料,包括日平均氣溫、相對(duì)濕度、氣壓等.同期的大氣污染物O3-8h濃度和PM2.5、SO2、NO2日均濃度監(jiān)測(cè)資料,來(lái)源于河北省環(huán)境保護(hù)局.

    1.2 研究方法

    COPD每日住院人數(shù)按時(shí)間序列統(tǒng)計(jì),通過(guò)泊松分布的廣義相加模型(GAM),分析氣溫與O3協(xié)同作用對(duì)COPD住院人數(shù)的影響.利用R3.4.3軟件和EmpowerStats(http://empowerstats.com/en/)軟件實(shí)現(xiàn)[2].

    在研究氣溫與O3協(xié)同對(duì)COPD住院人數(shù)影響之前,首先,分別研究日平均氣溫、O3濃度對(duì)COPD住院人數(shù)的影響關(guān)系,確定日平均氣溫、O3對(duì)COPD的影響是否存在曲線關(guān)系.以COPD日住院人數(shù)為響應(yīng)變量,以氣溫為暴露變量,控制了時(shí)間長(zhǎng)期趨勢(shì)(Time)、季節(jié)(Season)、節(jié)假日(Holiday)、相對(duì)濕度(RH)、污染物PM2.5、SO2、NO2濃度混雜因素的影響后,得到的影響關(guān)系式:

    式中:Y為第日COPD住院人數(shù);(Y|)為第日COPD住院人數(shù)的期望值;為截距,為回歸系數(shù);為非線性樣條函數(shù),df是自由度;使用自相關(guān)函數(shù)(PACF),通過(guò)最小信息準(zhǔn)則法(Akaike Information Criterion, AIC),確定各因子自由度.

    其次,按照平滑曲線閾值效應(yīng)和飽和效應(yīng)方法,分析日平均氣溫、O3濃度的分層閾值點(diǎn),用對(duì)數(shù)似然比進(jìn)行檢驗(yàn).第三,分析氣溫與O3對(duì)COPD住院人數(shù)的協(xié)同效應(yīng),以<0.05作為是否具有統(tǒng)計(jì)學(xué)意義進(jìn)行檢驗(yàn).

    1.3 質(zhì)量控制

    COPD住院病例資料來(lái)源于石家莊市醫(yī)療保險(xiǎn)管理中心,所有資料都有醫(yī)院診斷病種代碼記錄.按照國(guó)家疾病分類(lèi)標(biāo)準(zhǔn)第10版(IDC-10)編碼,篩選疾病編碼為J43-44的信息.空氣污染數(shù)據(jù)和氣象資料均來(lái)自國(guó)家認(rèn)證的大氣自動(dòng)監(jiān)測(cè)系統(tǒng)和氣象觀測(cè)系統(tǒng),沒(méi)有缺測(cè)記錄.

    2 結(jié)果

    2.1 石家莊市COPD住院人數(shù)與氣象環(huán)境要素描述性分析

    表1 2013~2016年石家莊市COPD日住院人數(shù)、氣象要素及空氣污染物描述性分析

    注:O3-8h為8h平均;PM2.5、SO2、NO2均為24h平均;P25、P50、P75分別表示第25、50、75百分位數(shù).

    如表1所示,COPD日均住院總?cè)藬?shù)4.4人次,日最多住院人數(shù)為16人.對(duì)應(yīng)期間的日平均氣溫14.7℃,日平均氣溫范圍為-9.4~35.5℃,季節(jié)變化差異顯著.2013~2016年O3-8h平均濃度為81.5mg/m3, O3-8h濃度范圍為2.0~262.0mg/m3,年內(nèi)變化差異顯著.PM2.5濃度平均為114.8mg/m3,SO2和NO2日均濃度分別是63.8,56.9mg/m3.

    2.2 日平均氣溫、O3與COPD住院人數(shù)的暴露—反應(yīng)關(guān)系

    實(shí)線代表住院人數(shù)隨暴露變量的變化,虛線代表其95%可信區(qū)間

    在控制了時(shí)間長(zhǎng)期變化趨勢(shì)、季節(jié)、節(jié)假日效應(yīng)、相對(duì)濕度和空氣污染物PM2.5、SO2、NO2濃度混雜影響后,分析石家莊市日平均氣溫、O3與COPD住院人數(shù)的暴露反應(yīng)關(guān)系(圖1).由圖1可以看出,日平均氣溫與COPD日均住院人數(shù)呈寬口的“U”型變化,當(dāng)日平均氣溫在0.3℃£<19℃之間時(shí),COPD日均住院人數(shù)相對(duì)較少,呈現(xiàn)在寬“U”型底部;當(dāng)日平均氣溫319℃,隨著日平均氣溫每上升1℃, COPD住院人數(shù)將增加4.3%(95%:1.024, 1.062);當(dāng)日平均氣溫<0.3℃時(shí),隨著日平均氣溫每下降1℃, COPD住院人數(shù)將增加1.5%(95%:0.952, 1.019).可見(jiàn),日平均氣溫對(duì)COPD住院人數(shù)的影響存在冷效應(yīng)(<0.3℃)、熱效應(yīng)(319℃)和相對(duì)適宜的氣溫區(qū)間(0.3£<19℃).而O3濃度與COPD日均住院人數(shù)總體上呈波動(dòng)變化趨勢(shì),隨著臭氧O3濃度每升高10mg/m3,COPD日均住院人數(shù)增加1.1%(95%:1.005, 1.017).

    2.3 不同氣溫環(huán)境下O3濃度的變化

    在控制了時(shí)間長(zhǎng)期變化趨勢(shì)、季節(jié)變化影響后,圖2是2013~2016年石家莊O3濃度隨日平均氣溫變化關(guān)系圖.由圖2可見(jiàn),隨著氣溫的升高,O3濃度總體呈增長(zhǎng)趨勢(shì),但其增長(zhǎng)速度不同.根據(jù)平滑曲線閾值效應(yīng)和飽和效應(yīng)分析,日平均氣溫以19℃分層,得到O3濃度隨日平均氣溫的關(guān)系:當(dāng)<19℃,隨著氣溫每升高1℃,O3濃度增加的相對(duì)危險(xiǎn)度RR(95%)為2.897(95%: 2.119,3.674);當(dāng)319℃時(shí),隨著日平均氣溫每上升1℃,O3濃度增加的相對(duì)危險(xiǎn)度RR(95%)為9.071(95%: 7.928,10.215),對(duì)數(shù)似然比檢驗(yàn),< 0.001,具有統(tǒng)計(jì)意義.

    圖2 石家莊市日平均氣溫與O3濃度的變化關(guān)系

    實(shí)線是O3隨氣溫變化,虛線是其95%可信區(qū)間

    在上述分析的基礎(chǔ)上,對(duì)2013~2016年石家莊市日平均氣溫和O3濃度進(jìn)行統(tǒng)計(jì),結(jié)果發(fā)現(xiàn),石家莊市日平均氣溫<19℃時(shí),89.1%的日數(shù)O3平均濃度<100mg/m3;石家莊市319℃,主要出現(xiàn)在3~10月,其中5~9月出現(xiàn)日數(shù)占到91%;O3濃度O33100mg/m3也出現(xiàn)3~10月,其中5~9月出現(xiàn)日數(shù)占到81.9%.按照中國(guó)空氣質(zhì)量標(biāo)準(zhǔn)(GB 3095-2012)[18]和環(huán)境空氣質(zhì)量指數(shù)技術(shù)規(guī)定(HJ 633-2012)[19],當(dāng)O3濃度小于100mg/m3,空氣質(zhì)量分指數(shù)為一級(jí),空氣質(zhì)量?jī)?yōu),對(duì)人體健康無(wú)影響,因此,按照O3濃度100mg/m3進(jìn)行分層分析.

    2.4 不同O3污染環(huán)境下氣溫對(duì)COPD住院人數(shù)的影響

    由表3可以看出,兩個(gè)模型不同的O3分層中,日平均氣溫與COPD住院人數(shù)的影響存在穩(wěn)定一致的影響關(guān)系.模型Ⅱ中,對(duì)于石家莊市COPD住院總?cè)藬?shù),當(dāng)O3<100mg/m3,日平均氣溫每上升1 ℃,COPD住院總?cè)藬?shù)的相對(duì)危險(xiǎn)度為1.014 (95%CI:1.002, 1.026),當(dāng)O33100mg/m3時(shí),日平均氣溫每上升1℃, COPD住院總?cè)藬?shù)的相對(duì)危險(xiǎn)度為1.042(95%:1.020, 1.065).可見(jiàn),在O33100mg/m3時(shí),對(duì)石家莊COPD住院人數(shù)增加的風(fēng)險(xiǎn)是O3< 100mg/m3時(shí)的3倍,交互作用檢驗(yàn)為0.042 (<0.05),說(shuō)明兩層間的差異具有統(tǒng)計(jì)意義.也就是說(shuō),在較高O3濃度下,氣溫對(duì)COPD住院人數(shù)增加的風(fēng)險(xiǎn)相對(duì)較低的O3濃度更顯著.從性別和年齡分人群看,同樣表現(xiàn)出高O3(O33100mg/ m3)下對(duì)COPD住院人數(shù)的影響風(fēng)險(xiǎn)大于較低O3濃度(O3<100mg/m3)下的風(fēng)險(xiǎn),只是男性和60歲以上的老人,分層間的值都小于0.05,具有統(tǒng)計(jì)意義,而女性和18~60歲的成年人,由于樣本人數(shù)少,分層間的顯著性檢驗(yàn),不具有統(tǒng)計(jì)意義.

    表3 不同O3濃度下日平均氣溫對(duì)COPD住院人數(shù)影響的相對(duì)危險(xiǎn)度及95%可信區(qū)間(95%CI)

    注:為交互作用檢驗(yàn)值;模型Ⅰ是單因素模型,模型Ⅱ是調(diào)整了時(shí)間變化趨勢(shì)、季節(jié)效應(yīng)、節(jié)假日效應(yīng),控制了相對(duì)濕度、PM2.5,NO2和SO2污染物混雜因素影響后的模型;*代表<0.05.

    2.5 氣溫與O3協(xié)同作用對(duì)COPD住院人數(shù)的影響

    在模型Ⅱ的基礎(chǔ)上,按照O3不同濃度水平分層,計(jì)算了氣溫冷、熱效應(yīng)及相對(duì)舒適區(qū)間內(nèi),日平均氣溫每上升1℃對(duì)石家莊分人群COPD住院人數(shù)的定量關(guān)系(表4).在O3<100mg/m3分層,不同氣溫區(qū)間的效應(yīng)有所差異,其中當(dāng)<0.3℃時(shí),隨著氣溫每上升1℃對(duì)COPD住院人數(shù)影響的相對(duì)危險(xiǎn)度小于1,說(shuō)明發(fā)病風(fēng)險(xiǎn)隨著氣溫的升高而降低,對(duì)總?cè)藬?shù)其相對(duì)危險(xiǎn)度0.986(95%:0.946, 1.028),男性1.003 (95%: 0.953, 1.054)、女性 0.946 (95%:0.878, 1.018)、成年人0.994 (95%:0.902, 1.096)、老人0.982 (95%:0.938, 1.028).在O33100mg/m3分層,除<0.3℃時(shí)樣本數(shù)只有5個(gè)沒(méi)有計(jì)算RR(95%)外,其他2個(gè)氣溫閾值區(qū)間,平均氣溫每上升1℃對(duì)COPD住院人數(shù)影響的相對(duì)危險(xiǎn)度都大于1,說(shuō)明發(fā)病有增加的風(fēng)險(xiǎn),當(dāng)319℃時(shí),無(wú)論是總?cè)藬?shù)還是男女,成人和老人,都表現(xiàn)出最大的效應(yīng)值,即隨著氣溫每上升1℃,COPD疾病住院總?cè)藬?shù)增加的風(fēng)險(xiǎn)為8.2%、男性 7.9%、女性8.1%、成人3.8%、老人8.4%,60歲以上的老人是最敏感人群,男女間差異不顯著.

    表4 不同O3濃度、不同氣溫區(qū)間下,日平均氣溫每上升1℃對(duì)石家莊市COPD住院人數(shù)影響的相對(duì)危險(xiǎn)度及95%可信區(qū)間(95%CI)

    注:NA表示樣本量少,沒(méi)有計(jì)算RR(95%).*代表<0.05.

    3 討論

    本研究發(fā)現(xiàn),日平均氣溫對(duì)COPD住院人數(shù)的影響呈非線性關(guān)系,日平均氣溫的閾值為0.3,19℃,存在冷效應(yīng)、熱效應(yīng)和相對(duì)適宜的氣溫區(qū)間.當(dāng)氣溫在0.3~19℃之間時(shí),COPD住院人數(shù)相對(duì)最少,為相對(duì)適宜的氣溫區(qū)間,當(dāng)<0.3℃時(shí),日平均氣溫每下降1℃,COPD住院人數(shù)增加的風(fēng)險(xiǎn)為1.5%,當(dāng)319℃時(shí),隨著日平均氣溫每上升1℃,COPD住院人數(shù)增加的風(fēng)險(xiǎn)為4.3%.研究表明高溫、低溫對(duì)哮喘呼吸疾病的急性發(fā)作都有影響[20].Gómez等[21]發(fā)現(xiàn)低溫對(duì)人群死亡的影響明顯大于高溫效應(yīng),其比值比OR值為3.40(95%:2.95,3.93).高溫能引起呼吸系統(tǒng)疾病入院率升高,中東和北歐城市的氣溫從臨界值開(kāi)始每升高1℃,呼吸系統(tǒng)疾病入院率增加4.5%[22],本研究結(jié)果與前期相關(guān)報(bào)道相一致.但不同地區(qū)、不同呼吸疾病敏感人群其最佳氣溫閾值有所不同[13,23].高溫條件下,COPD患者可能出現(xiàn)過(guò)度呼吸,引起肺過(guò)度充氣,導(dǎo)致疾病復(fù)發(fā)[24].有關(guān)低溫與呼吸系統(tǒng)疾病關(guān)系方面,有研究表明冷空氣對(duì)呼吸道黏膜有刺激作用,低溫環(huán)境,肺功能降低,炎癥細(xì)胞數(shù)量容易增加誘發(fā)疾病發(fā)生[25-26].

    與氣溫相比,國(guó)內(nèi)外關(guān)于O3污染對(duì)呼吸疾病的影響研究相對(duì)較少.本研究發(fā)現(xiàn)氣溫與O3呈非線性關(guān)系,O3對(duì)COPD住院人數(shù)的影響呈波動(dòng)變化.研究表明,長(zhǎng)期的臭氧暴露與呼吸原因引起的死亡風(fēng)險(xiǎn)相關(guān),O3濃度增加10′10-3mg/m3相關(guān)的呼吸原因?qū)е滤劳龅南鄬?duì)風(fēng)險(xiǎn)為1.04%[27],與本研究結(jié)果相一致.Faridi等[28]發(fā)現(xiàn)長(zhǎng)期接觸O3導(dǎo)致呼吸系統(tǒng)疾病死亡率增加0.9%~2.3%.Wong等[29]發(fā)現(xiàn)中國(guó)城市NO2,SO2,PM10和O34種污染物對(duì)人群健康的綜合影響等于或大于西方工業(yè)國(guó)家.Tao等[6]研究珠江三角洲O3污染的急性影響發(fā)現(xiàn),O3濃度增加10mg/m3,總死亡率將增加0.81%,對(duì)呼吸系統(tǒng)死亡率的估計(jì)影響更強(qiáng).以上研究均表明O3與呼吸原因?qū)е碌募膊 ⑺劳鲲L(fēng)險(xiǎn)有關(guān),與本研究結(jié)果一致.

    進(jìn)一步分析發(fā)現(xiàn), O3與氣溫對(duì)COPD疾病的影響具有協(xié)同加強(qiáng)作用.當(dāng)O33100mg/m3時(shí),氣溫不分層時(shí),氣溫每上升1℃對(duì)COPD日均住院總?cè)藬?shù)增加的相對(duì)危險(xiǎn)度1.042 (95%:1.020, 1.065);當(dāng)319℃時(shí),其效應(yīng)值最大為1.084 (95%:1.048, 1.121).說(shuō)明較高氣溫、較高O3污染環(huán)境下,對(duì)COPD發(fā)病風(fēng)險(xiǎn)有增加效應(yīng).有關(guān)氣象環(huán)境因素對(duì)呼吸疾病的協(xié)同作用分析方面,王敏診[23]分析了氣溫與濕度的交互作用對(duì)呼吸疾病的影響,發(fā)現(xiàn)低溫低濕的效應(yīng)最強(qiáng);張瑩等[16]分析了氣溫與PM2.5的協(xié)同作用,表明高溫與PM2.5對(duì)人群健康的影響存在協(xié)同加強(qiáng)效應(yīng).Ren等[30]研究氣溫與O3對(duì)美國(guó)東部60個(gè)社區(qū)總死亡率的短期影響,指出高溫可以調(diào)節(jié)O3對(duì)死亡的風(fēng)險(xiǎn),并且出現(xiàn)一定的區(qū)域差異.胡悅等[7]研究表明夏秋季大氣O3污染對(duì)增加居民呼吸系統(tǒng)疾病的風(fēng)險(xiǎn)更大,從季節(jié)變化的角度印證了氣溫與O3污染存在協(xié)同效應(yīng),與本文的結(jié)論相一致.而Zhang等[31]研究上海O3與包括自然原因死亡率和因心腦血管、呼吸系統(tǒng)疾病死亡率之間的關(guān)系時(shí),發(fā)現(xiàn)在我國(guó)上海,O3污染對(duì)呼吸疾病死亡率的影響在寒冷季比暖季更強(qiáng),可見(jiàn)氣溫與O3的協(xié)同性具有地域差異.

    本文分析了氣溫與O3對(duì)COPD人群的協(xié)同作用,開(kāi)展了男女性別、不同年齡段分人群分析及分層交互作用檢驗(yàn).研究方法上,與王敏診等[23]、張瑩等[16]相同,都采用的廣義相加模型.本研究的優(yōu)點(diǎn)是利用平滑曲線閾值效應(yīng)和飽和效應(yīng)方法,得到的閾值點(diǎn)相對(duì)與百分位數(shù)[32]、中位數(shù)[33]、多氣溫分層[16]等得到的閾值點(diǎn)更精確;O3不同分層間,進(jìn)行了似然比顯著性檢驗(yàn);另外本研究選擇的是我國(guó)北方污染嚴(yán)重的石家莊市,典型的呼吸類(lèi)疾病——慢性阻塞性肺病,且是市醫(yī)保住院病例資料,相對(duì)于從某個(gè)醫(yī)院得到的資料更具有代表性.所不足的是氣象和大氣污染資料是本市監(jiān)測(cè)均值,在個(gè)體暴露方面可能有一定的局限性.

    4 結(jié)論

    4.1 平均氣溫和臭氧污染對(duì)石家莊COPD住院人數(shù)均有影響.日平均氣溫的影響具有冷效應(yīng)(<0.3℃)、熱效應(yīng)(319℃)和相對(duì)適宜的氣溫環(huán)境(0.3~19℃);O3的影響沒(méi)有最低閾值.

    4.2 當(dāng)O3濃度大于100mg/m3,氣溫大于等于19℃時(shí),氣溫與O3協(xié)同加強(qiáng)效應(yīng)最強(qiáng),且對(duì)老年人影響更敏感;而當(dāng)氣溫小于0.3℃時(shí),沒(méi)發(fā)現(xiàn)氣溫冷效應(yīng)與臭氧濃度對(duì)COPD住院人數(shù)的聯(lián)合作用.

    [1] 張 瑩,邵 毅,王式功,等.北京市空氣污染物對(duì)呼吸系統(tǒng)疾病門(mén)診人數(shù)的影響 [J]. 中國(guó)環(huán)境科學(xué), 2014,34(9):2401-2407. Zhang Y, Shao Y, Wang S G, et al. Relationship between air pollutant and respiratory diseases hospital outpatient visits in Beijing [J]. China Environmental Science, 2014,34(9):2401-2407.

    [2] Fu G Q, Jiang Y F, LiuL P, et al. Effects of PM2.5exposure in different air quality grades on daily outpatient visits for childhood asthma in Shijiazhuang, China [J]. Biomed Environ Sci, 2018,31(12):888-892.

    [3] 陳仁杰,陳秉衡,闞海東.我國(guó)113個(gè)城市大氣顆粒物的健康經(jīng)濟(jì)學(xué)評(píng)價(jià) [J]. 中國(guó)環(huán)境科學(xué), 2010,30(3):410-415. Chen R J, Chen B H, Kan H D. A health-based economic assessment of particulate air pollution in 113chinese cities [J]. China Environmental Science, 2010,30(3):410-415.

    [4] Gryparis A, Forsberg B, Katsouyanni K, et al. Acute effects of ozone on mortality from the “air pollution and health [J]. American Journal of Respiratory and Critical Care Medicine, 2004,170(10):1080-1087.

    [5] Eremenko M, Dufour G, Foret G, et al. Tropospheric ozone distributions over Europe during the heat wave in July 2007 observed from infrared nadir spectra recorded by IASI [J]. Geophysical Research Letters, 2008,35(18):60-74.

    [6] Tao Y B, Huang W, Huang X L, et al. Estimated acute effects of ambient ozone and nitrogen dioxide on mortality in the pearl river delta of southern China [J]. Environmental Health Perspectives, 2011, 120(3):393-398.

    [7] 胡 悅,郭 昳,關(guān)茗洋,等.石家莊市大氣臭氧濃度與居民急救關(guān)系的時(shí)間序列分析 [J]. 環(huán)境與健康雜志, 2016,33(10):872-875. Hu Y, Guo Y, Guan M Y, et al. Time-series analysis of association between ozone concentration and dailyemergency ambulance dispatches in Shijiazhuang [J]. Environ Health, 2016,33(10):872-875.

    [8] Anderson B G, Bell M L. Weather-related mortality [J]. Epidemiology, 2009,20(2):205-213.

    [9] Baccini M, Biggeri A, Accetta G, et al. Heat effects on mortality in 15European cities [J]. Epidemiology, 2008,19(5):711-719.

    [10] Basu R. High ambient temperature and mortality: a review of epidemiologic studies from 2001 to 2008 [J]. Environmental Health, 2009,8(1):1-13.

    [11] Filleul L, Cassadou S, Medina S, et al. The relation between temperature, ozone, and mortality in nine French cities during the heat wave of 2003 [J]. Environ Health Perspect, 2006,114(9):1344-1347.

    [12] Chung J Y, Honda Y, Hong Y C, et al. Ambient temperature andmortality: an international study in four capital cities of East Asia [J]. Sci Total Environ, 2009,408(2):390-396.

    [13] Kan H D, Jia J, Chen B H. Temperature and daily mortality inshanghai: a time-series study [J]. Biomed Environ Sci, 2003,16(2):133-139.

    [14] Liu Y, Kan H, Xu J, et al. Temporal relationship between hospital admissions for pneumonia and weather conditions in Shanghai, China: a time-series analysis [J]. BMJ Open, 2014,4(7):e004961.

    [15] Knowlton K, Rotkin-Ellman M, King G, et al. The 2006california heat wave: impacts on hospitalizations and emergency department visits [J]. Environmental Health Perspectives, 2009,117(1):61-67.

    [16] 張 瑩,王式功,賈旭偉,等.氣溫與PM2.5協(xié)同作用對(duì)疾病急診就診人數(shù)的影響[J]. 中國(guó)環(huán)境科學(xué), 2017,37(8):3175-3182. Zhang Y, Wang S G, Jia X W, et al. Synergetic effect of mean temperature and PM2.5on emergency room visits for different diseases [J]. China Environmental Science, 2017,37(8):3175-3182.

    [17] 蒲 茜.長(zhǎng)三角高溫?zé)崂似陂g臭氧的形成機(jī)理及輻射效應(yīng)研究[D]. 南京:南京大學(xué), 2017. Pu X. Analysis of elevated ozone formation mechanism and radiative effect during heat waves in Yangtze river delta [D]. Nanjing: Nanjing University, 2017.

    [18] GB 3095-2012 中國(guó)空氣質(zhì)量標(biāo)準(zhǔn)[S].GB 3095-2012 Ambient air quality standards [S].

    [19] HJ 633-2012 環(huán)境空氣質(zhì)量指數(shù)(AQI)技術(shù)規(guī)定(試行) [S].HJ 633-2012 Technical regulation on ambient quality index (on trial) [S].

    [20] Xu Z, Huang C, Hu W, et al. Extreme temperatures and emergency department admissions for childhood asthma in Brisbane, Australia [J] .Occup Environ Med, 2013,70(10):730-5.

    [21] Gómez-Acebo I, Dierssen-Sotos T, Llorca J. Effect of coldtemperatures on mortality in cantabria (northern Spain): acase– crossover study [J]. Public health, 2010,124(7):398-403.

    [22] Michelozzi P, Accetta G, Sario M, et al. High temperature andhospitalizations for cardiovascular and respiratory causes in 12europeancities [J]. Am J Respir Crit Care Med, 2009,179:383-389.

    [23] 王敏珍,鄭 山,王式功,等.氣溫與濕度的交互作用對(duì)呼吸系統(tǒng)疾病的影響 [J]. 中國(guó)環(huán)境科學(xué), 2016,36(2):581-588. Wang M Z, Zhen S, Wang S G, et al. Interaction of temperature and relative humidity on emergency room visits for respiratory diseases [J]. China Environmental Science, 2016,36(2):581-588.

    [24] Viegi G, Pistelli F, Sherrill DL, et al. Definition, epidemiology and natural history of COPD [J]. EurRespir J, 2007,30:993-1013.

    [25] Larsson K, Tornling G, Gavhed D, et al. Inhalation of cold airincreases the number of inflammatory cells in the lungs in healthysubjects [J]. European Respiratory Journal, 1998,12(4):825-830.

    [26] Koskela H, Tukiainen H, Kononoff A, et al. Effect ofwhole-body exposure to cold and wind on lung function in asthmaaticpatients [J]. Chest, 1994,105(6):1728–1731.

    [27] Turner M C, Jerrett M, Pope C, et al. Long-term ozone exposure and mortality in a large prospective study [J]. Am J Respir Crit Care Med, 2016,193(10):1134-1142.

    [28] Faridi S, Shamsipour M, Krzyzanowski M, et al. Long-term trends and health impact of PM and O in Tehran, Iran, 2006~2015. [J]. Environ Int, 2018,114:37-49.

    [29] Wong C M, Vichit-Vadakan N, Kan H, et al. Public health and air pollution in Asia (PAPA): a multicity study of short term effects of air pollution on mortality [J]. Environmental Health Perspectives, 2008, 116(9):1195-1201.

    [30] Ren C, Williams G M, Mengersen K, et al. Does temperaturemodify short-term effects of ozone on total mortality in 60largeeastern US communities?-an assessment using the NMMAPSdata [J]. Environment International, 2008,34(4):451-458.

    [31] Zhang Y H, Huang W, London S J, et al. Ozone and daily mortality in shanghai, China [J]. Environmental Health Perspectives, 2006,114(8): 1227-1232.

    [32] Li G, Zhou M, Cai Y, et al. Does temperature enhance acutemortality effects of ambient particle pollution in tianjincity, China [J]. Sci. Total Environ., 2011,409:1811-1817.

    [33] Ebi K L, Exuzides K A, Lau E, et al. Weather changes associatedwith hospitalizations for cardiovascular diseases and stroke in California, 1983~1998 [J]. International Journalof Biometeorology, 2004,49(1): 48-58.

    Synergistic effect of temperature and O3on the numbers of COPD hospitalizations.

    FU Gui-qin*, TIAN Ya-qin, ZHANG Cheng-wei

    (Hebei Meteorological Service Center, Key Laboratory for Meteorology and Ecological Environment of Hebei Province, Shijiazhuang 050021, China)., 2019,39(12):5312~5318

    In this study, we applied the generalized additive models (GAMs) including threshold and saturation effect of smooth curve analysis model to evaluate the synergistic effects of air temperature and ozoneon the numbers of chronic obstructive pulmonary disease (COPD) hospitalizations in Shijiazhuang, China. After controlling the potential confounding factors, such as long-term time trends, seasonal patterns, and the effect of holidays and day of the week of time series data, we estimated the synergistic effect of temperature and ozone on COPD hospitalizations, controlling forrelative humidity, PM2.5, SO2, and NO2. In addition, we also evaluated differential associations on samples stratified by sex and age. The results showed that, the thresholds of temperature and ozone in the dose-response relationship with the numbers of COPD hospitalizations were 0.3℃, 19℃, and 100μg/m3, respectively. When the concentration of ozone was above 100μg/m3, the relative risk (RR) of COPD hospitalization was 1.042 (95% confidence interval [CI]: 1.020, 1.065) per1℃ increase in temperature, and achieved the maximum of 1.082 (95% CI: 1.049, 1.116) when temperature was above 19℃. There were synergistic effects of high temperature and high ozone concentration on the COPD hospitalizations, to which the elderly were more sensitive. However, no significant difference was found regarding gender.

    chronic obstructive pulmonary disease (COPD);temperature;ozone;synergistic effect;generalized additive model (GAM)

    X503.1

    A

    1000-6923(2019)12-5312-07

    付桂琴(1968-),女,河北蔚縣人,正高級(jí)工程師,碩士,主要從事環(huán)境健康等應(yīng)用氣象研究.發(fā)表論文30余篇.

    2019-05-21

    河北省重點(diǎn)研發(fā)項(xiàng)目(18273705D);上海健康氣象重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金(QXJK201606)

    *責(zé)任作者, 正高級(jí)工程師, fgq84@tom.com

    猜你喜歡
    平均氣溫石家莊市氣溫
    氣溫的神奇“魔法”
    石家莊市第二中學(xué)
    石家莊市井陘礦區(qū)賈莊村
    公民與法治(2022年1期)2022-07-26 05:57:40
    基于FY-3D和FY-4A的氣溫時(shí)空融合
    石家莊市
    石家莊市
    烏蘭縣近38年氣溫變化特征分析
    從全球氣候變暖大背景看萊州市30a氣溫變化
    1981—2010年拐子湖地區(qū)氣溫變化特征及趨勢(shì)分析
    近50年來(lái)全球背景下青藏高原氣候變化特征分析
    欧美激情久久久久久爽电影| 亚洲人成伊人成综合网2020| 亚洲精华国产精华精| 18禁裸乳无遮挡免费网站照片| 最近在线观看免费完整版| 成人亚洲精品av一区二区| 小蜜桃在线观看免费完整版高清| 三级男女做爰猛烈吃奶摸视频| 国产精品1区2区在线观看.| 在线播放国产精品三级| 在线a可以看的网站| 亚洲男人的天堂狠狠| 老司机福利观看| 精品久久久久久久久久久久久| 久久精品国产亚洲av天美| 国产综合懂色| 中文在线观看免费www的网站| 国产av麻豆久久久久久久| 欧美+日韩+精品| 国产探花极品一区二区| 看免费av毛片| 在线播放无遮挡| 黄色日韩在线| 成年女人看的毛片在线观看| 亚洲精品在线观看二区| 国产精品久久久久久精品电影| 亚洲综合色惰| 午夜a级毛片| 久久精品国产自在天天线| 亚洲美女视频黄频| 精品一区二区三区视频在线观看免费| 听说在线观看完整版免费高清| 亚洲国产欧洲综合997久久,| 深夜a级毛片| 日韩欧美国产一区二区入口| 一进一出抽搐gif免费好疼| 免费av观看视频| 丁香六月欧美| 亚洲国产精品合色在线| 免费在线观看影片大全网站| 国产av在哪里看| 我的女老师完整版在线观看| 九九热线精品视视频播放| 全区人妻精品视频| 日日摸夜夜添夜夜添av毛片 | 欧美乱色亚洲激情| 日韩大尺度精品在线看网址| 亚洲人成网站高清观看| 国内精品久久久久久久电影| 一进一出好大好爽视频| 国产野战对白在线观看| 又粗又爽又猛毛片免费看| 欧美最新免费一区二区三区 | 亚洲精品在线观看二区| av天堂在线播放| eeuss影院久久| 国产成人福利小说| 夜夜爽天天搞| 最近最新中文字幕大全电影3| 伦理电影大哥的女人| 少妇被粗大猛烈的视频| www.www免费av| 级片在线观看| 身体一侧抽搐| 少妇的逼好多水| 男人和女人高潮做爰伦理| av在线老鸭窝| 亚洲专区中文字幕在线| 99国产综合亚洲精品| 久久久久性生活片| 欧美激情久久久久久爽电影| 色噜噜av男人的天堂激情| 欧美成人性av电影在线观看| 国产v大片淫在线免费观看| 日韩国内少妇激情av| 亚洲美女黄片视频| 露出奶头的视频| 非洲黑人性xxxx精品又粗又长| 脱女人内裤的视频| 亚洲乱码一区二区免费版| 亚洲最大成人av| 日本撒尿小便嘘嘘汇集6| 91九色精品人成在线观看| 亚洲精品在线观看二区| 亚洲狠狠婷婷综合久久图片| 免费人成视频x8x8入口观看| 女生性感内裤真人,穿戴方法视频| 亚洲狠狠婷婷综合久久图片| 国内毛片毛片毛片毛片毛片| 欧美一区二区亚洲| 黄色视频,在线免费观看| 婷婷精品国产亚洲av在线| 国产精品不卡视频一区二区 | 成人特级黄色片久久久久久久| 婷婷亚洲欧美| 亚洲电影在线观看av| 亚洲欧美日韩无卡精品| 最近视频中文字幕2019在线8| 成人无遮挡网站| 一级作爱视频免费观看| 如何舔出高潮| 亚洲av成人精品一区久久| 久久久久久久久中文| 伦理电影大哥的女人| 亚洲精品在线观看二区| 国内揄拍国产精品人妻在线| 国产精品一区二区三区四区久久| 久久久久久国产a免费观看| 老鸭窝网址在线观看| 国产成人aa在线观看| 亚洲男人的天堂狠狠| 国产精品人妻久久久久久| 国产精品亚洲一级av第二区| 丁香欧美五月| 99riav亚洲国产免费| 久久久精品大字幕| 自拍偷自拍亚洲精品老妇| 观看美女的网站| 婷婷丁香在线五月| 亚洲国产高清在线一区二区三| 五月玫瑰六月丁香| 男人狂女人下面高潮的视频| 欧美不卡视频在线免费观看| 他把我摸到了高潮在线观看| 好男人在线观看高清免费视频| 中文字幕久久专区| 三级男女做爰猛烈吃奶摸视频| 国产白丝娇喘喷水9色精品| or卡值多少钱| av在线天堂中文字幕| 99久久无色码亚洲精品果冻| 18禁裸乳无遮挡免费网站照片| 两人在一起打扑克的视频| 中文字幕人妻熟人妻熟丝袜美| 午夜日韩欧美国产| 黄色配什么色好看| 亚洲美女黄片视频| АⅤ资源中文在线天堂| 久久久久久久久久黄片| 国产黄片美女视频| 国产精品一区二区性色av| 国内精品美女久久久久久| 美女大奶头视频| 女生性感内裤真人,穿戴方法视频| 国产精品爽爽va在线观看网站| а√天堂www在线а√下载| 国产精品久久久久久人妻精品电影| 亚洲欧美日韩高清专用| 又爽又黄无遮挡网站| 午夜免费男女啪啪视频观看 | 国产成人啪精品午夜网站| .国产精品久久| 精品午夜福利在线看| 国产精华一区二区三区| 国产大屁股一区二区在线视频| 69人妻影院| 国产在视频线在精品| 日韩欧美 国产精品| 琪琪午夜伦伦电影理论片6080| 亚洲成a人片在线一区二区| 看十八女毛片水多多多| 成人特级av手机在线观看| 极品教师在线视频| 12—13女人毛片做爰片一| 听说在线观看完整版免费高清| 精品久久久久久成人av| 中文字幕av成人在线电影| 精品免费久久久久久久清纯| 黄色丝袜av网址大全| 国产高清视频在线播放一区| 一个人免费在线观看电影| 99久久成人亚洲精品观看| 国产高潮美女av| 亚洲国产欧美人成| 国产av一区在线观看免费| 日本免费一区二区三区高清不卡| 日韩欧美免费精品| 啦啦啦观看免费观看视频高清| av在线观看视频网站免费| 欧美成人a在线观看| 国产探花在线观看一区二区| 久久国产乱子免费精品| 熟女人妻精品中文字幕| 在线a可以看的网站| 久久久久久国产a免费观看| 色av中文字幕| 国产精品综合久久久久久久免费| 国产伦人伦偷精品视频| av专区在线播放| 五月伊人婷婷丁香| 久9热在线精品视频| 欧美日韩国产亚洲二区| 亚洲av成人av| 欧美丝袜亚洲另类 | avwww免费| 99在线视频只有这里精品首页| 国产伦在线观看视频一区| 搡女人真爽免费视频火全软件 | www.999成人在线观看| 日韩高清综合在线| 日本熟妇午夜| 亚洲激情在线av| 欧美xxxx性猛交bbbb| 黄色一级大片看看| 精品久久国产蜜桃| 国产高清三级在线| 最新中文字幕久久久久| 一进一出抽搐动态| 欧美日韩黄片免| 搡老妇女老女人老熟妇| 成人欧美大片| 色吧在线观看| 如何舔出高潮| 亚洲真实伦在线观看| 在线观看午夜福利视频| 一个人看视频在线观看www免费| 日本精品一区二区三区蜜桃| 女人被狂操c到高潮| 成年女人看的毛片在线观看| 69人妻影院| 窝窝影院91人妻| 一区福利在线观看| 热99re8久久精品国产| 国产黄片美女视频| 偷拍熟女少妇极品色| 免费看日本二区| 国产又黄又爽又无遮挡在线| 国产av麻豆久久久久久久| а√天堂www在线а√下载| 国产三级黄色录像| 丰满乱子伦码专区| 两个人视频免费观看高清| 亚洲成a人片在线一区二区| 国产在线精品亚洲第一网站| 国产麻豆成人av免费视频| 啪啪无遮挡十八禁网站| 99久国产av精品| 欧美激情久久久久久爽电影| 热99在线观看视频| www.www免费av| 最近最新中文字幕大全电影3| 成年女人看的毛片在线观看| 舔av片在线| 一卡2卡三卡四卡精品乱码亚洲| 国产精品1区2区在线观看.| 久久精品久久久久久噜噜老黄 | av在线蜜桃| 午夜福利免费观看在线| 在线看三级毛片| 啦啦啦观看免费观看视频高清| 亚洲 国产 在线| 精品免费久久久久久久清纯| 午夜激情福利司机影院| 亚洲av中文字字幕乱码综合| 国内少妇人妻偷人精品xxx网站| 国产三级黄色录像| 午夜激情福利司机影院| 成人永久免费在线观看视频| 嫩草影视91久久| or卡值多少钱| 首页视频小说图片口味搜索| 成人高潮视频无遮挡免费网站| 少妇丰满av| 日日干狠狠操夜夜爽| 亚洲av中文字字幕乱码综合| 国产精品,欧美在线| 男人狂女人下面高潮的视频| 亚洲专区中文字幕在线| 俺也久久电影网| 免费看光身美女| 国产一区二区激情短视频| 尤物成人国产欧美一区二区三区| 成人三级黄色视频| 国产视频内射| 嫩草影视91久久| 国产亚洲精品久久久com| 99久久九九国产精品国产免费| 中文字幕熟女人妻在线| 夜夜躁狠狠躁天天躁| 久久精品久久久久久噜噜老黄 | 色综合婷婷激情| 欧美黑人欧美精品刺激| 中文字幕免费在线视频6| 波多野结衣高清无吗| 午夜福利视频1000在线观看| 国产精品亚洲美女久久久| 国产国拍精品亚洲av在线观看| 内地一区二区视频在线| 亚洲精品成人久久久久久| 欧美xxxx性猛交bbbb| 国产一级毛片七仙女欲春2| 最近最新免费中文字幕在线| 色综合婷婷激情| 十八禁网站免费在线| 中文字幕免费在线视频6| 2021天堂中文幕一二区在线观| 麻豆国产av国片精品| 亚洲精华国产精华精| 在线观看舔阴道视频| 变态另类丝袜制服| 啦啦啦韩国在线观看视频| 国产私拍福利视频在线观看| 欧美极品一区二区三区四区| 欧美又色又爽又黄视频| 日韩 亚洲 欧美在线| 日本成人三级电影网站| 精品免费久久久久久久清纯| 18+在线观看网站| 亚洲av成人av| 亚洲三级黄色毛片| 日本精品一区二区三区蜜桃| 又爽又黄a免费视频| 男女下面进入的视频免费午夜| 色视频www国产| 精品人妻偷拍中文字幕| 青草久久国产| 一本综合久久免费| 深夜精品福利| 丰满人妻熟妇乱又伦精品不卡| 欧美在线黄色| 18禁黄网站禁片午夜丰满| 亚洲无线观看免费| 亚洲国产高清在线一区二区三| 久久久精品欧美日韩精品| 高潮久久久久久久久久久不卡| 欧美午夜高清在线| 亚洲天堂国产精品一区在线| 男女做爰动态图高潮gif福利片| 一级黄片播放器| 悠悠久久av| 欧美乱妇无乱码| 国产亚洲精品久久久久久毛片| 丝袜美腿在线中文| 白带黄色成豆腐渣| 亚洲国产精品成人综合色| 18+在线观看网站| 久久久久久国产a免费观看| 人妻丰满熟妇av一区二区三区| www.色视频.com| 亚洲av电影不卡..在线观看| 村上凉子中文字幕在线| 男人狂女人下面高潮的视频| 亚洲乱码一区二区免费版| x7x7x7水蜜桃| 久久久久亚洲av毛片大全| 国产成人aa在线观看| 国产高潮美女av| 久久久国产成人免费| 全区人妻精品视频| 欧美乱妇无乱码| 日本黄色视频三级网站网址| 俺也久久电影网| 小说图片视频综合网站| 午夜免费激情av| 国产亚洲精品久久久久久毛片| 免费人成视频x8x8入口观看| 伦理电影大哥的女人| 成人午夜高清在线视频| 亚洲最大成人中文| 3wmmmm亚洲av在线观看| 欧美激情在线99| 麻豆国产97在线/欧美| 国产日本99.免费观看| 欧美精品啪啪一区二区三区| 精品一区二区免费观看| 免费看美女性在线毛片视频| 日韩欧美在线二视频| 国产精品综合久久久久久久免费| 少妇裸体淫交视频免费看高清| 国产在线精品亚洲第一网站| 久久久久亚洲av毛片大全| 丝袜美腿在线中文| 国产精品电影一区二区三区| 桃色一区二区三区在线观看| 少妇的逼水好多| 国内少妇人妻偷人精品xxx网站| 一个人免费在线观看的高清视频| 内地一区二区视频在线| 精品人妻偷拍中文字幕| 岛国在线免费视频观看| 国产在视频线在精品| 国产男靠女视频免费网站| 18美女黄网站色大片免费观看| 麻豆一二三区av精品| 精品人妻视频免费看| 性色avwww在线观看| 99久久精品一区二区三区| 久久久成人免费电影| 51国产日韩欧美| 又紧又爽又黄一区二区| 深爱激情五月婷婷| 天天一区二区日本电影三级| 两个人视频免费观看高清| 黄色一级大片看看| 亚洲精华国产精华精| 国产老妇女一区| 免费观看的影片在线观看| 精品一区二区三区视频在线观看免费| 国产人妻一区二区三区在| 日韩免费av在线播放| 国产美女午夜福利| av中文乱码字幕在线| 亚洲一区高清亚洲精品| eeuss影院久久| 国产久久久一区二区三区| avwww免费| 日本三级黄在线观看| 国产真实伦视频高清在线观看 | 国产麻豆成人av免费视频| 亚洲精品在线观看二区| 婷婷丁香在线五月| 一进一出好大好爽视频| 国产综合懂色| 久久久久久国产a免费观看| 在线播放无遮挡| 国产伦人伦偷精品视频| 欧美最新免费一区二区三区 | 激情在线观看视频在线高清| 久久天躁狠狠躁夜夜2o2o| 国产视频一区二区在线看| 欧美成狂野欧美在线观看| 麻豆久久精品国产亚洲av| 51国产日韩欧美| 午夜影院日韩av| 如何舔出高潮| 我的女老师完整版在线观看| 99热精品在线国产| 偷拍熟女少妇极品色| 757午夜福利合集在线观看| 欧美黑人巨大hd| 日韩中文字幕欧美一区二区| 午夜精品久久久久久毛片777| 亚洲av成人不卡在线观看播放网| 女人被狂操c到高潮| 午夜a级毛片| 最近最新中文字幕大全电影3| 亚洲熟妇熟女久久| 国产精品不卡视频一区二区 | 成人特级黄色片久久久久久久| 久久久国产成人免费| 神马国产精品三级电影在线观看| 中出人妻视频一区二区| 亚洲av.av天堂| 天堂网av新在线| 露出奶头的视频| 伊人久久精品亚洲午夜| 俺也久久电影网| 又紧又爽又黄一区二区| av女优亚洲男人天堂| av中文乱码字幕在线| 亚洲av电影不卡..在线观看| 久久中文看片网| 国产精品野战在线观看| 一个人免费在线观看电影| 黄色日韩在线| 成人性生交大片免费视频hd| 人人妻人人澡欧美一区二区| 亚洲美女搞黄在线观看 | av欧美777| 久久国产精品人妻蜜桃| 久久久久久国产a免费观看| 欧美+日韩+精品| 久久天躁狠狠躁夜夜2o2o| 三级国产精品欧美在线观看| 国产精品一区二区三区四区免费观看 | 欧美色视频一区免费| av天堂中文字幕网| 精品欧美国产一区二区三| 亚洲精品在线观看二区| 亚洲久久久久久中文字幕| 97超视频在线观看视频| 人人妻人人看人人澡| 狠狠狠狠99中文字幕| 熟女电影av网| 免费在线观看成人毛片| 欧美bdsm另类| 国产色婷婷99| 制服丝袜大香蕉在线| 亚洲经典国产精华液单 | 亚洲乱码一区二区免费版| 日本三级黄在线观看| 国产国拍精品亚洲av在线观看| 国产午夜精品论理片| 精品久久久久久久末码| 欧美一区二区精品小视频在线| aaaaa片日本免费| 无遮挡黄片免费观看| 欧美黄色淫秽网站| 又粗又爽又猛毛片免费看| av女优亚洲男人天堂| 91久久精品电影网| 午夜激情欧美在线| 一a级毛片在线观看| 亚洲在线观看片| 国产欧美日韩精品亚洲av| 男女那种视频在线观看| 亚洲av成人av| 伊人久久精品亚洲午夜| 久久久久九九精品影院| 不卡一级毛片| 女人被狂操c到高潮| 国产私拍福利视频在线观看| netflix在线观看网站| 一个人观看的视频www高清免费观看| 国产亚洲精品久久久久久毛片| 午夜激情福利司机影院| 国产伦精品一区二区三区视频9| 国产精品电影一区二区三区| 欧美成人性av电影在线观看| 午夜福利在线观看免费完整高清在 | 国产精品1区2区在线观看.| 午夜精品久久久久久毛片777| 亚洲精品日韩av片在线观看| 欧美一级a爱片免费观看看| 亚洲av一区综合| 亚洲va日本ⅴa欧美va伊人久久| 一级毛片久久久久久久久女| 99国产精品一区二区三区| 99国产精品一区二区蜜桃av| 在线观看66精品国产| 亚洲欧美清纯卡通| 亚洲一区二区三区色噜噜| 别揉我奶头~嗯~啊~动态视频| 中文字幕av成人在线电影| 午夜免费成人在线视频| 麻豆一二三区av精品| 国产精品女同一区二区软件 | 国内久久婷婷六月综合欲色啪| 高清日韩中文字幕在线| 一夜夜www| 成年女人看的毛片在线观看| 美女cb高潮喷水在线观看| 亚洲国产高清在线一区二区三| 久久久国产成人免费| av在线观看视频网站免费| 亚洲av电影不卡..在线观看| 免费看a级黄色片| 极品教师在线视频| eeuss影院久久| 亚洲精品乱码久久久v下载方式| 久久99热这里只有精品18| 午夜福利成人在线免费观看| 国产午夜精品久久久久久一区二区三区 | 亚洲第一区二区三区不卡| 中文字幕熟女人妻在线| 真人做人爱边吃奶动态| 精品午夜福利在线看| 成人欧美大片| 露出奶头的视频| 国产免费男女视频| 免费人成视频x8x8入口观看| av福利片在线观看| 国产亚洲精品久久久com| 老熟妇仑乱视频hdxx| 久99久视频精品免费| 午夜免费成人在线视频| 99久久精品国产亚洲精品| 成人三级黄色视频| 亚洲美女搞黄在线观看 | 日韩大尺度精品在线看网址| 亚洲国产欧美人成| 在线观看舔阴道视频| 天美传媒精品一区二区| 在线观看66精品国产| 欧美乱妇无乱码| 老熟妇仑乱视频hdxx| 夜夜夜夜夜久久久久| 少妇被粗大猛烈的视频| 免费人成在线观看视频色| 老鸭窝网址在线观看| 久久久久亚洲av毛片大全| 真实男女啪啪啪动态图| 午夜亚洲福利在线播放| 成人毛片a级毛片在线播放| 欧美极品一区二区三区四区| 色哟哟·www| 一级黄色大片毛片| 在线看三级毛片| 欧美黑人欧美精品刺激| 成人国产综合亚洲| bbb黄色大片| 宅男免费午夜| 国产精品乱码一区二三区的特点| 精品人妻熟女av久视频| 香蕉av资源在线| 一个人看视频在线观看www免费| 别揉我奶头~嗯~啊~动态视频| 婷婷六月久久综合丁香| 免费人成在线观看视频色| 嫩草影院入口| 午夜影院日韩av| 中文字幕久久专区| 国产精品久久久久久人妻精品电影| 成人国产一区最新在线观看| 色视频www国产| 99在线人妻在线中文字幕| 国产av一区在线观看免费| 深夜a级毛片| 成人一区二区视频在线观看| 色精品久久人妻99蜜桃| 自拍偷自拍亚洲精品老妇| 免费观看精品视频网站| 天天躁日日操中文字幕| 日本五十路高清| 男女那种视频在线观看| 69av精品久久久久久| 午夜福利在线在线| 欧美一级a爱片免费观看看| 欧美成人a在线观看| 久久国产乱子伦精品免费另类| 国产欧美日韩精品一区二区| 99riav亚洲国产免费| eeuss影院久久| 国产精品不卡视频一区二区 | 国产午夜福利久久久久久|