• <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ì)分析
    成人一区二区视频在线观看| 亚洲中文字幕一区二区三区有码在线看 | 1024香蕉在线观看| 亚洲精华国产精华精| 两人在一起打扑克的视频| 成人欧美大片| 12—13女人毛片做爰片一| 国产精品av视频在线免费观看| 国产一区二区三区在线臀色熟女| av片东京热男人的天堂| 国产精品亚洲美女久久久| 999久久久国产精品视频| 丰满人妻一区二区三区视频av | 中国美女看黄片| 性欧美人与动物交配| 中文资源天堂在线| 久久精品人妻少妇| 九色成人免费人妻av| 男插女下体视频免费在线播放| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av中文字字幕乱码综合| 欧美+亚洲+日韩+国产| 19禁男女啪啪无遮挡网站| 日本一二三区视频观看| 淫妇啪啪啪对白视频| 变态另类成人亚洲欧美熟女| 美女午夜性视频免费| 国产精品女同一区二区软件 | 女人高潮潮喷娇喘18禁视频| 久久久国产精品麻豆| 久久久国产成人免费| 久久久久久久久免费视频了| 亚洲精品在线美女| 99久久精品热视频| 国产成人精品无人区| 夜夜夜夜夜久久久久| 亚洲 欧美一区二区三区| 小说图片视频综合网站| 亚洲成人中文字幕在线播放| 午夜免费观看网址| 亚洲精品色激情综合| 久久这里只有精品中国| 亚洲成人中文字幕在线播放| or卡值多少钱| 免费在线观看亚洲国产| 久久久久久久久中文| 免费在线观看影片大全网站| 亚洲精品美女久久久久99蜜臀| 成人无遮挡网站| 成年免费大片在线观看| 熟女少妇亚洲综合色aaa.| 成熟少妇高潮喷水视频| 欧美丝袜亚洲另类 | 国产毛片a区久久久久| 精品人妻1区二区| 国产一区二区在线av高清观看| 两性夫妻黄色片| svipshipincom国产片| 欧美一级a爱片免费观看看| 亚洲人成网站高清观看| 欧美中文综合在线视频| 特大巨黑吊av在线直播| 欧美乱码精品一区二区三区| 欧美在线黄色| 亚洲第一欧美日韩一区二区三区| 看片在线看免费视频| 特级一级黄色大片| netflix在线观看网站| 久久午夜综合久久蜜桃| 免费观看精品视频网站| 好男人电影高清在线观看| 国产午夜精品久久久久久| 亚洲欧美日韩高清专用| 久久久久久久久中文| 男人舔奶头视频| 国产精品九九99| 亚洲色图 男人天堂 中文字幕| 久久草成人影院| av女优亚洲男人天堂 | 亚洲国产日韩欧美精品在线观看 | 免费高清视频大片| 热99re8久久精品国产| 成人欧美大片| 男人舔女人的私密视频| 听说在线观看完整版免费高清| 美女 人体艺术 gogo| 欧美乱码精品一区二区三区| 99re在线观看精品视频| 午夜福利在线观看吧| 久久久久精品国产欧美久久久| 国产精品一区二区精品视频观看| 日韩三级视频一区二区三区| 免费在线观看成人毛片| www日本在线高清视频| 悠悠久久av| 国产高清视频在线观看网站| 精品无人区乱码1区二区| 91在线观看av| 成人国产一区最新在线观看| 在线观看美女被高潮喷水网站 | 90打野战视频偷拍视频| 国产精品亚洲av一区麻豆| 不卡av一区二区三区| 亚洲成人精品中文字幕电影| 又大又爽又粗| 两个人看的免费小视频| 亚洲 国产 在线| 九九热线精品视视频播放| 精品无人区乱码1区二区| 久久人妻av系列| 成人特级黄色片久久久久久久| 亚洲国产欧洲综合997久久,| av女优亚洲男人天堂 | 99久久无色码亚洲精品果冻| 亚洲国产中文字幕在线视频| 亚洲精品在线美女| 免费大片18禁| 成年女人毛片免费观看观看9| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲片人在线观看| 国产亚洲av嫩草精品影院| 看免费av毛片| 久久久久久九九精品二区国产| 香蕉av资源在线| 久久性视频一级片| 国产熟女xx| 久久国产乱子伦精品免费另类| 岛国在线免费视频观看| e午夜精品久久久久久久| 亚洲av成人一区二区三| 亚洲国产精品sss在线观看| 两个人视频免费观看高清| 熟女少妇亚洲综合色aaa.| 午夜两性在线视频| 嫩草影院精品99| 欧美黄色淫秽网站| av天堂中文字幕网| 亚洲精品国产精品久久久不卡| 99在线视频只有这里精品首页| 99国产精品99久久久久| 午夜a级毛片| 中亚洲国语对白在线视频| 天堂网av新在线| 亚洲专区中文字幕在线| 国产高清有码在线观看视频| 国产成人av激情在线播放| 男人舔奶头视频| 少妇裸体淫交视频免费看高清| 欧美成人免费av一区二区三区| 成人一区二区视频在线观看| 亚洲最大成人中文| 曰老女人黄片| 久久久久久人人人人人| 最新中文字幕久久久久 | 国产精品99久久久久久久久| 精品国产美女av久久久久小说| 毛片女人毛片| 午夜日韩欧美国产| 国产免费av片在线观看野外av| 午夜精品在线福利| 午夜精品一区二区三区免费看| 国产精品亚洲美女久久久| 特级一级黄色大片| 国产精品乱码一区二三区的特点| 一区二区三区激情视频| 欧美黑人欧美精品刺激| 桃色一区二区三区在线观看| 精品不卡国产一区二区三区| 欧美性猛交╳xxx乱大交人| 午夜福利成人在线免费观看| 老汉色∧v一级毛片| 成人特级av手机在线观看| 欧美中文日本在线观看视频| 亚洲色图av天堂| 五月玫瑰六月丁香| 最新在线观看一区二区三区| 亚洲国产高清在线一区二区三| 久久香蕉国产精品| 免费看日本二区| 亚洲av五月六月丁香网| 国产一区二区在线av高清观看| 狂野欧美激情性xxxx| 三级男女做爰猛烈吃奶摸视频| 长腿黑丝高跟| 女人被狂操c到高潮| 亚洲专区国产一区二区| 99国产综合亚洲精品| 成人三级做爰电影| 男人舔女人的私密视频| 久久久久九九精品影院| 老司机深夜福利视频在线观看| 亚洲美女视频黄频| 久久久久性生活片| 欧美三级亚洲精品| 国产高清有码在线观看视频| 国产激情偷乱视频一区二区| 舔av片在线| 成年人黄色毛片网站| 精品一区二区三区四区五区乱码| 国内久久婷婷六月综合欲色啪| 最近最新中文字幕大全电影3| 久久久久国产精品人妻aⅴ院| 国产精品一及| 成人性生交大片免费视频hd| 国产精品99久久99久久久不卡| 国产高清视频在线观看网站| 亚洲avbb在线观看| 亚洲精品一区av在线观看| 亚洲国产精品成人综合色| svipshipincom国产片| av天堂中文字幕网| 嫩草影视91久久| 黄色女人牲交| 亚洲精品粉嫩美女一区| 免费在线观看亚洲国产| 亚洲专区国产一区二区| 欧美成人性av电影在线观看| 亚洲av成人精品一区久久| av天堂在线播放| 亚洲国产欧美一区二区综合| 韩国av一区二区三区四区| 精品国产三级普通话版| 三级毛片av免费| 亚洲欧美日韩东京热| 国产野战对白在线观看| 国产91精品成人一区二区三区| 国产精品爽爽va在线观看网站| 色播亚洲综合网| 日本一二三区视频观看| 欧美国产日韩亚洲一区| 手机成人av网站| 人人妻,人人澡人人爽秒播| 国产精品一区二区三区四区久久| 中文字幕熟女人妻在线| 十八禁人妻一区二区| 免费在线观看成人毛片| 桃色一区二区三区在线观看| 在线播放国产精品三级| 香蕉久久夜色| 亚洲一区二区三区色噜噜| 五月玫瑰六月丁香| 午夜免费激情av| 男女那种视频在线观看| 身体一侧抽搐| 一进一出好大好爽视频| 亚洲中文av在线| 亚洲中文字幕一区二区三区有码在线看 | 亚洲中文字幕日韩| 少妇人妻一区二区三区视频| 19禁男女啪啪无遮挡网站| 国产精华一区二区三区| 久久久久国产一级毛片高清牌| 香蕉久久夜色| 午夜影院日韩av| 91av网站免费观看| 亚洲av成人精品一区久久| 成人三级黄色视频| 久久午夜亚洲精品久久| 又紧又爽又黄一区二区| 性色av乱码一区二区三区2| 黑人欧美特级aaaaaa片| 国内揄拍国产精品人妻在线| 最新在线观看一区二区三区| 国产99白浆流出| 亚洲七黄色美女视频| 人妻丰满熟妇av一区二区三区| 亚洲国产看品久久| 成人鲁丝片一二三区免费| av在线蜜桃| 一个人观看的视频www高清免费观看 | 欧美一区二区精品小视频在线| 99久久无色码亚洲精品果冻| 亚洲中文日韩欧美视频| 99在线人妻在线中文字幕| 欧美中文综合在线视频| 97人妻精品一区二区三区麻豆| 久99久视频精品免费| av天堂中文字幕网| 欧美在线一区亚洲| 两性午夜刺激爽爽歪歪视频在线观看| 久久久精品欧美日韩精品| 两个人看的免费小视频| 国产蜜桃级精品一区二区三区| x7x7x7水蜜桃| 亚洲一区二区三区色噜噜| 亚洲国产精品成人综合色| 黑人操中国人逼视频| 男女床上黄色一级片免费看| 久久午夜综合久久蜜桃| 黄色女人牲交| 欧美成人免费av一区二区三区| 国产三级中文精品| 亚洲人成电影免费在线| 国产一区二区三区视频了| 亚洲黑人精品在线| 国产精品99久久久久久久久| 亚洲狠狠婷婷综合久久图片| 一级毛片高清免费大全| 亚洲精品粉嫩美女一区| 可以在线观看毛片的网站| 久久婷婷人人爽人人干人人爱| 首页视频小说图片口味搜索| 午夜免费成人在线视频| 久久久色成人| 国产精品野战在线观看| 国产极品精品免费视频能看的| 757午夜福利合集在线观看| 午夜免费成人在线视频| 亚洲av日韩精品久久久久久密| 国产麻豆成人av免费视频| 精品国内亚洲2022精品成人| 国产蜜桃级精品一区二区三区| 欧美丝袜亚洲另类 | 国产真实乱freesex| 黄色日韩在线| 国产精品精品国产色婷婷| 成人一区二区视频在线观看| 身体一侧抽搐| 噜噜噜噜噜久久久久久91| 国产麻豆成人av免费视频| 免费一级毛片在线播放高清视频| 动漫黄色视频在线观看| 男女之事视频高清在线观看| 亚洲成av人片在线播放无| АⅤ资源中文在线天堂| 日韩欧美国产在线观看| 国产又色又爽无遮挡免费看| 国产成人欧美在线观看| 搡老岳熟女国产| 免费在线观看亚洲国产| 亚洲真实伦在线观看| 亚洲人成电影免费在线| 麻豆成人av在线观看| 深夜精品福利| 亚洲欧美日韩卡通动漫| 亚洲人成电影免费在线| 欧美乱色亚洲激情| 国产亚洲精品久久久com| 亚洲va日本ⅴa欧美va伊人久久| 女生性感内裤真人,穿戴方法视频| 国产在线精品亚洲第一网站| 亚洲熟妇中文字幕五十中出| 久久久精品大字幕| 国内揄拍国产精品人妻在线| 亚洲av电影不卡..在线观看| 亚洲国产看品久久| 99久久综合精品五月天人人| 精品久久久久久,| 高清毛片免费观看视频网站| av天堂中文字幕网| 麻豆av在线久日| 手机成人av网站| 亚洲七黄色美女视频| 亚洲自偷自拍图片 自拍| 亚洲性夜色夜夜综合| 99精品在免费线老司机午夜| 亚洲av日韩精品久久久久久密| 免费观看精品视频网站| 国产野战对白在线观看| 国产成人影院久久av| 小蜜桃在线观看免费完整版高清| 亚洲专区中文字幕在线| 午夜福利免费观看在线| 午夜激情欧美在线| 欧美精品啪啪一区二区三区| 美女午夜性视频免费| 免费人成视频x8x8入口观看| 真实男女啪啪啪动态图| 精品不卡国产一区二区三区| 亚洲熟妇中文字幕五十中出| 精品99又大又爽又粗少妇毛片 | 一进一出好大好爽视频| 免费看十八禁软件| 午夜精品一区二区三区免费看| 1000部很黄的大片| 黄色视频,在线免费观看| 国产v大片淫在线免费观看| 亚洲色图av天堂| 久久中文看片网| 99在线视频只有这里精品首页| 欧美不卡视频在线免费观看| АⅤ资源中文在线天堂| 成年女人毛片免费观看观看9| ponron亚洲| 麻豆成人午夜福利视频| 在线观看舔阴道视频| 美女免费视频网站| 午夜福利在线观看吧| 黑人欧美特级aaaaaa片| 日本黄大片高清| 偷拍熟女少妇极品色| 久久热在线av| 欧美一区二区精品小视频在线| 51午夜福利影视在线观看| 国产精品一区二区免费欧美| 亚洲真实伦在线观看| 亚洲精华国产精华精| 特级一级黄色大片| 欧美性猛交黑人性爽| 啦啦啦韩国在线观看视频| 天堂√8在线中文| 久久久精品大字幕| 一二三四社区在线视频社区8| 国产高清三级在线| 天堂动漫精品| 精品国内亚洲2022精品成人| 视频区欧美日本亚洲| 亚洲中文日韩欧美视频| 黄频高清免费视频| 国产精品98久久久久久宅男小说| 久久精品91无色码中文字幕| 午夜福利在线观看免费完整高清在 | 久久久久久九九精品二区国产| 精品熟女少妇八av免费久了| 婷婷精品国产亚洲av在线| 夜夜夜夜夜久久久久| 亚洲成人中文字幕在线播放| 欧美国产日韩亚洲一区| 久久久久久久久久黄片| 一区二区三区国产精品乱码| 免费大片18禁| 天天一区二区日本电影三级| a级毛片在线看网站| 亚洲国产高清在线一区二区三| 日韩高清综合在线| 色老头精品视频在线观看| 色播亚洲综合网| 久久久久国产精品人妻aⅴ院| 狠狠狠狠99中文字幕| 日韩精品中文字幕看吧| www.www免费av| 久久香蕉精品热| 成年人黄色毛片网站| 999久久久国产精品视频| 精品一区二区三区四区五区乱码| 日韩欧美精品v在线| 国产精品久久电影中文字幕| 一级毛片高清免费大全| 亚洲精品粉嫩美女一区| 久久国产精品人妻蜜桃| 一a级毛片在线观看| 午夜成年电影在线免费观看| 亚洲成人精品中文字幕电影| 99久久国产精品久久久| 一区二区三区高清视频在线| e午夜精品久久久久久久| 国产成人系列免费观看| 亚洲天堂国产精品一区在线| 变态另类丝袜制服| 欧美乱码精品一区二区三区| 亚洲国产欧美一区二区综合| 精品国产亚洲在线| 一级作爱视频免费观看| 精品久久久久久,| 成人亚洲精品av一区二区| 99久久精品热视频| 日日夜夜操网爽| 最近视频中文字幕2019在线8| 久久久久久久久中文| 精品免费久久久久久久清纯| 三级男女做爰猛烈吃奶摸视频| 给我免费播放毛片高清在线观看| 好男人电影高清在线观看| 亚洲一区高清亚洲精品| 在线观看一区二区三区| 国产高清videossex| 亚洲av成人av| 波多野结衣巨乳人妻| 在线观看美女被高潮喷水网站 | 精品午夜福利视频在线观看一区| 久久这里只有精品中国| 亚洲欧美日韩高清在线视频| 中文字幕熟女人妻在线| 成人国产综合亚洲| 青草久久国产| 日本与韩国留学比较| 精品99又大又爽又粗少妇毛片 | 中文资源天堂在线| 啪啪无遮挡十八禁网站| 精品一区二区三区视频在线观看免费| 在线观看免费午夜福利视频| 天天躁狠狠躁夜夜躁狠狠躁| 国产高清videossex| 国产精品九九99| 欧美乱色亚洲激情| 国模一区二区三区四区视频 | 狂野欧美激情性xxxx| 国产蜜桃级精品一区二区三区| 啦啦啦韩国在线观看视频| 俄罗斯特黄特色一大片| 一个人看的www免费观看视频| 久久九九热精品免费| 97碰自拍视频| 精品99又大又爽又粗少妇毛片 | 国产伦在线观看视频一区| 99久久无色码亚洲精品果冻| 老汉色∧v一级毛片| av片东京热男人的天堂| 婷婷六月久久综合丁香| 窝窝影院91人妻| 黄色视频,在线免费观看| 国产精品98久久久久久宅男小说| 国产欧美日韩一区二区精品| 国内精品美女久久久久久| 99久久国产精品久久久| 色噜噜av男人的天堂激情| 国产综合懂色| 黄色视频,在线免费观看| 国产黄片美女视频| 在线观看日韩欧美| 一a级毛片在线观看| 国产成人av激情在线播放| 黄色视频,在线免费观看| 在线观看66精品国产| 日韩有码中文字幕| 国产一区二区在线观看日韩 | 欧美乱码精品一区二区三区| 中亚洲国语对白在线视频| 久久久久久人人人人人| 国产97色在线日韩免费| 日韩欧美三级三区| 一级毛片高清免费大全| 免费人成视频x8x8入口观看| 日本免费一区二区三区高清不卡| 色综合欧美亚洲国产小说| 一二三四社区在线视频社区8| 欧美日韩国产亚洲二区| 最近最新免费中文字幕在线| 亚洲专区中文字幕在线| 欧美另类亚洲清纯唯美| 亚洲中文av在线| 757午夜福利合集在线观看| 毛片女人毛片| 美女免费视频网站| 天堂网av新在线| 精品免费久久久久久久清纯| 午夜福利在线观看免费完整高清在 | 国产精品av久久久久免费| 变态另类成人亚洲欧美熟女| 狂野欧美白嫩少妇大欣赏| 一a级毛片在线观看| 啦啦啦观看免费观看视频高清| 极品教师在线免费播放| 亚洲国产高清在线一区二区三| 在线看三级毛片| 国产主播在线观看一区二区| 久久久久久久久中文| 日本 欧美在线| 亚洲专区中文字幕在线| 日韩欧美在线二视频| 亚洲精品美女久久av网站| 免费看光身美女| 久久久色成人| 巨乳人妻的诱惑在线观看| 日日夜夜操网爽| 国产精品自产拍在线观看55亚洲| 亚洲一区二区三区不卡视频| 日本免费a在线| 欧美日韩瑟瑟在线播放| 欧美黑人巨大hd| 麻豆成人av在线观看| 香蕉国产在线看| 老汉色∧v一级毛片| 国产淫片久久久久久久久 | 亚洲国产欧美一区二区综合| 亚洲欧美精品综合一区二区三区| 欧美黑人巨大hd| 婷婷亚洲欧美| 91久久精品国产一区二区成人 | 日韩人妻高清精品专区| 欧美日韩亚洲国产一区二区在线观看| 1024香蕉在线观看| 国内精品久久久久精免费| 变态另类丝袜制服| 亚洲欧美日韩高清在线视频| 久久中文看片网| 精品久久久久久久久久免费视频| 99久久精品热视频| 最近最新中文字幕大全电影3| 99热只有精品国产| 欧美最黄视频在线播放免费| 亚洲精品国产精品久久久不卡| 亚洲av片天天在线观看| 亚洲性夜色夜夜综合| 亚洲欧美日韩东京热| 国模一区二区三区四区视频 | 亚洲av熟女| 夜夜看夜夜爽夜夜摸| 久久久国产成人精品二区| 中文字幕人妻丝袜一区二区| 又黄又粗又硬又大视频| 给我免费播放毛片高清在线观看| 国产不卡一卡二| 一区二区三区国产精品乱码| 成人无遮挡网站| 国产成人系列免费观看| 成人一区二区视频在线观看| 午夜福利在线观看吧| 欧美日本亚洲视频在线播放| 淫妇啪啪啪对白视频| 色播亚洲综合网| 波多野结衣高清无吗| 国产精品一区二区免费欧美| 91在线精品国自产拍蜜月 | 一个人免费在线观看的高清视频| 久久香蕉精品热| 免费看日本二区| 久久香蕉国产精品| 亚洲人成网站在线播放欧美日韩| 99久久精品国产亚洲精品| 18禁国产床啪视频网站| www日本在线高清视频| 无人区码免费观看不卡| 美女高潮喷水抽搐中文字幕| 久久久久国产精品人妻aⅴ院|