• <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è)廢水的研究進展
    玉樹依然美麗
    国产三级黄色录像| 国产欧美日韩一区二区精品| 欧美bdsm另类| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美黑人欧美精品刺激| 亚洲va日本ⅴa欧美va伊人久久| 高清在线国产一区| 午夜精品久久久久久毛片777| 日韩欧美免费精品| 一a级毛片在线观看| 女同久久另类99精品国产91| 国产综合懂色| 成人性生交大片免费视频hd| 有码 亚洲区| 免费看a级黄色片| 亚洲成av人片在线播放无| 免费av观看视频| 亚洲va日本ⅴa欧美va伊人久久| 成人特级黄色片久久久久久久| 网址你懂的国产日韩在线| 国产精品乱码一区二三区的特点| 国产激情偷乱视频一区二区| 日韩欧美在线二视频| 日韩欧美精品免费久久 | 亚洲欧美日韩高清专用| 午夜福利免费观看在线| 自拍偷自拍亚洲精品老妇| 国产精品国产高清国产av| 级片在线观看| 久久久国产成人精品二区| 精品午夜福利在线看| 成人国产一区最新在线观看| 午夜日韩欧美国产| 亚洲成人久久爱视频| 丁香欧美五月| 国产午夜精品论理片| 成人亚洲精品av一区二区| 老鸭窝网址在线观看| 少妇丰满av| 亚洲第一区二区三区不卡| 国产精品久久久久久久久免 | 99视频精品全部免费 在线| 99热这里只有精品一区| 欧美日韩乱码在线| 在线观看66精品国产| 亚洲av不卡在线观看| 免费人成视频x8x8入口观看| 内射极品少妇av片p| 国产一区二区亚洲精品在线观看| 麻豆成人av在线观看| 男女那种视频在线观看| bbb黄色大片| 欧美性感艳星| 色综合亚洲欧美另类图片| 嫩草影视91久久| 国产亚洲精品久久久久久毛片| 国产在视频线在精品| 美女 人体艺术 gogo| 男插女下体视频免费在线播放| 国产精品女同一区二区软件 | 人妻夜夜爽99麻豆av| 久久这里只有精品中国| 精品欧美国产一区二区三| 欧美成人性av电影在线观看| 国产精品av视频在线免费观看| 亚洲色图av天堂| 国产欧美日韩精品一区二区| 蜜桃亚洲精品一区二区三区| 久久性视频一级片| 亚洲性夜色夜夜综合| 亚洲国产精品999在线| 一个人看的www免费观看视频| 成人欧美大片| 国产精品永久免费网站| 欧美xxxx性猛交bbbb| 免费在线观看亚洲国产| 亚洲一区高清亚洲精品| 精品99又大又爽又粗少妇毛片 | 国产日本99.免费观看| 老鸭窝网址在线观看| 欧美一区二区亚洲| 夜夜爽天天搞| 好男人在线观看高清免费视频| 一级黄片播放器| 亚洲成av人片免费观看| 搡老熟女国产l中国老女人| 老鸭窝网址在线观看| 综合色av麻豆| 亚洲成人免费电影在线观看| 国产激情偷乱视频一区二区| av在线老鸭窝| 国产精品伦人一区二区| 无人区码免费观看不卡| 国产成人影院久久av| av视频在线观看入口| 美女大奶头视频| 国产亚洲欧美在线一区二区| 亚洲片人在线观看| 亚洲中文字幕一区二区三区有码在线看| 午夜免费男女啪啪视频观看 | 九九在线视频观看精品| 黄色丝袜av网址大全| 欧美成狂野欧美在线观看| 国产精品久久视频播放| 中文字幕熟女人妻在线| 亚洲国产欧洲综合997久久,| 精品人妻1区二区| 成人永久免费在线观看视频| 直男gayav资源| 精品久久久久久久末码| 91久久精品电影网| 怎么达到女性高潮| 成年女人毛片免费观看观看9| 女生性感内裤真人,穿戴方法视频| 美女免费视频网站| 亚洲专区中文字幕在线| 搡老熟女国产l中国老女人| 老司机福利观看| 亚洲av免费高清在线观看| 亚洲av美国av| 日韩人妻高清精品专区| 国产精品亚洲美女久久久| 免费在线观看日本一区| 成年免费大片在线观看| 最近最新免费中文字幕在线| 我的女老师完整版在线观看| 1000部很黄的大片| 一个人看的www免费观看视频| 国产成人a区在线观看| 好男人在线观看高清免费视频| 国产一区二区激情短视频| 在线观看av片永久免费下载| 日韩欧美国产在线观看| 亚洲人与动物交配视频| 亚洲国产精品成人综合色| 国产一区二区三区在线臀色熟女| 我的女老师完整版在线观看| 男女下面进入的视频免费午夜| eeuss影院久久| 亚洲av.av天堂| 精品乱码久久久久久99久播| 日韩欧美三级三区| 色吧在线观看| 五月伊人婷婷丁香| 午夜福利在线观看免费完整高清在 | 久久午夜福利片| 嫩草影院入口| 黄色视频,在线免费观看| 男女做爰动态图高潮gif福利片| 悠悠久久av| 波多野结衣巨乳人妻| 最近最新免费中文字幕在线| 综合色av麻豆| 国产91精品成人一区二区三区| 国产爱豆传媒在线观看| 18禁黄网站禁片免费观看直播| 亚洲精品日韩av片在线观看| 亚洲乱码一区二区免费版| 一个人看视频在线观看www免费| 久久精品国产亚洲av天美| 极品教师在线视频| 国产视频内射| 国产成人啪精品午夜网站| 国产视频内射| 亚洲国产色片| 欧美zozozo另类| 一级a爱片免费观看的视频| 午夜亚洲福利在线播放| 精品久久久久久久久av| 看十八女毛片水多多多| 老女人水多毛片| 国产精品久久电影中文字幕| 1024手机看黄色片| 天天躁日日操中文字幕| 亚洲18禁久久av| 亚洲av成人不卡在线观看播放网| 国产老妇女一区| 国产免费男女视频| 99热这里只有精品一区| 在线观看美女被高潮喷水网站 | 可以在线观看毛片的网站| 国内精品久久久久久久电影| 制服丝袜大香蕉在线| 国产乱人伦免费视频| 婷婷丁香在线五月| 熟女人妻精品中文字幕| 一级毛片久久久久久久久女| 免费看a级黄色片| 欧美成人免费av一区二区三区| 丝袜美腿在线中文| 在线观看av片永久免费下载| 啦啦啦韩国在线观看视频| 久久精品国产99精品国产亚洲性色| 国产高清有码在线观看视频| 麻豆久久精品国产亚洲av| 少妇的逼水好多| www.熟女人妻精品国产| 桃色一区二区三区在线观看| 又爽又黄无遮挡网站| 色噜噜av男人的天堂激情| 精品日产1卡2卡| 欧美xxxx黑人xx丫x性爽| 国产高清视频在线观看网站| 宅男免费午夜| 少妇丰满av| 男女下面进入的视频免费午夜| av中文乱码字幕在线| 国产乱人视频| 999久久久精品免费观看国产| 免费看光身美女| 午夜日韩欧美国产| 国内精品久久久久精免费| 成人美女网站在线观看视频| 一卡2卡三卡四卡精品乱码亚洲| 亚洲人成伊人成综合网2020| 性色avwww在线观看| www.www免费av| 国模一区二区三区四区视频| 性插视频无遮挡在线免费观看| 99久久无色码亚洲精品果冻| 久久久成人免费电影| 日本黄色视频三级网站网址| 亚洲 欧美 日韩 在线 免费| 91久久精品电影网| 午夜福利高清视频| 日韩欧美精品v在线| 日本五十路高清| 色哟哟哟哟哟哟| 国产一区二区三区视频了| 又黄又爽又免费观看的视频| 在线观看舔阴道视频| 激情在线观看视频在线高清| 亚洲最大成人中文| 国产精品一区二区免费欧美| 午夜精品一区二区三区免费看| 夜夜躁狠狠躁天天躁| 国产大屁股一区二区在线视频| 亚洲成av人片在线播放无| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩综合久久久久久 | 伊人久久精品亚洲午夜| 啦啦啦韩国在线观看视频| 日日夜夜操网爽| 国产成+人综合+亚洲专区| 又紧又爽又黄一区二区| 不卡一级毛片| 成人av一区二区三区在线看| 日韩高清综合在线| 久久天躁狠狠躁夜夜2o2o| 天天躁日日操中文字幕| 亚洲最大成人av| 淫秽高清视频在线观看| 亚洲自偷自拍三级| 男女之事视频高清在线观看| www.www免费av| 在线观看午夜福利视频| 桃红色精品国产亚洲av| 1024手机看黄色片| 在线播放国产精品三级| 欧美绝顶高潮抽搐喷水| 国产精品久久久久久久电影| 亚洲精品日韩av片在线观看| 国产三级黄色录像| 热99re8久久精品国产| 日日夜夜操网爽| 人妻夜夜爽99麻豆av| www.999成人在线观看| 草草在线视频免费看| 亚洲国产欧美人成| 观看美女的网站| 亚洲精品久久国产高清桃花| 欧美zozozo另类| 麻豆国产av国片精品| 免费在线观看亚洲国产| 露出奶头的视频| 一级毛片久久久久久久久女| 听说在线观看完整版免费高清| 人人妻,人人澡人人爽秒播| ponron亚洲| 欧美+日韩+精品| 国产熟女xx| 亚洲国产日韩欧美精品在线观看| 看十八女毛片水多多多| 国产精品综合久久久久久久免费| 大型黄色视频在线免费观看| 日韩中字成人| 亚洲欧美日韩无卡精品| 国产成年人精品一区二区| av视频在线观看入口| 一级作爱视频免费观看| 亚洲成a人片在线一区二区| av黄色大香蕉| 免费av观看视频| 亚洲专区国产一区二区| 男女做爰动态图高潮gif福利片| 亚洲欧美日韩卡通动漫| 日本在线视频免费播放| 脱女人内裤的视频| 麻豆久久精品国产亚洲av| 国产蜜桃级精品一区二区三区| 一个人免费在线观看的高清视频| 国产欧美日韩一区二区精品| 天堂影院成人在线观看| 黄色日韩在线| 欧美三级亚洲精品| 深夜精品福利| 亚洲最大成人手机在线| 久久久久精品国产欧美久久久| 成人一区二区视频在线观看| 亚洲色图av天堂| 中文字幕av成人在线电影| 搡老岳熟女国产| 欧美成人性av电影在线观看| 一卡2卡三卡四卡精品乱码亚洲| 乱码一卡2卡4卡精品| 中文字幕av在线有码专区| 成人av一区二区三区在线看| 高清日韩中文字幕在线| 十八禁人妻一区二区| 99热6这里只有精品| 国产视频内射| 久久国产精品人妻蜜桃| 91狼人影院| 免费av不卡在线播放| 日韩欧美国产一区二区入口| 3wmmmm亚洲av在线观看| 久久久久久久午夜电影| 天天躁日日操中文字幕| 在线免费观看不下载黄p国产 | 看十八女毛片水多多多| 精品一区二区三区人妻视频| 一a级毛片在线观看| 88av欧美| 日韩高清综合在线| 亚洲第一电影网av| 日韩精品中文字幕看吧| www日本黄色视频网| 黄色视频,在线免费观看| 1024手机看黄色片| 午夜福利成人在线免费观看| 夜夜躁狠狠躁天天躁| 国内精品一区二区在线观看| 国产一区二区三区视频了| 午夜a级毛片| 一进一出好大好爽视频| 超碰av人人做人人爽久久| 亚洲精华国产精华精| 国产亚洲欧美98| 中文字幕久久专区| 在线a可以看的网站| 亚洲成人精品中文字幕电影| 精品久久国产蜜桃| 亚洲狠狠婷婷综合久久图片| 熟女人妻精品中文字幕| 无人区码免费观看不卡| 成人欧美大片| 人妻制服诱惑在线中文字幕| 两个人的视频大全免费| 国产精品久久久久久亚洲av鲁大| 国产黄a三级三级三级人| 特大巨黑吊av在线直播| 最近最新中文字幕大全电影3| 少妇的逼水好多| 成人高潮视频无遮挡免费网站| 国产亚洲精品av在线| 熟女电影av网| 午夜a级毛片| 亚洲无线在线观看| 国产精品嫩草影院av在线观看 | 夜夜夜夜夜久久久久| 亚洲精品日韩av片在线观看| 欧美精品国产亚洲| 国产黄色小视频在线观看| 亚洲欧美激情综合另类| 精品久久久久久成人av| 又黄又爽又免费观看的视频| 亚洲av成人不卡在线观看播放网| 变态另类丝袜制服| 美女xxoo啪啪120秒动态图 | 直男gayav资源| 一本综合久久免费| 久久香蕉精品热| 男人舔女人下体高潮全视频| 搡老妇女老女人老熟妇| 中文字幕精品亚洲无线码一区| 99热这里只有是精品50| 亚洲成人免费电影在线观看| www日本黄色视频网| 国产黄色小视频在线观看| 三级毛片av免费| 国产精品一区二区免费欧美| 欧美区成人在线视频| 亚洲片人在线观看| 少妇裸体淫交视频免费看高清| 国产亚洲欧美在线一区二区| 免费观看的影片在线观看| 能在线免费观看的黄片| 乱码一卡2卡4卡精品| 婷婷精品国产亚洲av| 国产不卡一卡二| or卡值多少钱| 熟女人妻精品中文字幕| 人妻丰满熟妇av一区二区三区| 国产精品不卡视频一区二区 | ponron亚洲| 五月玫瑰六月丁香| 欧美黑人巨大hd| 成人鲁丝片一二三区免费| netflix在线观看网站| 波野结衣二区三区在线| 婷婷六月久久综合丁香| 免费人成视频x8x8入口观看| 欧美成人免费av一区二区三区| 亚洲无线观看免费| 欧美+亚洲+日韩+国产| 嫩草影院精品99| 成人无遮挡网站| 麻豆国产av国片精品| 9191精品国产免费久久| 亚洲美女视频黄频| 男人和女人高潮做爰伦理| 人妻丰满熟妇av一区二区三区| av天堂在线播放| 久久久久国内视频| 午夜日韩欧美国产| 久久精品夜夜夜夜夜久久蜜豆| 欧美乱妇无乱码| 欧美日韩国产亚洲二区| 久久久久久久久大av| 狂野欧美白嫩少妇大欣赏| 欧美精品啪啪一区二区三区| 国产精品,欧美在线| 99riav亚洲国产免费| 亚洲精品456在线播放app | 在线观看一区二区三区| а√天堂www在线а√下载| 欧美极品一区二区三区四区| 日本熟妇午夜| 少妇丰满av| 一边摸一边抽搐一进一小说| 可以在线观看毛片的网站| 老熟妇仑乱视频hdxx| 宅男免费午夜| 日韩亚洲欧美综合| 中文字幕高清在线视频| 国产三级中文精品| 欧美日韩乱码在线| 97超级碰碰碰精品色视频在线观看| 日韩国内少妇激情av| 国产蜜桃级精品一区二区三区| 一边摸一边抽搐一进一小说| 免费看日本二区| 一级黄色大片毛片| 亚洲18禁久久av| 51午夜福利影视在线观看| 亚洲在线自拍视频| 亚洲精品一区av在线观看| 日本a在线网址| 免费黄网站久久成人精品 | 美女大奶头视频| 91久久精品国产一区二区成人| 蜜桃亚洲精品一区二区三区| 日韩中文字幕欧美一区二区| 少妇裸体淫交视频免费看高清| 久久久精品欧美日韩精品| 啦啦啦观看免费观看视频高清| 亚洲av成人不卡在线观看播放网| 51午夜福利影视在线观看| 国产综合懂色| 亚洲成人久久性| 99国产综合亚洲精品| 内地一区二区视频在线| 国产免费一级a男人的天堂| 露出奶头的视频| 成熟少妇高潮喷水视频| 精品一区二区三区视频在线观看免费| 色在线成人网| 亚洲精华国产精华精| 国产精品一区二区免费欧美| 日韩免费av在线播放| 在线观看美女被高潮喷水网站 | 欧美日韩亚洲国产一区二区在线观看| 日本黄色视频三级网站网址| 国产国拍精品亚洲av在线观看| 日韩av在线大香蕉| 国产v大片淫在线免费观看| 日本黄色视频三级网站网址| 淫秽高清视频在线观看| 校园春色视频在线观看| 午夜福利在线观看免费完整高清在 | 日韩欧美免费精品| 老女人水多毛片| 村上凉子中文字幕在线| 免费观看的影片在线观看| 国产成+人综合+亚洲专区| 无遮挡黄片免费观看| 久久午夜福利片| 国产av麻豆久久久久久久| 色噜噜av男人的天堂激情| 51国产日韩欧美| 日韩精品中文字幕看吧| 国产精品日韩av在线免费观看| 日本三级黄在线观看| 欧美一区二区精品小视频在线| 亚洲aⅴ乱码一区二区在线播放| 成人一区二区视频在线观看| 亚洲国产精品sss在线观看| 久久性视频一级片| 国产精品亚洲美女久久久| 黄片小视频在线播放| 又爽又黄无遮挡网站| 亚洲中文字幕一区二区三区有码在线看| 一级毛片久久久久久久久女| 久久久久久久亚洲中文字幕 | 深夜a级毛片| 91在线精品国自产拍蜜月| 怎么达到女性高潮| 男人舔女人下体高潮全视频| 一级作爱视频免费观看| 午夜日韩欧美国产| 亚洲综合色惰| 亚洲成人久久爱视频| 老熟妇乱子伦视频在线观看| 琪琪午夜伦伦电影理论片6080| 特大巨黑吊av在线直播| 精品人妻偷拍中文字幕| 国产高清有码在线观看视频| 亚洲人成网站高清观看| www日本黄色视频网| 三级毛片av免费| 国产视频一区二区在线看| 在线观看美女被高潮喷水网站 | 国产色婷婷99| 看十八女毛片水多多多| 国产一区二区亚洲精品在线观看| 国产成人啪精品午夜网站| 大型黄色视频在线免费观看| 最近最新免费中文字幕在线| 国产熟女xx| 深夜精品福利| 亚洲av不卡在线观看| 午夜福利高清视频| 成人午夜高清在线视频| 看黄色毛片网站| 99热只有精品国产| 精品无人区乱码1区二区| 中文亚洲av片在线观看爽| 国产精品自产拍在线观看55亚洲| 一本精品99久久精品77| 日韩欧美精品v在线| 别揉我奶头~嗯~啊~动态视频| 欧美激情久久久久久爽电影| 欧美性猛交黑人性爽| 又黄又爽又刺激的免费视频.| 久久精品夜夜夜夜夜久久蜜豆| aaaaa片日本免费| 亚洲成人免费电影在线观看| 91久久精品国产一区二区成人| 国产精品不卡视频一区二区 | 国产亚洲精品久久久com| 日韩国内少妇激情av| 欧美一区二区国产精品久久精品| 国产av在哪里看| 久久精品久久久久久噜噜老黄 | 一二三四社区在线视频社区8| 丁香六月欧美| 亚洲最大成人av| 熟妇人妻久久中文字幕3abv| 国内毛片毛片毛片毛片毛片| 亚洲,欧美,日韩| 色精品久久人妻99蜜桃| av天堂中文字幕网| 毛片女人毛片| 在线观看av片永久免费下载| xxxwww97欧美| 搡老岳熟女国产| 欧美极品一区二区三区四区| 色综合亚洲欧美另类图片| 久久国产精品人妻蜜桃| 欧美黑人欧美精品刺激| 99久久精品一区二区三区| 我要看日韩黄色一级片| 欧美精品啪啪一区二区三区| 亚洲欧美日韩卡通动漫| 亚洲精品一卡2卡三卡4卡5卡| 免费观看的影片在线观看| 亚洲专区中文字幕在线| 亚洲成a人片在线一区二区| av国产免费在线观看| 尤物成人国产欧美一区二区三区| 91麻豆精品激情在线观看国产| 免费观看精品视频网站| 亚洲狠狠婷婷综合久久图片| 欧美激情久久久久久爽电影| 亚洲欧美清纯卡通| 久久欧美精品欧美久久欧美| 特大巨黑吊av在线直播| 欧美黄色片欧美黄色片| 日韩欧美在线乱码| 久久6这里有精品| 欧美三级亚洲精品| 精品无人区乱码1区二区| 一a级毛片在线观看| 夜夜躁狠狠躁天天躁| 亚洲性夜色夜夜综合| 精品福利观看| 久久午夜亚洲精品久久| 亚洲真实伦在线观看| 精品99又大又爽又粗少妇毛片 | 国产一区二区在线av高清观看| 国产精品亚洲美女久久久| 欧美成人a在线观看| 两个人的视频大全免费| 男人狂女人下面高潮的视频| 高清日韩中文字幕在线|