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

    天津大氣穩(wěn)定度和逆溫特征對PM2.5污染的影響

    2018-08-23 05:42:12劉敬樂蔡子穎韓素芹馬志強
    中國環(huán)境科學(xué) 2018年8期
    關(guān)鍵詞:逆溫穩(wěn)定度均值

    姚 青,劉敬樂,蔡子穎,韓素芹,馬志強

    ?

    天津大氣穩(wěn)定度和逆溫特征對PM2.5污染的影響

    姚 青1*,劉敬樂2,蔡子穎1,韓素芹2,馬志強3

    (1.天津市環(huán)境氣象中心,天津 300074;2.天津市氣象科學(xué)研究所,天津 300074;3.京津冀環(huán)境氣象預(yù)報預(yù)警中心,北京 100089)

    利用氣象鐵塔資料分析了逆溫頻率和強度,采用溫差-風(fēng)速法計算了天津地區(qū)大氣穩(wěn)定度,探討了其相互關(guān)系及對PM2.5濃度月均值和超標(biāo)日的影響,并對一次重污染過程中大氣穩(wěn)定度和逆溫分布特征進行了分析.結(jié)果表明, 2015年9月~2017年8月A,B,C,D,E和F類大氣穩(wěn)定度發(fā)生頻率依次為6.7%,11.4%,22.4%,46.1%,11.1%和2.2%,秋冬季節(jié)穩(wěn)定類天氣(E,F類)較多,全年白天各時段均以不穩(wěn)定類大氣為主,夜間大氣穩(wěn)定度以中性為主,秋冬季夜間穩(wěn)定類高達30%~40%.觀測期內(nèi)冬季逆溫頻率最高, 5:00~8:00和21:00~23:00超過90%,冬季逆溫強度也最高.隨著穩(wěn)定類大氣層結(jié)日數(shù)的增多,PM2.5月均值和污染日數(shù)均有所增大,同時逐月PM2.5均值、污染日發(fā)生頻率均與逆溫發(fā)生頻率呈正相關(guān)關(guān)系.2016年12月16~21日的一次重污染天氣過程顯示,PM2.5受到大氣穩(wěn)定度和逆溫發(fā)展的影響,霾形成、霧-霾交替和消散等階段大氣穩(wěn)定度和逆溫特征具有顯著的不同.大氣持續(xù)趨于穩(wěn)定及逆溫強度的逐漸增大,對污染生成和維持起了非常重要的作用,污染過程中大氣穩(wěn)定度和逆溫特征的精細化分析有助于提升重污染天氣預(yù)報預(yù)警水平.

    大氣穩(wěn)定度;逆溫;重污染;PM2.5;天津

    近年來伴隨著工業(yè)和城市交通快速發(fā)展,大量化石能源消耗,京津冀、長三角、珠三角等城市群呈現(xiàn)區(qū)域性復(fù)合型大氣污染特征,突出表現(xiàn)為秋冬季節(jié)中東部地區(qū)霧霾天氣頻發(fā),尤其是京津冀地區(qū),由于受山地、渤海等特殊地形和特大城市群構(gòu)成的復(fù)雜下墊面影響,區(qū)域和局地氣象條件、近地層環(huán)流特征特別是邊界層結(jié)構(gòu)具有獨特性,造成了污染物在邊界層內(nèi)的局地擴散以及區(qū)域輸送的變化規(guī)律呈現(xiàn)復(fù)雜多變態(tài)勢.大量研究[1-5]表明,京津冀區(qū)域重污染天氣常伴有大氣靜穩(wěn)天氣,較強的逆溫層,以及較強的近地層水汽輸送,此類不利氣象條件易導(dǎo)致污染物的累積,從而形成重污染過程,并突出性的表現(xiàn)為持續(xù)霧霾天氣.氣象條件是造成2013年初我國中東部地區(qū)持續(xù)嚴重霧霾天氣的重要外因[1-2],氣溶膠分布特征[3-6],不同尺度環(huán)流形勢和邊界層結(jié)構(gòu)分析[7-9]及其對氣象條件的反饋作用[10-13]等被廣泛于霧霾天氣過程研究中.

    大氣穩(wěn)定度可作為衡量大氣污染嚴重程度的重要指標(biāo),污染物在大氣中擴散的速度、距離和范圍受到大氣穩(wěn)定度的影響,當(dāng)大氣層結(jié)處于穩(wěn)定狀態(tài)時,污染物不易擴散.大氣穩(wěn)定度有多種分類方法,康凌等[14]對針對戈壁沙漠、低矮丘陵山地和沿海地區(qū)等不同下墊面,分別采用采用P-T法、ΔT法、ΔT/U法、ΔT/U2法、輻射法和風(fēng)向標(biāo)準差法計算大氣穩(wěn)定度,分析了不同下墊面大氣穩(wěn)定度的適用方法.金莉莉等[15]利用烏魯木齊市4座100m梯度氣象塔的氣象資料,認為溫差-風(fēng)速法更加適合烏魯木齊的大氣穩(wěn)定度分類.楊靜等[16]的研究顯示烏魯木齊F類穩(wěn)定度的出現(xiàn)頻率與API相關(guān)系數(shù)0.80,是一個比較好的空氣污染潛勢預(yù)報因子.北京[17]、哈爾濱[18]等地的研究表明,逆溫是導(dǎo)致大氣層結(jié)穩(wěn)定和形成大氣污染的重要原因.

    天津地區(qū)有關(guān)學(xué)者利用氣象鐵塔資料,開展了有關(guān)邊界層低層風(fēng)溫資料[19]、粗糙度[20]、霧[21-22]和污染物[23-30]的初步研究,獲得了大量一手觀測資料,對于天津地區(qū)邊界層結(jié)構(gòu)和污染物垂直分布特征有了一定程度的了解.深入研究天津城市氣象條件,如大氣穩(wěn)定度、逆溫等的變化特征,以及對PM2.5等主要大氣污染物的影響,厘清重污染天氣形成機理,對于改進空氣質(zhì)量模式,提升本區(qū)域霧霾天氣等環(huán)境氣象業(yè)務(wù)預(yù)報預(yù)警水平,制定切實可行的大氣污染防控措施,具有重要的科學(xué)意義和技術(shù)價值.

    1 資料與方法

    本文采用的氣象數(shù)據(jù)來自中國氣象局天津大氣邊界層觀測站(北緯39°06¢,東經(jīng)117°10¢,海拔高度2.2m,臺站編號:54517)氣象鐵塔梯度氣象觀測系統(tǒng),該站點位于天津市城區(qū)南部,其北距快速路約100m,東臨友誼路-友誼南路,西面和南面主要為住宅區(qū),已開展逾30a的梯度氣象觀測,觀測資料已在天津地區(qū)邊界層氣象和大氣環(huán)境研究中得到諸多應(yīng)用[19-22].該系統(tǒng)在距地面5, 10, 20, 30, 40, 60, 80, 100, 120, 140, 160, 180, 200, 220, 250m觀測平臺上均開展風(fēng)速、風(fēng)向、溫度和相對濕度的觀測,采用中環(huán)天儀(天津)氣象儀器有限公司生產(chǎn)的DZZ6型自動氣象站,氣象資料經(jīng)過嚴格質(zhì)控后取小時平均值,質(zhì)控方法見參考文獻[31].PM2.5質(zhì)量濃度數(shù)據(jù)來源于天津市環(huán)境監(jiān)測中心對外發(fā)布的全市平均值,依據(jù)《環(huán)境空氣質(zhì)量指數(shù)(AQI)技術(shù)規(guī)定(試行)》[32]判定污染日數(shù),考慮到天津地區(qū)主要污染物為PM2.5,本文僅考慮PM2.5的污染狀況,以PM2.5質(zhì)量濃度日均值超過75μg/m3記做一個PM2.5超標(biāo)日.觀測期為2015年9月1日~2017年8月31日,數(shù)據(jù)經(jīng)過質(zhì)控,小時樣本數(shù)為17304,有效率為98.6%.能見度數(shù)據(jù)采用美國Belfort公司生產(chǎn)的MODLE6000型前向散射能見度儀測量,儀器測量上限為20km,下限為6m,精確度10%,能見度儀安裝于觀測站院內(nèi)的觀測場內(nèi),探頭離地高約2.2m.

    蔡子穎等[33]采用溫差-風(fēng)速法和凈輻射法分別計算了本站2015年10月~2017年2月的大氣穩(wěn)定度,并詳細分析了各種計算方法的優(yōu)缺點,結(jié)果表明采用溫差-風(fēng)速法可以充分利用氣象鐵塔的客觀觀測資料,并且較好的反映熱力湍流和機械湍流的影響.本文采用該方法計算大氣穩(wěn)定度,具體方法為根據(jù)鐵塔40m處的風(fēng)速,以及100m和30m處的溫度差,按照不同風(fēng)速區(qū)間,將溫度差插值進入表1,即可得到相應(yīng)的大氣穩(wěn)定度,其中A為強不穩(wěn)定,B為不穩(wěn)定,C為弱不穩(wěn)定,D為中性,E為較穩(wěn)定,F為穩(wěn)定.

    考慮到氣象鐵塔的基座為水泥地面,5m處溫度因受到地面輻射影響,與10m及以上平臺測定的溫度差異較大,本文中逆溫的判定以10m平臺的溫度開始,逐層計算,在10~250m共14層高度的溫度廓線中,如有上層溫度高于下層溫度,即記做一次逆溫,定義逆溫強度(Temperature Inversion,TI,單位為℃/100m)如下:

    TI=100*(2-1)/(2-1)

    式中:1<2,1為高度1下的溫度,2為高度2下的溫度,如2>1則記做逆溫,如一次逆溫事件中1=10m,則記做貼地逆溫,如1>10m,則記做脫地逆溫,本文不考慮多層逆溫的情況.考慮到溫度傳感器準確性為±0.1℃,如兩層間溫度差小于0.2℃,可能代表較弱的逆溫,也可能是儀器測量誤差造成.經(jīng)計算,逆溫差0.2℃時逆溫強度平均值為0.36℃,低于全體逆溫樣本時的逆溫強度平均值(0.82℃),表明仍以弱逆溫為主,因而本文保留這些數(shù)據(jù).

    表1 溫差-風(fēng)速法確定穩(wěn)定度的分類標(biāo)準

    2 結(jié)果與討論

    2.1 大氣穩(wěn)定度特征分析

    圖1 大氣穩(wěn)定度發(fā)生頻率的逐月分布特征

    觀測期內(nèi)A~F類大氣穩(wěn)定度的發(fā)生頻率依次為6.7%,11.4%,22.4%,46.1%,11.1%和2.2%,D類穩(wěn)定度(中性)最多,穩(wěn)定類天氣(E,F類)為13.3%,各類穩(wěn)定度發(fā)生頻率的分布規(guī)律較為合理.圖1給出了大氣穩(wěn)定度發(fā)生頻率的逐月分布情況,A~F類穩(wěn)定度發(fā)生頻率的范圍依次為1.8%~12.9%,6.8%~16.4%, 14.2%~45.7%,32.9%~58.1%,2.3%~24.3%和0~6.2%,其中2016年7月天津地區(qū)多降水天氣,強對流天氣頻發(fā),未出現(xiàn)F類穩(wěn)定度.逐月分布顯示,大氣穩(wěn)定度的季節(jié)變化差異較大,秋冬季節(jié)穩(wěn)定類天氣較多, 2015年10月E+F類發(fā)生頻率為30.5%,2015年12月、2016年2月、3月、10月、11月、12月,2017年2月發(fā)生頻率均在20%左右,較為穩(wěn)定的大氣是造成天津地區(qū)秋冬季節(jié)大氣擴散條件變差、污染物濃度較高和霧霾天氣頻發(fā)的重要氣象條件[34].

    受到晝夜間熱力和動力因子差異的影響,大氣穩(wěn)定度分布具有明顯的日變化(圖2).各季節(jié)白天大氣穩(wěn)定度均以不穩(wěn)定類為主,春季和夏季持續(xù)時間明顯高于秋季和冬季,春夏季7:00~17:00不穩(wěn)定類頻率超過50%,春季10:00~14:00這一比例甚至接近100%;秋季和冬季不穩(wěn)定類占優(yōu)勢地位的時間有所縮短,分別為8:00~ 16:00和9:00~16:00,可能與太陽輻射時間縮短、熱力對流強度減弱有關(guān).夜間大氣穩(wěn)定度以中性為主,穩(wěn)定類次之,其中春季中性類在60%以上,夏季則高達80%,秋季和冬季也維持在50%左右,穩(wěn)定類夏季最低,僅為10%左右,秋季和冬季高達30%~40%,根據(jù)溫差-風(fēng)速法的計算原理,可能與秋冬季頻繁發(fā)生的逆溫現(xiàn)象有關(guān).

    圖2 不同季節(jié)下不穩(wěn)定類、中性和穩(wěn)定類的日變化特征

    2.2 逆溫特征分析

    由圖3可見,冬季逆溫頻率最高,兩個冬季均在60%以上,其中2016年12月和2017年1月逆溫頻率超過80%,夏季逆溫頻率在40%以下,貼地逆溫頻率也有類似的季節(jié)分布規(guī)律,但冬季頻率沒有逆溫頻率突出,可能與冬季逆溫存在大量的脫地逆溫有關(guān).逆溫和貼地逆溫強度分布規(guī)律較為復(fù)雜,大體上呈現(xiàn)秋冬高、春夏低的特點,除了冬季逆溫強度較高外,2015年10月、2016年3月和2017年4月等月份的突出高值,可能與單次較強的逆溫事件有關(guān),需要結(jié)合邊界層氣象資料和污染事件進一步分析.

    圖3 逆溫和貼地逆溫發(fā)生頻率及強度的逐月分布特征

    由圖4可見,冬季在各個時段逆溫頻率均高于其他季節(jié),春、夏、秋三季則總體差別不大,冬季凌晨和夜間受地面輻射降溫影響,逆溫頻發(fā),5:00~8:00和21:00~23:00逆溫頻率超過90%,秋季次之,春季和夏季的白天偶有逆溫發(fā)生.貼地逆溫沒有逆溫的季節(jié)差異明顯,凌晨各時段貼地逆溫頻率均在50%左右,白天8:00~15:00各季節(jié)貼地逆溫發(fā)生頻率接近0,這與日出后,地面升溫,破壞了貼地逆溫層有關(guān),日落后貼地逆溫頻率逐漸增大,冬季>秋季>春季>夏季,可能與不同季節(jié)傍晚地面降溫時間先后和降溫幅度不同有關(guān).圖5給出了各季節(jié)逆溫和貼地逆溫強度的日變化狀況,考慮到白天部分時段逆溫和貼地逆溫發(fā)生頻率較低(10%),樣本數(shù)不足,在統(tǒng)計逆溫強度時未統(tǒng)計該部分數(shù)據(jù),主要是春夏秋三季11:00~16:00的逆溫強度和8:00~16:00的貼地逆溫強度,從圖5可見,一般而言秋冬季逆溫強度高于春夏季,且夜間(19:00~次日4:00)冬季逆溫強度最高,秋季次之,春季和夏季最低,逆溫強度和逆溫頻率的季節(jié)差異類似,反映了太陽輻射的季節(jié)變化特征.

    圖5 各季節(jié)逆溫和貼地逆溫強度的日變化特征

    2.3 大氣穩(wěn)定度、逆溫特征與污染的關(guān)系

    觀測期間天津全市PM2.5質(zhì)量濃度均值為(71±60)μg/m3,小時濃度中位值為53μg/m3,最大小時值為509μg/m3,發(fā)生在2016年2月8日1:00.圖6給出了觀測期間PM2.5質(zhì)量濃度的逐月均值和逐月超標(biāo)日數(shù),可見大體上呈現(xiàn)冬季高、夏季低的特征,2a的峰值均出現(xiàn)在12月,正值天津地區(qū)2015年和2016年發(fā)布的重污染紅色預(yù)警期間,持續(xù)嚴重的重污染天氣對推高PM2.5月均值的峰值起到了重要作用,此外2016年3月出現(xiàn)了一個較小的峰值,這與當(dāng)月我國北方出現(xiàn)的幾次沙塵天氣過程有關(guān)[35].按照《環(huán)境空氣質(zhì)量指數(shù)(AQI)技術(shù)規(guī)定(試行)》[32]統(tǒng)計各月的超標(biāo)日數(shù),可以發(fā)現(xiàn)2015年11月、12月,2016年11月、12月和2017年1月PM2.5超標(biāo)日數(shù)較高,均超過全月日數(shù)的50%,其中2016年12月最高,超標(biāo)日數(shù)為22d,2015年12月和2016年11月次之,為20d,夏季超標(biāo)日數(shù)較少,一般在6d以內(nèi).超標(biāo)日數(shù)的月分布規(guī)律與PM2.5均值的月分布相一致,這表明污染天氣頻發(fā)對PM2.5濃度均值的拉升作用,對PM2.5的治理策略應(yīng)首先集中在降低PM2.5污染程度和頻度上.

    圖6 觀測期間PM2.5質(zhì)量濃度的逐月均值和逐月超標(biāo)日數(shù)

    采用溫差-風(fēng)速法計算的大氣穩(wěn)定度既反映了動力特征(風(fēng)速),也反映了大氣的熱力特征(溫度層結(jié)),從表1可以看出,如近地層存在逆溫(100>30),則易出現(xiàn)E和F類穩(wěn)定度,在大風(fēng)天氣下(37.0m/s)則均為D類穩(wěn)定度.圖7給出了不同大氣穩(wěn)定度下逆溫頻率、貼地逆溫頻率和貼地逆溫占逆溫比的關(guān)系,可見隨著逆溫和貼地逆溫發(fā)生概率增加,大氣趨于穩(wěn)定時,逆溫發(fā)生概率不足20%,幾乎沒有貼地逆溫發(fā)生,大氣多為不穩(wěn)定層結(jié);當(dāng)逆溫和貼地逆溫發(fā)生頻率增高至為60%和29%時,表現(xiàn)為中性層結(jié);逆溫發(fā)生概率接近100%,貼地逆溫在80%以上時,大氣為較穩(wěn)定層結(jié)(E類);而大氣進一步穩(wěn)定達到F類穩(wěn)定度時,逆溫和貼地逆溫頻率均有不同程度下降,這可能與逆溫層深厚,逆溫頂和逆溫底均向上發(fā)展導(dǎo)致大氣更加穩(wěn)定,而鐵塔資料僅為地面至250m高,高于這一高度的逆溫層底無法在統(tǒng)計資料中體現(xiàn),因而部分脫地逆溫可能被遺漏有關(guān).

    圖7 不同大氣穩(wěn)定度下逆溫頻率、貼地逆溫頻率和貼地逆溫占逆溫比

    由圖8可見,隨著穩(wěn)定類大氣層結(jié)日數(shù)的增多,PM2.5均值和污染日數(shù)均有所增大.由圖9可見,隨著逆溫發(fā)生頻率的增大,PM2.5均值和污染日數(shù)均有所增大,這與何永晴等[36]在西寧的研究結(jié)果一致.逆溫頻發(fā)月份PM2.5月均值和污染日也相應(yīng)較高,這可能與逆溫天氣下,大氣邊界層高度降低,湍流減弱,地面污染物的垂直擴散受到了直接的限制,更容易堆積在近地面有關(guān).相關(guān)性分析表明逆溫對空氣污染的影響強于穩(wěn)定度,這可能與大氣穩(wěn)定度的計算方法有關(guān).本文采用的溫差-風(fēng)速法使用40m高度的風(fēng)速,當(dāng)40m處風(fēng)速為3.4~7.9m/s時如30~100m存在較弱的逆溫則有可能被記做E類穩(wěn)定度,但40m高度風(fēng)速較大時,地面可能難以維持靜小風(fēng)(1.5m/s)狀態(tài),大氣擴散條件可能并不差,因而相應(yīng)的PM2.5月均值和污染日發(fā)生頻率并不高(圖8方框內(nèi)所示),計算方法可能夸大了穩(wěn)定類層結(jié)的統(tǒng)計結(jié)果,這可能是造成PM2.5均值、污染日發(fā)生頻率與穩(wěn)定類(E+F)發(fā)生頻率相關(guān)性較弱的主要原因.鄧雪嬌等[37]對廣州地區(qū)清潔過程和污染過程污染物與穩(wěn)定度Z/L的相關(guān)性分析表明,清潔過程二者存在正相關(guān),污染期間則為一定的負相關(guān),是某種程度上造成污染物與穩(wěn)定度弱相關(guān)性的重要原因.這需要結(jié)合具體污染過程進行下一步深入分析.

    圖8 逐月PM2.5均值、污染日發(fā)生頻率與穩(wěn)定類(E+F)發(fā)生頻率的關(guān)系

    圖9 逐月PM2.5均值、污染日發(fā)生頻率與逆溫發(fā)生頻率的關(guān)系

    2.4 一次重污染過程的個例分析

    2016年12月我國中東部大部分地區(qū)霾日數(shù)整體偏多,月內(nèi)共出現(xiàn)3次大范圍霧-霾天氣過程,分別為12月2~4日、10~13日和16~21日,其中12月16~21日是2016年范圍最廣、持續(xù)時間最長、強度最強的霾天氣,袁東敏等[38]的氣象診斷分析表明,地面小風(fēng)速、高濕度以及淺薄的邊界層是促使本次重度霾發(fā)生的重要因子.從環(huán)流形勢上看,這段時間我國中東部大部分地區(qū)位于冷高壓前部的均壓場控制,地面氣壓梯度弱、風(fēng)速小,以下沉氣流為主,有利于污染物累積,環(huán)流形勢有利于霧-霾天氣的發(fā)生和維持[39].圖10給出了此次過程PM2.5質(zhì)量濃度、相對濕度、大氣能見度、大氣穩(wěn)定度和逆溫強度的逐時分布狀況,本次重污染過程大體上經(jīng)歷了污染物累積、維持和快速消散,分別記作P1,P2和P3 3個階段.

    圖10 典型重污染過程大氣穩(wěn)定度、逆溫強度和PM2.5質(zhì)量濃度的時間分布特征

    12月15日午后~18日午后PM2.5逐漸累積,至18日12:00PM2.5質(zhì)量濃度達到383μg/m3,大氣能見度也由10km左右下降至1km左右,相對濕度具有明顯的日變化特征,且逐日升高,這一階段可以判斷為霾的增強階段,記作P1階段.大氣穩(wěn)定度以較穩(wěn)定(E)為主,逆溫強度也逐日增大,穩(wěn)定層結(jié)和逆溫發(fā)展對污染物的聚集起到了重要作用,PM2.5質(zhì)量濃度達到第一個峰值(383μg/m3,18日12:00),這同時也是霾形成和增強的階段.

    其后隨著RH長時間維持在90%以上,受到霧的濕清除作用[40],PM2.5濃度有所降低,維持在200~ 300μg/m3間,能見度維持在1km以下,最低時僅為0.03km,接近儀器測量極限,表現(xiàn)出典型的霧天氣特征,記作P2階段.受地面輻射降溫影響,近地面氣象條件有利于形成逆溫層結(jié),逆溫強度突然大幅度升高,最高時達到7.69℃/100m,強逆溫嚴重抑制了污染物和水汽的垂直擴散,導(dǎo)致PM2.5質(zhì)量濃度出現(xiàn)了又一個峰值(361μg/m3,19日16:00),能見度隨著RH的小幅下降短暫恢復(fù)到1km左右,霧天氣下大氣層結(jié)以中性為主[22],隨著霧趨于成熟,大氣穩(wěn)定度以中性為主,不穩(wěn)定層結(jié)時有發(fā)生,同時伴有逆溫發(fā)生頻率降低,逆溫強度減弱.吳彬貴等[21]的研究表明,華北地形槽和地面弱高壓場下的輻射霧中淺薄逆溫層夜生日消,強度從下而上迅速減弱,霧形成后逆溫底抬升到霧頂高度上下,霧體中逐漸演變呈弱不穩(wěn)定層結(jié),霧形成后逆溫層底抬升,霧體中逐漸演變呈不穩(wěn)定層結(jié).逆溫層整體抬升,超出了氣象鐵塔觀測高度,可能是本次重污染過程中霧成熟階段逆溫頻率降低的主要原因.這一階段霧-霾天氣交替發(fā)展.

    22日白天起受冷空氣影響,隨著地面偏北風(fēng)逐漸加大,氣污染擴散條件逐漸增強,大霧快速消散,能見度好轉(zhuǎn),PM2.5濃度迅速下降,從嚴重污染減弱至優(yōu)良級別,本次重污染天氣結(jié)束,記作P3階段.這一階段大氣穩(wěn)定度以不穩(wěn)定狀態(tài)為主,逆溫頻率降低,逆溫強度顯著減弱,溫度層結(jié)有利于污染物擴散.

    本次重污染過程中大氣穩(wěn)定度和逆溫特征具有明顯的階段性,大氣持續(xù)趨于穩(wěn)定及逆溫強度的逐漸增大,對污染生成和維持起了非常重要的作用,對其的精細分析有助于厘清重污染過程中的大氣擴散條件和邊界層結(jié)構(gòu)特征的演變規(guī)律,對于提升精細化重污染天氣預(yù)報預(yù)警水平也有所裨益.

    3 結(jié)論

    3.1 采用溫差-風(fēng)速法計算得到的2015年9月~ 2017年8月大氣穩(wěn)定度資料顯示,A~F類大氣穩(wěn)定度的發(fā)生頻率依次為6.7%,11.4%,22.4%,46.1%, 11.1%和2.2%,秋冬季節(jié)穩(wěn)定類天氣(E,F類)較多,全年白天各時段均以不穩(wěn)定類大氣為主,夜間大氣穩(wěn)定度以中性為主,秋冬季穩(wěn)定類高達30~40%.

    3.2 觀測期內(nèi)冬季逆溫頻率最高,5:00~8:00和21:00~23:00逆溫頻率超過90%,春、夏、秋三季則總體差別不大,冬季逆溫強度最高,秋季次之,春季和夏季最低.

    3.3 觀測期間天津全市PM2.5質(zhì)量濃度均值為(71±60)μg/m3, PM2.5質(zhì)量濃度月均值和污染日峰值均出現(xiàn)在12月.隨著穩(wěn)定類大氣層結(jié)日數(shù)的增多,PM2.5均值和污染日數(shù)均有所增大,同時逐月PM2.5均值、污染日發(fā)生頻率與逆溫發(fā)生頻率均呈正相關(guān)關(guān)系;2016年12月16~21日的一次重污染天氣過程顯示PM2.5質(zhì)量濃度受到大氣穩(wěn)定度和逆溫發(fā)展的影響,霾形成、霧-霾交替和消散等階段大氣穩(wěn)定度和逆溫特征具有顯著的不同,大氣持續(xù)趨于穩(wěn)定及逆溫強度的逐漸增大,對污染生成和維持起了非常重要的作用.

    [1] Wang Y S, Yao L, Wang L L, et al. Mechanism for the formation of the January 2013 heavy haze pollution episode over central and eastern China [J]. Science China: Earth Sciences, 2014,57:14-25.

    [2] 張人禾,李 強,張若楠.2013年1月中國東部持續(xù)性強霧霾天氣產(chǎn)生的氣象條件分析[J]. 中國科學(xué):地球科學(xué), 2014,44(1):27-36.

    [3] 楊 欣,陳義珍,劉厚鳳,等.北京2013年1月連續(xù)強霾過程的污染特征及成因分析[J]. 中國環(huán)境科學(xué), 2014,34(2):282-288.

    [4] 花 叢,張碧輝,張恒德.2013年1~2月華北霧、霾天氣邊界層特征對比分析[J]. 氣象, 2015,41(9):1144-1151.

    [5] Han S Q, Wu J H, Zhang Y F, et al. Characteristics and formation mechanism of a winter haze-fog episode in Tianjin, China [J]. Atmospheric Enviroment, 2014,98:323-330.

    [6] 姚 青,蔡子穎,韓素芹,等.天津冬季相對濕度對氣溶膠濃度譜分布和大氣能見度的影響[J]. 中國環(huán)境科學(xué), 2014,34(3):596-603.

    [7] 王 躍,王莉莉,趙廣娜,等.北京冬季 PM2.5重污染時段不同尺度環(huán)流形勢及邊界層結(jié)構(gòu)分析[J]. 氣候與環(huán)境研究, 2014,19(2):173– 84.

    [8] 孫兆斌,李梓銘,廖曉農(nóng),等.北京大氣熱力和動力結(jié)構(gòu)對污染物輸送和擴散條件的影響[J]. 中國環(huán)境科學(xué), 2017,37(5):1693-1705.

    [9] 孫永亮,趙天良,邱玉珺,等.成都一次霾過程中顆粒物消光作用的垂直變化[J]. 中國環(huán)境科學(xué), 2018,38(5):1629-1636.

    [10] 王自發(fā),李 杰,王 哲,等.2013年1月我國中東部強霾污染的數(shù)值模擬和防控對策[J]. 中國科學(xué):地球科學(xué), 2014,44(1):3-14.

    [11] 楊雨靈,譚吉華,孫家仁,等.華北地區(qū)一次強灰霾污染的天氣學(xué)效應(yīng) [J]. 氣候與環(huán)境研究, 2015,20(5):555?570.

    [12] 程良曉,范 萌,陳良富,等.秋季秸稈焚燒對京津冀地區(qū)污染過程的影響分析[J]. 中國環(huán)境科學(xué), 2017,37(8):2801-2812.

    [13] 張晗宇,溫 維,程水源,等.京津冀區(qū)域典型重污染過程與反饋效應(yīng)研究[J]. 中國環(huán)境科學(xué), 2018,38(4):1209-1220.

    [14] 康 凌,張宏升,王志遠,等.不同下墊面大氣穩(wěn)定度分類方法的對比研究[J]. 北京大學(xué)學(xué)報(自然科學(xué)版), 2011,47(1):66-70.

    [15] 金莉莉,何 清,李振杰,等.烏魯木齊市低層大氣穩(wěn)定度分布特征的統(tǒng)計分析[J]. 氣象, 2016,42(8):996-1006.

    [16] 楊 靜,李 霞,李 秦,等.烏魯木齊近30a大氣穩(wěn)定度和混合層高度變化特征及與空氣污染的關(guān)系[J]. 干旱區(qū)地理, 2011,34(5):747-752.

    [17] 楊孝文,周 穎,程水源,等.北京冬季一次重污染過程的污染特征及成因分析[J]. 中國環(huán)境科學(xué), 2016,36(3):679-686.

    [18] 龐 博,汪喜江.哈爾濱市冬末逆溫、大氣穩(wěn)定度與大氣污染概率關(guān)系[J]. 黑龍江氣象, 2011,28(2):7-11.

    [19] 解以揚,劉學(xué)軍.天津氣象塔風(fēng)溫梯度觀測資料的統(tǒng)計特征[J]. 氣象, 2003,29(1):12-16.

    [20] 韓素芹,劉彬賢,解以揚,等.利用255m鐵塔研究城市化對地面粗糙度的影響[J]. 氣象, 2008,34(1):54-58.

    [21] 吳彬貴,張宏升,張長春,等.天津城區(qū)大霧過程近地面層特征研究[J]. 氣候與環(huán)境研究, 2010,15(2):179-190.

    [22] 蔡子穎,韓素芹,吳彬貴,等.天津一次霧過程的邊界層特征研究[J]. 氣象, 2012,38(9):1103-1109.

    [23] 姚 青,蔡子穎,韓素芹,等.一次沙塵過程對天津氣溶膠濃度分布的影響[J]. 中國沙漠, 2013,33(4):1138-1143.

    [24] 姚 青,蔡子穎,韓素芹,等.2009年秋冬季天津低能見度天氣下氣溶膠污染特征[J]. 氣象, 2012,38(9):1096-1102.

    [25] 高文康,唐貴謙,姚 青,等.天津重污染期間大氣污染物濃度垂直分布特征[J]. 環(huán)境科學(xué)研究, 2012,25(7):731-738.

    [26] 黃 鶴,姚 青,張文煜,等.天津夏季低層大氣O3,NO2濃度垂直觀測研究[J]. 環(huán)境科學(xué)研究, 2009,22(3):315-320.

    [27] Hao T Y, Han S Q, Chen S C, et al. The role of fog in haze episode in Tianjin, China:A case study for November 2015. Atmospheric Research. 2017,194:235-244.

    [28] Wu H, Zhang Y F, Han S Q, et al. Vertical characteristics of PM2.5during the heating season in Tianjin, China [J]. Science of the Total Environment, 2015,523:152-160.

    [29] Han S Q, Cai Z Y, Zhang Y F, et al. Long-term trends in fog and boundary layer characteristics in Tianjin, China [J]. Particuology, 2015,20:61-68.

    [30] Han S Q, Zhang Y F, Wu J, et al. Evaluation of regional background particulate matter concentration based on vertical distribution characteristics. Atmos. Chem. Phys. 2015,15:11165–11177.

    [31] 金莉莉,何 清,李振杰,等.烏魯木齊市氣象塔梯度觀測資料質(zhì)量控制與處理[J]. 氣象, 2016,42(6):732-742.

    [32] 蔡子穎,韓素芹,張 敏,等.天津地區(qū)污染天氣分析中垂直擴散指標(biāo)構(gòu)建及運用 [J]. 環(huán)境科學(xué), 2018,39(6):2548-2556.

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

    [34] 劉麗麗,王莉莉.天津冬季重霾污染過程及氣象和邊界層特征分析 [J]. 氣候與環(huán)境研究, 2015,20(2):129-140.

    [35] 胡 寧,孫 軍. 2016年3月大氣環(huán)流和天氣分析[J]. 氣象, 2016, 42(6):770-776.

    [36] 何永晴,趙玉成,馮蜀青.西寧市逆溫特征及其與空氣污染物濃度的相關(guān)分析[J]. 冰川凍土, 2014,36(3):608-613.

    [37] 鄧雪嬌,李 菲,吳 兌,等.廣州地區(qū)典型清潔與污染過程的大氣湍流與物質(zhì)交換特征[J]. 中國環(huán)境科學(xué), 2011,31(9):1424-1430.

    [38] 袁東敏,馬小會.2016年12月16~21日重度霾過程及大氣環(huán)流異常[J]. 氣候與環(huán)境研究, 2017,22(6):757?764.

    [39] 劉 超,馬學(xué)款.2016年12月大氣環(huán)流和天氣分析[J]. 氣象, 2017, 43(3):378-384.

    [40] 康漢青,朱 彬,樊曙先.南京北郊冬季大氣氣溶膠及其濕清除特征研究 [J]. 氣候與環(huán)境研究, 2009,14(5):523-530.

    Atmospheric stability and inverse temperature characteristics of Tianjin and its effect on PM2.5pollution.

    YAO Qing1*, LIU Jin-le2, CAI Zi-ying1, HAN Su-qin2, MA Zhi-qiang3

    (1.Tianjin Environmental Meteorology Center, Tianjin 300074, China;2.Tianjin Institute of Meteorological Science, Tianjin 300074, China;3.Environmental Meteorology Forecast Center of Beijing-Tianjin-Hebei, Beijing 100089, China)., 2018,38(8):2865~2872

    Temperature-velocity method was used to calculate the atmospheric stability in Tianjin and the frequency and intensity of the inversion of temperature were analysed based on the meteorological tower. In addition, we investigated the relationship between atmospheric stability and the temperature inversion, which had an influence on the monthly averaged mass concentration of PM2.5. The characteristics of atmospheric stability and temperature inversion were analysed during a heavy haze period. The results showed that the frequencies of atmospheric stability were 6.7%, 11.4%, 22.4%, 46.1%, 11.1% and 2.2%, respectively, from September 2015 to August 2017. The stagnant meteorological conditions occurred frequently during fall and winter. The unstable and neutral atmospheric stability dominated the daytime and nighttime of all year, respectively. The stable atmosphere contributed 30%~40% during nighttime in fall and winter. During the observation period, the frequency of temperature inversion was the highest in winter, which up to more than 90% during 5:00~8:00and 21:00~23:00. The intensity of temperature inversion was also the highest in winter. With the increase of the number of days of stable atmosphere, the monthly averaged mass concentration of PM2.5and the number of polluted days increased. Meanwhile, the monthly averaged mass concentration of PM2.5and the frequency of polluted day were all positively correlated with the frequency of inversion temperature. A heavy pollution process from December 16, 2016 to December 21, 2016 showed that the mass concentration of PM2.5was effected by the atmospheric stability and the evolution of inversion temperature. The characteristics of atmospheric stability and temperature inversion during the formation of haze, fog-haze alternation and dissipation of haze had obvious difference. As the atmosphere continues to stabilize and the inversion temperature increases, it plays a very important role in the generation and maintenance of pollution. Fine analysis of atmospheric stability and inversion characteristics in the process of pollution will help to improve the prediction and warning level of heavy pollution weather.

    atmospheric stability;temperature inversion;heavily polluted weather;PM2.5;Tianjin

    X131.1

    A

    1000-6923(2018)08-2865-08

    姚 青(1980-),男,湖北宜昌人,高級工程師,碩士,主要從事大氣環(huán)境與大氣化學(xué)方面的科研業(yè)務(wù)工作.發(fā)表論文50余篇.

    2018-01-02

    國家重點研發(fā)計劃課題(2016YFC0203302);國家自然科學(xué)基金資助項目(41771242);國家國際科技合作專項項目(2015DFA20870-02);中國氣象局預(yù)報員專項(CMAYBY2016-005);中國氣象局第二批青年英才培養(yǎng)項目

    * 責(zé)任作者, 高級工程師, yao.qing@163.com

    猜你喜歡
    逆溫穩(wěn)定度均值
    小諸葛漫“話”天氣之逆溫
    濟南市低空逆溫特征及其對空氣質(zhì)量影響分析
    酒泉市大氣邊界層逆溫特征及其與沙塵天氣的關(guān)系
    高穩(wěn)晶振短期頻率穩(wěn)定度的仿真分析
    2011—2013年長春市逆溫特征分析
    均值不等式失效時的解決方法
    多MOSFET并聯(lián)均流的高穩(wěn)定度恒流源研究
    均值與方差在生活中的應(yīng)用
    工藝參數(shù)對橡膠球鉸徑向剛度穩(wěn)定度的影響
    關(guān)于均值有界變差函數(shù)的重要不等式
    亚洲精品成人av观看孕妇| 国产精品久久视频播放| 女同久久另类99精品国产91| 岛国在线观看网站| www日本在线高清视频| 中文字幕精品免费在线观看视频| 亚洲 欧美一区二区三区| 国产精品一区二区在线不卡| 最近最新中文字幕大全电影3 | 99re6热这里在线精品视频| 亚洲专区中文字幕在线| ponron亚洲| 老汉色∧v一级毛片| 亚洲专区中文字幕在线| 精品电影一区二区在线| 国产在线一区二区三区精| 人人妻,人人澡人人爽秒播| 国产日韩欧美亚洲二区| 欧美日韩精品网址| 不卡一级毛片| 人人妻,人人澡人人爽秒播| 国产在线观看jvid| 亚洲国产看品久久| 国产精品免费一区二区三区在线 | 久久国产精品大桥未久av| 中文字幕人妻丝袜制服| 少妇裸体淫交视频免费看高清 | 国产精品免费视频内射| 午夜两性在线视频| 丰满饥渴人妻一区二区三| 精品电影一区二区在线| 国产精品久久久av美女十八| 日韩 欧美 亚洲 中文字幕| 国产在线一区二区三区精| 国产深夜福利视频在线观看| 少妇 在线观看| 欧美国产精品va在线观看不卡| 天天躁狠狠躁夜夜躁狠狠躁| 欧美大码av| 两人在一起打扑克的视频| 久久精品91无色码中文字幕| 亚洲精品美女久久av网站| 他把我摸到了高潮在线观看| 国产xxxxx性猛交| 丝袜在线中文字幕| 一边摸一边做爽爽视频免费| 在线十欧美十亚洲十日本专区| 啦啦啦视频在线资源免费观看| 国产精品欧美亚洲77777| 国产高清videossex| av电影中文网址| 女人爽到高潮嗷嗷叫在线视频| 久久精品国产综合久久久| 一边摸一边抽搐一进一出视频| 亚洲精品一卡2卡三卡4卡5卡| 水蜜桃什么品种好| 建设人人有责人人尽责人人享有的| 久久久久久人人人人人| 色在线成人网| 亚洲国产看品久久| 国产91精品成人一区二区三区| 亚洲精品美女久久av网站| 视频区欧美日本亚洲| 免费观看人在逋| 亚洲一区二区三区欧美精品| 亚洲色图 男人天堂 中文字幕| 国产精品秋霞免费鲁丝片| 亚洲成a人片在线一区二区| 亚洲精品国产一区二区精华液| 欧美一级毛片孕妇| 国产精品免费大片| 欧美另类亚洲清纯唯美| 国产精品免费大片| 精品少妇久久久久久888优播| 丁香六月欧美| 国产精品 欧美亚洲| 国产精品 欧美亚洲| 中文字幕最新亚洲高清| 亚洲午夜理论影院| 色婷婷久久久亚洲欧美| 欧美亚洲日本最大视频资源| 久久国产乱子伦精品免费另类| 亚洲成人免费av在线播放| 国产精品国产av在线观看| 91麻豆精品激情在线观看国产 | 精品国产一区二区三区四区第35| 天堂动漫精品| 国产有黄有色有爽视频| 51午夜福利影视在线观看| 最近最新免费中文字幕在线| 成年人黄色毛片网站| 亚洲精品乱久久久久久| 男女高潮啪啪啪动态图| 日日夜夜操网爽| 国产精品久久久久成人av| 欧美亚洲日本最大视频资源| 在线视频色国产色| 国产精品久久久人人做人人爽| 国产一卡二卡三卡精品| 亚洲精品国产精品久久久不卡| 午夜亚洲福利在线播放| av网站免费在线观看视频| 久久香蕉激情| 亚洲国产中文字幕在线视频| 大陆偷拍与自拍| 久久亚洲精品不卡| 午夜福利,免费看| 国产男女超爽视频在线观看| 国产有黄有色有爽视频| 嫁个100分男人电影在线观看| 美女 人体艺术 gogo| 国产男女内射视频| 精品第一国产精品| 国精品久久久久久国模美| 亚洲午夜理论影院| 老司机在亚洲福利影院| 国产精品1区2区在线观看. | 啦啦啦 在线观看视频| 女人久久www免费人成看片| 99精国产麻豆久久婷婷| 母亲3免费完整高清在线观看| 两性夫妻黄色片| 国产97色在线日韩免费| 久久久国产一区二区| 久久国产精品人妻蜜桃| 久久精品国产综合久久久| 国产精品免费视频内射| 18禁观看日本| 美女午夜性视频免费| 一个人免费在线观看的高清视频| 交换朋友夫妻互换小说| 精品无人区乱码1区二区| 在线播放国产精品三级| 欧美乱码精品一区二区三区| 热re99久久精品国产66热6| tube8黄色片| 国产一卡二卡三卡精品| 国产免费现黄频在线看| 亚洲av日韩在线播放| 波多野结衣一区麻豆| 法律面前人人平等表现在哪些方面| 中文字幕高清在线视频| 一边摸一边做爽爽视频免费| 日本欧美视频一区| 午夜福利在线免费观看网站| 高清黄色对白视频在线免费看| 成人黄色视频免费在线看| av网站在线播放免费| 午夜日韩欧美国产| 久久久久久免费高清国产稀缺| 国产真人三级小视频在线观看| a级片在线免费高清观看视频| 久久中文字幕一级| 青草久久国产| 美女扒开内裤让男人捅视频| 午夜日韩欧美国产| 女人被躁到高潮嗷嗷叫费观| 久久亚洲精品不卡| 俄罗斯特黄特色一大片| 欧美日韩视频精品一区| 亚洲国产中文字幕在线视频| 中文字幕精品免费在线观看视频| 亚洲精品在线观看二区| 国产一区二区三区综合在线观看| 99riav亚洲国产免费| 亚洲精品国产一区二区精华液| 国精品久久久久久国模美| 国产精品一区二区精品视频观看| 午夜福利视频在线观看免费| 岛国毛片在线播放| 91国产中文字幕| 午夜成年电影在线免费观看| 精品卡一卡二卡四卡免费| 999精品在线视频| a级片在线免费高清观看视频| 久久精品国产亚洲av高清一级| 日本vs欧美在线观看视频| 国产成人免费无遮挡视频| 亚洲精品中文字幕在线视频| 国产又爽黄色视频| 日韩视频一区二区在线观看| 成年人黄色毛片网站| 亚洲精品国产一区二区精华液| 王馨瑶露胸无遮挡在线观看| 亚洲成人免费电影在线观看| 亚洲第一av免费看| 人人妻人人澡人人爽人人夜夜| 欧美乱妇无乱码| 国产精品永久免费网站| 欧美性长视频在线观看| 乱人伦中国视频| 日本黄色日本黄色录像| 男女之事视频高清在线观看| bbb黄色大片| 日韩欧美在线二视频 | 天堂动漫精品| 久久久久久久久久久久大奶| 久久精品国产99精品国产亚洲性色 | 热99re8久久精品国产| 亚洲欧美日韩高清在线视频| 亚洲欧美精品综合一区二区三区| 在线视频色国产色| 精品国产一区二区三区久久久樱花| 中文字幕人妻熟女乱码| 久久人妻av系列| 成年人午夜在线观看视频| 国产欧美日韩一区二区三区在线| 国产一卡二卡三卡精品| 中亚洲国语对白在线视频| e午夜精品久久久久久久| 9191精品国产免费久久| 久久人人97超碰香蕉20202| 黄色视频不卡| 桃红色精品国产亚洲av| 伦理电影免费视频| 国精品久久久久久国模美| 色婷婷久久久亚洲欧美| 男人舔女人的私密视频| 在线视频色国产色| 久久久精品区二区三区| 麻豆乱淫一区二区| 午夜福利一区二区在线看| 久久久久久久精品吃奶| 99久久综合精品五月天人人| 久久精品亚洲精品国产色婷小说| 国产深夜福利视频在线观看| 国产精品九九99| 国产熟女午夜一区二区三区| 国产精品一区二区免费欧美| 亚洲精品成人av观看孕妇| 少妇被粗大的猛进出69影院| 欧美精品人与动牲交sv欧美| 老司机福利观看| 高清欧美精品videossex| 97人妻天天添夜夜摸| avwww免费| 国产成人精品久久二区二区91| 亚洲免费av在线视频| 丁香六月欧美| 国产欧美日韩一区二区三区在线| 美女高潮喷水抽搐中文字幕| 桃红色精品国产亚洲av| 国产高清国产精品国产三级| 欧美不卡视频在线免费观看 | 国产免费男女视频| 国产在线观看jvid| 日韩欧美免费精品| 首页视频小说图片口味搜索| 无限看片的www在线观看| 在线观看一区二区三区激情| 男女下面插进去视频免费观看| 国产亚洲欧美精品永久| 激情在线观看视频在线高清 | 中文欧美无线码| 人人澡人人妻人| ponron亚洲| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲aⅴ乱码一区二区在线播放 | 69av精品久久久久久| 国产精品免费大片| 久热爱精品视频在线9| 他把我摸到了高潮在线观看| 又大又爽又粗| 中文字幕最新亚洲高清| 久久久久视频综合| 国产精品久久久久成人av| 国产xxxxx性猛交| 大型av网站在线播放| 精品午夜福利视频在线观看一区| 两性午夜刺激爽爽歪歪视频在线观看 | a级毛片在线看网站| 国产激情久久老熟女| 美女国产高潮福利片在线看| 一区福利在线观看| 国产成人影院久久av| 成人18禁高潮啪啪吃奶动态图| 欧美日韩av久久| 另类亚洲欧美激情| 久久精品成人免费网站| 成人影院久久| 亚洲欧美精品综合一区二区三区| 91麻豆精品激情在线观看国产 | 国产极品粉嫩免费观看在线| 啦啦啦视频在线资源免费观看| 国产野战对白在线观看| 亚洲一码二码三码区别大吗| 法律面前人人平等表现在哪些方面| 中文亚洲av片在线观看爽 | 久久ye,这里只有精品| 精品一区二区三区视频在线观看免费 | 欧美精品av麻豆av| 黄色成人免费大全| 久久中文看片网| 波多野结衣一区麻豆| 9色porny在线观看| 精品一区二区三区视频在线观看免费 | 啦啦啦视频在线资源免费观看| 久久人人97超碰香蕉20202| 黄色视频,在线免费观看| 老司机午夜福利在线观看视频| 日日夜夜操网爽| 精品亚洲成a人片在线观看| 亚洲人成77777在线视频| 中出人妻视频一区二区| 亚洲精品一卡2卡三卡4卡5卡| 在线免费观看的www视频| 一区在线观看完整版| 一级毛片高清免费大全| av中文乱码字幕在线| 高清黄色对白视频在线免费看| 国产精品一区二区在线不卡| 午夜视频精品福利| 午夜福利在线免费观看网站| 久久草成人影院| 一级黄色大片毛片| 亚洲精品一二三| 如日韩欧美国产精品一区二区三区| 啦啦啦在线免费观看视频4| 大型黄色视频在线免费观看| 精品国产一区二区三区久久久樱花| 美女扒开内裤让男人捅视频| 婷婷成人精品国产| 免费高清在线观看日韩| 熟女少妇亚洲综合色aaa.| 亚洲五月色婷婷综合| 亚洲熟妇中文字幕五十中出 | 18在线观看网站| 亚洲情色 制服丝袜| 亚洲全国av大片| 最近最新中文字幕大全免费视频| 亚洲色图 男人天堂 中文字幕| 精品少妇久久久久久888优播| 人人妻,人人澡人人爽秒播| cao死你这个sao货| 国产精品免费大片| 热re99久久精品国产66热6| 国产精品一区二区在线观看99| 亚洲精品乱久久久久久| 777久久人妻少妇嫩草av网站| 国产一卡二卡三卡精品| 一进一出抽搐动态| 欧美激情久久久久久爽电影 | 久久国产精品人妻蜜桃| 好男人电影高清在线观看| 国产成人精品无人区| 精品第一国产精品| 国产无遮挡羞羞视频在线观看| 成在线人永久免费视频| 高清av免费在线| 国产成人欧美| 欧美日韩乱码在线| 欧美成人免费av一区二区三区 | 一进一出好大好爽视频| 欧美日韩一级在线毛片| 久久精品亚洲熟妇少妇任你| 热99久久久久精品小说推荐| 午夜精品久久久久久毛片777| 欧美人与性动交α欧美软件| 80岁老熟妇乱子伦牲交| 在线观看免费午夜福利视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲欧美色中文字幕在线| 91av网站免费观看| 国产成人av教育| 成人精品一区二区免费| 日韩欧美国产一区二区入口| 国产av精品麻豆| 成人三级做爰电影| 欧美日韩精品网址| 久久中文看片网| bbb黄色大片| 王馨瑶露胸无遮挡在线观看| 两个人免费观看高清视频| 国产欧美日韩一区二区精品| 热re99久久国产66热| 丁香欧美五月| av网站免费在线观看视频| 国产精品 欧美亚洲| 亚洲五月色婷婷综合| 啦啦啦 在线观看视频| 极品人妻少妇av视频| 不卡av一区二区三区| 久久久精品免费免费高清| 亚洲一区二区三区欧美精品| 亚洲在线自拍视频| 国产日韩欧美亚洲二区| 免费看十八禁软件| 欧美乱码精品一区二区三区| 在线永久观看黄色视频| 国产欧美日韩综合在线一区二区| 亚洲一区中文字幕在线| 欧美日韩中文字幕国产精品一区二区三区 | 久久香蕉国产精品| 国产国语露脸激情在线看| 曰老女人黄片| 天堂动漫精品| 亚洲精品国产一区二区精华液| 成人av一区二区三区在线看| а√天堂www在线а√下载 | 99精品欧美一区二区三区四区| 首页视频小说图片口味搜索| 亚洲欧美激情在线| 久久国产精品大桥未久av| 久久人妻av系列| 在线观看一区二区三区激情| 欧美日韩一级在线毛片| 亚洲欧美精品综合一区二区三区| 俄罗斯特黄特色一大片| 黄色a级毛片大全视频| 搡老岳熟女国产| av福利片在线| 手机成人av网站| 亚洲性夜色夜夜综合| 国产高清国产精品国产三级| 美女午夜性视频免费| 女警被强在线播放| 建设人人有责人人尽责人人享有的| 满18在线观看网站| 亚洲欧美日韩高清在线视频| 国产淫语在线视频| 国产成人av教育| 999久久久精品免费观看国产| 国产精品 国内视频| 精品卡一卡二卡四卡免费| 精品人妻熟女毛片av久久网站| 精品久久蜜臀av无| 三上悠亚av全集在线观看| 亚洲人成电影观看| 高清欧美精品videossex| 亚洲男人天堂网一区| 熟女少妇亚洲综合色aaa.| 嫩草影视91久久| 一夜夜www| 黑人操中国人逼视频| 亚洲av美国av| 一区在线观看完整版| 精品欧美一区二区三区在线| 男人舔女人的私密视频| 欧美日韩成人在线一区二区| 欧美 亚洲 国产 日韩一| 成人亚洲精品一区在线观看| 搡老乐熟女国产| 午夜福利乱码中文字幕| 日本vs欧美在线观看视频| 国产99白浆流出| 最近最新免费中文字幕在线| 国产精品久久久久成人av| 亚洲av日韩精品久久久久久密| 国产高清国产精品国产三级| 欧美av亚洲av综合av国产av| 一进一出抽搐动态| 国产亚洲一区二区精品| av网站免费在线观看视频| 热99久久久久精品小说推荐| 可以免费在线观看a视频的电影网站| 国产成人欧美| 亚洲五月婷婷丁香| 欧美日韩av久久| 午夜亚洲福利在线播放| 亚洲成人手机| 一本一本久久a久久精品综合妖精| www.精华液| 啦啦啦免费观看视频1| 在线av久久热| av线在线观看网站| 波多野结衣一区麻豆| 久久香蕉精品热| 99国产精品一区二区蜜桃av | 精品国产乱码久久久久久男人| 91精品国产国语对白视频| 捣出白浆h1v1| 亚洲在线自拍视频| 免费在线观看完整版高清| 国产精品 欧美亚洲| 国产精品香港三级国产av潘金莲| 狂野欧美激情性xxxx| 精品人妻1区二区| 无人区码免费观看不卡| 欧美亚洲日本最大视频资源| 18在线观看网站| 亚洲精品成人av观看孕妇| 69av精品久久久久久| 99国产精品一区二区三区| 露出奶头的视频| 黄色毛片三级朝国网站| 亚洲av美国av| 精品第一国产精品| 91麻豆精品激情在线观看国产 | 最新美女视频免费是黄的| 午夜福利一区二区在线看| 久久久国产一区二区| 亚洲国产毛片av蜜桃av| 精品久久久久久电影网| 三上悠亚av全集在线观看| 一级作爱视频免费观看| 久久香蕉精品热| 1024视频免费在线观看| 9热在线视频观看99| 大型av网站在线播放| 欧美性长视频在线观看| 黄频高清免费视频| 欧美乱码精品一区二区三区| 极品教师在线免费播放| 亚洲精品国产一区二区精华液| av中文乱码字幕在线| 免费在线观看黄色视频的| 悠悠久久av| 99香蕉大伊视频| 久久久久久人人人人人| 欧美日韩福利视频一区二区| 成人av一区二区三区在线看| 久久香蕉激情| 国产无遮挡羞羞视频在线观看| 高清毛片免费观看视频网站 | 午夜日韩欧美国产| 90打野战视频偷拍视频| 亚洲av片天天在线观看| 露出奶头的视频| 精品高清国产在线一区| 精品人妻熟女毛片av久久网站| 50天的宝宝边吃奶边哭怎么回事| 久久精品国产清高在天天线| √禁漫天堂资源中文www| 欧美精品高潮呻吟av久久| 久久久久视频综合| 久久国产乱子伦精品免费另类| aaaaa片日本免费| 日本五十路高清| 女人精品久久久久毛片| 国产亚洲精品第一综合不卡| 亚洲一区高清亚洲精品| 91成人精品电影| 亚洲第一欧美日韩一区二区三区| 亚洲国产精品sss在线观看 | 嫁个100分男人电影在线观看| 亚洲五月色婷婷综合| 嫁个100分男人电影在线观看| 麻豆成人av在线观看| 最新的欧美精品一区二区| 侵犯人妻中文字幕一二三四区| 欧美精品人与动牲交sv欧美| 国产乱人伦免费视频| 亚洲精品自拍成人| 久热爱精品视频在线9| 一边摸一边抽搐一进一小说 | 丝袜人妻中文字幕| 久久久精品区二区三区| 国产av又大| 成人av一区二区三区在线看| 操出白浆在线播放| 精品人妻1区二区| 亚洲精品乱久久久久久| 日本欧美视频一区| 亚洲精品国产一区二区精华液| 亚洲精品久久成人aⅴ小说| 国内久久婷婷六月综合欲色啪| 欧美老熟妇乱子伦牲交| 黄片播放在线免费| 亚洲中文av在线| 精品卡一卡二卡四卡免费| 极品教师在线免费播放| 欧美人与性动交α欧美软件| 男女下面插进去视频免费观看| 老熟妇乱子伦视频在线观看| 久久天躁狠狠躁夜夜2o2o| 动漫黄色视频在线观看| 欧美色视频一区免费| 亚洲 国产 在线| 亚洲人成伊人成综合网2020| 日韩成人在线观看一区二区三区| 色综合欧美亚洲国产小说| 久久久精品免费免费高清| 亚洲情色 制服丝袜| 国产精品九九99| 美国免费a级毛片| 久久久久久久精品吃奶| 麻豆成人av在线观看| 久久久国产成人免费| 日韩制服丝袜自拍偷拍| 老司机午夜福利在线观看视频| 精品一品国产午夜福利视频| 国内久久婷婷六月综合欲色啪| 欧美黄色片欧美黄色片| 91精品国产国语对白视频| 叶爱在线成人免费视频播放| 不卡av一区二区三区| 久久天躁狠狠躁夜夜2o2o| 超碰97精品在线观看| 国产精品美女特级片免费视频播放器 | 人人妻,人人澡人人爽秒播| 久久精品国产99精品国产亚洲性色 | 亚洲,欧美精品.| 丁香六月欧美| 国产精品一区二区在线观看99| 久久精品亚洲熟妇少妇任你| 母亲3免费完整高清在线观看| 成人特级黄色片久久久久久久| 99精品久久久久人妻精品| 涩涩av久久男人的天堂| 美女扒开内裤让男人捅视频| 亚洲男人天堂网一区| 99国产精品99久久久久| 51午夜福利影视在线观看| 一边摸一边抽搐一进一出视频| 亚洲人成电影观看| 成人国产一区最新在线观看| 久久久久久亚洲精品国产蜜桃av| 午夜福利乱码中文字幕| 亚洲精品国产色婷婷电影| 国产一区二区三区综合在线观看| 99久久99久久久精品蜜桃| 交换朋友夫妻互换小说| 国产亚洲av高清不卡| 人人妻人人澡人人看| 国产一区在线观看成人免费|