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

    青藏高原東部玉樹地區(qū)雨滴譜特征及云微物理結(jié)構(gòu)分析

    2021-10-07 05:21:04舒磊李茂善華生索朗江措呂釗伏薇王靈芝陰蜀城
    關(guān)鍵詞:積雨云玉樹層狀

    舒磊 李茂善 華生 索朗江措 呂釗 伏薇 王靈芝 陰蜀城

    (1 成都信息工程大學(xué)大氣科學(xué)學(xué)院/高原大氣與環(huán)境四川省重點實驗室/氣候環(huán)境變化聯(lián)合實驗室,成都 610225;2 昌都市氣象局大氣探測中心,昌都 854000)

    0 引言

    云-降水物理過程的觀測在研究天氣變化和參數(shù)化模擬等方面都具有重要意義。云中微物理過程復(fù)雜并且多種因素綜合作用下產(chǎn)生了雨滴,它在不同經(jīng)緯度、不同海拔高度、不同降雨類型下變化很大。青藏高原位于中國西南部,平均海拔高度超過4000 m,覆蓋面積超過300萬 km,是世界上海拔最高、面積最大的高原。從2014年第三次青藏高原大氣科學(xué)試驗外場觀測全面展開以來,研究人員取得了階段性進展。該項目目前已完成2014年云降水物理過程外場觀測試驗,并同步開展地面雷達和加密探空觀測,獲得高原地區(qū)多種云降水綜合觀測數(shù)據(jù)。

    在研究云-降水的微物理機制以及宏觀結(jié)構(gòu)特性時,架構(gòu)簡單、數(shù)據(jù)穩(wěn)定的雨滴譜分布,是人們重點獲取的微物理參量之一。雨滴譜分布描述了單位體積空氣、單位尺度間隔內(nèi)的雨滴個數(shù),是反射率因子、降雨率及液態(tài)含水量的計算中最為關(guān)鍵的信息。雨滴譜直接反映了在云降水中碰撞、合并、破碎的降水粒子在下落過程中的增長和演變,揭示了云滴向雨滴轉(zhuǎn)化形成降水過程的轉(zhuǎn)化過程的機制,為云降水的微物理機制和動力過程提供參照,對云-降水物理過程研究和提高參數(shù)化模擬水平都有重要意義。宮福久和陳寶君等利用3個典型的積雨云、層狀云和積層混合云降水的雨滴譜資料研究三類云降水的不同,發(fā)現(xiàn)層狀云雨滴譜比較窄,呈指數(shù)分布,積雨云和積層混合云降水譜寬,且在大水滴一方呈多峰結(jié)構(gòu);王可法等利用 Parsivel激光雨滴譜儀連續(xù)觀測太湖地區(qū)自然降水,發(fā)現(xiàn)對流云降水采樣的粒子譜數(shù)據(jù)常常包含有少數(shù)異常高、低數(shù)據(jù)。指出實際觀測中的少數(shù)小粒子有極大下落速度,主要原因是對流云降水容易產(chǎn)生部分過小雨滴粒子,其散射吸收能力較弱并且入射激光可能衍射或透射造成,而個別較低速度的大粒子主要是由儀器測量中雨滴重疊所造成。

    云物理研究指出,雨滴譜分布一般呈指數(shù)分布。Marshall和Palmer得出的表達式為:

    其中:

    N

    (

    D

    )為雨滴密度分布函數(shù),

    N

    λ

    為兩個參數(shù),并且

    λ

    (cm)和雨強

    R

    (mm/h)有如下關(guān)系:

    λ

    =41

    R

    N

    為常數(shù),為8000 m/mm。由于M-P分布具有一般雨滴譜的共同特點,從而得到廣泛應(yīng)用;Takeuchi、Ulbrich和Willis分別用Γ分布擬合實際譜,通過和M-P分布對比分析,發(fā)現(xiàn)Γ分布大大提高了在小滴和大滴區(qū)段的擬合精度;熊飛麟利用激光雨滴譜資料,研究了估算雨滴譜Gamma分布三參數(shù)的階矩法,提出一種新矩量組合M矩量組合分別代表雨滴數(shù)濃度,液態(tài)水含量和雷達反射率因子三個物理量,使得計算方便物理意義清晰。并且在實際應(yīng)用中若能通過探測手段獲得三個矩量,就能通過這一矩量組合反演得到云中雨滴譜;霍朝陽等通過對雨滴譜Gamma擬合多種方法的比較并結(jié)合C-FMCW雷達對降水系統(tǒng)垂直結(jié)構(gòu)的精細反映,獲得了對復(fù)雜降水過程進行精確到分鐘的降水類型分辨能力,精確反演降水方面得到很大的提高;張洪勝等,根據(jù)南京地區(qū)雨滴譜觀測資料,利用改進Gamma分布研究了南京地區(qū)降水雨滴譜擬合模型的參數(shù)特征發(fā)現(xiàn)

    μ

    λ

    間存在較好的二項式關(guān)系;于建宇等,利用那曲地區(qū)型地面激光雨滴譜儀的降水資料運用M-P擬合和Gamma擬合的方法,得到三類云降水的雨強貢獻均來自于大雨滴譜(即直徑大于1 mm的雨滴);Mehta Shyam利用中距離微雨雷達(MRR-2)觀測得到雨滴尺寸分布(DSD)隨高度和雨速的數(shù)據(jù)。并且提出了一種新的Gamma雨滴譜分布(DSD)模型。運用指數(shù)分布、對數(shù)正態(tài)分布和Gamma分布等不同的DSD模型對觀測數(shù)據(jù)進行擬合。結(jié)果發(fā)現(xiàn)無論是低階矩還是高階矩,擬合結(jié)果都具有Gamma DSD的特性。

    降水云的種類繁多,并且云降水發(fā)生以及分布具有不確定性,導(dǎo)致對降水的研究變得復(fù)雜且龐大。目前主要的研究內(nèi)容可以分為兩個部分:宏觀的結(jié)構(gòu)特征、微觀的物理特征。潘曉等夏季高原降水主要以深厚對流降水為主(67.8%),淺薄降水次之(26.4%),深厚強對流降水出現(xiàn)最少(5.8%);常祎等利用2014年7月1日—8月31日期間在西藏那曲的觀測數(shù)據(jù)發(fā)現(xiàn)高原雨滴譜分布相對于同緯度和季節(jié)的平原地區(qū)較寬,導(dǎo)致高原對流產(chǎn)生降水;楊俊梅等對山西汾陽地區(qū)2次降水進行分析,通過太原站多普勒雷達回波的強度大小,把雷達回波<25 dBz判斷為層狀云降水,回波強度>35 dBz判斷為對流云降水。

    就目前國內(nèi)外研究進展來看,高原地區(qū)的云降水物理過程主要是從天氣分析統(tǒng)計還有數(shù)值模擬試驗兩方面進行研究,典型的天氣個例用常規(guī)觀測資料、再分析資料還有衛(wèi)星資料等。高原地區(qū)高海拔、自然條件惡劣,對高原云降水物理過程的觀測資料匱乏,缺少有效的高分辨率資料,對高原云降水微觀參數(shù)特征和精細化結(jié)構(gòu)的研究較少,對高原云降水的分類的研究有限,遠遠不能滿足科學(xué)研究和實際應(yīng)用的要求。本文使用數(shù)據(jù)來自“第三次青藏高原氣象科學(xué)試驗”和青藏高原云-降水物理過程及大氣水循環(huán)的觀測和機理研究中獲取的觀測資料,由于加入激光雨滴譜儀器、Ka毫米波云雷達、微波輻射計和飛機觀測使得觀測技術(shù)手段與以往相比更全面觀測技術(shù)手段與以往相比更全面、可靠,更能有效的揭示青藏高原地區(qū)大氣動力過程對云微物理過程的影響,深入了解高原降水造成的微物理參數(shù)特征。對改進模式云物理過程參數(shù)化,提高數(shù)值預(yù)報水平具有重要意義。

    1 資料與方法

    本文所使用的雷達資料是2019年7月2日—8月31日在青海玉樹市氣象局觀測到的Ka波段毫米波云雷達數(shù)據(jù)。Ka波段毫米波云雷達探測高度為30~15 km,探測時間分辨率為5 s,垂直分辨率為30 m。雷達體質(zhì)采用全固態(tài)、全相參多普勒、脈沖壓縮,能夠在長時間的觀測中穩(wěn)定獲得探測資料。

    本文所使用的雨滴譜資料是2019年7月12日—8月31日在青海玉樹市氣象局的HSC-OTT Parsivel 2激光降水粒子譜儀連續(xù)觀測得到。HSC-OTT Parsivel 2激光降水粒子譜儀是一種由北京華創(chuàng)風(fēng)云集團聯(lián)合德國OTT北京公司生產(chǎn)的以激光測量為基礎(chǔ)的光學(xué)傳感器。其測量原理是:當(dāng)有降水粒子通過一維光束帶下落時,光信號轉(zhuǎn)換為電子信號,通過測量遮擋物的寬度和經(jīng)過的時間計算降水粒子的尺度和速度,并對地面雨滴譜信息自動分類。雨滴譜儀將降水粒子按照等效體積直徑和落速大小劃分為32個等級,有效測量范圍分別為0.2~25 mm和0.2~20 m/s。由于直徑大于6 mm的雨滴已經(jīng)在下落過程中破碎,因此剔除雨滴譜儀觀測記錄中直徑大于6 mm的樣本。

    雨滴譜儀實測數(shù)據(jù)的物理意義是指平均直徑

    D

    (1≤

    i

    ≤32)和下落末速度

    V

    (1≤

    j

    ≤32)對應(yīng)的降水粒子個數(shù)

    A

    ,該雨滴譜儀取樣面積54 cm,采樣間隔為1 min。那么與直徑通道

    D

    對應(yīng)單位空間體積的降水粒子總數(shù)

    N

    (

    D

    )可以表示為:

    式中:

    N

    (

    D

    )是每單位體積每單位尺寸間隔的雨滴數(shù)量濃度。

    T

    S

    分別表示雨滴譜的采樣周期和采樣面積,

    S

    =54 cm,雨滴譜儀觀測數(shù)據(jù)的時間分辨率

    T

    為60 s。由此可以進一步計算得到粒子數(shù)總密度

    N

    、雷達反射率因子

    Z

    (mm/m)和降雨強度

    R

    (mm/h),計算公式如下:

    雷達反射率因子

    Z

    降雨強度

    R

    式中:

    V

    (

    D

    )(m/s)表示直徑

    D

    的雨滴譜對應(yīng)下落末速度。分析雨滴譜的粒徑變化時,中值體積直徑

    D

    是指從小滴起累加降水量達到總降水量一半處的雨滴直徑,它代表了半數(shù)雨滴的直徑水平,因此用

    D

    說明降水過程中雨滴大小的變化。

    雨滴譜的參數(shù)化是運用數(shù)學(xué)模型來描述自然界降水粒子濃度隨尺度的變化。有關(guān)雨滴譜分布的模型很多,如M-P分布、Gamma分布。其中,Gamma分布模型的表達式為公式(1)。

    2 降水類型分類

    有關(guān)降水云的劃分有許多種辦法,傅云飛等根據(jù)雷達回波特征和雨滴譜雨強進行降水分類,將云降水分為層狀云降水、混合云降水和積雨云降水三種類型;本文參照柳臣中等的方法,按照天氣學(xué)云的定義,采用雨滴譜雷達回波和Ka波段毫米云雷達圖像將雨滴譜儀記錄到降水?dāng)?shù)據(jù)的云進行分類。人工篩選統(tǒng)計了共計6136份雨滴譜資料,把連續(xù)降水時間大于15 min定為一次有效降水,將篩選出的雨滴譜資料劃分為70次降水。由于雷達回波強度<10 dBz的資料對降水的貢獻微乎其微,因此將雷達回波強度<10 dBz的數(shù)據(jù)資料視為非降水回波。明顯降水的回波因子在30~50 dBz,本文中雷達反射率因子

    Z

    >35 dBz為積雨云降水,

    Z

    <30 dBz為層狀云降水,處于30~35 dBz為混合云降水。

    降水類型統(tǒng)計如表1和表2所示。從表1可以看出,按照頻次來看玉樹地區(qū)主要降水類型為層狀云。層狀云降水占到總降水次數(shù)的57.14%,積雨云降水次之,混合云降水次數(shù)最少。7月中下旬(20 d)層狀云降水和積雨云降水的次數(shù)與8月(31 d)相比差距不大。主要是在混合云降水上有浮動,7月中下旬(20 d)混合云降水頻次要低于8月(31 d)。從表2來看,積雨云降水的樣本數(shù)最多占比40.67%,層狀云次之37.00%,混合云降水樣本數(shù)最少22.31%。與潘曉等統(tǒng)計的夏季高原降水主要以深厚對流降水為主(67.8%),淺薄降水次之(26.4%),深厚強對流降水出現(xiàn)最少(5.8%)有出入。統(tǒng)計出的積雨云降水的樣本數(shù)略高于層狀云的樣本數(shù),原因是層狀云降水的次數(shù)雖然多,但單次有效降水的持續(xù)時間并沒有積雨云降水時間長。

    表1 玉樹地區(qū)三類降水的次數(shù)統(tǒng)計Table 1 The statistics of frequency for three types of rainfall in Yushu region

    表2 玉樹地區(qū)三類降水樣本統(tǒng)計Table 2 The statistics of precipitation samples for three types of rainfall in Yushu region

    3 結(jié)果分析

    3.1 三類降水微物理特征參量

    青海玉樹地區(qū)7月中下旬—8月的觀測期間內(nèi)雨滴譜測量總降水量為125 mm。通過對70次樣本降水的雨滴譜進行平均,得到每個降水過程的1 min的平均數(shù)據(jù)(圖1)。青藏高原東部玉樹地區(qū)雨滴譜儀測量的單次降雨量平均值為1.80 mm,其中7月中下旬單次降雨平均降雨量1.78 mm,8月單次降雨平均降雨量為1.82 mm,單次最大降雨量為11.81 mm。與常祎等的研究結(jié)果略有差別(2014年那曲地區(qū)7月日平均降水量3.22 mm,8月為5.24 mm)。單次雨強極大值的平均值為3.5 mm/h,其中7月中下旬雨強平均值為3.32 mm/h,8月雨強平均值為3.72 mm/h。表明玉樹地區(qū)降雨量小,持續(xù)時間短,降雨頻次多。

    圖1 玉樹地區(qū)70次雨強極大值和累計降雨量Fig.1 70 times rainfall intensity and accumulated rainfall in Yushu region

    特征直徑包括

    D

    平均粒徑、

    D

    平均體積直徑、

    D

    眾數(shù)直徑、

    D

    優(yōu)勢直徑、

    D

    中數(shù)直徑、

    D

    中數(shù)體積直徑。對表3中特征直徑進行綜合分析,平均粒徑

    D

    中積雨云降水的粒徑是最大的,層狀云降水的粒徑是最小的,混合云降水的數(shù)值介于兩種降雨特征之間;三類云降水的眾數(shù)直徑

    D

    數(shù)值相同,表明玉樹地區(qū)降水中最大頻率直徑在0.562 mm左右是最多的;平均體積直徑

    D

    和眾數(shù)體積直徑

    D

    有著良好的正相關(guān)關(guān)系,也有助于區(qū)分三類云降水;對含水量貢獻最多的直徑

    D

    混合云降水和積雨云降水的值相同為1.375 mm,而層狀云明顯低于前者。

    表3 三類云降水微物理特征參數(shù)Table 3 Microphysical parameters for three types of rainfall

    為了獲得觀測期間玉樹地區(qū)的雨滴譜分布特征,在處理數(shù)據(jù)的時候?qū)y(tǒng)計出的70次降水進行標(biāo)準(zhǔn)化處理:對單次降水的雨滴譜進行平均,得到每個降水過程的1 min的平均數(shù)據(jù)。歸類到三類云降水中,得到的不同降雨類型的平均譜圖數(shù)據(jù)。圖2三條線各代表三類云降水,積雨云降水Cb在圖像的最上層,層狀云降水St在最下層,混合云降水Sc處于二者之間。表明相同粒徑下,雨滴譜濃度大小為Cb>Sc>St,并且濃度越大譜寬越寬。三類云降水的譜圖都呈現(xiàn)單峰形,峰值位于0.562 mm左右。隨著雨滴粒徑的增大,滴譜濃度數(shù)隨之減弱,變化趨勢大致相同。在2.5 mm附近濃度隨著粒徑增大而減小的趨勢變小,這可能是下落時小雨滴之間的碰撞合并或者大雨滴破碎有關(guān)。積雨云、混合云和層狀云三類云降水最大雨滴粒徑分別為5.2 mm、3.7 mm和3.5 mm。

    圖2 玉樹地區(qū)三類云降水雨滴譜濃度分布(灰色點劃線為總體、紅色實線為積雨云降水、黃色劃線為混合云降水、藍色點線為層狀云降水)Fig.2 Raindrop size distributions for three types of rainfall in Yushu region(the gray dotted line is overall,the red solid line is cumulonimbus precipitation,the yellow dotted line is stratocumulus rainfall,the blue dotted line is stratiform rainfall)

    3.2 Gamma 分布的擬合特性

    運用階矩法估算雨滴譜參數(shù)參照熊飛麟等的算法,對于Gamma分布的雨滴譜來說定義第

    n

    階距為:

    積分后得到:

    在Gamma分布中,每個矩量由

    N

    、

    λ

    μ

    這三個分布參數(shù)決定,求Gamma分布三個參數(shù)就需要這三個矩量,通常把階數(shù)0、1稱為低階,2、3、4稱為中階,5、6稱為高階。本文采用0、3和6階距來對Gamma分布模型中的3個參數(shù)進行估算(表4)。

    表4 三類降水Gamma分布擬合參數(shù)Table 4 Average values of standardized Gamma distribution parameters for three types of rainfall

    通過Gamma擬合所得到的參數(shù)可以看出,積雨云濃度參數(shù)

    N

    的值最小,lg

    N

    的值為3.68 m/mm,形狀因子

    μ

    為0.347、斜率因子

    λ

    為3.280 mm;層狀云濃度參數(shù)最大,lg

    N

    的值為4.20 m/mm,形狀因子

    μ

    為1.16、斜率因子

    λ

    為5.830 mm;混合云處于兩者之間,lg

    N

    的值為3.957 m/mm,形狀因子

    μ

    為0.667、斜率因子

    λ

    為4.360 mm。與張洪勝等南京地區(qū)lg

    N

    的值為7.2 m/mm,形狀因子

    μ

    為11.6、斜率因子

    λ

    為9.2 mm相對比,3個參數(shù)都偏小,表明雨滴譜分布地域性差異較大。與李山山高原東坡瀘定站lg

    N

    為3.70 m/mm,石渠站lgN為3.16 m/mm相比,玉樹地區(qū)整體的數(shù)值略大于高原東坡地區(qū)。三類降水的Gamma分布擬合與M-P分布擬合如圖3所示。兩種擬合之間的差距不大,積雨云Cb的譜寬最寬,最大粒徑為4 mm;混合云Sc其次,最大粒徑為3.3 mm;層狀云St譜寬最窄,最大粒徑為2.5 mm。對比擬合曲線與實測數(shù)據(jù)的差異性可以更加客觀地對比模擬的適應(yīng)度。通過數(shù)據(jù)對比,記錄0.312 mm以下的數(shù)據(jù),3類云降雨的差異均超過200%。模擬較好的粒徑區(qū)間為0.562~5.5 mm,積雨云降水的差異為28.47%、混合云的差異為11.18%、層狀云的差異為9.9%。因此,玉樹地區(qū)雨滴譜在適當(dāng)范圍內(nèi)是能夠運用Gamma模擬的。

    圖3 三類云降水的平均雨滴譜分布與M-P、Gamma擬合曲線(灰色點劃線代表總體、紅色實線為實測雨滴譜分布、綠色點劃線為M-P擬合分布、黑色點線為Gamma擬合分布)(a)積雨云降水Cb;(b)混合云降水Sc;(c)層狀云降水StFig.3 Raindrop size distributions and M-P,Gamma-distribution fitting curves for three types of rainfall(the gray dotted line represents the whole,the red solid line represents the measured raindrop size distribution,the green dotted line represents the M-P distribution,and the black dotted line represents the Gamma distribution)(a) cumulonimbus rainfall,(b) stratocumulus rainfall,(c) stratiform rainfall

    3.3 擬合參數(shù)斜率因子λ 與形狀因子μ 之間的關(guān)系

    μ

    -λ的經(jīng)驗關(guān)系已經(jīng)被廣泛證明能夠更好的描述降雨云中的雨滴譜變化,表5所示大量研究表明在不同降水類型、不同氣候條件、不同地理環(huán)境之間

    μ

    -λ的經(jīng)驗關(guān)系各不相同。圖4中紅色實線積雨云關(guān)系線,黃色點劃線為混合云關(guān)系線,藍色點線為層狀云關(guān)系線。從劃分出來的三類云降水?dāng)?shù)據(jù)通過二次多項式進行擬合得到

    μ

    -λ的經(jīng)驗關(guān)系:

    圖4 實測μ-λ散點圖以及三類降水?dāng)M合曲線(紅色實線積雨云關(guān)系線,黃色點劃線為混合云關(guān)系線,藍色點線為層狀云關(guān)系線)Fig.4 Scatter plots μ-λ relations and fitting curves for three types of rainfall(red solid line:cumulonimbus relationship line,yellow dotted line:stratocumulus cloud relationship line,blue dotted line:stratiform cloud relationship line)

    表5 各地區(qū)μ-λ的經(jīng)驗關(guān)系Table 5 The μ-λ relations of other parts of the world

    可以看出:任意降水云的模擬都具有較好的二項式關(guān)系,隨著形狀因子

    μ

    的增大,

    λ

    也相應(yīng)增加;當(dāng)

    μ

    小于1時,層狀云的斜率因子

    λ

    最大,積雨云最小,混合云處于兩者之間。而當(dāng)

    μ

    大于1時,情況出現(xiàn)翻轉(zhuǎn),積雨云的斜率因子

    λ

    最大,層狀云其次,混合云最小。原因是層狀云雨滴譜濃度分布主要貢獻在偏小的小粒子上,譜寬最窄,使得其斜率的變化值最小。而混合云與積雨云的雨滴譜濃度分布主要貢獻在除了較小的小粒子外,還有較大粒子的貢獻,譜寬較寬,斜率變化率較大。

    3.4 降雨強度與雷達反射率因子之間的關(guān)系

    雷達定量估計降水是通過雷達反射率因子

    Z

    和降雨強度

    R

    的關(guān)系來完成的,而雷達反射率因子是由降水粒子譜分布決定的,不同的粒子譜分布會得到不同的雷達反射率因子。雷達反射率因子

    Z

    和降雨強度

    R

    的計算表達式如下所示

    帶入公式(4)和(5)根據(jù)對前面降水的分類,對三類降水的

    Z

    -

    R

    關(guān)系進行計算,擬合了三類降水的

    Z

    -

    R

    分布情況。計算出的

    R

    值為0~2.8,表明高原降水的雨強遠遠低于同類型的東西部地區(qū)。隨著雨強的增大,反射率因子也逐漸變大,二者有良好的線性關(guān)系。圖5所示,紅色實線積雨云關(guān)系線,黃色點劃線為混合云關(guān)系線,藍色點線為層狀云關(guān)系線,灰色點劃線為天氣雷達關(guān)系線。散點圖中層狀云降水主要出現(xiàn)在0~0.75 mm/h內(nèi),相應(yīng)的雷達回波因子是最小的;混合云集中在0~1.75 mm/h內(nèi);積雨云在0~2.8 mm/h內(nèi)都有分布。在不同類型降水中,

    B

    值不斷減小,

    b

    值也對應(yīng)減小,積雨云降水B數(shù)值為混合云與層狀云的2倍。與天氣雷達相比,玉樹地區(qū)三類云降水過程的

    Z

    -

    R

    關(guān)系參數(shù)

    B

    值大于天氣雷達(

    B

    =300,

    b

    =1.4)。原因在于玉樹地區(qū)的降雨強度低于一般地區(qū),導(dǎo)致擬合

    B

    值較大,還有雷達自身的波長、脈寬、波束寬度不同,這些因素都會影響到

    Z

    -

    R

    關(guān)系。表6中列舉了我國南部(廣東)和東部地區(qū)(南京)的

    Z

    -

    R

    分布,與青藏高原東部地區(qū)相比較,除廣東地區(qū)層狀云降水的

    B

    值為899.45遠大于各類參數(shù)值外,其他降水都在標(biāo)準(zhǔn)天氣雷達

    B

    值為300左右。對于

    b

    值來說,高原東部與我國東部相差不大,都小于天氣雷達1.4;我國南部地區(qū)的

    B

    值大于天氣雷達1.4。從降雨類型來看,我國南部和東部地區(qū)的

    B

    值從積雨云、混合云到層狀云是逐漸遞增的;高原東部地區(qū)的

    B

    值從積雨云、混合云到層狀云是逐漸遞減的,出現(xiàn)這類現(xiàn)象的原因可能是不同地理環(huán)境以及不同氣候條件。

    圖5 三類云降水的雷達反射率因子Z與雨強R關(guān)系散點圖及擬合結(jié)果(紅色實線積雨云關(guān)系線,黃色點劃線為混合云關(guān)系線,藍色點線為層狀云關(guān)系線,灰色點劃線為天氣雷達關(guān)系線)Fig.5 Scatter plot of radar reflectivity factor Z versus rain rate R and the fitted relationship for three types of rainfall(red solid line:cumulonimbus relationship line,yellow dotted line:stratocumulus cloud relationship line,blue dotted line:stratiform cloud relationship line,gray dotted line:weather radar relationship line)

    表6 玉樹地區(qū)以及其他地區(qū)降水的Z-R關(guān)系Table 6 The statistics of Z-R relationship in precipitation in Yushu region and others

    3.5 降水云的微物理結(jié)構(gòu)

    利用玉樹地區(qū)Ka波段雷達資料,通過分類后的三類降水資料中選取典型的三次降水過程分析降水云的微物理結(jié)構(gòu)。分別是8月23日積雨云、8月19日混合云、8月12日層狀云。7月中下旬—8月玉樹地區(qū)整個Ka波段毫米云雷達數(shù)據(jù)測得最大雷達回波為27 dBz。

    積雨云8月23日23時開始,云雷達上空出現(xiàn)了積雨云,此時回波反射因子不超過10 dBz,反射率因子<20 dBz,徑向速度<3 m/s,速度譜寬在0.6 m/s左右,在23:20左右地面開始觀測到有雨滴降落,云底部距離地面1 km或已經(jīng)接地。23:30左右,云頂高度上升到接近11 km左右,出現(xiàn)了強回波區(qū),強回波域在離地2~3 km,圖6a中顯示為紅色與暗紅色位置,回波強度達到23 dBz,徑向速度大于9 m/s,表明有較大粒子出現(xiàn)。譜寬值出現(xiàn)大于1.8 m/s的區(qū)域,說明伴隨雨滴下落的還有冰相粒子存在。結(jié)合圖7a,雨滴譜分布為單峰形粒子數(shù)濃度在124 m/mm,粒子直徑從0.4~2 mm不等,降雨強度為1.63 mm/h,最大粒徑位于2.37 mm位置。隨后云中回波強度減弱但速度譜寬大于1.2 m/s的區(qū)域不斷增加,粒子數(shù)濃度積雨不斷增加。到23:45左右,反射率因子再次出現(xiàn)大值區(qū),徑向速度-6 m/s區(qū)域升高到4 km位置,地面降水增大,圖7b所示粒子數(shù)濃度在252 m/mm,降雨強度達到2.88 mm/h,地面累計降水為0.68 mm。在00:00和00:30中又有云團發(fā)生發(fā)展,但雷達回波逐漸減弱,云層逐漸消散出現(xiàn)了亮帶,位于2.2 km高度上,積雨云降水過程基本結(jié)束。此處過程累計降水為8.86 mm,平均雨強為3.7 mm/h,結(jié)合圖7的演變我們可以看出,積雨云降水隨著降水強度的增加,雨滴譜濃度增大緩慢,主體逐漸向著大雨滴方向移動,大粒徑不斷出現(xiàn),并且小粒徑的數(shù)濃度不斷增加。

    圖6 Ka波段毫米云雷達觀測到的2019年8月23日22:45—次日00:45(北京時,下同)玉樹地區(qū)的積雨云降水(a)反射率因子;(b)徑向速度;(c)速度譜寬;(d)降雨量Fig.6 Time-height distributions of reflectivity (a),radial velocity (b),velocity spectral width (c),rainfall (d) by the Ka-band millimeter cloud radar during 22:45 BT 23 to next day 00:45 BT August 2019 in Yushu region

    圖7 2019年8月23—24日15 min內(nèi)積雨云雨滴譜分布演變規(guī)律(a)23:45;(b)00:00;(c)00:15;(d)00:30;(e)00:45;(f)01:00Fig.7 Raindrop size distributions of cumulonimbus rainfall in 15 min(a) 23:45,(b) 00:00,(c) 00:15,(d) 00:30,(e) 00:45,(f) 01:00

    8月19日19:00,雷達上空存在兩層云,高層云回波較弱,均小于0 dBz,在底層對流云消亡后隨之消失。底層為發(fā)展中的對流云回波此時回波反射因子在15 dBz左右,徑向速度小于0 m/s,速度譜寬小于0.6 m/s。在19:40左右云頂高度在7 km左右,離地面2 km左右出現(xiàn)強回波區(qū)域,回波強度達到23.5 dBz,徑向速度大于9 m/s,譜寬值出現(xiàn)大于1.5 m/s的區(qū)域。地面也開始觀測到有雨滴降落,圖8a中粒子數(shù)濃度為49.06 m/mm,粒子直徑在0.4~2 mm,最大粒徑為2.12 mm,此時降雨強度為0.32 mm/h。在19:45后發(fā)展最為旺盛,最大雨強為5.99 mm/h 粒子數(shù)達到峰值,最大粒徑為2.75 mm,隨后雷達回波逐漸減弱降水過程結(jié)束。此次過程累計降水為0.94 mm,平均雨強為1.37 mm/h,結(jié)合圖8混合云降水雨滴譜分布形態(tài)變化明顯,三個時次的雨滴譜濃度分別為49.06 m/mm,406.85 m/mm,250.67 m/mm,粒子增長速度較快。

    圖8 2019年8月19日15 min內(nèi)混合云雨滴譜分布演變規(guī)律(a)19:30;(b)19:45;(c)20:00Fig.8 Raindrop size distributions of stratocumulus rainfall in 15min(a) 19:30,(b) 19:45,(c) 20:00

    8月12日23:00時,雷達上空出現(xiàn)了較弱的對流云,雷達回波最大不超過15 dBz,徑向速度數(shù)值偏小,零星值達到6 m/s;譜寬值在0.6 m/s左右,云頂高度為6~7 km。并且有明顯的零度層亮帶位于離地面高度2 km左右的位置,雷達反射率在其上部隨高度遞減,在其下部隨著高度的降低。說明含有冰相粒子在下落過程中,上層的冰相粒子下落至在零度層后有融化過程,冰相粒子融化為水滴,水滴碰并增長變成大水滴,導(dǎo)致雷達反射率變強。到了零度層以下,由于雨滴的下落速度較快,降水粒子的濃度反而變小,雷達反射率強度反而變低。雨滴譜濃度分布為單峰形整體變化不大,平均雨強為2.6 mm/h,累計降水為4.28 mm,平均數(shù)濃度為409.56 m/mm。

    綜上所述,玉樹地區(qū)不同云降水所對應(yīng)的雷達回波和雨滴譜分布演變規(guī)律有著明顯的差別。積雨云回波強度比較強,最大值過20 dBz。云頂發(fā)展高度較高,并且零度層亮帶難以觀測到;強回波區(qū)域在離地2~3 km左右?;旌辖邓陌l(fā)生發(fā)展高度位于積雨云回波強度類似,但最大值在20 dBz左右;云頂發(fā)展高度較高,接近13 km左右,沒有零度層亮帶出現(xiàn),強回波區(qū)域在離地2~3 km左右;持續(xù)時間較短,垂直尺度大于水平尺度。層狀云降水中,回波強度比較弱,最大值低于15 dBz。云頂發(fā)展高度在6~7 km左右較為平整;能夠直接看見零度層亮帶,位于2 km左右,強回波區(qū)域與零度層亮帶同高度。積雨云中雨滴譜濃度增大緩慢,主體逐漸向著大雨滴方向移動,大粒徑不斷出現(xiàn),并且小粒徑的數(shù)濃度不斷增加,持續(xù)時間較長?;旌显浦杏甑巫V濃度增大和減弱的時間較短,持續(xù)時間弱于積雨云。此次層狀云降水時間較長,各項雨滴譜參數(shù)都較為穩(wěn)定整體變化較小。三類云降水的譜圖都呈現(xiàn)單峰形,峰值位于0.562 mm粒徑左右。隨著雨滴粒徑的增大,滴譜濃度數(shù)隨之減弱,變化趨勢大致相同。連續(xù)降水使得不斷增加的小雨滴之間的碰撞合并在1.375 mm粒徑位置第一次出現(xiàn)了濃度隨著粒徑的增大而減小的趨勢變小的情況;第二次出現(xiàn)這種情況是在2.5 mm粒徑左右的附近,原因為隨著降水過程持續(xù)時間變長或降雨強度增加,較大雨滴破碎成大雨滴使得濃度變小的趨勢變小。與霍朝陽年分析廣東層狀云降水雨滴譜濃度圖相比,廣東地區(qū)層狀云譜圖只出現(xiàn)了第二次減緩趨勢。與于建宇年分析青藏高原那曲地區(qū)積雨云降水,譜圖變化趨勢相同,但那曲地區(qū)積雨云第二次減緩趨勢更加明顯。與Wang等研究高原東坡稻城地區(qū)雨滴譜相比,玉樹地區(qū)的譜寬更窄。

    4 結(jié)論與討論

    本文利用青藏高原玉樹氣象局的HSC-OTT Parsivel 2激光降水粒子譜儀在2019年7月中下旬至8月觀測得到的70次降水過程的雨滴譜資料和Ka波段云雷達觀測資料,研究青藏高原東部玉樹地區(qū)降雨的雨滴譜特征,主要結(jié)論如下:

    (1)高原地區(qū)7月中下旬至8月總降水量為125 mm,兩個月降雨強度相當(dāng),玉樹地區(qū)主要降水類型為層狀云,降水占到總降水次數(shù)的57.14%;積雨云降水次之,混合云降水次數(shù)最少。三類云降水最大雨滴粒徑分別為5.2 mm、3.7 mm和3.5 mm,平均粒徑為0.801 mm。三類云降水的譜圖都呈現(xiàn)單峰形,峰值位于0.562 mm左右。相同粒徑下,雨滴譜濃度大小為Cb>Sc>St,并且濃度越大譜寬越寬。

    (2)三類云在玉樹地區(qū)的Gamma分布擬合模擬較好的粒徑為0.562~5.5 mm,但各個粒徑區(qū)間內(nèi)的擬合值都有高估的情況,特別是對小雨滴濃度值高估明顯;任意降水云的模擬都具有較好的二項式關(guān)系,隨著形狀因子

    μ

    的增大,

    λ

    也相應(yīng)增加。層狀云雨滴譜濃度分布主要貢獻在小粒徑上且譜寬最窄,使得其斜率的變化值最小。而混合云與積雨云的雨滴譜濃度分布主要貢獻分布均勻且譜寬較寬所以,變化率較大;積雨云、混合云和層狀云的

    Z

    -

    R

    關(guān)系擬合,隨著雨強的增大,反射率因子也逐漸變大,二者有良好的線性關(guān)系,擬合的

    B

    值從積雨云、混合云到層狀云是逐漸遞減的。

    (3)選取個例中三類降水云的結(jié)構(gòu)相同點為零度層位于2 km左右,降雨時云底高度很低或接地。積雨云和混合云發(fā)生發(fā)展旺盛,最高云高可達13 km,但混合云持續(xù)時間短;層狀云發(fā)展較為平穩(wěn),云高在6~7 km。連續(xù)降水使得不斷增加的小雨滴之間的碰撞合并在1.375 mm粒徑位置第一次出現(xiàn)了濃度隨著粒徑增大而減小的趨勢變小的情況;第二次出現(xiàn)這種情況是在2.5 mm粒徑左右的附近,原因為隨著降水過程持續(xù)時間變長或降雨強度增加,較大雨滴破碎成大雨滴使得濃度變小的趨勢變小。

    本次雨滴譜分布受到地理環(huán)境的局限性,僅適用于玉樹地區(qū)。按照降雨次數(shù)來劃分使得數(shù)據(jù)更具有代表性,層狀云的估值較為準(zhǔn)確,但在統(tǒng)計中數(shù)據(jù)對于積雨云的低估使得積雨云降水整體數(shù)據(jù)偏小,將在今后的研究中進一步深入研究。

    猜你喜歡
    積雨云玉樹層狀
    玉樹留芳
    成都積雨云形成的層結(jié)條件分析
    軋制復(fù)合制備TA1/AZ31B/TA1層狀復(fù)合材料組織與性能研究
    重型機械(2019年3期)2019-08-27 00:58:44
    玉樹花開、雞鳴東白
    寶藏(2018年12期)2019-01-29 01:51:40
    穿越積雨云
    蔡俊的詩
    青春(2018年1期)2018-01-09 18:34:16
    兩級結(jié)構(gòu)層狀Ti-TiBw/Ti復(fù)合材料擴散焊及其拉伸行為
    焊接(2016年9期)2016-02-27 13:05:22
    高韌性抗層狀撕裂Q345FTE-Z35鋼板開發(fā)
    新疆鋼鐵(2015年2期)2015-11-07 03:27:52
    層狀雙氫氧化物處理工業(yè)廢水的研究進展
    玉樹依然美麗
    美女免费视频网站| 桃色一区二区三区在线观看| 免费观看精品视频网站| 国产免费男女视频| 亚洲精品一区av在线观看| 国内揄拍国产精品人妻在线| 少妇猛男粗大的猛烈进出视频 | a级毛片a级免费在线| 精品久久久久久久末码| 日本撒尿小便嘘嘘汇集6| 精品久久久噜噜| 成年版毛片免费区| 在线观看一区二区三区| 一个人观看的视频www高清免费观看| 国产亚洲欧美98| 99热精品在线国产| 日本黄色视频三级网站网址| 美女黄网站色视频| 久久九九热精品免费| 直男gayav资源| 色综合色国产| 成人国产麻豆网| 天美传媒精品一区二区| 特大巨黑吊av在线直播| 免费无遮挡裸体视频| 精华霜和精华液先用哪个| 黄色视频,在线免费观看| 久99久视频精品免费| 天堂网av新在线| 日本-黄色视频高清免费观看| 一进一出好大好爽视频| 禁无遮挡网站| 国产在线男女| 2021天堂中文幕一二区在线观| 男女边吃奶边做爰视频| x7x7x7水蜜桃| 日韩在线高清观看一区二区三区 | videossex国产| 蜜桃亚洲精品一区二区三区| 色综合色国产| 在线免费观看的www视频| 久久午夜福利片| 亚洲男人的天堂狠狠| 欧美精品啪啪一区二区三区| 成人国产麻豆网| 欧美zozozo另类| 极品教师在线视频| 婷婷亚洲欧美| 美女高潮喷水抽搐中文字幕| 可以在线观看的亚洲视频| av专区在线播放| 国产免费一级a男人的天堂| 他把我摸到了高潮在线观看| 一边摸一边抽搐一进一小说| 免费观看精品视频网站| 人人妻人人澡欧美一区二区| 免费观看人在逋| xxxwww97欧美| 成人高潮视频无遮挡免费网站| 亚洲成人中文字幕在线播放| 亚洲av熟女| 99久国产av精品| 亚洲欧美日韩东京热| 五月伊人婷婷丁香| 成人无遮挡网站| 尤物成人国产欧美一区二区三区| 在线天堂最新版资源| 国产91精品成人一区二区三区| 久久久久国内视频| 99精品在免费线老司机午夜| 日韩一本色道免费dvd| 男女边吃奶边做爰视频| 午夜激情欧美在线| 熟妇人妻久久中文字幕3abv| 久久亚洲精品不卡| 亚洲五月天丁香| 色尼玛亚洲综合影院| 在线观看66精品国产| 香蕉av资源在线| 啪啪无遮挡十八禁网站| 联通29元200g的流量卡| 国产男靠女视频免费网站| 黄色视频,在线免费观看| 欧美中文日本在线观看视频| 国产精品1区2区在线观看.| 久久精品国产亚洲av香蕉五月| 欧美黑人巨大hd| 又爽又黄a免费视频| 日本黄色视频三级网站网址| 国产一区二区三区在线臀色熟女| 免费看av在线观看网站| 国产精品1区2区在线观看.| 91麻豆精品激情在线观看国产| 高清毛片免费观看视频网站| 看黄色毛片网站| 国产男靠女视频免费网站| 日本免费a在线| 精品午夜福利视频在线观看一区| 熟女电影av网| 日韩强制内射视频| 午夜福利成人在线免费观看| .国产精品久久| 国产精品久久久久久精品电影| 一本久久中文字幕| 亚洲成av人片在线播放无| 国产真实伦视频高清在线观看 | 国产午夜福利久久久久久| 免费高清视频大片| 欧美3d第一页| 日本黄色视频三级网站网址| 九九热线精品视视频播放| 两人在一起打扑克的视频| 欧美+日韩+精品| 日韩中字成人| 久久久久九九精品影院| 国产精品不卡视频一区二区| 搡女人真爽免费视频火全软件 | 亚洲成a人片在线一区二区| 欧美人与善性xxx| 欧美bdsm另类| 麻豆成人午夜福利视频| 欧美人与善性xxx| 精品国产三级普通话版| 国产亚洲精品久久久com| 我要搜黄色片| 俺也久久电影网| 女的被弄到高潮叫床怎么办 | 国产精品三级大全| 国模一区二区三区四区视频| 欧美中文日本在线观看视频| 欧美最黄视频在线播放免费| 麻豆成人av在线观看| av国产免费在线观看| 中文字幕熟女人妻在线| 俺也久久电影网| 日本三级黄在线观看| 成人永久免费在线观看视频| 国产精品久久久久久亚洲av鲁大| 亚洲内射少妇av| 亚洲av免费高清在线观看| 亚洲第一电影网av| 成人高潮视频无遮挡免费网站| 无遮挡黄片免费观看| 舔av片在线| 国产成人av教育| 真人做人爱边吃奶动态| 22中文网久久字幕| 国产av不卡久久| 男人和女人高潮做爰伦理| 国产在视频线在精品| 波野结衣二区三区在线| 久久精品久久久久久噜噜老黄 | 欧美日本视频| 免费电影在线观看免费观看| 一进一出抽搐动态| 免费无遮挡裸体视频| 五月玫瑰六月丁香| 色av中文字幕| 久久久久久九九精品二区国产| 我要看日韩黄色一级片| 蜜桃久久精品国产亚洲av| 精品不卡国产一区二区三区| 男女边吃奶边做爰视频| 精品久久久久久久久久久久久| 国产精品一区二区性色av| 国产精品自产拍在线观看55亚洲| 一本精品99久久精品77| 亚洲精品日韩av片在线观看| 麻豆国产av国片精品| 中国美白少妇内射xxxbb| 亚洲美女视频黄频| 国产淫片久久久久久久久| 中文资源天堂在线| 久久99热这里只有精品18| 真人一进一出gif抽搐免费| 国产亚洲精品av在线| 99久久精品热视频| 国内精品久久久久久久电影| 久久香蕉精品热| 欧美性猛交黑人性爽| 色av中文字幕| 精品人妻1区二区| 久久午夜福利片| 精品久久久久久久人妻蜜臀av| 深夜a级毛片| 18禁裸乳无遮挡免费网站照片| 国产激情偷乱视频一区二区| 黄色女人牲交| 国产69精品久久久久777片| 国产精品亚洲一级av第二区| 超碰av人人做人人爽久久| 内地一区二区视频在线| 日日啪夜夜撸| 禁无遮挡网站| 国产白丝娇喘喷水9色精品| 亚洲熟妇熟女久久| 精品人妻一区二区三区麻豆 | eeuss影院久久| 亚洲av免费高清在线观看| 搞女人的毛片| 女生性感内裤真人,穿戴方法视频| 特大巨黑吊av在线直播| 午夜福利在线在线| 亚洲最大成人手机在线| 国产精品国产高清国产av| 国产成人福利小说| 欧美极品一区二区三区四区| 在线观看av片永久免费下载| 嫁个100分男人电影在线观看| 日韩欧美在线乱码| 天堂动漫精品| 久久久久精品国产欧美久久久| 国产三级中文精品| 国产亚洲精品久久久com| 欧美+日韩+精品| 五月玫瑰六月丁香| 欧美极品一区二区三区四区| 在线观看av片永久免费下载| 国产亚洲精品久久久com| 老司机深夜福利视频在线观看| 国产精品,欧美在线| 美女高潮的动态| 全区人妻精品视频| 午夜亚洲福利在线播放| 在线国产一区二区在线| 国产精品,欧美在线| 性欧美人与动物交配| 999久久久精品免费观看国产| 最近中文字幕高清免费大全6 | 99热网站在线观看| 蜜桃亚洲精品一区二区三区| 国产高潮美女av| 天天躁日日操中文字幕| 国产免费男女视频| 亚洲熟妇熟女久久| 国产在线精品亚洲第一网站| 国产伦精品一区二区三区视频9| 国产精品美女特级片免费视频播放器| 日韩中文字幕欧美一区二区| 亚洲欧美日韩卡通动漫| 窝窝影院91人妻| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕精品亚洲无线码一区| 女同久久另类99精品国产91| 一级毛片久久久久久久久女| a级一级毛片免费在线观看| 国产主播在线观看一区二区| 如何舔出高潮| 男人舔奶头视频| 高清毛片免费观看视频网站| 亚洲av二区三区四区| 国产一区二区三区在线臀色熟女| 成人欧美大片| 国产男人的电影天堂91| 在线播放国产精品三级| 国产不卡一卡二| 午夜激情福利司机影院| 我的老师免费观看完整版| 欧美成人a在线观看| 亚洲无线观看免费| 久久久久久久午夜电影| 欧美最黄视频在线播放免费| 88av欧美| 一本精品99久久精品77| 久久精品国产亚洲av香蕉五月| 桃色一区二区三区在线观看| 亚洲av中文av极速乱 | 99国产极品粉嫩在线观看| 国内精品宾馆在线| 成年版毛片免费区| 一区二区三区四区激情视频 | 99热精品在线国产| 亚洲18禁久久av| 日本三级黄在线观看| 亚洲自拍偷在线| 好男人在线观看高清免费视频| 日日摸夜夜添夜夜添av毛片 | 欧美高清成人免费视频www| 亚洲av免费在线观看| 国产v大片淫在线免费观看| 亚洲精品亚洲一区二区| 国产真实乱freesex| 婷婷精品国产亚洲av| 偷拍熟女少妇极品色| 国产av麻豆久久久久久久| 我的女老师完整版在线观看| 日本三级黄在线观看| 精品人妻视频免费看| 国产精品美女特级片免费视频播放器| 亚洲成人久久性| 久久久久九九精品影院| 国产视频一区二区在线看| 女人十人毛片免费观看3o分钟| 成人精品一区二区免费| 黄色女人牲交| 直男gayav资源| www日本黄色视频网| 国产中年淑女户外野战色| 日本-黄色视频高清免费观看| 日韩av在线大香蕉| 成人国产综合亚洲| 欧美在线一区亚洲| 亚洲国产精品合色在线| 亚洲av第一区精品v没综合| 亚洲成av人片在线播放无| 少妇的逼好多水| 十八禁国产超污无遮挡网站| 免费av毛片视频| 亚洲欧美精品综合久久99| 午夜免费成人在线视频| 久久精品国产亚洲网站| 午夜福利在线在线| xxxwww97欧美| 久久人人精品亚洲av| 熟女人妻精品中文字幕| ponron亚洲| 真人一进一出gif抽搐免费| 在线观看免费视频日本深夜| 久久天躁狠狠躁夜夜2o2o| 国产精品乱码一区二三区的特点| 亚洲七黄色美女视频| 中文字幕av在线有码专区| 长腿黑丝高跟| 日韩大尺度精品在线看网址| 久久人人爽人人爽人人片va| 亚洲欧美日韩高清在线视频| 国产高清视频在线观看网站| 免费在线观看日本一区| 国产一区二区在线av高清观看| 十八禁国产超污无遮挡网站| 国产成年人精品一区二区| 直男gayav资源| 一a级毛片在线观看| 欧美色欧美亚洲另类二区| 又粗又爽又猛毛片免费看| 久久久色成人| 一进一出抽搐gif免费好疼| 午夜福利视频1000在线观看| 国产一级毛片七仙女欲春2| 国模一区二区三区四区视频| 中文亚洲av片在线观看爽| 日本 欧美在线| 看十八女毛片水多多多| 性插视频无遮挡在线免费观看| 欧美激情久久久久久爽电影| 观看美女的网站| 亚洲国产日韩欧美精品在线观看| 极品教师在线免费播放| 欧美绝顶高潮抽搐喷水| 嫩草影院新地址| 国产精品99久久久久久久久| 亚洲七黄色美女视频| 伊人久久精品亚洲午夜| 91久久精品国产一区二区成人| 黄色一级大片看看| 免费搜索国产男女视频| 婷婷精品国产亚洲av| 夜夜爽天天搞| 九九在线视频观看精品| 久久国内精品自在自线图片| 国产免费一级a男人的天堂| 精品欧美国产一区二区三| 精品一区二区三区视频在线观看免费| 在现免费观看毛片| 国产精品一区二区性色av| 日韩亚洲欧美综合| 国内少妇人妻偷人精品xxx网站| 久久久国产成人免费| 国产一区二区三区视频了| 日韩亚洲欧美综合| 久久精品人妻少妇| 最新在线观看一区二区三区| 深夜a级毛片| 91精品国产九色| 日韩中字成人| 日本在线视频免费播放| 国产精品,欧美在线| 十八禁国产超污无遮挡网站| АⅤ资源中文在线天堂| 日日摸夜夜添夜夜添小说| 国产一区二区激情短视频| 真人做人爱边吃奶动态| 18禁在线播放成人免费| 亚洲欧美日韩东京热| 久久6这里有精品| 亚洲av五月六月丁香网| 最后的刺客免费高清国语| 一级黄色大片毛片| 亚洲欧美清纯卡通| 精品一区二区三区人妻视频| 国产精品美女特级片免费视频播放器| 亚洲真实伦在线观看| 韩国av在线不卡| 久久人人精品亚洲av| 成年女人毛片免费观看观看9| 精品久久久久久,| 性欧美人与动物交配| 婷婷色综合大香蕉| 又爽又黄a免费视频| 日韩欧美三级三区| 成人一区二区视频在线观看| 国产精品爽爽va在线观看网站| 国产亚洲精品av在线| 制服丝袜大香蕉在线| 久久精品国产亚洲av涩爱 | 99九九线精品视频在线观看视频| 国产免费男女视频| 最近最新免费中文字幕在线| 永久网站在线| 三级毛片av免费| 自拍偷自拍亚洲精品老妇| 日本黄色片子视频| 日本 av在线| 久久精品国产亚洲av天美| 91麻豆精品激情在线观看国产| 亚洲av成人精品一区久久| 国产国拍精品亚洲av在线观看| 日韩一区二区视频免费看| 亚洲最大成人av| 尾随美女入室| av视频在线观看入口| .国产精品久久| 毛片女人毛片| 亚洲av中文av极速乱 | 在线免费观看的www视频| 中文字幕久久专区| av在线天堂中文字幕| 熟妇人妻久久中文字幕3abv| 午夜福利在线观看吧| 九九在线视频观看精品| 能在线免费观看的黄片| 国产一区二区三区av在线 | 老司机午夜福利在线观看视频| eeuss影院久久| 99久久成人亚洲精品观看| 亚洲av免费在线观看| 美女被艹到高潮喷水动态| 亚洲性夜色夜夜综合| 国产精品99久久久久久久久| 久久午夜福利片| 国产一区二区三区视频了| 国产麻豆成人av免费视频| 午夜福利18| 成熟少妇高潮喷水视频| 亚洲av美国av| 午夜日韩欧美国产| 两个人视频免费观看高清| 在线观看免费视频日本深夜| 亚洲成人久久爱视频| 无遮挡黄片免费观看| 别揉我奶头 嗯啊视频| 深夜精品福利| 韩国av一区二区三区四区| 久久久久久久精品吃奶| 999久久久精品免费观看国产| 51国产日韩欧美| 久久人人精品亚洲av| 亚洲不卡免费看| 波多野结衣高清作品| 日韩在线高清观看一区二区三区 | 韩国av一区二区三区四区| 欧美在线一区亚洲| 国产精品一区二区性色av| 欧美色欧美亚洲另类二区| 久久精品久久久久久噜噜老黄 | 欧美日韩乱码在线| 在线观看免费视频日本深夜| 狂野欧美白嫩少妇大欣赏| 18禁黄网站禁片免费观看直播| 日本在线视频免费播放| 色综合站精品国产| 国产伦一二天堂av在线观看| 午夜福利在线观看吧| 国产高潮美女av| 女人被狂操c到高潮| 免费av毛片视频| 免费一级毛片在线播放高清视频| 观看美女的网站| 小说图片视频综合网站| 少妇人妻精品综合一区二区 | 免费av不卡在线播放| 欧美色视频一区免费| 亚洲乱码一区二区免费版| 国产精品久久电影中文字幕| 精品午夜福利在线看| 久久久久免费精品人妻一区二区| 欧美黑人巨大hd| 日韩欧美 国产精品| 国产日本99.免费观看| 3wmmmm亚洲av在线观看| 淫妇啪啪啪对白视频| 99九九线精品视频在线观看视频| 变态另类成人亚洲欧美熟女| 少妇人妻精品综合一区二区 | 在线观看66精品国产| 能在线免费观看的黄片| 国产一区二区在线观看日韩| 免费在线观看日本一区| 69人妻影院| 一区二区三区激情视频| 成年免费大片在线观看| 他把我摸到了高潮在线观看| 18+在线观看网站| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产高清在线一区二区三| 一级黄片播放器| 最新中文字幕久久久久| 国产欧美日韩精品一区二区| 欧美激情在线99| 国产免费一级a男人的天堂| 久久精品国产99精品国产亚洲性色| 香蕉av资源在线| 亚洲精华国产精华液的使用体验 | 亚洲av熟女| 精品一区二区三区人妻视频| 如何舔出高潮| АⅤ资源中文在线天堂| 又黄又爽又免费观看的视频| 男人和女人高潮做爰伦理| 深夜a级毛片| 99热这里只有是精品50| 制服丝袜大香蕉在线| 日韩中文字幕欧美一区二区| 欧美3d第一页| 五月玫瑰六月丁香| 日韩,欧美,国产一区二区三区 | 国产 一区精品| 精品久久久久久久久久免费视频| 午夜亚洲福利在线播放| 精品一区二区三区人妻视频| 精品午夜福利视频在线观看一区| 亚洲熟妇熟女久久| 老熟妇乱子伦视频在线观看| 国产精品美女特级片免费视频播放器| 精品国产三级普通话版| 波多野结衣巨乳人妻| 亚洲五月天丁香| 神马国产精品三级电影在线观看| 99久久成人亚洲精品观看| 久久久精品欧美日韩精品| 国产视频一区二区在线看| 国产亚洲欧美98| 亚洲精品色激情综合| 熟妇人妻久久中文字幕3abv| 亚洲欧美日韩高清专用| 91av网一区二区| 大又大粗又爽又黄少妇毛片口| 日日干狠狠操夜夜爽| 国产女主播在线喷水免费视频网站 | 极品教师在线视频| 哪里可以看免费的av片| 伦理电影大哥的女人| 在线免费十八禁| 欧美bdsm另类| 中文字幕av在线有码专区| 国模一区二区三区四区视频| 2021天堂中文幕一二区在线观| 十八禁国产超污无遮挡网站| 最近中文字幕高清免费大全6 | 俺也久久电影网| 国产午夜精品久久久久久一区二区三区 | 美女高潮的动态| 亚洲18禁久久av| 黄色一级大片看看| 午夜激情福利司机影院| 美女黄网站色视频| 日本熟妇午夜| 亚洲国产色片| 国产蜜桃级精品一区二区三区| 欧美日本视频| 精品欧美国产一区二区三| 少妇丰满av| 免费av观看视频| 毛片女人毛片| 在线观看一区二区三区| 国产精品人妻久久久久久| 午夜老司机福利剧场| 深夜a级毛片| 日韩一区二区视频免费看| 中亚洲国语对白在线视频| 最好的美女福利视频网| 日本一本二区三区精品| 亚洲av熟女| 久久欧美精品欧美久久欧美| 国产色爽女视频免费观看| 午夜亚洲福利在线播放| 国产爱豆传媒在线观看| 日韩欧美国产在线观看| 精品免费久久久久久久清纯| 国内精品美女久久久久久| 内地一区二区视频在线| av在线天堂中文字幕| 女同久久另类99精品国产91| 精品人妻熟女av久视频| 亚洲欧美日韩无卡精品| 露出奶头的视频| av黄色大香蕉| 免费不卡的大黄色大毛片视频在线观看 | 久久99热这里只有精品18| 久久久久久久亚洲中文字幕| 国产精品日韩av在线免费观看| 欧美激情国产日韩精品一区| 啪啪无遮挡十八禁网站| 久久精品91蜜桃| 国产中年淑女户外野战色| 少妇丰满av| 午夜精品一区二区三区免费看| 我要看日韩黄色一级片| h日本视频在线播放| 在线免费十八禁| 美女大奶头视频| 观看免费一级毛片| 久久精品综合一区二区三区| 黄色配什么色好看| 真人一进一出gif抽搐免费|