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

    基于星地同步觀測(cè)的華北平原中部背景地區(qū)冬季霾污染過程

    2016-10-13 19:45:20王新輝陶明輝王子峰陳良富李莘莘
    中國(guó)環(huán)境科學(xué) 2016年6期
    關(guān)鍵詞:沙塵能見度氣溶膠

    王新輝,蘇 林,陶明輝,王子峰,陳良富,李莘莘,汪 洋

    ?

    基于星地同步觀測(cè)的華北平原中部背景地區(qū)冬季霾污染過程

    王新輝1,2,蘇 林*,陶明輝1,王子峰1,陳良富1,李莘莘1,汪 洋1,2

    (1.中國(guó)科學(xué)院遙感與數(shù)字地球研究所,遙感科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100101;2.中國(guó)科學(xué)院大學(xué),北京 100049)

    利用地基觀測(cè)結(jié)果和多源衛(wèi)星遙感觀測(cè),結(jié)合氣象數(shù)據(jù)及HYSPLIT4后向軌跡模式,對(duì)華北平原中部背景地區(qū)(河南省鄭州市中牟縣東南郊)冬季霾事件的污染物特征和形成過程進(jìn)行分析.綜合觀測(cè)時(shí)間為2014年12月13日~2015年1月16日,共有5次霾過程,占觀測(cè)總天數(shù)的82%.地面監(jiān)測(cè)結(jié)果顯示,不同的污染過程污染物濃度變化曲線相似,O3濃度在清潔天濃度較高;NO、SO2、PM10、PM2.5呈較強(qiáng)正相關(guān)性,NOx、SO2與 PM10相關(guān)系數(shù)0.64、0.57,與PM2.5相關(guān)系數(shù)0.56、0.45;近地面污染物以細(xì)粒子污染物為主,其中又以氣態(tài)污染物二次生成的細(xì)粒子為主.AMPLE地基激光雷達(dá)和CALIPSO數(shù)據(jù)表明,華北平原霾層中上部受浮塵影響,以粗粒子污染物為主.氣象探空數(shù)據(jù)表明該地區(qū)冬季大氣對(duì)流層穩(wěn)定利于霾的維持,假相當(dāng)位溫垂直差、指數(shù)、露點(diǎn)差與能見度相關(guān)系數(shù)分別為0.52、0.56和 0.38.分析近地面風(fēng)速風(fēng)向?qū)策^程的影響表明,該地區(qū)冬季以南方向靜小風(fēng)為主,風(fēng)速與能見度相關(guān)系數(shù)為0.32 ,PM1受東北方向污染源影響,PM1~2.5及PM2.5~10受西北方向污染源影響;結(jié)合高空風(fēng)場(chǎng)分析,霾過程1受西北浮塵影響,霾過程5受南來水汽影響.通過后向軌跡分析,該地區(qū)冬季的低空污染傳輸主要來自東北和西北方向,其中東北方向區(qū)域傳輸來自河北與山東,占來源比例的14%,近距離污染傳輸主要來自站點(diǎn)以西的鄭州、洛陽,約占來源比例的33%.

    霾過程;污染物特征;氣象因素;區(qū)域輸送

    霾是由于空氣中大量微小氣溶膠粒子積聚在大氣中,使得空氣變渾濁,致使水平能見度下降至10km以下的一種天氣現(xiàn)象[1-2].近年來,我國(guó)生活與工業(yè)的污染排放增加迅速,城市光化學(xué)煙霧事件變多,霾天氣急劇增多,尤其在長(zhǎng)三角、珠三角、京津冀等城市密集區(qū),區(qū)域性污染日益嚴(yán)重,極端污染事件頻發(fā)[3-5].國(guó)外學(xué)者對(duì)霾污染的研究較早,Andreae等[6]1988年在亞馬遜盆地進(jìn)行綜合觀測(cè),對(duì)1~4km霾層結(jié)構(gòu)分析發(fā)現(xiàn),霾層中污染氣體濃度高.早期對(duì)霾的研究主要通過消光系數(shù)、化學(xué)組分、粒徑與顆粒物濃度等分析,研究顆粒物特性以及與氣溶膠混合狀態(tài)的相互作用,指出細(xì)粒子污染容易形成灰霾天氣[7-8].我國(guó)灰霾的概念于2002年提出后,開始系統(tǒng)研究[9]. Zhang等[10]根據(jù)49年681個(gè)站點(diǎn)的能見度,將我國(guó)霾區(qū)劃分為華北平原、長(zhǎng)三角、珠三角和四川盆地.目前國(guó)內(nèi)對(duì)霾的研究主要有如下方向:通過數(shù)值分析模擬與大氣模式研究并預(yù)測(cè)霾的發(fā)生發(fā)展與污染程度[11-12],通過地基觀測(cè)獲得顆粒物濃度、氣溶膠光學(xué)特性、化學(xué)組分研究霾的時(shí)空變化、污染特性或揭示霾的形成機(jī)制[13-15],以及通過氣象數(shù)據(jù)與參數(shù)研究氣象條件對(duì)霾形成發(fā)展的影響[16-17].隨衛(wèi)星遙感技術(shù)發(fā)展,AQUA、CALIPSO、AURA、PARASOL等系列衛(wèi)星的聯(lián)合觀測(cè)為霾的研究提供新視角[18].如H?nel等[19]結(jié)合CALIPSO數(shù)據(jù)對(duì)GAW站點(diǎn)霾層分析,指出東亞地區(qū)低層霾多為人為細(xì)顆粒氣溶膠,高層霾是沙塵、生物質(zhì)燃煙以及工業(yè)污染物的混合物.Tao等[4]結(jié)合MODIS火點(diǎn)數(shù)據(jù)、氣溶膠產(chǎn)品、CALIPSO數(shù)據(jù)分析東亞地區(qū)夏季棕色霾,指出沙塵傳輸、生物質(zhì)燃燒與工業(yè)污染以及夏季高濕度的天氣條件是其成因.

    現(xiàn)今有關(guān)霾的研究多針對(duì)京津冀、珠三角、長(zhǎng)三角等污染嚴(yán)重地區(qū)統(tǒng)計(jì)分析,或針對(duì)某次重霾事件進(jìn)行分析與追因,對(duì)背景地區(qū)的霾過程分析較少,且依據(jù)地基觀測(cè)數(shù)據(jù)的研究占多數(shù).通過站點(diǎn)進(jìn)行地基觀測(cè)具有數(shù)據(jù)精度高、時(shí)間分辨率高等優(yōu)勢(shì),但地基站點(diǎn)數(shù)量有限,無法進(jìn)行連續(xù)空間分析,且觀測(cè)結(jié)果容易受周圍源影響,衛(wèi)星觀測(cè)恰好彌補(bǔ)這一不足,可實(shí)現(xiàn)霾的區(qū)域分析以及污染傳輸?shù)谋O(jiān)控.華北平原是我國(guó)四大霾區(qū)之一,本文選取華北平原中部背景地區(qū)進(jìn)行了為期一個(gè)月的冬季外場(chǎng)地基觀測(cè)試驗(yàn),并同時(shí)結(jié)合氣象數(shù)據(jù)、衛(wèi)星數(shù)據(jù)和后向軌跡模型對(duì)該背景地區(qū)的霾污染物特征、氣象條件影響、形成過程以及污染傳輸源頭進(jìn)行綜合的分析與研究.

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

    1.1 站點(diǎn)描述

    背景觀測(cè)站點(diǎn)位于河南鄭州市的中牟縣東南郊.站點(diǎn)經(jīng)緯度坐標(biāo)為114.04°E,34.71°N.站點(diǎn)所處地勢(shì)平緩、四周是大面積麥田,觀測(cè)數(shù)據(jù)受人為源影響少.東偏北方向30km處為開封市區(qū),西北方向35km處為鄭州市區(qū).聯(lián)合觀測(cè)時(shí)間段為2014年12月13日~2015年1月16日.

    1.2 地基觀測(cè)與氣象數(shù)據(jù)

    此次地基觀測(cè)數(shù)據(jù)包括近地表的氣象五參數(shù)數(shù)據(jù)(溫度、濕度、風(fēng)速、風(fēng)向、氣壓)、近地面能見度數(shù)據(jù)、氣溶膠光學(xué)厚度數(shù)據(jù)、近地表顆粒物濃度數(shù)據(jù)和AMPLE(多波長(zhǎng)拉曼偏振大氣探測(cè)激光雷達(dá))探測(cè)數(shù)據(jù).氣溶膠光學(xué)厚度由CE318太陽光度計(jì)測(cè)得,由于云或水汽的影響,有效數(shù)據(jù)反演結(jié)果存在一定時(shí)間段的缺失. AMPLE地基探測(cè)激光雷達(dá)2012年由中國(guó)航天科技集團(tuán)公司九院北京遙測(cè)技術(shù)研究所研制,空間分辨率為3.75m,探測(cè)距離達(dá)到30km.

    污染氣體濃度的觀測(cè)數(shù)據(jù)來自中牟市環(huán)境監(jiān)測(cè)中心提供的空氣質(zhì)量數(shù)據(jù).溫度、濕度以及風(fēng)速風(fēng)向等的探空資料來自于懷俄明大學(xué)大氣科學(xué)學(xué)院提供的無線探測(cè)資料,時(shí)間分辨率為1次/12h,對(duì)應(yīng)中國(guó)地區(qū)時(shí)間為8:00與20:00.風(fēng)場(chǎng)數(shù)據(jù)來自NCEP (Nationgal Centers for Environmental Prediction),分辨率為1°,用于分析大范圍的氣象條件[20],可分析潛在的水汽傳輸以及沙塵傳輸?shù)挠绊?

    1.3 衛(wèi)星數(shù)據(jù)集

    本文采用MODIS(中分辨率成像光譜儀)550nm波段10km AOD(氣溶膠光學(xué)厚度)產(chǎn)品、OMI(臭氧監(jiān)測(cè)傳感器)的UVAI(紫外氣溶膠指數(shù))以及CALIPSO產(chǎn)品.其中AOD與真實(shí)值有(±0.05±15%)的誤差,可用來分析區(qū)域污染的變化[21].UVAI可用于探測(cè)沙塵和生物燃燒,監(jiān)測(cè)沙塵和煙羽的傳輸[22].UVAI對(duì)紫外吸收型氣溶膠的高度非常敏感(>2km),對(duì)表面類型不夠敏感[23].本文采用OMI L2格網(wǎng)數(shù)據(jù),UVAI分辨率為0.25°.CALIPSO產(chǎn)品可提供氣溶膠和云層的垂直廓線、后向散射和退偏比[24],其中532nm波段的線性退偏比對(duì)沙塵等非球形粒子尤其敏感[25],同時(shí)可以提供氣溶膠類型產(chǎn)品,并將氣溶膠分成6類,包括煤煙、沙塵、污染沙塵、清潔和污染的大陸型,以及清潔的海洋型[26].

    1.4 軌跡模式HYSPLIT4

    采用美國(guó)NOAA研制的軌跡模式HYSPLIT4定性分析此次霾天氣污染物的輸送過程.HYSPLIT4 是一種歐拉和拉格朗日混合型的計(jì)算模式,包含多種物理過程,可針對(duì)不同類型排放源,進(jìn)行較完整的輸送擴(kuò)散和沉降過程模擬,廣泛應(yīng)用于大氣污染物輸送研究[27-28].

    2 污染物濃度特征分析

    2.1 霾污染過程

    根據(jù)霾的定義[1-2],結(jié)合日均能見度小于10km的日期,霾過程可分為5個(gè),分別是12月13~15日,12月18~20日,12月23~30日,1月2~5日,1月7~15日.霾污染天數(shù)為27d,占整個(gè)觀測(cè)期間天數(shù)的82%.霧霾天的定義是相對(duì)濕度日均值大于80%的霾天[29],反之為灰霾天.根據(jù)圖1所示,觀測(cè)期間日均相對(duì)濕度均低于80%,為灰霾天氣.依據(jù)AQI(空氣質(zhì)量指數(shù))對(duì)中國(guó)大陸空氣質(zhì)量等級(jí)的劃分[30],觀測(cè)期間重污染天數(shù)較多,占觀測(cè)天數(shù)的46%.其中,12月18日,以及1月的3~5日、8日、10日和15日出現(xiàn)六級(jí)重度污染(AQI值大于300),五級(jí)重度污染(AQI值為200~300)有9d.達(dá)到優(yōu)良級(jí)別(AQI值小于100)的有7d,占觀測(cè)天數(shù)的18%.

    從圖1可以看出,過程2和過程3霾的發(fā)生變化過程具有明顯的突發(fā)性,可在2~3h內(nèi)形成,并經(jīng)過進(jìn)一步的積聚形成重污染.過程2在17日16:00~19:00能見度迅速從14km降低至9km,次日10時(shí)達(dá)最低值.過程3在22日16:00~19:00能見度從14.3km降低至7.5km次日10:00達(dá)到能見度最低.與過程2和過程3霾發(fā)生的突發(fā)性不同,過程4和過程5霾形成前從能見度開始下降到霾形成經(jīng)過相對(duì)長(zhǎng)的時(shí)間,霾形成后較為穩(wěn)定,一直維持低能見度狀態(tài).該地區(qū)能見度存在著規(guī)律的日內(nèi)變化,由于邊界層高度升高、水平對(duì)流和垂直交換增強(qiáng)以及溫度升高、濕度降低等原因[31-32],能見度白天升高,相反在夜晚降低;幾個(gè)過程中霾發(fā)生時(shí)間均在下午或夜間,契合能見度的日變化.

    2.2 污染物特征及變化分析

    由圖2可見,顆粒物PM1、PM2.5、PM10的平均濃度分別為90,101,128μg/m3,SO2、NO、O3、CO平均濃度分別為60,109,20μg/m3、2.74mg/m3.O3日均值峰值主要出現(xiàn)在污染過程之前,如12月15日、12月21日、12月26日、12月31日、1月6日,日均值最高值在12月26日達(dá)到54μg/m3.O3濃度在霾過程開始時(shí)下降,原因可能是O3生成與前體物NO和VOC有關(guān)[33],而當(dāng)某些還原性污染物如 NO濃度較大時(shí)[34],O3易與其發(fā)生反應(yīng),加之霧霾天太陽光照弱,光化學(xué)反應(yīng)減弱[35],使得O3濃度在霾天下降.AOD與PM10、PM2.5的變化有較好的相關(guān),相關(guān)系數(shù)分別為0.4212、0.5237.PM10、PM2.5、SO2、NO變化過程具有一致性.CO與其他5種成分也具有相關(guān)變化,但其變化和其他成分的變化相比大多有時(shí)間上的延遲,日均值的最高值在1月15日達(dá)到5.93mg/m3.通過對(duì)污染物的小時(shí)均值計(jì)算,在此過程中NO與SO2濃度相關(guān)系數(shù)為0.6312, NO、SO2與PM10濃度的相關(guān)系數(shù)分別為0.6426、0.5666,與PM2.5濃度的相關(guān)系數(shù)分別為0.5592、0.4489,4種污染物呈現(xiàn)明顯正相關(guān)性,說明四者具有較高的同源性,SO2和NO的爆發(fā)長(zhǎng)伴隨著PM的爆發(fā)[36],同時(shí)SO2、NO、O3、CO與污染氣體二次生成細(xì)顆粒相關(guān)性高[37],說明顆粒物中二次生成細(xì)顆粒物貢獻(xiàn)較大.

    圖3是觀測(cè)站點(diǎn)7種大氣成分濃度在不同天氣條件下(清潔天和5個(gè)霾污染過程)24h的變化.PM1、PM2.5與PM10的日變化相關(guān)性較好,呈現(xiàn)雙峰變化,白天在9:00左右達(dá)到峰值,夜間由于大氣邊界層變低,氣流速度變緩也出現(xiàn)一個(gè)峰值.O3日變化呈現(xiàn)單峰正態(tài)變化,隨著日照變化,光化學(xué)反應(yīng)變強(qiáng),峰值出現(xiàn)在15:00.SO2主要來自燃煤[35],圖3顯示其峰值出現(xiàn)在11:00~12:00左右,另有兩個(gè)污染過程在7點(diǎn)左右也出現(xiàn)一個(gè)峰值.NO濃度峰值在9:00~10:00左右,次高峰出現(xiàn)在19:00.CO主要來源是燃煤和機(jī)動(dòng)車尾氣[38], CO日變化出現(xiàn)多個(gè)峰值,最大的峰值出現(xiàn)在7:00~8:00上班高峰期,次高峰出現(xiàn)在15:00~ 17:00幼兒園與小學(xué)放學(xué)期間以及19:00~23:00城市夜生活間, 與汽車尾氣排放有較大關(guān)系.上述分析表明,觀測(cè)背景點(diǎn)空氣質(zhì)量與污染物變化一定程度上受城市影響,這也是由于華北平原地勢(shì)平坦開闊一定程度上利于污染的擴(kuò)散和傳輸.另外不同天氣狀況下污染物的變化曲線差別不大,O3的變化曲線最為穩(wěn)定.除了O3濃度在清潔天相對(duì)霾污染過程濃度高,CO區(qū)別較小,其他污染物清潔天濃度低于霾天.

    結(jié)合地基測(cè)得污染物濃度數(shù)據(jù)、地基雷達(dá)數(shù)據(jù)以及CALIPSO衛(wèi)星影像,對(duì)顆粒物濃度時(shí)序變化和垂直分布進(jìn)一步分析.從圖4可以看出,該地區(qū)1月份顆粒物污染程度遠(yuǎn)大于12月份;除了霾過程1,其他4個(gè)霾過程中日均值PM2.5/PM10在0.7以上,清潔天的值均小于0.7,說明該地區(qū)霾過程中近地表的顆粒物污染貢獻(xiàn)主要來自細(xì)粒子;同時(shí)PM1/PM2.5值大多在0.8以上,PM1主要是由污染氣體二次生成顆粒污染物[39],說明該地區(qū)污染氣體二次生成細(xì)顆粒污染物占較大比重.霾過程1粗顆粒貢獻(xiàn)較大,其中13日、15日PM2.5~10/PM10分別為48.7%、40%,可能存在揚(yáng)塵或沙塵傳輸,進(jìn)一步通過雷達(dá)影像(圖5)和UVAI紫外氣溶膠指數(shù)(圖9)驗(yàn)證.

    過程1AMPLE地基雷達(dá)探測(cè)532nm波段退偏比時(shí)間序列圖如圖5所示.退偏比越高顆粒物越不規(guī)則,粒徑越大[40],反之說明粒徑越小.根據(jù)退偏比信息可以看到,整層的顆粒物分布情況和濃度變化.13日近地面的顆粒物退偏比在0.2左右,粗粒子較多,1km上空有較厚的粗顆粒層;14日近地面退偏比降低,粗顆粒減少同時(shí)霾層變厚;15日整層退偏比都有所增加.從整段期間來看,13~14日4km左右的高空有明顯的浮塵沉降過程,14日1km處浮塵有上升的過程,與沉降的浮塵結(jié)合在一起在2~3km處積聚后在15日沉降,導(dǎo)致近地面粗顆粒物濃度增加.

    地基雷達(dá)數(shù)據(jù)獲得點(diǎn)的整層的時(shí)序變化,為了獲取華北霾天區(qū)域性垂直霾層顆粒物污染情況,需用衛(wèi)星數(shù)據(jù)CALIPSO.選擇衛(wèi)星過境華北平原中部且污染嚴(yán)重的1月4日、8日、10日、16日分析華北平原垂直霾層特征.圖6顯示,盡管地面監(jiān)測(cè)顯示顆粒物以細(xì)粒子為主,但光學(xué)觀測(cè)整層中粗粒子占主導(dǎo).霾層中上部受浮塵主導(dǎo),后向散射系數(shù)較大,粗粒子較多;霾層層頂高度最高達(dá)到8~9km,沙塵層厚度大于5km;上層氣溶膠類型主要是沙塵與污染沙塵型氣溶膠,邊界層內(nèi)存在煤煙型和污染大陸型以及污染沙塵氣溶膠組成的混合氣溶膠.

    3 氣象影響與霾過程分析

    3.1 近地表風(fēng)

    霾分布區(qū)域內(nèi)的表面風(fēng)速作為水平動(dòng)力因素影響污染物被大氣稀釋的程度和擴(kuò)散的范圍[41-42].表面風(fēng)速偏大時(shí),污染物以及水汽等向區(qū)域外輸送加強(qiáng),反之有利于污染物的積聚以及霾的發(fā)展;另一方面,風(fēng)速較大時(shí)可吹起地面塵土等顆粒,在局地范圍內(nèi)易形成浮塵或揚(yáng)沙天氣.風(fēng)向能反映出污染源周圍受影響的方位及污染物的來向和去向.能見度大小是霾的表征,從圖7可見,5個(gè)霾過程發(fā)生前期以靜風(fēng)(<0.5m/s)為主,靜風(fēng)風(fēng)頻為10.7%.霾的消散伴隨著較大風(fēng)速,霾過程中能見度的變化均伴隨風(fēng)速的變化,如12月19日、12月14日、12月18日和1月1日.靜穩(wěn)無風(fēng)天氣是霧霾積聚維持的重要環(huán)境條件.

    由圖8可知,PM1、PM1~2.5、PM2.5~10主要集中在靜小風(fēng)2m/s范圍內(nèi),約占67.9%,且風(fēng)速小于2m/s多為南風(fēng)向,約占79%;風(fēng)速越小,PM1質(zhì)量濃度越大,當(dāng)風(fēng)速大于4m/s 時(shí),PM1~2.5和PM2.5~10濃度有增加趨勢(shì).來自東北風(fēng)向的顆粒物數(shù)濃度比例較大,約為26.6%;PM1高值多出現(xiàn)在東北方向,PM1~2.5、PM2.5~10高值多來自西和西北方向.綜上說明,觀測(cè)期間霾污染過程中,顆粒物多為本地積聚,此外PM1細(xì)顆粒物主要受東北方向污染源影響,而PM1~2.5、PM2.5~10主要受西和西北方向的影響.

    3.2 高空風(fēng)場(chǎng)

    結(jié)合衛(wèi)星數(shù)據(jù)分析高空風(fēng)場(chǎng)水平傳輸對(duì)觀測(cè)站點(diǎn)及華北平原地區(qū)霾發(fā)展過程中的影響,以12月12~15日(包含霾過程1)與1月6日~16日(包含霾過程5)為例.

    圖9顯示的700hPa(約海拔3km)高度的風(fēng)場(chǎng)與紫外氣溶膠指數(shù)UVAI的疊加圖,星號(hào)為站點(diǎn)位置.UVAI反映沙塵氣溶膠的程度與分布,觀測(cè)站點(diǎn)處的UVAI與圖5顯示 3km上空退偏比變化情況相同,13日浮塵剛過境站點(diǎn)位置測(cè)得退偏比較小,14與15日再一次沙塵過境情況下,站點(diǎn)上空懸浮著沙塵層,導(dǎo)致退偏比變高.高空風(fēng)場(chǎng)顯示,12~15日中國(guó)北部地區(qū)持續(xù)為西北風(fēng),UVAI分布變化顯示在此期間西部的浮塵一直向東輸送,華北平原上空沙塵污染嚴(yán)重,而漂浮的沙塵顆粒沿途會(huì)發(fā)生圖6顯示的沉降與人為污染混合在一起加重華北平原的霾污染程度.

    3.3 大氣垂直穩(wěn)定性

    分析該站點(diǎn)與溫度和濕度相聯(lián)系的垂直向大氣熱力因素在幾次霾過程演變中的作用,用到了表征大氣熱力不穩(wěn)定度的常用指數(shù)的指數(shù)、近地面溫度露點(diǎn)差(-d)以及假相當(dāng)垂直位溫差(THV).本文用的指數(shù)計(jì)算公式為,=(850-500)+d850-(-d)700[43].式中,和d分別為氣溫和露點(diǎn)溫度,下標(biāo)850,700和500表示和d所在的高度分別為850,700,500hPa.-d代表的是近地面空氣的飽和程度,THV代表的是逆溫的情況.

    由圖11可看出,THV、-d與能見度呈正相關(guān)性,指數(shù)則與能見度具負(fù)相關(guān)性.THV、指數(shù)、-d的逐日變化與能見度相關(guān)系數(shù)分別為0.52、0.56和 0.38.-d越小,空氣中的水汽越接近于飽和,利于水汽凝結(jié),造成霾天氣增強(qiáng).THV、指數(shù)代表了大氣對(duì)流層中低層層結(jié)的穩(wěn)定程度,THV越小、指數(shù)越大,大氣層結(jié)穩(wěn)定性增強(qiáng),霧霾天氣增強(qiáng).

    3.4 氣流軌跡分析

    采用美國(guó)NOAA研制的軌跡模式HYSPLIT4分析觀測(cè)期間該地區(qū)霾天氣的低空污染輸送過程.由圖12可見,清潔天氣團(tuán)來自西北方向,較長(zhǎng)的路徑說明風(fēng)速較大,相對(duì)有利于污染物的擴(kuò)散;有部分污染天的24h氣團(tuán)也來自西北方向,占污染天比例的33%,但軌跡明顯短于清潔天;此外嚴(yán)重污染天氣軌跡來自西、南和東北方向,后向軌跡較短,較小的風(fēng)速有助于污染物的積聚;東北方向的污染氣團(tuán)主要來自河北、山東地區(qū),約占來源比例為14%;近距離傳輸污染氣團(tuán)多來自西面的鄭州和洛陽市區(qū),同時(shí)河南省東面的開封以及南部的城區(qū)也對(duì)該地區(qū)有污染傳輸.

    4 結(jié)論

    4.1 觀測(cè)期間共有5次霾過程,霾天占82%,五六級(jí)重度污染天占污染天數(shù)的46%,PM2.5質(zhì)量濃度達(dá)到環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)(GB3095-2012)一級(jí)(35μg/m3)標(biāo)準(zhǔn)天數(shù)比例為21%,達(dá)到環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)二級(jí)(75μg/m3)的天數(shù)為65%.

    4.2 背景站點(diǎn)近地面污染物以細(xì)粒子污染物為主,約占顆粒物的79%,其中又以氣態(tài)污染物二次生成的細(xì)粒子為主.雖然地基觀測(cè)顯示地面顆粒物以細(xì)粒子為主,光學(xué)觀測(cè)顯示,霾層中上部受浮塵主導(dǎo),粗粒子污染在整層中占據(jù)主導(dǎo).同時(shí)結(jié)合衛(wèi)星數(shù)據(jù)分析,冬季華北平原上層氣溶膠類型主要是沙塵與污染沙塵型氣溶膠,邊界層內(nèi)存在煤煙型和污染大陸型以及污染沙塵氣溶膠組成的混合氣溶膠.

    4.3 該地區(qū)冬季近地面以南向的靜小風(fēng)為主,顆粒物多為本地積聚,另外PM1受東北方向污染源影響,PM1~2.5及PM2.5~10受到西北方向污染源影響.高空風(fēng)場(chǎng)通過遠(yuǎn)距離傳輸水汽和沙塵影響該地區(qū)霾的生成和發(fā)展.冬季大氣對(duì)流層層節(jié)較為穩(wěn)定,假相當(dāng)位溫垂直差、近地面溫度露點(diǎn)差、能見度呈明顯正相關(guān),相關(guān)系數(shù)分別為0.52和0.38,指數(shù)與能見度呈負(fù)相關(guān),相關(guān)系數(shù)為0.56,利于霾的維持.

    4.4 綜合5個(gè)霾過程分析可知,冬季華北平原中部地區(qū)霾污染受到西北浮塵沉降以及南來水汽的影響.對(duì)觀測(cè)站點(diǎn)地區(qū)進(jìn)行后向軌跡分析,該地區(qū)冬季的低空污染傳輸主要來自東北和西北方向,其中東北方向區(qū)域傳輸來自河北與山東,占來源比例的14%,近距離污染傳輸主要來自站點(diǎn)以西的鄭州、洛陽,約占來源比例的33%.

    吳 兌.關(guān)于霾與霧的區(qū)別和灰霾天氣預(yù)警的討論 [J]. 氣象, 2005,31(4):3-7.

    201中華人民共和國(guó)氣象行業(yè)標(biāo)準(zhǔn) [S].

    Ma J, ZWang W, Chen Y, et al. The IPAC-NC field campaign: a pollution and oxidization pool in the lower atmosphere over Huabei, China [J]. Atmospheric Chemistry and Physics, 2012, 12(9):3883-3908.

    Tao M, Chen L, Su L, et al. Satellite observation of regional haze pollution over the North China Plain [J]. Journal of Geophysical Research: Atmospheres (1984~2012), 2012, 117(D12), DOI:10. 1029/2012JD017915.

    Ji D, Wang Y, Wang L, et al. Analysis of heavy pollution episodes in selected cities of northern China [J]. Atmospheric Environment, 2012(50):338-348.

    Andreae M O, Browell E V, Garstang M, et al. Biomass-burning emissions and associated haze layers over Amazonia [J]. Journal of Geophysical Research, 1988,93(D2):1509-1527.

    Okada K, Ikegami M, Zaizen Y, et al. The mixture state of individual aerosol particles in the 1997 Indonesian haze episode [J]. J Aerosol Sci., 2001(32):1269-1279.

    Whiteaker J R, Suess D T, Prather K A. Effects of meteorological conditions on aerosol composition and mixing state in Bakersfield, CA [J]. Environmental Science & Technology, 2002,36(11):2345- 2353.

    吳 兌.近十年中國(guó)灰霾天氣研究綜述 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2012,32(2):257-269.

    Zhang X Y, Wang Y Q, Niu T, et al. Atmospheric aerosol compositions in China: spatial/temporal variability, chemical signature, regional haze distribution and comparisons with global aerosols [J]. Atmospheric Chemistry and Physics, 2012,12(2): 779-799.

    李 鋒,朱 彬,安俊嶺,等.2013年12月初長(zhǎng)江三角洲及周邊地區(qū)重霾污染的數(shù)值模擬 [J]. 中國(guó)環(huán)境科學(xué), 2015,35(7):1965- 1974.

    呂夢(mèng)瑤,劉紅年,張 寧,等.南京市灰霾影響因子的數(shù)值模擬 [J]. 高原氣象, 2011,30(4):929-941.

    Huang K, Zhuang G, Fu J S et al., Shanghai [J]. Atmospheric Chemistry and Physics, 2012,12(1):105-124.

    Jung J, Lee H, Kim Y J, et al. Optical properties of atmospheric aerosols obtained by in situ and remote measurements during 2006 Campaign of Air Quality Research in Beijing (CAREBeijing-2006) [J]. Journal of Geophysical Research: Atmospheres (1984–2012), 2009,114(D2), DOI:10.1029/ 2008JD010337.

    吳 蒙,范紹佳,吳 兌,等.廣州地區(qū)灰霾與清潔天氣變化特征及影響因素分析 [J]. 中國(guó)環(huán)境科學(xué), 2012,32(8):1409-1415.

    王躍思,姚 利,王莉莉,等.2013年元月我國(guó)中東部地區(qū)強(qiáng)霾污染成因分析 [J]. 中國(guó)科學(xué):地球科學(xué), 2014,44(1):15-26.

    廖曉農(nóng),張小玲,王迎春,等.北京地區(qū)冬夏季持續(xù)性霧-霾發(fā)生的環(huán)境氣象條件對(duì)比分析 [J]. 環(huán)境科學(xué), 2014,35(6):2031-2044.

    Tao M, Chen L, Wang Z, et al. A study of urban pollution and haze clouds over northern China during the dusty season based on satellite and surface observations [J]. Atmospheric Environment, 2014,82:183-192.

    H?nel A, Baars H, Althausen D, et al. One-year aerosol profiling with EUCAARI Raman lidar at Shangdianzi GAW station: Beijing plume and seasonal variations [J]. Journal of Geophysical Research: Atmospheres (1984–2012), 2012, 117(D13), DOI:10. 1029/2012JD017577.

    Kalnay E, Kanamitsu M, Kistler R, et al. The NCEP/NCAR 40-year reanalysis project [J]. Bulletin of the American Meteorological Society, 1996,77(3):437-471.

    Levy R C, Remer L A, Kleidman R G, et al. Global evaluation of the Collection 5MODIS dark-target aerosol products over land [J]. Atmospheric Chemistry and Physics, 2010,10(21):10399-10420.

    Torres O, Tanskanen A, Veihelmann B, et al. Aerosols and surface UV products from Ozone Monitoring Instrument observations: An overview [J]. Journal of Geophysical Research: Atmospheres (1984–2012), 2007,112(D24), DOI:10.1029/2007JD008809.

    Ahn C, Torres O, Bhartia P K. Comparison of ozone monitoring instrument UV aerosol products with Aqua/Moderate Resolution Imaging Spectroradiometer and Multiangle Imaging Spectroradiometer observations in 2006 [J]. Journal of Geophysical Research: Atmospheres (1984–2012), 2008, 113(D16),DOI:10.1029/2007JD008832.

    Liu Z, Liu D, Huang J, et al. Airborne dust distributions over the Tibetan Plateau and surrounding areas derived from the first year of CALIPSO lidar observations [J]. Atmospheric Chemistry and Physics, 2008,8(16):5045-5060.

    Mielonen T, Arola A, Komppula M, et al. Comparison of CALIOP level 2aerosol subtypes to aerosol types derived from AERONET inversion data [J]. Geophysical Research Letters, 2009,36(18),DOI:10.1029/2009GL039609.

    Omar A H, Winker D M, Vaughan M A, et al. The CALIPSO automated aerosol classification and lidar ratio selection algorithm [J]. Journal of Atmospheric and Oceanic Technology, 2009,26(10):1994-2014.

    趙 恒,王體健,江 飛,等.利用后向軌跡模式研究TRACE-P期間香港大氣污染物的來源 [J]. 熱帶氣象學(xué)報(bào), 2009,25(2): 181-186.

    吳 兌.再論都市覆與霧的區(qū)別 [J]. 廣東氣象, 2006,32:9-15.

    2012H J.環(huán)境空氣質(zhì)量指數(shù)(AQI)技術(shù)規(guī)定(試行) [S].

    張 劍,劉紅年,唐麗娟,等.蘇州城區(qū)能見度與顆粒物濃度和氣象要素的相關(guān)分析 [J]. 環(huán)境科學(xué)研究, 2011,24(9):982-987.

    鄭小波,羅宇翔,段長(zhǎng)春,等.云貴高原近45年來日照及能見度變化及其成因初步分析 [J]. 高原氣象, 2010,29(4):992-998.

    Nishanth T, Praseed K M, Satheesh Kumar M K, et al. Observational study of surface O3, NO, CH4and total NMHCs at Kannur, India [J]. Aerosol. Air. Qual. Res, 2014(14):1074-1088.

    Han S, Bian H, Feng Y, et al. Analysis of the Relationship between O3, NO and NO2in Tianjin, China [J]. Aerosol Air Qual. Res., 2011(11):128-139.

    Lu Z, Streets D G, de Foy B, et al. Ozone Monitoring Instrument observations of interannual increases in SO2emissions from Indian coal-fired power plants during 2005~2012 [J]. Environmental Science & Technology, 2013,47(24):13993- 14000.

    黃 鶴,蔡子穎,韓素芹,等.天津市PM10, PM2.5和PM1連續(xù)在線觀測(cè)分析 [J]. 環(huán)境科學(xué)研究, 2011,24(8):897-903.

    Xie Y Y, Zhao B, Zhang L, et al. Spatiotemporal variations of PM2.5and PM10concentrations between 31Chinese cities and their relationships with SO2, NO2, CO and O3[J]. Particuology, 2015,20:141-149.

    曹國(guó)良,安心琴,周春紅,等.中國(guó)區(qū)域反應(yīng)性氣體排放源清單 [J]. 中國(guó)環(huán)境科學(xué), 2010,30(7):900-906.

    江 琪,孫業(yè)樂,王自發(fā),等.應(yīng)用顆粒物化學(xué)組分監(jiān)測(cè)儀(ACSM)實(shí)時(shí)在線測(cè)定致霾細(xì)粒子無機(jī)和有機(jī)組分 [J]. 科學(xué)通報(bào), 2013,58(36):3818-3828.

    Althausen D, Engelmann R, Baars H, et al. Portable Raman lidar PollyXT for automated profiling of aerosol backscatter extinction and depolarization [J]. Journal of Atmospheric and Oceanic Technology, 2009,26(11):2366-2378.

    程敘耕,何金海,車慧正,等.1980~2010年中國(guó)區(qū)域地面風(fēng)速對(duì)能見度影響的地理分布特征 [J]. 中國(guó)沙漠, 2013,33(6):1832- 1839.

    朱 焱,楊金彪,朱蓮芳,等.蘇州市能見度與影響因子關(guān)系研究 [J]. 氣象科學(xué), 2011,31(5):626-631.

    章國(guó)材,矯梅燕,李延香.[M]. 北京:氣象出版社, 2007:89-95.

    * 責(zé)任作者, 研究員, sulin@radi.ac.cn

    a rural area of central Northern China Plain based on satellite and surface observations

    WANG Xin-hui1,2, SU Lin1*, TAO Ming-hui1, WANG Zi-feng1, CHEN Liang-fu1, LI Shen-shen1, WANG Yang1,2

    (1.State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, China Academy of Sciences, Beijing 100101, China;2.University of China Academy of Sciences, Beijing 100049, China)., 2016,36(6):1610~1620

    This paper employed ground-based measurements, multiple satellite observations, NCEP reanalysis data, meteorological sounding data and Modeled HYSPLIT4backward trajectory results to ae the pollution characteristics and forming process of haze pollution in Zhongmou, a rural site located in the central North China Plain. The Comprehensive observation campaign lasted from December 13th, 2014 to January 16th, 2015 and a total of five haze pollution episodes were captured Haze days accounted for 82% of the total number of days, while heavily polluted days (Class V and VI based on AQI level) accounted for 46%. The ground-based observations showed that: 1) the concentrations of PM and SO2, NOwere highly correlated, 2) fine particulate matters were dominant among the ground-level pollutants, 3) and the secondary aerosols transformed from gaseous pollutants accounted for a predominant fraction of fine particles. CALIPSO and AMPLE vertical detection data showed that thick haze clouds were dominated by floating dust in the middle and upper layer, and coarse dust particles were prevalent over North China Plain, which accounted for a predominant fraction in the columnar optical volume. Based on meteorological sounding data, the atmospheric stratification of the planetary boundary layer (PBL) was stable, suppressing air turbulence in PBL and favoring formation and maintenance of haze. Meteorological data also showed that near-ground wind speed and direction had a great influence on haze maintenance and dissipation; High-level wind field data indicated that the first haze was influenced by floating dust from the northwest, and the fifth episode was influenced by water vapor from the south. By tracking the 100m air mass using HYSPLIY4model, the low-altitude pollutants were transported from northwest and northeast. The contribution from northeast, accounting for almost 14%, passed Hebei and Shandong province, while short distance transportation traces were mostly from Zhengzhou and Luoyang, accounting for almost 33%.

    haze episode;pollutant character;meteorological factor;regional transport

    X513

    A

    1000-6923(2016)06-1610-11

    王新輝(1990-),女,山東泰安人,碩士研究生,主要從事大氣遙感方向的研究.

    2015-12-02

    中國(guó)科學(xué)院先導(dǎo)專項(xiàng)課題(XDB05020100);國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目(41130528)

    猜你喜歡
    沙塵能見度氣溶膠
    什么是沙塵天氣?
    氣溶膠傳播之謎
    2005—2017年白云機(jī)場(chǎng)能見度變化特征及其與影響因子關(guān)系研究
    氣溶膠中210Po測(cè)定的不確定度評(píng)定
    可怕的沙塵天氣
    四川盆地秋季氣溶膠與云的相關(guān)分析
    低能見度下高速公路主動(dòng)誘導(dǎo)技術(shù)的應(yīng)用
    前向散射能見度儀的常見異?,F(xiàn)象處理及日常維護(hù)
    前向散射能見度儀故障實(shí)例分析
    大氣氣溶膠成核監(jiān)測(cè)
    亚洲国产高清在线一区二区三| 日产精品乱码卡一卡2卡三| 色哟哟哟哟哟哟| 热99在线观看视频| 亚洲七黄色美女视频| 国产黄色视频一区二区在线观看 | 成人性生交大片免费视频hd| 一本精品99久久精品77| 插逼视频在线观看| 日韩中字成人| 丝袜美腿在线中文| 亚洲精品乱码久久久v下载方式| a级一级毛片免费在线观看| 人人妻人人看人人澡| 大又大粗又爽又黄少妇毛片口| 大又大粗又爽又黄少妇毛片口| 日韩大尺度精品在线看网址| 国产在线男女| 日韩欧美在线乱码| 内射极品少妇av片p| 日韩精品中文字幕看吧| 在线观看免费视频日本深夜| 我的老师免费观看完整版| 国产三级中文精品| 国产 一区精品| 九九久久精品国产亚洲av麻豆| 99久久中文字幕三级久久日本| 99在线人妻在线中文字幕| 亚洲av一区综合| 亚洲七黄色美女视频| 午夜福利视频1000在线观看| 在线播放国产精品三级| 黄色一级大片看看| 亚洲成人中文字幕在线播放| 亚洲七黄色美女视频| 国产单亲对白刺激| 欧美激情久久久久久爽电影| 国内揄拍国产精品人妻在线| 中国美女看黄片| 91精品国产九色| 国产av麻豆久久久久久久| 中文亚洲av片在线观看爽| 九九久久精品国产亚洲av麻豆| 丝袜美腿在线中文| 精品99又大又爽又粗少妇毛片| 亚洲性久久影院| 久久欧美精品欧美久久欧美| 欧美性猛交╳xxx乱大交人| 三级毛片av免费| 日韩一区二区视频免费看| 日本成人三级电影网站| 亚洲av中文av极速乱| 五月伊人婷婷丁香| 日韩中字成人| 亚洲美女视频黄频| 男女边吃奶边做爰视频| 日本精品一区二区三区蜜桃| 九色成人免费人妻av| 日韩欧美免费精品| 精品久久久久久久末码| 午夜激情欧美在线| 成年av动漫网址| 啦啦啦韩国在线观看视频| 小说图片视频综合网站| 国产午夜精品论理片| 国产高清不卡午夜福利| 波多野结衣高清无吗| 久久精品国产亚洲av香蕉五月| 亚洲国产欧美人成| 男女视频在线观看网站免费| 99热6这里只有精品| 久久久a久久爽久久v久久| 久久精品国产亚洲av涩爱 | 97在线视频观看| 亚洲成人久久爱视频| 大又大粗又爽又黄少妇毛片口| 免费在线观看影片大全网站| 国产精品爽爽va在线观看网站| 午夜精品在线福利| 成人精品一区二区免费| 国产真实伦视频高清在线观看| 国产男人的电影天堂91| 三级经典国产精品| 久久99热这里只有精品18| 精品国产三级普通话版| 精品午夜福利在线看| 中文字幕av成人在线电影| 久久精品人妻少妇| 欧美日本亚洲视频在线播放| 亚洲成人中文字幕在线播放| 久久精品国产自在天天线| 亚州av有码| 少妇被粗大猛烈的视频| 熟女电影av网| 黄色配什么色好看| 91久久精品国产一区二区成人| 99在线人妻在线中文字幕| 长腿黑丝高跟| 欧美性猛交╳xxx乱大交人| 最后的刺客免费高清国语| 午夜激情欧美在线| 久久精品国产亚洲av香蕉五月| 美女大奶头视频| 在线免费观看的www视频| 狂野欧美激情性xxxx在线观看| 亚洲精品影视一区二区三区av| 嫩草影院入口| 秋霞在线观看毛片| 高清毛片免费看| 亚洲欧美日韩高清在线视频| 看免费成人av毛片| 日韩一本色道免费dvd| 丰满人妻一区二区三区视频av| 男人舔女人下体高潮全视频| 精品久久国产蜜桃| 97热精品久久久久久| 日本黄色片子视频| 国产高清三级在线| 波野结衣二区三区在线| 波多野结衣高清无吗| 亚洲av.av天堂| 色综合站精品国产| 国产精品久久久久久av不卡| 美女cb高潮喷水在线观看| 在线免费观看不下载黄p国产| 国内久久婷婷六月综合欲色啪| 男女那种视频在线观看| 伦理电影大哥的女人| 国产真实伦视频高清在线观看| 亚洲av免费在线观看| 免费人成视频x8x8入口观看| 午夜激情欧美在线| 亚洲人成网站高清观看| 又黄又爽又免费观看的视频| 久久精品人妻少妇| 亚洲自偷自拍三级| 欧美在线一区亚洲| 在线看三级毛片| 大又大粗又爽又黄少妇毛片口| 欧美最新免费一区二区三区| 2021天堂中文幕一二区在线观| 亚洲内射少妇av| 国产亚洲av嫩草精品影院| 蜜桃亚洲精品一区二区三区| 午夜日韩欧美国产| 永久网站在线| 国产成人91sexporn| 日韩一本色道免费dvd| 国产真实伦视频高清在线观看| 高清日韩中文字幕在线| 日本免费a在线| 国产欧美日韩一区二区精品| 九九爱精品视频在线观看| 国产午夜精品论理片| 高清午夜精品一区二区三区 | 国产高清视频在线播放一区| 免费不卡的大黄色大毛片视频在线观看 | 在线a可以看的网站| av中文乱码字幕在线| 日本a在线网址| 观看免费一级毛片| 人人妻人人澡欧美一区二区| 一本久久中文字幕| 久久久久免费精品人妻一区二区| 99热网站在线观看| 在线免费观看不下载黄p国产| 又黄又爽又免费观看的视频| 久久午夜亚洲精品久久| 亚洲人与动物交配视频| 伦精品一区二区三区| 国产 一区精品| 日韩高清综合在线| 欧美潮喷喷水| 欧美日本亚洲视频在线播放| 久久婷婷人人爽人人干人人爱| 人妻制服诱惑在线中文字幕| 不卡视频在线观看欧美| 在线天堂最新版资源| 此物有八面人人有两片| 国产成人91sexporn| av视频在线观看入口| 国产真实乱freesex| 内地一区二区视频在线| 成人欧美大片| 少妇熟女aⅴ在线视频| 国产精品一区二区性色av| av国产免费在线观看| 有码 亚洲区| 男人的好看免费观看在线视频| 国产精品一区二区免费欧美| 看十八女毛片水多多多| 欧美一区二区精品小视频在线| 少妇被粗大猛烈的视频| 天堂动漫精品| 国产亚洲91精品色在线| 国产精品爽爽va在线观看网站| 在线观看午夜福利视频| 麻豆乱淫一区二区| 亚洲人成网站高清观看| 久久久久国产网址| 精品人妻偷拍中文字幕| 亚洲18禁久久av| 在线观看美女被高潮喷水网站| 久久国内精品自在自线图片| 亚洲精品456在线播放app| 亚洲丝袜综合中文字幕| 日韩,欧美,国产一区二区三区 | 一级av片app| 国产精品久久视频播放| 免费黄网站久久成人精品| 网址你懂的国产日韩在线| 国产高清激情床上av| 淫秽高清视频在线观看| 久久天躁狠狠躁夜夜2o2o| 校园春色视频在线观看| 看免费成人av毛片| 国产高清激情床上av| 久久韩国三级中文字幕| 欧美成人精品欧美一级黄| 99久久精品一区二区三区| 国产高潮美女av| 久久久久久久午夜电影| 欧美xxxx性猛交bbbb| 亚洲成av人片在线播放无| 免费观看在线日韩| 欧美日韩精品成人综合77777| 99国产精品一区二区蜜桃av| 免费观看精品视频网站| 色综合亚洲欧美另类图片| 亚洲最大成人手机在线| 搡老熟女国产l中国老女人| 免费不卡的大黄色大毛片视频在线观看 | 久久午夜亚洲精品久久| 搡老熟女国产l中国老女人| 最新中文字幕久久久久| 亚洲国产精品成人久久小说 | av福利片在线观看| а√天堂www在线а√下载| 日本在线视频免费播放| 久久精品人妻少妇| 亚洲精品影视一区二区三区av| 国产综合懂色| 日本五十路高清| av在线播放精品| 尤物成人国产欧美一区二区三区| 不卡一级毛片| 久久久久久国产a免费观看| 免费在线观看成人毛片| 精品福利观看| 日韩欧美一区二区三区在线观看| 亚洲精品456在线播放app| 黄色日韩在线| 丰满乱子伦码专区| 国产av一区在线观看免费| 日产精品乱码卡一卡2卡三| 丝袜美腿在线中文| 日韩大尺度精品在线看网址| 美女免费视频网站| 插阴视频在线观看视频| 亚洲人成网站在线播| 高清午夜精品一区二区三区 | 免费观看在线日韩| 亚洲最大成人av| 99热这里只有是精品在线观看| 晚上一个人看的免费电影| 国产aⅴ精品一区二区三区波| 观看免费一级毛片| 国内精品一区二区在线观看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲,欧美,日韩| 国产高清不卡午夜福利| 91午夜精品亚洲一区二区三区| 亚洲国产欧洲综合997久久,| 毛片女人毛片| 青春草视频在线免费观看| 在线观看av片永久免费下载| 免费大片18禁| 亚洲欧美日韩高清在线视频| 婷婷精品国产亚洲av在线| 日日干狠狠操夜夜爽| 成熟少妇高潮喷水视频| 日韩人妻高清精品专区| 99久久精品国产国产毛片| 免费观看精品视频网站| 少妇被粗大猛烈的视频| 别揉我奶头~嗯~啊~动态视频| 亚洲专区国产一区二区| 久久久久久国产a免费观看| 老司机福利观看| 中文亚洲av片在线观看爽| 九九热线精品视视频播放| 最近最新中文字幕大全电影3| 欧美+亚洲+日韩+国产| 69人妻影院| 老师上课跳d突然被开到最大视频| 国产高潮美女av| 日本a在线网址| 蜜桃久久精品国产亚洲av| 午夜福利18| 一边摸一边抽搐一进一小说| 成人永久免费在线观看视频| 色吧在线观看| 一卡2卡三卡四卡精品乱码亚洲| 不卡一级毛片| 观看美女的网站| 观看免费一级毛片| 国产真实伦视频高清在线观看| 狂野欧美白嫩少妇大欣赏| 三级经典国产精品| 激情 狠狠 欧美| 亚洲av不卡在线观看| 亚洲欧美成人精品一区二区| 一区二区三区高清视频在线| 国产熟女欧美一区二区| 日本黄大片高清| 国产蜜桃级精品一区二区三区| 亚洲av熟女| 国产久久久一区二区三区| 国产精品乱码一区二三区的特点| 午夜福利18| 超碰av人人做人人爽久久| 国产精品一区www在线观看| 最近最新中文字幕大全电影3| 九九在线视频观看精品| 久久久成人免费电影| 插阴视频在线观看视频| 国产 一区 欧美 日韩| 欧美日韩综合久久久久久| 国产精品一及| 毛片女人毛片| 啦啦啦韩国在线观看视频| 在线a可以看的网站| 搞女人的毛片| 亚洲av成人精品一区久久| 欧美一级a爱片免费观看看| 色综合色国产| 国产一级毛片七仙女欲春2| 久久久久久久久久黄片| 床上黄色一级片| 最后的刺客免费高清国语| 国产av一区在线观看免费| 亚洲av二区三区四区| 日韩精品有码人妻一区| 国产成人a∨麻豆精品| 精品午夜福利视频在线观看一区| 精品一区二区三区av网在线观看| 欧美区成人在线视频| 久久亚洲精品不卡| 午夜福利在线观看吧| 又爽又黄a免费视频| av.在线天堂| 国产精品国产三级国产av玫瑰| 一本一本综合久久| 欧美国产日韩亚洲一区| 亚洲av二区三区四区| 亚洲精品一区av在线观看| 不卡视频在线观看欧美| 免费在线观看影片大全网站| 一级毛片久久久久久久久女| 麻豆av噜噜一区二区三区| 真人做人爱边吃奶动态| 国产乱人视频| a级毛片a级免费在线| a级毛色黄片| 真实男女啪啪啪动态图| 在线观看66精品国产| 国产一区二区在线av高清观看| 国产精品美女特级片免费视频播放器| av专区在线播放| av中文乱码字幕在线| 免费大片18禁| 免费在线观看成人毛片| 国产精品人妻久久久久久| 亚洲人成网站在线播放欧美日韩| 亚洲av.av天堂| 精品人妻熟女av久视频| 校园春色视频在线观看| 国产高清不卡午夜福利| 99热网站在线观看| 99久国产av精品国产电影| 美女被艹到高潮喷水动态| 成人av一区二区三区在线看| 欧美一区二区国产精品久久精品| 在线看三级毛片| 白带黄色成豆腐渣| a级毛片a级免费在线| 老司机福利观看| 精品一区二区三区视频在线观看免费| 香蕉av资源在线| 18禁在线无遮挡免费观看视频 | 久久人人精品亚洲av| 欧美性感艳星| 国产91av在线免费观看| 在线天堂最新版资源| 国产黄a三级三级三级人| 日本撒尿小便嘘嘘汇集6| 亚洲欧美精品自产自拍| 成人午夜高清在线视频| 欧美高清成人免费视频www| 国产精品免费一区二区三区在线| 午夜精品在线福利| 亚洲,欧美,日韩| 久久亚洲精品不卡| 亚洲欧美清纯卡通| 国产精品野战在线观看| 午夜影院日韩av| 欧美在线一区亚洲| 久久精品国产鲁丝片午夜精品| 亚洲中文字幕一区二区三区有码在线看| 亚洲av免费在线观看| 性欧美人与动物交配| 免费av不卡在线播放| 国产精品国产三级国产av玫瑰| 可以在线观看的亚洲视频| 尾随美女入室| 在线观看一区二区三区| 天天躁日日操中文字幕| 观看免费一级毛片| av.在线天堂| 国产一区二区亚洲精品在线观看| 夜夜夜夜夜久久久久| 老熟妇乱子伦视频在线观看| 在线免费观看的www视频| 国产中年淑女户外野战色| 国产精品乱码一区二三区的特点| 菩萨蛮人人尽说江南好唐韦庄 | 国产人妻一区二区三区在| 午夜免费男女啪啪视频观看 | 亚洲欧美日韩高清专用| 午夜免费男女啪啪视频观看 | 国产高清有码在线观看视频| 午夜久久久久精精品| 日本黄大片高清| 中国国产av一级| 午夜久久久久精精品| 久久精品人妻少妇| а√天堂www在线а√下载| 免费在线观看成人毛片| h日本视频在线播放| 亚洲av.av天堂| 日本欧美国产在线视频| 亚洲天堂国产精品一区在线| 身体一侧抽搐| 在线免费观看的www视频| 午夜免费激情av| 啦啦啦观看免费观看视频高清| 91久久精品电影网| 少妇人妻一区二区三区视频| 18禁在线播放成人免费| 高清毛片免费看| 久久热精品热| 夜夜爽天天搞| 国产av麻豆久久久久久久| 午夜福利在线在线| 校园人妻丝袜中文字幕| 一进一出抽搐gif免费好疼| 亚洲精品一区av在线观看| 成人无遮挡网站| 国产精品电影一区二区三区| 日本五十路高清| 校园人妻丝袜中文字幕| 国产色爽女视频免费观看| av在线观看视频网站免费| 性插视频无遮挡在线免费观看| 亚洲不卡免费看| av天堂中文字幕网| 精品午夜福利在线看| 国产黄色视频一区二区在线观看 | 午夜爱爱视频在线播放| 国产极品精品免费视频能看的| 欧美在线一区亚洲| 国产探花极品一区二区| 少妇裸体淫交视频免费看高清| 搡老熟女国产l中国老女人| 亚洲最大成人av| 精品国内亚洲2022精品成人| 国产精品精品国产色婷婷| 看黄色毛片网站| 欧美性猛交黑人性爽| 久久久欧美国产精品| 一级毛片久久久久久久久女| 欧美日韩国产亚洲二区| 亚洲精品影视一区二区三区av| 麻豆久久精品国产亚洲av| 国产高清有码在线观看视频| 国产精品不卡视频一区二区| 综合色av麻豆| 亚洲,欧美,日韩| av卡一久久| 欧美+日韩+精品| 观看美女的网站| 熟女电影av网| 亚洲国产色片| 中文字幕av成人在线电影| 国产人妻一区二区三区在| 精品人妻视频免费看| 久久99热6这里只有精品| 免费看a级黄色片| 国产视频内射| 99热网站在线观看| 亚洲在线自拍视频| 哪里可以看免费的av片| 又黄又爽又刺激的免费视频.| 噜噜噜噜噜久久久久久91| 日韩欧美在线乱码| 亚洲av免费在线观看| 亚洲一区高清亚洲精品| 色综合色国产| 国产亚洲91精品色在线| 晚上一个人看的免费电影| 亚洲专区国产一区二区| 日韩成人av中文字幕在线观看 | 免费看a级黄色片| 亚洲欧美日韩高清专用| 少妇熟女aⅴ在线视频| 国产精品日韩av在线免费观看| 欧美日韩在线观看h| 欧美在线一区亚洲| av专区在线播放| 亚洲精品成人久久久久久| 亚洲成人久久爱视频| 女同久久另类99精品国产91| 欧美日韩乱码在线| av在线老鸭窝| 看非洲黑人一级黄片| 在线免费观看不下载黄p国产| 变态另类成人亚洲欧美熟女| 一级毛片久久久久久久久女| 国产精品一区二区免费欧美| 深夜a级毛片| 精品人妻熟女av久视频| 女同久久另类99精品国产91| 日韩成人av中文字幕在线观看 | 欧美又色又爽又黄视频| 长腿黑丝高跟| 国产真实乱freesex| 少妇熟女aⅴ在线视频| 国产成人91sexporn| 免费看av在线观看网站| 国产免费一级a男人的天堂| 色播亚洲综合网| 可以在线观看毛片的网站| 一区福利在线观看| 日本黄色片子视频| 国产爱豆传媒在线观看| 18禁裸乳无遮挡免费网站照片| 成人高潮视频无遮挡免费网站| 国产黄a三级三级三级人| 九九爱精品视频在线观看| 国产乱人视频| 欧美成人一区二区免费高清观看| 亚洲无线观看免费| 久久精品国产亚洲av香蕉五月| 97热精品久久久久久| 18+在线观看网站| 97超视频在线观看视频| 亚洲不卡免费看| 亚洲av中文av极速乱| 欧美色视频一区免费| 亚洲经典国产精华液单| 九九热线精品视视频播放| 尤物成人国产欧美一区二区三区| 日韩精品青青久久久久久| 一级毛片我不卡| 国语自产精品视频在线第100页| 九九热线精品视视频播放| 老师上课跳d突然被开到最大视频| 深夜a级毛片| 一级黄片播放器| 寂寞人妻少妇视频99o| 丝袜美腿在线中文| 搞女人的毛片| 噜噜噜噜噜久久久久久91| 国产精品电影一区二区三区| 黑人高潮一二区| 色尼玛亚洲综合影院| 成人三级黄色视频| 国产美女午夜福利| 成人高潮视频无遮挡免费网站| 亚洲精品国产成人久久av| 观看免费一级毛片| 美女免费视频网站| 欧美色欧美亚洲另类二区| 又粗又爽又猛毛片免费看| 男人和女人高潮做爰伦理| 亚洲第一区二区三区不卡| 插阴视频在线观看视频| 又黄又爽又免费观看的视频| 听说在线观看完整版免费高清| 成人精品一区二区免费| 少妇的逼好多水| 国产高清三级在线| 国产精品av视频在线免费观看| 校园春色视频在线观看| 搡老岳熟女国产| 乱系列少妇在线播放| 国产精品亚洲一级av第二区| 欧美激情在线99| 亚洲国产精品合色在线| 日韩成人av中文字幕在线观看 | 日日啪夜夜撸| 久久久久国内视频| 欧美激情国产日韩精品一区| 欧美性猛交╳xxx乱大交人| 看黄色毛片网站| 欧美性猛交黑人性爽| 99热精品在线国产| 国产精品一区二区三区四区久久| 中文在线观看免费www的网站| 午夜a级毛片| 午夜精品国产一区二区电影 | 最近中文字幕高清免费大全6| 久久精品国产亚洲av涩爱 | 久久热精品热| 最近在线观看免费完整版|