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

    一次積層混合云降水不同尺度結(jié)構(gòu)的數(shù)值模擬

    2014-09-09 07:30:12張微周毓荃
    大氣科學學報 2014年4期

    張微,周毓荃

    (1.南京信息工程大學,江蘇 南京 210044;2.中國氣象科學研究院,北京 100081)

    一次積層混合云降水不同尺度結(jié)構(gòu)的數(shù)值模擬

    張微1,周毓荃2

    (1.南京信息工程大學,江蘇 南京 210044;2.中國氣象科學研究院,北京 100081)

    利用中尺度數(shù)值模式WRF-ARW(V3.2)對2009年4月18—19日發(fā)生在張家口地區(qū)的一次積層混合云降水進行了模擬,并結(jié)合觀測資料從不同尺度對這次降水過程進行了對比分析。結(jié)果表明:700 hPa西風槽、850 hPa低渦是影響這次降水的主要天氣系統(tǒng),來自南方的暖濕空氣和西北內(nèi)蒙古低渦帶來的水汽是這次降水的主要水汽來源,兩股水汽在張家口附近低層出現(xiàn)了大尺度輻合,有利于該地區(qū)云系的發(fā)展、降水的形成;降水云系呈東北—西南向帶狀分布,帶長約1 000 km,帶寬300 km,在大片的云帶中分布著很多個小的高值中心,中心區(qū)域一般在幾十千米;結(jié)合雷達回波可以看到在均勻的回波層中鑲嵌著柱狀對流回波,具有典型的積層混合云降水回波特征;沿著雷達回波做剖面,發(fā)現(xiàn)云中云水含量分布無論是水平方向還是垂直方向都是不均勻的,雨水的大值中心與上層的霰、雪的大值中心相對應,中心水平范圍在10~20 km。

    WRF模式;積層混合云;不同尺度結(jié)構(gòu);數(shù)值模擬

    0 引言

    我國全境年平均降水量約為630 mm,人均水資源占有量不到世界平均水平的四分之一,是世界上13個貧水國家之一。其中華北、西北部分地區(qū)降水量偏少,分布很不均勻。積層混合云是一種重要的降水云型,在梅雨鋒、冷鋒、地形云降水等多種天氣過程中都會出現(xiàn)(許梓秀和王鵬云,1989;Fabry et al.,1993;劉黎平等,2004;Anagnostou,2004;宮福久等,2005;Fuhrer and Schaer,2005)。很多大范圍的降水都是由較為深厚的層狀云和嵌入其中的對流云組成的混合云產(chǎn)生的。因此,研究積層混合云自然降水過程的宏、微觀特征,對進一步開發(fā)利用空中云水資源,提高人工增雨作業(yè)的科學水平都有重要意義。

    數(shù)值模擬是研究云降水系統(tǒng)的重要手段。20世紀80年代后中尺度數(shù)值模式得到很大發(fā)展,國內(nèi)外許多學者采用數(shù)值模式對積層混合云降水進行研究。Yuter et al.(2005)利用衛(wèi)星資料等多種資料,結(jié)合數(shù)值模擬結(jié)果,分析了熱帶大洋上的積層混合云系中的層狀云區(qū)與對流云區(qū)的熱力、動力結(jié)構(gòu);結(jié)果表明,熱帶洋面上的積層混合云中,對流泡附近高層輻散、低層輻合。Orville et al.(1987)利用數(shù)值模式,模擬了穩(wěn)定層狀云中催化產(chǎn)生嵌入對流的可能性。Fuhrer and Schaer(2002)在觀測的基礎(chǔ)上,利用云數(shù)值模式,模擬了從均勻的層云降水發(fā)展至嵌有對流單體的層云降水的轉(zhuǎn)變過程,結(jié)果表明,由于對流云的動力作用,雨強、降水量明顯增加,在均勻上升氣流的條件下,嵌入的對流泡表現(xiàn)為多單體形式。洪延超(1996a,1996b)用積云對流速度場疊加輻合場的方法,建立了一個二維平面對稱積層混合云數(shù)值模式,用于模擬研究層狀云和嵌入其內(nèi)的對流云組成的混合云,并用該模式及暴雨云的平均大氣層結(jié)模擬研究了暴雨積層混合云的演變過程、兩種云的相互作用、云體結(jié)構(gòu)及降水特征,分析了暴雨產(chǎn)生的物理原因。王維佳和陶遐齡(2006,2007)建立了一維半地形積層混合云模式,用于對積層混合云降水過程的研究。鄒倩等(2008)對三維非靜力中尺度模式ARPS的云微物理方案進行了改進,利用改進后的模式模擬了華北地區(qū)一次積層混合云降水個例,通過對模擬結(jié)果的分析并結(jié)合實況資料研究了積層混合云的降水特征、云物理結(jié)構(gòu)特征和微物理過程。于翡和姚展予(2009)利用WRF-ARW中尺度數(shù)值模式對發(fā)生在華北地區(qū)一次積層混合云降水過程進行實例模擬,并結(jié)合實測資料,從動力學、熱力學和云微物理的角度研究積層混合云的發(fā)展機制。

    本文利用中尺度數(shù)值模式WRF的較新版本W(wǎng)RF-ARW(V3.2),對2009年4月18—19日發(fā)生在張家口地區(qū)的一次積層混合云降水個例進行實例模擬,并結(jié)合觀測資料,從天氣系統(tǒng)、中尺度云雨帶、云中微物理量等不同尺度分析了該次降水過程,有利于加深對積層混合云降水的認識,對播云催化位置選取、人工增雨作業(yè)具有一定的指導意義。

    1 資料與模擬方案介紹

    觀測資料主要有:地面雨量資料(1、6、24 h);FY-2C衛(wèi)星的觀測資料及其反演的產(chǎn)品;張家口站多普勒雷達(CB)的觀測資料;飛機觀測資料。模式的初始場采用2009年4月17—19日的NCEP再分析資料(FNL from GFS),空間分辨率為1°×1°,時間分辨率為6 h。

    圖1 模擬區(qū)域Fig.1 Simulation domains

    本文使用WRF的較新版本W(wǎng)RF-ARW(V3.2)進行模擬。模式采用三重雙向嵌套方案(圖1),分辨率分別為27、9、3 km,分別用于天氣系統(tǒng)、局地中尺度系統(tǒng)和對流系統(tǒng)的研究;網(wǎng)格中心位于張家口市(115°E,41°N),垂直方向采用σ坐標,共分27層。對于可分辨尺度降水,云微物理過程均采用Morrison雙參數(shù)方案,該方案包含6種水成物(雨水、云水、冰晶、霰(雹)、水汽、雪)以及多種液態(tài)、固態(tài)和混合態(tài)過程,是WRF中較為復雜的云微物理過程方案,可以很好地反映云、雨過程。對于不可分辨尺度降水,粗、中網(wǎng)格采用Grell-Devenyi積云對流參數(shù)化集合方案,細網(wǎng)格已能分辨一般的對流降水,故本文沒有采用積云參數(shù)化方案。邊界物理過程參數(shù)化使用了Yonsei University(YSU)方案。

    模擬時段為2009年4月17日20時—19日20時(北京時間),共48 h。粗網(wǎng)格與中網(wǎng)格每小時輸出一次結(jié)果,細網(wǎng)格每半小時輸出一次,以此來研究這次降水過程中云的結(jié)構(gòu)及其發(fā)展過程。

    圖2 NCEP(a,c,e)和模擬(b,d,f)的2009年4月18日08時500 hPa(a,b)、700 hPa(c,d)、850 hPa(e,f)位勢高度(單位:gpm)Fig.2 (a,c,e)NCEP and (b,d,f)simulated geopotential height at (a,b)500 hPa,(c,d)700 hPa,and (e,f)850 hPa at 08:00 BST 18 April 2009(units:gpm)

    2 天氣形勢對比分析

    2.1 天氣形勢變化情況

    圖2為2009年4月18日08時NCEP和模擬的高度場??梢钥吹?模式較好模擬了500 hPa高度上河北南部至河南東部一帶西北走向的槽,及700 hPa高度上內(nèi)蒙西部的橫槽和內(nèi)蒙中部到山西一帶的豎槽。由于地形原因,850 hPa(約1 500 m)部分地區(qū)出現(xiàn)了斷續(xù),可以看到,模式較好模擬了800 hPa低壓中心數(shù)值及分布,中心值均為1 400 gpm,低壓中心分布在111°E、43°N附近??偟膩碚f,模式模擬的位勢高度場與NCEP基本一致,模擬效果較好。

    對比分析2009年4月18日02時—19日08時各高度場的變化可知,模式對高度場演變的模擬較好,模擬結(jié)果與NCEP基本一致,可以很好地反映這次降水過程中天氣系統(tǒng)的演變情況。分析各高度場的變化(圖略)可知:在這次過程中,500 hPa上張家口附近較平滑,沒有明顯的天氣形勢變化;700 hPa上有西風槽經(jīng)過張家口地區(qū),18日02時槽線位于內(nèi)蒙古—陜西境內(nèi),呈南北向,02—14時槽線東移,14時槽線移至河北正前方,并與蒙古地區(qū)下來的低壓結(jié)合,槽線加深,此時張家口正處于槽前,14—20時槽線繼續(xù)東移并轉(zhuǎn)向成東北—西南走向,20時槽線位于張家口境內(nèi),20時后槽線繼續(xù)東移并減弱,次日08時槽線移出張家口地區(qū);850 hPa上18日02時蒙古地區(qū)有一低渦,02—14時低渦東移,14時低渦移到內(nèi)蒙古—河北地區(qū),14—20時低渦向東偏南移動,影響張家口地區(qū),20時后低渦消散。

    結(jié)合18日08時至19日08時張家口地區(qū)的小時加密雨量分布來看,從18日09時起張家口北部地區(qū)開始出現(xiàn)零星小雨,16時雨量開始逐漸增大并持續(xù)向南移動,18日18時至19日02時降雨量最大,19日03時雨量開始逐漸減小并漸止。由此可知,700 hPa西風槽、850 hPa低渦是影響這次張家口地區(qū)降水的主要天氣系統(tǒng)。

    2.2 水汽場分析

    為了對此次過程的水汽條件進行分析,模式模擬了18日08時—19日08時800 hPa水汽通量和風場(圖略)。通過對比可知,模擬的水汽通量大小、分布及演變與NCEP基本一致,模擬的風場與NCEP也基本一致,模擬效果較好。

    分析18日08時—19日08時水汽通量變化可知,張家口地區(qū)的水汽輸送主要有兩個:一個是來自于南方的暖濕空氣,在南風急流的作用下隨著時間推移源源不斷地向北輸送;另一個水汽輸送是來自于西北蒙古地區(qū),08—14時水汽隨著氣旋環(huán)流向東移動,14—20時向南移動,14—16時到達張家口地區(qū),受這兩股水汽輸送的影響,從16時起張家口地區(qū)地面降水明顯增加。配合800 hPa的高度場可知,14—20時張家口地區(qū)有一個低壓中心,即在該時段內(nèi),張家口地區(qū)低層出現(xiàn)了大尺度輻合,有利于云系發(fā)展和降水形成。

    針對14—20時張家口地區(qū)水汽輻合輻散情況做進一步分析。圖3為18日16—20時各層水汽通量散度分布??梢?在800 hPa上,18日16時張家口附近已形成一條明顯的水汽輻合帶,輻合帶呈東北—西南走向的帶狀分布,與衛(wèi)星觀測云帶位置對應較好,輻合帶長約700~800 km,寬約200~300 km,且隨時間向東向張家口地區(qū)移動,20時開始減弱,并隨時間逐漸縮小、消散。在800 hPa的水汽輻合帶位置,500 hPa上對應為輻散,輻散區(qū)域也呈東北—西南走向的帶狀分布。通過分析各層水汽通量散度可知,14—22時張家口附近地區(qū)出現(xiàn)了低層水汽輻合、高層水汽輻散的配置,有利于水汽的垂直輸送、云系的發(fā)展,進而產(chǎn)生降水。

    3 中尺度云、降水場對比分析

    3.1 云場對比分析

    把模式輸出的水成物(云水比含量、冰晶比含量)進行垂直積分得到云帶的分布,與FY-2C數(shù)據(jù)反演的光學厚度進行對比。由于FY-2C的反演受可見光通道限制,反演的產(chǎn)品只有在正午前后效果最佳,因此選取離降水時段較近且光學厚度較好的時刻進行比對。圖4是18日14時FY-2C反演光學厚度和模擬云帶的對比(其中模擬云帶采用了模式輸出的第二層數(shù)據(jù),分辨率為9 km)。比較發(fā)現(xiàn),FY-2C數(shù)據(jù)反演出的光學厚度與模擬云帶的相關(guān)性較好,尤其是大值區(qū)的分布對應較好。對比多個時次的結(jié)果可看到,模式模擬的云帶可以較好地反映此次過程中云帶的移動、變化,模式對中尺度系統(tǒng)的模擬效果較好。

    分析云帶變化情況(圖略)可知:10時內(nèi)蒙古中東部與河北北部的光學厚度較厚,且分布著幾個高值中心,中心值為30~36 m,山西南部和河南大部的光學厚度也較厚,最大值超過36 m;同時,在模擬云帶上,內(nèi)蒙古中東部及河北北部的水成物較豐富,且水成物總體分布不均勻,存在多個小的高值中心。14時云帶發(fā)展、東移,此時云帶呈東北—西南向帶狀分布于河北北面,帶長約1 000 km,帶寬約300 km,在大片的云帶中分布著很多個小的高值中心,中心區(qū)域一般為幾十千米,中心數(shù)值為3.0~4.0 g/m2。14—17時云帶向東南方向移動,17時移至河北西北部,且覆蓋了張家口部分地區(qū);17—22時云帶繼續(xù)東移、發(fā)展,22時幾乎覆蓋了整個張家口地區(qū),17—22時云帶中仍分布著很多小的高值中心,22時至次日04時云帶繼續(xù)東移,19日04時云帶高值中心移出張家口地區(qū),08時云帶移出張家口地區(qū),張家口上空幾乎無云帶覆蓋。

    圖3 模擬的2009年4月18日16—20時800 hPa(a,c,e)與500 hPa(b,d,f)水汽通量散度分布(單位:g/(s·cm2·hPa)) a,b.16時;c,d.18時;e,f.20時Fig.3 Simulated (a,c,e)800 hPa and (b,d,f)500 hPa moisture flux divergence from 16:00 BST to 20:00 BST 18 April 2009(units:g/(s·cm2·hPa)) a,b.16:00 BST;c,d.18:00 BST;e,f.20:00 BST

    圖4 2009年4月18日14時FY-2C反演光學厚度(a;單位:m)和模擬云帶(b;單位:g/m2)Fig.4 (a)Optical thickness(units:m) retrieved from FY-2C and (b)simulated cloud band(units:g/m2) at 14:00 BST 18 April 2009

    3.2 地面降水對比分析

    圖5為觀測與模擬的24 h降水量分布(其中圖5b為模式第二層的輸出結(jié)果)。通過對比可知:模擬的24 h降水量分布與觀測結(jié)果基本一致,尤其是5 mm以上的降水分布位置、走向與觀測結(jié)果非常接近;降水中心的分布及強度的模擬效果也較好;從雨帶的分布與強降水落區(qū)來看,模式具有較好的模擬能力。

    分析各時段降水分布(圖略)可知:18日08—14時6 h降水呈東北—西南向,帶狀分布在河北西北地區(qū),且在較大的雨帶中分布著三個小的高值中心,中心降水量為5~10 mm,張家口地區(qū)降水量較小。14—20時降水量增大,幾乎覆蓋整個河北地區(qū),且降水量中心擴大南移,模式模擬的降水較連續(xù),此時幾個中心區(qū)連成一片,形成帶狀。18日20時—19日08時整個華北地區(qū)普遍降水,降水高值分布在河南南部地區(qū),張家口地區(qū)降水量較14—20時有所減少,其降水中心也移至張家口南部。對比6、12、24 h降水量可知,張家口地區(qū)降水的模擬值比實測值大,且在較大的雨帶中分布著若干個小的高值中心,中心區(qū)域水平范圍在幾十千米,中心降水量一般是大范圍降水量的3~5倍,降水總體分布不均勻。

    圖5 觀測(a)和模擬(b)的2009年4月18日08時—19日08時24 h降水量分布(單位:mm)Fig.5 (a)Observed and (b)simulated 24 h rainfall distributions from 08:00 BST 18 to 08:00 BST 19 April 2009(units:mm)

    圖6 觀測(a)和模擬(b)的2009年4月18日17時雷達回波反射率(單位:dBz)Fig.6 (a)Observed and (b)simulated radar echo reflectivity at 17:00 BST 18April 2009(units:dBz)

    3.3 雷達回波對比分析

    圖6為觀測和模擬的18日17時雷達回波反射率。圖6a是張家口站多普勒雷達(CB)的雷達組合反射率因子圖,可以看到雷達回波整體邊緣呈支離破碎狀,強度小于35 dBz的為層狀云降水回波,在較均勻的層狀云回波中鑲嵌著幾個柱狀對流回波(已標出),回波值超過50 dBz,具有典型的積層混合云降水回波特征。圖6b為由模式輸出的第三層數(shù)據(jù)繪制的雷達回波反射率圖,分辨率為3 km,用以研究小尺度系統(tǒng)。為了更直觀地與觀測結(jié)果進行對比,根據(jù)張家口站中心位置及雷達覆蓋范圍,在圖6b中選取大致相同的區(qū)域(紅色虛線內(nèi))。

    通過比較可知,模擬與實測的雷達回波反射率分布位置、形狀及最大值較一致,都是在較均勻的回波層中鑲嵌著幾個柱狀的對流回波。對流回波的位置在圖6中已標出,模擬位置與實測位置稍偏一點;對流回波呈柱狀分布,水平范圍一般在10~20 km;回波中心值超過50 dBz??傮w而言,模式對小尺度系統(tǒng)的模擬效果較好。

    4 小尺度微物理分析

    4.1 云中水成物分布特征

    圖7是模擬的18日17時雷達回波圖及沿剖線AB的各物理量的剖面圖。圖7a給出了模擬的雷達回波平面圖及剖線AB的位置,方向是由A指向B,且穿過兩個對流云區(qū)。圖7b是沿直線AB的雷達回波剖面圖,可以看到沿著剖線的對流回波較強,最大值超過50 dBz,回波頂高在11 km左右。圖7c是云水的垂直分布,可以看到無論是沿著水平方向還是垂直方向都有多個含水量高值中心,中心區(qū)域在10 km左右,最大值達到0.5 g/kg,云水含量總體分布不均勻;同時可以看到云水的大值中心分布在零度層以上,即云中存在較豐沛的過冷云水。圖7d是冰晶的垂直分布,可以看到冰晶水平分布范圍較廣,存在多個高值中心,最大值達到0.1 g/kg,垂直分布較窄,主要位于400~200 hPa。圖7e—g分別是霰、雪和雨水的垂直分布,可以看到霰、雪和下層的雨水分布是相對應的,霰、雪的大值中心對應著下層雨水的大值中心,水平分布上都存在兩個高值中心,水平范圍為10~20 km,其位置與圖7b中雷達回波相對應,可見霰、雪和雨水是主要的降水粒子,是形成雷達回波的主要影響因子。圖7h是總水凝物含量剖面圖,可以看到沿直線AB上總的水凝物分布很廣,水平方向存在兩個大值中心,最大值為5.0 g/kg。

    4.2 積層混合云降水討論

    為了進一步研究這次降水過程,沿著云帶移動方向做剖面分析。圖8是沿云帶移動方向做的一組剖面圖。圖8a給出了模擬的雷達回波平面圖及剖線的位置,剖線平行于云帶的移動方向;選取剖線穿過的一個對流云區(qū)(圖中虛線所示),可以看到該對流回波呈橢圓狀,水平尺度為10~20 km,內(nèi)嵌在大片的層狀云回波中。圖8b是沿著剖線的垂直速度分布,相應的對流云區(qū)在圖中用虛線標示,可以看到17時在對流云區(qū)內(nèi)以上升氣流為主,最大上升速度超過2.5 m/s,但在對流云區(qū)的前上部和后下方存在與上升速度同量級的下沉氣流,表明此時的對流云正處于成熟階段。圖8c是沿著剖線的云水、冰晶數(shù)濃度分布,可以看到云系前部上層存在大量的冰晶,但沒有云水,推測云系前部主要由高卷云、高層云組成,隨著經(jīng)緯度變化云系開始出現(xiàn)云水、霰和雪,云系開始出現(xiàn)層云,繼續(xù)西移后云中開始出現(xiàn)大量的霰、雪和雨水,此時以積雨云、高積云為主,在云系后部又出現(xiàn)大量冰晶但沒有云水,也是以卷云為主。圖8d是沿剖線的霰、雪和雨水的分布。

    結(jié)合圖8b、c、d可以看到,在對流云區(qū)(虛線)內(nèi),霰、雪比含量很大,并與下層的雨水分布相對應,而云水比含量并不大,尤其是在霰、雪比含量的大值區(qū)幾乎沒有云水,存在較典型的冷云降水機制;霰、雪比含量大值區(qū)的上升速度較強,最大值超過2.5 m/s,霰、雪在上升過程中不斷消耗過冷云水,生成較大的降水粒子,大降水粒子下落,落到零度層以下融化成雨水,雨水降落至地面形成降水,因此在對流云區(qū)的后部零度層以下由于降水粒子的拖曳出現(xiàn)了下沉氣流。結(jié)合圖8c、d可知,在層狀云區(qū)內(nèi),選定區(qū)域(點線內(nèi))的下層有雨水分布,而與之對應的上層卻沒有霰、雪存在,只有一些過冷云水分布,該區(qū)域的雨水主要是由過冷云水直接轉(zhuǎn)化而成的,存在著暖云降水過程。

    綜合所述,此次降水過程中的積層混合云不僅有積狀云和層狀云形態(tài)的混合,在微物理方面還存在冷云過程和暖云過程相態(tài)的混合,降水機制較復雜。本研究只針對了一次過程,選取多個個例對積層混合云進行深入研究,對播云催化位置選取、人工增雨作業(yè)具有一定的指導意義。

    圖7 模擬的17:00雷達回波平面圖及沿剖線AB各物理量的剖面(圖中陰影為雷達回波,單位:dBz;實線為各物理量分布,單位:g/kg;虛線為溫度,單位:℃;箭矢為風場,單位:m/s) a.模擬雷達回波平面圖及剖線AB;b.沿剖線AB的雷達回波剖面圖;c.云水比含量;d.冰晶比含量;e.霰比含水量;f.雨比含量;g.雪比含量;h.總水凝物Fig.7 Horizontal distribution of simulated radar echo and vertical distribution of simulated hydrometeors along line AB at 17:00 BST(shadings are radar echo(units:dBz);solid lines denote distribution of hydrometeors(units:k/kg);dashed lines show temperature(units:℃);arrows are wind field(units:m/s)) a.horizontal distribution of simulated radar echo and position of line AB;b.vertical distribution of simulated radar echo along line AB;c.cloud water mixing ratio;d.ice mixing ratio;e.graupel mixing ratio;f.rain mixing ratio;g.snow mixing ratio;h.total hydrometeors

    圖8 模擬的17:00雷達回波平面圖及沿剖線的各物理量剖面 a.模擬的雷達回波(單位:dBz)平面圖及剖線CD;b.垂直速度(陰影,單位:m/s)分布;c.云水比含量(陰影,單位:g/kg)、冰晶數(shù)濃度(黑色實線,單位:kg-1)、溫度(紅色虛線,單位:℃)分布;d.霰、雪和比含量(陰影,單位:g/kg)、雨水比含量(黑色實線,單位:g/kg)、溫度(紅色虛線,單位:℃)分布Fig.8 Horizontal distribution of simulated radar echo and vertical distribution of simulated hydrometeors at 17:00 BST a.horizontal distribution of simulated radar echo(units:dBz) and position of line CD;b.vertical velocity profile(shadings;units:m/s);c.cloud water mixing ratio(shadings;units:g/kg) and ice number concentration(black solid lines;units:kg-1) and temperature(red dashed lines;units:℃);d.the sum of graupel and snow mixing ratio(shadings;units:g/kg),rain mixing ratio(black solid lines;units:g/kg) and temperature(red dashed lines;units:℃)

    5 結(jié)論

    1)比較分析實測和模擬資料可知:700 hPa西風槽、850 hPa低渦是形成這次降水的主要天氣系統(tǒng),來自南方的暖濕空氣和西北內(nèi)蒙古低渦帶來的水汽是這次降水的主要水汽來源,14—20時兩股水汽在張家口附近出現(xiàn)低層輻合、中高空輻散的配置,有利于水汽的垂直輸送、云系的發(fā)展。

    2)比較分析中尺度系統(tǒng)可知:張家口附近云系呈東北—西南向帶狀分布,帶長約1 000 km,帶寬300 km,在大片的云帶中分布著很多個小的高值中心,中心區(qū)域一般在幾十千米;降水也呈東北—西南向帶狀分布,在較大的雨帶中分布著幾個小的雨核,降水總體分布不均勻;結(jié)合雷達回波可見整體回波邊緣呈支離破碎狀,在較均勻的層狀云回波中鑲嵌著柱狀的對流回波,具有典型的積層混合云降水回波特征。

    3)比較分析小尺度系統(tǒng)可知:對流云波呈柱狀分布,水平范圍在10~20 km,回波中心值超過50 dBz;沿對流回波做剖面發(fā)現(xiàn),云系中云水含量在水平方向和垂直方向上都分布不均勻,雨水的大值中心與上層的霰、雪的大值中心相對應,其位置與剖線處雷達回波相對應,表明霰、雪和雨水是主要的降水粒子,是形成雷達回波的主要影響因子。此次積層混合云降水不僅有積狀云和層狀云形態(tài)的混合,在微物理方面還存在著冷云過程和暖云過程相態(tài)的混合,降水機制較復雜。

    4)WRF模式對此次降水過程的模擬效果較好。對大尺度高度場、水汽場,中尺度云場、降水場,小尺度對流回波及云中各水凝物的分布模擬效果較好。模式對實況的模擬在時間上沒有滯后,落區(qū)也較準確,具有較好的模擬能力。

    宮福久,周德平,陳寶君,等.2005.東北冷渦云物理及導變技術(shù)[M].北京:氣象出版社:11.

    洪延超.1996a.積層混合云數(shù)值模擬研究(Ⅰ)——模式及其微物理過程參數(shù)化[J].氣象學報,54(5):544-557.

    洪延超.1996b.積層混合云數(shù)值模擬研究(Ⅱ)——云相互作用及暴雨產(chǎn)生機制[J].氣象學報,54(6):661-674.

    劉黎平,阮征,覃丹宇.2004.長江流域梅雨鋒暴雨過程的中尺度結(jié)構(gòu)個例分析[J].中國科學D輯,34(12):1193-1201.

    王維佳,陶遐齡.2006.一維半地形積層混合云模式簡介[J].四川氣象,26(1):43-45.

    王維佳,陶遐齡.2007.地形積層混合云模式數(shù)值試驗研究[J].四川氣象,27(1):29-32.

    許梓秀,王鵬云.1989.冷鋒前部中尺度雨帶特征及其機制分析[J].氣象學報,47(2):198-206.

    于翡,姚展予.2009.一次積層混合云降水實例的數(shù)值模擬分析[J].氣象,35(12):3-11.

    鄒倩,劉琦俊,劉衛(wèi)國.2008.積層混合云結(jié)構(gòu)和云微物理的數(shù)值模擬[J].干旱氣象,26(2):74-83.

    Anagnostou E N.2004.A convective/stratiform precipitation classification algorithm for volume scanning weather radar observations[J].Meteorological Applications,11(4):291-300.

    Fabry 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.

    Fuhrer O,Schaer C.2002.Effects of moist convection on flow past 2D-topography (2002-10MntMet)[C]//Proceedings of the 10th AMS Conference on Mountain Meteorology:17-21.

    Fuhrer O,Schaer C.2005.Embedded cellular convection in moist flow past topography[J].J Atmos Sci,62(8):2810-2828.

    Orville H D,Hirsch J H,Farley R D.1987.Further results on numerical cloud seeding simulations of stratiform-type clouds[J].The Journal of Weather Modification,19(1):57-61.

    Yuter S E,Houze R A,Smith E A,et al.2005.Physical characterization of tropical oceanic convection observed in KWAJEX[J].J Appl Meteor,44(4):385-415.

    (責任編輯:倪東鴻)

    Numericalsimulationonaconvective-stratiformmixedcloudprecipitationindifferentscales

    ZHANG Wei1,ZHOU Yu-quan2

    (1.Nanjing University of Information Science and Technology,Nanjing 210044,China;2.Chinese Academy of Meteorological Sciences,Beijing 100081,China)

    The precipitation of a convective-stratiform mixed cloud in Zhangjiakou during 18—19 April 2009 is simulated using the mesoscale numerical model WRF-ARW(V3.2) and comparatively analyzed with the observation data in different scales.Results show that the west wind trough at 700 hPa and the low vortex at 850 hPa are the main weather systems responsible for convective-stratiform mixed cloud.The warm wet air from the southern and the low vortex from northwest Inner Mongolia are main water vapor suppliers,which converge in Zhangjiakou region and thus are conducive to cloud system development and precipitation formation.The cloud system shows a NE-SW banded distribution,which is 1 000 km long and 300 km wide,and has lots of cloud water centers,which are about tens of kilometers.The cloud system has the echo characteristics of typical convective-stratiform mixed cloud precipitation,such as columnar echo of cumulus cloud embedded in uniform echo layer.The cross-section of radar echo shows that the cloud water content is heterogeneous either horizontally or vertically.The centers of rainwater are corresponding with the centers of graupel and snow in the upper layer,with horizontal range of 10—20 km.

    WRF model;convective-stratiform mixed cloud;structure in different scales;numerical simulation

    2012-03-27;改回日期2014-04-22

    公益性行業(yè)(氣象)科研專項(GYHY201206025)

    張微,碩士,研究方向為云降水物理和人工影響天氣,zhangwei198633@126.com.

    10.13878/j.cnki.dqkxxb.20120327002.

    1674-7097(2014)04-0459-10

    P426.51;P426.611

    A

    10.13878/j.cnki.dqkxxb.20120327002

    張微,周毓荃.2014.一次積層混合云降水不同尺度結(jié)構(gòu)的數(shù)值模擬[J].大氣科學學報,37(4):459-468.

    Zhang Wei,Zhou Yu-quan.2014.Numerical simulation on a convective-stratiform mixed cloud precipitation in different scales[J].Trans Atmos Sci,37(4):459-468.(in Chinese)

    18禁观看日本| 久久久精品欧美日韩精品| 亚洲国产欧洲综合997久久,| 成年女人永久免费观看视频| 欧美性猛交黑人性爽| 最近最新免费中文字幕在线| 首页视频小说图片口味搜索| 99国产精品99久久久久| 国产极品精品免费视频能看的| 国产高清激情床上av| 日本a在线网址| 在线看三级毛片| 看片在线看免费视频| 国产精品美女特级片免费视频播放器 | 免费看a级黄色片| 巨乳人妻的诱惑在线观看| 国产精品自产拍在线观看55亚洲| 一进一出抽搐动态| 免费观看精品视频网站| 一区二区三区激情视频| 草草在线视频免费看| 亚洲av第一区精品v没综合| 美女高潮的动态| 国产激情偷乱视频一区二区| 999精品在线视频| 草草在线视频免费看| 日韩三级视频一区二区三区| 国产三级黄色录像| 国产成人福利小说| 欧美黑人巨大hd| 99久国产av精品| 国产免费男女视频| 91av网一区二区| 亚洲九九香蕉| 最新在线观看一区二区三区| 久久精品国产综合久久久| 久久精品人妻少妇| 床上黄色一级片| 国产欧美日韩一区二区精品| 天堂网av新在线| 国产亚洲精品久久久久久毛片| 黄色女人牲交| 18美女黄网站色大片免费观看| 91麻豆av在线| 久久久久国产精品人妻aⅴ院| 亚洲自拍偷在线| 国产精品精品国产色婷婷| 国内精品久久久久久久电影| 亚洲av成人一区二区三| 91字幕亚洲| 色综合站精品国产| 99久久国产精品久久久| 成人精品一区二区免费| 日本五十路高清| 露出奶头的视频| 99久久国产精品久久久| 91在线观看av| 精品久久久久久,| 亚洲九九香蕉| 国产v大片淫在线免费观看| 超碰成人久久| 男插女下体视频免费在线播放| xxxwww97欧美| 欧美高清成人免费视频www| 国产爱豆传媒在线观看| 国产三级黄色录像| 久久这里只有精品19| 麻豆一二三区av精品| 999久久久国产精品视频| 1024香蕉在线观看| 亚洲精品美女久久久久99蜜臀| 日本在线视频免费播放| 一级毛片女人18水好多| 欧美极品一区二区三区四区| 国产 一区 欧美 日韩| 可以在线观看毛片的网站| 亚洲一区二区三区不卡视频| 又紧又爽又黄一区二区| 婷婷精品国产亚洲av在线| 久久精品aⅴ一区二区三区四区| 日韩欧美免费精品| 在线观看日韩欧美| 操出白浆在线播放| 国产精品久久电影中文字幕| 亚洲欧美精品综合一区二区三区| 国产私拍福利视频在线观看| 欧美成狂野欧美在线观看| 精品一区二区三区av网在线观看| 一级黄色大片毛片| 亚洲自拍偷在线| 色综合亚洲欧美另类图片| 午夜福利免费观看在线| 美女cb高潮喷水在线观看 | 国产伦精品一区二区三区视频9 | 国产真实乱freesex| 美女被艹到高潮喷水动态| 国内少妇人妻偷人精品xxx网站 | 久久这里只有精品19| 日本五十路高清| 18禁黄网站禁片免费观看直播| 国产精品,欧美在线| 久久久久久国产a免费观看| 国产成人精品久久二区二区91| 午夜福利在线观看吧| 夜夜夜夜夜久久久久| 欧美极品一区二区三区四区| 一级毛片女人18水好多| 亚洲专区国产一区二区| 我的老师免费观看完整版| 亚洲 国产 在线| 中文字幕人妻丝袜一区二区| 免费人成视频x8x8入口观看| 精品一区二区三区四区五区乱码| 美女扒开内裤让男人捅视频| 国产激情偷乱视频一区二区| 国内精品久久久久精免费| 国产视频一区二区在线看| 伊人久久大香线蕉亚洲五| 国产精品九九99| 母亲3免费完整高清在线观看| 国产成人一区二区三区免费视频网站| 国产激情久久老熟女| 成人永久免费在线观看视频| 嫩草影视91久久| 国产午夜精品久久久久久| 国模一区二区三区四区视频 | 午夜视频精品福利| 成人特级av手机在线观看| 一进一出抽搐动态| 亚洲avbb在线观看| 亚洲欧美日韩东京热| 亚洲中文日韩欧美视频| 国产精品国产高清国产av| 黄色 视频免费看| 亚洲成人久久性| 男人舔女人下体高潮全视频| 老鸭窝网址在线观看| 国产精品永久免费网站| 丁香六月欧美| 少妇熟女aⅴ在线视频| 久久国产精品影院| 免费看a级黄色片| 一个人免费在线观看电影 | 亚洲精品在线美女| 亚洲18禁久久av| 日韩有码中文字幕| 精品久久久久久成人av| 波多野结衣巨乳人妻| 日本黄大片高清| 老汉色av国产亚洲站长工具| 亚洲国产欧美一区二区综合| 久久这里只有精品19| 日日夜夜操网爽| 两个人的视频大全免费| 19禁男女啪啪无遮挡网站| a在线观看视频网站| 黄色丝袜av网址大全| 精品国产亚洲在线| 免费高清视频大片| 久99久视频精品免费| 中文在线观看免费www的网站| 看免费av毛片| 欧美日韩综合久久久久久 | 91麻豆精品激情在线观看国产| av在线天堂中文字幕| 麻豆国产av国片精品| netflix在线观看网站| 久久精品夜夜夜夜夜久久蜜豆| 香蕉国产在线看| 日日摸夜夜添夜夜添小说| 麻豆久久精品国产亚洲av| 亚洲欧美激情综合另类| 亚洲人与动物交配视频| 久久久精品大字幕| 一级毛片女人18水好多| 国产综合懂色| 黑人巨大精品欧美一区二区mp4| 国产精品,欧美在线| 十八禁人妻一区二区| 久久午夜综合久久蜜桃| 婷婷精品国产亚洲av| 五月伊人婷婷丁香| 成人午夜高清在线视频| 又爽又黄无遮挡网站| 久久午夜综合久久蜜桃| 高清毛片免费观看视频网站| 老司机福利观看| 白带黄色成豆腐渣| 男女做爰动态图高潮gif福利片| 免费搜索国产男女视频| 99久久国产精品久久久| 99精品在免费线老司机午夜| ponron亚洲| 国产精品日韩av在线免费观看| 亚洲av第一区精品v没综合| 欧美黄色淫秽网站| 在线观看舔阴道视频| 色在线成人网| 18禁美女被吸乳视频| 欧美一级毛片孕妇| 久久精品亚洲精品国产色婷小说| 又粗又爽又猛毛片免费看| 每晚都被弄得嗷嗷叫到高潮| 久久久久国产精品人妻aⅴ院| 午夜福利18| 欧美性猛交黑人性爽| 国模一区二区三区四区视频 | 天天一区二区日本电影三级| 午夜福利在线观看免费完整高清在 | 国产精品女同一区二区软件 | www.熟女人妻精品国产| 成人高潮视频无遮挡免费网站| 999精品在线视频| 在线十欧美十亚洲十日本专区| 热99在线观看视频| 嫁个100分男人电影在线观看| 成人欧美大片| 亚洲人成网站在线播放欧美日韩| 变态另类丝袜制服| 女警被强在线播放| 美女被艹到高潮喷水动态| 大型黄色视频在线免费观看| 亚洲av日韩精品久久久久久密| bbb黄色大片| 小蜜桃在线观看免费完整版高清| 一区二区三区高清视频在线| 十八禁网站免费在线| 久99久视频精品免费| 国产成人精品久久二区二区免费| 波多野结衣高清无吗| 亚洲性夜色夜夜综合| 国产亚洲欧美98| 搞女人的毛片| avwww免费| 免费人成视频x8x8入口观看| 丰满人妻熟妇乱又伦精品不卡| av天堂中文字幕网| 美女高潮喷水抽搐中文字幕| 日韩有码中文字幕| 99久久99久久久精品蜜桃| 日本与韩国留学比较| 国产精品精品国产色婷婷| 精品久久蜜臀av无| 成人高潮视频无遮挡免费网站| 亚洲电影在线观看av| 18禁国产床啪视频网站| 国产精品九九99| 香蕉国产在线看| 免费一级毛片在线播放高清视频| 日韩欧美国产一区二区入口| av在线天堂中文字幕| 特大巨黑吊av在线直播| 三级男女做爰猛烈吃奶摸视频| 国产精品免费一区二区三区在线| 久久九九热精品免费| 一区二区三区国产精品乱码| 日韩欧美精品v在线| 亚洲国产精品成人综合色| 老司机午夜福利在线观看视频| 中文资源天堂在线| 日韩 欧美 亚洲 中文字幕| 亚洲国产欧美一区二区综合| 最近最新中文字幕大全电影3| 每晚都被弄得嗷嗷叫到高潮| 99国产综合亚洲精品| 欧美一级毛片孕妇| 成人高潮视频无遮挡免费网站| 美女午夜性视频免费| 99国产精品一区二区蜜桃av| 老司机午夜十八禁免费视频| 一级毛片高清免费大全| 国产高潮美女av| 在线免费观看的www视频| 悠悠久久av| 麻豆国产97在线/欧美| 黄色片一级片一级黄色片| 国产激情偷乱视频一区二区| 婷婷精品国产亚洲av在线| 国产精品精品国产色婷婷| 无人区码免费观看不卡| 国产伦精品一区二区三区四那| 国产精品香港三级国产av潘金莲| 女人被狂操c到高潮| 亚洲成人久久性| 这个男人来自地球电影免费观看| 欧美不卡视频在线免费观看| 国产伦人伦偷精品视频| 国产熟女xx| 一二三四社区在线视频社区8| 亚洲第一电影网av| 国产精品影院久久| 亚洲七黄色美女视频| 99久久精品国产亚洲精品| 久久久国产精品麻豆| 在线观看免费午夜福利视频| 中文在线观看免费www的网站| 小说图片视频综合网站| 老汉色av国产亚洲站长工具| 亚洲真实伦在线观看| 男女床上黄色一级片免费看| 99热6这里只有精品| 日韩精品中文字幕看吧| 宅男免费午夜| 国产成人影院久久av| 99国产精品99久久久久| h日本视频在线播放| 一个人看视频在线观看www免费 | 成人永久免费在线观看视频| 最好的美女福利视频网| а√天堂www在线а√下载| 日本一本二区三区精品| 午夜免费成人在线视频| 国产一区二区在线av高清观看| 日韩欧美一区二区三区在线观看| 亚洲九九香蕉| 成人18禁在线播放| 又紧又爽又黄一区二区| 草草在线视频免费看| 国产真实乱freesex| 欧美乱色亚洲激情| 亚洲精品粉嫩美女一区| 日日夜夜操网爽| 五月玫瑰六月丁香| 久久国产精品人妻蜜桃| 久久天堂一区二区三区四区| 成人18禁在线播放| 国产免费av片在线观看野外av| 一卡2卡三卡四卡精品乱码亚洲| 天天躁狠狠躁夜夜躁狠狠躁| www国产在线视频色| 精品人妻1区二区| 狂野欧美白嫩少妇大欣赏| 日韩欧美 国产精品| 久久午夜亚洲精品久久| 亚洲精品在线美女| 久久午夜亚洲精品久久| 老汉色∧v一级毛片| 青草久久国产| 身体一侧抽搐| 国产真实乱freesex| 色哟哟哟哟哟哟| 国产亚洲av嫩草精品影院| 99视频精品全部免费 在线 | 亚洲专区国产一区二区| 欧美成人一区二区免费高清观看 | 一级作爱视频免费观看| ponron亚洲| 一级作爱视频免费观看| 大型黄色视频在线免费观看| 精品无人区乱码1区二区| 久久精品人妻少妇| 免费在线观看成人毛片| 男人舔奶头视频| 热99re8久久精品国产| 日韩人妻高清精品专区| 黄色片一级片一级黄色片| 亚洲专区中文字幕在线| 99在线视频只有这里精品首页| 国产一区二区三区视频了| tocl精华| 国产精品综合久久久久久久免费| 国产成人av教育| 午夜福利高清视频| 亚洲欧美激情综合另类| 亚洲美女黄片视频| 少妇的丰满在线观看| 国产精品久久久av美女十八| 欧美成人免费av一区二区三区| 精华霜和精华液先用哪个| 久久精品国产99精品国产亚洲性色| 精品久久久久久,| 欧美中文日本在线观看视频| 最好的美女福利视频网| 国产淫片久久久久久久久 | 99久久综合精品五月天人人| 成人特级黄色片久久久久久久| 久久久久久大精品| 在线播放国产精品三级| 2021天堂中文幕一二区在线观| 国产欧美日韩一区二区精品| 国产精品亚洲一级av第二区| а√天堂www在线а√下载| 国产综合懂色| 日韩成人在线观看一区二区三区| 成年人黄色毛片网站| 999久久久精品免费观看国产| 久久午夜亚洲精品久久| 男人舔女人的私密视频| 久久久久性生活片| 一个人看视频在线观看www免费 | 久久精品亚洲精品国产色婷小说| 88av欧美| 日韩欧美国产在线观看| 88av欧美| 99国产精品99久久久久| 免费观看人在逋| 亚洲va日本ⅴa欧美va伊人久久| 亚洲欧美日韩卡通动漫| 精品不卡国产一区二区三区| 欧美成狂野欧美在线观看| 黄色 视频免费看| 国产成人aa在线观看| 神马国产精品三级电影在线观看| 国产单亲对白刺激| 最新中文字幕久久久久 | 亚洲欧美日韩东京热| 国产一区二区在线av高清观看| 欧美极品一区二区三区四区| 男人舔奶头视频| 免费观看的影片在线观看| 亚洲午夜精品一区,二区,三区| 亚洲av成人av| 亚洲av电影不卡..在线观看| 欧美日韩黄片免| 一个人免费在线观看电影 | 母亲3免费完整高清在线观看| 在线a可以看的网站| 日韩高清综合在线| 每晚都被弄得嗷嗷叫到高潮| e午夜精品久久久久久久| 亚洲成人中文字幕在线播放| 欧美国产日韩亚洲一区| 99热只有精品国产| 日韩大尺度精品在线看网址| a级毛片a级免费在线| 中文资源天堂在线| 免费人成视频x8x8入口观看| 久久精品91蜜桃| 99在线视频只有这里精品首页| 嫩草影视91久久| 又爽又黄无遮挡网站| 久久精品夜夜夜夜夜久久蜜豆| 午夜久久久久精精品| 免费高清视频大片| 香蕉久久夜色| 免费搜索国产男女视频| 不卡一级毛片| av片东京热男人的天堂| 国产一区在线观看成人免费| 国产综合懂色| 久久香蕉国产精品| 国产伦在线观看视频一区| 18禁美女被吸乳视频| 亚洲欧美精品综合一区二区三区| 最近在线观看免费完整版| 伊人久久大香线蕉亚洲五| 性欧美人与动物交配| 国产又色又爽无遮挡免费看| 中文在线观看免费www的网站| АⅤ资源中文在线天堂| 欧美一级a爱片免费观看看| 亚洲av免费在线观看| 97人妻精品一区二区三区麻豆| 非洲黑人性xxxx精品又粗又长| 中文字幕人成人乱码亚洲影| 国产高清videossex| 国产欧美日韩精品亚洲av| 日韩欧美在线乱码| 啪啪无遮挡十八禁网站| 久久精品91蜜桃| 午夜激情欧美在线| 国产精品香港三级国产av潘金莲| 91老司机精品| 国产精品九九99| 国产成+人综合+亚洲专区| 91久久精品国产一区二区成人 | 亚洲精品一卡2卡三卡4卡5卡| 一本一本综合久久| 最好的美女福利视频网| 国内久久婷婷六月综合欲色啪| 亚洲精品色激情综合| 日韩欧美在线乱码| 亚洲aⅴ乱码一区二区在线播放| 日韩中文字幕欧美一区二区| 国内精品久久久久久久电影| 草草在线视频免费看| 观看免费一级毛片| 舔av片在线| 精品福利观看| 欧美日韩福利视频一区二区| 日韩成人在线观看一区二区三区| 亚洲天堂国产精品一区在线| 亚洲av日韩精品久久久久久密| 天天一区二区日本电影三级| 青草久久国产| 日本与韩国留学比较| 国内久久婷婷六月综合欲色啪| 日韩三级视频一区二区三区| 亚洲欧美日韩东京热| 亚洲五月婷婷丁香| 男女午夜视频在线观看| 中文字幕最新亚洲高清| 亚洲18禁久久av| 久久久久久久久免费视频了| 国产成人系列免费观看| 精品国产美女av久久久久小说| 亚洲一区高清亚洲精品| 又粗又爽又猛毛片免费看| 成人性生交大片免费视频hd| 国产亚洲精品久久久久久毛片| 性色avwww在线观看| 757午夜福利合集在线观看| 久久久国产精品麻豆| 9191精品国产免费久久| 女生性感内裤真人,穿戴方法视频| 观看免费一级毛片| 欧美午夜高清在线| 精品久久久久久成人av| 亚洲美女黄片视频| 97超级碰碰碰精品色视频在线观看| www.自偷自拍.com| 美女 人体艺术 gogo| 波多野结衣高清无吗| x7x7x7水蜜桃| 男人舔奶头视频| 国产亚洲欧美98| 香蕉av资源在线| 亚洲中文日韩欧美视频| 亚洲精华国产精华精| 久久精品人妻少妇| 久久伊人香网站| 亚洲无线在线观看| 久久国产精品影院| 老司机午夜十八禁免费视频| 91麻豆精品激情在线观看国产| 在线国产一区二区在线| 又大又爽又粗| 欧美又色又爽又黄视频| 国产人伦9x9x在线观看| 欧美性猛交╳xxx乱大交人| 国产亚洲av高清不卡| 日本撒尿小便嘘嘘汇集6| 国产蜜桃级精品一区二区三区| 小蜜桃在线观看免费完整版高清| 色精品久久人妻99蜜桃| 欧美+亚洲+日韩+国产| 日本在线视频免费播放| 中文字幕熟女人妻在线| 亚洲成人免费电影在线观看| 免费在线观看成人毛片| 免费电影在线观看免费观看| 亚洲一区二区三区不卡视频| 亚洲成av人片免费观看| 精华霜和精华液先用哪个| 国产探花在线观看一区二区| 91久久精品国产一区二区成人 | 免费在线观看影片大全网站| 91麻豆av在线| 国产成人福利小说| 一个人看的www免费观看视频| 热99re8久久精品国产| 丰满的人妻完整版| 国产午夜精品论理片| 久久九九热精品免费| 97人妻精品一区二区三区麻豆| 啦啦啦观看免费观看视频高清| 两个人视频免费观看高清| 欧美一区二区精品小视频在线| 欧美极品一区二区三区四区| 久久久久久久午夜电影| 女同久久另类99精品国产91| 国产精品亚洲av一区麻豆| 亚洲熟妇中文字幕五十中出| 曰老女人黄片| 免费在线观看成人毛片| 美女高潮喷水抽搐中文字幕| 国内精品一区二区在线观看| 99精品久久久久人妻精品| 草草在线视频免费看| 欧美日本视频| 亚洲色图 男人天堂 中文字幕| 五月伊人婷婷丁香| 天堂av国产一区二区熟女人妻| 久久精品91无色码中文字幕| 亚洲美女黄片视频| 最近最新中文字幕大全电影3| 一二三四社区在线视频社区8| 国产激情久久老熟女| netflix在线观看网站| 九色成人免费人妻av| 亚洲欧美日韩高清在线视频| 成年女人永久免费观看视频| 不卡av一区二区三区| 久久久久国内视频| 一区二区三区高清视频在线| 狂野欧美白嫩少妇大欣赏| 此物有八面人人有两片| aaaaa片日本免费| 丁香欧美五月| 亚洲男人的天堂狠狠| 人妻夜夜爽99麻豆av| 三级毛片av免费| av天堂在线播放| 欧美另类亚洲清纯唯美| 精品国产美女av久久久久小说| 亚洲av成人精品一区久久| 国产91精品成人一区二区三区| 最新美女视频免费是黄的| 国产av不卡久久| 欧美日韩国产亚洲二区| 欧美成人性av电影在线观看| 国产精品久久久久久久电影 | 天堂√8在线中文| 窝窝影院91人妻| 亚洲狠狠婷婷综合久久图片| 在线观看舔阴道视频| 亚洲欧美日韩卡通动漫| 欧美不卡视频在线免费观看| 久久久久国内视频| 麻豆久久精品国产亚洲av| 成年免费大片在线观看| 亚洲片人在线观看| 欧美最黄视频在线播放免费| 少妇裸体淫交视频免费看高清| 成年版毛片免费区|