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

    太原市大氣PM2.5季節(jié)傳輸路徑和潛在源分析

    2019-08-28 03:08:12任浦慧解靜芳姜洪進王淑楠劉瑞卿
    中國環(huán)境科學(xué) 2019年8期
    關(guān)鍵詞:太原市氣流顆粒物

    任浦慧,解靜芳,姜洪進,王淑楠,劉瑞卿

    太原市大氣PM2.5季節(jié)傳輸路徑和潛在源分析

    任浦慧,解靜芳*,姜洪進,王淑楠,劉瑞卿

    (山西大學(xué)環(huán)境與資源學(xué)院,山西 太原 030006)

    為了研究太原市大氣PM2.5不同季節(jié)的傳輸路徑和污染源區(qū),利用HYSPLIT后向軌跡模型和NCEP 的GDAS 全球氣象要素數(shù)據(jù),對2017~2018年不同季節(jié)太原市逐日48h氣流后向軌跡進行聚類分析,同時結(jié)合小時污染物質(zhì)量濃度數(shù)據(jù),分析不同季節(jié)太原市PM2.5的潛在源貢獻因子(WPSCF)和濃度權(quán)重軌跡(WCWT).結(jié)果表明,太原市PM2.5的質(zhì)量濃度在季節(jié)上呈現(xiàn)冬季(77.56μg/m3)>秋季(69.89μg/m3)>春季(63.78μg/m3)>夏季(45.51μg/m3)的變化趨勢.PM2.5與SO2、NO2和CO之間存在明顯的同源性和二次轉(zhuǎn)化過程.春、秋和冬季大氣傳輸路徑主要以西和西北方向近距離、慢移速的軌跡為主,夏季以南和東方向軌跡為主.PM2.5潛在源區(qū)季節(jié)變化明顯:夏季主要受太原本地和晉中地區(qū)的影響;春、秋和冬季主要受陜西中北部、呂梁、臨汾和晉中等地的影響.

    PM2.5;后向軌跡;聚類分析;潛在源貢獻因子法;濃度權(quán)重軌跡法

    目前,太原市關(guān)于 PM2.5的研究主要集中于冬季采暖期間 PM2.5和氣態(tài)污染物質(zhì)量濃度的變化,研究時段相對較短[1-3];研究內(nèi)容主要集中在PM2.5的組分特征、來源解析和健康風(fēng)險評估等方面[4-6],而關(guān)于傳輸路徑和潛在來源區(qū)域的研究尚不多見.后向軌跡模型(HYSPLIT)、聚類分析、潛在源貢獻因子方法(PSCF)和濃度權(quán)重軌跡方法(CWT)均是通過氣流軌跡分析、識別大氣污染物潛在來源,進而研究污染物輸送擴散的有效手段.主要針對區(qū)域污染物分布、氣流輸送途徑、不同方向氣流軌跡對研究區(qū)域的傳輸以及影響空氣質(zhì)量濃度的潛在源區(qū)等進行研究[7-10].此外,有研究者發(fā)現(xiàn)氣流輸送途徑、軌跡特征及其污染物的分布存在明顯的季節(jié)變化[11-13].

    太原市的地形特征為三面環(huán)山,北高南低的簸箕型,大氣PM2.5不僅受本地排放源影響,也受周圍地域排放源的傳輸影響,而不同季節(jié)主導(dǎo)風(fēng)向不同、傳輸路徑各異,不同地域污染源強度變化不同,因此研究PM2.5的時空變化特征、傳輸途徑和潛在來源對解決太原市PM2.5污染問題有著重要的意義.由此,本研究在分析2017年6月1日~2018年5月31日太原市PM2.5、PM10、SO2、NO2、CO和O3濃度時間變化特征的基礎(chǔ)上,使用Spearman相關(guān)系數(shù)分析PM2.5與其他污染物之間的相關(guān)性.并使用后向軌跡聚類分析、PSCF和CWT方法,研究了太原市不同季節(jié)PM2.5的傳輸途徑和潛在來源.

    1 數(shù)據(jù)來源與方法

    1.1 數(shù)據(jù)來源

    太原市大氣PM2.5、PM10、SO2、NO2、CO和O3的監(jiān)測數(shù)據(jù),來自太原市空氣質(zhì)量實時數(shù)據(jù)發(fā)布平臺(http://183.203.223.83:85/aqi/);氣象數(shù)據(jù)(溫度、濕度、風(fēng)速)源自中國天氣網(wǎng)(http://www.weather.com.cn/);后向軌跡聚類、PSCF和CWT分析所需的氣象場數(shù)據(jù)來自美國國家環(huán)境預(yù)報中心(NCEP)的全球資料同化系統(tǒng)(GDAS) (https://ready.arl.noaa. gov/HYSPLIT. php).研究時段為2017年6月1日~ 2018年5月31日.

    1.2 方法

    HYSPLIT方法是將太原市市中心(37.54°N112.33°E)作為后向軌跡的起始點,模擬起始高度設(shè)置為500m,該高度能夠準確反映邊界層平均流場特征[14].以每日12:00前推48h進行軌跡計算.同時結(jié)合GIS 技術(shù)開發(fā)的TrajStat軟件,將每個季節(jié)的氣流軌跡進行聚類分析[15],再結(jié)合對應(yīng)的污染物質(zhì)量濃度特征進行統(tǒng)計分析.

    PSCF分析法[16]主要根據(jù)經(jīng)過某一網(wǎng)格的污染軌跡數(shù)()占經(jīng)過此網(wǎng)格的所有軌跡數(shù)()的比值來確定PSCF值的大小,PSCF的值越大代表經(jīng)過此網(wǎng)格的氣團對接受點的空氣質(zhì)量影響的概率越大.具體公式如下:

    由于PSCF是一種條件概率,當(dāng)分母n較小時,計算的PSCF值會有很大不確定性.為了減少不確定性,故引入權(quán)重因子,將PSCF值乘以一個權(quán)重函數(shù)(W)[17],即WPSCF=W×PSCF.

    W定義如下:

    PSCF值是一種條件概率,在一定程度上可反映網(wǎng)格對接受點污染程度的貢獻大小,但無法確定研究區(qū)域的污染程度.因此,本文進一步引用CWT分析法[18-19],結(jié)合()計算每個網(wǎng)格中軌跡的污染權(quán)重指數(shù)來反映不同軌跡的污染程度.計算公式如下:

    式中: CWT是網(wǎng)格的平均權(quán)重污染物濃度;C是軌跡經(jīng)過網(wǎng)格對應(yīng)的污染物濃度; τ是軌跡在網(wǎng)格的停留時間.同樣,為減少m值較小時所引起的不確定性,將PSCF分析法中的權(quán)重函數(shù)W適用于CWT分析法,即:

    WCWT =W×CWT (4)

    2 結(jié)果與討論

    2.1 PM2.5和PM10質(zhì)量濃度季節(jié)變化特征

    研究期間,太原市PM2.5和PM10的年均質(zhì)量濃度分別為(64±37) μg/m3和(140±64) μg/m3,二者分別為《環(huán)境空氣質(zhì)量標(biāo)準》(GB 3095-2012)[20]二級標(biāo)準(35和70μg/m3)的1.83和2.00倍.太原市PM2.5較2016年相比有所下降(約15%),空氣質(zhì)量有明顯的改善,但仍高于同期的北京和天津地區(qū)(https://www. aqistudy.cn).由圖1可知,研究期間31%的天數(shù)超過了PM2.5的二級日均標(biāo)準限值75μg/m3,42%的天數(shù)超過了PM10二級標(biāo)準限值150μg/m3.且PM2.5與PM10有顯著的相關(guān)性,春、夏、秋和冬的相關(guān)系數(shù)分別為0.766**、0.805**、0.894**和0.916**, ρ(PM2.5)/ ρ(PM10)可以反映可吸入顆粒物中細顆粒物的含量,研究期間(PM2.5)/(PM10)均值為0.45±0.13,且細顆粒物PM2.5的質(zhì)量濃度與PM10日變化情況一致.

    由圖2可知,太原市PM2.5的質(zhì)量濃度在季節(jié)上呈現(xiàn)冬季(77.56μg/m3)>秋季(69.89μg/m3)>春季(63.78μg/m3)>夏季(45.51μg/m3的變化趨勢.其中冬季PM2.5的質(zhì)量濃度是夏季的1.70倍,類似的季節(jié)變化可能是由于各季節(jié)主要的污染源、排放量以及氣象條件存在差異.夏季,太原市大氣垂直擴散條件良好,抑制顆粒物的積累,還受東南季風(fēng)的影響,帶來大量的降水,加速細顆粒物的濕沉降,從而使PM2.5質(zhì)量濃度較低;冬季,煤炭燃燒等人類活動使PM2.5的排放量增加,且邊界層高度低,低濕、無風(fēng)的靜穩(wěn)天氣不利于污染物的擴散[21].此外,夏季PM2.5的質(zhì)量濃度相對集中,而秋、冬季的差異較大;根據(jù)中位數(shù)及均值,可以看出夏、冬季的PM2.5質(zhì)量濃度基本呈正態(tài)分布.

    然而,與PM2.5不同的是:PM10在四季的平均質(zhì)量濃度從高到低依次為春季(162.65μg/m3)>秋季(155.99μg/m3)>冬季(142.43μg/m3)>夏季(99.22μg/m3).除了工業(yè)來源和不利的天氣條件對冬季空氣污染稀釋和擴散的綜合影響外,春季還會受到來自北方沙塵的嚴重影響.生物質(zhì)燃燒(例如秸稈露天焚燒)可能是導(dǎo)致秋季PM10增多的原因.此外,除夏季外,秋、冬和春季的PM10均值與中位數(shù)基本一致,而冬季PM10的均值明顯高于中位數(shù),這可能是由于冬季發(fā)生幾次以PM10為主要污染物的重污染過程.

    圖1 2017年6月~2018年5月太原市PM2.5和PM10的日均質(zhì)量濃度

    圖2 太原市不同季節(jié)PM2.5和PM10的質(zhì)量濃度

    2.2 不同季節(jié)PM2.5與氣態(tài)前體物及氣象因素之間的關(guān)系

    由表1可知,PM2.5與SO2、NO2和CO在4個季節(jié)均呈不同程度的顯著正相關(guān);而PM2.5與O3只在夏季和冬季分別呈現(xiàn)相反的正、負相關(guān)關(guān)系.大氣中PM2.5不僅來自于污染源的直接排放,也來自于SO2和NO等氣態(tài)前體物在大氣中通過均相或非均相(在顆粒物表面)反應(yīng)生成硫酸鹽、硝酸鹽和有機氣溶膠等二次污染物[22].上述太原市大氣PM2.5與SO2、NO2和CO之間的良好正相關(guān)性,表明其相互間存在明顯的二次轉(zhuǎn)化過程,大氣中的SO2和NO2通過均相或非均相反應(yīng)生成硫酸鹽和硝酸鹽,對PM2.5質(zhì)量濃度變化具有重要影響.夏季高溫、低濕和強的太陽輻射易產(chǎn)生大量的·OH,繼而與VOCs反應(yīng)形成O3.O3和·OH均可加速顆粒物表面的均相或非均相化學(xué)反應(yīng),生成二次顆粒物;而冬季PM2.5形成的氣溶膠消光作用對O3生成產(chǎn)生抑制作用,所以,PM2.5與O3在夏季和冬季分別呈現(xiàn)相反的正、負相關(guān)性.

    表1 不同季節(jié)PM2.5與氣態(tài)前體物的Spearman 秩相關(guān)系數(shù)

    注:*表示<0.05, **表示<0.01.

    此外,由表1可知,PM2.5與SO2、NO2和CO具有較強的同源性,SO2主要來源為化石燃料的燃燒,NO2主要來自機動車排放和燃煤,CO主要來源于冶金工業(yè)、內(nèi)燃機排氣和化石燃料的不完全燃燒.因此,煤炭燃燒和汽車尾氣的排放對太原市大氣PM2.5具有較高的貢獻.

    由表2可知,春季PM2.5與溫度、濕度和風(fēng)速均無明顯的相關(guān)關(guān)系;夏季PM2.5只與溫度呈顯著正相關(guān); 秋季PM2.5與溫度無明顯的相關(guān)關(guān)系,而分別與濕度和風(fēng)速呈顯著正負相關(guān)性;冬季PM2.5與溫度和濕度均呈顯著正相關(guān),與風(fēng)速呈顯著負相關(guān).

    夏季高溫,太陽輻射增強,可促進氣態(tài)前體物氣相和顆粒物表面氧化反應(yīng)的進行,加速二次污染的形成,導(dǎo)致夏季溫度升高PM2.5質(zhì)量濃度增大;秋季較高的濕度能促使大氣硫氧化物和氮氧化物被氧化成二次硫酸鹽和硝酸鹽[23],加速顆粒物的吸濕增長,同時較低的風(fēng)速減弱了對顆粒物的擴散作用,易導(dǎo)致PM2.5的積累;冬季取暖化石燃料燃燒細顆粒物排放增大,低風(fēng)速、高濕和逆溫天氣都是導(dǎo)致顆粒物污染嚴重的重要因素.

    表2 不同季節(jié)PM2.5與氣象因素的Spearman 秩相關(guān)系數(shù)

    注:*表示<0.05, **表示<0.01.

    2.3 后向軌跡模型對PM2.5、PM10、SO2和NO2濃度的季節(jié)變化分析

    根據(jù)不同季節(jié)PM2.5與其他污染物之間相關(guān)性分析結(jié)果,選擇相互關(guān)系較為密切的PM10、SO2、NO2進行后向軌跡聚類分析,結(jié)果見圖3.

    由圖3可知,研究期間,不同季節(jié)氣流軌跡最終聚類的總類數(shù)不同,春、夏和冬均為4類,只有秋季為5類.由不同季節(jié)氣流軌跡來看,太原市不同季節(jié)氣流軌跡與季風(fēng)氣候關(guān)系密切.春、秋和冬季受冷空氣南下的影響,盛行西北風(fēng),西北方向軌跡較長,但占比最高的氣流軌跡主要集中于太原市西和西北方向,且以近距離、移速慢的軌跡為主;而夏季氣流軌跡較為分散,以南和東方向軌跡為主.

    由圖4可知,除夏季以外,其他3個季節(jié)基本遵循氣流軌跡占比越大攜帶的污染物質(zhì)量濃度越高. 3個季節(jié)中來自西和西北方向近距離傳輸?shù)臍饬鬈壽E占比最高,攜帶污染物(PM2.5、PM10、SO2和NO2)質(zhì)量濃度亦最高.但亦存在氣流軌跡占比最小,攜帶的污染物質(zhì)量濃度較大的情況.如春季來自西北方向的遠距離傳輸占比最小,卻攜帶了較多的污染物,且可以看出以PM10貢獻為主,可能是傳輸過程中攜帶了部分揚沙所致;秋季氣流軌跡4和5雖然占比相同,但所攜帶的污染物質(zhì)量濃度存在明顯差異,主要以SO2和NO2的差異最為明顯,說明氣流軌跡占比和其攜帶污染物質(zhì)量濃度大小之間的關(guān)系與軌跡的長度有關(guān).相同占比情況下,軌跡越短所攜帶的污染物質(zhì)量濃度越高.冬季氣流軌跡2、3和4占比大小差異不大,但氣流軌跡3攜帶的污染物質(zhì)量濃度最小,可能是由于太原市的地形特征是三面環(huán)山,北高南低的簸箕型,來自北方的氣流能吹散集聚的污染物,因此攜帶的污染物質(zhì)量濃度較低.夏季不同氣流軌跡長短不同,占比各異,但不同軌跡攜帶的污染物質(zhì)量濃度差異不大,只有來自東邊的氣流軌跡攜帶污染物(PM2.5、PM10、SO2和NO2)質(zhì)量濃度最小.而來自西北方向的氣流軌跡攜帶了較高的SO2和NO2.除此之外,唯有夏季有來自南面山西長治和晉中等地的氣流軌跡1,且攜帶的PM2.5和PM10質(zhì)量濃度相對較高.

    綜上所述,春、秋和冬季,污染物質(zhì)量濃度較高的輸送路徑均來自西和西北方向的近距離輸送,氣流受太行山的阻擋,在太原盆地聚集,顆粒物質(zhì)量濃度增大;夏季主要為南和東方向氣流軌跡,其中來自南邊的軌跡污染物質(zhì)量濃度較高,受地形影響不易擴散;而來自偏東的軌跡經(jīng)過地區(qū)污染源少,攜帶污染物質(zhì)量濃度較低,有利于本地污染物的稀釋和擴散.

    2.4 PSCF和CWT對不同季節(jié)PM2.5潛在源分析

    為了進一步厘清PM2.5潛在來源,采用PSCF進行計算分析,將HYSPLIT 計算的氣流軌跡所覆蓋的區(qū)域網(wǎng)格化,網(wǎng)格大小為0.25°′0.25°, PM2.5標(biāo)準值設(shè)為《環(huán)境空氣質(zhì)量標(biāo)準》(GB 3095-2012)[20]中的二級日均標(biāo)準限值75μg/m3, WPSCF計算結(jié)果見圖5.

    圖5 太原市四季PM2.5的WPSCF 分布特征

    由圖5可知,不同季節(jié)潛在源區(qū)分布不同.夏季PM2.5質(zhì)量濃度最低,其WPSCF值大多小于0.1,其相對較大值集中在太原和呂梁地區(qū),太原本地源貢獻較高;春季太原市及周邊的WPSCF值介于0.1~0.4,大于0.7的區(qū)域主要集中于經(jīng)陜西東部、山西呂梁、臨汾等地向太原市傳輸路徑的汾渭平原一帶;秋季太原市及周邊的WPSCF值介于0.2~0.5,WPSCF大于0.6的區(qū)域主要在山西呂梁、臨汾等地和陜西的東北部;冬季W(wǎng)PSCF的整體值較其他3個季節(jié)是最高的,且潛在污染源區(qū)集中,WPSCF值大于0.8的區(qū)域主要分布在太原市西部和山西呂梁地區(qū).通過PSCF結(jié)果可以看到,潛在源區(qū)存在明顯的季節(jié)變化,但也存在長時間影響的潛在源區(qū),如陜西的東北部在秋冬季W(wǎng)PSCF值均較高;此外,山西呂梁地區(qū)在秋冬季的WPSCF值也較高.

    為了進一步驗證WPSCF的分析結(jié)果,采用CWT方法,同樣將計算的氣流軌跡所覆蓋的區(qū)域網(wǎng)格化,網(wǎng)格大小亦為0.25°×0.25°, WCWT分析結(jié)果見圖6.

    圖6 太原市四季PM2.5的WCWT分布特征

    由圖6可見,與WPSCF結(jié)果類似,夏季W(wǎng)CWT值最小,潛在源區(qū)的WCWT數(shù)值基本小于75μg/m3,較高WCWT值的貢獻源區(qū)主要集中在太原市、晉中等地,這與夏季盛行東南風(fēng)有關(guān).春季主要是遠距離運輸,WCWT值在75~100mg/m3集中在由陜西經(jīng)運城、呂梁、晉中向太原運輸?shù)穆窂?同樣秋季W(wǎng)CWT值大于100μg/m3集中在山西的太原市南部、呂梁地區(qū)和陜西北部.冬季W(wǎng)CWT最大,貢獻源區(qū)也有所變化,向西和向南延伸,除山西中部、陜西中北部外,山西南部也有貢獻,這與各地取暖增加污染排放有關(guān),太原市西、北、東三面環(huán)山,整個地形北高南低,且冬季多均壓場的天氣形勢,本地污染物不易擴散,加之地面風(fēng)場輻合,易導(dǎo)致周邊污染匯聚.總之,春季和秋季主要受陜西和呂梁等地的影響,夏季以太原市本地源的貢獻為主,冬季主要來自太原市本地及周邊城市如晉中、呂梁和長治等地的影響.

    3 結(jié)論

    3.1 研究時間段內(nèi),太原市PM2.5質(zhì)量濃度在季節(jié)變化上呈現(xiàn)冬季(77.56μg/m3)>秋季(69.89μg/m3)>春季(63.78μg/m3)>夏季(45.51μg/m3)的趨勢.

    3.2 不同季節(jié)PM2.5與PM10、SO2、NO2和CO均呈顯著正相關(guān)關(guān)系,表明其存在明顯的同源性;而PM2.5與O3在夏季和冬季分別呈現(xiàn)相反的正、負相關(guān)關(guān)系.分析氣象因素,夏季高溫、秋冬季低風(fēng)速和較高濕度是導(dǎo)致細顆粒物污染嚴重的重要氣象因素.

    3.3 太原市大氣污染輸送路徑以近距離、移速慢的軌跡為主.其中春、秋和冬季以西和西北方向軌跡為主;夏季軌跡較分散,以南和東方向軌跡為主,各軌跡攜帶高質(zhì)量濃度污染物的軌跡與該軌跡占比基本一致.

    3.4 對太原市PM2.5潛在源貢獻因子和濃度權(quán)重軌跡分析表明,太原市PM2.5主要污染源區(qū)范圍較小,夏季主要受太原本地和呂梁地區(qū)的短距離輸送的影響,春、秋和冬季主要集中于陜西中北部、呂梁、臨汾和晉中等地區(qū).

    [1] 賈小花,解靜芳,馬 翔,等.太原市冬季PM2.5水溶性組分污染特征分析 [J]. 中國環(huán)境科學(xué), 2013,33(4):599-604. Jia X H, Xie J F, Ma X, et al. Analysis of water-soluble constituents in winter of PM2.5in Taiyuan City [J]. China Environmental Science (in Chinese), 2013,33(4):599-604.

    [2] 李義宇,楊鴻儒,王 楠,等.太原市一次重污染天氣過程的成因分析 [J]. 氣象與環(huán)境學(xué)報, 2018,34(2):11-18. Li Y Y, Yang H R, Wang N, et al. Formation analysis of a severe air pollution event in Taiyuan [J]. Journal of Meteorology and Environment(in Chinese), 2018,34(2):11-18.

    [3] 耿 紅,宣瑩瑩,蔡夏童.太原市2014年春節(jié)期間常規(guī)大氣污染物濃度變化及聚類分析[J]. 環(huán)境科學(xué)學(xué)報, 2015,35(4):965-974. Geng H, Xuan Y Y, Cai X T. Mass concentration variation and cluster analysis of urban air pollutants in Taiyuan, Shanxi Province during Chinese New Year of 2014 [J]. Acta Scientiae Circumstantiae, 2015, 35(4):965-974.

    [4] Zhang M, Xie J F, Wang Z T, et al. Determination and source identification of priority polycyclic aromatic hydrocarbons in PM2.5in Taiyuan, China [J]. Atmospheric Research, 2016,178:401-414.

    [5] He Q S, Yan Y L, Guo L L, et al. Characterization and source analysis of water-soluble inorganic ionic species in PM2.5in Taiyuan city, China [J]. Atmospheric Research, 2017,184:48-55.

    [6] 李麗娟,溫彥平,彭 林,等.太原市采暖季PM2.5中元素特征及重金屬健康風(fēng)險評價 [J]. 環(huán)境科學(xué), 2014,35(12):0250-3301.Li L J, Wen Y P, Peng L, et al. Characteristic of elements in PM2.5and health risk assessment of heavy metals during heating season in Taiyuan [J]. Environmental Science (in Chinese), 2014,35(12):0250- 3301.

    [7] 鄧發(fā)榮,康 娜,Kanike R K,等.長江三角洲地區(qū)大氣污染過程分析[J]. 中國環(huán)境科學(xué), 2018,38(2):401-411. Deng F R, Kang N, Kanike R K, et al. Analysis of air pollution episodes over different cites in the Yangtze River Delta [J]. China Environmental Science(in Chinese), 2018,38(2):401-411.

    [8] 翟 華,朱 彬,趙雪婷,等.長江三角洲初冬一次重污染天氣成因分析 [J]. 中國環(huán)境科學(xué), 2018,38(11):4001-4009.Zhai H, Zhu B, Zhao X T, et al. Analysis of a heavy air pollution event in early winter in the Yangtze River Delta [J]. China Environmental Science(in Chinese), 2018,38(11):4001-4009.

    [9] Zhao Y Z, Man S W, Kwon H L.Estimation of potential source regions of PM2.5in Beijing using backward trajectories [J]. Atmospheric Pollution Research, 2015,6(5):173-177.

    [10] Liu B, Wu J, Zhang J, et al. Characterization and source apportionment of PM2.5based on error estimation from EPA PMF 5.0 model at a medium city in China [J]. Environmental Pollution, 2017,222:10-22.

    [11] 任傳斌,吳立新,張媛媛,等.北京城區(qū)PM2.5輸送途徑與潛在源區(qū)貢獻的四季差異分析 [J]. 中國環(huán)境科學(xué), 2016,36(9):2591-2598. Ren C B, Wu L X, Zhang Y Y, et al. Analyze to the seasonal differences of transport pathways and potential source-zones of Beijing Urban PM2.5[J]. China Environmental Science(in Chinese), 2016,36(9):2591-2598.

    [12] Sun K, Liu X G, Gu J W, et al. Chemical characterization of size resolved aerosols in four seasons and haze days in the megacity Beijingof China [J]. Journal of environmental Sciences, 2015,32: 155-167.

    [13] 史旭榮,溫 杰,田瑛澤,等.天津不同氣團來向PM2.5中組分和污染源貢獻的季節(jié)變化 [J]. 中國環(huán)境科學(xué), 2018,38(7):2406-2144. Shi X R, Wen J, Tian Y Z, et al. Seasonal patterns of PM2.5sources and chemical composition from different air mass directions in Tianjin [J]. China Environmental Science (in Chinese), 2018,38(7):2406-2144.

    [14] 趙 恒,王體健,江 飛,等.利用后向軌跡模式研究TRACE-P期間香港大氣污染物的來源 [J]. 熱帶氣象學(xué)報, 2009,25(2):181-186. Zhao H, Wang T J, Jiang F et al. Investigation into the source of air pollutants to Hong Kong by using backward trajectory method during the TRACE-P campaign [J]. Journal of Tropical Meteorology(in Chinese), 2009,25(2):181-186.

    [15] Wang Y Q, Zhang X Y, Draxler R R. TrajStat: GIS-based software that uses various trajectory statistical analysis methods to identify potential sources from long-term air pollution measurement data [J]. Environmental Modelling & Software, 2009,24(8):938-939.

    [16] Begum B A, Kim E, Jeong C H, et al. Evaluation of the potential source contribution function usng the 2002 Quebee forest fire episode [J]. Atmospheric Environment, 2005,39(20):3719-3724.

    [17] Polissar A V, Hopke P K, Paatero P, et al. The aerosol at Barrow, Alaska: long-term trends and source locations [J]. Atmospheric Environment, 1999,33(16):2441-2458.

    [18] Seibert P, Kromp-Kolb H, Baltensperger U, et al. Trajectory analysis of high-alpine air pollution data [M]. Air Pollution Modeling and Its Application X. Springer US, 1994:253-69.

    [19] Hsu Y K, Holsen T M, Hopke P K. Comparison of hybrid receptor models to locate PCB sources in Chicago [J]. Atmospheric Environment, 2003,37(4):545-562.

    [20] GB3095-2012 環(huán)境空氣質(zhì)量標(biāo)準 [S].

    [21] Tao J, Gao J, Zhang L, et al. PM2.5pollution in a megacity of southwest China: Source apportionment and implication [J]. Atmospheric Chemistry and Physics., 2014,14,8679–8699.

    [22] 丁萌萌,周健楠,劉保獻,等.2015年北京城區(qū)大氣PM2.5中NH4+、NO3-、SO42-及前體氣體的污染特征[J]. 環(huán)境科學(xué), 2017,38(4): 1307-1316. Ding M M, Zhou J N, Liu B X,et al. Pollution characteristics of NH4+, NO3-, SO42-in PM2.5and their precursor gases during 2015 in an urban area of Beijing [J]. Environmental Science(in Chinese), 2017, 38(4):1307-1316.

    [23] Sun Y L, Zhuang G S, Tang A, et al. Chemical characteristics of PM2.5and PM10in haze-fog episodes in Beijing [J]. Environmental Science & Technology, 2006,40(10):3148-3155.

    Transport pathways and potential sources of PM2.5in different seasons in Taiyuan city.

    REN Pu-hui, XIE Jing-fang*, JIANG Hong-jin, WANG Shu-nan, LIU Rui-qing

    (College of Environmental & Resource Sciences, Shanxi University, Taiyuan 030006, China)., 2019,39(8):3144~3151

    In order to study the transmission paths and pollution source areas of atmospheric PM2.5in different seasons in Taiyuan City, the cluster analysis of daily 48-hour airflow backward trajectories during 2017 to 2018 was conducted, based on the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model and global data assimilation system (GDAS) meteorological element data from National Centers for Environmental Prediction (NCEP), and weighted potential source contribution function (WPSCF) and weighted concentration-weighted trajectory (WCWT) were calculated by combining with hourly pollutant concentration data. The results showed that the sort orders of the average mass concentrations of PM2.5in different seasons of Taiyuan city wereas follows: winter (77.56μg/m3)>autumn (69.89μg/m3)>spring (63.78μg/m3)>summer (45.51μg/m3). It seemed to exist the same source and the secondary conversion process between PM2.5and SO2, NO2and CO. The dominant pollutanttransport trajectories of spring, autumn and winter were short distance and slow moving speed of airflow in the west and northwest directions; and the dominant transport trajectories of summer was in the south and east directions. The potential source areas of PM2.5presented the obvious seasonal variations: in the summer, Taiyuan local and Jinzhong areas contributed greatly to the PM2.5pollutions in Taiyuan City, while the PM2.5pollutions of spring, autumn and winter were mainly affected by the central and northern Shanxi, Luliang, Linfen and Jinzhong areas.

    PM2.5;backward trajectories;cluster analysis;potential source contribution function (PSCF);concentration-weighted trajectory (CWT)

    X513

    A

    1000-6923(2019)08-3144-08

    任浦慧(1992-),女,山西晉城人,山西大學(xué)碩士研究生,主要從事大氣顆粒物的研究.發(fā)表論文1篇.

    2018-12-26

    國家大氣重污染成因與治理攻關(guān)項目(DQGG-05-11);國家公益性行業(yè)(農(nóng)業(yè))科研專項(201103024)

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

    猜你喜歡
    太原市氣流顆粒物
    氣流的威力
    太原市61個村要建污水處理設(shè)施嚴禁直排入河
    太原市園林綠化植物選擇及分析
    太原市六家藥茶企業(yè)獲省級授權(quán)
    太原市第四實驗小學(xué)
    山西電教(2017年3期)2017-10-19 03:39:16
    南平市細顆粒物潛在來源分析
    固體運載火箭變軌發(fā)動機噴管氣流分離研究
    飛片下的空氣形成的“超強高速氣流刀”
    錯流旋轉(zhuǎn)填料床脫除細顆粒物研究
    化工進展(2015年3期)2015-11-11 09:18:15
    基于停留時間分布的氣流床氣化爐通用網(wǎng)絡(luò)模型
    高潮久久久久久久久久久不卡| 国产野战对白在线观看| 日本精品一区二区三区蜜桃| 成年人黄色毛片网站| 校园春色视频在线观看| 777久久人妻少妇嫩草av网站| 极品少妇高潮喷水抽搐| 99国产综合亚洲精品| 一级作爱视频免费观看| 成人精品一区二区免费| 一进一出抽搐gif免费好疼 | 久久天躁狠狠躁夜夜2o2o| 久久久国产成人精品二区 | 亚洲专区中文字幕在线| 亚洲精品粉嫩美女一区| 日韩有码中文字幕| 99国产精品免费福利视频| 久久天躁狠狠躁夜夜2o2o| 人妻久久中文字幕网| 最新在线观看一区二区三区| 久久青草综合色| 一边摸一边做爽爽视频免费| 一本一本久久a久久精品综合妖精| 欧美精品人与动牲交sv欧美| 成人三级做爰电影| 国产精品98久久久久久宅男小说| 另类亚洲欧美激情| 最近最新免费中文字幕在线| 美女扒开内裤让男人捅视频| 黄色成人免费大全| 亚洲精品一二三| 老鸭窝网址在线观看| 99久久综合精品五月天人人| 国产免费男女视频| 啦啦啦免费观看视频1| 精品国产乱码久久久久久男人| 美女 人体艺术 gogo| 久久久久久久久免费视频了| 欧美亚洲日本最大视频资源| 亚洲精品美女久久久久99蜜臀| 桃红色精品国产亚洲av| 脱女人内裤的视频| 国产在线一区二区三区精| 女性生殖器流出的白浆| 亚洲欧美激情综合另类| 99国产精品一区二区三区| 亚洲成人免费av在线播放| 国产精品av久久久久免费| 免费在线观看亚洲国产| 老司机深夜福利视频在线观看| 自线自在国产av| 国产三级黄色录像| 亚洲免费av在线视频| 亚洲av成人一区二区三| 美女午夜性视频免费| 亚洲av第一区精品v没综合| 在线av久久热| 9热在线视频观看99| 一a级毛片在线观看| 狠狠狠狠99中文字幕| 欧美成人免费av一区二区三区 | 天堂√8在线中文| 如日韩欧美国产精品一区二区三区| 9色porny在线观看| 天天躁夜夜躁狠狠躁躁| 国产精品九九99| 欧美日韩一级在线毛片| 美女福利国产在线| 婷婷精品国产亚洲av在线 | 精品久久久精品久久久| 日韩 欧美 亚洲 中文字幕| 国产成人精品久久二区二区91| 成熟少妇高潮喷水视频| 97人妻天天添夜夜摸| 久久久精品免费免费高清| 亚洲五月天丁香| 夫妻午夜视频| 亚洲三区欧美一区| 五月开心婷婷网| 别揉我奶头~嗯~啊~动态视频| 精品一区二区三区av网在线观看| 国产99白浆流出| 国产欧美日韩精品亚洲av| 国产欧美日韩一区二区精品| 国产在视频线精品| 国产一区在线观看成人免费| 亚洲黑人精品在线| 精品久久久久久久久久免费视频 | 国产精品综合久久久久久久免费 | 黑人巨大精品欧美一区二区蜜桃| 久久国产精品大桥未久av| 99久久精品国产亚洲精品| 天天躁日日躁夜夜躁夜夜| 91成人精品电影| 天堂中文最新版在线下载| 精品久久蜜臀av无| 日本黄色视频三级网站网址 | 国产成+人综合+亚洲专区| 亚洲成人免费av在线播放| 亚洲av第一区精品v没综合| 99久久综合精品五月天人人| 亚洲精品在线美女| 国产一区二区三区在线臀色熟女 | 日韩免费高清中文字幕av| 亚洲成人免费av在线播放| 狠狠狠狠99中文字幕| 午夜久久久在线观看| 久久国产精品影院| 天天躁狠狠躁夜夜躁狠狠躁| 巨乳人妻的诱惑在线观看| 黑人欧美特级aaaaaa片| 亚洲欧美精品综合一区二区三区| 国产精品一区二区精品视频观看| 国产真人三级小视频在线观看| av在线播放免费不卡| 午夜亚洲福利在线播放| 国产成人精品无人区| 如日韩欧美国产精品一区二区三区| 黄色丝袜av网址大全| 亚洲自偷自拍图片 自拍| 免费在线观看影片大全网站| 热99久久久久精品小说推荐| 久久久久久久久久久久大奶| 女人被狂操c到高潮| 一级毛片女人18水好多| 国产精品一区二区精品视频观看| 巨乳人妻的诱惑在线观看| 天天添夜夜摸| 波多野结衣一区麻豆| 久久久久久亚洲精品国产蜜桃av| 亚洲在线自拍视频| 不卡av一区二区三区| 久久久久久人人人人人| 一进一出好大好爽视频| 中文字幕人妻丝袜制服| 中文欧美无线码| 久久草成人影院| 丁香六月欧美| 黄片小视频在线播放| 免费观看人在逋| 久久久国产成人免费| 国产高清国产精品国产三级| 亚洲色图av天堂| 国产男靠女视频免费网站| 国产一区二区激情短视频| 久久性视频一级片| 无限看片的www在线观看| 两个人免费观看高清视频| 18禁裸乳无遮挡免费网站照片 | 亚洲中文字幕日韩| 久久久国产精品麻豆| 女人精品久久久久毛片| 黄网站色视频无遮挡免费观看| 99re6热这里在线精品视频| 无人区码免费观看不卡| 国产精品1区2区在线观看. | 久久国产精品人妻蜜桃| 日韩欧美国产一区二区入口| 1024视频免费在线观看| 99久久99久久久精品蜜桃| 黄网站色视频无遮挡免费观看| 成人特级黄色片久久久久久久| 美女扒开内裤让男人捅视频| 免费观看精品视频网站| 欧美亚洲日本最大视频资源| 精品一品国产午夜福利视频| 美女高潮喷水抽搐中文字幕| 露出奶头的视频| 免费高清在线观看日韩| 亚洲熟女精品中文字幕| 亚洲精品国产色婷婷电影| 欧美乱色亚洲激情| 午夜精品国产一区二区电影| 一进一出好大好爽视频| 亚洲av日韩在线播放| 操美女的视频在线观看| 人人妻人人添人人爽欧美一区卜| 欧美日韩精品网址| 涩涩av久久男人的天堂| 女人久久www免费人成看片| 国产av又大| 成人18禁在线播放| www日本在线高清视频| 亚洲国产精品sss在线观看 | 男人舔女人的私密视频| 男女下面插进去视频免费观看| 久久人妻av系列| 久久久久久人人人人人| 亚洲欧美一区二区三区黑人| 脱女人内裤的视频| 国产av一区二区精品久久| a在线观看视频网站| 王馨瑶露胸无遮挡在线观看| 黄色女人牲交| 久9热在线精品视频| 欧美性长视频在线观看| 日韩欧美一区二区三区在线观看 | 亚洲精品成人av观看孕妇| 99精品欧美一区二区三区四区| 精品无人区乱码1区二区| 色老头精品视频在线观看| 国产成人系列免费观看| 日韩有码中文字幕| 水蜜桃什么品种好| 女性生殖器流出的白浆| 黑人欧美特级aaaaaa片| 日韩 欧美 亚洲 中文字幕| 香蕉国产在线看| 亚洲一区二区三区不卡视频| 一本一本久久a久久精品综合妖精| 欧美黄色淫秽网站| 久久婷婷成人综合色麻豆| 少妇裸体淫交视频免费看高清 | 我的亚洲天堂| 91成人精品电影| 国产精华一区二区三区| 久久热在线av| 国产99白浆流出| 老熟妇乱子伦视频在线观看| 午夜精品久久久久久毛片777| 人人澡人人妻人| 午夜视频精品福利| 久久香蕉国产精品| 99精品久久久久人妻精品| 国产成人免费无遮挡视频| 18禁观看日本| 91成年电影在线观看| 免费在线观看黄色视频的| 一级毛片女人18水好多| 免费在线观看影片大全网站| 国产欧美日韩精品亚洲av| 黄色片一级片一级黄色片| 丁香六月欧美| 好男人电影高清在线观看| 91麻豆精品激情在线观看国产 | 欧美日韩中文字幕国产精品一区二区三区 | 久久久久精品人妻al黑| 亚洲色图 男人天堂 中文字幕| 亚洲avbb在线观看| a级毛片黄视频| 久久久精品免费免费高清| 久久精品91无色码中文字幕| av一本久久久久| 亚洲五月色婷婷综合| 国产精品免费一区二区三区在线 | 国产高清videossex| 好男人电影高清在线观看| 最近最新免费中文字幕在线| 老司机亚洲免费影院| 国产一区二区三区综合在线观看| svipshipincom国产片| 极品少妇高潮喷水抽搐| 中文字幕人妻熟女乱码| 黄网站色视频无遮挡免费观看| 51午夜福利影视在线观看| 亚洲精品久久午夜乱码| 国产淫语在线视频| 在线av久久热| 欧美+亚洲+日韩+国产| a在线观看视频网站| 王馨瑶露胸无遮挡在线观看| www.999成人在线观看| 欧美一级毛片孕妇| 亚洲少妇的诱惑av| 国产极品粉嫩免费观看在线| 国产无遮挡羞羞视频在线观看| 90打野战视频偷拍视频| 精品少妇久久久久久888优播| 国产主播在线观看一区二区| 精品久久蜜臀av无| 亚洲成人免费av在线播放| 制服诱惑二区| 国产亚洲一区二区精品| 亚洲一码二码三码区别大吗| 国产免费男女视频| 日韩视频一区二区在线观看| 免费观看精品视频网站| 男女床上黄色一级片免费看| 高清在线国产一区| 在线观看免费视频网站a站| 国产麻豆69| 999久久久精品免费观看国产| 国产成人精品久久二区二区91| 色播在线永久视频| av有码第一页| 精品国产乱子伦一区二区三区| 国产高清激情床上av| 日本wwww免费看| 999精品在线视频| 欧美精品人与动牲交sv欧美| 激情在线观看视频在线高清 | 老熟妇乱子伦视频在线观看| 天堂中文最新版在线下载| 男男h啪啪无遮挡| 日韩欧美三级三区| 正在播放国产对白刺激| 国产蜜桃级精品一区二区三区 | av视频免费观看在线观看| 欧美黄色淫秽网站| 亚洲成人免费av在线播放| 中文字幕人妻熟女乱码| 欧美 日韩 精品 国产| 久久亚洲真实| 在线十欧美十亚洲十日本专区| 免费高清在线观看日韩| 一区二区三区国产精品乱码| av有码第一页| 国产精品电影一区二区三区 | 亚洲av成人一区二区三| 女性生殖器流出的白浆| 日日夜夜操网爽| 黄片小视频在线播放| 欧美另类亚洲清纯唯美| 国产真人三级小视频在线观看| 夜夜夜夜夜久久久久| 久久久久久人人人人人| av福利片在线| 国产成人欧美| 亚洲色图 男人天堂 中文字幕| 国产精品98久久久久久宅男小说| av网站免费在线观看视频| netflix在线观看网站| 日本黄色日本黄色录像| 亚洲 欧美一区二区三区| 久久天堂一区二区三区四区| 国产高清国产精品国产三级| 日韩精品免费视频一区二区三区| 校园春色视频在线观看| 精品国产一区二区三区四区第35| 亚洲免费av在线视频| 亚洲全国av大片| 午夜精品在线福利| 妹子高潮喷水视频| 国产精品偷伦视频观看了| 精品少妇一区二区三区视频日本电影| 一级,二级,三级黄色视频| 精品久久久精品久久久| 极品教师在线免费播放| 一本大道久久a久久精品| 90打野战视频偷拍视频| 午夜91福利影院| 久久久水蜜桃国产精品网| 我的亚洲天堂| 亚洲av成人av| 久久 成人 亚洲| 天堂动漫精品| 亚洲 欧美一区二区三区| 在线免费观看的www视频| 国产成人啪精品午夜网站| 日韩 欧美 亚洲 中文字幕| 国产97色在线日韩免费| 国产在线一区二区三区精| 国产黄色免费在线视频| 美女 人体艺术 gogo| 香蕉国产在线看| 久久精品亚洲精品国产色婷小说| 巨乳人妻的诱惑在线观看| 91国产中文字幕| 99久久99久久久精品蜜桃| 欧美日韩黄片免| 99国产综合亚洲精品| 国产亚洲欧美精品永久| 亚洲午夜精品一区,二区,三区| 成年动漫av网址| 18在线观看网站| 看免费av毛片| 精品视频人人做人人爽| 99香蕉大伊视频| 超碰成人久久| 国产一区二区三区综合在线观看| 亚洲成国产人片在线观看| 中文亚洲av片在线观看爽 | 国产亚洲欧美精品永久| 飞空精品影院首页| 亚洲av欧美aⅴ国产| 日日夜夜操网爽| 一a级毛片在线观看| 视频区图区小说| av免费在线观看网站| 午夜福利免费观看在线| 正在播放国产对白刺激| 涩涩av久久男人的天堂| 成人国语在线视频| 国产成人欧美| 国产成人av教育| 男人舔女人的私密视频| 亚洲一卡2卡3卡4卡5卡精品中文| 超碰97精品在线观看| 国产熟女午夜一区二区三区| 三上悠亚av全集在线观看| 日本一区二区免费在线视频| 99久久人妻综合| 日韩一卡2卡3卡4卡2021年| 老熟女久久久| 精品久久蜜臀av无| avwww免费| 免费av中文字幕在线| 久久精品国产亚洲av香蕉五月 | 国产欧美日韩精品亚洲av| 我的亚洲天堂| 国产精品.久久久| 一区二区三区国产精品乱码| 少妇裸体淫交视频免费看高清 | 久久人人爽av亚洲精品天堂| 亚洲片人在线观看| 成人18禁在线播放| 巨乳人妻的诱惑在线观看| 国产欧美日韩综合在线一区二区| 超碰97精品在线观看| 午夜精品久久久久久毛片777| 女人被躁到高潮嗷嗷叫费观| 午夜老司机福利片| 新久久久久国产一级毛片| 久久久国产成人精品二区 | 亚洲av日韩在线播放| 亚洲国产精品合色在线| 中亚洲国语对白在线视频| 中文亚洲av片在线观看爽 | 国产男靠女视频免费网站| 国产区一区二久久| 精品一品国产午夜福利视频| av天堂久久9| 天堂动漫精品| 美女福利国产在线| 性色av乱码一区二区三区2| 欧美精品啪啪一区二区三区| 丝袜人妻中文字幕| 久久婷婷成人综合色麻豆| 女警被强在线播放| 村上凉子中文字幕在线| 午夜成年电影在线免费观看| 免费在线观看完整版高清| 日日夜夜操网爽| 欧美日韩中文字幕国产精品一区二区三区 | 在线天堂中文资源库| 国产激情欧美一区二区| 一区二区三区精品91| 国产视频一区二区在线看| tocl精华| ponron亚洲| 精品久久久久久电影网| 搡老岳熟女国产| 美女扒开内裤让男人捅视频| 人妻一区二区av| av国产精品久久久久影院| 老司机影院毛片| 中文字幕精品免费在线观看视频| 久久精品91无色码中文字幕| 亚洲精品在线美女| 亚洲人成电影免费在线| 亚洲熟女毛片儿| 欧美精品人与动牲交sv欧美| 色精品久久人妻99蜜桃| 人人澡人人妻人| 亚洲国产精品sss在线观看 | 欧美老熟妇乱子伦牲交| 高潮久久久久久久久久久不卡| 五月开心婷婷网| 无限看片的www在线观看| 99国产精品一区二区蜜桃av | 黄色片一级片一级黄色片| av天堂在线播放| 高清av免费在线| 国产不卡av网站在线观看| 亚洲色图综合在线观看| 国产精品av久久久久免费| 欧美日韩av久久| 一区二区三区精品91| 91字幕亚洲| 两人在一起打扑克的视频| 国产主播在线观看一区二区| 黄色a级毛片大全视频| 亚洲人成77777在线视频| 免费在线观看视频国产中文字幕亚洲| 又大又爽又粗| 老熟妇仑乱视频hdxx| 国产亚洲欧美98| 精品视频人人做人人爽| 亚洲精品久久午夜乱码| 国产精品av久久久久免费| 国产精品一区二区精品视频观看| 亚洲色图 男人天堂 中文字幕| 精品人妻1区二区| 大片电影免费在线观看免费| 天堂动漫精品| 精品久久久精品久久久| 一边摸一边抽搐一进一小说 | 国产单亲对白刺激| 久久天堂一区二区三区四区| 中文字幕高清在线视频| 99国产精品一区二区蜜桃av | 欧美丝袜亚洲另类 | 一进一出抽搐gif免费好疼 | 欧美国产精品va在线观看不卡| 中亚洲国语对白在线视频| 亚洲精品国产区一区二| 亚洲成国产人片在线观看| 国产免费av片在线观看野外av| 日日夜夜操网爽| 国产亚洲精品久久久久5区| 精品福利观看| 大香蕉久久成人网| 村上凉子中文字幕在线| 亚洲欧洲精品一区二区精品久久久| 国产欧美亚洲国产| 久久精品亚洲精品国产色婷小说| 超色免费av| 欧美+亚洲+日韩+国产| 美女扒开内裤让男人捅视频| 亚洲av成人av| 欧美国产精品va在线观看不卡| 日韩大码丰满熟妇| av中文乱码字幕在线| 午夜成年电影在线免费观看| 国产亚洲欧美在线一区二区| www日本在线高清视频| 制服诱惑二区| 日韩成人在线观看一区二区三区| 国产精品一区二区免费欧美| 久久中文字幕人妻熟女| 日韩中文字幕欧美一区二区| 久久久久久亚洲精品国产蜜桃av| 久久天堂一区二区三区四区| 在线观看舔阴道视频| 久久精品91无色码中文字幕| 国产精品影院久久| www.熟女人妻精品国产| videosex国产| 亚洲色图av天堂| 国产极品粉嫩免费观看在线| 国产淫语在线视频| 中文字幕av电影在线播放| 狠狠狠狠99中文字幕| 欧美亚洲日本最大视频资源| 中文字幕色久视频| 久久久久精品人妻al黑| 不卡av一区二区三区| 国产av又大| 欧美在线黄色| 免费在线观看日本一区| 亚洲av熟女| 黄色 视频免费看| 亚洲综合色网址| 日本精品一区二区三区蜜桃| 热99国产精品久久久久久7| 欧美精品一区二区免费开放| 亚洲国产精品sss在线观看 | 午夜91福利影院| 1024视频免费在线观看| 老熟女久久久| 久久人妻av系列| 视频区欧美日本亚洲| 国产精品1区2区在线观看. | 成人特级黄色片久久久久久久| 性少妇av在线| 人妻一区二区av| 午夜福利影视在线免费观看| 欧美成狂野欧美在线观看| 亚洲伊人色综图| 成人18禁高潮啪啪吃奶动态图| 黄色视频,在线免费观看| 一区二区三区激情视频| 国产区一区二久久| 99国产精品免费福利视频| 亚洲国产精品合色在线| 99香蕉大伊视频| 一区二区日韩欧美中文字幕| 男男h啪啪无遮挡| 最新的欧美精品一区二区| 十八禁高潮呻吟视频| 丰满迷人的少妇在线观看| 国产精品一区二区在线观看99| 午夜两性在线视频| 成年女人毛片免费观看观看9 | 午夜久久久在线观看| 国产午夜精品久久久久久| 亚洲欧美一区二区三区久久| 久久性视频一级片| 亚洲精品在线观看二区| 人人澡人人妻人| 69av精品久久久久久| 一进一出抽搐gif免费好疼 | 久久香蕉精品热| 最新在线观看一区二区三区| 老汉色av国产亚洲站长工具| 手机成人av网站| 欧美日韩福利视频一区二区| 精品国产美女av久久久久小说| 精品卡一卡二卡四卡免费| 91国产中文字幕| 天天影视国产精品| 中文字幕av电影在线播放| 视频在线观看一区二区三区| videos熟女内射| 国产99久久九九免费精品| 欧美精品人与动牲交sv欧美| 极品少妇高潮喷水抽搐| 久久久久久久久免费视频了| 国产精品久久久久久人妻精品电影| 亚洲精品在线观看二区| www.熟女人妻精品国产| 欧美日韩乱码在线| 女人被躁到高潮嗷嗷叫费观| 在线播放国产精品三级| 俄罗斯特黄特色一大片| 国产高清视频在线播放一区| 国产精品欧美亚洲77777| 精品国产乱码久久久久久男人| 女人被躁到高潮嗷嗷叫费观| 国产片内射在线| 人人妻人人澡人人爽人人夜夜| 女性被躁到高潮视频| 久久性视频一级片| 国产精品一区二区精品视频观看| 久久天堂一区二区三区四区|