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

    積層混合云結(jié)構(gòu)特征及降水機理的個例模擬研究

    2015-12-06 03:05:57何暉高茜劉香娥周嵬賈星燦
    大氣科學(xué) 2015年2期

    何暉 高茜 劉香娥 周嵬 賈星燦

    ?

    積層混合云結(jié)構(gòu)特征及降水機理的個例模擬研究

    何暉1, 2高茜1, 2劉香娥1, 2周嵬1, 2賈星燦1, 2

    1北京市人工影響天氣辦公室,北京100089;2云降水物理研究與云水資源開發(fā)北京市重點實驗室,北京100089

    積層混合云是我國一種重要的降水系統(tǒng),其降水既有對流云又有層狀云特征?;诜e層混合云的重要性,本文利用中尺度數(shù)值模式WRF(Weather Research and Forecasting Model),結(jié)合三維粒子運行增長模式對2012年5月29日北京地區(qū)的一次積層混合云降過程進行了模擬研究。模擬的降水與雷達回波與實測結(jié)果基本一致。在此基礎(chǔ)上,重點分析了混合云系中積狀云與層狀云各自的微物理結(jié)構(gòu)特征與降水的發(fā)生機理等。結(jié)果表明:降水過程云內(nèi)存在著明顯的“播種—供給”機制,層狀云中“播種—供給”機制相對簡單。而對流云區(qū)中由于降水粒子可以發(fā)生上下多次的循環(huán)增長,“播種—供給”機制可在云的上下層間雙向進行,云中粒子群可以增長得更大。在積層混合云中,在低層,層狀云中已有的水凝物粒子進入內(nèi)嵌的積云塊中,而在高層水成物粒子又從積云中落到層云中,積層混合云系充分發(fā)揮了積云和層云各自的優(yōu)勢,從而降水效率較高。

    積層混合云 降水機制 播種—供給 循環(huán)增長

    1 引言

    積層混合云是由層狀云與嵌入其內(nèi)的對流云組成的云系統(tǒng)。積層混合云是我國一種重要的降水系統(tǒng),其降水既具有層狀云又兼有對流云特征:在持續(xù)低強度降水中帶有較高強度的陣性,其降水時間比較長,常常帶來大范圍的持續(xù)性或間歇性降水,可以有效地緩解干旱(李艷偉等,2009a)。大范圍的暴雨和特大暴雨基本上是由較為深厚的層狀云和嵌入內(nèi)的對流云組成的混合云產(chǎn)生的(洪延超,1996a;黃美元等,1987b)。積層混合云不僅對降水的發(fā)展過程有明顯影響,而且其動力、熱力場的結(jié)構(gòu)比層云或積云都要復(fù)雜(于翡和姚展予,2009),這是由于其積層混合的特殊結(jié)構(gòu),主要體現(xiàn)在層狀云和積狀云之間的相互影響,相互反饋以及對總體降水的影響等。

    針對積層混合云,國內(nèi)外學(xué)者開展了大量的觀測和數(shù)值模擬研究。通過觀測發(fā)現(xiàn),積層混合云中由于嵌入式對流單體中液態(tài)水含量比周圍的層云高,并且具有一定的上升氣流速度,可以產(chǎn)生更多的冰晶(Matejka et al., 1980;Hobbs and Rangno,1990;Evans et al., 2005),根據(jù)云中不同的垂直速度,往往可以增加20%~35%的降水量(Herzegh and Hobbs,1980;Houze et al. 1981;Rutledge and Hobbs,1983)。Hobbs and Locatelli(1978)指出積層混合云降水往往分布不均勻,嵌入式對流單體下方降水量一般比周圍區(qū)域高,會出現(xiàn)“雨核”的現(xiàn)象。Frederic et al.(1993)利用多種雷達資料和標準氣象儀器資料分析了層云降水對淺對流泡的影響;Sandra et al.(2005)分析熱帶洋面的積層混合云過程中的物理特征,并指出對流區(qū)是低層輻合、高層輻散。國內(nèi)的學(xué)者也對積層混合云開展了大量的觀測研究。黃美元等(1987a)利用雷達觀測資料等分析了暴雨積層混合云的基本特征及其與降水的關(guān)系;李艷偉(2009a)利用雷達、衛(wèi)星、地面、高空觀測資料,對48次貴州地區(qū)大范圍的積層混合云降水過程進行了分析歸納,并認為積層混合云系往往有2種途徑生成,分別是:積云并合擴展層化型與層云積化型積層混合云,而貴州省的48次積層混合云降水過程中,絕大多數(shù)都屬于積云并合擴展層化型積層混合云降水。

    在數(shù)值模擬研究方面,F(xiàn)rederic et al.(1993)利用數(shù)值模擬研究發(fā)現(xiàn)鑲嵌在層狀云中的對流 云發(fā)展更旺盛,生命史更長,降水效率更高,洪延超(1996a,1996b)、黃美元等(1987b)利用二維時變積云降水模式討論了積層混合云中層狀云對對流云發(fā)展和降水的影響,并得出積層混合云是一種非常有效的降水系統(tǒng)。李艷偉等(2009b)、呂玉環(huán)等(2012)利用中尺度數(shù)值模式WRF研究了積云并合的過程、產(chǎn)生并合的機理及其對降水的影響等。林磊和姚展予(2011)、于翡和姚展予(2009)等也利用中尺度模式分別對積層混合云的微物理特征、動力特征、云物理過程等進行了研究。

    上述有關(guān)積層混合云的觀測和數(shù)值模擬工作已經(jīng)取得了一些非常有意義的成果。但是上述的研究主要是集中在積層混合云的形成過程、云的熱力和動力學(xué)特征的分析方面。而在積層混合云降水的機制,尤其是積層混合云中積狀云與層狀云形成降水機制差異的分析以及兩者的相互影響等方面的研究還有所欠缺。

    本文利用中尺度模式WRF(Weather Research and Forecasting Model)結(jié)合三維粒子運行增長模式對2012年5月29日北京地區(qū)的一次積層混合云降水個例進行模擬研究。在模擬結(jié)果與實況較接近的基礎(chǔ)之上,分析了此次積層混合云形成的原因、降水效率的分布、對比分析了混合云系中積狀云與層狀云各自的微物理結(jié)構(gòu)特征與降水的發(fā)生機理等,并對云系內(nèi)對流云和層狀云的相互作用進行了研究。

    2 模擬方案和結(jié)果的初步分析

    2.1 模擬方案

    本文利用的WRF模式是由NCAR(National Center for Atmospheric Research)/NCEP(National Centers for Environmental Prediction)和FSL/NOAA(Forecast Systems Laboratory /National Oceanic and Atmospheric Administration)等聯(lián)合開發(fā)研制的中尺度數(shù)值模式。本文采用Kain-Fritsch積云對流參數(shù)化方案(Khain,2004)、Morrison 2-moment云微物理顯式方案(Morrison et al., 2009)、RRTM長波輻射計算方案(Mlawer et al., 1997)、Dudhia短波輻射計算方案(Dudhia, 1989)、the Yonsei University (YSU)PBL scheme邊界層方案(Noh et al., 2001)。WRF中尺度模式中的Morrison 2-moment方案包括了水汽、云水、雨水、云冰、雪、霰等的比含水量以及雨水、云冰、雪、霰等的比數(shù)濃度預(yù)報量。模式采用四重嵌套,格距分別是27、9、3和1 km。模式區(qū)域的中心位置是(40°N,115°E),模式在 垂直方向分27層,水平格點數(shù)分別為151×151、184×142、199×175、397×445。

    采用分辨率為1°×1°,時間間隔為6小時的NCEP再分析資料作為模式初始場和側(cè)邊界條件。模擬選用歐拉質(zhì)量坐標和Runge-Kutta 3rd時間積分方案,在格距為27和9 km的區(qū)域采用Kain- Fritsch積云對流參數(shù)化方案和Morrison 2-moment云微物理顯式方案,在格距3和1 km的區(qū)域僅采用顯式方案。模擬的開始時間為2012年5月28日20:00(北京時,下同),共模擬24小時。

    模擬云粒子在云中的運動增長行為采用的是三維粒子運行增長模式(許煥斌和段英,2001)。該模式采用全拉格朗日算法對目標粒子的運行增長進行全程追蹤,由粒子所在位置的變量值內(nèi)插決定流場、溫度、氣壓、水汽和背景含水量場,再利用這些值來計算粒子的增長,而運動軌跡則由流場和粒子的現(xiàn)時末速度來決定。微物理過程考慮了冰粒子與過冷云水、雨水的干/濕碰凍增長,以及凇附對粒子體密度和下落末速度的影響,較完整地描述了粒子的增長及密度和下落末速的變化。

    根據(jù)WRF模擬云的發(fā)展情況,選擇在5月29日12:18分別針對其中的一塊對流云區(qū)域和一塊層狀云區(qū)域運行三維粒子運行增長模式。具體的操作方法是:分別將WRF模擬的2012年5月29日12:18這兩塊區(qū)域的氣流場(水平氣流、垂直氣流)、氣壓場、高度場以及云中不同水成物含量場等輸出。按照粒子軌跡模式的要求進行不同高度層插值處理,最后形成方向為31個格點(格距為1 km),方向47個格點(格距為1 km),方向100個格點(起始高度為0.5 km,格距為0.1 km)的數(shù)據(jù)背景場,利用形成的數(shù)據(jù)背景場運行三維粒子運行增長模式120 min,對粒子運行軌跡特點進行對比研究。需要說明的是在三維粒子運行增長模式運行的120 min里,模式的初始場并沒有變化是靜態(tài)的,而實際大氣中的環(huán)境場是動態(tài)變化的,相關(guān)研究顯示在變化著的宏微觀場中粒子增長軌跡的規(guī)律性與靜態(tài)模擬的結(jié)果是相似的(許煥斌和段英,2001),因此本文中采取了計算量偏少的靜態(tài)模擬方法。

    2.2 模擬結(jié)果和觀測資料的對比分析

    2.2.1 實測和模擬的降水場

    圖1為觀測和模擬的北京地區(qū)2012年5月29日11:00~15:00和15:00~19:00降水量的對比圖。觀測的降水量來源于北京區(qū)域自動站的實測資料??梢钥闯觯M的降水位置和實測大致相同。觀測和模擬均顯示這次影響北京的降雨系統(tǒng)自西北向東南方向移動。在11:00~15:00降水基本都處于北京中部—西北部區(qū)域,最強降水量出現(xiàn)在北京中部,約為18~20 mm左右,降水帶為東北—西南向,模式模擬的強降水中心、降水量級及雨帶走向等均與實測較接近,但也可以看出模式模擬的降雨區(qū)較實測有所偏大;在15:00~19:00,觀測和模擬均顯示在北京的北部和北京的東南角出現(xiàn)了兩個強降雨中心,但模擬的北京東南角的強降水中心的降水量較實測偏大,北京東南角觀測的降水量大約為8 mm,而模擬的降水量超過了12 mm。

    2.2.2 實測和模擬的雷達回波組合反射率

    圖2為觀測和模擬的2012年5月29日13:00、14:00、15:00、16:00雷達回波組合反射率的對比圖。從實測的雷達回波圖可以看出,此次影響北京的降雨系統(tǒng)自西南向東北向移動。模擬的雷達回波也體現(xiàn)出了這個特點,并且就強雷達回波的位置來說,模擬和實測的也比較接近。在13:00,模擬和實測的強回波均出現(xiàn)在北京的西南部,到14:00模擬和實測的雷達回波均移到了北京的中部,15:00模擬和實測的雷達回波強中心都移到了北京的東北部,但模擬的回波在北京的東南部出現(xiàn)了一個強中心,而實測結(jié)果沒有。到29日16:00,模擬和實測的雷達回波均顯示以前的強回波在北京的西北部已變?nèi)酰谔旖虻奈鞅辈烤霈F(xiàn)了一個較強的回波中心,從實測與模擬的降雨量與雷達回波的對比可以看出,WRF模式對此次北京地區(qū)積層混合云的降水個例模擬的比較成功。

    圖1 觀測(左列)和模擬(右列)的2012年5月29日降水量(單位:mm):(a、b)11:00~15:00(北京時,下同);(c、d)15:00~19:00

    圖2 (a、b、c、d)觀測和(e、f、g、h)模擬的2012年5月29日雷達回波組合反射率(單位:dBZ):(a、e)13:00;(b、f)14:00;(c、g)15:00;(d、h)16:00

    3 環(huán)境對積層云系降水發(fā)展的作用

    3.1 合適的天氣條件

    圖3為模擬的2012年5月29日02:00,500、700、850 hPa高度場與風(fēng)場以及雷達回波組合反射率與地面風(fēng)場圖,由圖可見,在500、700和850 hPa均有一深厚的高空槽,高空槽后有較冷的偏北氣流,而槽前是暖濕的偏南氣流。冷暖空氣交匯,暖空氣沿著冷空氣爬升,冷空氣將暖空氣整層抬升,暖濕空氣凝結(jié)形成層狀云,導(dǎo)致在該時刻,在空中槽線對應(yīng)的位置出現(xiàn)了雷達回波(圖3d)。

    圖3 2012年5月29日02:00模擬的高度場(藍線,單位:dagpm)和風(fēng)矢量:(a)500 hPa;(b)700 hPa;(c)850 hPa;(d)雷達回波組合反射率(單位:dBZ)及地面風(fēng)場

    圖4 2012年5月29日02:00模擬的(a)800 hPa和(b)500 hPa散度場(紅色橢圓框區(qū)域代表云團所在的位置,單位:10?5 s?1)

    3.2 散度場

    5月29日02:00是云系初生的時刻,圖4為該時刻800 hPa和500 hPa散度場的分布(圖4)。從圖4可以看出,在當時的天氣條件下,云(團)附近的中低層(800 hPa附近)存在著輻合,而高層(500 hPa附近)存在著輻散,這種配置有利于環(huán)境對云(系)的發(fā)展輸送水汽和能量.并且輻合場由西向東移動,云系的移動方向與之一致。

    3.3 空氣穩(wěn)定度

    對29日02:00的不穩(wěn)定能量進行分析,圖5分別是模擬的29日02:00與03:00雷達回波組合反射率、沿02:00雷達回波組合反射率的假相當位溫垂直圖以及02:00最大對流有效位能分布圖。可以看出,在02:00對應(yīng)于此時的雷達回波區(qū)域出現(xiàn)了一定大小的對流有效位能,最強的對流有效位能為290 J kg?1。從沿02:00雷達回波組合反射率的假相當位溫垂直圖可以看出,在對應(yīng)于02:00雷達回波的A、B、C處,在假相當位溫垂直剖面圖上均出現(xiàn)了假相當位溫隨高度增加而減少的情況,表明這些區(qū)域存在對流不穩(wěn)定。從而在一小時之后的03:00雷達回波組合反射率圖上,在層狀云系的A、B、C處均出現(xiàn)了對流云塊,形成了積層混合云系,并逐漸西移至北京區(qū)域內(nèi)。

    圖5 2012年5月29日模擬的雷達回波組合反射率(單位:dBZ):(a)02:00;(b)03:00;(c)沿02:00雷達回波組合反射率DD'的假相當位溫垂直剖面(A、B、C標注的是出現(xiàn)不穩(wěn)定的區(qū)域,單位:K);(d)02:00最大對流有效位能分布(單位:J kg?1)

    4 云系不同部位降水效率的分布

    考慮到輸入計算域內(nèi)的總水汽量,一部分仍以水汽形式滯留于空中,一部分是通過凝結(jié)和凝華等過程轉(zhuǎn)化為云中各種濕物質(zhì),而這其中的一部分濕物質(zhì)會以降水的形式降落到地面。將地面降水量與凈凝結(jié)量凈凝華量之和的比定義為降水效率(Khain et al.2008; 洪延超和李宏宇,2011),即:

    =SP/(cond+sub), (1)

    其中代表降水效率,SP代表地面降水量,cond代表凈凝結(jié)量,sub代表凈凝華量。

    圖6為5月29日11:48~12:18地面降水效率的分布,在此段時間內(nèi)大部分云系的降水效率為10%~40%左右。而部分區(qū)域云系的降水效率達到50%~70%,結(jié)合12:18雷達回波組合反射率(圖15a)可以發(fā)現(xiàn),降水效率高的云系對應(yīng)的往往是積層混合云系。

    圖6 5月29日11:48~12:18地面降水效率的分布

    5 云系的微物理結(jié)構(gòu)特征及降水機理分析

    5.1 單站云微結(jié)構(gòu)的變化

    對位置為(40.1°N,116.1°E)的站點,結(jié)合模擬的雷達回波組合反射率分別選擇系統(tǒng)到來之前的12:00,處于系統(tǒng)最強時的12:30以及系統(tǒng)結(jié)束 時的13:30對該位置模擬的各種水成物含量的垂直分布進行分析(圖7)。

    從圖7可以看出,在系統(tǒng)來臨之前的12:00(圖7 a、d),該點的上空主要存在著冰晶和雪晶,冰晶的分布高度為3000~8500 m,峰值位于7000 m的高度,而雪晶的分布范圍為2500~7000 m,峰值位于4000 m的高度,雪晶位于冰晶層的下方,主要來源于冰晶的轉(zhuǎn)化,此時存在一定的雨水含量,由于此時云水含量很少,因此雨水主要來自冰雪晶的融化,此時云中主要存在冷云過程。

    到系統(tǒng)最強的12:30(圖7b、e),各種水成物含量急劇增加,云系內(nèi)雪晶、霰、云水分布的高度都非常廣,表明云中存在較強的冷暖云過程,雪晶的上下方分別為少量的冰晶和霰,霰的峰值達到了4 g kg?1,霰峰值所在的高度正好對應(yīng)雪晶的峰值所處的高度并接近于云水峰值所對應(yīng)的高度,可以看出,大量的雪經(jīng)高層降落的冰晶轉(zhuǎn)化生成后,凇附過冷水生長,并不斷向霰轉(zhuǎn)化。此時的雨水一方面來自雪晶、霰降落到暖區(qū)的融化,另一方面來自云水的轉(zhuǎn)化。

    圖7 5月29日(a、d)12:00、(b、e)12:30、(c、f)13:30雷達回波組合反射率(左列,●代表站點的位置)及相應(yīng)時刻的站點(40.1°N,116.1°E)處各種水成物含量的垂直分布(右列,qi:冰晶,qg:霰,qs:雪晶,qr:雨滴,qc:云滴;單位:g kg?1)

    到系統(tǒng)結(jié)束時的13:30(圖7c、f),各種水成物含量大幅下降,云水與霰的分布范圍也大為減小,雖然仍然有雨水的存在,但含量較系統(tǒng)最強時也大幅減小,此時云系中的冷暖云過程均較弱。

    由該站云內(nèi)各種粒子微物理結(jié)構(gòu)變化可以看出(圖8),29日12:00到13:30,該站云體分為三層,高層位于7000~9000 m,主要包括冰晶和少許的雪晶;中層位于3000~6000 m,主要包括雪晶、云水以及霰等。而低層位于3000 m以下,主要包括雨水和云水等。

    圖8 站點(40.1°N,116.1°E)處(a)冰晶、(b)雪晶、(c)云水、(d)霰、(e)雨水含量(單位:g kg?1)與該點的(f)主要云微物理轉(zhuǎn)化過程(單位:g kg?1 s?1)及(g)上升速度(單位:m s?1)和(h)雷達回波(單位:dBZ)自29日12:00到13:30的時間—高度剖面

    冰晶主要出現(xiàn)在高層,雪晶的位置出現(xiàn)在冰晶層偏下方,并且冰雪晶的分布高度隨著時間的推移有向下方移動的趨勢,可以看出云內(nèi)存在著冰、雪粒子從高層向下的“播種”過程(圖8a、b)。在12:30與13:15,由于云中上升氣流的增大(圖8g),云中水汽凝結(jié)為云水過程(Pcc)加強(圖8f)導(dǎo)致中高層云水含量較高,中心含水量值分別達到2 g kg?1與1 g kg?1。同時還可以看到在高云水區(qū)的下方雨水含量也較高,并且還出現(xiàn)了霰(圖8c、d、e)。根據(jù)云內(nèi)冰雪晶、過冷云水的配置以及云內(nèi)主要云微物理轉(zhuǎn)化過程(圖8f)可以看出,大量云水的存在也為高層“播種”下來的冰、雪粒子提供了優(yōu)越“供給”增長條件,霰主要是降落的冰晶撞凍雨水(Piacr)、雪晶撞凍云水(Psacw)、雪晶撞凍雨水(Psacr)增長轉(zhuǎn)化的結(jié)果(圖8f)。并且從圖8可以看出,在降雨最強的12:30與13:10左右,也正好是高層的冰晶和雪晶粒子,中層的霰粒子、雪晶和云水以及低層的雨水含量最多和云中上升氣流最強的時段,降雨的發(fā)生主要是高層雪晶降落到暖區(qū)的融化(Psmlt)、霰降落到暖區(qū)的融化(Pgmt)以及云水轉(zhuǎn)化為雨水(Prc)與雨水對云水的收集(Pra)等過程??梢娸^強的冷暖雨過程以及云內(nèi)存在明顯的“播種—供給”云機制導(dǎo)致了降雨量的增大。結(jié)合該站點的雷達回波時間高度剖面圖可以看到,在降雨最強的12:30與13:10左右也正是雷達回波偏強的時刻。

    隨著時間的推移,到13:30,高中空的冰晶、雪晶、霰和云水的含量已較少,可見此時冷云過程已有所減弱,“播種—供給”云機制也逐漸消失,從而導(dǎo)致了地面雨量的減少。

    5.2 對流云與層狀云中粒子運行軌跡的分析

    5.2.1 對流云粒子運行軌跡

    粒子的運行是由流場和粒子落速控制的,對流云中流場相對復(fù)雜,垂直氣流與水平氣流分布都非常不均勻,粒子在對流云場中的運行相對較為復(fù)雜。為研究對流云中的粒子運行軌跡特點,利用三維粒子運行增長模式對該次積層混合云降水個例中對流云區(qū)中的粒子運行軌跡進行模擬研究。

    圖9為WRF模擬的2012年5月29日12:18疊加風(fēng)場的組合雷達回波圖以及沿紅線的疊加流場上的云水含量垂直剖面圖。選擇的對流云區(qū)域為方框區(qū)域A(圖9a,39.9°~40.3°N,115.879°~116.243°E),將該區(qū)域模擬的相關(guān)數(shù)據(jù)輸出并進行插值處理后運行三維軌跡模式。在三維軌跡模式中,在格點[(方向)(15~25 km),(方向)(20~30 km),(方向)(6~7 km)]處(圖9b中的圓圈區(qū)域)放置示蹤粒子,一共放置1000個示蹤粒子,并對每個粒子進行編號(1~1000)。軌跡模式共運行120分鐘,步長為5秒,每五步輸出一次軌跡。

    圖9 5月29日12:18疊加風(fēng)場的(a)組合雷達回波圖以及(b)沿區(qū)域A紅線的疊加流場和溫度場(單位:°C)上的雷達回波垂直剖面,單位:dBZ

    圖10 編號為310、473以及739粒子的三維運行軌跡:(a)正面圖;(b)側(cè)面圖。紅線代表水平速度為零的零線,“×”代表粒子的起始位置

    圖10為對流云區(qū)中部分粒子的三維運行軌跡的正面和側(cè)面圖,圖11為其中編號為739粒子的三維運行軌跡的正面和側(cè)面圖,“×”代表粒子的起始位置,紅線代表水平速度為零的零線。圖12為編號為739粒子的運行軌跡在圖9紅線剖面上的投影。

    從圖10到圖12可以看出,在初始時刻由于粒子處于上升氣流的邊緣位置,上升氣流速度小于粒子落速,粒子下落,并隨著水平氣流向外移動。當粒子下降到一定高度時,由于水平氣流的方向與上層時的相反,粒子隨著水平氣流逐步進入主上升氣流區(qū),進入主上升氣流區(qū)后,由于對流云中復(fù)雜的流場,垂直氣流上下層分布不均勻,對流云中心的上升氣流旁往往有下沉氣流,有些區(qū)域上升氣流速度小于粒子落速,而有些區(qū)域上升氣流速度又大于粒子落速。在水平氣流方面,上下層的水平氣流方向也不一致,低層的水平氣流與高層水平氣流方向不一致或相反,同時隨著粒子的逐步增長,粒子的落速也在增大,粒子在云場中運行軌跡也較復(fù)雜,粒子在主上升氣流區(qū)內(nèi)呈現(xiàn)出圍繞水平速度為零的零線,上下循環(huán)往復(fù)的運行軌跡(圖10,圖11),隨著粒子的進一步增大,最終粒子從主上升氣流區(qū)的另一側(cè)下落(圖12)。

    圖11 編號為739粒子的三維運行軌跡:(a)正面圖;(b)側(cè)面圖。紅線代表水平速度為零的零線,“×”代表粒子的起始位置

    圖12 編號為739粒子的運行軌跡在圖9剖面上的投影(右圖為左圖方框區(qū)域的放大)

    5.2.2 層狀云粒子運行軌跡

    WRF模擬的2012年5月29日12:18方框區(qū)域B(圖9a,39.9°~40.3°N,115.879°~116.243°E)中,回波強度較為均勻,回波強度為20~30 dB左右,可以認為是層狀云均勻。將WRF模式的該塊區(qū)域的相關(guān)數(shù)據(jù)輸出并進行插值處理后運行軌跡模式以對層狀云中粒子運行軌跡特點進行對比分析。在三維運行增長模式中,在格點[(10~15 km),(10~15 km),(3~4 km)]處共放置250個示蹤粒子。

    圖13為層狀云中部分粒子疊加在一起的三維運行軌跡圖。圖14為標號為207粒子的運行軌跡。可以清楚得看出,相對于對流云的流場,層狀云流場要簡單得多,粒子大多呈現(xiàn)簡單的上升和下沉運動。

    通過以上分析,可以看出在一般的層狀云中,層狀云上升氣流相對較小且較為均勻,云中下部過冷水的運行軌跡是簡單的先上升再下降的曲線。在上升中,水粒子可凝結(jié)增長,在低于0°C時,一些粒子可轉(zhuǎn)化為冰晶,在高于冰面飽和的情況下,可凝華增長,在存在過冷水的云區(qū),還通過貝吉隆過程以及凇附過程增長,長大的粒子下落時,對下層云播撒,再通過碰并和凝結(jié)增長,繼續(xù)長大,掉落出云形成降水。這樣的“播撒—供給”機制相對簡單。而對流云區(qū)的氣流上升不均勻,且在對流泡旁也有下沉氣流,降水粒子的增長軌跡不再是簡單的上升,可以發(fā)生上下多次的循環(huán)增長,“播撒—供給”機制可在云的上下層間雙向進行,云中粒子群可以長得更大。

    圖13 層狀云中部分粒子的三維運行軌跡:(a)正面圖;(b)側(cè)面圖?!啊痢贝砹W拥钠鹗嘉恢?/p>

    圖14 層狀云中編號為207粒子的三維運行軌跡:(a)正面圖;(b)側(cè)面圖。“×”代表粒子的起始位置

    5.3 對流云和層狀云之間含水量的交換

    云是水凝物粒子群,它的含水量變化反映著云的變化。下面來分析一下在積層混合云中,對流云與層狀云中的水成物是如何交換的。圖15分別是12:18與12:30雷達回波組合反射率圖,在圖中方框區(qū)域,可以看到出現(xiàn)兩個對流云塊,兩個對流云塊通過層狀云相連接。我們沿線做雷達回波的垂直剖面,并同時疊加云中水成物總含量total=cloud+ice+rain+snow+graup水平通量散度剖面圖。從圖15可以看出,在兩個對流云塊的下部為水成物的輻合區(qū),對流云塊的上部為水成物的輻散區(qū),而中間的層狀云情況正好相反,層狀云的下部為水成物的輻散區(qū),層狀云的上部為水成物的輻合區(qū)。所以可以得出,在低層,層狀云中已有的水凝物粒子進入水內(nèi)嵌的積云塊中,能增長得更為充分。而在高層水成物粒子又從積云中落到層云中,此時積云充當了層云的播撒云。因此,積層混合云中發(fā)生層云—積云的粒子群交換。積層混合云系充分發(fā)揮了積云和層云各自的優(yōu)勢,降水效率較高。

    圖15 2012年5月29日(a、c)12:18、(b、d)12:30組合雷達回波圖(左)以及沿線段AB疊加在雷達回波上的云中水成物總含量的水平通量散度(c、d, 單位:10?4 kg cm?2 hPa s)垂直剖面(右)

    6 結(jié)論

    本文利用中尺度模式WRF結(jié)合三維粒子運行增長模式對2012年5月29日北京地區(qū)的一次積層混合云降水過程進行了模擬,研究了這次降水過程形成的原因、降水效率的分布、降水形成機理以及積狀云與層狀云之間含水量的交換等。研究表明:

    (1)合適的天氣條件、低層輻合高層輻散的散度場配置以及云系中存在局地對流不穩(wěn)定導(dǎo)致積層混合云系的出現(xiàn);

    (2)此次降水過程云內(nèi)存在著明顯的“播種—供給”機制。降水一方面來自于云雨的轉(zhuǎn)化,另一方面來自冰相粒子的融化。降雨的最強時段正是云中“播種—供給”機制最強的時段。

    (3)在降水過程云內(nèi)存在著“播種—供給”機制是常見的,在層積混合云中“播種—供給”機制的具體表現(xiàn)如何呢?運用粒子運行增長模式模擬可以發(fā)現(xiàn):首先,層狀云上升氣流相對較小且較為均勻,云粒子的運行軌跡是簡單的先上升再下降的曲線,“播種—供給”機制相對簡單;其次,對流云區(qū)的氣流上升不均勻,降水粒子的可以發(fā)生上下多次的循環(huán)增長,“播種—供給”機制可在云的上下層間雙向進行,云中粒子群可以長得更大。

    (4)在積層混合云中,在低層,層狀云中已有的水凝物粒子進入內(nèi)嵌的積云塊中,能增長得更為充分。而在高層,水成物粒子又從積云中落到層云中,此時積云充當了層云的播撒云。因此,積層混合云中發(fā)生層云—積云的粒子群交換。積層混合云系充分發(fā)揮了積云和層云各自的優(yōu)勢,降水效率較高。

    (References:)

    Dudhia, J. 1989. Numerical study of convection observed during the winter monsoon experiment using a mesoscale two-dimensional model [J], J. Atmos. Sci., 46 (2): 3077–3107.

    Evans A G, Locatelli J D, Stoelinga M T, et al. 2005. The IMPROVE-1 storm of 1–2 February 2001. Part II: Cloud structures and the growth of precipitation [J]. J. Atmos. Sci., 62 (10): 3456–3473.

    Frederic F, Zawadzki I, Cohn S. 1993. The influence of stratiform precipitation on shallow convective rain: A case study [J]. Mon. Wea. Rev., 121 (12): 3312–3325.

    黃美元, 洪延超, 吳玉霞. 1987a. 梅雨鋒云系和降水的若干研究 [J]. 大氣科學(xué), 11 (1): 23–30. Huang M Y,Hong Y C, Wu Y X. 1987a. Some studies on meiyu frontal cloud system and its precipitation [J]. Chinese Journal of Atmospheric Sciences (in Chinese) , 11 (1): 23?30.

    黃美元, 洪延超, 徐華英, 等. 1987b. 層狀云對積云發(fā)展和降水的影響——一種云與云之間影響的數(shù)值模擬 [J]. 氣象學(xué)報, 45 (1): 72–77. Huang M Y,Hong Y C, Xu H Y.1987b. The effects of the existence of stratiform cloud on the development of cumulus cloud and its precipitation [J]. Acta Meteorologica Sinica (in Chinese), 45 (1):72~ 77.

    Khain A P, Benmoshe N, Pokrovsky A. 2008. Factors determining the impact of aerosols on surface precipitation from clouds: An attempt at classi?cation [J]. J. Atmos. Sci., 65: 1721–1748.

    Herzegh P H, Hobbs P V. 1980. The mesoscale and microscale structure and organization of clouds and precipitation in midlatitude cyclones. II: Warm-Frontal clouds [J]. J. Atmos. Sci., 37: 597–611.

    Hobbs P V, Locatelli J D. 1978. Rainbands, precipitation cores and generating cells in a cyclonic storm [J]. J. Atmos. Sci., 35: 230–241.

    Hobbs P V, Rangno A L. 1990. Rapid development of high ice particle concentrations in small polar maritime cumuliform clouds [J]. J. Atmos. Sci., 47: 2710–2722.

    洪延超. 1996a. 積層混合云數(shù)值模擬研究 (I)——模式及其微物理過程參數(shù)化 [J]. 氣象學(xué)報, 54 (5): 544–557. Hong YC. 1996a. The numerical simulation on study of convective and stratiform mixed clouds, part(I) ─The model and paramterization of microphysical processes [J]. Acta Meteorologica Sinica (in Chinese), 54 (5): 544?557.

    洪延超. 1996b. 積層混合云數(shù)值模擬研究 (Ⅱ)──云相互作用及暴雨產(chǎn)生機制 [J]. 氣象學(xué)報, 54 (6): 661–674. Hong Y C. 1996b. The numerical simulation on study of convective and stratiform mixed clouds, part (Ⅱ) ─Intercation of clouds and formative mechanism of the heavy rain [J]. Acta Meteorologica Sinica (in Chinese), 54 (6): 661?674

    洪延超, 李宏宇. 2011. 一次鋒面層狀云云系結(jié)構(gòu)降水機制及人工增雨條件研究 [J]. 高原氣象, 30 (5): 1308–1323. Hong Yanchao, Li Hongyu. 2011. Cloud structure, precipitation Mechanism and artifical enhancement precipitation condition frontal stratiform cloud system [J]. Plateau Meteorology (in Chinese), 30 (5):1308?1323.

    Houze R A Jr, Rutledge S A, Matejka T J, et al. 1981. The mesoscale and microscale structure and organization of clouds and precipitation in midlatitude cyclones. III: Air motions and precipitation growth in a warm-frontal rainband [J]. J. Atmos. Sci., 38: 639–649.

    李艷偉, 牛生杰, 羅寧, 等. 2009a. 積云并合擴展層化型積層混合云的數(shù)值模擬分析 [J]. 地球物理學(xué)報, 52 (5): 1165–1175. Li Y W, Niu S J, Luo N, et al. 2009a. Numerical simulation about mixture of convective and stratiform clouds formed ny convection merger [J]. Chinese J. Geophys. (in Chinese), 52 (5): 1165?1175.

    李艷偉, 牛生杰, 姚展予, 等. 2009b. 云并合的初始位置探討 [J]. 大氣科學(xué), 33 (5): 1015–1026. Li Yanwei, Niu Shengjie, Yao Zhanyu , et al. 2009b. Initial part discussion of cl oud merger [J]. Chinese Journal of Atmosph eric Sciences (in Chinese), 33 (5): 1015?1026.

    林磊, 姚展予. 2011. 華北地區(qū)一次積層混合云降水的數(shù)值模擬研究 [J]. 氣象, 37 (12): 1473–1480. Lin Lei, Yao Zhanyu. 2011. Numerical study on the complex of the stratiform and embedded convective cloud precipitataion of North China [J]. Meteorological Monthly (in Chinese), 37 (12): 1473?1480.

    呂玉環(huán), 李艷偉, 金蓮姬, 等. 2012. 云并合過程中物理特征演變的模擬研究 [J]. 大氣科學(xué), 36 (3): 471–486. Lü Yuhuan, Li Yanwei, Jin Lianji, et al. 2012.Simulation of physical characteristics evolution in cloud merger [J].Chinese Journal of Atmospheric Sciences (in Chinese) 36 (3): 471?486.

    Matejka T J. Houze R A, Hobbs P V. 1980. Microphysics and dynamics of clouds associated with mesoscale rainbands in extratropical cyclones [J]. Quart. J. Roy. Meteor. Soc., 106 (447): 29–56.

    Morrison, H., J. Curry, et al. 2005. A new double-moment microphysics parameterization for application in cloud and climate models. Part I: Description [J]. J. Atmos. Sci., 62(6): 1665?1677.

    Noh Y, Cheon W G, Raasch, S. 2001. The improvement of the K-profile model for the PBL using LES. Preprints, Int. Workshop of Next Generation NWP Models [J]. Laboratory for Atmospheric Modeling Research, 12 (6): 65–66.

    Rutledge S A, Hobbs P V. 1983. The mesoscale and microscale structure and organization of clouds and precipitation in midlatitude cyclones. VIII: A model for the “Seeder-Feeder” process in warm-frontal rainbands [J]. J. Atmos. Sci., 40: 1185–1206.

    Sandra E Y, Robert A, Houze J R, et al. 2005. Physical characterization of tropical oceanic convection observed in KWAJEX [J]. J. Appl. Meteor., 44 (4): 385–415.

    許煥斌, 段英. 2001. 冰雹形成機制的研究并論人工雹胚與自然雹胚的“利益競爭”防雹假說 [J]. 大氣科學(xué), 25 (2): 277–288. Xu Huanbin, Duan Ying. 2001. The mechanism of hailstone’s formation and the hail-suppression hypothesis: “Beneficial competition” [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 25 (2): 277?288.

    于翡, 姚展予. 2009. 一次積層混合云降水實例的數(shù)值模擬分析 [J]. 氣象, 35 (12): 3?13. Yu Fei, Yao Zhanyu. 2009. Numerical study on the complex of the stratiform and embedded convective cloud precipitataion: A case study [J]. Meteorological Monthly(in Chinese), 35 (12): 3?13.

    何暉,高茜,劉香娥,等. 2015. 積層混合云結(jié)構(gòu)特征及降水機理的個例模擬研究[J]. 大氣科學(xué), 39 (2): 315?328, doi: 10.3878/j.issn.1006-9895.1404. 14102. He Hui, Gao Qian, Liu Xiang’e, et al. 2015. Numerical simulation of the structural characteristics and precipitation mechanism of stratiform clouds with embedded convections[J]. Chinese Journal of Atmospheric Sciences (in Chinese), 39 (2): 315?328.

    Numerical Simulation of the Structural Characteristics and Precipitation Mechanism of Stratiform Clouds with Embedded Convections

    HE Hui1, 2, GAO Qian1, 2, LIU Xiang’e1, 2, ZHOU Wei1, 2, and JIA Xingcan1, 2

    1,100089;2,,100089

    Stratiform clouds with embedded convections represent an important precipitation system in China. Precipitation from this system has characteristics of both convective and stratiform cloud. We studied this important process on May 29, 2012 using a mesoscale Weather Research and Forecasting (WRF) model and a three-dimension growth model. Since the resulting simulation of rainfall and radar echo were similar to observation results, we then analyzed the microphysical structural characteristics and precipitation mechanism of stratiform and convective clouds, respectively, in stratiform clouds with embedded convections. Results show that there is a “seeder-feeder” mechanism in this precipitation system, and this mechanism is simple in stratiform cloud. In convective cloud, however, because of the frequent circular particle movement between the upper and lower levels, this mechanism can occur in both the upper and lower cloud levels and so the particles can grow bigger. At low levels of the cloud system, the particles move into the embedded convective cloud, while at high levels, the particles fall from the convective cloud to the stratiform cloud. Stratiform clouds with embedded convections give full play to the advantages of convective and stratiform cloud, so their precipitation efficiency is high.

    Embedded convective cloud, Precipitation mechanism, Seeder-feeder, Circular growth

    1006-9895(2015)02-0315-14

    P401

    A

    10.3878/j.issn.1006-9895.1404.14102

    2014-01-01;網(wǎng)絡(luò)預(yù)出版日期2014-04-28

    公益性行業(yè)(氣象)科研專項GYHY201306065,北京市自然科學(xué)基金項目6154027,國家自然科學(xué)基金項目41205100、41175007

    何暉,男,1977年出生,博士,主要從事云降水物理研究。E-mail: hehui@bjmb.gov.cn

    亚洲欧美成人综合另类久久久| 国产黄色免费在线视频| 午夜老司机福利片| 国产1区2区3区精品| 免费黄色在线免费观看| 精品国产一区二区三区四区第35| 欧美日韩亚洲综合一区二区三区_| 亚洲人成电影观看| 一级片'在线观看视频| 好男人视频免费观看在线| 侵犯人妻中文字幕一二三四区| 欧美精品高潮呻吟av久久| 亚洲在久久综合| 日韩av在线免费看完整版不卡| 七月丁香在线播放| 女人爽到高潮嗷嗷叫在线视频| 国产无遮挡羞羞视频在线观看| 亚洲精品久久成人aⅴ小说| 女人久久www免费人成看片| 国产一区二区三区av在线| 欧美少妇被猛烈插入视频| 热99国产精品久久久久久7| 亚洲精品一区蜜桃| 五月开心婷婷网| 汤姆久久久久久久影院中文字幕| 国产片内射在线| 亚洲一区中文字幕在线| 精品亚洲成国产av| 亚洲欧美一区二区三区黑人| 日本猛色少妇xxxxx猛交久久| 99久国产av精品国产电影| 人人妻,人人澡人人爽秒播 | 18禁国产床啪视频网站| 丝袜人妻中文字幕| av一本久久久久| 国产精品99久久99久久久不卡 | 9色porny在线观看| 亚洲av综合色区一区| 午夜福利乱码中文字幕| 这个男人来自地球电影免费观看 | 九草在线视频观看| 国产麻豆69| 国精品久久久久久国模美| 韩国高清视频一区二区三区| 男女床上黄色一级片免费看| 婷婷色综合大香蕉| 国产免费视频播放在线视频| 久久久久久久国产电影| 精品酒店卫生间| 欧美人与性动交α欧美精品济南到| 国产精品一区二区在线不卡| 蜜桃国产av成人99| 欧美激情极品国产一区二区三区| 街头女战士在线观看网站| 欧美亚洲日本最大视频资源| 成年人免费黄色播放视频| 青草久久国产| 一区二区三区四区激情视频| 国产爽快片一区二区三区| 老汉色av国产亚洲站长工具| 视频区图区小说| 亚洲,一卡二卡三卡| 黄色 视频免费看| 国产毛片在线视频| 操出白浆在线播放| 一区二区三区乱码不卡18| 亚洲av成人不卡在线观看播放网 | 欧美黄色片欧美黄色片| 亚洲一卡2卡3卡4卡5卡精品中文| xxxhd国产人妻xxx| 亚洲精品一区蜜桃| 高清黄色对白视频在线免费看| 爱豆传媒免费全集在线观看| 亚洲欧美精品综合一区二区三区| 高清视频免费观看一区二区| 国产麻豆69| av国产久精品久网站免费入址| 国产在线一区二区三区精| 欧美日本中文国产一区发布| 精品一区二区三区四区五区乱码 | av女优亚洲男人天堂| 男女无遮挡免费网站观看| 亚洲七黄色美女视频| 国产精品亚洲av一区麻豆 | 视频在线观看一区二区三区| 综合色丁香网| 黑人巨大精品欧美一区二区蜜桃| 人体艺术视频欧美日本| 免费av中文字幕在线| 国产精品免费视频内射| 午夜福利一区二区在线看| 男女边摸边吃奶| 免费久久久久久久精品成人欧美视频| 日韩一区二区三区影片| 欧美激情极品国产一区二区三区| 亚洲精品乱久久久久久| 一级毛片我不卡| 91成人精品电影| 在线天堂最新版资源| 国产精品亚洲av一区麻豆 | 国产精品偷伦视频观看了| 国产一区二区 视频在线| 麻豆精品久久久久久蜜桃| 久久av网站| 免费观看a级毛片全部| 天美传媒精品一区二区| 熟妇人妻不卡中文字幕| 国产有黄有色有爽视频| svipshipincom国产片| 97精品久久久久久久久久精品| 蜜桃在线观看..| 欧美日韩成人在线一区二区| 在线观看免费高清a一片| 九九爱精品视频在线观看| 久久久国产精品麻豆| 精品一区二区三区四区五区乱码 | 一级黄片播放器| 男女无遮挡免费网站观看| 国产成人精品在线电影| 这个男人来自地球电影免费观看 | 免费在线观看黄色视频的| 亚洲国产欧美在线一区| 中文字幕精品免费在线观看视频| 精品亚洲乱码少妇综合久久| 日本黄色日本黄色录像| 亚洲激情五月婷婷啪啪| 亚洲精品国产一区二区精华液| 人妻人人澡人人爽人人| 极品少妇高潮喷水抽搐| 亚洲中文av在线| 日本一区二区免费在线视频| av国产精品久久久久影院| 亚洲欧美激情在线| 国产av一区二区精品久久| a级毛片在线看网站| 一边亲一边摸免费视频| 啦啦啦视频在线资源免费观看| 秋霞在线观看毛片| tube8黄色片| 亚洲久久久国产精品| 性高湖久久久久久久久免费观看| 啦啦啦视频在线资源免费观看| 亚洲欧美成人精品一区二区| 日本一区二区免费在线视频| 国产 精品1| 国产不卡av网站在线观看| 别揉我奶头~嗯~啊~动态视频 | 中文字幕色久视频| 超碰97精品在线观看| 国产成人精品福利久久| 精品午夜福利在线看| 少妇人妻精品综合一区二区| 在线精品无人区一区二区三| 久久这里只有精品19| 天天躁日日躁夜夜躁夜夜| 搡老岳熟女国产| 黄网站色视频无遮挡免费观看| 在线天堂中文资源库| 操美女的视频在线观看| 麻豆乱淫一区二区| 午夜福利在线免费观看网站| 欧美少妇被猛烈插入视频| 美女视频免费永久观看网站| 男女之事视频高清在线观看 | 中文字幕av电影在线播放| 国产乱人偷精品视频| 满18在线观看网站| 日本av手机在线免费观看| 国产在线一区二区三区精| 成人毛片60女人毛片免费| 久久久久人妻精品一区果冻| 国产日韩欧美亚洲二区| 久久免费观看电影| 亚洲婷婷狠狠爱综合网| 国产又色又爽无遮挡免| 99热网站在线观看| 国产精品一二三区在线看| 免费高清在线观看日韩| 国产精品免费大片| 啦啦啦啦在线视频资源| 免费高清在线观看日韩| 9色porny在线观看| 亚洲欧美色中文字幕在线| 久久精品亚洲av国产电影网| 中文字幕精品免费在线观看视频| 看免费成人av毛片| 熟女少妇亚洲综合色aaa.| 亚洲精品国产区一区二| 最近最新中文字幕大全免费视频 | 欧美成人精品欧美一级黄| 在线天堂最新版资源| 国产亚洲最大av| 久久久久久久国产电影| 欧美精品人与动牲交sv欧美| 欧美日韩亚洲国产一区二区在线观看 | 啦啦啦啦在线视频资源| 亚洲av在线观看美女高潮| 亚洲国产精品999| 交换朋友夫妻互换小说| 久久久久精品久久久久真实原创| 欧美久久黑人一区二区| 成年人午夜在线观看视频| 亚洲图色成人| 在线精品无人区一区二区三| 久久免费观看电影| 女人精品久久久久毛片| 久久人人爽av亚洲精品天堂| 国产精品久久久久久精品古装| 午夜福利乱码中文字幕| 亚洲av电影在线进入| 青草久久国产| 亚洲国产中文字幕在线视频| 午夜日韩欧美国产| 三上悠亚av全集在线观看| 蜜桃在线观看..| 日韩一本色道免费dvd| 婷婷色综合大香蕉| 日本av手机在线免费观看| 欧美日韩av久久| 黑人巨大精品欧美一区二区蜜桃| 亚洲人成77777在线视频| av.在线天堂| 国产欧美日韩综合在线一区二区| 国产精品三级大全| 人人妻人人澡人人看| 青青草视频在线视频观看| 免费在线观看视频国产中文字幕亚洲 | 在线免费观看不下载黄p国产| 考比视频在线观看| 91国产中文字幕| www.熟女人妻精品国产| 丝袜脚勾引网站| 亚洲精品久久成人aⅴ小说| 精品国产一区二区三区久久久樱花| 欧美日韩视频精品一区| 久久99热这里只频精品6学生| 亚洲三区欧美一区| 亚洲人成网站在线观看播放| 亚洲,欧美,日韩| 亚洲美女黄色视频免费看| 久久久久精品国产欧美久久久 | 岛国毛片在线播放| 欧美日韩视频精品一区| 日韩成人av中文字幕在线观看| 99久久99久久久精品蜜桃| 女性被躁到高潮视频| 丰满饥渴人妻一区二区三| 国产淫语在线视频| 亚洲视频免费观看视频| av有码第一页| 99精国产麻豆久久婷婷| 麻豆av在线久日| 国产精品av久久久久免费| 久久人妻熟女aⅴ| 只有这里有精品99| 亚洲欧美成人综合另类久久久| 色吧在线观看| svipshipincom国产片| 国产视频首页在线观看| 亚洲欧洲精品一区二区精品久久久 | 国产av码专区亚洲av| 欧美在线一区亚洲| 国产黄频视频在线观看| 国产亚洲精品第一综合不卡| av.在线天堂| 爱豆传媒免费全集在线观看| 精品少妇内射三级| 新久久久久国产一级毛片| 国产高清不卡午夜福利| av国产精品久久久久影院| av电影中文网址| 一二三四在线观看免费中文在| 搡老岳熟女国产| 国产成人精品在线电影| 大香蕉久久成人网| www.熟女人妻精品国产| 黄网站色视频无遮挡免费观看| 亚洲第一av免费看| 欧美日韩一区二区视频在线观看视频在线| 国产熟女午夜一区二区三区| av网站在线播放免费| 一本久久精品| 丝袜在线中文字幕| 免费在线观看完整版高清| 久久久久久久久久久久大奶| 日本猛色少妇xxxxx猛交久久| 亚洲精品日本国产第一区| 亚洲美女搞黄在线观看| 亚洲,欧美,日韩| av网站免费在线观看视频| 亚洲伊人久久精品综合| 中文字幕av电影在线播放| 韩国高清视频一区二区三区| 国产在线一区二区三区精| 久久久精品免费免费高清| 国产成人精品福利久久| 高清黄色对白视频在线免费看| 国产精品成人在线| 极品人妻少妇av视频| av网站在线播放免费| 黄色怎么调成土黄色| 色婷婷久久久亚洲欧美| 欧美成人精品欧美一级黄| 巨乳人妻的诱惑在线观看| 久久av网站| 国产av码专区亚洲av| 最新的欧美精品一区二区| 日韩精品有码人妻一区| 男女边吃奶边做爰视频| 亚洲国产最新在线播放| 三上悠亚av全集在线观看| 亚洲视频免费观看视频| 满18在线观看网站| 日韩大片免费观看网站| 少妇猛男粗大的猛烈进出视频| 丰满饥渴人妻一区二区三| 人人妻人人爽人人添夜夜欢视频| 午夜av观看不卡| 国产高清国产精品国产三级| 老司机在亚洲福利影院| 美女午夜性视频免费| 中文字幕人妻丝袜制服| 色吧在线观看| 国产1区2区3区精品| 少妇人妻精品综合一区二区| 亚洲国产看品久久| 精品国产一区二区三区久久久樱花| 观看av在线不卡| 成人漫画全彩无遮挡| 大香蕉久久成人网| 亚洲国产av新网站| 国精品久久久久久国模美| 18禁动态无遮挡网站| 看免费成人av毛片| 如日韩欧美国产精品一区二区三区| 国产成人91sexporn| 色视频在线一区二区三区| 狂野欧美激情性bbbbbb| 久久久国产欧美日韩av| 国产精品.久久久| 91国产中文字幕| 久久久国产一区二区| 国语对白做爰xxxⅹ性视频网站| 在线亚洲精品国产二区图片欧美| 国产伦理片在线播放av一区| 母亲3免费完整高清在线观看| 七月丁香在线播放| 嫩草影院入口| 欧美成人精品欧美一级黄| 尾随美女入室| 亚洲精品视频女| www.精华液| 9热在线视频观看99| 国产不卡av网站在线观看| 亚洲精品,欧美精品| 欧美日韩视频高清一区二区三区二| 婷婷色av中文字幕| 精品酒店卫生间| 汤姆久久久久久久影院中文字幕| www.熟女人妻精品国产| 丝袜美足系列| 一区二区三区激情视频| 一级爰片在线观看| 欧美人与善性xxx| 操出白浆在线播放| 2021少妇久久久久久久久久久| 久久精品熟女亚洲av麻豆精品| 亚洲精品,欧美精品| 国产99久久九九免费精品| 18禁观看日本| 成人亚洲精品一区在线观看| 日韩一卡2卡3卡4卡2021年| 超碰成人久久| avwww免费| 亚洲精品日韩在线中文字幕| 一级片'在线观看视频| 国产无遮挡羞羞视频在线观看| 黄色一级大片看看| 亚洲av男天堂| 777久久人妻少妇嫩草av网站| 亚洲,欧美精品.| 亚洲av成人精品一二三区| 亚洲精品第二区| 亚洲国产精品999| av免费观看日本| 国产1区2区3区精品| 亚洲图色成人| 丝袜美腿诱惑在线| 成年女人毛片免费观看观看9 | 黄色 视频免费看| 亚洲国产日韩一区二区| 嫩草影院入口| 亚洲国产精品成人久久小说| 日韩欧美一区视频在线观看| 免费人妻精品一区二区三区视频| 天堂8中文在线网| 蜜桃国产av成人99| 国产成人精品福利久久| 亚洲国产av影院在线观看| av电影中文网址| 青春草视频在线免费观看| 老鸭窝网址在线观看| 亚洲精品国产色婷婷电影| 午夜91福利影院| 欧美亚洲日本最大视频资源| 日韩,欧美,国产一区二区三区| 国产极品粉嫩免费观看在线| 精品一区二区三区av网在线观看 | 黄色一级大片看看| 日韩制服骚丝袜av| 国产免费一区二区三区四区乱码| 亚洲国产毛片av蜜桃av| 亚洲一码二码三码区别大吗| 夫妻午夜视频| 亚洲精品中文字幕在线视频| 亚洲国产日韩一区二区| 欧美日韩一级在线毛片| 亚洲色图综合在线观看| av国产精品久久久久影院| 伊人久久大香线蕉亚洲五| 国产av码专区亚洲av| 日韩不卡一区二区三区视频在线| 午夜久久久在线观看| 妹子高潮喷水视频| 国产淫语在线视频| 中文字幕人妻丝袜制服| 一级片'在线观看视频| 夜夜骑夜夜射夜夜干| 国产免费又黄又爽又色| 欧美激情高清一区二区三区 | 精品久久久精品久久久| 国产欧美日韩一区二区三区在线| 国产午夜精品一二区理论片| 永久免费av网站大全| 国产高清国产精品国产三级| 少妇人妻精品综合一区二区| 深夜精品福利| 婷婷色麻豆天堂久久| 亚洲精品国产色婷婷电影| 国产一级毛片在线| 深夜精品福利| 男人操女人黄网站| 精品少妇内射三级| 亚洲美女搞黄在线观看| 欧美黑人精品巨大| 国产又色又爽无遮挡免| 成人手机av| 国产在视频线精品| 1024香蕉在线观看| 最新在线观看一区二区三区 | 赤兔流量卡办理| 国产黄频视频在线观看| 高清不卡的av网站| 18禁动态无遮挡网站| 亚洲成人av在线免费| 亚洲国产中文字幕在线视频| 狂野欧美激情性bbbbbb| 国产又色又爽无遮挡免| 午夜福利,免费看| 日本色播在线视频| 丰满乱子伦码专区| 777久久人妻少妇嫩草av网站| 日韩欧美一区视频在线观看| 制服诱惑二区| 一级毛片电影观看| 亚洲五月色婷婷综合| 亚洲精品一二三| 亚洲欧洲国产日韩| 99久久99久久久精品蜜桃| 精品视频人人做人人爽| 母亲3免费完整高清在线观看| 一级毛片电影观看| 在线亚洲精品国产二区图片欧美| 亚洲国产日韩一区二区| 男女边摸边吃奶| 久久久久久久久免费视频了| 黄色 视频免费看| 亚洲欧美精品自产自拍| 伦理电影免费视频| 丝袜脚勾引网站| 伦理电影大哥的女人| 免费久久久久久久精品成人欧美视频| 国产精品久久久久成人av| 日本欧美视频一区| 视频区图区小说| 午夜免费鲁丝| 国产野战对白在线观看| 国产精品 国内视频| 热re99久久国产66热| 日本欧美视频一区| 少妇被粗大猛烈的视频| 亚洲婷婷狠狠爱综合网| 亚洲伊人色综图| 亚洲精品日韩在线中文字幕| 校园人妻丝袜中文字幕| 亚洲精品一区蜜桃| 成年美女黄网站色视频大全免费| 日韩av免费高清视频| 伊人久久国产一区二区| 国产精品成人在线| 亚洲精品第二区| 国产成人午夜福利电影在线观看| 1024视频免费在线观看| 午夜日韩欧美国产| 亚洲国产看品久久| 国产极品粉嫩免费观看在线| 日韩免费高清中文字幕av| 午夜福利乱码中文字幕| 亚洲欧洲国产日韩| 色婷婷久久久亚洲欧美| 在线天堂中文资源库| 成人国产麻豆网| 亚洲四区av| 九草在线视频观看| 精品久久蜜臀av无| 纯流量卡能插随身wifi吗| 人体艺术视频欧美日本| 中国三级夫妇交换| 亚洲精品中文字幕在线视频| 最近中文字幕2019免费版| 两个人看的免费小视频| 亚洲国产成人一精品久久久| 欧美日韩亚洲综合一区二区三区_| 久久青草综合色| 国产淫语在线视频| 99热网站在线观看| 亚洲欧美一区二区三区国产| 国产不卡av网站在线观看| 最近最新中文字幕免费大全7| 亚洲熟女毛片儿| 国产成人精品久久久久久| 亚洲色图综合在线观看| 另类亚洲欧美激情| av福利片在线| 久久99热这里只频精品6学生| 最近中文字幕2019免费版| 亚洲免费av在线视频| 老司机靠b影院| 男女下面插进去视频免费观看| 热re99久久国产66热| 国产精品久久久人人做人人爽| 国产97色在线日韩免费| 免费黄色在线免费观看| 男女之事视频高清在线观看 | 国产精品香港三级国产av潘金莲 | 一本一本久久a久久精品综合妖精| 精品一区二区三区四区五区乱码 | 国产精品免费视频内射| 欧美激情高清一区二区三区 | 国产成人系列免费观看| av在线观看视频网站免费| 五月天丁香电影| av有码第一页| 成人18禁高潮啪啪吃奶动态图| 亚洲图色成人| 香蕉国产在线看| 久久久久网色| 久久精品人人爽人人爽视色| 午夜免费鲁丝| 国产亚洲av片在线观看秒播厂| 日本av免费视频播放| 欧美日韩av久久| 久久久久久久久久久免费av| 少妇人妻 视频| 如何舔出高潮| 交换朋友夫妻互换小说| 我要看黄色一级片免费的| www.av在线官网国产| 国产日韩一区二区三区精品不卡| 国产色婷婷99| 国产一区亚洲一区在线观看| 丰满饥渴人妻一区二区三| 伊人亚洲综合成人网| 最近中文字幕高清免费大全6| tube8黄色片| 成人国产av品久久久| 电影成人av| 黄色毛片三级朝国网站| 捣出白浆h1v1| 欧美中文综合在线视频| 欧美在线黄色| 又大又黄又爽视频免费| 国产一级毛片在线| av网站免费在线观看视频| 亚洲欧美清纯卡通| 在线天堂中文资源库| 青青草视频在线视频观看| 97在线人人人人妻| 亚洲av成人精品一二三区| 99热全是精品| 欧美日韩综合久久久久久| 亚洲av电影在线观看一区二区三区| 亚洲成人国产一区在线观看 | 欧美激情极品国产一区二区三区| 国产高清国产精品国产三级| 国产淫语在线视频| 中文字幕色久视频| 丝瓜视频免费看黄片| 久久99精品国语久久久| 午夜免费鲁丝| 国产毛片在线视频| 伦理电影免费视频| 两个人免费观看高清视频| 大片电影免费在线观看免费| 美女午夜性视频免费| 日本色播在线视频| 悠悠久久av| 国产伦人伦偷精品视频| 国产免费视频播放在线视频| 人妻 亚洲 视频| 男女边吃奶边做爰视频| 亚洲欧美色中文字幕在线| av又黄又爽大尺度在线免费看| 老熟女久久久| av片东京热男人的天堂|