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

    我國(guó)東部夏季一次強(qiáng)對(duì)流活動(dòng)過(guò)程中對(duì)流層上部大氣成分變化的分析

    2012-12-14 11:19:50易明建傅云飛劉鵬衡志煒高越洪星園
    大氣科學(xué) 2012年5期
    關(guān)鍵詞:平流層對(duì)流層強(qiáng)對(duì)流

    易明建 傅云飛 劉鵬 衡志煒 高越 洪星園

    中國(guó)科學(xué)技術(shù)大學(xué)地球和空間科學(xué)學(xué)院,合肥 230026

    我國(guó)東部夏季一次強(qiáng)對(duì)流活動(dòng)過(guò)程中對(duì)流層上部大氣成分變化的分析

    易明建 傅云飛 劉鵬 衡志煒 高越 洪星園

    中國(guó)科學(xué)技術(shù)大學(xué)地球和空間科學(xué)學(xué)院,合肥 230026

    利用星載微波臨邊遙感探測(cè)結(jié)果,對(duì)2006年6月28~29日江淮地區(qū)的一次強(qiáng)對(duì)流天氣過(guò)程中對(duì)流層上部一氧化碳(CO)、臭氧(O3)、水汽(H2O)和冰水含量(IWC)的分布特點(diǎn)進(jìn)行了研究。強(qiáng)對(duì)流天氣過(guò)程前后的對(duì)比分析表明,CO混合比增大,在200 hPa處增加了0.12 ppm (1 ppm=10-6);O3混合比減小,在70 hPa 處減少了0.30 ppm;H2O混合比在250 hPa處增加了400 ppm;IWC在強(qiáng)降水發(fā)生之前有大幅增長(zhǎng),在200 hPa處最大含量可達(dá)0.03 g/m3。CO和O3含量與垂直運(yùn)動(dòng)速度兩者的相關(guān)變化表明,對(duì)流垂直輸送作用可能是造成對(duì)流層上層和平流層低層大氣成分變化的機(jī)制之一。而H2O和IWC含量的增加主要局限于對(duì)流層頂以下,這表明對(duì)流層上部水物質(zhì)的質(zhì)量和形態(tài)是由垂直輸送作用和對(duì)流系統(tǒng)內(nèi)部的微物理過(guò)程共同決定的。

    對(duì)流 大氣成分 微波臨邊遙感

    1 引言

    對(duì)流層頂是地球大氣對(duì)流層與平流層的過(guò)渡層,也是對(duì)流層大氣與平流層大氣進(jìn)行物質(zhì)和能量直接交換的關(guān)鍵區(qū)域,這里的大氣微量成分不僅是對(duì)流層與平流層之間動(dòng)力輸送的示蹤物,還是影響光化學(xué)過(guò)程和大氣輻射平衡的關(guān)鍵因子(Sherwood and Dessler, 2000; Salby et al., 2003; 陳洪濱等, 2006;卞建春, 2009)。在熱帶地區(qū),對(duì)流層頂高度變化在14 km至18 km之間,該地區(qū)大部分的對(duì)流活動(dòng)位于對(duì)流層頂以下,但也有一些強(qiáng)對(duì)流活動(dòng)中的上升氣流可以穿過(guò)對(duì)流層頂,進(jìn)入到平流層低層(Kuhn et al., 1971; Kuhn and Stearns, 1973; Alcala and Dessler, 2002)。如,Liu and Zipser(2005)通過(guò)對(duì)5年時(shí)間 TRMM(Tropical Rainfall Measuring Mission)資料的統(tǒng)計(jì)發(fā)現(xiàn),約有 1.3%的熱帶對(duì)流系統(tǒng)可以達(dá)到14 km高度,其中陸地上深對(duì)流發(fā)生的頻率要比洋面高,特別是在中部非洲等地區(qū),深對(duì)流發(fā)生的可能性更大。

    近年來(lái)的研究結(jié)果表明,對(duì)流活動(dòng)中的向上輸送及其混合過(guò)程是影響對(duì)流層頂水汽變化的重要過(guò)程,特別是冰相粒子的向上輸送對(duì)于對(duì)流層頂水汽變化有著更重要的作用(Sassi et al., 2001; Hassim and Lane, 2010)。其實(shí),對(duì)流活動(dòng)不只是引起對(duì)流層頂水汽的變化,還引起對(duì)流層頂其它的一些微量氣體(如CO、O3等)分布的變化。已有研究結(jié)果表明,對(duì)流層上部(即370 K等位溫面以上)多達(dá)60%的空氣可以通過(guò)平流層底部(即380 K等位溫面)向上逸出(Dessler, 2002),而在對(duì)流層低層排放的微量氣體(如 CO和CH4等),也會(huì)由于對(duì)流活動(dòng)的垂直輸送,上升至對(duì)流層頂附近甚至進(jìn)入到平流層(Ricaud et al., 2007; 卞建春等, 2011)。

    迄今為止,夏季熱帶對(duì)流及其影響仍是研究的熱點(diǎn),而對(duì)熱帶外地區(qū)強(qiáng)對(duì)流活動(dòng)和相應(yīng)的大氣成分變化的研究比較少。但是,對(duì) TRMM(Tropical Rainfall Measuring Mission satellite)測(cè)雨雷達(dá)十年探測(cè)結(jié)果的統(tǒng)計(jì)表明,夏季亞洲副熱帶地區(qū)的對(duì)流降水頂部高度高于熱帶地區(qū),尤其在陸面上更是如此,由此推測(cè),夏季整個(gè)副熱帶地區(qū),特別是副熱帶陸面,對(duì)流將可能強(qiáng)于熱帶地區(qū)(傅云飛等, 2008;劉鵬和傅云飛, 2010)。實(shí)際上,利用機(jī)載設(shè)備進(jìn)行的觀測(cè)實(shí)驗(yàn)已經(jīng)證實(shí),在熱帶外地區(qū)發(fā)生的強(qiáng)對(duì)流活動(dòng)中,對(duì)流層低層空氣可以隨強(qiáng)烈對(duì)流運(yùn)動(dòng)上升到接近平流層的高度(Fischer et al., 2003),在對(duì)流云頂部也發(fā)現(xiàn)CO的濃度要比對(duì)流云外部高出近一倍(Poulida et al., 1996)。由此可知,熱帶外地區(qū)強(qiáng)對(duì)流活動(dòng)對(duì)大氣成分變化的影響也非常顯著。本文將利用多平臺(tái)高分辨率的衛(wèi)星資料,就 2006年 6月28日~29日發(fā)生在江淮之間的一次強(qiáng)對(duì)流活動(dòng)前后對(duì)流層上層大氣成分變化進(jìn)行研究,并討論對(duì)流活動(dòng)對(duì)大氣成分垂直輸送的作用。

    2 資料介紹

    本文主要使用了TRMM PR、CloudSat CPR、Aura MLS等衛(wèi)星資料和NCEP資料集。

    TRMM是一顆非太陽(yáng)同步極軌衛(wèi)星,軌道傾角約35°,飛行高度350 km(2001 年8月7日后改為400 km)(Kummerow et al., 1998, 2000)。其環(huán)繞地球一周約需96 min,每天在38°S至38°N之間約有16條軌道。文中使用的標(biāo)準(zhǔn)資料2A25是根據(jù)PR回波反演的三維降水率(mm/h)(Iguchi et al., 2000;Meneghini et al., 2000)。PR的星下點(diǎn)分辨率是4.3 km(水平)、0.25 km(垂直),垂直探測(cè)高度自地表至20 km高,掃描寬度約為220 km。

    CloudSat是太陽(yáng)同步衛(wèi)星,其軌道高度為705 km,軌道面和平均太陽(yáng)子午線(xiàn)之間的固定傾角為98.2°(Stephens et al., 2002; Haynes and Stephens,2007; Stephens and Kummerow, 2007; Stephens et al.,2008)。CloudSat衛(wèi)星搭載的CPR(Cloud Profiling Radar)每 0.16秒對(duì)大氣發(fā)射一次探測(cè)脈沖,衛(wèi)星在軌道上向前飛行1.1 km后接收大氣返回脈沖,這樣就對(duì)大氣完成一次垂直探測(cè)。CPR在垂直方向共探測(cè)有125層,每層分辨率為240 m,共對(duì)30 km厚度的大氣層進(jìn)行探測(cè)。衛(wèi)星繞軌道 1周約需 99分鐘,每天共有 14條軌道,每條軌道上可以產(chǎn)生37088個(gè)云剖面(馬占山等, 2008)。文中所使用的雷達(dá)反射率因子有效數(shù)值范圍為-40~50 dBZ。

    Aura衛(wèi)星的軌道高度為705 km,與CloudSat同屬于A-Train系列,由于兩者之間相互距離非常近,它們通過(guò)探測(cè)地的時(shí)間差幾乎可以忽略,極大的方便了衛(wèi)星資料的綜合應(yīng)用。MLS(Micro wave Limb Sounder)(Waters et al., 2006)是 Aura衛(wèi)星的重要探測(cè)器之一,主要用于測(cè)量地球大氣臨邊的微波熱輻射,通過(guò)臨邊掃描來(lái)觀測(cè)從平流層到對(duì)流層頂?shù)?18、190、240、640 GHz和2.5 THz光譜范圍的微波散射。MLS每24.7秒完成一次從地表到約90 km高度的掃描,每條廓線(xiàn)沿衛(wèi)星軌跡的空間間隔為1.5°(約165 km),垂直分辨率是3 km左右,準(zhǔn)全球覆蓋(南北緯 82°之間)。MLS資料既可用于反演 O3、H2O、BrO、ClO、HCl、HOCl、OH、HO2、HCN、CO、HNO3、N2O和SO2混合比的垂直廓線(xiàn),也可用于反演冰面相對(duì)濕度、云含冰量、位勢(shì)高度和溫度。MLS資料早期版本(v1.5)與氣球探空資料的差別一般在5%~10%(Froidevaux et al., 2006),而本文使用的v2.2版資料已經(jīng)針對(duì)前一版本存在的問(wèn)題作了一些修正和改進(jìn),其準(zhǔn)確性和可信度進(jìn)一步得到提高(Livesey et al., 2007)。

    NCEP(National Centers for Environmental Prediction)資料為 17層氣壓坐標(biāo)系再分析資料(Kalnay et al., 1996; Kistler et al., 2001),水平分辨率為 2.5°×2.5°。本文使用了每日4次的位勢(shì)高度、溫度、相對(duì)濕度、以及水平和垂直速度物理量來(lái)計(jì)算和分析大氣環(huán)流形勢(shì)。

    3 天氣背景

    3.1 降水和對(duì)流云特征

    本文研究的強(qiáng)對(duì)流過(guò)程發(fā)生于2006年6月28日~29日江淮之間。這次過(guò)程中安徽北部和中部均出現(xiàn)了強(qiáng)雷暴、大風(fēng)(極值 32 m/s)、龍卷和短時(shí)暴雨。圖1給出了雨量計(jì)觀測(cè)的六小時(shí)累計(jì)降水,降水主要集中在2006年6月29日00時(shí)到06時(shí)(協(xié)調(diào)世界時(shí),下同)之間,降水發(fā)生的主要范圍在(31°N~35°N, 115°E~121°E)(圖 1 中方框所示,稱(chēng)之為I區(qū),下文均與此相同),六小時(shí)累計(jì)降水最大可達(dá)40 mm以上。TRMM PR反演的地表降水分布(如圖 2)顯示,從安徽到江蘇北部為一狹長(zhǎng)雨帶,寬約100 km,長(zhǎng)達(dá)數(shù)百公里。

    通過(guò)對(duì)地面觀測(cè)結(jié)果和衛(wèi)星反演結(jié)果的比較,已經(jīng)很清楚地描述了強(qiáng)對(duì)流天氣過(guò)程的降水分布和降水強(qiáng)度變化。CloudSat衛(wèi)星在29日16:42經(jīng)過(guò)Ⅰ區(qū)上空,從圖3給出的CloudSat軌道和探測(cè)的雷達(dá)反射率因子可知,有一深對(duì)流云團(tuán)中心位于(33.7°N, 118.2°E)處,云頂高度約14 km,云砧明顯。該對(duì)流云出現(xiàn)在地面最強(qiáng)降水過(guò)程之后,但是其位置和垂直高度仍然代表了該地區(qū)發(fā)生的強(qiáng)對(duì)流云團(tuán)的一般特點(diǎn)。按照對(duì)流云頂高度劃分,這是一次深對(duì)流過(guò)程(云頂高度大于 10 km)。由于CloudSat相鄰軌道在地球表面的東西向間距約為24°,且軌道不會(huì)連續(xù)經(jīng)過(guò)同一地點(diǎn),因此在 6月28~29日這段時(shí)間里,只能得到這一例對(duì)流云圖像。在江淮地區(qū)的這次強(qiáng)對(duì)流天氣過(guò)程中,是否還存在有更強(qiáng)更高的對(duì)流云,單靠目前所掌握的衛(wèi)星資料還不能予以證實(shí),但是如果以降水強(qiáng)度的變化推測(cè),在29日凌晨極有可能出現(xiàn)比圖3b中更強(qiáng)更高的對(duì)流云。

    3.2 大氣環(huán)流背景

    強(qiáng)對(duì)流和暴雨天氣的發(fā)生一般滿(mǎn)足以下物理?xiàng)l件:位勢(shì)不穩(wěn)定層結(jié),低層水汽輻合,具有使不穩(wěn)定能量釋放的機(jī)制。下面將對(duì)這次強(qiáng)對(duì)流天氣過(guò)程中的一些物理?xiàng)l件進(jìn)行簡(jiǎn)單的分析,以了解強(qiáng)對(duì)流天氣的大氣環(huán)流特點(diǎn)和垂直對(duì)流的成因。

    在江淮地區(qū)強(qiáng)對(duì)流天氣之前,500 hPa高度場(chǎng)上可以看到高緯度的低壓槽分裂出短波槽東移南下,槽后冷空氣與西南氣流在江淮地區(qū)交匯。圖 4給出了6月29日00時(shí)850 hPa、500 hPa位勢(shì)高度和風(fēng)場(chǎng)。從850 hPa(圖4a)可以看到,低空有一支偏南方向的低空急流,這支急流在低空輸送暖濕空氣;從500 hPa(圖4b)可以看到,高空為一支西北方向急流和冷平流區(qū)。高低空大氣環(huán)流形勢(shì)有利于強(qiáng)對(duì)流發(fā)展:其一,產(chǎn)生上下層溫度平流差,形成大氣不穩(wěn)定層結(jié),提供了強(qiáng)對(duì)流和暴雨所需的必要條件;其二,低空急流中心的左前方是低空輻合區(qū),高空急流中心的右后方和左前方是高空輻散區(qū)。當(dāng)高低空急流交匯,低空急流的輻合區(qū)與高空急流的輻散區(qū)重合時(shí),造成上升運(yùn)動(dòng)發(fā)展,觸發(fā)潛在不穩(wěn)定的釋放。

    首先,強(qiáng)對(duì)流天氣發(fā)生前后大氣不穩(wěn)定層結(jié)發(fā)生了明顯變化。大氣穩(wěn)定程度常用沙氏指數(shù)(Is)來(lái)表示,它定義為850 hPa等壓面上的濕空氣塊沿干絕熱線(xiàn)抬升到達(dá)凝結(jié)高度后,再沿濕絕熱線(xiàn)上升至500 hPa時(shí)該氣塊具有的溫度T'與500 hPa等壓面上的環(huán)境溫度T500的差值(Showalter, 1953)。Is單位為℃,Is數(shù)值大小與對(duì)流天氣有密切關(guān)系。通常,當(dāng)Is<0時(shí),大氣層結(jié)不穩(wěn)定,且負(fù)值越大,不穩(wěn)定程度越大,反之,則表示氣層是穩(wěn)定的。當(dāng)-3<Is<0時(shí),表示有發(fā)生雷暴的可能性。從表1可以看到,從28日12時(shí)至29日06時(shí),沙氏指數(shù)連續(xù)為負(fù)值,它反映了在此期間,對(duì)流層大氣處于對(duì)流不穩(wěn)定狀態(tài),若有觸發(fā)機(jī)制,極其容易使不穩(wěn)定能量爆發(fā)。在29日00時(shí),沙氏指數(shù)達(dá)到負(fù)的最大。隨后,發(fā)生強(qiáng)對(duì)流和大量降水使得大氣不穩(wěn)定能量逐漸釋放,到29日12時(shí),沙氏指數(shù)回到正值,大氣恢復(fù)穩(wěn)定狀態(tài)。

    圖1 2006年6月雨量計(jì)觀測(cè)的六小時(shí)累計(jì)降水:(a) 28日12~18時(shí);(b) 28日18時(shí)~29日00時(shí);(c) 29日00~06時(shí);(d) 29日06~12時(shí)Fig.1 Accumulative six-hour precipitation by hyetometer in Jun 2006: (a) 1200-1800 UTC 28 Jun; (b) 1800 UTC 28 Jun-0000 UTC 29 Jun; (c) 0000 UTC-0600 UTC 29 Jun; (d) 0600 UTC-1200 UTC 29 Jun

    圖2 TRMM探測(cè)的地表降水Fig.2 Surface rain distribution measured by TRMM

    此外,大氣強(qiáng)迫抬升運(yùn)動(dòng)是大氣不穩(wěn)定得以釋放并造成強(qiáng)對(duì)流和暴雨天氣的直接機(jī)制。圖5給出了6月29日00時(shí)850 hPa、200 hPa的散度場(chǎng)。與圖 4中位勢(shì)高度和風(fēng)場(chǎng)相對(duì)應(yīng),在江淮地區(qū),850 hPa(圖5a)低空為大氣輻合運(yùn)動(dòng),200 hPa(圖5b)高空為大氣輻散運(yùn)動(dòng)。低空輻合區(qū)與高空輻散區(qū)重合,空氣受迫抬升引發(fā)垂直運(yùn)動(dòng)。圖6給出了29日00時(shí)500 hPa垂直速度ω以及I區(qū)平均的ω隨時(shí)間和高度的變化。從圖6a可以看到,江淮地區(qū)500 hPa等壓面上大氣有較強(qiáng)的上升運(yùn)動(dòng),垂直速率為 0.2 Pa/s。事實(shí)上,江淮地區(qū)對(duì)流層整層大氣上升運(yùn)動(dòng)從28日就已經(jīng)開(kāi)始,到29日凌晨時(shí)從900 hPa一直到200 hPa高度上均為上升運(yùn)動(dòng),最強(qiáng)上升運(yùn)動(dòng)在400 hPa附近,速率大于0.2 Pa/s(圖6b)。此時(shí)正對(duì)應(yīng)于江淮地區(qū)發(fā)生最強(qiáng)降水過(guò)程。從29日06時(shí)開(kāi)始,整個(gè)對(duì)流層大氣轉(zhuǎn)為下沉運(yùn)動(dòng),強(qiáng)對(duì)流天氣過(guò)程也由此開(kāi)始減弱。

    圖3 (a)CloudSat CPR探測(cè)軌道和(b)雷達(dá)反射率因子剖面圖Fig.3 (a) CloudSat CPR orbit and (b) profile of radar reflectivity factor

    圖4 2006年6月29日00時(shí)位勢(shì)高度(單位:gpm)和風(fēng)場(chǎng)(單位:m/s):(a)850 hPa;(b)500 hPaFig.4 Geopotential height fileds (units: gpm) and horizontal wind at 0000 UTC Jun 29 2006 at (a) 850 hPa and (b) 500 hPa

    4 對(duì)流層上層大氣成分分布

    在這一次強(qiáng)對(duì)流事件中,搭載微波臨邊遙感(MLS) 的 Aura衛(wèi)星分別于 28日 05:29和 29日18:15經(jīng)過(guò)I區(qū)上空(圖7),從地面觀測(cè)的累計(jì)降水隨時(shí)間變化(圖1)可知,這兩個(gè)時(shí)刻分別是在強(qiáng)對(duì)流降水發(fā)生之前和之后。圖8和圖9分別給出了沿這兩條軌道A和B的溫度(T)、一氧化碳(CO)、臭氧(O3)、對(duì)流層水汽(H2O)、平流層水汽(H2O)和冰水含量(IWC, ice water content)的分布。選擇這幾種大氣成分來(lái)分析它們?cè)趶?qiáng)對(duì)流過(guò)程中的變化,有著各自不同的代表意義。CO和O3的變化直接受動(dòng)力輸送的影響,而與對(duì)流過(guò)程中溫度變化等物理過(guò)程無(wú)關(guān),因此可以作為大氣垂直運(yùn)動(dòng)的示蹤氣體。CO源地是在對(duì)流層低層,隨著高度增加到對(duì)流層頂附近混合比逐漸減小,對(duì)流上升運(yùn)動(dòng)將會(huì)導(dǎo)致下層高混合比 CO進(jìn)入上層,可以使得上層混合比增加。O3最大混合比一般是在10 hPa附近, 從10 hPa以下直到300 hPa的對(duì)流層上部,越向下混合比逐漸減小,對(duì)流上升運(yùn)動(dòng)將會(huì)導(dǎo)致下層低混合比的O3進(jìn)入上層,可以使得上層混合比降低。而H2O除了同樣受到垂直運(yùn)動(dòng)輸送的影響以外,還與強(qiáng)對(duì)流過(guò)程中隨溫度變化的微物理過(guò)程密切相關(guān)。大部分的 H2O在經(jīng)過(guò)冷的對(duì)流層頂前已經(jīng)凝結(jié)沉降,平流層的H2O混合比比對(duì)流層要小得多。IWC反映了大氣中冰相水的含量,是對(duì)流層上層大氣中水物質(zhì)組成和變化另一個(gè)重要因子。

    圖5 2006年6月29日00時(shí)的散度(單位:10-5 s-1):(a) 850 hPa ; (b) 200 hPaFig.5 Divergence fields (units: 10-5 s-1) at 0000 UTC Jun 29 2006 at (a) 850 hPa and (b) 200 hPa

    圖6 (a)2006年6月29日00時(shí)500 hPa垂直速度(ω)分布(單位:Pa/s);(b)I區(qū)平均的垂直速度(ω)隨時(shí)間和高度的變化(單位:Pa/s)Fig.6 (a) Vertical velocity (ω) at 500 hPa at 0000 UTC Jun 29 2006; (b) time series of average vertical velocity in region I (ω) from 1000 hPa to 100 hPa.(units:Pa/s)

    圖7 Aura衛(wèi)星軌道Fig.7 Orbit of Aura satellite

    表1 沙氏指數(shù)(Is)變化Table 1 Variation of Showalter index (Is)

    從圖8中可以看到,對(duì)流層上層有較強(qiáng)的上升運(yùn)動(dòng),I區(qū)上空垂直速度(ω,黑色箭頭)比周邊地區(qū)要大,300 hPa附近速度最大可達(dá)0.2 Pa/s。這個(gè)時(shí)候雖然強(qiáng)降水還沒(méi)有發(fā)生,但是空氣受迫抬升運(yùn)動(dòng)已經(jīng)發(fā)展到一定程度。表明這一次強(qiáng)對(duì)流天氣過(guò)程中,強(qiáng)降水的發(fā)生對(duì)應(yīng)著對(duì)流活動(dòng)發(fā)展達(dá)到頂峰,而對(duì)流層大氣的垂直上升運(yùn)動(dòng)在此前十幾個(gè)小時(shí)就已經(jīng)開(kāi)始了。圖8還給出了此時(shí)的溫度和這幾種大氣成分在不同高度的混合比分布特點(diǎn),CO混合比在I區(qū)出現(xiàn)明顯擾動(dòng),在南部邊緣CO混合比出現(xiàn)了顯著的增加。O3混合比從200 hPa到50 hPa逐漸增加,水平方向從南往北混合比增加,100 hPa以下O3混合比都非常小。在150 hPa附近,I區(qū)偏北的位置上有一個(gè)很小的區(qū)域,CO混合比減小并且O3混合比增加,這與上升氣流的作用結(jié)果相反。有可能是對(duì)流層頂附近中小尺度的瞬變擾動(dòng)造成了這種現(xiàn)象,但是需要采用更高分辨率的探空結(jié)果作進(jìn)一步的分析才可能給予證實(shí)。因?yàn)檫@并不是本文的重點(diǎn)和關(guān)鍵所在,這里不再詳述。H2O混合比在對(duì)流層上層變化也非常明顯,250 hPa以下混合比大于10-4量級(jí),I區(qū)附近H2O混合比值最大可達(dá)0.005,在25 0hPa以上混合比變得很小。平流層低層H2O混合比小于10-5量級(jí),I區(qū)上空相比周邊地區(qū)差別很小。IWC分布的特點(diǎn)就很明顯,在I區(qū)要比周邊區(qū)域大得多,200 hPa高度處大約為 0.035 g/m3。對(duì)流系統(tǒng)頂部大量冰晶粒子的出現(xiàn),意味著降水粒子的大量形成和增長(zhǎng),這正是地面強(qiáng)降水啟動(dòng)和維持所必需的前提條件。

    圖9是沿軌道B的大氣成分剖面圖。此時(shí)地面強(qiáng)降水過(guò)程已經(jīng)結(jié)束,圖中也可以看到對(duì)流層上層垂直速度減小了很多。與圖8中的幾種大氣成分混合比進(jìn)行比較發(fā)現(xiàn):I區(qū)上空CO混合比增加,200 hPa高度上從 0.06×10-6增加到了 0.18×10-6,增加了一倍多,這種增加一直到100 hPa以上都是非常明顯的。O3混合比濃度由于在對(duì)流層上部的值很小,與之前相比差別不明顯,但是水平方向的分布出現(xiàn)了明顯擾動(dòng),I區(qū)上空及其周邊O3混合比要比同一高度上的其他地區(qū)略小。H2O混合比在300 hPa高度上有些減少,但往上卻有明顯增加,在250 hPa處混合比從0.0001增加到0.0005以上。在150 hPa的平流層低層,低緯地區(qū)H2O混合比大,I區(qū)H2O混合比變得很小,甚至不到 1.0×10-6。IWC在 I區(qū)上空仍然比其他緯度上要大,不過(guò)比起之前已經(jīng)減少了很多,200 hPa高度處不到0.01 g/m3,下降了一半以上。

    通過(guò)上述分析可知,對(duì)流層上層大氣成分在強(qiáng)對(duì)流前后發(fā)生了一些變化。通過(guò)對(duì)圖7中Ⅰ區(qū)內(nèi)的衛(wèi)星觀測(cè)數(shù)據(jù)逐軌平均,并按照時(shí)間坐標(biāo)插值到每日00時(shí),得到了6月26日至7月3日之間這幾種大氣成分區(qū)域平均的混合比隨時(shí)間的變化,同時(shí)配以垂直速度的變化(圖10)。從中可以看到,在27日和 28日已經(jīng)開(kāi)始零星出現(xiàn)有上升運(yùn)動(dòng),但持續(xù)時(shí)間很短。到 29日凌晨前后,持續(xù)的上升運(yùn)動(dòng)開(kāi)始了。之后的30日和7月2日垂直速度也較大,是因?yàn)楫?dāng)?shù)赜职l(fā)生了兩次對(duì)流過(guò)程(石瑩等,2009)。在28日之后對(duì)流運(yùn)動(dòng)發(fā)展的同時(shí),對(duì)流層上部與平流層低層之間的過(guò)渡層也變得更冷,對(duì)流層上層大氣成分分布開(kāi)始出現(xiàn)明顯變化。雖然這比地面強(qiáng)降水發(fā)生的時(shí)間要略早一些,但是與垂直上升運(yùn)動(dòng)開(kāi)始的時(shí)間是一致的。大氣成分分布的變化也很好地印證了上升運(yùn)動(dòng)的作用結(jié)果:在 200 hPa一直到50 hPa的高度上,CO混合比都表現(xiàn)出了顯著的增加;O3混合比在100 hPa以上出現(xiàn)了一定的減少;H2O混合比在對(duì)流層上部(250 hPa)增加的同時(shí)在對(duì)流層頂附近(150 hPa)附近減少;IWC的變化尤其明顯,在地面強(qiáng)降水發(fā)生之前不斷增加,之后便迅速減少。

    5 結(jié)論

    本文利用TRMM PR,CloudSat CPR,Aura MLS和NCEP資料,對(duì)2006年6月28~29日 (UTC) 發(fā)生在江淮地區(qū)的一次強(qiáng)對(duì)流天氣中的降水分布和對(duì)流云結(jié)構(gòu),大氣環(huán)流背景,以及對(duì)流層上層大氣成分分布進(jìn)行了分析。目的在于揭示夏季熱帶外地區(qū)強(qiáng)對(duì)流特征,強(qiáng)對(duì)流形成和發(fā)展的天氣背景,以及對(duì)流垂直運(yùn)動(dòng)與對(duì)流層上層大氣成分分布變化的關(guān)系,為進(jìn)一步計(jì)算強(qiáng)對(duì)流對(duì)大氣成分在對(duì)流層與平流層之間的直接輸送和研究對(duì)流層與平流層物質(zhì)交換的物理機(jī)制提供前期基礎(chǔ)。

    圖8 6月28日05時(shí)29分溫度、大氣成分分布與垂直速度:(a)T;(b)CO;(c)O3;(d)對(duì)流層H2O;(e)平流層H2O;(f)IWC。矢量箭頭:垂直速度ω (單位:Pa/s);虛線(xiàn):I區(qū)范圍Fig.8 The profiles of temperature, atmospheric composition, and vertical velocity ω (arrows) at 0529 UTC 28 Jun: (a) T; (b) CO; (c) O3; (d) H2O in the troposphere; (e) H2O in the stratosphere; (f) IWC. Region between the two white dashed lines: region I

    圖9 同圖8,但為6月29日18:15Fig.9 Same as Fig. 8, but for 1815 UTC 29 Jun

    圖10 對(duì)流發(fā)生前后溫度、大氣成分與垂直速度的變化。兩條白虛線(xiàn)代表6月28日05時(shí)和6月29日18時(shí)Fig.10 The variation of temperature, atmospheric composition, and vertical velocity in and near the convection process. The two white dashed lines indicate the time: at 0500 UTC 28 Jun and at 1800 UTC 29 Jun

    研究結(jié)果表明,最強(qiáng)地面降水出現(xiàn)在 29日凌晨前后,此時(shí)對(duì)應(yīng)有最強(qiáng)的大氣上升運(yùn)動(dòng)。測(cè)云雷達(dá)探測(cè)到了強(qiáng)對(duì)流減弱階段的一個(gè)對(duì)流云剖面,雷達(dá)反射率剖面顯示云頂高度達(dá)到14 km,深入到了對(duì)流層上層。對(duì)天氣背景的分析表明,高低空大氣環(huán)流形勢(shì)提供了有利于強(qiáng)對(duì)流天氣發(fā)展的條件,對(duì)流層大氣不穩(wěn)定狀態(tài)持續(xù)發(fā)展,低空輻合區(qū)與高空輻散區(qū)相重合促使大氣受迫抬升,觸發(fā)不穩(wěn)定能量釋放,最終造成大氣的強(qiáng)烈對(duì)流運(yùn)動(dòng)和強(qiáng)降水的發(fā)生。

    MLS探測(cè)結(jié)果表明,強(qiáng)對(duì)流天氣過(guò)程空氣上升運(yùn)動(dòng)速率變化與對(duì)流層上層和平流層低層的大氣成分分布變化有很好的對(duì)應(yīng)關(guān)系,對(duì)于CO而言,上升氣流將下層高混合比空氣引入到上層,增加了上層的CO混合比濃度;對(duì)O3而言,上升可以將下層低混合比空氣引入到上層,稀釋了上層的O3混合比濃度。大部分的H2O只能停留在250 hPa以下,但是在250 hPa以上的對(duì)流層上部,H2O混合比有較明顯的增加;而在原本相對(duì)就比較干燥的對(duì)流層頂附近,H2O混合比減小。IWC的大量增加出現(xiàn)在地面強(qiáng)降水發(fā)生之前,在強(qiáng)降水發(fā)生之后IWC含量迅速減小。

    上述CO和O3的分布發(fā)生變化的原因可以從對(duì)流垂直輸送的角度得到合理的解釋?zhuān)?H2O和IWC等水物質(zhì)在對(duì)流層上部和平流層的分布變化則不盡然。在受到垂直輸送的作用的同時(shí),水分子在不同相態(tài)間的運(yùn)動(dòng)和轉(zhuǎn)化更主要的是由對(duì)流系統(tǒng)上部的微物理過(guò)程決定的。因此,盡管對(duì)流層上部H2O和IWC含量都有一定程度的增加,但是這些增加的不同相態(tài)水的質(zhì)量,帶給對(duì)流層頂以上的水汽收支的直接貢獻(xiàn)并不顯著。相反,強(qiáng)對(duì)流過(guò)程中大量降水粒子的形成并沉降為地面降水,甚至使得對(duì)流層頂附近的 H2O混合比在短時(shí)間內(nèi)出現(xiàn)下降。

    本文對(duì)強(qiáng)對(duì)流天氣過(guò)程和對(duì)流層上層大氣成分分布變化的分析,定性地給出了兩者之間存在的相關(guān)變化關(guān)系,指出了對(duì)流垂直輸送對(duì)對(duì)流層上層大氣成分分布的作用結(jié)果,同時(shí)還發(fā)現(xiàn)了對(duì)流上部冰水含量在地面強(qiáng)降水發(fā)生前的急劇增長(zhǎng),這也許可以作為一個(gè)前兆因子來(lái)預(yù)判強(qiáng)降水的發(fā)生。雖然本文只是對(duì)單一典型個(gè)例進(jìn)行的分析,其普遍性還有待更充實(shí)的證據(jù)進(jìn)一步予以證實(shí),但是其結(jié)果仍然很有意義。在此基礎(chǔ)之上,可以對(duì)更多類(lèi)似個(gè)例進(jìn)行定量的計(jì)算和統(tǒng)計(jì),還可以借助天氣尺度模式進(jìn)行數(shù)值模擬,對(duì)赤道外地區(qū)強(qiáng)對(duì)流過(guò)程對(duì)大氣成分的垂直輸送作用給出更翔實(shí)的例證和更加客觀和量化的評(píng)估。

    致謝感謝美國(guó)國(guó)家航空航天局(NASA)提供的資料數(shù)據(jù)使得本文的工作得以完成。對(duì)審稿人提出的寶貴意見(jiàn)表示感謝!

    References)

    Alcala C M, Dessler A E. 2002. Observations of deep convection in the tropics using the Tropical Rainfall Measuring Mission (TRMM)precipitation radar [J]. J. Geophys. Res., 107 (D24): 4792.

    Dessler A E. 2002. The effect of deep, tropical convection on the tropical tropopause layer [J]. J. Geophys. Res., 107 (D3): 4033.

    Fischer H, de Reus M, Traub M, et al. 2003. Deep convective injection of boundary layer air into the lowermost stratosphere at midlatitudes [J].Atmospheric Chemistry and Physics, 3 (3): 739–745.

    Froidevaux L, Livesey N J, Read W G, et al. 2006. Early validation analyses of atmospheric profiles from EOS MLS on the Aura satellite [J]. IEEE Trans. Geosci. Remote Sens., 44 (5): 1106–1121.

    Haynes J M, Stephens G L. 2007. Tropical oceanic cloudiness and the incidence of precipitation: Early results from CloudSat [J]. Geophys. Res.Lett., 34: L09811.

    Hassim M E E, Lane T P. 2010. A model study on the influence of overshooting convection on TTL water vapour [J]. Atmospheric Chemistry and Physics, 10 (20): 9833–9849.

    Iguchi T, Kozu T, Meneghini R, et al. 2000. Rain-profiling algorithm for the TRMM precipitation radar [J]. J. Appl. Meteor., 39 (12): 2038–2052.

    Kalnay E, Kanamitsu M, Kistler R, et al. 1996. The NCEP/NCAR 40-year reanalysis project [J]. Bull. Amer. Meteor. Soc., 77 (3): 437–472.

    Kistler R, Kalnay E, Collins W, et al. 2001. The NCEP-NCAR 50-year reanalysis: Monthly means CD-ROM and documentation [J]. Bull. Amer.Meteor. Soc., 82 (2): 247–267.

    Kuhn P M, Lojko M S, Petersen E V. 1971. Water vapor: Stratospheric injection by thunderstorms [J]. Science, 174 (4016): 1319–1321.

    Kuhn P M, Stearns L P. 1973. Radiometric observations of atmospheric water vapor injection by thunderstorms [J]. J. Atmos. Sci., 30 (3): 507–509.

    Kummerow C, Barnes W, Kozu T, et al. 1998. The Tropical Rainfall Measuring Mission (TRMM) sensor package [J]. J. Atmos. Oceanic Tech.,15 (3): 809–817.

    Kummerow C, Simpson J, Thiele O, et al. 2000. The status of the Tropical Rainfall Measuring Mission (TRMM) after two years in orbit [J]. J. Appl.Meteor., 39 (12): 1965–1982.

    Livesey N J, Read W G, Lambert A, et al. 2007. Earth Observing System(EOS) Microwave Limb Sounder (MLS) Version 2.2 Level 2 data quality and description document [R]. JPL Technical Report D-33509. Available from http://mls.jpl.nasa.gov/.

    Liu C T, Zipser E J. 2005. Global distribution of convection penetrating the tropical tropopause [J]. J. Geophys. Res., 110 (D23): D23104.

    Meneghini R, Iguchi T, Kozu T, et al. 2000. Use of the surface reference technique for path attenuation estimates from the TRMM precipitation radar [J]. J. Appl. Meteor., 39 (12): 2053–2070.

    Poulida O, Dickerson R R, Heymsfield A. 1996. Stratosphere-troposphere exchange in a midlatitude mesoscale convective complex. 1. Observations[J]. J. Geophys. Res., 101 (D3): 6823–6836.

    Ricaud P, Barret B, Attié J L, et al. 2007. Impact of land convection on troposphere-stratosphere exchange in the tropics [J]. Atmospheric Chemistry and Physics, 7 (21): 5639–5657.

    Salby M, Sassi F, Callaghan P, et al. 2003. Fluctuations of cloud, humidity,and thermal structure near the tropical tropopause [J]. J. Climate, 16 (21):3428–3446.

    Sassi F, Salby M, Read W G. 2001. Relationship between upper tropospheric humidity and deep convection [J]. J. Geophys. Res., 106 (D15):17133–17146.

    Sherwood S C, Dessler A E. 2000. On the control of stratospheric humidity[J]. Geophys. Res. Lett., 27 (16): 2513–2516.

    Showalter A K. 1953. A stability index for thunderstorms forecasting [J].Bull. Amer. Meteor. Soc., 34 (6): 250–252.

    Stephens G L, Kummerow C D. 2007. The remote sensing of clouds and precipitation from space: A review [J]. J. Atmos. Sci., 64 (11): 3742–3765.

    Stephens G L, Vane D G, Boain R J, et al. 2002. The cloudsat mission and the a-train—A new dimension of space-based observations of clouds and precipitation [J]. Bull. Amer. Meteor. Soc., 83 (12): 1771–1790.

    Stephens G L, Vane D G, Tanelli S, et al. 2008. CloudSat mission:Performance and early science after the first year of operation [J]. J.Geophys. Res., 113 (D23): D00A18.

    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 Sens., 44 (5): 1075–1092.

    卞建春. 2009. 上對(duì)流層/下平流層大氣垂直結(jié)構(gòu)研究進(jìn)展 [J]. 地球科學(xué)進(jìn)展, 24 (3): 262–271. Bian Jianchun. 2009. Recent advances in the study of atmospheric vertical structures in upper troposphere and lower stratosphere [J]. Advances in Earth Science (in Chinese), 24 (3):262–271.

    卞建春, 嚴(yán)仁嫦, 陳洪濱. 2011. 亞洲夏季風(fēng)是低層污染物進(jìn)入平流層的重要途徑[J]. 大氣科學(xué), 35 (5): 897–902. Bian Jianchun, Yan Renchang, Chen Hongbin. 2011. Tropospheric pollutant transport to the stratosphere by Asian summer monsoon [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 35 (5): 897–902.

    陳洪濱, 卞建春, 呂達(dá)仁. 2006. 上對(duì)流層—下平流層交換過(guò)程研究的進(jìn)展與展望 [J]. 大氣科學(xué), 30 (5): 813–820. Chen Hongbin, Bian Jianchun, Lü Daren. 2006. Advances and prospects in the study of stratosphere-troposphere exchange [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 2006, 30 (5): 813–820.

    傅云飛, 張愛(ài)民, 劉勇, 等. 2008. 基于星載測(cè)雨雷達(dá)探測(cè)的亞洲對(duì)流和層云降水季尺度特征分析 [J]. 氣象學(xué)報(bào), 66 (5): 730–746. Fu Yunfei,Zhang Aimin, Liu Yong, et al. 2008. Characteristics of seasonal scale convective and stratiform precipitation in Asia based on measurements by TRMM Precipitation Radar [J]. Acta Meteorologica Sinica (in Chinese),66 (5): 730–746.

    劉鵬, 傅云飛. 2010. 利用星載測(cè)雨雷達(dá)探測(cè)結(jié)果對(duì)夏季中國(guó)南方對(duì)流和層云降水氣候特征的分析 [J]. 大氣科學(xué), 34 (4): 802–814. Liu Peng, Fu Yunfei. 2010. Climatic characteristics of summer convective and stratiform precipitation in southern China based on measurements by TRMM precipitation radar [J]. Chinese Journal of Atmospheric Sciences(in Chinese), 34 (4): 802–814.

    馬占山, 劉奇俊, 秦琰琰, 等. 2008. 云探測(cè)衛(wèi)星CloudSat [J]. 氣象, 34(8): 104–111. Ma Zhanshan, Liu Qijun, Qin Yanyan, et al. 2008.Introductions to a new type cloud detecting satellite—CloudSat [J].Meteorological Monthly (in Chinese), 34 (8): 104–111.

    石瑩, 孫成武, 朱時(shí)良, 等. 2009. 2006年7月3日黃淮地區(qū)大暴雨過(guò)程分析[J]. 氣象科學(xué), 29 (1): 110–114. Shi Ying, Sun Chengwu, Zhu Shiliang, et al. 2009. The analysis of a heavy rainfall in Huanghuai area on 3 July, 2006 [J]. Scientia Meteorologica Sinica (in Chinese), 29 (1):110–114.

    Analysis of the Variation of Atmospheric Composition in the Upper Troposphere during a Strong Convection in Eastern China in Summer

    YI Mingjian, FU Yunfei, LIU Peng, HENG Zhiwei, GAO Yue, and HONG Xingyuan
    School of Earth and Space Sciences, University of Science and Technology of China, Hefei230026

    The distribution of atmospheric composition, including carbon monoxide (CO), ozone (O3), water vapor(H2O), and ice water content (IWC), during a strong convective weather occurring in the Jianghuai areas is investigated based on measurements by the MLS (Microwave Limb Sounder) on the Aura satellite. Results show that the increase of CO mixing ratio at 200 hPa is 0.12 ppm, the decrease of O3mixing ratio at 70 hPa is about 0.30 ppm, and the increase of H2O mixing ratio at 250 hPa is 400 ppm. The value of IWC has a sharp growth before the surface heavy rain occurs, and its maximum at 200 hPa is about 0.03 g/m3. The consistency variation between vertical velocity and mixing ratio of CO and O3indicates that the vertical transport is a conceivable mechanism changing the distribution of atmospheric composition in the upper troposphere and the lower stratosphere. Besides, H2O mixing ratio and IWC also increase obviously when the convection is triggered, but the growth is limited below the tropopause. This suggests that not only the vertical transport but also the microphysical processes in the convection system decide the quality and phase of water in the upper troposphere.

    convection, atmospheric composition, microwave limb sounder

    1006-9895(2012)05-0901-11

    P421

    A

    10.3878/j.issn.1006-9895.2012.11124

    2011-07-07, 2012-03-16收修定稿

    國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目2010CB428601,中國(guó)科學(xué)院知識(shí)創(chuàng)新工程重要方向項(xiàng)目KZCX2-YW-Q11-04,國(guó)家自然科學(xué)基金資助項(xiàng)目41075041、40730950、41105031,中央高?;究蒲袠I(yè)務(wù)費(fèi)專(zhuān)項(xiàng)資金資助WK2080000030,中國(guó)博士后科學(xué)基金面上資助 20100470859

    易明建,男,1980年出生,博士后,從事大氣動(dòng)力學(xué)和衛(wèi)星遙感應(yīng)用研究,E-mail: mjyi@ustc.edu.cn

    傅云飛,教授,E-mail: fyf@ustc.edu.cn

    猜你喜歡
    平流層對(duì)流層強(qiáng)對(duì)流
    青藏高原上空平流層水汽的時(shí)空演變特征
    基于副氣囊的平流層浮空器高度控制
    哈爾濱2020年一次局地強(qiáng)對(duì)流天氣分析
    黑龍江氣象(2021年2期)2021-11-05 07:06:54
    郴州地區(qū)對(duì)流層頂氣候概況
    突發(fā)性強(qiáng)對(duì)流天氣快速識(shí)別預(yù)警改進(jìn)方法
    青藏高原東北部地區(qū)一次強(qiáng)對(duì)流天氣特征分析
    實(shí)時(shí)干涉測(cè)量中對(duì)流層延遲與鐘差精修正建模
    載人航天(2016年4期)2016-12-01 06:56:24
    1979~2011年間平流層溫度及平流層水汽的演變趨勢(shì)
    成都地區(qū)2005~2015年對(duì)流層NO2柱濃度趨勢(shì)與時(shí)空分布
    遼寧強(qiáng)對(duì)流天氣物理量閾值探索統(tǒng)計(jì)分析
    午夜爱爱视频在线播放| 久99久视频精品免费| 亚洲熟妇熟女久久| 日韩成人av中文字幕在线观看 | 观看免费一级毛片| 精品无人区乱码1区二区| 亚洲精品一区av在线观看| 可以在线观看毛片的网站| 老司机午夜福利在线观看视频| 少妇的逼好多水| 国产午夜精品久久久久久一区二区三区 | 两个人视频免费观看高清| 免费看a级黄色片| 久久鲁丝午夜福利片| 国产精品一二三区在线看| 国产成人a∨麻豆精品| 舔av片在线| 舔av片在线| 国产片特级美女逼逼视频| 校园春色视频在线观看| 搞女人的毛片| 自拍偷自拍亚洲精品老妇| 在线免费观看的www视频| 国产国拍精品亚洲av在线观看| 国产精品女同一区二区软件| 亚洲四区av| 在现免费观看毛片| 久99久视频精品免费| 国产av一区在线观看免费| 大香蕉久久网| 丰满的人妻完整版| 99久久精品国产国产毛片| 在线播放国产精品三级| 国产午夜福利久久久久久| 久久这里只有精品中国| 18禁裸乳无遮挡免费网站照片| 高清毛片免费看| 亚洲国产精品合色在线| 婷婷色综合大香蕉| 亚洲婷婷狠狠爱综合网| 日本一本二区三区精品| 不卡一级毛片| 一级毛片电影观看 | 久久久久久久久大av| 毛片一级片免费看久久久久| 一本精品99久久精品77| 午夜福利在线观看免费完整高清在 | 丰满的人妻完整版| 好男人在线观看高清免费视频| 男女那种视频在线观看| 国产白丝娇喘喷水9色精品| 黄色一级大片看看| 一边摸一边抽搐一进一小说| 日本欧美国产在线视频| 成人av在线播放网站| 国产麻豆成人av免费视频| 一进一出抽搐动态| 麻豆乱淫一区二区| 免费看av在线观看网站| 欧美极品一区二区三区四区| 校园春色视频在线观看| 青春草视频在线免费观看| 亚洲av二区三区四区| 午夜a级毛片| 性插视频无遮挡在线免费观看| 色av中文字幕| 亚洲中文字幕日韩| 国产精品日韩av在线免费观看| 熟女人妻精品中文字幕| 欧美人与善性xxx| 国产一区二区三区av在线 | 日本一二三区视频观看| 美女被艹到高潮喷水动态| 一进一出抽搐gif免费好疼| 亚洲国产欧美人成| 欧美中文日本在线观看视频| 国产精品1区2区在线观看.| 一级a爱片免费观看的视频| 国产乱人视频| av黄色大香蕉| 99热这里只有精品一区| 天堂av国产一区二区熟女人妻| 91在线观看av| 久久久久国产精品人妻aⅴ院| 最好的美女福利视频网| 最好的美女福利视频网| 99久久中文字幕三级久久日本| 久久精品综合一区二区三区| 嫩草影院精品99| 中国国产av一级| 一个人免费在线观看电影| 少妇猛男粗大的猛烈进出视频 | 欧美zozozo另类| 午夜精品国产一区二区电影 | 少妇人妻精品综合一区二区 | 小说图片视频综合网站| 哪里可以看免费的av片| 欧美丝袜亚洲另类| 老熟妇乱子伦视频在线观看| 久久人人爽人人爽人人片va| 麻豆精品久久久久久蜜桃| 尾随美女入室| 99久久精品热视频| 18禁在线播放成人免费| 卡戴珊不雅视频在线播放| 热99re8久久精品国产| 精品日产1卡2卡| 别揉我奶头 嗯啊视频| 亚洲国产精品久久男人天堂| 成人综合一区亚洲| av天堂在线播放| 一级毛片久久久久久久久女| 欧美极品一区二区三区四区| 我的老师免费观看完整版| 床上黄色一级片| 精品无人区乱码1区二区| 在线国产一区二区在线| 成人国产麻豆网| 天天躁日日操中文字幕| 女的被弄到高潮叫床怎么办| 国产亚洲精品综合一区在线观看| 国产精品亚洲美女久久久| 亚洲欧美成人精品一区二区| 国产精品乱码一区二三区的特点| 免费无遮挡裸体视频| 国产爱豆传媒在线观看| 午夜a级毛片| 日本a在线网址| 久久精品国产亚洲av天美| 日韩,欧美,国产一区二区三区 | 天美传媒精品一区二区| 国产精品久久久久久av不卡| 国内精品宾馆在线| 热99re8久久精品国产| 免费黄网站久久成人精品| 精品日产1卡2卡| 成人精品一区二区免费| 精品人妻偷拍中文字幕| 无遮挡黄片免费观看| 国产精品电影一区二区三区| 国产av麻豆久久久久久久| 国产免费男女视频| 精品一区二区三区视频在线观看免费| 五月伊人婷婷丁香| 99热只有精品国产| 有码 亚洲区| 日本免费a在线| 午夜亚洲福利在线播放| 俺也久久电影网| 桃色一区二区三区在线观看| 日韩欧美国产在线观看| 久久久久国产精品人妻aⅴ院| 身体一侧抽搐| 亚洲真实伦在线观看| 国产蜜桃级精品一区二区三区| 丰满乱子伦码专区| 日日摸夜夜添夜夜爱| 久久人人爽人人片av| 国产精品乱码一区二三区的特点| 此物有八面人人有两片| 天堂√8在线中文| a级毛片a级免费在线| 日本五十路高清| avwww免费| 欧美一区二区国产精品久久精品| 免费人成视频x8x8入口观看| 少妇人妻精品综合一区二区 | 久久久久久久久大av| 女生性感内裤真人,穿戴方法视频| 一个人看的www免费观看视频| 国产精品综合久久久久久久免费| 欧美bdsm另类| 欧美性猛交黑人性爽| 简卡轻食公司| 午夜精品一区二区三区免费看| 国产午夜福利久久久久久| 五月玫瑰六月丁香| 给我免费播放毛片高清在线观看| 亚洲第一电影网av| 免费av不卡在线播放| 蜜臀久久99精品久久宅男| 国产精品日韩av在线免费观看| 99在线视频只有这里精品首页| 尤物成人国产欧美一区二区三区| 黄色配什么色好看| 亚洲精品一区av在线观看| 亚洲无线在线观看| 国产成人aa在线观看| 国产大屁股一区二区在线视频| 免费看光身美女| 亚洲美女黄片视频| 欧美一区二区亚洲| 97超级碰碰碰精品色视频在线观看| 日日摸夜夜添夜夜添小说| 身体一侧抽搐| 国产精品一区二区三区四区久久| 男插女下体视频免费在线播放| 色尼玛亚洲综合影院| 国产免费男女视频| 欧美日韩在线观看h| 超碰av人人做人人爽久久| 99热全是精品| 国产伦精品一区二区三区四那| 高清午夜精品一区二区三区 | 一级黄色大片毛片| 一级a爱片免费观看的视频| 嫩草影院入口| 日日摸夜夜添夜夜添av毛片| 国产精品一区二区三区四区免费观看 | 国产精品美女特级片免费视频播放器| 亚洲专区国产一区二区| 亚洲av一区综合| 久久人妻av系列| 欧美日本亚洲视频在线播放| 成人av在线播放网站| 国产精品国产高清国产av| 亚洲熟妇中文字幕五十中出| 亚洲久久久久久中文字幕| 91久久精品电影网| 日韩欧美精品免费久久| 国产伦在线观看视频一区| 我的女老师完整版在线观看| 老司机福利观看| 老熟妇仑乱视频hdxx| 亚洲aⅴ乱码一区二区在线播放| 天堂网av新在线| 国产一区二区三区av在线 | 精品久久久噜噜| 婷婷色综合大香蕉| 夜夜夜夜夜久久久久| 亚洲人成网站在线播放欧美日韩| 日韩av在线大香蕉| 日日摸夜夜添夜夜添av毛片| 国产精品国产三级国产av玫瑰| 美女 人体艺术 gogo| 在线观看免费视频日本深夜| 69人妻影院| 成人综合一区亚洲| 亚洲美女视频黄频| 99久国产av精品国产电影| 久久精品国产鲁丝片午夜精品| 卡戴珊不雅视频在线播放| 成人亚洲精品av一区二区| 亚洲成av人片在线播放无| av在线蜜桃| 国产精品日韩av在线免费观看| 国产精品嫩草影院av在线观看| 18+在线观看网站| av在线亚洲专区| 在线a可以看的网站| 精品一区二区三区视频在线| 亚洲在线自拍视频| 寂寞人妻少妇视频99o| 亚洲无线在线观看| 亚洲欧美精品综合久久99| 白带黄色成豆腐渣| 色噜噜av男人的天堂激情| 性色avwww在线观看| 日本在线视频免费播放| 精品不卡国产一区二区三区| 欧美精品国产亚洲| 国产美女午夜福利| 老司机午夜福利在线观看视频| 成人欧美大片| 熟女人妻精品中文字幕| 99热这里只有精品一区| 天天躁夜夜躁狠狠久久av| 国产av在哪里看| 激情 狠狠 欧美| 非洲黑人性xxxx精品又粗又长| 亚洲欧美日韩无卡精品| 国产真实乱freesex| 免费黄网站久久成人精品| 日韩国内少妇激情av| 亚洲成人精品中文字幕电影| 日本爱情动作片www.在线观看 | 亚洲精品乱码久久久v下载方式| 亚洲精品一卡2卡三卡4卡5卡| 精品国产三级普通话版| 日韩成人av中文字幕在线观看 | 日韩国内少妇激情av| 欧美高清性xxxxhd video| 中国美白少妇内射xxxbb| 亚洲欧美日韩无卡精品| 丰满乱子伦码专区| 六月丁香七月| 联通29元200g的流量卡| 国产黄片美女视频| 看非洲黑人一级黄片| 露出奶头的视频| 国产91av在线免费观看| 欧美激情国产日韩精品一区| 波多野结衣高清无吗| 观看美女的网站| 最新在线观看一区二区三区| 亚洲最大成人中文| 在线免费观看的www视频| 久99久视频精品免费| 精品少妇黑人巨大在线播放 | 晚上一个人看的免费电影| 国语自产精品视频在线第100页| 波多野结衣高清作品| 久久久久性生活片| 国产精品一区二区三区四区免费观看 | 九色成人免费人妻av| 色视频www国产| 亚洲色图av天堂| 亚洲熟妇熟女久久| 亚洲一区高清亚洲精品| 国产老妇女一区| 有码 亚洲区| av女优亚洲男人天堂| 日韩成人伦理影院| 午夜a级毛片| 能在线免费观看的黄片| 国产av麻豆久久久久久久| 国产乱人视频| 亚洲18禁久久av| 99在线视频只有这里精品首页| 美女被艹到高潮喷水动态| 日日啪夜夜撸| 欧美性感艳星| 中出人妻视频一区二区| 大香蕉久久网| 欧美日韩精品成人综合77777| 99国产精品一区二区蜜桃av| 99久国产av精品| h日本视频在线播放| 日韩强制内射视频| av专区在线播放| 国产精品人妻久久久影院| 在线播放国产精品三级| 久久久成人免费电影| 亚洲欧美精品综合久久99| 久久欧美精品欧美久久欧美| 麻豆乱淫一区二区| 永久网站在线| 欧美色视频一区免费| 久久精品91蜜桃| av卡一久久| 国产探花极品一区二区| 狠狠狠狠99中文字幕| 日日摸夜夜添夜夜添av毛片| 99九九线精品视频在线观看视频| 亚洲性久久影院| 九九热线精品视视频播放| av专区在线播放| 亚洲熟妇熟女久久| 精品久久久久久久人妻蜜臀av| 国产av不卡久久| 亚洲欧美日韩高清在线视频| 亚洲内射少妇av| 欧美一级a爱片免费观看看| 色吧在线观看| 亚洲av成人精品一区久久| 99热只有精品国产| 男女那种视频在线观看| 国产一区亚洲一区在线观看| 淫妇啪啪啪对白视频| 亚洲av免费在线观看| 又爽又黄a免费视频| 日韩欧美一区二区三区在线观看| 亚洲精品乱码久久久v下载方式| avwww免费| 精品午夜福利视频在线观看一区| 精品国产三级普通话版| 精品人妻一区二区三区麻豆 | 精品无人区乱码1区二区| 中文字幕av成人在线电影| 久99久视频精品免费| 22中文网久久字幕| 麻豆乱淫一区二区| 国产 一区精品| 啦啦啦啦在线视频资源| 3wmmmm亚洲av在线观看| 直男gayav资源| 国产成人91sexporn| 直男gayav资源| 丰满的人妻完整版| 老女人水多毛片| 日韩欧美三级三区| 国产色婷婷99| 麻豆国产97在线/欧美| 成人综合一区亚洲| 欧美性猛交黑人性爽| 国产亚洲精品综合一区在线观看| 听说在线观看完整版免费高清| 人妻丰满熟妇av一区二区三区| 少妇被粗大猛烈的视频| 最近在线观看免费完整版| 色尼玛亚洲综合影院| 观看免费一级毛片| 婷婷亚洲欧美| 99热只有精品国产| 一级黄片播放器| 三级经典国产精品| 97在线视频观看| 国产日本99.免费观看| 国产麻豆成人av免费视频| 国产精品电影一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 日韩制服骚丝袜av| 亚洲国产精品国产精品| 国产伦一二天堂av在线观看| 别揉我奶头 嗯啊视频| av.在线天堂| 色综合亚洲欧美另类图片| videossex国产| 午夜日韩欧美国产| 丰满的人妻完整版| 久久九九热精品免费| 在线国产一区二区在线| 国产精品久久久久久精品电影| 亚洲国产精品久久男人天堂| 国产精品免费一区二区三区在线| 亚洲丝袜综合中文字幕| 中文字幕熟女人妻在线| 特级一级黄色大片| 免费搜索国产男女视频| 免费观看在线日韩| 少妇人妻精品综合一区二区 | 在线天堂最新版资源| 婷婷六月久久综合丁香| 国内精品宾馆在线| 国产午夜精品论理片| 国产 一区精品| 俄罗斯特黄特色一大片| 亚洲人成网站在线播| 一个人看的www免费观看视频| 搡女人真爽免费视频火全软件 | 午夜视频国产福利| 成人漫画全彩无遮挡| 少妇熟女aⅴ在线视频| 69av精品久久久久久| 三级经典国产精品| 久久久色成人| 亚洲精品在线观看二区| 欧美成人免费av一区二区三区| 少妇人妻一区二区三区视频| 亚洲国产精品久久男人天堂| 十八禁网站免费在线| 亚洲中文日韩欧美视频| 一区二区三区四区激情视频 | 日日撸夜夜添| 亚洲av成人精品一区久久| 国产精品三级大全| 一个人看视频在线观看www免费| 性插视频无遮挡在线免费观看| 女同久久另类99精品国产91| 别揉我奶头~嗯~啊~动态视频| 久久欧美精品欧美久久欧美| 午夜精品国产一区二区电影 | av在线观看视频网站免费| 亚州av有码| 免费av毛片视频| 久久人人爽人人片av| 日韩精品中文字幕看吧| 久久久精品欧美日韩精品| 国产成人a区在线观看| 欧美日本亚洲视频在线播放| 最近在线观看免费完整版| 国内精品美女久久久久久| 1024手机看黄色片| 成年版毛片免费区| 俄罗斯特黄特色一大片| 国产精品日韩av在线免费观看| 最近2019中文字幕mv第一页| 国产精品久久久久久久电影| av卡一久久| 午夜福利在线观看免费完整高清在 | 亚洲av熟女| 国产精品1区2区在线观看.| 俺也久久电影网| 精品人妻偷拍中文字幕| 一级av片app| 婷婷六月久久综合丁香| 联通29元200g的流量卡| 一个人免费在线观看电影| 尾随美女入室| 麻豆国产97在线/欧美| 国产淫片久久久久久久久| 午夜精品一区二区三区免费看| 亚洲中文字幕一区二区三区有码在线看| av卡一久久| 老熟妇乱子伦视频在线观看| 国产精品久久久久久av不卡| 三级经典国产精品| 美女内射精品一级片tv| 亚洲欧美日韩东京热| 欧美日韩乱码在线| 特级一级黄色大片| 在线观看美女被高潮喷水网站| 国产精品亚洲美女久久久| 18禁在线播放成人免费| 91在线观看av| 深爱激情五月婷婷| 亚洲最大成人中文| 国产精品久久视频播放| 你懂的网址亚洲精品在线观看 | 男人的好看免费观看在线视频| 伦精品一区二区三区| 成人漫画全彩无遮挡| 亚洲精华国产精华液的使用体验 | 亚洲精品粉嫩美女一区| 色尼玛亚洲综合影院| 男女做爰动态图高潮gif福利片| 成熟少妇高潮喷水视频| 少妇人妻一区二区三区视频| 亚洲av电影不卡..在线观看| 免费av毛片视频| 国产色婷婷99| 99热只有精品国产| 午夜精品一区二区三区免费看| 亚洲18禁久久av| 精品久久久久久久久久免费视频| 日日撸夜夜添| 亚洲在线自拍视频| 十八禁网站免费在线| 欧美一级a爱片免费观看看| 亚洲精品色激情综合| 男女那种视频在线观看| 免费电影在线观看免费观看| 亚洲人成网站在线播| 国产不卡一卡二| 麻豆精品久久久久久蜜桃| av中文乱码字幕在线| 国产精品美女特级片免费视频播放器| 久久热精品热| 欧美色欧美亚洲另类二区| 国产成人91sexporn| 午夜日韩欧美国产| 97热精品久久久久久| 欧美又色又爽又黄视频| 午夜久久久久精精品| 中国国产av一级| 国产亚洲精品久久久久久毛片| 最近的中文字幕免费完整| 草草在线视频免费看| 少妇的逼水好多| 一区二区三区四区激情视频 | 十八禁网站免费在线| 天天一区二区日本电影三级| 午夜精品在线福利| 国产一区二区在线观看日韩| 日韩一区二区视频免费看| 成人午夜高清在线视频| 一进一出抽搐gif免费好疼| 如何舔出高潮| 久久人人爽人人爽人人片va| 日韩,欧美,国产一区二区三区 | videossex国产| 日本五十路高清| 精品久久久久久久久久久久久| 亚洲无线在线观看| 久久久久国产精品人妻aⅴ院| 简卡轻食公司| 男女做爰动态图高潮gif福利片| 看黄色毛片网站| 成年女人看的毛片在线观看| 韩国av在线不卡| 99久久中文字幕三级久久日本| 成人永久免费在线观看视频| 亚洲av免费高清在线观看| 日韩中字成人| 日日摸夜夜添夜夜添小说| 亚洲aⅴ乱码一区二区在线播放| 成人av一区二区三区在线看| 一卡2卡三卡四卡精品乱码亚洲| 欧美一区二区国产精品久久精品| 在线观看免费视频日本深夜| 桃色一区二区三区在线观看| 99久久中文字幕三级久久日本| 特级一级黄色大片| 国产真实乱freesex| 国产伦精品一区二区三区视频9| 日日摸夜夜添夜夜添小说| 精品久久久久久久末码| 日日啪夜夜撸| 国产 一区 欧美 日韩| 久久综合国产亚洲精品| 国产极品精品免费视频能看的| 极品教师在线视频| 亚洲精品粉嫩美女一区| 在线观看66精品国产| 啦啦啦啦在线视频资源| 变态另类丝袜制服| 三级男女做爰猛烈吃奶摸视频| aaaaa片日本免费| 亚洲18禁久久av| 免费无遮挡裸体视频| 麻豆精品久久久久久蜜桃| 五月玫瑰六月丁香| 大香蕉久久网| 亚洲av免费高清在线观看| 91午夜精品亚洲一区二区三区| 啦啦啦啦在线视频资源| 国产大屁股一区二区在线视频| 一a级毛片在线观看| 国产乱人偷精品视频| 亚洲av美国av| 成人一区二区视频在线观看| 91久久精品国产一区二区三区| 国产免费男女视频| 一个人观看的视频www高清免费观看| 亚洲人成网站在线播放欧美日韩| 国产精品嫩草影院av在线观看| 欧美+亚洲+日韩+国产| 国产单亲对白刺激| 男女做爰动态图高潮gif福利片| 精品久久久久久久久久免费视频| 久久久久久国产a免费观看| 色综合色国产| 久久99热6这里只有精品| 天堂动漫精品| 久久久久国产精品人妻aⅴ院|