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

    上甸子本底站鹵代溫室氣體大氣濃度短期波動(dòng)的個(gè)例分析

    2014-09-09 07:30:15馬林安興琴朱彬周凌晞姚波車勝利李東宇
    大氣科學(xué)學(xué)報(bào) 2014年4期
    關(guān)鍵詞:大氣區(qū)域

    馬林,安興琴,朱彬,周凌晞,姚波,車勝利,李東宇

    (1.遼寧省氣象裝備保障中心,遼寧 沈陽(yáng) 110166;2.南京信息工程大學(xué) 大氣物理學(xué)院,江蘇 南京 210044;3.中國(guó)氣象科學(xué)研究院,北京 100081)

    上甸子本底站鹵代溫室氣體大氣濃度短期波動(dòng)的個(gè)例分析

    馬林1,2,安興琴3,朱彬2,周凌晞3,姚波3,車勝利1,李東宇1

    (1.遼寧省氣象裝備保障中心,遼寧 沈陽(yáng) 110166;2.南京信息工程大學(xué) 大氣物理學(xué)院,江蘇 南京 210044;3.中國(guó)氣象科學(xué)研究院,北京 100081)

    利用軌跡分析法、印痕分析和流場(chǎng)分析法,結(jié)合上甸子站鹵代溫室氣體H-1301、HCFC-22、CFC-11和SF6在線濃度觀測(cè)數(shù)據(jù),選取2012年9月7—12日上甸子測(cè)站鹵代溫室氣體濃度短期波動(dòng)典型個(gè)例進(jìn)行分析。軌跡分析結(jié)果表明:7日12時(shí),污染發(fā)生前,氣團(tuán)主要來(lái)自較遠(yuǎn)的偏西北、偏北方向,水平輸送距離長(zhǎng),移動(dòng)迅速,垂直高度高,對(duì)應(yīng)的鹵代溫室氣體濃度偏低,H-1301、HCFC-22、CFC-11和SF6的體積分?jǐn)?shù)分別為4×10-12、350×10-12、260×10-12、10×10-12;9、10日有一定比例的氣團(tuán)在測(cè)站的偏南區(qū)域近地面回旋打轉(zhuǎn),水平輸送距離短,垂直高度低,在邊界層內(nèi)緩慢移動(dòng),不利于污染物在邊界層內(nèi)擴(kuò)散,導(dǎo)致鹵代溫室氣體濃度偏高,對(duì)測(cè)站濃度的短期抬升貢獻(xiàn)較大,9日12時(shí)H-1301、HCFC-22、CFC-11的峰值體積分?jǐn)?shù)分別達(dá)到45×10-12、1 200×10-12、310×10-12,10日03時(shí)SF6的峰值體積分?jǐn)?shù)達(dá)到28×10-12;11日西南方向回旋氣團(tuán)消失;12日氣團(tuán)完全來(lái)自較遠(yuǎn)的西北方向且軌跡移動(dòng)較快。印痕分析與軌跡分析結(jié)果一致:7、8日敏感性系數(shù)較高區(qū)域主要分布在測(cè)站以北,9、10日敏感性系數(shù)較高區(qū)域分布在測(cè)站偏南,11、12日測(cè)站偏南的敏感性系數(shù)較高區(qū)域消失。流場(chǎng)分析結(jié)果表明:9、10日環(huán)流形勢(shì)有利于污染物在測(cè)站區(qū)域累積,造成測(cè)站濃度的短期抬升。

    鹵代溫室氣體;上甸子區(qū)域本底站;個(gè)例分析

    0 引言

    鹵代溫室氣體是指含氟、氯、溴、碘等鹵素的溫室氣體,其化學(xué)性質(zhì)穩(wěn)定,無(wú)毒、無(wú)臭、不可燃,多用于制冷劑、推進(jìn)劑、發(fā)泡劑、噴霧劑等。由于鹵代溫室氣體會(huì)損耗臭氧層(莫天麟和蔣龍海,1984),并導(dǎo)致全球變暖,這類物質(zhì)已被列入《蒙特利爾議定書(shū)》和《京都議定書(shū)》的減排清單。根據(jù)國(guó)際公約要求,2010—2015年中國(guó)全面禁止大部分鹵代溫室氣體(姚波等,2007)。與此同時(shí),隨著經(jīng)濟(jì)活動(dòng)和工業(yè)化的迅速發(fā)展,中國(guó)已成為生產(chǎn)和使用氟氯烴的大國(guó),其排放和淘汰CFCs、HCFCs以及其他鹵代烴的情況受到全球的廣泛關(guān)注(張芳,2006)。目前在全球范圍內(nèi)已有相當(dāng)數(shù)量的鹵代溫室氣體常規(guī)觀測(cè)系統(tǒng),如美國(guó)國(guó)家海洋和大氣管理局/地球系統(tǒng)研究實(shí)驗(yàn)室(NOAA/ESRL)、改進(jìn)的全球大氣實(shí)驗(yàn)網(wǎng)(AGAGE)、歐盟鹵代溫室氣體觀測(cè)網(wǎng)(SOGE)等在南北半球的一些大氣本底站開(kāi)展了多種鹵代溫室氣體的在線觀測(cè)(姚波等,2010a)。為了解我國(guó)典型區(qū)域鹵代溫室氣體的本底濃度水平及排放分布,評(píng)估我國(guó)減排進(jìn)程和履約效果,開(kāi)展我國(guó)鹵代溫室氣體濃度的長(zhǎng)期連續(xù)監(jiān)測(cè)與研究具有重要意義(周凌晞等,2006)。姚波等(2010b)利用自組裝的氣相色譜儀—電子捕獲檢測(cè)器系統(tǒng)(GC-ECD)測(cè)定大氣中12種鹵代溫室氣體;陳立民等(1999)通過(guò)低溫脫水和低溫濃縮建立了氣相色譜(ECD檢測(cè)器)測(cè)定大氣中超痕量氯氟烴類物質(zhì)(CECs)的分析方法,對(duì)上海市清潔大氣中CFCs本底濃度進(jìn)行監(jiān)測(cè);孫楊等(2004)基于氣相色譜儀(GC)和電子捕獲檢測(cè)器(ECD),設(shè)計(jì)了一種連續(xù)測(cè)量大氣中痕量氯氟烴的自動(dòng)系統(tǒng),實(shí)現(xiàn)了對(duì)大氣中CFCs的連續(xù)測(cè)定;楊越等(2001)通過(guò)對(duì)北半球中緯度到北極圈地區(qū)大規(guī)模的氣體樣品的采集,運(yùn)用GC-ECD的檢測(cè)方法,得到了隨緯度變化的CCl3F、CCl2F和CClF2的特征濃度;張芳等(2006)采用預(yù)濃縮-GC/MS方法,研究了珠江三角洲大氣中的CCl3F、CFC-12、CFC-113和CFC-114等4種痕量氟氯烷烴氣體的濃度水平和變化特征。

    高體積分?jǐn)?shù)的污染氣體主要?dú)w因于局地排放,同時(shí)也會(huì)受到周邊地區(qū)或者遠(yuǎn)距離輸送的影響(金祺等,2012)。對(duì)鹵代溫室氣體大氣濃度短期波動(dòng)進(jìn)行個(gè)例分析,可以考察不同來(lái)向的輸送對(duì)測(cè)站鹵代溫室氣體大氣濃度變化的影響,有利于分析測(cè)站大氣鹵代溫室氣體本底濃度變化及源匯機(jī)制。本文利用2009年上甸子本底站GC-ECD觀測(cè)系統(tǒng)在線連續(xù)觀測(cè)的大氣H-1301、HCFC-22、CFC-11和SF6數(shù)據(jù),選取2012年9月7—12日鹵代溫室氣體濃度短期波動(dòng)的典型個(gè)例,采用大氣擴(kuò)散模式、后向軌跡方法,分析不同來(lái)向的輸送對(duì)上甸子測(cè)站觀測(cè)濃度的影響,探尋引起測(cè)站污染物濃度變化的氣團(tuán)源地。

    1 研究區(qū)域、數(shù)據(jù)來(lái)源及研究方法

    1.1 研究區(qū)域

    圖1 上甸子大氣本底站及周邊主要城市Fig.1 Shangdianzi regional background monitoring station and major cities around it

    上甸子區(qū)域大氣本底站是世界氣象組織/全球大氣觀測(cè)網(wǎng)的區(qū)域本底站之一(圖1),建立于1983年,測(cè)站(117°07’E,40°39’N)位于北京市東北方向,距離北京市區(qū)120 km,海拔高度為293.3 m(蒲維維等,2012)。該站地處暖溫帶半濕潤(rùn)氣候區(qū),夏季炎熱多雨,冬季寒冷干燥,春、秋季短暫,具有典型的華北區(qū)域氣候特征(孟昭陽(yáng)等,2007)。在一定程度上,上甸子本底監(jiān)測(cè)站的長(zhǎng)期觀測(cè)資料反映了京、津、唐乃至整個(gè)華北地區(qū)的大氣狀況。因此,對(duì)上甸子站鹵代溫室氣體短期波動(dòng)進(jìn)行個(gè)例分析具有重要意義。

    1.2 數(shù)據(jù)來(lái)源

    本文所用的鹵代溫室氣體觀測(cè)數(shù)據(jù)來(lái)自于上甸子測(cè)站的GC-ECD在線觀測(cè)系統(tǒng)(姚波等,2010a),該系統(tǒng)位于上甸子站實(shí)驗(yàn)室內(nèi),采樣塔(高度10 m)位于實(shí)驗(yàn)室正北約60 m處,采樣口位于采樣塔8 m高處。

    1.3 研究方法

    1.3.1 軌跡分析

    大氣擴(kuò)散模式HYSPLIT(Hybrid Single Particle Lagrangian Integrated) (http://www.arl.noaa.gov/ready/hysplit4.html)是由美國(guó)國(guó)家海洋和大氣管理局(National Oceanic and Atmospheric Administration,NOAA)以及澳大利亞氣象局(Bureau of Meteorology Australia,BOM)聯(lián)合研制的一種基于Eulerian-Lagrangian的混合計(jì)算模式,其平流和擴(kuò)散計(jì)算采用Lagrangian方法,濃度計(jì)算采用Eulerian方法(Draxler and Hess,1997)。本文利用HYSPLIT模式(Wang et al.,2009)計(jì)算了2009年9月7—12日到達(dá)上甸子測(cè)站的逐日8個(gè)時(shí)次(北京時(shí)00時(shí)、03時(shí)、06時(shí)、09時(shí)、12時(shí)、15時(shí)、18時(shí)、21時(shí),下同)120 h的后向氣流軌跡,起始高度選擇為距地500 m(劉娜等,2012)。模式所用氣象資料采用NCEP模式(National Centers for Environmental Prediction)輸出的經(jīng)全球資料同化系統(tǒng)GDAS(Global Data Assimilation System)處理的FNL全球分析資料(http://www.arl.noaa.gov/ss/transport/archives.html)(魏文華等,2012)。

    1.3.2 印痕分析

    采用FLEXPART模式模擬“印痕分布”,即“潛在源區(qū)”,表征格點(diǎn)區(qū)域?qū)y(cè)站觀測(cè)濃度的可能影響。FLEXPART模式是一個(gè)拉格朗日粒子擴(kuò)散模式(Stohl et al.,2005),主要研究大氣示蹤物的長(zhǎng)期中尺度輸送、擴(kuò)散、干濕沉降及其輻射衰減等過(guò)程,通過(guò)后向軌跡模擬確定源對(duì)特定接受點(diǎn)的可能貢獻(xiàn)。本文選取2009年9月7—12日的時(shí)間序列,利用FLEXPART模式模擬潛在源區(qū)的分布。

    1.3.3 流場(chǎng)分析

    排放到大氣中的污染物,由于垂直湍流的擴(kuò)散作用,在經(jīng)過(guò)數(shù)十千米的水平傳輸后,在邊界層內(nèi)基本混合均勻,此時(shí),邊界層頂?shù)牡剞D(zhuǎn)風(fēng)氣流方向大體代表了大氣污染物的傳輸方向,即邊界層頂?shù)牧鲌?chǎng)對(duì)污染物的輸送起主導(dǎo)作用(時(shí)進(jìn)剛等,2010)。

    利用NCEP再分析資料(沈新勇等,2012)繪制9月7—12日各典型時(shí)刻所對(duì)應(yīng)的850 hPa平均高度場(chǎng)并疊加風(fēng)場(chǎng),反映污染事件發(fā)生時(shí)所對(duì)應(yīng)的大氣環(huán)流形勢(shì)特征。

    圖2 2009年9月7—12日上甸子測(cè)站大氣H-1301(a)、HCFC-22(b)、CFC-11(c)和SF6(d)的體積分?jǐn)?shù)變化Fig.2 Volume fraction changes of (a)H-1301,(b)HCFC-22,(c)CFC-11 and (d)SF6 at Shangdianzi regional background monitoring station during 7—12 September 2009

    2 研究結(jié)果

    2.1 軌跡分析

    圖2描述的是2009年9月7—12日上甸子測(cè)站大氣H-1301、HCFC-22、CFC-11和SF6的體積分?jǐn)?shù)變化??梢?jiàn),9月7—8日,上甸子大氣H-1301、HCFC-22、CFC-11和SF6的體積分?jǐn)?shù)較低,分別維持在(4~5)×10-12、(200~350)×10-12、(240~270)×10-12、(6~10)×10-12,9—10日急劇抬升,H-1301、HCFC-22、CFC-11和SF6的體積分?jǐn)?shù)峰值分別達(dá)到45×10-12、1 200×10-12、310×10-12、28×10-12,11日開(kāi)始下降,12日回落到4×10-12、355×10-12、265×10-12、8×10-12。7—12日是上甸子測(cè)站大氣H-1301、HCFC-22、CFC-11和SF6的體積分?jǐn)?shù)短期波動(dòng)的典型個(gè)例。

    為了深入分析污染發(fā)生過(guò)程,選取7日12時(shí)(污染發(fā)生前)、9日12時(shí)(污染發(fā)生時(shí))、10日00時(shí)(污染高峰時(shí))、11日03時(shí)(污染減退時(shí))和12日00時(shí)(污染發(fā)生后)5個(gè)典型時(shí)刻,分析其后向氣流軌跡。圖3為上甸子測(cè)站2009年9月7—12日5個(gè)典型時(shí)刻的后向氣流軌跡??梢钥闯?7日12時(shí),氣團(tuán)來(lái)自較遠(yuǎn)的偏西北地區(qū),軌跡較長(zhǎng),移動(dòng)速度較快;9月12日,氣團(tuán)軌跡變短,移動(dòng)速度減慢,大部分來(lái)自較近的西北方向,但在測(cè)站的西南方向出現(xiàn)了一個(gè)小回旋;10日00時(shí),氣團(tuán)軌跡較短,移動(dòng)較慢,測(cè)站西南方向的回旋達(dá)到最大;11日03時(shí),氣團(tuán)軌跡較短,在西南方向的回旋消失;12日00時(shí),氣團(tuán)完全來(lái)自較遠(yuǎn)的西北方向,且軌跡變長(zhǎng),移動(dòng)較快。

    圖3 上甸子測(cè)站2009年9月7—12日5個(gè)典型時(shí)刻的后向軌跡(圖中數(shù)字7代表7日12時(shí)、數(shù)字9代表9日12時(shí)、數(shù)字10代表10日00時(shí)、數(shù)字11代表11日03時(shí)、數(shù)字12代表12日00時(shí))Fig.3 Backward trajectory of the five typical moments at Shangdianzi regional background monitoring station during 7—12 September 2009(figures 7,9,10,11 and 12 represent 1200 UTC 7,1200 UTC 9,0000 UTC 10,0300 UTC 11,and 0000 UTC 12,respectively)

    圖4為5個(gè)典型時(shí)刻后向軌跡的垂直變化??梢钥闯?5個(gè)時(shí)刻垂直高度在500~900 hPa之間;在到達(dá)測(cè)站的前24 h,以10日00時(shí)的高度較低,位于800~900 hPa,且軌跡較平穩(wěn),而其余時(shí)刻的高度則較高且軌跡波動(dòng)較大。

    圖4 上甸子測(cè)站2009年9月7—12日5個(gè)典型時(shí)刻后向軌跡的垂直變化Fig.4 Vertical variations of back trajectory of the five typical moments at Shangdianzi regional background monitoring station during 7—12 September 2009

    圖5 2009年9月7—12日5個(gè)典型時(shí)刻上甸子區(qū)域大氣本底站印痕分布(黑色圓點(diǎn)代表上甸子站) a.7日12時(shí);b.9日12時(shí);c.10日00時(shí);d.11日03時(shí);e.12日00時(shí)Fig.5 Footprint distribution of the five typical moments at Shangdianzi regional background monitoring station during 7—12 September 2009(the black dots denote Shangdianzi station) a.1200 UTC 7;b.1200 UTC 9;c.0000 UTC 10;d.0300 UTC 11;e.0000 UTC 12

    2.2 印痕分析

    圖5給出了利用FLEXPART模式模擬的5個(gè)典型時(shí)刻上甸子測(cè)站的印痕分布。可見(jiàn),7日12時(shí),敏感性較高的區(qū)域分布在以上甸子測(cè)站為中心的偏西北方向,覆蓋范圍較小;9日12時(shí),高值區(qū)域在西北方向擴(kuò)大,敏感性較高的區(qū)域在測(cè)站西南方開(kāi)始出現(xiàn);10月00時(shí),敏感性較高區(qū)域在測(cè)站西南方的覆蓋范圍達(dá)到峰值;11日03時(shí),敏感性較高區(qū)域在測(cè)站西南方的分布范圍開(kāi)始減小;12月00日,敏感性較高的區(qū)域在測(cè)站西南方向的覆蓋消失,只在測(cè)站北部分布。

    2.3 流場(chǎng)分析

    圖6 2009年9月7—12日5個(gè)典型時(shí)刻上甸子區(qū)域大氣本底站850 hPa風(fēng)場(chǎng)(箭矢;單位:m·s-1)和位勢(shì)高度場(chǎng)(陰影;單位gpm) a.7日12時(shí);b.9日12時(shí);c.10日00時(shí);d.11日03時(shí);e.12日00時(shí)Fig.6 850 hPa wind(arrows;units:m·s-1) and geopotential height(shadings;units:gpm) fields of the five typical moments at Shangdianzi regional background monitoring station during 7—12 September 2009 a.1200 UTC 7;b.1200 UTC 9;c.0000 UTC 10;d.0300 UTC 11;e.0000 UTC 12

    圖6為上甸子大氣本底站2009年9月7—12日850 hPa風(fēng)場(chǎng)??梢?jiàn),7日12時(shí),我國(guó)東北以東地區(qū)存在一低壓系統(tǒng),貝加爾湖以南地區(qū)為高壓脊控制,上甸子站處于低壓槽后和高壓脊前的偏北氣流中,風(fēng)速較大。9日12時(shí),貝加爾湖以南高壓向東南方向移動(dòng),污染發(fā)生時(shí)上甸子站位于南下高壓后部的偏南氣流中;同時(shí)貝加爾湖以西又有新的高壓東進(jìn),在我國(guó)內(nèi)蒙古中東部存在一個(gè)西南氣流和西北氣流的匯合區(qū)。10日00時(shí),西北氣流和西南氣流的匯合區(qū)移到北京地區(qū)上空,風(fēng)速較小,有利于污染物積累,上甸子測(cè)站污染物濃度達(dá)到峰值。11日03時(shí),貝加爾湖以西的高壓已經(jīng)推進(jìn)到貝加爾湖以南地區(qū),上甸子站上空為脊前的西北氣流控制,風(fēng)速較大,有利于污染物擴(kuò)散,污染開(kāi)始消退。12日00時(shí),貝加爾湖以南的高壓主體位于我國(guó)內(nèi)蒙古中西部以及山西和陜西等地,我國(guó)華北地區(qū)上空為脊前偏北氣流控制,污染過(guò)程結(jié)束。

    2.4 三種方法的結(jié)果比較

    應(yīng)該上述3種方法得到如下相同結(jié)論:上甸子測(cè)站7、8日氣團(tuán)的來(lái)向主要為偏西北方向,垂直高度較高,印痕結(jié)果表明敏感性系數(shù)較高區(qū)域分布以測(cè)站北向?yàn)橹?9、10日雖然也有來(lái)向西北向氣團(tuán),但在測(cè)站西南方向氣團(tuán)達(dá)到頂峰,且移動(dòng)速度慢,垂直高度相對(duì)低,基本在近地面緩慢移動(dòng),敏感性系數(shù)較高區(qū)域在測(cè)站偏南方向覆蓋面積最大;11、12日南向氣團(tuán)逐漸消失,氣團(tuán)完全來(lái)自于偏北向,移動(dòng)較快,垂直高度較高,敏感性系數(shù)較高區(qū)域在測(cè)站偏南部分布消失,全部分布在測(cè)站偏北向。從污染發(fā)生過(guò)程的5個(gè)典型時(shí)刻來(lái)看,7日12時(shí),污染發(fā)生前,上甸子站處于低壓槽后和高壓脊前的偏北氣流中,風(fēng)速較大,氣團(tuán)來(lái)自較遠(yuǎn)的偏西北地區(qū),移動(dòng)速度較快,敏感性較高的潛在源區(qū)以上甸子本底站為中心向偏西北方向分布,覆蓋范圍較小;9月12日,我國(guó)內(nèi)蒙古中東部存在一個(gè)西南氣流和西北氣流的匯合區(qū),上甸子站位于南下高壓后部的偏南氣流中,氣團(tuán)雖大部分來(lái)自于較近的西北方向,但在測(cè)站的西南方向開(kāi)始出現(xiàn)氣團(tuán),且移動(dòng)速度減慢,敏感性較高的區(qū)域在測(cè)站的西南方開(kāi)始出現(xiàn);10日00時(shí),西北氣流和西南氣流的匯合區(qū)移到北京地區(qū)上空,風(fēng)速較小,測(cè)站西南方向的氣團(tuán)達(dá)到最大,氣團(tuán)移動(dòng)較慢,敏感性較高的區(qū)域在測(cè)站西南方的覆蓋范圍達(dá)到峰值;11日03時(shí),上甸子站上空為脊前的西北氣流控制,風(fēng)速較大,來(lái)自測(cè)站西南方向氣團(tuán)開(kāi)始減少,敏感性較高的區(qū)域在測(cè)站西南方的分布范圍開(kāi)始減小;12日00時(shí),我國(guó)華北地區(qū)上空為脊前偏北氣流控制,氣團(tuán)完全來(lái)自較遠(yuǎn)的西北方向,移動(dòng)速度較快,敏感性較高區(qū)域在測(cè)站西南方向的覆蓋消失,只在測(cè)站北部分布。

    圖7 2008年?yáng)|亞地區(qū)HCFC-22排放量分布(單位:10-12 g/(m2·s))(引自Stohl et al.(2010))Fig.7 The emission intensity of HCFC-22 over East Asia in 2008(units:10-12 g/(m2·s))(from Stohl et al.(2010))

    此外,根據(jù)Stohl et al.(2010)對(duì)東亞2008年HCFC-22排放量的研究結(jié)果(圖7),再結(jié)合本文研究結(jié)論,可知上甸子測(cè)站北部地區(qū)的HCFC-22排放量很低,由于鹵代溫室氣體來(lái)自人為源排放,同源性好,HCFC-22的排放水平能大體能代表其他鹵代溫室氣體的排放水平,所以測(cè)站北部是人類活動(dòng)較少區(qū)域,來(lái)自偏北方向的氣團(tuán)對(duì)上甸子測(cè)站污染的貢獻(xiàn)較小。而上甸子測(cè)站的南部地區(qū)分布著北京、天津、石家莊、唐山等人口密集、工業(yè)化發(fā)達(dá)的大中城市,途經(jīng)該地區(qū)的偏南氣團(tuán)會(huì)攜帶高濃度污染物到上甸子站,造成測(cè)站污染物濃度升高。

    3 結(jié)論與討論

    2012年9月7—12日期間的典型污染個(gè)例分析表明:上甸子測(cè)站的西北區(qū)域?qū)儆谌丝谙鄬?duì)較少的非城市工業(yè)區(qū),鹵代物種的人為排放低,污染發(fā)生前,氣團(tuán)大部分來(lái)自測(cè)站偏北、偏西北方向,移動(dòng)快速,垂直混合及下沉運(yùn)動(dòng)對(duì)測(cè)站的采樣濃度影響較小,對(duì)應(yīng)的H-1301、HCFC-22、CFC-11和SF6濃度較低。9、10日出現(xiàn)了一定比例的偏南方向氣團(tuán),在測(cè)站的偏南區(qū)域近地面回旋打轉(zhuǎn),水平軌跡短,垂直高度低,水平移動(dòng)和垂直運(yùn)動(dòng)小,在近低層緩慢移動(dòng),有利于污染物在邊界層內(nèi)累積;加上上甸子測(cè)站的偏東南區(qū)域是人口密集、高排放的城市工業(yè)區(qū),分布著一些排放合成鹵代物種的企業(yè)。因此,在此區(qū)域近地層回旋打轉(zhuǎn)的偏南氣團(tuán)將高排放的污染氣團(tuán)攜帶到上甸子測(cè)站,造成這類軌跡對(duì)應(yīng)的H-1301、HCFC-22、CFC-11和SF6的濃度一致增高,與偏北軌跡輸送相比,對(duì)上甸子本底站鹵代溫室氣體濃度水平短期抬升的貢獻(xiàn)較大。FLEXPART模式的印痕結(jié)果與軌跡分析結(jié)果一致:7、8日污染發(fā)生前,敏感性系數(shù)較高區(qū)域主要分布在測(cè)站以北,9、10日污染發(fā)生時(shí),敏感性系數(shù)較高區(qū)域分布在測(cè)站偏南位置,11、12日污染消散時(shí),測(cè)站偏南分布的敏感性系數(shù)較高區(qū)域消失。850 hPa流場(chǎng)的分析結(jié)果更深入揭示了造成這次污染個(gè)例的環(huán)流形勢(shì),9、10日測(cè)站區(qū)域的偏南氣流和輻合下沉氣流,有利于污染物的累積,造成上甸子測(cè)站濃度的短期抬升。

    致謝:周凌晞研究員的溫室氣體團(tuán)隊(duì)在資料獲取、論文撰寫(xiě)過(guò)程中給予幫助,挪威大氣研究所Andre博士提供FLEXPART模式,上甸子測(cè)站工作人員采集、提供數(shù)據(jù),在此一并感謝!

    陳立民,段楊,樂(lè)志威,等.1999.大氣中氯氟烴類物質(zhì)濃度變化的研究[J].環(huán)境科學(xué),20(1):27-29.

    金祺,銀燕,譚穩(wěn).2012.黃山地區(qū)復(fù)雜地形下污染氣體輸送過(guò)程的數(shù)值模擬[J].大氣科學(xué)學(xué)報(bào),35(6):680-688.

    劉娜,余曄,陳晉北,等.2012.蘭州春季沙塵過(guò)程PM10輸送路徑及其潛在源區(qū)[J].大學(xué)科學(xué)學(xué)報(bào),35(4):477-486.

    孟昭陽(yáng),丁國(guó)安,湯潔,等.2007.北京上甸子本底站2003年秋冬季痕量氣體濃度變化特征[J].氣象科技,35(4):551-552.

    莫天麟,蔣龍海.1984.大氣中鹵代化合物對(duì)氣候的影響[J].氣象科學(xué),4(1):37-43.

    蒲維維,于波,趙秀娟,等.2012.北京地區(qū)酸雨的天氣影響因素及降水化學(xué)特征分析[J].氣象科學(xué),32(5):565-572.

    沈新勇,劉佳,秦南南,等.2012.臺(tái)風(fēng)麥莎的正壓特征波動(dòng)結(jié)構(gòu)及其穩(wěn)定性[J].大氣科學(xué)學(xué)報(bào),35(3):257-271.

    時(shí)進(jìn)剛,王亞男,陳曉麗.2010.中尺度氣象流場(chǎng)特征及其對(duì)大氣污染物跨界環(huán)境影響的作用——以黃河中上游能源化工區(qū)為例[C]//中國(guó)環(huán)境科學(xué)學(xué)會(huì)環(huán)境影響評(píng)價(jià)專業(yè)委員會(huì)2010年學(xué)術(shù)年會(huì)暨國(guó)家環(huán)境保護(hù)環(huán)境影響評(píng)價(jià)重點(diǎn)實(shí)驗(yàn)室(籌)2010年學(xué)術(shù)年會(huì):43-47.

    孫楊,王躍思,盛國(guó)英,等.2004.改進(jìn)GC/ECD法連續(xù)測(cè)定大氣中的CFCs[J].環(huán)境污染治理技術(shù)與設(shè)備,5(8):37-39.

    魏文華,王體健,石春娥,等.2012.合肥市霧日氣象條件分析[J].氣象科學(xué),32(4):437-442.

    楊越,陳立民,孫國(guó)忠,等.2001.北半球大氣中CCl3F和CCl2FCClF2隨緯度變化的特征和趨勢(shì)研究[J].復(fù)旦學(xué)報(bào):自然科學(xué)版,40(4):55-60.

    姚波,周凌晞,Vollmer M K,等.2007.北京上甸子區(qū)域本底站鹵代溫室氣體在線測(cè)[C]//中國(guó)氣象學(xué)會(huì)年會(huì)2007年年會(huì)大氣成分觀測(cè)、研究與預(yù)報(bào)分會(huì)場(chǎng)論文集.廣州:中國(guó)氣象學(xué)會(huì).

    姚波,周凌晞,張芳,等.2010a.北京上甸子區(qū)域大氣本底站HCFC-22在線觀測(cè)研究[J].環(huán)境科學(xué),31(8):1749-1754.

    姚波,周凌晞,張芳,等.2010b.氣相色譜—電子捕獲檢測(cè)器法在線觀測(cè)12種鹵代溫室氣體[J].環(huán)境化學(xué),29(6):1158-1162.

    張芳.2006.珠江三角洲地區(qū)大氣中與全球變化相關(guān)的痕量鹵代烴的初步研究[D].廣州:中國(guó)科學(xué)院廣州地球化學(xué)研究所.

    張芳,王新明,李龍鳳,等.2006.近年來(lái)珠三角地區(qū)大氣中痕量氟氯烴(CFCs)的濃度水平與變化特征[J].地球與環(huán)境,34(4):19-24.

    周凌晞,張曉春,郝慶菊,等.2006.溫室氣體本底觀測(cè)研究[J].氣候變化研究進(jìn)展,2(2):63-67.

    Draxler R R,Hess G D.1997.Description of the HYSPLIT24 modeling system[C]//NOAA Technical memorandum ERL ARL:220-224.

    Stohl A,Forster C,Frank A,et al.2005.Technical Note:The Lagrangian particle dispersion model FLEXPART version 6.2[J].Atmospheric Chemistry and Physics,5:2461-2474.

    Stohl A,Kim J,Li S,et al.2010.Hydrochlorofluorocarbon and hydrofluorocarbon emissions in East Asia determined by inverse modeling[J].Atmos Chem Phys,10:3545-3560.

    Wang Y Q,Zhang X Y,Draxler R.2009.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,24(8):938-939.

    (責(zé)任編輯:倪東鴻)

    Acasestudyofshort-termfluctuationinatmosphericconcentrationofhalogenatedgreenhousegasesatShangdianziregionalbackgroundmonitoringstation

    MA Lin1,2,AN Xing-qin3,ZHU Bin2,ZHOU Ling-xi3,YAO Bo3,CHE Sheng-li1,LI Dong-yu1

    (1.Liaoning Meteorological Equipment Support Center,Shenyang 110166,China;2.School of Atmospheric Physics,NUIST,Nanjing 210044,China;3.Chinese Academy of Meteorological Sciences,Beijing 100081,China)

    Based on the online observed concentration of H-1301,HCFC-22,CFC-11 and SF6 at Shangdianzi background station during 7—12 September 2012,a typical case about short-term fluctuation of halogenated greenhouse gases is analyzed by means of trajectory,footprint and flow field analysis methods.The trajectory analysis results show that,before occurrence of pollution,air masses mainly come from the further north and northwest at 1200 UTC 7 September 2012,with a long horizontal transport distance,rapid moving speed and high vertical height.Meanwhile,the corresponding volume fractions of halogenated greenhouse gases are relatively low,which are 4×10-12,350×10-12,260×10-12and 10×10-12for H-1301,HCFC-22,CFC-11 and SF6,respectively.However,in 9—10 September 2012,a certain percentage of air masses convolute over the south of the station with a short horizontal transport distance and low vertical height,moving slowly in the boundary layer,which hinders the spread of pollutants in the boundary layer and results in the relatively high concentrations of halogenated greenhouse gases.Such situation contributes a lot to the short-term rapid growth of concentration levels.The peak volume fractions of H-1301,HCFC-22 and CFC-11 reach 45×10-12,1 200×10-12and 310×10-12at 1200 UTC 9 September,respectively,and that of SF6 does 28×10-12at 0300 UTC 10 September.On 11 September,the convoluting air mass in the southwest direction disappears.On 12 September,air masses all come from the far northwest and have rapid moving trajectories.The footprint analysis results are the same as the trajectory analysis results,which are that the regions with high sensitivity coefficient are mainly located in the north of the station in 7—8 September,and in the south of the station in 9—10 September.Then the regions with high sensitivity coefficients in the south of the station disappear in 11—12 September.In addition,the flow field analysis results show that the circulation pattern in 9—10 September is favorable for the accumulation of pollutants in the observed region,resulting in the short-term lift of pollutants’ concentrations.

    halogenated greenhouse gases;Shangdianzi regional background monitoring station;case study

    2012-03-12;改回日期2014-03-25

    國(guó)家自然科學(xué)基金資助項(xiàng)目(41030107);國(guó)際科技合作計(jì)劃項(xiàng)目(1015);中國(guó)氣象科學(xué)研究院基本科研業(yè)務(wù)費(fèi)專項(xiàng)(2010Y003)

    安興琴,博士,研究員,研究方向?yàn)閰^(qū)域和城市大氣污染物分布和傳輸規(guī)律的數(shù)值模擬、污染控制和污染源評(píng)估,anxq@cams.cma.gov.cn.

    10.13878/j.cnki.dqkxxb.20120312002.

    1674-7097(2014)04-0493-09

    P402

    A

    10.13878/j.cnki.dqkxxb.20120312002

    馬林,安興琴,朱彬,等.2014.上甸子本底站鹵代溫室氣體大氣濃度短期波動(dòng)的個(gè)例分析[J].大氣科學(xué)學(xué)報(bào),37(4):493-501.

    Ma Lin,An Xing-qin,Zhu Bin,et al.2014.A case study of short-term fluctuation in atmospheric concentration of halogenated greenhouse gases at Shangdianzi regional background monitoring station[J].Trans Atmos Sci,37(4):493-501.(in Chinese)

    猜你喜歡
    大氣區(qū)域
    大氣的呵護(hù)
    軍事文摘(2023年10期)2023-06-09 09:15:06
    永久基本農(nóng)田集中區(qū)域“禁廢”
    太赫茲大氣臨邊探測(cè)儀遙感中高層大氣風(fēng)仿真
    分割區(qū)域
    大氣古樸揮灑自如
    大氣、水之后,土十條來(lái)了
    關(guān)于四色猜想
    分區(qū)域
    基于嚴(yán)重區(qū)域的多PCC點(diǎn)暫降頻次估計(jì)
    區(qū)域
    民生周刊(2012年10期)2012-10-14 09:06:46
    一区福利在线观看| 女人被狂操c到高潮| 黄色丝袜av网址大全| 成人永久免费在线观看视频| 欧美日本亚洲视频在线播放| 极品教师在线免费播放| 在线免费观看不下载黄p国产 | 亚洲成人久久爱视频| 一本精品99久久精品77| 亚洲中文字幕日韩| 亚洲国产精品sss在线观看| 国产国拍精品亚洲av在线观看 | 国产亚洲精品av在线| 人妻夜夜爽99麻豆av| 国产蜜桃级精品一区二区三区| 搡老妇女老女人老熟妇| 女人被狂操c到高潮| 9191精品国产免费久久| 国产精品久久久久久亚洲av鲁大| 最近视频中文字幕2019在线8| 一级作爱视频免费观看| 久久久久久久亚洲中文字幕 | 黄色女人牲交| 午夜福利在线观看吧| 熟女少妇亚洲综合色aaa.| 无限看片的www在线观看| 成人国产一区最新在线观看| 天天添夜夜摸| 亚洲一区高清亚洲精品| 内射极品少妇av片p| 女人高潮潮喷娇喘18禁视频| 久久久国产精品麻豆| 麻豆国产av国片精品| 两个人看的免费小视频| 我的老师免费观看完整版| 国产91精品成人一区二区三区| 母亲3免费完整高清在线观看| 特大巨黑吊av在线直播| 一个人看视频在线观看www免费 | 国产色婷婷99| 动漫黄色视频在线观看| 精品乱码久久久久久99久播| 桃红色精品国产亚洲av| 国产 一区 欧美 日韩| av福利片在线观看| 久久婷婷人人爽人人干人人爱| 岛国在线免费视频观看| 色尼玛亚洲综合影院| 99久久九九国产精品国产免费| 色噜噜av男人的天堂激情| 一夜夜www| 欧美丝袜亚洲另类 | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国内精品美女久久久久久| www.色视频.com| 中亚洲国语对白在线视频| 一卡2卡三卡四卡精品乱码亚洲| 伊人久久精品亚洲午夜| 午夜福利欧美成人| 欧美黑人巨大hd| 日韩欧美在线二视频| 成年版毛片免费区| 九九久久精品国产亚洲av麻豆| 日韩av在线大香蕉| 高清在线国产一区| 国产午夜福利久久久久久| 亚洲av免费在线观看| 黄片大片在线免费观看| 国产成人av激情在线播放| 精品一区二区三区视频在线观看免费| 女生性感内裤真人,穿戴方法视频| 1024手机看黄色片| 免费观看精品视频网站| 国产av在哪里看| 黄色视频,在线免费观看| 日本撒尿小便嘘嘘汇集6| 午夜免费男女啪啪视频观看 | 亚洲精品456在线播放app | 波多野结衣巨乳人妻| 日本免费一区二区三区高清不卡| 叶爱在线成人免费视频播放| 亚洲国产中文字幕在线视频| av天堂中文字幕网| 久久久久久久久久黄片| 好男人在线观看高清免费视频| 国内精品久久久久精免费| av在线蜜桃| 久久精品影院6| 一级作爱视频免费观看| 国产精品99久久久久久久久| 国产伦精品一区二区三区视频9 | 免费av不卡在线播放| 99久久精品热视频| 国产淫片久久久久久久久 | 搞女人的毛片| 久久久久久久久大av| 18禁裸乳无遮挡免费网站照片| 亚洲天堂国产精品一区在线| 国产极品精品免费视频能看的| 天堂动漫精品| e午夜精品久久久久久久| 国产精品亚洲美女久久久| 国产午夜精品久久久久久一区二区三区 | 国产亚洲av嫩草精品影院| 亚洲在线自拍视频| 午夜免费成人在线视频| 中亚洲国语对白在线视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 此物有八面人人有两片| 欧美激情在线99| avwww免费| 18美女黄网站色大片免费观看| 九色国产91popny在线| 亚洲精品在线美女| 国产一区二区三区在线臀色熟女| 亚洲国产中文字幕在线视频| 国产精品久久久久久久电影 | 午夜免费男女啪啪视频观看 | 亚洲国产精品久久男人天堂| 国产免费一级a男人的天堂| 国产一区二区在线av高清观看| 三级毛片av免费| 香蕉久久夜色| 久久中文看片网| 好男人在线观看高清免费视频| 一区二区三区国产精品乱码| 色噜噜av男人的天堂激情| 国模一区二区三区四区视频| 麻豆一二三区av精品| 波多野结衣高清无吗| 欧美一区二区亚洲| 国产久久久一区二区三区| 久久精品91蜜桃| 麻豆国产97在线/欧美| 国产日本99.免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 欧美黄色片欧美黄色片| 久久久久性生活片| 三级男女做爰猛烈吃奶摸视频| 99精品欧美一区二区三区四区| 午夜精品在线福利| 无限看片的www在线观看| 国产精品香港三级国产av潘金莲| 99精品久久久久人妻精品| 日韩欧美在线乱码| 国产精品久久久人人做人人爽| 一本一本综合久久| 日韩国内少妇激情av| 九九久久精品国产亚洲av麻豆| 一本一本综合久久| 欧美日韩瑟瑟在线播放| 免费电影在线观看免费观看| 最新在线观看一区二区三区| 日本黄色视频三级网站网址| 90打野战视频偷拍视频| 高清在线国产一区| 又粗又爽又猛毛片免费看| 少妇裸体淫交视频免费看高清| 成人18禁在线播放| 日韩欧美国产一区二区入口| 午夜福利在线观看吧| 99精品在免费线老司机午夜| 日韩成人在线观看一区二区三区| 日韩欧美精品免费久久 | 免费在线观看影片大全网站| 日韩av在线大香蕉| 人人妻人人看人人澡| 一级a爱片免费观看的视频| 免费看十八禁软件| 无限看片的www在线观看| 日韩欧美一区二区三区在线观看| 在线免费观看不下载黄p国产 | 一区二区三区国产精品乱码| 91麻豆精品激情在线观看国产| 欧美日韩福利视频一区二区| 免费无遮挡裸体视频| 日韩欧美三级三区| www国产在线视频色| 亚洲人成网站在线播放欧美日韩| 我要搜黄色片| 啦啦啦韩国在线观看视频| 国产真实伦视频高清在线观看 | 少妇裸体淫交视频免费看高清| 在线播放国产精品三级| 1000部很黄的大片| 性色av乱码一区二区三区2| 国产在视频线在精品| 日韩欧美 国产精品| 观看美女的网站| 中文在线观看免费www的网站| 又粗又爽又猛毛片免费看| 丝袜美腿在线中文| 18+在线观看网站| 亚洲片人在线观看| 在线免费观看不下载黄p国产 | 国产亚洲精品av在线| 亚洲美女视频黄频| 两人在一起打扑克的视频| 国产精品98久久久久久宅男小说| 国产一区二区三区在线臀色熟女| 淫妇啪啪啪对白视频| 无遮挡黄片免费观看| 在线观看午夜福利视频| 伊人久久精品亚洲午夜| 亚洲无线在线观看| 国产精品一区二区免费欧美| 一a级毛片在线观看| 午夜精品在线福利| 黄色日韩在线| 亚洲av一区综合| 男女那种视频在线观看| 国产不卡一卡二| tocl精华| 黄色日韩在线| 搡老妇女老女人老熟妇| 成人一区二区视频在线观看| 色在线成人网| 一进一出抽搐动态| 日韩中文字幕欧美一区二区| av视频在线观看入口| 国产激情欧美一区二区| 精品国产美女av久久久久小说| 18禁黄网站禁片午夜丰满| 日韩精品中文字幕看吧| 亚洲性夜色夜夜综合| 国产在视频线在精品| 精品人妻1区二区| 在线观看午夜福利视频| 久久久久国产精品人妻aⅴ院| 国产高清视频在线观看网站| 成人18禁在线播放| 啦啦啦免费观看视频1| 久9热在线精品视频| 岛国在线观看网站| 亚洲av电影在线进入| 午夜久久久久精精品| 国产三级在线视频| 国产欧美日韩一区二区精品| 国产精品亚洲一级av第二区| 国产高清有码在线观看视频| 中文资源天堂在线| 日本精品一区二区三区蜜桃| 亚洲不卡免费看| av天堂在线播放| 十八禁人妻一区二区| 久久精品亚洲精品国产色婷小说| 欧美日韩中文字幕国产精品一区二区三区| 国内久久婷婷六月综合欲色啪| 此物有八面人人有两片| 欧美色视频一区免费| 亚洲最大成人中文| 听说在线观看完整版免费高清| 国产高清videossex| 国语自产精品视频在线第100页| 岛国视频午夜一区免费看| 国产成人aa在线观看| av专区在线播放| 国产高清激情床上av| 一级黄片播放器| 国产熟女xx| 欧美bdsm另类| 国产97色在线日韩免费| 久久精品人妻少妇| 麻豆国产av国片精品| 成年版毛片免费区| 国产高清有码在线观看视频| 在线国产一区二区在线| 深夜精品福利| 国产亚洲欧美98| 1000部很黄的大片| 毛片女人毛片| 精品不卡国产一区二区三区| 丝袜美腿在线中文| 国产精品乱码一区二三区的特点| 亚洲国产欧美人成| 人妻丰满熟妇av一区二区三区| 午夜福利欧美成人| 99久久综合精品五月天人人| 日韩欧美精品v在线| 精品熟女少妇八av免费久了| 欧美三级亚洲精品| 亚洲精品久久国产高清桃花| 国产精品国产高清国产av| 成人特级av手机在线观看| 天堂影院成人在线观看| 男女下面进入的视频免费午夜| aaaaa片日本免费| 欧美日韩黄片免| 亚洲aⅴ乱码一区二区在线播放| 在线国产一区二区在线| 久久久久久久午夜电影| 国产成人a区在线观看| 久久亚洲精品不卡| 色av中文字幕| 人人妻,人人澡人人爽秒播| 午夜福利视频1000在线观看| 99视频精品全部免费 在线| 69人妻影院| 国产亚洲精品综合一区在线观看| 国产一区二区亚洲精品在线观看| 淫妇啪啪啪对白视频| 婷婷六月久久综合丁香| 精品99又大又爽又粗少妇毛片 | 三级国产精品欧美在线观看| 午夜免费激情av| 少妇熟女aⅴ在线视频| 国产精品爽爽va在线观看网站| av黄色大香蕉| av欧美777| 性欧美人与动物交配| 欧美一区二区精品小视频在线| 男插女下体视频免费在线播放| 免费无遮挡裸体视频| 最近视频中文字幕2019在线8| 欧美一级毛片孕妇| 国产毛片a区久久久久| 美女大奶头视频| 一个人看视频在线观看www免费 | 美女黄网站色视频| 最近最新中文字幕大全免费视频| 天堂网av新在线| 在线视频色国产色| 欧美在线一区亚洲| 亚洲欧美日韩高清在线视频| 999久久久精品免费观看国产| 亚洲av第一区精品v没综合| 日本在线视频免费播放| 久久久精品大字幕| 午夜影院日韩av| 欧美xxxx黑人xx丫x性爽| 成人18禁在线播放| 我的老师免费观看完整版| 国产精品一区二区三区四区久久| 欧美色欧美亚洲另类二区| 亚洲色图av天堂| 哪里可以看免费的av片| 成年女人毛片免费观看观看9| 黄色日韩在线| 五月玫瑰六月丁香| 国产日本99.免费观看| av福利片在线观看| 成人特级av手机在线观看| 亚洲国产中文字幕在线视频| 国产高潮美女av| 99国产精品一区二区蜜桃av| 怎么达到女性高潮| 一级a爱片免费观看的视频| 日本一二三区视频观看| 久久精品亚洲精品国产色婷小说| 少妇丰满av| 99久国产av精品| 在线视频色国产色| 亚洲av不卡在线观看| 久久国产乱子伦精品免费另类| aaaaa片日本免费| 久久欧美精品欧美久久欧美| 国产亚洲欧美98| 女人高潮潮喷娇喘18禁视频| 性欧美人与动物交配| 亚洲美女黄片视频| 成年女人毛片免费观看观看9| 亚洲自拍偷在线| 99久久精品国产亚洲精品| 嫩草影视91久久| 操出白浆在线播放| 亚洲精品在线观看二区| 国产日本99.免费观看| 久久久久亚洲av毛片大全| 国产精品美女特级片免费视频播放器| 日韩精品青青久久久久久| 国产美女午夜福利| 人人妻,人人澡人人爽秒播| 一级黄片播放器| 波多野结衣高清无吗| 国产美女午夜福利| 亚洲精品一区av在线观看| 97超级碰碰碰精品色视频在线观看| 午夜免费成人在线视频| 狂野欧美激情性xxxx| 亚洲美女黄片视频| 成人亚洲精品av一区二区| av国产免费在线观看| 国产精品久久久久久久久免 | 色噜噜av男人的天堂激情| 熟女少妇亚洲综合色aaa.| 国内揄拍国产精品人妻在线| 国内久久婷婷六月综合欲色啪| 色综合亚洲欧美另类图片| 日韩欧美国产在线观看| 精品乱码久久久久久99久播| 啦啦啦韩国在线观看视频| 国产一区二区在线av高清观看| 亚洲av中文字字幕乱码综合| 精品久久久久久久毛片微露脸| 精品久久久久久久人妻蜜臀av| 有码 亚洲区| 成人精品一区二区免费| 亚洲狠狠婷婷综合久久图片| 两性午夜刺激爽爽歪歪视频在线观看| 老鸭窝网址在线观看| av天堂中文字幕网| 国产亚洲av嫩草精品影院| 精品一区二区三区人妻视频| 热99re8久久精品国产| 人人妻,人人澡人人爽秒播| 桃色一区二区三区在线观看| 欧美在线黄色| 国内久久婷婷六月综合欲色啪| 18禁国产床啪视频网站| 中文字幕久久专区| 亚洲无线在线观看| 日韩国内少妇激情av| 女人十人毛片免费观看3o分钟| 亚洲最大成人中文| 久久天躁狠狠躁夜夜2o2o| 女同久久另类99精品国产91| 亚洲成av人片免费观看| 两个人的视频大全免费| 亚洲欧美激情综合另类| 18禁黄网站禁片免费观看直播| 国产精品久久电影中文字幕| 久久婷婷人人爽人人干人人爱| 最新美女视频免费是黄的| 亚洲成人中文字幕在线播放| 99精品久久久久人妻精品| 此物有八面人人有两片| 小蜜桃在线观看免费完整版高清| 黄色女人牲交| 高清在线国产一区| 99在线人妻在线中文字幕| 成人精品一区二区免费| 国产黄色小视频在线观看| 一个人看视频在线观看www免费 | 午夜免费男女啪啪视频观看 | 国产久久久一区二区三区| 91麻豆精品激情在线观看国产| 国产成人av激情在线播放| 亚洲精品国产精品久久久不卡| netflix在线观看网站| 一夜夜www| 精品福利观看| 国产亚洲av嫩草精品影院| 国产爱豆传媒在线观看| 又爽又黄无遮挡网站| 国产真人三级小视频在线观看| 蜜桃久久精品国产亚洲av| 国产私拍福利视频在线观看| 国产激情欧美一区二区| 我要搜黄色片| 一级黄色大片毛片| av女优亚洲男人天堂| 男人舔奶头视频| 久久草成人影院| 亚洲片人在线观看| 天堂网av新在线| 亚洲无线在线观看| 国产精品国产高清国产av| 99riav亚洲国产免费| 麻豆久久精品国产亚洲av| 狂野欧美白嫩少妇大欣赏| 男女午夜视频在线观看| 亚洲国产精品久久男人天堂| АⅤ资源中文在线天堂| 国产成+人综合+亚洲专区| 欧美zozozo另类| 成人特级黄色片久久久久久久| 熟女人妻精品中文字幕| 美女免费视频网站| 国产aⅴ精品一区二区三区波| 久久精品影院6| 91九色精品人成在线观看| 亚洲黑人精品在线| 日本黄大片高清| 午夜激情福利司机影院| 欧美bdsm另类| 人妻丰满熟妇av一区二区三区| 久久久精品欧美日韩精品| 成人国产综合亚洲| 好男人在线观看高清免费视频| av天堂在线播放| 亚洲国产中文字幕在线视频| 长腿黑丝高跟| 久久久久久久久中文| 最近最新中文字幕大全免费视频| 亚洲电影在线观看av| 国产成人av教育| eeuss影院久久| 国产综合懂色| 舔av片在线| 国语自产精品视频在线第100页| 真人一进一出gif抽搐免费| 日韩欧美国产一区二区入口| 欧洲精品卡2卡3卡4卡5卡区| 国产成人福利小说| 搞女人的毛片| 久久性视频一级片| 日本与韩国留学比较| 日本在线视频免费播放| 女生性感内裤真人,穿戴方法视频| 欧美绝顶高潮抽搐喷水| 亚洲五月婷婷丁香| 2021天堂中文幕一二区在线观| 蜜桃久久精品国产亚洲av| 长腿黑丝高跟| 少妇熟女aⅴ在线视频| 久久性视频一级片| 亚洲美女黄片视频| 真实男女啪啪啪动态图| 噜噜噜噜噜久久久久久91| 亚洲av五月六月丁香网| 亚洲欧美日韩高清专用| 亚洲精品日韩av片在线观看 | 国产不卡一卡二| 夜夜躁狠狠躁天天躁| 叶爱在线成人免费视频播放| 麻豆成人av在线观看| 欧美日韩福利视频一区二区| 色噜噜av男人的天堂激情| 人妻丰满熟妇av一区二区三区| 免费观看人在逋| 日本黄色视频三级网站网址| 国产69精品久久久久777片| 午夜免费观看网址| 欧美日本视频| 性色avwww在线观看| 午夜福利在线在线| 精品久久久久久久久久久久久| 色综合站精品国产| h日本视频在线播放| 操出白浆在线播放| 精品国产三级普通话版| 看片在线看免费视频| 欧美zozozo另类| 成人av一区二区三区在线看| 特大巨黑吊av在线直播| 国产国拍精品亚洲av在线观看 | 一个人免费在线观看的高清视频| 免费无遮挡裸体视频| 久99久视频精品免费| avwww免费| 日日摸夜夜添夜夜添小说| 女人十人毛片免费观看3o分钟| 悠悠久久av| 成熟少妇高潮喷水视频| 床上黄色一级片| 少妇熟女aⅴ在线视频| 99在线视频只有这里精品首页| av国产免费在线观看| 免费观看精品视频网站| 99国产精品一区二区蜜桃av| 久久久久久九九精品二区国产| 欧美丝袜亚洲另类 | 亚洲成人精品中文字幕电影| 国产aⅴ精品一区二区三区波| 18禁在线播放成人免费| 午夜福利在线观看免费完整高清在 | 亚洲人成伊人成综合网2020| 久久久久久久久大av| 搡老熟女国产l中国老女人| 男插女下体视频免费在线播放| 午夜福利18| 黑人欧美特级aaaaaa片| 搡老妇女老女人老熟妇| 每晚都被弄得嗷嗷叫到高潮| 一卡2卡三卡四卡精品乱码亚洲| 女生性感内裤真人,穿戴方法视频| 亚洲av二区三区四区| 一级黄片播放器| 久久中文看片网| 欧美乱色亚洲激情| 欧美色视频一区免费| 一区福利在线观看| 18+在线观看网站| 男插女下体视频免费在线播放| 国产真实伦视频高清在线观看 | 国产精品永久免费网站| av在线天堂中文字幕| 国产一级毛片七仙女欲春2| 国产成人欧美在线观看| 少妇的逼好多水| 国语自产精品视频在线第100页| 很黄的视频免费| 国产私拍福利视频在线观看| 内射极品少妇av片p| 午夜福利视频1000在线观看| av国产免费在线观看| 丁香六月欧美| 最后的刺客免费高清国语| 日韩欧美精品v在线| 啦啦啦韩国在线观看视频| 日本 欧美在线| 欧美日本亚洲视频在线播放| 日韩大尺度精品在线看网址| 夜夜爽天天搞| 人人妻人人澡欧美一区二区| 人妻丰满熟妇av一区二区三区| 免费搜索国产男女视频| 精品人妻1区二区| 久久久精品大字幕| 夜夜爽天天搞| 国产三级在线视频| 精品国内亚洲2022精品成人| 中亚洲国语对白在线视频| 亚洲美女黄片视频| 成人一区二区视频在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲成人精品中文字幕电影| 成人欧美大片| 国产精品久久视频播放| 在线视频色国产色| a级一级毛片免费在线观看| 一本久久中文字幕| 午夜福利欧美成人| 久久精品国产综合久久久| tocl精华|