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

    天空輻射計(jì)觀測(cè)反演北京城區(qū)氣溶膠光學(xué)特性

    2020-05-13 08:39:08楊先逸車慧正陳權(quán)亮梁苑新
    應(yīng)用氣象學(xué)報(bào) 2020年3期
    關(guān)鍵詞:氣溶膠城區(qū)光學(xué)

    楊先逸 車慧正 陳權(quán)亮 梁苑新

    1)(成都信息工程大學(xué)大氣科學(xué)學(xué)院, 高原大氣與環(huán)境四川省重點(diǎn)實(shí)驗(yàn)室, 成都 610225)2)(中國(guó)氣象科學(xué)研究院災(zāi)害天氣國(guó)家重點(diǎn)實(shí)驗(yàn)室, 中國(guó)氣象局大氣化學(xué)重點(diǎn)開(kāi)放實(shí)驗(yàn)室, 北京 100081)

    引 言

    大氣氣溶膠可以通過(guò)散射和吸收太陽(yáng)輻射,影響地球的輻射收支,進(jìn)而對(duì)氣候變化產(chǎn)生影響[1]。此外氣溶膠粒子還可以作為云凝結(jié)核或冰核,不僅對(duì)云形成產(chǎn)生重要影響,也可以通過(guò)改變?cè)平Y(jié)構(gòu)和光學(xué)特性進(jìn)而影響云滴的大小和云的生命周期,間接影響氣候變化[2-4]。研究表明:大氣氣溶膠不僅影響全球和區(qū)域氣候變化,而且還會(huì)造成環(huán)境污染,對(duì)人類健康產(chǎn)生不良影響。盡管各國(guó)科學(xué)家針對(duì)氣溶膠進(jìn)行了多年研究,但氣溶膠的濃度變化和光學(xué)特性以及顯著的時(shí)空變化仍是當(dāng)前全球氣候變化評(píng)估和預(yù)測(cè)中不確定性的較大來(lái)源之一[5]。

    地基觀測(cè)是一種較為有用且準(zhǔn)確的測(cè)量方法,雖然地面觀測(cè)手段在獲取氣溶膠空間分布方面不占優(yōu)勢(shì),但其持續(xù)穩(wěn)定觀測(cè)獲得的原始數(shù)據(jù)和相對(duì)成熟的反演方法使其在氣溶膠觀測(cè)研究中占有重要地位[6]。AERONET (Aerosol Robotic Network) 是美國(guó)NASA (National Aeronautics and Space Administration) 在全球范圍內(nèi)建立的太陽(yáng)光度計(jì)觀測(cè)網(wǎng)絡(luò),采用法國(guó)CIMEL公司生產(chǎn)的CE318型太陽(yáng)光度計(jì),基于太陽(yáng)直接輻射和天空散射輻射反演氣溶膠光學(xué)特性參數(shù)[7-8]。CARSNET (China Aerosol Remote Sensing Network) 是中國(guó)氣象局建立的中國(guó)氣溶膠地基遙感監(jiān)測(cè)網(wǎng),最初用于監(jiān)測(cè)沙塵暴,經(jīng)逐步完善后用于分析中國(guó)區(qū)域大氣氣溶膠光學(xué)特性的分布和變化特征[9-10]。PHOTONS (Pho- tométrie pour le Traitement Opérationnel de Normalisation Satellitaire) 是法國(guó)在歐洲和非洲地區(qū)建立的氣溶膠光學(xué)特性觀測(cè)網(wǎng)[11]。SKYNET(Skyrad Network)[12]采用日本PREDE公司生產(chǎn)的POM型天空輻射計(jì),通過(guò)測(cè)量可見(jiàn)光到近紅外波段的太陽(yáng)直接輻射和天空散射輻射反演一系列氣溶膠特征參數(shù)[13-14], 該網(wǎng)絡(luò)主要位于亞洲和歐洲,致力于氣溶膠-云-輻射相互作用研究。歐洲地區(qū)也同樣建立了ESR(European Skynet Radiometer Net- work) 天空輻射計(jì)觀測(cè)站網(wǎng)用以研究歐洲地區(qū)氣溶膠光學(xué)特性[15-16]。

    隨著經(jīng)濟(jì)的飛速發(fā)展,工業(yè)化和城市化進(jìn)程加速,我國(guó)已成為人為氣溶膠的主要源地之一,特別是在人口密集、經(jīng)濟(jì)發(fā)達(dá)的京津冀地區(qū)[17-19]。大氣氣溶膠光學(xué)特性除了受人口密度和人為因素影響外,還與沙塵、霧霾等特殊天氣條件密切相關(guān)[20]。車慧正等[21]發(fā)現(xiàn)北京城區(qū)氣溶膠粒子體積譜分布具有明顯城市氣溶膠譜分布特征。孫文文等[22]指出北京城區(qū)氣溶膠主要受到內(nèi)蒙古地區(qū)的沙塵及近海的海鹽顆粒物遠(yuǎn)距離輸送的影響,局地排放的氣溶膠顆粒物主要來(lái)自機(jī)動(dòng)車尾氣和工業(yè)排放。楊溯等[23]通過(guò)分析北京城區(qū)氣溶膠光學(xué)特性變化指出沙塵氣溶膠是該地區(qū)春季重要的氣溶膠類型。夏祥鰲等[6]指出除沙塵事件外,華北地區(qū)人為氣溶膠排放也是導(dǎo)致該地區(qū)大氣氣溶膠濃度偏大的重要原因。王玲等[24]利用CE318型太陽(yáng)光度計(jì)觀測(cè)數(shù)據(jù)反演北京城區(qū)大氣氣溶膠成分,發(fā)現(xiàn)霾和清潔天氣條件下氣溶膠成分差異顯著。顏鵬等[25]發(fā)現(xiàn)在京津地區(qū)氣溶膠光學(xué)特性受天氣過(guò)程影響很大。

    本文利用2018年10月11日—2019年9月30日POM-02型天空輻射計(jì)觀測(cè)數(shù)據(jù),反演北京城區(qū)氣溶膠光學(xué)特性,對(duì)了解實(shí)施污染防治后的北京城區(qū)氣溶膠光學(xué)特性變化特征具有重要意義。利用天空輻射計(jì)反演結(jié)果,結(jié)合氣象觀測(cè)數(shù)據(jù)以及激光雷達(dá)數(shù)據(jù)多角度分析2019年1月北京城區(qū)一次典型污染事件,為進(jìn)一步認(rèn)識(shí)北京城區(qū)氣溶膠光學(xué)特性提供參考。

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

    1.1 觀測(cè)簡(jiǎn)介

    本文選取2018年10月11日—2019年9月30日為觀測(cè)時(shí)段(共355 d),觀測(cè)儀器為日本PREDE公司生產(chǎn)的POM-02天空輻射計(jì),架設(shè)在北京海淀區(qū)中國(guó)氣象科學(xué)研究院樓頂(39.933°N,116.317°E,105 m), 該站點(diǎn)位于北京西北三環(huán)以內(nèi),為典型城市站點(diǎn)。天空輻射計(jì)可以自動(dòng)觀測(cè)可見(jiàn)光和近紅外波段的太陽(yáng)直接輻射和散射輻射,可視化半角度為1°。根據(jù)Nakajima等[13]反演方法,5個(gè)通道(400,500,670,870 nm和1020 nm)反演得到的氣溶膠光學(xué)特性包括氣溶膠光學(xué)厚度τ、單次散射反照率ω、粒子譜分布、?ngstr?m波長(zhǎng)指數(shù)α。觀測(cè)期間各月觀測(cè)樣本容量能夠滿足統(tǒng)計(jì)學(xué)分析要求。我國(guó)的PM2.5標(biāo)準(zhǔn)值為24 h平均濃度小于75 μg·m-3,依據(jù)此標(biāo)準(zhǔn)將觀測(cè)期間日平均PM2.5濃度大于75 μg·m-3作為污染天,分別統(tǒng)計(jì)了儀器各月有效觀測(cè)日數(shù)和有效觀測(cè)次數(shù)(有效觀測(cè)次數(shù)為經(jīng)過(guò)質(zhì)量控制后反演得到的數(shù)據(jù)量), 以及有效觀測(cè)中污染天數(shù)量和污染天有效觀測(cè)次數(shù)。

    1.2 數(shù) 據(jù)

    中分辨率成像光譜儀(moderate-resolution imaging spectroradio-meter,MODIS)數(shù)據(jù)來(lái)自美國(guó)國(guó)家航空航天局(https:∥ladsweb.modaps.eosdis.nasa.gov/search/)。逐小時(shí)顆粒物質(zhì)量濃度(PM2.5和PM10)數(shù)據(jù)來(lái)自北京市環(huán)境保護(hù)監(jiān)測(cè)中心(http:∥www.bjmemc.com.cn/)官園站。地面氣象觀測(cè)數(shù)據(jù)來(lái)自國(guó)家氣象信息中心,包括風(fēng)向、風(fēng)速、溫度、降水量、水平能見(jiàn)度和相對(duì)濕度等要素。以上數(shù)據(jù)時(shí)段均為2018年10月—2019年9月。

    基于拉曼-米激光雷達(dá)觀測(cè)得到氣溶膠消光系數(shù)和退偏振比數(shù)據(jù),該儀器位于中國(guó)氣象局院內(nèi),用于分析2019年1月9日00:00(北京時(shí),下同)—14日23:00北京城區(qū)氣溶膠垂直分布變化。儀器存在160 m的探測(cè)盲區(qū),垂直分辨率為7.50 m,觀測(cè)周期為5 min。

    1.3 數(shù)據(jù)可靠性

    目前獲取氣溶膠光學(xué)厚度數(shù)據(jù)主要途徑包括地基觀測(cè)和衛(wèi)星遙感手段,其中衛(wèi)星遙感具有覆蓋面廣、信息獲取方便等特點(diǎn)。毛節(jié)泰等[26]對(duì)比了2001年北京城區(qū)太陽(yáng)光度計(jì)觀測(cè)結(jié)果和MODIS光學(xué)產(chǎn)品,證實(shí)了MODIS反演氣溶膠光學(xué)特征的可靠性。許多學(xué)者利用MODIS衛(wèi)星反演的氣溶膠光學(xué)厚度進(jìn)行研究[27-30]。本文選用Aqua 10 km 氣溶膠光學(xué)厚度產(chǎn)品與天空輻射計(jì)觀測(cè)得到的氣溶膠光學(xué)厚度(τs)進(jìn)行對(duì)比,期望誤差采用±(0.05+20%τs)表示。由圖1可以看到,衛(wèi)星反演結(jié)果和地基觀測(cè)結(jié)果相比普遍偏高,約有52.94%的匹配點(diǎn)落入期望誤差范圍以內(nèi),兩者相關(guān)系數(shù)為0.90(達(dá)到0.01顯著性水平),均方根誤差為0.17,具有較好的一致性。衛(wèi)星氣溶膠光學(xué)厚度反演算法與植被覆蓋密切相關(guān),而城市地區(qū)植被相對(duì)較少,導(dǎo)致反演結(jié)果出現(xiàn)偏差。

    圖1 北京城區(qū)(39.933°N,116.317°E,105 m)天空輻射計(jì)觀測(cè)氣溶膠光學(xué)厚度與MODIS衛(wèi)星反演氣溶膠光學(xué)厚度對(duì)比(實(shí)線和虛線分別表示1:1線和期望誤差分界線)Fig.1 Comparison of the aerosol optical depth observed by skyradiometer and those of MODIS in Beijing(39.933°N,116.317°E,105 m)(solid and dashed lines denote one-one line and the expected error line,respectively)

    2 氣溶膠光學(xué)特性

    2.1 基于所有觀測(cè)結(jié)果的氣溶膠光學(xué)特性逐月變化

    由圖2a可以看到,北京城區(qū)τ在不同季節(jié)存在一定差異。τ在春季、夏季、秋季、冬季平均值分別為0.55,0.58,0.47和0.49。春夏季τ較大,秋冬季τ較小,該結(jié)果與車慧正等[21]、楊溯等[23]研究結(jié)論一致。τ月平均值2—7月較大,其中全年最高值出現(xiàn)在6月,為0.71,其他月份τ相對(duì)較小,不超過(guò)0.50。本文利用天空輻射計(jì)得出北京城區(qū)500 nmω年平均值為0.93,與車慧正等[21]利用天空輻射計(jì)獲得的ω年平均值0.92相差0.01。同時(shí)略高于高中明等[31]利用太陽(yáng)光度計(jì)得到北京城區(qū)440 nm和675 nmω年平均值(分別為0.90和0.92)。由圖2b可以看到,北京城區(qū)ω平均值8月最高,為0.96;5月最低,為0.89。結(jié)果表明:北京城區(qū)夏季氣溶膠粒子散射能力較強(qiáng),春季較弱。夏季ω平均值為0.95,春季為0.90。春季空氣中具有吸收性的沙塵、浮塵粒子較多,會(huì)導(dǎo)致ω值偏低[6,32]。由圖2c可以看到,α整體上在夏季(1.11)、秋季(1.09)較高,春季(0.89)、冬季(0.95)較低,說(shuō)明春季和冬季氣溶膠粒子粒徑較大,粗粒子占總粒子比例較高,夏季氣溶膠粒子粒徑較小,細(xì)粒子占總粒子比例較高。

    圖2 2018年10月—2019年9月北京城區(qū)氣溶膠光學(xué)特性的逐月變化(a)500 nm氣溶膠光學(xué)厚度,(b)500 nm單次散射反照率,(c)440~870 nm ?ngstr?m波長(zhǎng)指數(shù)Fig.2 Monthly averaged variation in aerosol optical properties from skyradiometer measurements in Beijing from Oct 2018 to Sep 2019(a)aerosol optical depth at 500 nm, (b)single scattering albedo at 500 nm, (c)?ngstr?mm exponent at 440-870 nm

    由圖3可以看到,北京城區(qū)氣溶膠粒子體積譜分布在每個(gè)月都呈典型的雙峰型分布特征,符合典型的城市氣溶膠分布特征。細(xì)模態(tài)粒子粒徑主要集中在0.17 μm左右,粗模態(tài)粒子粒徑主要集中在3.62~5.29 μm。整體看,不同月份的氣溶膠粒子體積譜分布存在明顯差異。4月粗模態(tài)粒子濃度顯著高于其他月份,這可能與沙塵事件影響有關(guān) (統(tǒng)計(jì)發(fā)現(xiàn)2019年4月4—5日有沙塵天氣,其他月份沒(méi)有記錄)。細(xì)模態(tài)粒子濃度在6月最高,說(shuō)明氣溶膠細(xì)粒子數(shù)目明顯增多。細(xì)模態(tài)粒子粒徑在7月主要集中在0.25 μm,說(shuō)明夏季高溫高濕條件容易導(dǎo)致氣溶膠粒子吸濕增長(zhǎng)。8月細(xì)模態(tài)和粗模態(tài)氣溶膠粒子濃度均較小,這可能與降水過(guò)程帶來(lái)的濕清除效應(yīng)導(dǎo)致大氣中氣溶膠粒子減小有關(guān)。

    為進(jìn)一步探討不同類型氣溶膠粒子的體積譜分布特征,結(jié)合ω和α對(duì)氣溶膠粒子進(jìn)行分類:Ⅰ類,細(xì)模態(tài)強(qiáng)吸收性(α>1.2,ω≤0.85);Ⅱ類,混合型強(qiáng)吸收性(0.6<α≤1.2,ω≤0.85);Ⅲ類,粗模態(tài)強(qiáng)吸收性(α≤0.6,ω≤0.85);Ⅳ類,細(xì)模態(tài)中度吸收性(α>1.2,0.85<ω≤0.95);Ⅴ類,混合型中度吸收性(0.6<α≤1.2,0.85<ω≤0.95);Ⅵ類,粗模態(tài)中度吸收性(α≤0.6,0.85<ω≤0.95);Ⅶ類,細(xì)模態(tài)弱吸收性(α>1.2,ω>0.95);Ⅷ類,混合型弱吸收性(0.6<α≤1.2,ω>0.95);Ⅸ類,粗模態(tài)弱吸收性(α≤0.6,ω>0.95)。結(jié)合圖4a、圖4b 和圖4c,可以看到隨著粗粒子所占比例的增加,粗模態(tài)粒子濃度也明顯增大,對(duì)應(yīng)細(xì)模態(tài)粒子濃度減小。同時(shí),隨著氣溶膠粒子吸收性減弱,粗、細(xì)模態(tài)粒子濃度的變化幅度不大。該結(jié)果指出粗、細(xì)粒子比例會(huì)影響粒子體積譜分布,而粒子的吸收散射能力與其體積譜分布沒(méi)有必然聯(lián)系。北京城區(qū)細(xì)粒子占主導(dǎo)的粒子體積譜分布呈典型的雙峰型,粗、細(xì)模態(tài)粒子濃度大小相當(dāng);混合型粒子體積譜分布呈雙峰型,粗模態(tài)粒子濃度大于細(xì)模態(tài)粒子濃度;粗粒子占主導(dǎo)的粒子體積譜分布呈單峰型,粗模態(tài)粒子粒徑主要集中在3.62 μm。

    圖3 2018年10月—2019年9月北京城區(qū)粒子譜分布逐月變化Fig.3 Monthly averaged variation in volume size distribution from skyradiometer measurements in Beijing from Oct 2018 to Sep 2019

    圖4 2018年10月—2019年9月北京城區(qū)不同氣溶膠類型的粒子譜分布(a)強(qiáng)吸收性粒子,(b)中度吸收性粒子,(c)弱吸收性粒子Fig.4 Volume size distribution of different aerosol types in Beijing from Oct 2018 to Sep 2019 (a)highly absorbing particles,(b)moderately absorbing particles,(c)weakly absorbing particles

    2.2 污染與清潔天氣條件下氣溶膠光學(xué)特性對(duì)比

    根據(jù)國(guó)家PM2.5污染分級(jí)標(biāo)準(zhǔn),統(tǒng)計(jì)2018年10月11日—2019年9月30日各污染等級(jí)的實(shí)際日數(shù)和天空輻射計(jì)的有效觀測(cè)日數(shù)。如果出現(xiàn)云層較厚、降水、強(qiáng)沙塵暴等天氣現(xiàn)象,無(wú)法滿足天空輻射計(jì)的正常工作要求,導(dǎo)致缺測(cè)。由表1可以看到,基于PM2.5污染分級(jí)標(biāo)準(zhǔn),對(duì)比不同PM2.5污染等級(jí)下實(shí)際日數(shù)與儀器觀測(cè)日數(shù)的數(shù)目,發(fā)現(xiàn)天空輻射計(jì)在輕度、中度污染等級(jí)具有較好的觀測(cè)能力。2018年10月—2019年9月北京城區(qū)空氣質(zhì)量總體表現(xiàn)良好,實(shí)際日數(shù)中空氣質(zhì)量為優(yōu)的日數(shù)占總?cè)諗?shù)的48%,非污染日數(shù)比例超過(guò)83%,污染級(jí)別主要集中于輕度污染等級(jí),輕度污染日數(shù)占總污染日數(shù)的62%,但距離世界衛(wèi)生組織PM2.5標(biāo)準(zhǔn)值還存在差距。

    表1 2018年10月11日—2019年9月30日空氣質(zhì)量等級(jí)統(tǒng)計(jì)Table 1 Air quality rank statistics from 11 Oct 2018 to 30 Sep 2019

    在污染天氣和清潔天氣條件下,2018年10月11日—2019年9月30日τ平均值分別為0.85和0.49,即污染天氣條件下氣溶膠粒子對(duì)光的削減能力強(qiáng);ω平均值分別為0.96和0.92,即污染天氣條件下氣溶膠粒子的散射能力較強(qiáng),污染天氣條件下ω月平均最高值出現(xiàn)在1月,清潔天氣條件下ω月平均最低值出現(xiàn)在3月;α平均值分別為0.97和1.01,冬季污染天氣和清潔天氣條件下,α分別為1.02 和0.91,春季污染天氣和清潔天氣條件下,α分別為0.87和0.90。上述分析表明:在污染天氣條件下,冬季北京城區(qū)以細(xì)粒子為主,春季以粗粒子為主,其原因主要是該地區(qū)冬季污染天氣類型以霾為主,而春季主要為沙塵所導(dǎo)致;PM2.5濃度平均值分別為107.22 μg·m-3和47.16 μg·m-3,污染天氣條件下PM2.5濃度月平均最高值出現(xiàn)在1月,為146.53 μg·m-3。

    2.3 2019年1月北京重污染過(guò)程氣溶膠光學(xué)特性

    根據(jù)2.2節(jié)分析可知,北京城區(qū)冬季污染天氣條件下PM2.5濃度最大,污染最為嚴(yán)重,很多學(xué)者也針對(duì)北京城區(qū)冬季污染事件進(jìn)行了研究[33-34]。本文將研究2019年1月北京城區(qū)冬季大氣污染狀況。由圖5可以看到,大氣污染狀況越嚴(yán)重,τ越大,PM2.5在空氣中含量越高。有效觀測(cè)日的τ日平均值與PM2.5濃度日變化趨勢(shì)大體一致。根據(jù)國(guó)家PM2.5超標(biāo)分級(jí)標(biāo)準(zhǔn)發(fā)現(xiàn)2019年1月10—14日北京城區(qū)PM2.5濃度持續(xù)處于高值期,最高日平均值達(dá)227.00 μg·m-3,τ均高于0.40,選取2019年1月8—15日(圖5陰影部分)的一次完整污染事件進(jìn)行分析。

    圖6a為2019年1月8—15日北京大氣顆粒物濃度逐小時(shí)變化,PM2.5和PM10濃度變化趨勢(shì)基本相同。在污染天氣條件下,PM2.5和PM10濃度大小較為接近,說(shuō)明氣溶膠粒子中PM2.5占絕大部分,細(xì)粒子為主要大氣污染物。11日和12日晚,顆粒物濃度出現(xiàn)爆發(fā)性增長(zhǎng),2019年1月12日21:00 PM2.5濃度高達(dá)589 μg·m-3,次日凌晨出現(xiàn)降低趨勢(shì)。由圖6b可見(jiàn),能見(jiàn)度與相對(duì)濕度變化趨勢(shì)相反,相對(duì)濕度是能見(jiàn)度的主要影響因子之一[35]。污染天氣下相對(duì)濕度可達(dá)90%,對(duì)應(yīng)1月9—14日能見(jiàn)度日平均值較低。1月13日04:00能見(jiàn)度最低達(dá)到0.10 km,此時(shí)相對(duì)濕度較大,溫度偏低,有利于霧的形成,導(dǎo)致能見(jiàn)度進(jìn)一步降低。根據(jù)地表面風(fēng)場(chǎng)數(shù)據(jù)可知,1月11日風(fēng)速小于2 m·s-1, 污染物累積不易擴(kuò)散,顆粒物濃度增大,污染事件加重。13日風(fēng)速增大,擴(kuò)散條件得到改善,污染事件有所好轉(zhuǎn)。14—15日地面風(fēng)速較大,有利于污染物擴(kuò)散,污染事件結(jié)束。

    表2為2019年1月8—15日北京城區(qū)不同波段τ日平均值,以500 nm為例分析污染事件中τ變化。由表2可以看到,τ隨波長(zhǎng)的增大而減小,表明氣溶膠粒子的消光性具有波長(zhǎng)選擇性。1月8日τ為0.10,天氣狀況良好,而9日τ上升到0.68,表明北京城區(qū)整體大氣消光性開(kāi)始增加,污染事件開(kāi)始。10日τ有所下降,又迅速上升至11日的0.89,近地面風(fēng)速較小,污染物不易擴(kuò)散導(dǎo)致污染加重。12日北京城區(qū)污染持續(xù),τ高達(dá)0.71,相對(duì)濕度接近90%。由圖6a可知,1月12日18:00—13日07:00 PM2.5濃度均大于300 μg·m-3,是整個(gè)污染事件中最嚴(yán)重的時(shí)段。13日10:00—17:00,PM2.5濃度迅速下降至75 μg·m-3以下,受西北冷空氣的滲透,擴(kuò)散條件改善。14日由于相對(duì)濕度較大,導(dǎo)致氣溶膠吸濕增長(zhǎng),消光能力增強(qiáng),τ增加至0.59。15日擴(kuò)散條件增強(qiáng),τ低至0.15,污染事件結(jié)束。

    圖5 2019年1月北京城區(qū)500 nm氣溶膠光學(xué)厚度和PM2.5濃度日平均值變化(陰影部分代表污染事件)Fig.5 Daily averaged aerosol optical depth at 500 nm and PM2.5 in Beijing from 1 Jan to 31 Jan in 2019(the shaded denotes pollution period)

    表3為2019年1月8—15日北京城區(qū)不同波段ω的日平均值,以500 nm為例分析污染事件中ω變化。1月8—10日ω呈上升趨勢(shì),由0.93上升至0.99,表明氣溶膠粒子散射能力增強(qiáng)。1月11日ω下降至0.98,是由于當(dāng)日擴(kuò)散條件極差,吸收性氣溶膠(如黑碳?xì)馊苣z)累積所致。1月12日ω上升至0.99,是由于相對(duì)濕度偏高,氣溶膠粒子吸濕增長(zhǎng)使得半徑增大,散射能力增強(qiáng)。1月14日ω較高值同樣也受到氣溶膠粒子吸濕增長(zhǎng)影響。1月15日ω下降至0.84,表明大氣中具有強(qiáng)吸收能力的氣溶膠粒子增加。

    圖7為2019年1月8—15日細(xì)、粗粒子占比的變化情況,表示α及細(xì)粒子與粗粒子的比例γ(PM2.5與PM10濃度之比)隨時(shí)間的變化趨勢(shì)。1月8—13日γ從0.37上升至0.92,處于平穩(wěn)上升階段,表明近地面氣溶膠細(xì)粒子逐漸增多。1月13—15日γ迅速下降, 15日γ低至0.17。污染期間α處于波動(dòng)變化趨勢(shì),1月9日α降低至0.82,粗粒子數(shù)量增多。1月10日α增加至1.09,說(shuō)明此時(shí)大氣中細(xì)粒子比例增加。1月11日α為0.91,是由于低風(fēng)速產(chǎn)生的不利擴(kuò)散條件導(dǎo)致氣溶膠粒子累積,使粗粒子數(shù)量也同時(shí)增加。1月12日α為1.16,此時(shí)相對(duì)濕度迅速升高,細(xì)粒子開(kāi)始吸濕增長(zhǎng),從污染累積階段逐漸進(jìn)入污染爆發(fā)階段。1月14—15日α由1.00減小至0.45,該時(shí)期氣溶膠細(xì)粒子占比顯著降低,粗粒子占比較大。

    圖6 2019年1月8—15日顆粒物濃度和氣象要素變化(a)PM2.5與PM10濃度變化,(b)能見(jiàn)度與相對(duì)濕度變化Fig.6 Temporal variation of particulate concentration and meteorological elements from 8 Jan to 15 Jan in 2019(a)PM2.5 and PM10,(b)visibility and relative humidity

    表2 2019年1月8—15日氣溶膠光學(xué)厚度統(tǒng)計(jì)Table 2 Statistics of aerosol optical depth from 8 Jan to 15 Jan in 2019

    表3 2019年1月8—15日單次散射反照率統(tǒng)計(jì)Table 3 Statistics of single scattering albedo from 8 Jan to 15 Jan in 2019

    圖7 2019年1月8—15日北京城區(qū)?ngstr?m波長(zhǎng)指數(shù)逐日變化Fig.7 Daily averaged variation of ?ngstr?m exponent from skyradiometer in Beijing from 8 Jan to 15 Jan in 2019

    由圖8可以看到,1月8日氣溶膠粒子體積濃度都非常低。1月9日污染天氣條件下,體積譜呈雙峰分布,其中粗模態(tài)粒子體積濃度大約是細(xì)模態(tài)粒子體積濃度的4倍。1月10日仍為雙峰分布,但粗、細(xì)模態(tài)粒子體積濃度相當(dāng),細(xì)模態(tài)氣溶膠粒子半徑主要集中在0.17 μm,粗模態(tài)粒子半徑主要集中在3.62 μm。1月11日為雙峰分布,而粗模態(tài)粒子體積濃度爆發(fā)性增長(zhǎng)到0.23 μm3·μm-2,細(xì)模態(tài)粒子體積濃度也增長(zhǎng)到0.07 μm3·μm-2,這是由于污染累積階段不利的擴(kuò)散條件導(dǎo)致。1月12日呈三峰型分布特征,其中粒子半徑主要集中在0.17 μm、1.69 μm 和5.29 μm。該結(jié)果表明在高相對(duì)濕度的條件下,氣溶膠粒子易吸濕增長(zhǎng),使得半徑逐漸增大。1月14日呈雙峰分布特征,吸濕增長(zhǎng)同樣導(dǎo)致氣溶膠粒子半徑增大,細(xì)模態(tài)粒子半徑主要集中在0.17~0.25 μm,粗模態(tài)粒子半徑主要集中在5.29 μm。1月15日呈單峰型分布,細(xì)模態(tài)粒子占比非常少,粗模態(tài)粒子半徑集中在5.29 μm左右,峰值半徑對(duì)應(yīng)的體積濃度高達(dá)0.15 μm3·μm-2。該結(jié)論表明:1月15日氣溶膠粗粒子和細(xì)粒子比例嚴(yán)重失調(diào),與1月15日α異常偏小相對(duì)應(yīng)。

    圖9為2019年1月9—14日一次典型北京冬季污染事件過(guò)程中基于激光雷達(dá)觀測(cè)得到的532 nm 氣溶膠消光系數(shù)及退偏振比分布,可以看到,1月9日00:00在0.5 km高度附近出現(xiàn)氣溶膠消光層,消光系數(shù)增大,粒子退偏振比也有所增加。11:00 在1.5 km 高度以下消光系數(shù)大于0.40 km-1, 退偏振比為0.05~0.10;在3.5~5.0 km 高度消光系數(shù)大于0.75 km-1,退偏振比小于0.10,可能是云層的影響。1月9日18:00—10日04:00,1.5~3.0 km高度消光系數(shù)變化范圍為0.30~0.90 km-1,退偏振比小于0.10,該區(qū)域可能為云和氣溶膠共存。此后,氣溶膠消光層高度出現(xiàn)下降并擴(kuò)散,消光系數(shù)高值區(qū)(大于0.60 km-1)呈帶狀緊貼近地面,近地面污染事件發(fā)生。1月11日00:00在2.5 km高度以上區(qū)域的消光系數(shù)大于0.45 km-1,退偏振比大于0.25,可能是不規(guī)則冰晶粒子組成的冰云影響。1月11日消光系數(shù)逐漸增大,達(dá)0.99 km-1以上,退偏振比無(wú)明顯變化,氣溶膠粒子濃度不斷增大,處于污染物累積階段。1月11日19:00氣溶膠消光層開(kāi)始降低,近地面消光系數(shù)可達(dá)0.90 km-1以上,氣溶膠粒子濃度顯著增加,進(jìn)入污染爆發(fā)階段。1月12日氣溶膠消光系數(shù)維持在高值范圍,退偏振比穩(wěn)定在1.0 km以內(nèi)。直到1月13日凌晨,消光系數(shù)迅速減小。在1月13日10:00在0.5 km高度附近出現(xiàn)氣溶膠消光帶(退偏振比大于0.20,消光系數(shù)小于0.30 km-1),由于近地面風(fēng)速較大導(dǎo)致地面揚(yáng)塵。1月13日18:00消光系數(shù)在0.5 km高度附近迅速增大至0.90 km-1以上,退偏振比小于0.10,說(shuō)明北京城區(qū)再次受到氣溶膠污染。1月14日19:00受到強(qiáng)西北風(fēng)的影響,風(fēng)速增大有利于污染物擴(kuò)散,消光系數(shù)迅速減小,另外西北風(fēng)有利于將西北地區(qū)的沙塵粒子輸送到北京,造成非球形粒子增多,因此,在1.5 km高度處退偏振比明顯增加。

    圖8 2019年1月8—15日北京城區(qū)粒子譜分布變化Fig.8 Averaged volume size from skyradiometer in Beijing from 8 Jan to 15 Jan in 2019

    圖9 2019年1月9—14日北京城區(qū)532 nm氣溶膠消光系數(shù)(a)和退偏振比(b)分布Fig.9 Temporal and spatial distribution of extinction coefficient(a) and the depolarization ratio(b) at 532 nm in Beijing from 9 Jan to 14 Jan in 2019

    3 結(jié) 論

    通過(guò)2018年10月—2019年9月天空輻射計(jì)連續(xù)觀測(cè)數(shù)據(jù),反演得到北京城區(qū)氣溶膠光學(xué)特性參數(shù)包括τ,ω,α和粒子譜分布特征,結(jié)合激光雷達(dá)重點(diǎn)分析了2019年1月一次典型污染事件過(guò)程中氣溶膠光學(xué)特性變化及其與氣象條件的相互聯(lián)系。分析結(jié)果表明:

    1) 北京城區(qū)氣溶膠光學(xué)厚度春夏季較大,秋冬季節(jié)較小;夏季氣溶膠散射能力強(qiáng),春季較弱;春冬季粗粒子較多,夏季細(xì)粒子較多。氣溶膠粒子體積譜呈雙峰型分布,粗、細(xì)粒子比會(huì)影響體積譜分布,而粒子的吸收、散射能力對(duì)體積譜分布影響不顯著。

    2) 污染天氣條件下PM2.5濃度是清潔天氣下PM2.5濃度平均值的2.27倍;污染天氣和清潔天氣條件下氣溶膠光學(xué)厚度平均值分別為0.85和0.49;氣溶膠粒子在污染天氣下具有更強(qiáng)的散射能力。受冬季霾污染和春季沙塵事件影響,冬季污染天氣條件下細(xì)粒子更多,春季污染天氣條件下粗粒子更多。

    3) 分析2019年1月8—15日一次典型污染事件,低風(fēng)速、高濕度等不利氣象條件、氣溶膠吸濕增長(zhǎng)和二次轉(zhuǎn)化、污染物局地排放及區(qū)域輸送共同導(dǎo)致嚴(yán)重污染事件的發(fā)生。

    猜你喜歡
    氣溶膠城區(qū)光學(xué)
    長(zhǎng)沙市望城區(qū)金地三千府幼兒園
    滑輪組的裝配
    光學(xué)常見(jiàn)考題逐個(gè)擊破
    氣溶膠傳播之謎
    金霞早油蟠在保定滿城區(qū)的表現(xiàn)及栽培技術(shù)
    氣溶膠中210Po測(cè)定的不確定度評(píng)定
    福田要建健康城區(qū)
    四川盆地秋季氣溶膠與云的相關(guān)分析
    對(duì)淮安市城區(qū)河道生態(tài)建設(shè)的思考
    光學(xué)遙感壓縮成像技術(shù)
    一级毛片我不卡| 国产在线一区二区三区精| 伦精品一区二区三区| 高清午夜精品一区二区三区| 搡女人真爽免费视频火全软件| 一级毛片电影观看| a 毛片基地| 国产在线一区二区三区精| 九草在线视频观看| 91aial.com中文字幕在线观看| 亚洲自偷自拍三级| 春色校园在线视频观看| 乱码一卡2卡4卡精品| h日本视频在线播放| 深夜a级毛片| 97在线视频观看| 五月天丁香电影| 丝袜喷水一区| 大话2 男鬼变身卡| 久久久久久伊人网av| 91精品伊人久久大香线蕉| 丝袜在线中文字幕| kizo精华| 久久精品国产亚洲av涩爱| 久久精品熟女亚洲av麻豆精品| 80岁老熟妇乱子伦牲交| 亚洲性久久影院| 午夜日本视频在线| 日韩中文字幕视频在线看片| 97在线视频观看| 国产成人a∨麻豆精品| 国产成人免费无遮挡视频| 一区二区三区乱码不卡18| 一级黄片播放器| 欧美精品高潮呻吟av久久| 色网站视频免费| 亚洲欧美成人综合另类久久久| 亚洲,一卡二卡三卡| 国产精品伦人一区二区| 久久精品国产亚洲av涩爱| 久久久久精品性色| 免费黄色在线免费观看| tube8黄色片| 欧美区成人在线视频| 国产伦精品一区二区三区视频9| 久久 成人 亚洲| 亚洲av不卡在线观看| 97超视频在线观看视频| 午夜91福利影院| 亚洲国产成人一精品久久久| 中文字幕久久专区| 一个人看视频在线观看www免费| 黄色日韩在线| av在线观看视频网站免费| 国产亚洲午夜精品一区二区久久| 少妇熟女欧美另类| 亚洲久久久国产精品| 如日韩欧美国产精品一区二区三区 | 中文字幕精品免费在线观看视频 | 国产精品久久久久久精品古装| 亚洲性久久影院| 亚洲欧美精品专区久久| 日本爱情动作片www.在线观看| 国产精品伦人一区二区| av有码第一页| 十分钟在线观看高清视频www | 人妻制服诱惑在线中文字幕| 久久久久久久久大av| 性色avwww在线观看| 新久久久久国产一级毛片| 99精国产麻豆久久婷婷| 七月丁香在线播放| 国产精品嫩草影院av在线观看| 国产一区有黄有色的免费视频| 免费观看在线日韩| 亚洲欧洲国产日韩| 中国国产av一级| 国产精品女同一区二区软件| 人人妻人人爽人人添夜夜欢视频 | 男的添女的下面高潮视频| av播播在线观看一区| 亚洲av电影在线观看一区二区三区| 国产精品99久久99久久久不卡 | 日韩av在线免费看完整版不卡| 寂寞人妻少妇视频99o| 少妇精品久久久久久久| 日韩一区二区三区影片| videos熟女内射| 18+在线观看网站| 九九在线视频观看精品| 人体艺术视频欧美日本| 国产精品熟女久久久久浪| 日本欧美国产在线视频| 欧美精品亚洲一区二区| 极品教师在线视频| 免费黄网站久久成人精品| 国产精品人妻久久久久久| 成人免费观看视频高清| 国产永久视频网站| 国产成人免费观看mmmm| 成人漫画全彩无遮挡| 六月丁香七月| 欧美日韩av久久| 91成人精品电影| 成人免费观看视频高清| 国产精品人妻久久久影院| 成年人午夜在线观看视频| 久久久精品免费免费高清| 日日啪夜夜撸| 国模一区二区三区四区视频| 国产伦精品一区二区三区四那| 亚洲电影在线观看av| 亚洲精品日韩av片在线观看| 日本与韩国留学比较| 一本一本综合久久| 久久精品国产亚洲av天美| 亚洲成人av在线免费| 亚洲一级一片aⅴ在线观看| 女性被躁到高潮视频| 国产av一区二区精品久久| 国产女主播在线喷水免费视频网站| 三级经典国产精品| 国产精品国产三级国产专区5o| 国产亚洲一区二区精品| 三上悠亚av全集在线观看 | 亚洲国产精品一区二区三区在线| 亚洲国产精品一区三区| 亚洲精品456在线播放app| 女性生殖器流出的白浆| 在线观看三级黄色| a级毛色黄片| 国产欧美亚洲国产| 各种免费的搞黄视频| 狂野欧美激情性bbbbbb| 亚洲国产日韩一区二区| 十分钟在线观看高清视频www | 欧美日韩综合久久久久久| 日韩三级伦理在线观看| 秋霞在线观看毛片| 精品酒店卫生间| 中国三级夫妇交换| 在线免费观看不下载黄p国产| 国产毛片在线视频| 亚洲国产精品专区欧美| 大香蕉97超碰在线| 91精品一卡2卡3卡4卡| 美女主播在线视频| 2021少妇久久久久久久久久久| 久久这里有精品视频免费| 黄片无遮挡物在线观看| 男男h啪啪无遮挡| 亚洲国产毛片av蜜桃av| 免费高清在线观看视频在线观看| 精品久久久精品久久久| kizo精华| 18禁裸乳无遮挡动漫免费视频| 精品人妻熟女av久视频| 国产一区二区三区综合在线观看 | 午夜福利影视在线免费观看| 中国美白少妇内射xxxbb| 少妇丰满av| 少妇 在线观看| 国产精品麻豆人妻色哟哟久久| 男的添女的下面高潮视频| 综合色丁香网| 免费人成在线观看视频色| 日本与韩国留学比较| videossex国产| 男人添女人高潮全过程视频| 观看av在线不卡| 精品久久久精品久久久| 亚洲国产精品999| www.色视频.com| a 毛片基地| 亚洲国产日韩一区二区| 中文精品一卡2卡3卡4更新| 国产日韩欧美视频二区| 国产男女超爽视频在线观看| 久久综合国产亚洲精品| 又黄又爽又刺激的免费视频.| 日韩av在线免费看完整版不卡| 日日摸夜夜添夜夜爱| 欧美一级a爱片免费观看看| 黄片无遮挡物在线观看| 男人和女人高潮做爰伦理| 亚洲欧美成人精品一区二区| 亚洲精品乱久久久久久| 国产69精品久久久久777片| 国产男女超爽视频在线观看| 18禁在线播放成人免费| 自线自在国产av| 91久久精品电影网| 国产成人a∨麻豆精品| 大又大粗又爽又黄少妇毛片口| 久久99蜜桃精品久久| 成人无遮挡网站| 久热这里只有精品99| 美女主播在线视频| 精品一区二区三区视频在线| 亚洲欧美日韩东京热| 五月玫瑰六月丁香| 亚洲国产成人一精品久久久| 国产亚洲5aaaaa淫片| 丰满少妇做爰视频| 国产精品99久久久久久久久| 欧美精品一区二区免费开放| 亚洲情色 制服丝袜| 国产成人aa在线观看| 黑人巨大精品欧美一区二区蜜桃 | 六月丁香七月| 亚洲国产成人一精品久久久| 日日摸夜夜添夜夜添av毛片| 男女免费视频国产| 九九久久精品国产亚洲av麻豆| 偷拍熟女少妇极品色| 国产中年淑女户外野战色| 观看免费一级毛片| 国产黄片美女视频| 精华霜和精华液先用哪个| 一级毛片久久久久久久久女| 麻豆精品久久久久久蜜桃| 两个人免费观看高清视频 | 亚洲欧美一区二区三区国产| 国产黄片美女视频| 美女xxoo啪啪120秒动态图| 97在线人人人人妻| 欧美日韩视频精品一区| av女优亚洲男人天堂| 国产一区二区三区综合在线观看 | 少妇人妻 视频| 最近2019中文字幕mv第一页| 久久人妻熟女aⅴ| 午夜免费鲁丝| 80岁老熟妇乱子伦牲交| 最新中文字幕久久久久| 亚洲婷婷狠狠爱综合网| 丁香六月天网| 波野结衣二区三区在线| 男人舔奶头视频| 亚洲欧洲日产国产| 在现免费观看毛片| 一级爰片在线观看| 精品国产露脸久久av麻豆| 婷婷色综合大香蕉| 综合色丁香网| 国产69精品久久久久777片| 成人午夜精彩视频在线观看| 亚洲电影在线观看av| 九九久久精品国产亚洲av麻豆| 精品少妇久久久久久888优播| 美女视频免费永久观看网站| 2022亚洲国产成人精品| 五月天丁香电影| 亚洲av二区三区四区| 成人18禁高潮啪啪吃奶动态图 | 欧美高清成人免费视频www| 久久久久久久久大av| 国产精品伦人一区二区| 日产精品乱码卡一卡2卡三| 国产精品久久久久成人av| 国产一区二区三区综合在线观看 | 久久久久视频综合| 亚洲精品乱码久久久v下载方式| 亚洲精品亚洲一区二区| 国产在视频线精品| 亚洲av二区三区四区| 久久人人爽av亚洲精品天堂| 蜜桃在线观看..| 国产一区亚洲一区在线观看| 伊人久久国产一区二区| 国产欧美日韩综合在线一区二区 | 亚洲丝袜综合中文字幕| 婷婷色麻豆天堂久久| 一级a做视频免费观看| 人妻人人澡人人爽人人| tube8黄色片| 国产av码专区亚洲av| 精品久久久久久电影网| 丰满少妇做爰视频| 乱码一卡2卡4卡精品| 午夜福利在线观看免费完整高清在| 女人精品久久久久毛片| 免费观看的影片在线观看| 久久狼人影院| 91午夜精品亚洲一区二区三区| 久久人人爽av亚洲精品天堂| 精品99又大又爽又粗少妇毛片| 波野结衣二区三区在线| 亚洲国产欧美在线一区| 免费黄网站久久成人精品| 亚洲欧美精品专区久久| 99热全是精品| 22中文网久久字幕| 国产乱人偷精品视频| 亚洲精品久久久久久婷婷小说| 亚洲欧美清纯卡通| 亚洲成人av在线免费| 中文字幕av电影在线播放| 国产精品不卡视频一区二区| 男女国产视频网站| 在线观看免费视频网站a站| 色吧在线观看| 久久久久久久久久久丰满| 观看av在线不卡| 麻豆精品久久久久久蜜桃| 欧美精品人与动牲交sv欧美| 午夜福利影视在线免费观看| 在线天堂最新版资源| 国产免费视频播放在线视频| 不卡视频在线观看欧美| 欧美性感艳星| 欧美+日韩+精品| a级毛片免费高清观看在线播放| 成人美女网站在线观看视频| 国产欧美另类精品又又久久亚洲欧美| 99久久人妻综合| 精品久久久久久久久亚洲| 日韩一区二区三区影片| 久久99精品国语久久久| 三级国产精品片| 久久久亚洲精品成人影院| 在现免费观看毛片| 免费观看av网站的网址| 中国国产av一级| 久久精品国产鲁丝片午夜精品| 高清黄色对白视频在线免费看 | 久久狼人影院| 欧美亚洲 丝袜 人妻 在线| videossex国产| 岛国毛片在线播放| 午夜福利视频精品| 久久女婷五月综合色啪小说| 99九九在线精品视频 | 色婷婷久久久亚洲欧美| 久久99一区二区三区| 婷婷色综合www| 人体艺术视频欧美日本| 春色校园在线视频观看| 午夜福利在线观看免费完整高清在| 精品久久久精品久久久| 一级毛片黄色毛片免费观看视频| 一本大道久久a久久精品| 国产av精品麻豆| videossex国产| 久久久国产精品麻豆| 91在线精品国自产拍蜜月| 亚洲国产精品999| 久久影院123| 美女主播在线视频| 热re99久久精品国产66热6| av在线观看视频网站免费| 国产中年淑女户外野战色| 91午夜精品亚洲一区二区三区| 简卡轻食公司| 久久热精品热| 最后的刺客免费高清国语| 亚洲精华国产精华液的使用体验| 亚洲图色成人| 久久久久久久久久久免费av| 国产男人的电影天堂91| 狂野欧美白嫩少妇大欣赏| 国产高清不卡午夜福利| 成人漫画全彩无遮挡| 观看免费一级毛片| 汤姆久久久久久久影院中文字幕| 午夜福利视频精品| 亚洲精品一区蜜桃| 成人黄色视频免费在线看| 欧美日韩国产mv在线观看视频| 久久婷婷青草| 少妇人妻 视频| 亚洲精品国产av成人精品| 尾随美女入室| 丰满人妻一区二区三区视频av| 少妇熟女欧美另类| 看免费成人av毛片| 久久人人爽人人爽人人片va| 香蕉精品网在线| 交换朋友夫妻互换小说| 校园人妻丝袜中文字幕| 久久久久久久久久久丰满| 国产精品国产三级专区第一集| 能在线免费看毛片的网站| 日日摸夜夜添夜夜添av毛片| av免费在线看不卡| 人妻系列 视频| 制服丝袜香蕉在线| 精品少妇久久久久久888优播| 欧美精品高潮呻吟av久久| 99久久中文字幕三级久久日本| 观看免费一级毛片| 一边亲一边摸免费视频| 亚洲精华国产精华液的使用体验| 久久久久久久精品精品| 久久久久久久久久久免费av| 春色校园在线视频观看| 22中文网久久字幕| 欧美丝袜亚洲另类| 国产欧美日韩综合在线一区二区 | 国产av精品麻豆| 99久久精品国产国产毛片| 大话2 男鬼变身卡| 国产精品一区www在线观看| 伊人亚洲综合成人网| 国产美女午夜福利| 国产精品偷伦视频观看了| 国产一级毛片在线| 精品一区二区免费观看| 一级爰片在线观看| 亚洲精品国产成人久久av| 色5月婷婷丁香| 天堂8中文在线网| 亚洲国产日韩一区二区| 最后的刺客免费高清国语| 国产成人精品久久久久久| 国产精品一区www在线观看| 久久久午夜欧美精品| 婷婷色av中文字幕| 国产一区二区在线观看av| tube8黄色片| 九色成人免费人妻av| 国产又色又爽无遮挡免| 国产黄片视频在线免费观看| 大陆偷拍与自拍| 大香蕉97超碰在线| 久久久午夜欧美精品| 一级av片app| 最新中文字幕久久久久| 特大巨黑吊av在线直播| a级毛片在线看网站| 精品一区在线观看国产| 一级毛片我不卡| 亚洲,欧美,日韩| 丰满迷人的少妇在线观看| 人人妻人人爽人人添夜夜欢视频 | 一本一本综合久久| 在线天堂最新版资源| 久久6这里有精品| 亚洲国产精品999| av一本久久久久| a级一级毛片免费在线观看| 黑人高潮一二区| 99热全是精品| 国产黄色免费在线视频| 韩国av在线不卡| 亚洲精品日韩在线中文字幕| 亚洲国产精品成人久久小说| 欧美日本中文国产一区发布| 日韩成人av中文字幕在线观看| 国产精品人妻久久久久久| 丰满乱子伦码专区| 国产淫语在线视频| 亚洲欧洲国产日韩| 少妇 在线观看| 欧美三级亚洲精品| 99久久综合免费| 国产欧美亚洲国产| 国产高清国产精品国产三级| 夜夜爽夜夜爽视频| 精品久久久久久久久亚洲| 久久久久国产精品人妻一区二区| av天堂中文字幕网| av国产精品久久久久影院| 在线播放无遮挡| 日韩中文字幕视频在线看片| 日韩中字成人| 国产成人aa在线观看| 国产真实伦视频高清在线观看| 老熟女久久久| 一级a做视频免费观看| 久久久久久久久久成人| 一区二区三区乱码不卡18| 国内少妇人妻偷人精品xxx网站| 青春草国产在线视频| 中文欧美无线码| 国产乱人偷精品视频| 一边亲一边摸免费视频| 不卡视频在线观看欧美| 免费观看的影片在线观看| 亚洲精品乱码久久久久久按摩| 亚洲精品中文字幕在线视频 | 成人二区视频| 久久久久久伊人网av| 精品一品国产午夜福利视频| 天美传媒精品一区二区| 日本免费在线观看一区| 国产日韩欧美视频二区| 亚洲av成人精品一二三区| 中文字幕久久专区| 视频中文字幕在线观看| 美女脱内裤让男人舔精品视频| 国产成人aa在线观看| 精品久久久久久久久亚洲| 午夜影院在线不卡| 偷拍熟女少妇极品色| 国产高清不卡午夜福利| 午夜激情福利司机影院| 免费播放大片免费观看视频在线观看| 男女无遮挡免费网站观看| 观看美女的网站| 最近中文字幕2019免费版| 最近的中文字幕免费完整| 国产一区二区在线观看日韩| 精品一区二区三卡| 校园人妻丝袜中文字幕| 欧美三级亚洲精品| 国产无遮挡羞羞视频在线观看| 亚洲伊人久久精品综合| 人妻 亚洲 视频| 久久人人爽av亚洲精品天堂| 全区人妻精品视频| 国产精品一二三区在线看| 日本vs欧美在线观看视频 | 久久6这里有精品| 三级国产精品欧美在线观看| 亚洲精品中文字幕在线视频 | 亚洲精品乱码久久久v下载方式| 免费久久久久久久精品成人欧美视频 | a级毛片免费高清观看在线播放| 国产成人精品福利久久| 国产淫片久久久久久久久| 久久国产乱子免费精品| 成人综合一区亚洲| 国精品久久久久久国模美| 高清黄色对白视频在线免费看 | 爱豆传媒免费全集在线观看| 亚洲,欧美,日韩| 99久久精品国产国产毛片| 街头女战士在线观看网站| 精品人妻一区二区三区麻豆| av天堂久久9| 免费av不卡在线播放| 婷婷色综合www| 午夜福利在线观看免费完整高清在| 欧美日韩视频高清一区二区三区二| 色婷婷久久久亚洲欧美| 日本免费在线观看一区| 国产成人免费无遮挡视频| av有码第一页| 我要看黄色一级片免费的| 欧美精品一区二区免费开放| 国产一区有黄有色的免费视频| 少妇裸体淫交视频免费看高清| 天堂中文最新版在线下载| 亚洲成人一二三区av| 在线观看三级黄色| 国产精品一区二区三区四区免费观看| 国产精品久久久久久精品电影小说| 日韩伦理黄色片| 两个人免费观看高清视频 | 女人久久www免费人成看片| 国产精品一区www在线观看| 亚洲欧洲国产日韩| 国产精品人妻久久久久久| 一级毛片我不卡| 国产午夜精品久久久久久一区二区三区| 在线观看人妻少妇| 精品少妇久久久久久888优播| 人妻夜夜爽99麻豆av| 在线亚洲精品国产二区图片欧美 | 午夜激情福利司机影院| 我的女老师完整版在线观看| 亚洲在久久综合| 亚洲欧美中文字幕日韩二区| 一区二区av电影网| 中文资源天堂在线| 在现免费观看毛片| 最新中文字幕久久久久| 欧美高清成人免费视频www| 夜夜骑夜夜射夜夜干| 91aial.com中文字幕在线观看| 日韩视频在线欧美| 国产精品久久久久久av不卡| 国产亚洲5aaaaa淫片| 国产日韩欧美在线精品| 国产成人免费无遮挡视频| 亚洲av国产av综合av卡| 久久久a久久爽久久v久久| 91在线精品国自产拍蜜月| 日韩人妻高清精品专区| 91成人精品电影| 久久精品久久精品一区二区三区| 在线观看av片永久免费下载| 久久久久久久久久久丰满| 国产综合精华液| 人妻夜夜爽99麻豆av| 黑人巨大精品欧美一区二区蜜桃 | 国产精品三级大全| 欧美日韩综合久久久久久| 男人狂女人下面高潮的视频| 日本猛色少妇xxxxx猛交久久| 国产精品一区二区在线观看99| 国产亚洲一区二区精品| 又粗又硬又长又爽又黄的视频| 国产黄片美女视频| 欧美丝袜亚洲另类| 亚洲欧洲日产国产| 人体艺术视频欧美日本| 国产精品蜜桃在线观看| 成人影院久久| 欧美日韩综合久久久久久| 亚洲国产欧美日韩在线播放 | 性高湖久久久久久久久免费观看| 亚洲精品成人av观看孕妇| 五月玫瑰六月丁香| 亚洲成色77777| 久久精品久久久久久久性| 有码 亚洲区| 插逼视频在线观看| 三上悠亚av全集在线观看 | 国产高清不卡午夜福利| 亚洲精品456在线播放app| 欧美另类一区| 肉色欧美久久久久久久蜜桃| 熟妇人妻不卡中文字幕|