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

    南北極區(qū)冬春季節(jié)大氣平流層CO次峰現(xiàn)象

    2014-09-22 08:02:12史華鋒李倩邵愛梅卞建春呂達(dá)仁
    大氣科學(xué) 2014年1期
    關(guān)鍵詞:極區(qū)廓線平流層

    史華鋒 李倩邵愛梅卞建春呂達(dá)仁

    1蘭州大學(xué)大氣科學(xué)學(xué)院,蘭州730000

    2中國(guó)科學(xué)院大氣物理研究所中層大氣和全球環(huán)境探測(cè)重點(diǎn)實(shí)驗(yàn)室,北京100029

    3中國(guó)人民解放軍93601部隊(duì),大同037006

    1 引言

    CO作為一種重要的大氣污染物,其在大氣中的濃度分布和時(shí)空變化受到了廣泛關(guān)注。在大氣對(duì)流層,通過相應(yīng)的化學(xué)反應(yīng)過程,CO對(duì)大氣中OH和臭氧等重要成分的濃度具有顯著的調(diào)制作用,從而影響了大氣化學(xué)系統(tǒng)進(jìn)而產(chǎn)生了相應(yīng)的氣候效應(yīng)。此外,較長(zhǎng)的化學(xué)壽命和較大的濃度梯度,也使CO常被用作研究大氣運(yùn)動(dòng)和污染物輸送的示蹤物,例如近年來在亞洲夏季風(fēng)區(qū)上對(duì)流層—下平流層衛(wèi)星觀測(cè)到的CO等污染物積聚形成的“污池”現(xiàn)象(Park et al.,2008),進(jìn)一步揭示了亞洲夏季風(fēng)在大氣平流層—對(duì)流層交換過程中的重要作用。

    大氣對(duì)流層CO主要來自于含碳物質(zhì)的不完全燃燒、生物活動(dòng)和CH4氧化,通過和OH的化學(xué)反應(yīng)被清除。大氣平流層CO主要來自于CH4與O、Cl和 OH等成分的氧化反應(yīng)。同樣,平流層中的CO也會(huì)與該區(qū)域較為富集的OH產(chǎn)生化學(xué)反應(yīng)。大氣平流層上部及大氣中間層的CO主要來自于該層次大氣中CO2的光解,因而中高層大氣中的CO濃度較高(Hays and Olivero, 1970)。圖1為Aura-MLS衛(wèi)星觀測(cè)(詳細(xì)介紹見 2.觀測(cè)資料)南北緯82°之間、垂直方向215~0.1 hPa的2005~2010年CO濃度多年緯向平均分布。該圖展示了地球大氣中CO的平均分布狀況。可以看出,大氣中CO在垂直方向存在明顯的濃度梯度分布。其中,熱帶地區(qū)大氣對(duì)流層頂附近CO濃度較高,得益于熱帶地區(qū)旺盛的垂直向上輸送過程。而相應(yīng)的平流層中下部(100~10 hPa)CO濃度較低。自平流層上部開始,CO2的光解則使CO濃度急劇增大,在0.1 hPa高度平均濃度(本文所指濃度均為體積混合比濃度)高達(dá)0.5~1.5 ppm(10?6,下同)。

    在中層大氣存在著極向極的大尺度環(huán)流(Shepherd,2007),其上升支位于夏半球,下沉支位于冬半球。通過對(duì)具有中間層“屬性”的大氣成分(CO、H、NOx)的觀測(cè)證明冬半球的下沉運(yùn)動(dòng)貫穿中間層—平流層,中間層的空氣可下沉至平流層中部(Fisher et al.,1993;Filipiak et al.,2005;Engel et al.,2006)。這種季節(jié)性的下沉運(yùn)動(dòng)在圖1中表現(xiàn)為在平流層上部自赤道向高緯地區(qū)CO濃度的顯著增加?;诘鼗鵆O的觀測(cè)研究,Solomon et al.(1985)很早就提出將CO作為高層大氣運(yùn)動(dòng)示蹤物,尤其是高緯地區(qū)夏季的上升運(yùn)動(dòng)和冬季的下沉運(yùn)動(dòng)。此后CO被廣泛用于研究中間層及以上大氣的大尺度環(huán)流和平流層—中間層的垂直傳輸(Manney et al.,1999;Allen et al.,2000;de Zafra et al.,2004;Lee et al.,2011)。

    圖1 2005~2010年Aura-MLS CO濃度(ppb,10?9,下同)多年緯向平均分布Fig.1 Zonal mean of Aura-MLS (Aura satellite,Microwave Limb Sounder) CO volume mixing ratio (VMR) during 2005–2010

    衛(wèi)星觀測(cè)是研究中層大氣化學(xué)系統(tǒng)和動(dòng)力結(jié)構(gòu)的重要手段(呂達(dá)仁和陳洪濱,2003;陳洪濱,2009)。搭載于加拿大SCISAT-1衛(wèi)星的傅里葉變換光譜儀(ACE-FTS)和美國(guó) Aura衛(wèi)星的微波臨邊探測(cè)器(MLS)是目前仍在運(yùn)行的兩種衛(wèi)星觀測(cè)儀器,并已積累了較長(zhǎng)時(shí)間的觀測(cè)資料(詳細(xì)介紹見2.觀測(cè)資料)。我們?cè)趯?duì)2004~2008年的ACE-FTS的CO廓線資料進(jìn)行分析時(shí)發(fā)現(xiàn),冬春季節(jié)南北半球高緯地區(qū)有相當(dāng)一部分廓線在平流層中部存在次峰值現(xiàn)象,采用 Manney et al.(2007)提出的ACE-FTS與Aura-MLS的采樣比較方法,使用Aura-MLS資料對(duì)這一現(xiàn)象進(jìn)行了確認(rèn)。此前Clerbaux et al.(2005)曾分析了 ACE-FTS衛(wèi)星觀測(cè)數(shù)據(jù)前八個(gè)月的資料,文中的廓線圖在平流層也展現(xiàn)了 CO次峰值現(xiàn)象,但該研究并未就此現(xiàn)象給予進(jìn)一步的詳細(xì)分析。本文將利用ACE-FTS和Aura-MLS衛(wèi)星觀測(cè)數(shù)據(jù),結(jié)合MERRA氣象再分析資料,深入研究極區(qū)大氣平流層CO次峰現(xiàn)象的時(shí)空特征及其形成機(jī)制。

    2 觀測(cè)資料

    2.1 ACE-FTS

    大氣化學(xué)實(shí)驗(yàn)(Atmospheric Chemistry Experiment,簡(jiǎn)稱 ACE)是由加拿大主持的一項(xiàng)地球大氣衛(wèi)星探測(cè)計(jì)劃(Bernath et al.,2005)。隸屬于ACE的衛(wèi)星SCISAT-1于2003年8月12日發(fā)射,圓形極地軌道,高度650 km,傾角74°。SCISAT-1搭載的主要儀器之一是一個(gè)高光譜分辨率(0.02 cm?1)的傅里葉變換光譜儀(Fourier Transform Spectrometer,簡(jiǎn)稱 FTS),工作波段 750~4400 cm?1,通過太陽掩星方式探測(cè)南北緯 85°之間、云頂至 110 km高空范圍內(nèi)的大氣溫、濕和幾十種大氣成分垂直分布。ACE-FTS 的 CO資料垂直分辨率為 3~4 km。初期的數(shù)據(jù)驗(yàn)證工作和隨后的伴隨驗(yàn)證均表明,F(xiàn)TS對(duì)CO的觀測(cè)誤差很小,上對(duì)流層至40 km優(yōu)于5%,40 km以上優(yōu)于10%(Clerbaux et al.,2008)。本文將誤差小于10%作為遴選數(shù)據(jù)的標(biāo)準(zhǔn)。

    本文采用的ACE-FTS CO資料包括v2.2(2004年1月至2007年8月)版本和v3.0(2008年1月至2008年6月)版本,共12402條大氣CO廓線。圖2為ACE-FTS觀測(cè)資料2004年1月至2008年6月所有廓線的全球分布圖。該圖顯示,ACE-FTS在中高緯度的觀測(cè)密度比在低緯度大的多,這是由于 ACE主要科學(xué)目標(biāo)之一是研究高緯地區(qū)的大氣化學(xué)、動(dòng)力過程。

    圖2 2004年1月至2008年6月ACE-FTS 所有廓線觀測(cè)點(diǎn)的全球分布Fig.2 Global scatter distribution of the geographical locations of ACE-FTS satellite observation from January 2004 to June 2008

    2.2 Aura-MLS

    微波臨邊探測(cè)器(Microwave Limb Sounder,簡(jiǎn)稱 MLS)是美國(guó)地球觀測(cè)系統(tǒng)(Earth Observing System,EOS)Aura衛(wèi)星載荷儀器之一(Waters et al.,2006)。Aura-MLS通過探測(cè)衛(wèi)星軌道前方地球大氣切線方向的熱輻射對(duì)大氣成分、溫度和冰云等進(jìn)行觀測(cè)。該衛(wèi)星于2004年7月15日發(fā)射,近極地太陽同步軌道,高度705 km。Aura-MLS每天可獲得多達(dá)3500條的CO廓線資料,范圍覆蓋南北緯82°之間,垂直可用范圍215~0.0046 hPa,垂直分辨率3.5~5.0 km。本文所用的 Aura-MLS CO衛(wèi)星資料為第二級(jí)(Level 2)v3.3版。表1所列為MLS CO在不同高度上的水平分辨率、觀測(cè)精度和系統(tǒng)偏差(Pumphrey et al.,2007; Livesey et al.,2008)。由于觀測(cè)存在較大的隨機(jī)誤差,因此在使用 Aura-MLS資料時(shí)要盡可能的采用多線平均以去除隨機(jī)誤差的影響。

    表1 Aura-MLS CO資料參數(shù)(引自http://disc.gsfc.nasa.gov/)Table 1 Data information summary for Aura–MLS CO(quoted from http://disc.gsfc.nasa.gov/)

    本文所用的OH資料也來自于Aura-MLS v3.3版,資料的可用范圍為32~0.0032 hPa,關(guān)于資料的分辨率、觀測(cè)精度和系統(tǒng)偏差及數(shù)據(jù)質(zhì)量控制詳情可參閱Aura-MLS v3.3資料手冊(cè)(http://disc.gsfc.nasa.gov/)和相關(guān)文獻(xiàn)(Pickett et al.,2006)。

    2.3 MERRA

    MERRA(the Modern-Era Retrospective analysis for Research and Applications)是 GMAO(the Global Modeling and Assimilation Office)利用三維變分資料同化系統(tǒng) GEOS-5 DAS生成的氣象再分析資料(Rienecker et al.,2011),其特點(diǎn)是吸收了大量的衛(wèi)星觀測(cè)資料和輻射觀測(cè)資料,并使用了一系列最新的物理過程方案和資料同化方法。MERRA資料采用同化模式水平分辨率 0.67°×0.5°,垂直方向72層,模式高度 0.01 hPa。其產(chǎn)品分辨率因要素而異,本文使用的風(fēng)場(chǎng)資料即來自 MERRA數(shù)據(jù)集,水平網(wǎng)格距1.25°×1.25°,垂直方向42層,垂直高度至0.1 hPa。

    3 研究方法和結(jié)果分析

    3.1 現(xiàn)象展示

    如前文所述,Aura-MLS的觀測(cè)頻次和空間覆蓋范圍均大于ACE-FTS,在對(duì)比兩種衛(wèi)星資料時(shí),往往需要考慮如何從 Aura-MLS數(shù)據(jù)中選擇與ACE-FTS數(shù)據(jù)在觀測(cè)時(shí)間和地點(diǎn)上較為一致的觀測(cè)廓線。Manney et al.(2007)提出了一種較為可行的辦法,即在ACE-FTS觀測(cè)時(shí)間前后12小時(shí),觀測(cè)點(diǎn)南北各1個(gè)緯度之間、東西各8個(gè)經(jīng)度范圍內(nèi)的 Aura-MLS平均廓線被近似認(rèn)為與該 ACE-FTS廓線具有“同性質(zhì)”。

    本文采用Manney et al.(2007)提出的方法,對(duì)Aura-MLS和ACE-FTS衛(wèi)星觀測(cè)資料進(jìn)行分析,并挑出兩類衛(wèi)星資料中具有“同性質(zhì)”的廓線進(jìn)行比較。圖3顯示了南北緯70°~80°之間的CO廓線在不同季節(jié)的分布特征。從該圖中可以看到,ACE-FTS和Aura-MLS衛(wèi)星觀測(cè)數(shù)據(jù)均顯示在南北半球冬春季節(jié)的平流層中部(約10 hPa)存在CO濃度次峰值現(xiàn)象,其中次峰值出現(xiàn)的位置和峰值濃度均有很好的一致性,而在夏秋季節(jié)則不具備這種特征。圖 3a、b顯示北半球冬季(DJF時(shí)段)CO濃度的峰值高度出現(xiàn)在約8 hPa附近,峰值濃度約300 ppb;春季(MAM時(shí)段)時(shí)峰值高度略有下降,出現(xiàn)在約10 hPa高度,峰值濃度降為約180 ppb。圖3c、d顯示南半球冬季(JJA時(shí)段)峰值高度出現(xiàn)在約8 hPa,峰值濃度約420 ppb;而春季(SON時(shí)段)時(shí)峰值出現(xiàn)的高度沒有出現(xiàn)明顯下降的現(xiàn)象,僅峰值濃度下降為300 ppb左右。綜合比較來看,相較北半球而言,南半球冬春季節(jié)平流層中部CO濃度廓線具有更顯著的次峰值結(jié)構(gòu)特征。

    本文研究發(fā)現(xiàn),冬春季節(jié)平流層中部CO濃度次峰現(xiàn)象不僅出現(xiàn)在高緯地區(qū),甚至可擴(kuò)展到南北緯 50°附近,但是相對(duì)高緯地區(qū),次峰現(xiàn)象在中緯地區(qū)出現(xiàn)的幾率較低。而這一現(xiàn)象分別開始于2~3月份(北半球)和8~9月份(南半球),即冬半球的晚冬初春季節(jié)。

    利用衛(wèi)星觀測(cè)資料對(duì)南北半球極區(qū)大氣平流層CO濃度分布進(jìn)行2005~2010多年單日平均(圖

    4),結(jié)果表明當(dāng)次峰值現(xiàn)象出現(xiàn)時(shí),CO濃度分布具有在垂直方向上呈懸球狀的空間分布結(jié)構(gòu),即在平流層的中部出現(xiàn)一個(gè)較大的CO富集區(qū),隨高度增加CO的濃度減少,在平流層頂附近時(shí)濃度又隨高度增加。

    圖3 Aura-MLS和ACE-FTS南北緯70°~80°不同季節(jié)CO濃度(ppb)垂直廓線族。MAM:3~5月;JJA:6~8月;SON:9~11月;DJF:12~2月Fig.3 Comparisons of the CO concentration vertical profiles (ppb) located in latitudes 70–80 between Aura-MLS and ACE-FTS.MAM (Mar–Apr–May), JJA(Jun–Jul–Aug), SON (Sep–Oct–Nov), and DJF (Dec–Jan–Feb), respectively

    圖4 Aura-MLS CO 濃度(ppb)2005~2010年多年單日平均在不同高度上的水平分布。(a)2月14日北半球;(b)8月28日南半球。左起第一列:1.0 hPa高度;第二列:2.1 hPa高度;第三列:4.6 hPa高度;第四列:10 hPa高度Fig.4 Distribution of daily mean CO mixing ratio (ppb) from Aura-MLS data (2005–2010) at different pressure levels for (a) the Northern Hemisphere on February 14 and (b) the Southern Hemisphere on August 28

    3.2 機(jī)制分析

    本節(jié)首先對(duì) 2005~2010年南北半球高緯地區(qū)大氣平流層CO濃度變化進(jìn)行了分析,并具體以較為典型的2007年南半球冬春季節(jié)極區(qū)平流層CO的濃度變化為例,探討CO次峰現(xiàn)象的形成機(jī)制。

    圖5展示了2005~2010年南北半球平流層中上部不同高度高緯至極地地區(qū)(緯度 60°~82°之間)Aura-MLS CO平均濃度隨時(shí)間的變化。在北半球(圖5a)1 hPa高度上 CO的濃度峰值出現(xiàn)于11月至來年1月;而在南半球(圖5b)1 hPa高度上 CO的濃度峰值則出現(xiàn)于4月至9月,持續(xù)時(shí)間幾乎是北半球的兩倍,表明了南北半球平流層大氣化學(xué)系統(tǒng)時(shí)空特征上的差異。相較于1 hPa高度,較低高度上的極區(qū)大氣中也同樣依次對(duì)應(yīng)出現(xiàn)了CO的濃度峰值,但隨高度下降該峰值發(fā)生時(shí)間相應(yīng)延遲,同時(shí)峰值濃度減小,持續(xù)時(shí)間縮短。通過圖5可以看出,南半球和北半球的大氣平流層中上層CO濃度峰值都存在著顯著的年際變化,其中北半球的峰值結(jié)構(gòu)年際變化過程更為復(fù)雜,例如 2005~2006年北半球極區(qū)冬季平流層頂CO濃度的表現(xiàn)為雙峰結(jié)構(gòu),值得關(guān)注的是在此期間出現(xiàn)了近年來北半球最強(qiáng)的平流層爆發(fā)性增溫過程(Manney et al.,2008)。平流層大氣CO濃度是極低的,水平方向的輸送不會(huì)引起極區(qū)CO濃度的顯著上升,只有通過平流層以上大氣的下沉過程才有機(jī)會(huì)造成平流層大氣中 CO濃度的上升,因此大氣中自上而下的CO峰值濃度的遞減和時(shí)間位相上的滯后是冬季中間層大氣成分向平流層大氣垂直輸送過程的反映。

    除了垂直輸送過程(下文以S1表示),高緯地區(qū)大氣平流層中CO濃度變化還取決于其他幾個(gè)因素:CO2的光解和CH4氧化,OH對(duì)CO的氧化清除,水平輸送過程(以S2表示),中括號(hào)表示該種大氣成分的濃度,即存在以下關(guān)系:

    圖5顯示在中間層—平流層垂直運(yùn)動(dòng)極弱的夏半年,平流層中的CO濃度極低且很穩(wěn)定,這一較低的濃度水平主要是由CO2光解、CH4氧化和OH清除之間的化學(xué)平衡來維持的。冬春季節(jié)光照減弱,溫度降低,極渦內(nèi)的大氣成分化學(xué)活性下降,來自中間層的下沉空氣使極區(qū)CH4濃度降低(劉傳熙等,2009),ACE-FTS衛(wèi)星資料分析結(jié)果顯示這種下沉運(yùn)動(dòng)對(duì)CO2影響不大(具體表現(xiàn)為CO2在中層大氣垂直濃度分布基本一致,圖略),因此可忽略CO2光解、CH4氧化對(duì)CO濃度變化的影響。OH的濃度很低,但是化學(xué)活性很強(qiáng),其對(duì)CO的化學(xué)清除作用不可忽略。此外冬末春初時(shí)平流層極渦變?nèi)醪②呌诒罎?,為極渦內(nèi)外空氣的交換創(chuàng)造了條件(Manney et al., 2009)。接下來本文將進(jìn)一步分析S1、S2和OH濃度變化對(duì)極區(qū)平流層CO濃度的影響。此外,在大氣平流層中O3對(duì)CO的作用比較微弱(Solomon et al., 1985;Minschwaner et al., 2010),因此本文不針對(duì) O3與 CO的相互作用展開深入討論。

    圖6a顯示8月份平流層上部CO濃度基本保持不變,從9月初開始,平流層上部CO濃度迅速下降,并首次低于平流層中部的CO濃度,相應(yīng)的廓線分析顯示(圖未給出)CO垂直分布出現(xiàn)次峰值現(xiàn)象,此后僅一個(gè)月時(shí)間,平流層上部(1 hPa高度左右)CO濃度便迅速降低并恢復(fù)至夏半年水平,而平流層中部(10 hPa高度左右)CO濃度峰值則一直持續(xù)至11月底。圖6b顯示8~11月份高緯地區(qū)OH濃度整體處于穩(wěn)步回升階段。OH主要來自于光解反應(yīng),與太陽光照密切相關(guān),OH濃度的回升與太陽直射點(diǎn)的南移相對(duì)應(yīng)。但是平流層中部的OH濃度遠(yuǎn)低于上部(請(qǐng)注意圖中為指數(shù)色標(biāo)),即使在 11月末,平流層中部的 OH濃度也只有約0.01 ppb。圖6c顯示8~11月份環(huán)南緯60°的西風(fēng)大值中心隨著時(shí)間推移向下移動(dòng),與CO濃度峰值的變化特征相對(duì)應(yīng)。8~11月份平流層中下部始終維持西風(fēng)環(huán)流,但是在九月中旬,平流層上部的西風(fēng)環(huán)流逐漸轉(zhuǎn)變?yōu)闁|風(fēng)環(huán)流,結(jié)合溫度場(chǎng)的變化特征(圖未給出),該過程對(duì)應(yīng)了南半球高緯平流層的最后一次增溫過程。圖6d為南緯60°經(jīng)向風(fēng)V北風(fēng)分量的緯向平均隨時(shí)間的變化,代表了中緯地區(qū)大氣向極區(qū)大氣輸送的特征。與緯向風(fēng)U的變化特征相比,經(jīng)向風(fēng)V的北風(fēng)分量變化較為復(fù)雜,八月上旬和中旬平流層上部北風(fēng)分量的分別出現(xiàn)兩次異常加強(qiáng),自9月初開始,北風(fēng)分量持續(xù)增強(qiáng),并在緯向風(fēng)由西風(fēng)環(huán)流轉(zhuǎn)變?yōu)闁|風(fēng)環(huán)流的同時(shí)V風(fēng)速達(dá)到極大,這一變化與CO的濃度的變化密切對(duì)應(yīng)。圖6e為南緯60°經(jīng)向風(fēng)北風(fēng)分量的平均渦動(dòng)V′隨時(shí)間的變化,代表了中緯地區(qū)向極區(qū)大氣渦動(dòng)輸送的變化(負(fù)值代表指向極區(qū)的渦動(dòng)輸送,正值相反)。V′與經(jīng)向風(fēng)V的北風(fēng)分量的變化特征相似,但在平流層上部的變化幅度最大,而在平流層中下部變化幅度較小。綜合對(duì)比以上分析可以發(fā)現(xiàn),平流層上部CO濃度的急劇變化與緯向風(fēng)的減弱及逆轉(zhuǎn)、經(jīng)向輸送的增強(qiáng)和OH濃度的指數(shù)式增加相對(duì)應(yīng),而這些因素都能夠影響到大氣中CO的濃度并在冬春季節(jié)導(dǎo)致平流層上層CO濃度降低。與此同時(shí),平流層中下部CO濃度的緩慢變化則與始終維持的西風(fēng)環(huán)流、較弱的經(jīng)向輸送和極低的OH濃度相對(duì)應(yīng)。

    圖5 南北緯60°~82°之間Aura-MLS CO平均濃度(ppb)隨時(shí)間的變化:(a)北半球;(b)南半球Fig.5 Variation of the Aura-MLS CO mixing ratio (ppb) averaged over latitudes 60°~82°: (a) Northern Hemisphere; (b) Southern Hemisphere

    此外值得關(guān)注的是以CO的800 ppb(約0.5 hPa高度極區(qū)大氣的平均CO濃度)等值線的變化來看,圖 6a顯示至少在八月初 800 ppb等值線便開始收縮,意味著此時(shí)中間層大氣中的CO向下輸送被減弱,因而削弱乃至切斷了極區(qū)平流層CO的主要來源。

    4 結(jié)論和討論

    在 ACE-FTS 和 Aura-MLS 衛(wèi)星觀測(cè)數(shù)據(jù) CO廓線資料的分析中,我們發(fā)現(xiàn)冬春季節(jié)高緯地區(qū)大氣平流層中部存在著CO濃度的次峰值現(xiàn)象,并對(duì)

    這一現(xiàn)象出現(xiàn)的范圍和時(shí)間進(jìn)行了詳細(xì)的分析。分析結(jié)果表明平流層CO次峰值現(xiàn)象具有顯著的空間分布特性和季節(jié)變化特征,該現(xiàn)象主要出現(xiàn)在高緯地區(qū),偶爾可擴(kuò)展到南北緯 50°附近,但是相對(duì)高緯地區(qū),次峰現(xiàn)象在中緯地區(qū)出現(xiàn)的幾率較低。通常這一現(xiàn)象出現(xiàn)于晚冬初春季節(jié),其中北半球 CO次峰現(xiàn)象相較南北球 CO次峰現(xiàn)象持續(xù)時(shí)間較短、強(qiáng)度較弱。當(dāng)次峰現(xiàn)象出現(xiàn)時(shí),相應(yīng)時(shí)期大氣中 CO的緯向平均呈現(xiàn)為懸球狀的空間分布特征,即在極區(qū)大氣平流層中部CO的濃度較高,而在平流層上部CO的濃度較低。

    圖6 2007年8~11月南半球高緯地區(qū)各要素隨時(shí)間的變化:(a)南緯60°~82°之間Aura-MLS CO的平均濃度(ppb);(b)南緯60°~82°之間Aura-MLS OH的平均濃度(ppb);(c)南緯60°緯向風(fēng)U(m s?1)的緯向平均;(d)南緯60°經(jīng)向風(fēng)V(m s?1)北風(fēng)分量的緯向平均;(e)南緯60°經(jīng)向風(fēng)V(m s?1)北風(fēng)分量的平均渦動(dòng)V'Fig.6 (a) The evolution of the mean mixing ratio of Aura-MLS CO (ppb) averaged in 60°S–82°S; (b) same as in (a), but for the Aura-MLS OH mixing ratio(ppb); (c) the MERRA U-wind (m s?1) averaged at 60°S; (d) same as in (c), but for the southward MERRA V-wind (m s?1); (e) same as in (c), but for the transient southward MERRA V-wind (m s?1).All figures are plotted from Aug to Nov, 2007

    基于Aura-MLS 衛(wèi)星觀測(cè)CO資料,在綜合分析了2005~2010年間平流層—中間層CO濃度變化的基礎(chǔ)上,本文以 2007年南半球冬春季節(jié)的 CO濃度變化為例,重點(diǎn)研究了在CO次峰結(jié)構(gòu)形成過程中的相關(guān)因素,并深入探討了各類相關(guān)動(dòng)力傳輸過程和化學(xué)反應(yīng)過程的作用機(jī)制。結(jié)果表明,平流層上部 CO濃度的急劇減少與緯向風(fēng)的減弱及逆轉(zhuǎn)、經(jīng)向輸送的增強(qiáng)和大氣中OH濃度的指數(shù)式增加相對(duì)應(yīng),而平流層中部CO濃度的緩慢變化則與始終維持的西風(fēng)環(huán)流、較弱的經(jīng)向輸送和極低的OH濃度相對(duì)應(yīng)。

    在進(jìn)行CO次峰現(xiàn)象的成因機(jī)制分析中,本文也同時(shí)考慮到CO2光解、CH4氧化等化學(xué)反應(yīng)過程的影響。冬春季節(jié)光照的減弱和氣溫的下降,使極區(qū)的整個(gè)化學(xué)反應(yīng)體系活性下降,因而可忽略CO2光解、CH4氧化對(duì)CO濃度變化的影響。值得關(guān)注的是與來自中間層CO一并下傳的還有大量的氮氧化物 NOx(Engel et al.,2006; Siskind et al.,2007),而 NOx在 CH4氧化為 CO的過程中則是不可或缺的。此外,我們發(fā)現(xiàn)中緯大氣平流層偶爾也會(huì)出現(xiàn)CO次峰現(xiàn)象,而這一現(xiàn)象很可能與極渦的分裂、偏心有關(guān),同時(shí)一些相應(yīng)的動(dòng)力過程諸如行星波的破碎、平流層爆發(fā)性增溫等對(duì)CO分布變化存在著重要的影響。此外我們也發(fā)現(xiàn)南半球和北半球大氣中間層向大氣平流層傳輸過程持續(xù)時(shí)間及強(qiáng)度存在著顯著差異,以及南、北半球大氣平流層CO空間分布結(jié)構(gòu)也存在著不同,這些都需要進(jìn)行更深入的分析。結(jié)合更為細(xì)致和更長(zhǎng)期的衛(wèi)星觀測(cè)資料分析結(jié)果,同時(shí)結(jié)合目前得到廣泛應(yīng)用的中層大氣化學(xué)傳輸模式(田文壽等,2009),將非常有利于未來對(duì)大氣中CO濃度時(shí)空分布和變化特征與機(jī)制,展開進(jìn)一步的研究工作。

    致謝 本文所用衛(wèi)星資料分別來自于NASA的Aura-MLS(http://disc.sci.gsfc.nasa.gov/Aura)和加拿大 CSA 的 ACE-FTS(http://www.ace.uwaterloo.ca/index.html),MERRA同化資料由 NASA的 GMAO(https://gmao.gsfc.nasa.gov/products/)提供,在此對(duì)上述機(jī)構(gòu)一并表示感謝。兩位匿名審稿人提出的寶貴意見,使本文增色不少,謹(jǐn)致以衷心的感謝。

    References)

    Allen D R, Stanford J L, Nakamura N, et al.2000.Antarctic polar descent and planetary wave activity observed in ISAMS CO from April to July 1992 [J].Geophys.Res.Lett., 27: 665–668.

    Bernath P F, McElroy C T, Abrams M C, et al.2005.Atmospheric Chemistry Experiment (ACE): Mission overview [J].Geophys.Res.Lett.,32 (15): L15S01, doi:10.1029/2005GL022386.

    陳洪濱.2009.中高層大氣研究的空間探測(cè) [J].地球科學(xué)進(jìn)展, 24 (3):229–241.Chen Hongbin.2009.An overview of the space-based observations for upper atmospheric research [J].Advances in Earth Science (in Chinese), 24 (3): 229–241.

    Clerbaux C, Coheur P F, Hurtmans D, et al.2005.Carbon monoxide distribution from the ACE-FTS solar occultation measurements [J].Geophys.Res.Lett., 32 (16): L16S01, doi:10.1029/2005GL022394.

    Clerbaux C, George M, Turquety S, et al.2008.CO measurements from the ACE-FTS satellite instrument: Data analysis and validation using ground-based, airborne and spaceborne observations [J].Atmos.Chem.Phys., 8 (9): 2569–2594.

    de Zafra R L, Muscari G.2004.CO as an important high altitude tracer of dynamics in the polar stratosphere and mesosphere [J].J.Geophys.Res.,109: D06105, doi:10.1029/2003JD004099.

    Engel A, M?bius T, Haase H P, et al.2006.Observation of mesospheric air inside the arctic stratospheric polar vortex in early 2003 [J].Atmos.Chem.Phys., 6: 267–282.

    Fisher M, O’Neill A, Sutton R, et al.1993.Rapid descent of mesospheric air into the stratospheric polar vortex [J].Geophys.Res.Lett., 20: 1267–1270.

    Filipiak M J, Harwood R S, Jiang J H, et al.2005.Carbon monoxide measured by the EOS Microwave Limb Sounder on Aura: First results[J].Geophys.Res.Lett., 32: L14825.

    Hays P B, Olivero J J.1970.Carbon dioxide and monoxide above the troposphere [J].Planet.Space Sci., 18: 1729–1733.

    呂達(dá)仁, 陳洪濱.2003.平流層和中層大氣研究的進(jìn)展 [J].大氣科學(xué),27 (4): 750–769.Lü Daren, Chen Hongbin.2003.Advances in middle atmosphere physics research [J].Chinese Journal of Atmospheric Sciences (in Chinese), 27 (4): 750–769.

    劉傳熙, 王海平, 劉毅.2009.2003~2004年冬季北半球爆發(fā)性增溫期間極地平流層甲烷、水汽的演變特征 [J].自然科學(xué)進(jìn)展, 19 (3): 425–433.Liu Chuanxi, Wang Haiping, Liu Yi.2009.Evolution of the methane and water vapor during SSW in 2003–2004 Northern Hemisphere winter [J].Progress in Natural Science (in Chinese), 19 (3):425–433.

    Livesey N J, Filipiak M J, Froidevaux L, et al.2008.Validation of AuraMicrowave Limb Sounder O3and CO observations in the upper troposphere and lower stratosphere [J].J.Geophys.Res., 113: D15S02,doi: 10.1029/2007JD008805.

    Lee J N, Wu D L, Manney G L, et al.2011.Aura Microwave Limb Sounder observations of the polar middle atmosphere: Dynamics and transport of CO and H2O [J].J.Geophys.Res., 116: D05110, doi:10.1029/2010JD014608.

    Manney G L, Michelson H A, Santee M L, et al.1999.Polar vortex dynamics during spring and fall diagnosed using trace gas observations from the Atmospheric Trace Molecule Spectroscopy instrument [J].J.Geophys.Res., 104: 18841–18866.

    Manney G L, Dafferet W H, Zawodny J M, et al.2007.Solar occultation satellite data and derived meteorological products: Sampling issues and comparisons with Aura Microwave Limb Sounder [J].J.Geophys.Res.,112: D24S50, doi:10.1029/2007JD008709.

    Manney G L, Daffer W H, Strawbridge K B, et al.2008.The high Arctic in extreme winters: Vortex, temperature, and MLS and ACE-FTS trace gas evolution [J].Atmos.Chem.Phys., 8: 505–522.

    Manney G L, Harwood R S, MacKenzie I A, et al.2009.Satellite observations and modeling of transport in the upper troposphere through the lower mesosphere during the 2006 major stratospheric sudden warming [J].Atmos.Chem.Phys., 9: 4775–4795.

    Minschwaner K, Manney G L, Livesey N J, et al.2010.The photochemistry of carbon monoxide in the stratosphere and mesosphere evaluated from observations by the Microwave Limb Sounder on the Aura satellite [J].J.Geophys.Res., 115: D13303, doi:10.1029/2009JD012654.

    Pickett H M, Drouin B J, Canty T, et al.2006.Validation of Aura MLS HOxmeasurements with remote-sensing balloon instruments [J].Geophys.Res.Lett., 33 (1): L01808, doi:10.1029/2005GL024048.

    Pumphrey H C, Filipiak M J, Livesey N J, et al.2007.Validation of middle atmosphere carbon monoxide retrievals from the Microwave Limb Sounder on Aura [J].J.Geophys.Res., 112: D24S38, doi:10.1029/2007JD008723.

    Park M, Randel W J, Emmons L K, et al.2008.Chemical isolation in the Asian monsoon anticyclone observed in Atmospheric Chemistry Experiment (ACE-FTS) data [J].Atmos.Chem.Phys., 8: 757–764.

    Rienecker M, Suarez M J, Gelaro R, et al.2011.MERRA: NASA’s Modern-Era Retrospective Analysis for Research and Applications [J].J.Climate, 24: 3624–3648.

    Solomon S, Garcia R R, Olivero J J, et al.1985.Photochemistry and transport of carbon monoxide in the middle atmosphere [J].J.Atmos.Sci.,42 (10): 1072–1083.

    Shepherd T G.2007.Transport in the middle atmosphere [J].J.Meteor.Soc.Japan., 85: 165–191.

    Siskind D E, Eckermann S D, Coy L, et al.2007.On recent interannual variability of the Arctic winter mesosphere: Implications for tracer descent [J].Geophys.Res.Lett., 34: L09806, doi:10.1029/2007GL029293.

    田文壽, 張敏, 舒建川.2009.中層大氣模式的應(yīng)用及發(fā)展前景 [J].地球科學(xué)進(jìn)展, 24 (3): 252–261.Tian Wenshou, Zhang Min, Shu Jianchuan.2009.The applications and future development of middle atmosphere models [J].Advances in Earth Science (in Chinese), 24 (3):252–261.

    Waters J W, Froidevaux L, Harwood R S, et al.2006.The Earth Observing System Microwave Limb Sounder (EOS MLS) on the Aura satellite [J].IEEE Trans.Geosci.Remote, 44 (5): 1075–1092.

    猜你喜歡
    極區(qū)廓線平流層
    青藏高原上空平流層水汽的時(shí)空演變特征
    基于副氣囊的平流層浮空器高度控制
    直升機(jī)極區(qū)海洋環(huán)境加速試驗(yàn)環(huán)境譜研究
    不同降水強(qiáng)度下風(fēng)廓線雷達(dá)譜矩特征與測(cè)風(fēng)準(zhǔn)確性分析
    一種極區(qū)統(tǒng)一坐標(biāo)系及其導(dǎo)航參數(shù)轉(zhuǎn)換方法*
    同化Argo海洋廓線觀測(cè)對(duì)CAS-ESM-C的上層海洋溫度模擬的改進(jìn)
    1979~2011年間平流層溫度及平流層水汽的演變趨勢(shì)
    基于快速局域線性回歸的IRAS/FY-3B大氣溫濕廓線反演
    GRAPES模式頂外部背景廓線構(gòu)造方法初步研究
    極區(qū)間接橫向慣性導(dǎo)航方法
    亚洲在线自拍视频| 中文亚洲av片在线观看爽| 欧美+亚洲+日韩+国产| 2021天堂中文幕一二区在线观| 日本黄色视频三级网站网址| 日韩欧美精品v在线| 亚洲自拍偷在线| 色综合站精品国产| 国产亚洲欧美在线一区二区| 免费人成视频x8x8入口观看| 18美女黄网站色大片免费观看| 国产精品亚洲av一区麻豆| 亚洲欧美精品综合久久99| 人人妻人人看人人澡| 不卡一级毛片| www.999成人在线观看| 精品久久久久久久末码| 精品久久久久久久久久久久久| xxx96com| 日本黄大片高清| 一进一出好大好爽视频| 久久久久久久精品吃奶| 国产精品亚洲av一区麻豆| 色综合婷婷激情| 少妇粗大呻吟视频| 亚洲,欧美精品.| 久久精品91无色码中文字幕| 美女扒开内裤让男人捅视频| 国产乱人伦免费视频| 在线视频色国产色| 亚洲男人天堂网一区| 国模一区二区三区四区视频 | 久久精品国产99精品国产亚洲性色| 人妻丰满熟妇av一区二区三区| 日本成人三级电影网站| 国产精品,欧美在线| 母亲3免费完整高清在线观看| 国产成人欧美在线观看| 国产伦在线观看视频一区| 午夜福利在线在线| 久久精品综合一区二区三区| 在线观看免费午夜福利视频| 成人精品一区二区免费| 午夜视频精品福利| 亚洲无线在线观看| 中文字幕精品亚洲无线码一区| 免费在线观看日本一区| 神马国产精品三级电影在线观看 | 亚洲av成人一区二区三| 日韩精品青青久久久久久| 九色成人免费人妻av| 欧美成狂野欧美在线观看| 国产又黄又爽又无遮挡在线| 亚洲一区二区三区色噜噜| 99热只有精品国产| 久久香蕉激情| 久久久久久久久中文| 999久久久国产精品视频| 亚洲18禁久久av| 一边摸一边抽搐一进一小说| 国产精品永久免费网站| 国产精品免费视频内射| 老汉色∧v一级毛片| 99热只有精品国产| 操出白浆在线播放| 国产亚洲精品综合一区在线观看 | 国内精品久久久久久久电影| 老司机深夜福利视频在线观看| 女人高潮潮喷娇喘18禁视频| 日本a在线网址| 少妇的丰满在线观看| 欧美绝顶高潮抽搐喷水| 国产精品香港三级国产av潘金莲| 国产成人影院久久av| 久久久久精品国产欧美久久久| av免费在线观看网站| 男女做爰动态图高潮gif福利片| 亚洲aⅴ乱码一区二区在线播放 | 国产成年人精品一区二区| 禁无遮挡网站| 免费一级毛片在线播放高清视频| 国产三级在线视频| av在线天堂中文字幕| e午夜精品久久久久久久| 搡老熟女国产l中国老女人| 精品人妻1区二区| 亚洲人成77777在线视频| 无人区码免费观看不卡| 看片在线看免费视频| 国产精品久久久久久亚洲av鲁大| 51午夜福利影视在线观看| 老司机午夜福利在线观看视频| 特级一级黄色大片| 搡老岳熟女国产| 国产成人精品无人区| 精品久久久久久成人av| 一区二区三区激情视频| 亚洲av美国av| 欧美色视频一区免费| 欧美一级a爱片免费观看看 | 欧美极品一区二区三区四区| 国产精品国产高清国产av| 久久亚洲精品不卡| 十八禁人妻一区二区| 99热这里只有是精品50| 日韩欧美 国产精品| 亚洲国产欧洲综合997久久,| 久久久久性生活片| 欧美 亚洲 国产 日韩一| 一级a爱片免费观看的视频| av视频在线观看入口| 亚洲avbb在线观看| av福利片在线观看| 九色成人免费人妻av| 国产成人aa在线观看| 精品久久久久久久人妻蜜臀av| 亚洲av成人不卡在线观看播放网| 免费看美女性在线毛片视频| 国产单亲对白刺激| 国产精品亚洲一级av第二区| 亚洲精品粉嫩美女一区| 成人国语在线视频| 在线观看免费视频日本深夜| 亚洲人成77777在线视频| 少妇的丰满在线观看| 成年版毛片免费区| АⅤ资源中文在线天堂| 国产成人欧美在线观看| 精品国产超薄肉色丝袜足j| 亚洲国产欧洲综合997久久,| 国产单亲对白刺激| 欧美日韩黄片免| 国产一区二区在线av高清观看| 老司机午夜十八禁免费视频| 日韩国内少妇激情av| 亚洲片人在线观看| 麻豆成人午夜福利视频| 久久久久久久久中文| 少妇人妻一区二区三区视频| 男人的好看免费观看在线视频 | 欧美日韩乱码在线| 久久久久久久久免费视频了| 村上凉子中文字幕在线| 青草久久国产| 99久久99久久久精品蜜桃| 啦啦啦免费观看视频1| 老鸭窝网址在线观看| 一本精品99久久精品77| 老汉色∧v一级毛片| 欧美性猛交╳xxx乱大交人| 欧美成人免费av一区二区三区| 婷婷精品国产亚洲av| 亚洲av美国av| 久久久久久亚洲精品国产蜜桃av| 午夜两性在线视频| 成在线人永久免费视频| √禁漫天堂资源中文www| 全区人妻精品视频| 丝袜美腿诱惑在线| 亚洲欧洲精品一区二区精品久久久| 不卡一级毛片| 我要搜黄色片| 淫妇啪啪啪对白视频| 午夜免费成人在线视频| 俺也久久电影网| 亚洲av中文字字幕乱码综合| 午夜福利18| cao死你这个sao货| 免费一级毛片在线播放高清视频| 久久精品国产亚洲av香蕉五月| 欧美乱妇无乱码| 国产亚洲精品综合一区在线观看 | 免费观看精品视频网站| 18禁黄网站禁片免费观看直播| 国产精品久久久久久人妻精品电影| 最近最新中文字幕大全电影3| 一区二区三区高清视频在线| 日本 av在线| 丝袜美腿诱惑在线| 真人一进一出gif抽搐免费| 国产av一区在线观看免费| 成年版毛片免费区| 90打野战视频偷拍视频| xxx96com| 在线观看一区二区三区| 久9热在线精品视频| 国产精品亚洲美女久久久| av在线播放免费不卡| 午夜福利视频1000在线观看| 国内精品一区二区在线观看| 床上黄色一级片| 老司机午夜福利在线观看视频| 亚洲熟女毛片儿| 久久久国产精品麻豆| 特大巨黑吊av在线直播| 又粗又爽又猛毛片免费看| 男女之事视频高清在线观看| 久久 成人 亚洲| 天天躁狠狠躁夜夜躁狠狠躁| 久9热在线精品视频| 久久久久免费精品人妻一区二区| 欧美精品啪啪一区二区三区| 亚洲人与动物交配视频| 亚洲av电影不卡..在线观看| 校园春色视频在线观看| 精品高清国产在线一区| 一进一出抽搐gif免费好疼| 日韩欧美三级三区| a在线观看视频网站| 日本一本二区三区精品| 黑人巨大精品欧美一区二区mp4| 黄色女人牲交| a级毛片在线看网站| svipshipincom国产片| 国产乱人伦免费视频| 国产精品av视频在线免费观看| av视频在线观看入口| 校园春色视频在线观看| 老司机在亚洲福利影院| 两个人免费观看高清视频| 这个男人来自地球电影免费观看| 欧美中文综合在线视频| 免费无遮挡裸体视频| 国内揄拍国产精品人妻在线| 久久这里只有精品19| 久久国产精品影院| 一个人免费在线观看电影 | 一级毛片精品| 久久久久久久久久黄片| 成人亚洲精品av一区二区| 少妇粗大呻吟视频| 男人的好看免费观看在线视频 | 日韩av在线大香蕉| 欧美日韩国产亚洲二区| 久久久水蜜桃国产精品网| 亚洲欧美激情综合另类| 国产一级毛片七仙女欲春2| 成人三级黄色视频| 成人av在线播放网站| 亚洲一区二区三区不卡视频| 久久精品综合一区二区三区| 欧美大码av| 日韩欧美国产一区二区入口| 这个男人来自地球电影免费观看| 欧美激情久久久久久爽电影| 成在线人永久免费视频| 日日摸夜夜添夜夜添小说| 国产三级在线视频| 久久久国产成人免费| 精品福利观看| а√天堂www在线а√下载| 一级毛片女人18水好多| 丰满的人妻完整版| 91麻豆av在线| 精品人妻1区二区| 成人手机av| 中文字幕高清在线视频| av天堂在线播放| 免费看十八禁软件| 91在线观看av| 波多野结衣高清作品| 久久久久久亚洲精品国产蜜桃av| 亚洲国产精品合色在线| 亚洲av成人av| 国产aⅴ精品一区二区三区波| 女同久久另类99精品国产91| 亚洲国产欧美人成| 黑人欧美特级aaaaaa片| 国产精品香港三级国产av潘金莲| 少妇的丰满在线观看| 国内少妇人妻偷人精品xxx网站 | 人妻丰满熟妇av一区二区三区| 午夜福利欧美成人| 久久精品国产99精品国产亚洲性色| 国产高清视频在线播放一区| 欧美成狂野欧美在线观看| 国产精品国产高清国产av| 露出奶头的视频| 欧美国产日韩亚洲一区| 在线视频色国产色| 国产区一区二久久| 老熟妇仑乱视频hdxx| 麻豆av在线久日| 亚洲专区字幕在线| 久久久久久大精品| 亚洲av中文字字幕乱码综合| 国产黄a三级三级三级人| 日本三级黄在线观看| 国产男靠女视频免费网站| 欧美色视频一区免费| 亚洲专区国产一区二区| 欧美zozozo另类| 在线国产一区二区在线| 亚洲熟妇熟女久久| 日韩大尺度精品在线看网址| 丰满人妻一区二区三区视频av | 神马国产精品三级电影在线观看 | 中文字幕av在线有码专区| 日韩欧美 国产精品| 国产精品自产拍在线观看55亚洲| 国产精品久久电影中文字幕| 999久久久精品免费观看国产| 欧美日韩中文字幕国产精品一区二区三区| 中文字幕久久专区| 亚洲男人天堂网一区| 午夜成年电影在线免费观看| 亚洲自拍偷在线| 狂野欧美白嫩少妇大欣赏| 亚洲中文字幕日韩| 中国美女看黄片| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩中文字幕国产精品一区二区三区| 啦啦啦韩国在线观看视频| 国产私拍福利视频在线观看| svipshipincom国产片| 97人妻精品一区二区三区麻豆| 狂野欧美白嫩少妇大欣赏| 亚洲中文字幕一区二区三区有码在线看 | 91字幕亚洲| 高清毛片免费观看视频网站| 女同久久另类99精品国产91| 久久久久国内视频| 国产精品一区二区精品视频观看| 男女视频在线观看网站免费 | 国产伦一二天堂av在线观看| 欧美日韩中文字幕国产精品一区二区三区| 日本一本二区三区精品| 88av欧美| 人妻丰满熟妇av一区二区三区| 久久精品国产清高在天天线| 国产亚洲欧美98| 免费无遮挡裸体视频| 亚洲熟女毛片儿| 99精品在免费线老司机午夜| av免费在线观看网站| 亚洲国产看品久久| 日韩欧美免费精品| 女警被强在线播放| 9191精品国产免费久久| 午夜福利18| 最近最新免费中文字幕在线| 国产aⅴ精品一区二区三区波| 一区二区三区高清视频在线| 日韩欧美免费精品| 两个人看的免费小视频| av天堂在线播放| 哪里可以看免费的av片| 国产91精品成人一区二区三区| 老鸭窝网址在线观看| 中文字幕久久专区| 亚洲成人中文字幕在线播放| 麻豆成人午夜福利视频| 欧美日韩精品网址| 欧美黑人欧美精品刺激| 热99re8久久精品国产| 天天一区二区日本电影三级| 国产精品一及| 丝袜人妻中文字幕| 国内精品久久久久精免费| 国产精品乱码一区二三区的特点| 亚洲精品在线观看二区| 国产精品av久久久久免费| 亚洲国产精品久久男人天堂| 成人av一区二区三区在线看| 免费看美女性在线毛片视频| 国产av一区在线观看免费| 日本免费一区二区三区高清不卡| 日本a在线网址| 正在播放国产对白刺激| 国产三级中文精品| 夜夜爽天天搞| 国产99久久九九免费精品| 欧美日本亚洲视频在线播放| 18美女黄网站色大片免费观看| 看片在线看免费视频| 国产一区在线观看成人免费| 日韩有码中文字幕| 日韩欧美 国产精品| 亚洲一卡2卡3卡4卡5卡精品中文| 悠悠久久av| 久久午夜亚洲精品久久| ponron亚洲| 亚洲人成电影免费在线| 国产精品影院久久| 97超级碰碰碰精品色视频在线观看| av天堂在线播放| av片东京热男人的天堂| 18禁黄网站禁片免费观看直播| 老司机午夜福利在线观看视频| 国产在线观看jvid| 亚洲国产精品久久男人天堂| 久久精品影院6| 国产探花在线观看一区二区| 亚洲精品国产精品久久久不卡| 一级毛片精品| 舔av片在线| 久久精品国产综合久久久| or卡值多少钱| 国产成+人综合+亚洲专区| 午夜亚洲福利在线播放| 成年版毛片免费区| 高清毛片免费观看视频网站| 老汉色∧v一级毛片| 麻豆国产av国片精品| 免费无遮挡裸体视频| 美女黄网站色视频| 亚洲一区高清亚洲精品| 国产亚洲av高清不卡| 精品免费久久久久久久清纯| 搡老熟女国产l中国老女人| 国内毛片毛片毛片毛片毛片| 在线十欧美十亚洲十日本专区| 成人午夜高清在线视频| 色老头精品视频在线观看| 欧美黑人巨大hd| 女同久久另类99精品国产91| 麻豆国产av国片精品| 999久久久国产精品视频| 久久久精品国产亚洲av高清涩受| 日韩成人在线观看一区二区三区| 99热这里只有精品一区 | 精品人妻1区二区| 性欧美人与动物交配| 美女午夜性视频免费| 校园春色视频在线观看| 黄片小视频在线播放| 99热这里只有是精品50| 欧美性猛交黑人性爽| 亚洲精品av麻豆狂野| 亚洲av电影在线进入| 亚洲专区字幕在线| 老鸭窝网址在线观看| 国产亚洲精品久久久久久毛片| 制服诱惑二区| 99国产极品粉嫩在线观看| 精品无人区乱码1区二区| 午夜影院日韩av| 成年女人毛片免费观看观看9| 巨乳人妻的诱惑在线观看| 久久伊人香网站| 宅男免费午夜| 亚洲狠狠婷婷综合久久图片| 一本综合久久免费| 欧美黑人巨大hd| 久久精品91蜜桃| 操出白浆在线播放| 日韩欧美精品v在线| 99热这里只有精品一区 | 久久久久久国产a免费观看| 亚洲国产日韩欧美精品在线观看 | 国产三级黄色录像| 又黄又粗又硬又大视频| 一级片免费观看大全| 精品国产亚洲在线| 日本一区二区免费在线视频| 搡老岳熟女国产| www.熟女人妻精品国产| 天堂动漫精品| 看片在线看免费视频| 亚洲中文av在线| 韩国av一区二区三区四区| 亚洲一码二码三码区别大吗| 国产精品1区2区在线观看.| 99久久无色码亚洲精品果冻| 国产精品亚洲美女久久久| 欧美三级亚洲精品| 少妇人妻一区二区三区视频| 99在线视频只有这里精品首页| 一二三四在线观看免费中文在| 天天添夜夜摸| 日日爽夜夜爽网站| 国产欧美日韩精品亚洲av| 精品第一国产精品| 亚洲美女黄片视频| 亚洲欧美日韩高清专用| 精品少妇一区二区三区视频日本电影| 最近视频中文字幕2019在线8| 高潮久久久久久久久久久不卡| 婷婷六月久久综合丁香| 可以免费在线观看a视频的电影网站| 婷婷精品国产亚洲av在线| 国产精品美女特级片免费视频播放器 | 在线看三级毛片| av天堂在线播放| 亚洲欧美激情综合另类| 神马国产精品三级电影在线观看 | 老汉色∧v一级毛片| 国产成人影院久久av| 激情在线观看视频在线高清| 亚洲七黄色美女视频| 亚洲av片天天在线观看| 一个人免费在线观看的高清视频| 国产探花在线观看一区二区| 五月玫瑰六月丁香| 久久久久久久久免费视频了| 美女黄网站色视频| 一本精品99久久精品77| 亚洲激情在线av| 黄色丝袜av网址大全| 在线永久观看黄色视频| 亚洲第一电影网av| 香蕉av资源在线| 制服人妻中文乱码| 国产精品久久久人人做人人爽| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产看品久久| 久久精品国产99精品国产亚洲性色| 亚洲欧美精品综合一区二区三区| 一级片免费观看大全| 国产高清视频在线播放一区| 老司机福利观看| 波多野结衣高清作品| 亚洲专区字幕在线| av超薄肉色丝袜交足视频| 又紧又爽又黄一区二区| 露出奶头的视频| 1024视频免费在线观看| 日日爽夜夜爽网站| 岛国视频午夜一区免费看| 久久热在线av| 国产精品,欧美在线| 一本精品99久久精品77| 欧美黄色片欧美黄色片| 国产精品一区二区免费欧美| 18禁国产床啪视频网站| 精品久久久久久久久久久久久| 精品国产美女av久久久久小说| 国产人伦9x9x在线观看| 免费在线观看完整版高清| 日本在线视频免费播放| 国产欧美日韩一区二区三| 久久性视频一级片| 男男h啪啪无遮挡| 天天添夜夜摸| 婷婷亚洲欧美| 亚洲精品久久国产高清桃花| 成人精品一区二区免费| 麻豆国产97在线/欧美 | 特级一级黄色大片| 69av精品久久久久久| www日本在线高清视频| 欧美zozozo另类| 变态另类成人亚洲欧美熟女| 久久午夜亚洲精品久久| 亚洲全国av大片| 久久久久国产一级毛片高清牌| 丰满的人妻完整版| 免费在线观看日本一区| 床上黄色一级片| 国产亚洲精品av在线| 一个人免费在线观看电影 | 国产精品久久久久久久电影 | 日韩欧美国产一区二区入口| 中出人妻视频一区二区| www.999成人在线观看| 国内精品一区二区在线观看| 久99久视频精品免费| 亚洲美女视频黄频| 9191精品国产免费久久| 高清毛片免费观看视频网站| 日韩成人在线观看一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 成人午夜高清在线视频| 亚洲全国av大片| 黄色丝袜av网址大全| 99国产综合亚洲精品| av在线天堂中文字幕| 小说图片视频综合网站| 欧美3d第一页| 国产av不卡久久| 午夜激情av网站| 久久久久性生活片| 在线永久观看黄色视频| 国产伦一二天堂av在线观看| 好男人在线观看高清免费视频| 色噜噜av男人的天堂激情| 成人精品一区二区免费| 高清毛片免费观看视频网站| 美女大奶头视频| 亚洲精品av麻豆狂野| 国产精品久久视频播放| 色在线成人网| 脱女人内裤的视频| 久久精品91无色码中文字幕| 久久人人精品亚洲av| 免费在线观看日本一区| a级毛片在线看网站| 2021天堂中文幕一二区在线观| 午夜精品一区二区三区免费看| 免费人成视频x8x8入口观看| 久久久久亚洲av毛片大全| 亚洲人成网站高清观看| 欧美在线一区亚洲| 日本五十路高清| 午夜精品一区二区三区免费看| 色播亚洲综合网| 精品人妻1区二区| 欧美一级a爱片免费观看看 | 少妇粗大呻吟视频| www.熟女人妻精品国产| 欧美日韩乱码在线| 久久久久久九九精品二区国产 | 亚洲无线在线观看| 日本黄色视频三级网站网址| 老司机福利观看| 国产一区二区三区在线臀色熟女| 人妻夜夜爽99麻豆av| 亚洲一区高清亚洲精品| 国产黄a三级三级三级人| 久久国产精品影院| 国产亚洲欧美98| 欧美黑人精品巨大| 亚洲欧美日韩高清在线视频| 一进一出抽搐gif免费好疼| 久久中文字幕人妻熟女|