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

    風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)測(cè)風(fēng)資料的對(duì)比

    2016-02-15 02:18:02曲巧娜盛春巖車(chē)軍輝
    干旱氣象 2016年6期
    關(guān)鍵詞:高度層風(fēng)廓探空

    曲巧娜,盛春巖,孫 青,車(chē)軍輝

    (1.山東省氣象科學(xué)研究所,山東 濟(jì)南 250031;2.山東省氣象局,山東 濟(jì)南 250031)

    風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)測(cè)風(fēng)資料的對(duì)比

    曲巧娜1,盛春巖1,孫 青2,車(chē)軍輝1

    (1.山東省氣象科學(xué)研究所,山東 濟(jì)南 250031;2.山東省氣象局,山東 濟(jì)南 250031)

    為了探討風(fēng)廓線(xiàn)雷達(dá)資料的可用性,對(duì)2013年9月—2015年10月青島站和濟(jì)南站的風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)測(cè)風(fēng)數(shù)據(jù)進(jìn)行相關(guān)、誤差及有效樣本比率分析。結(jié)果表明:(1)濟(jì)南站和青島站絕大多數(shù)高度層00:00和12:00風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)的水平風(fēng)速顯著正相關(guān),通過(guò)α=0.05及以上信度檢驗(yàn);(2)濟(jì)南站00:00和12:00,晴天1.5 km以上及雨天0.64 km以上大多高度層風(fēng)廓線(xiàn)雷達(dá)的水平風(fēng)速比L波段探空雷達(dá)偏小約2 m·s-1,且當(dāng)風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)水平風(fēng)向差≤20°時(shí),有效樣本比率基本在70%以上,資料質(zhì)量很高;(3)青島站00:00和12:00,6.48 km以下大多高度層風(fēng)廓線(xiàn)雷達(dá)探測(cè)的水平風(fēng)速比L波段探空雷達(dá)偏小2~4 m·s-1,水平風(fēng)速資料可用,但當(dāng)2部雷達(dá)風(fēng)向差≤20°時(shí),有效樣本比率僅為20%,海陸風(fēng)及2種儀器的布設(shè)距離是水平風(fēng)向差異的主要原因。

    風(fēng)廓線(xiàn)雷達(dá);L波段探空雷達(dá);誤差分析;有效樣本比率

    曲巧娜,盛春巖,孫 青,等.風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)測(cè)風(fēng)資料的對(duì)比[J].干旱氣象,2016,34(6):1078-1086,[QU Qiaona,SHENG Chunyan,SUN Qing,et al.Comparison of HorizontalWind Between Wind Profile Radar and L Band Sounding Radar[J].Journal of Arid Meteorology,2016,34(6):1078-1086],DOI:10.11755/j.issn.1006-7639(2016)-06-1078

    引 言

    風(fēng)廓線(xiàn)雷達(dá)是21世紀(jì)高空探測(cè)系統(tǒng)的重要組成部分,已經(jīng)在全球進(jìn)入業(yè)務(wù)化應(yīng)用階段,其不間斷實(shí)時(shí)探測(cè)站點(diǎn)上空大氣的風(fēng)向風(fēng)速、大氣湍流等,與氣球軌跡法測(cè)風(fēng)原理差別較大[1-2]。近幾年我國(guó)也相繼開(kāi)展了風(fēng)廓線(xiàn)雷達(dá)資料與探空資料的對(duì)比工作。孫旭映等[3]指出自地面至高空200m左右的范圍內(nèi),廓線(xiàn)儀探測(cè)的風(fēng)向、風(fēng)速與氣球探測(cè)值有一定偏差,而200 m以上具有很好的相關(guān)性,但在風(fēng)廓線(xiàn)雷達(dá)的最大探測(cè)高度處偏差很大;研究表明[4-5],風(fēng)廓線(xiàn)雷達(dá)水平風(fēng)速和水平風(fēng)向測(cè)量值一般都比氣球軌跡法的結(jié)果偏大,分別偏大1.5~3.0 m·s-1和10°~20°;陳添宇等[6]指出,祁連山地區(qū)風(fēng)廓線(xiàn)雷達(dá)探測(cè)的風(fēng)速誤差較大,尤其是1 400 m以下誤差更大;楊梅等[7]對(duì)探空風(fēng)廓線(xiàn)與雷達(dá)風(fēng)廓線(xiàn)資料進(jìn)行了相關(guān)分析,指出兩者具有較好的相關(guān)性。董保舉[8-9]、胡明寶[10]、鄧闖[4]等分析 了不同天氣條件下風(fēng)廓線(xiàn)雷達(dá)的數(shù)據(jù)獲取率,得出夏季數(shù)據(jù)獲取率明顯高于冬季,且降水出現(xiàn)時(shí)環(huán)境風(fēng)場(chǎng)的不均勻性造成水平風(fēng)向、風(fēng)速的測(cè)量誤差較大。

    山東省已建成的風(fēng)廓線(xiàn)雷達(dá)有6部,分別是青島、煙臺(tái)、臨沂、東營(yíng)、蓬萊、濟(jì)南,其中兼具探空氣球觀測(cè)的站點(diǎn)僅有青島站和濟(jì)南站。山東省WRF模式已業(yè)務(wù)運(yùn)行多年,并不斷致力于模式的優(yōu)化,由于風(fēng)廓線(xiàn)雷達(dá)可以提供高時(shí)間密度的垂直風(fēng)探測(cè)資料,該資料在WRF模式中的同化應(yīng)用是改進(jìn) WRF模式預(yù)報(bào)質(zhì)量的重要方面。青島站、濟(jì)南站風(fēng)廓線(xiàn)雷達(dá)分別于2006年6月和2012年9月進(jìn)入運(yùn)行階段,資料較新,資料的可靠性如何?另外,降水對(duì)風(fēng)廓線(xiàn)雷達(dá)探測(cè)有一定的影響。為此,本文將晴空和降水天氣下風(fēng)廓線(xiàn)雷達(dá)的水平風(fēng)速、風(fēng)向與L波段探空資料進(jìn)行對(duì)比,檢驗(yàn)數(shù)據(jù)的可靠性,為模式同化風(fēng)廓線(xiàn)雷達(dá)資料的質(zhì)量控制方面做準(zhǔn)備。

    1 資料與方法

    1.1 資料及處理

    選用山東青島站(120.33°E,36.06°N)、濟(jì)南站(117.55°E,36.68°N)2013年9月—2015年10月每日00:00(世界時(shí),下同)和12:00的風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)平均水平風(fēng)向、風(fēng)速資料。青島站的風(fēng)廓線(xiàn)雷達(dá)為SCRTWO-01型固定式的對(duì)流層Ⅱ型雷達(dá),2006年6月安裝在城陽(yáng)區(qū)上馬鎮(zhèn)(120.22°E,36.27°N)運(yùn)行,2013年3月移至城陽(yáng)區(qū)河套鎮(zhèn)氣象探測(cè)保障中心(120.13°E,36.23°N)運(yùn)行。其探測(cè)高度為地面至6~8 km,時(shí)間分辨率為5 min (5波束)和3 min(3波束),生成的實(shí)時(shí)數(shù)據(jù)為5 min一組;雷達(dá)的垂直分辨率,地面至960 m為120 m,960 m以上為240 m。濟(jì)南站是CLC-11-D型固定式邊界層風(fēng)廓線(xiàn)雷達(dá),平均每6 min一次探測(cè)距地面10 km垂直高度的大氣風(fēng)向、風(fēng)速和低空水汽分布狀況,生成的實(shí)時(shí)數(shù)據(jù)為6 min一組;垂直分辨率,地面至1 060m為60m,1 060~2 980m為120m,2 980 m以上為240 m。

    L波段雷達(dá)探空資料和風(fēng)廓線(xiàn)雷達(dá)資料的時(shí)間和高度并不嚴(yán)格對(duì)應(yīng)。探空氣球施放時(shí)間為11:15 和23:15,觀測(cè)時(shí)間一般約為1 h;風(fēng)廓線(xiàn)雷達(dá)資料的時(shí)間間隔為逐5 min或逐6 min,經(jīng)過(guò)處理生成逐30 min、1 h平均產(chǎn)品。為便于2種資料的對(duì)比分析,選取風(fēng)廓線(xiàn)雷達(dá)逐1 h平均產(chǎn)品數(shù)據(jù)。風(fēng)廓線(xiàn)雷達(dá)的探測(cè)高度是固定的,而L波段雷達(dá)探測(cè)高度隨機(jī)且密度較大。其中,濟(jì)南站風(fēng)廓線(xiàn)雷達(dá)探測(cè)高度為地面至7 780 m(共47層);青島站探測(cè)高度為地面至7 920 m(共37層)。因此,首選風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)相同高度的探測(cè)數(shù)據(jù),其次是當(dāng)2種雷達(dá)探測(cè)高度差的絕對(duì)值≤3 m,認(rèn)定為同一高度,而當(dāng)絕對(duì)值>3 m時(shí),則認(rèn)為誤差較大,這一高度上的2種資料同時(shí)記為缺測(cè)(實(shí)際上差值>3 m的高度較少)。

    1.2 誤差分析

    為探討風(fēng)廓線(xiàn)雷達(dá)資料的可用性,利用實(shí)際業(yè)務(wù)中廣泛應(yīng)用的L波段探空雷達(dá)資料,采用平均絕對(duì)誤差、平均誤差進(jìn)行風(fēng)廓線(xiàn)雷達(dá)資料的誤差分析,其公式如下:

    其中,MAE是平均絕對(duì)誤差,ME是平均誤差;VP是指風(fēng)廓線(xiàn)雷達(dá)水平風(fēng)速,VL是指L波段探空雷達(dá)水平風(fēng)速,單位:m·s-1;N是總樣本數(shù)。

    1.3 有效樣本比率

    另外,采用有效樣本比率[5]對(duì)水平風(fēng)速和風(fēng)向的差異性進(jìn)行分析。當(dāng)風(fēng)廓線(xiàn)雷達(dá)和L波段探空雷達(dá)的水平風(fēng)速差≤3 m·s-1(風(fēng)向差≤20°)時(shí)定義為有效樣本(Ve),2類(lèi)資料均不缺測(cè)樣本為總樣本(Vt),那么有效樣本比率Vr定義為:

    2 水平風(fēng)速差異性對(duì)比

    2.1 水平風(fēng)速的相關(guān)性

    表1給出2013年9月—2015年10月期間青島站12:00風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)探測(cè)的水平風(fēng)速在不同高度上的相關(guān)性,考慮篇幅原因,文中僅列出了部分高度層??傮w看來(lái),青島站12:00風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)的水平風(fēng)速在不同高度上均呈顯著正相關(guān),但不同天氣條件下,相關(guān)性有較大差異。其中,晴天,480 m以下2種資料的相關(guān)系數(shù)基本在0.5左右,600 m以上相關(guān)系數(shù)大多在0.6 ~0.8之間,各高度層相關(guān)系數(shù)均通過(guò)α=0.001的信度檢驗(yàn);雨天有降水時(shí),720 m以下及5 760 m以上高度層2種資料的相關(guān)性較低,其他高度層相關(guān)性顯著,相關(guān)系數(shù)達(dá)0.6以上,通過(guò)α=0.05及以上信度檢驗(yàn)??梢?jiàn),晴天情況下低層、雨天低層和高層12:00的2種資料差異較大,且晴天的相關(guān)性較雨天大。

    由青島站00:00風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)探測(cè)的水平風(fēng)速在不同高度上的相關(guān)系數(shù)可知(表略),晴天時(shí),2種資料在各高度層上均顯著正相關(guān)(通過(guò)α=0.001的信度檢驗(yàn)),600 m以下相關(guān)系數(shù)<0.7,720~6 240 m高度層的相關(guān)系數(shù)較高為0.7~0.9,其中1 680 m高度的相關(guān)系數(shù)最高,而6 480 m以上相關(guān)系數(shù)大多<0.7,可見(jiàn)中層的相關(guān)性更大;雨天有降水時(shí),840 m以下相關(guān)系數(shù)<0.8,其中120 m 和240m高度未通過(guò)顯著性檢驗(yàn),而960~5 280m相關(guān)系數(shù)大多在0.8~0.9之間波動(dòng),通過(guò)了α=0.05及以上信度檢驗(yàn),5 280 m以上相關(guān)性較低,其中5 520 m和6 480 m未通過(guò)顯著性檢驗(yàn)。綜合對(duì)比青島站00:00和12:00的資料,發(fā)現(xiàn)風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)水平風(fēng)速的相關(guān)系數(shù)中層整體好于低層、高層,00:00好于12:00,晴天高于雨天。

    從濟(jì)南站12:00風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)探測(cè)的水平風(fēng)速在不同高度上的相關(guān)性看出(表2),晴天時(shí),2種資料在各高度層上均顯著正相關(guān)(通過(guò)α=0.001的信度檢驗(yàn)),除了280 m和340 m高度上的相關(guān)系數(shù)(<0.5)較低外,其他高度層上

    相關(guān)系數(shù)較高,在0.8~0.9之間;雨天有降水時(shí),460~940 m及6 340 m以上大多高度層相關(guān)系數(shù)(≤0.7)略低,其中460 m未通過(guò)信度檢驗(yàn),其他高度層的相關(guān)系數(shù)在0.9左右,均通過(guò)α=0.05及以上信度檢驗(yàn)。濟(jì)南站00:00(表略),晴天時(shí),各高度層上2種資料均顯著正相關(guān)(通過(guò)α=0.001信度檢驗(yàn)),除760~1 900 m高度層相關(guān)系數(shù)在0.7左右以外,其他高度上相關(guān)系數(shù)大多在0.8以上;雨天有降水時(shí),除340 m和1 060 m高度未通過(guò)信度檢驗(yàn)外,700~6 100 m高度層上的相關(guān)系數(shù)在0.9左右,其他高度上相關(guān)系數(shù)≤0.8,通過(guò)α=0.05及以上信度檢驗(yàn)??梢?jiàn),濟(jì)南站風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)絕大多數(shù)高度層上水平風(fēng)速的相關(guān)系數(shù)12:00好于00:00,晴天高于雨天,中層整體好于低層和高層。

    表1 2013年9月—2015年10月12:00青島站風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)水平風(fēng)速在不同高度上的相關(guān)性Tab.1 The correlations of horizontal wind speed at different heights between wind profile radar and L band sounding radar at12:00 UTC from September 2013 to October 2015 Qingdao station of Shandong

    表2 2013年9站風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)不同高度上水平風(fēng)速的相關(guān)性Tab.2 Th speed at different heights between wind profile radar and L band soundeptember 2013 to October 2015 Ji'nan station of Shandong

    綜上所述,濟(jì)南站風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)水平風(fēng)速的相關(guān)性整體上比青島站偏高;青島站的相關(guān)系數(shù)12:00比00:00偏低,而濟(jì)南站的相關(guān)系數(shù)12:00比00:00偏高;除雨天青島站12:00有7個(gè)高度層、00:00有4個(gè)高度層及濟(jì)南站12:00有1個(gè)高度層、00:00有2個(gè)高度層未通過(guò)信度檢驗(yàn)外,其他各高度層均通過(guò)α=0.05及以上信度檢驗(yàn)。說(shuō)明青島站和濟(jì)南站的風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)的水平風(fēng)速相關(guān)性較好,風(fēng)速具體相差多少,還需進(jìn)一步分析。

    圖1 2013年9月—2015年10月青島站(a、c、e)和濟(jì)南站(b、d、f)風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)不同天氣條件下12:00和00:00的水平風(fēng)速平均絕對(duì)誤差Fig.1 Themean absolute errors of horizonalwind speed with differentweather conditions between wind profile radar and L band sounding radar at12:00 UTC and 00:00 UTC from September 2013 to October 2015 Qingdao(a,c,e)and Ji'nan(b,d,f)stations of Shandong

    2.2 水平風(fēng)速的誤差

    圖1給出2013年9月—2015年10月青島站和濟(jì)南站風(fēng)廓線(xiàn)雷達(dá)相對(duì)于L波段探空雷達(dá)12:00和00:00的水平風(fēng)速平均絕對(duì)誤差??梢钥闯?,總體上青島站0.24—5.76 km高度上平均絕對(duì)誤差有著隨高度升高先減小后增大的特征,拐點(diǎn)約在3 km附近,6.48 km以下誤差為2~4 m·s-1,且00:00的誤差整體較12:00?。? km以上誤差及其波動(dòng)明顯增大(圖1a)。晴天時(shí)(圖1c),4.32 km以下平均絕對(duì)誤差隨高度升高先減小后增大,平均絕對(duì)誤差主要分布在1.4~3.7 m·s-1之間,拐點(diǎn)約在2 km附近,且00:00比12:00誤差略偏小,1.68 km高度上00:00的誤差最?。?.4 m·s-1);5 km以上誤差波動(dòng)較大,12:00和00:00的平均絕對(duì)誤差在3.0 m·s-1左右波動(dòng),且00:00比12:00誤差略偏大。雨天有降水時(shí)(圖1e),12:00比00:00誤差明顯偏大,6 km以下水平風(fēng)速的平均絕對(duì)誤差在1~5 m·s-1之間,6 km以上平均絕對(duì)誤差為1~8 m·s-1??梢?jiàn),無(wú)論是12:00還是00:00,風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)的水平風(fēng)速降水天氣比晴天的平均絕對(duì)誤差偏大,且0.6 km以下和4.32 km以上高度層誤差相對(duì)較大。另外,晴天,高層12:00的誤差比00:00小,低層正相反,12:00的誤差比00:00大;降水天氣,各高度層12:00均比00:00誤差明顯偏大。

    濟(jì)南站,1 km以上總體樣本誤差隨著高度先增大后減小再增大,誤差主要在2 m·s-1左右波動(dòng)(圖1b)。晴天時(shí)(圖1d),12:00水平風(fēng)速的平均絕對(duì)誤差從地面至 0.64 km先減小,0.64—2 km增大,2—4 km再次減小,4 km以上再次增大;00:00平均絕對(duì)誤差也表現(xiàn)出類(lèi)似的變化特征,只是拐點(diǎn)高度較12:00偏低。另外,0—1.9 km和3.7 km高度以上12:00的平均絕對(duì)誤差整體比00:00偏小,而2—3.7 km高度上12:00比00:00偏大。雨天有降水時(shí)(圖1f),5 km以下誤差隨高度升高變化較小,5 km以上隨著高度升高增長(zhǎng)較快,除少數(shù)高度層外,12:00平均絕對(duì)誤差均比00:00偏大。

    從青島站總體樣本平均誤差來(lái)看,7.44 km以下大多高度層風(fēng)廓線(xiàn)雷達(dá)測(cè)風(fēng)比L波段探空雷達(dá)水平風(fēng)速偏小,00:00平均誤差較12:00偏?。▓D2a)。晴天時(shí)(圖2c),2個(gè)時(shí)次的平均誤差均表現(xiàn)出先減小后增大再減小的特征,00:00平均誤差較12:00整體偏小。其中0.48—7.44 km風(fēng)廓線(xiàn)雷達(dá)測(cè)風(fēng)均比L波段探空雷達(dá)水平風(fēng)速偏小,平均誤差為0~3 m·s-1。雨天有降水時(shí)(圖 2e),6.72 km以下,除2—3 km高度上12:00及3.84 km高度上00:00風(fēng)廓線(xiàn)雷達(dá)探測(cè)風(fēng)速比L波段探空雷達(dá)風(fēng)速偏大外,其他時(shí)次高度上風(fēng)廓線(xiàn)雷達(dá)測(cè)風(fēng)均比L波段探空雷達(dá)風(fēng)速偏小,平均誤差為0~4 m·s-1,且表現(xiàn)出先減后增的變化特征,00:00平均誤差較12:00明顯偏小,且波動(dòng)也較小。

    圖2 2013年9月—2015年10月青島站(a、c、e)和濟(jì)南站(b、d、f)風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)不同天氣條件下12:00和00:00的水平風(fēng)速平均誤差Fig.2 Themean errors of horizonalwind speed with differentweather conditions between wind profile radar and L band sounding radar at12:00 UTC and 00:00 UTC from September 2013 to October 2015 in Qingdao(a,c,e)and Ji'nan(b,d,f)stations of Shandong

    從濟(jì)南站風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)12:00 和00:00的總體樣本水平風(fēng)速平均誤差看出(圖2b),平均誤差隨著高度升高逐漸減小,從地面至7.06 km平均誤差由1.5 m·s-1減小至-2.6 m·s-1,且12:00的平均誤差較00:00略小。其中,1 km以下2個(gè)時(shí)刻風(fēng)廓線(xiàn)雷達(dá)測(cè)風(fēng)均比L波段探空雷達(dá)水平風(fēng)速偏大,1 km以上風(fēng)廓線(xiàn)雷達(dá)測(cè)風(fēng)均比L波段探空雷達(dá)水平風(fēng)速偏小。晴天時(shí)(圖2d),2個(gè)時(shí)刻的水平風(fēng)速平均誤差變化特征與總體樣本一致。其中1.42 km以下2個(gè)時(shí)刻風(fēng)廓線(xiàn)雷達(dá)測(cè)風(fēng)比L波段探空雷達(dá)水平風(fēng)速偏大,平均誤差為1 m·s-1左右,1.5 km以上風(fēng)廓線(xiàn)雷達(dá)測(cè)風(fēng)比L波段探空雷達(dá)水平風(fēng)速偏小。雨天有降水時(shí)(圖2f),0.64 km以下2個(gè)時(shí)刻風(fēng)廓線(xiàn)雷達(dá)比L波段探空雷達(dá)水平風(fēng)速偏大,0.64 km以上,除00:00的1.3 km和2.74 km個(gè)別高度層風(fēng)廓線(xiàn)雷達(dá)比L波段探空雷達(dá)水平風(fēng)速偏大外,其他時(shí)次高度層風(fēng)廓線(xiàn)雷達(dá)比L波段探空雷達(dá)水平風(fēng)速偏小。

    2.3 水平風(fēng)速差的有效樣本比率

    對(duì)青島站和濟(jì)南站風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)12:00和00:00水平風(fēng)速進(jìn)行有效樣本比率計(jì)算。從青島站12:00和00:00總體有效樣本比率隨高度的變化來(lái)看,6.48 km高度以下有效樣本比率隨著高度升高先增大后減小,且2 km以下有效樣本比率00:00高于12:00,5.5 km以上正相反,12:00明顯高于00:00(圖3a)。晴天時(shí)(圖3c),2個(gè)時(shí)刻的有效樣本比率也表現(xiàn)出隨高度先增后減的變化特征,0.72 km以下12:00有效樣本比率為50% ~70%,00:00有效樣本比率約為70%,00:00略高于12:00;0.72—4 km比率較高,最高達(dá)91%;6 km以上00:00有效樣本比率為55% ~78%,12:00有效樣本比率為68% ~88%,12:00明顯高于00:00。雨天有降水時(shí)(圖3e),有效樣本比率不穩(wěn)定,整層高度在40%~90%之間浮動(dòng),12:00明顯較00:00波動(dòng)大。

    圖3 2013年9月—2015年10月青島站(a、c、e)、濟(jì)南站(b、d、f)不同天氣條件下12:00和00:00水平風(fēng)速有效樣本比率隨高度的變化Fig.3 The effective sample radios of horizonalwind speed with differentweather conditions at12:00 UTC and 00:00 UTC from September 2013 to October 2015 in Qingdao(a,c,e)and Ji'nan(b,d,f)stations of Shandong

    濟(jì)南站2個(gè)時(shí)刻的總體水平風(fēng)速有效樣本比率隨著高度升高先減小后增大再減小的變化特征,且2個(gè)時(shí)刻相差不大(圖3b)。晴天時(shí)(圖3d),2個(gè)時(shí)刻的水平風(fēng)速有效樣本比率也表現(xiàn)出隨高度先減小后增大再減小的變化特征。其中,1.5 km以下12:00(80% ~90%)有效樣本比率略高于00:00 (70% ~80%),1.66—4.42 km有效樣本比率增大至90%以上,12:00略低于00:00,4.66 km以上有效樣本比率逐漸減小到70%左右,12:00略高于00:00。雨天有降水時(shí)(圖3f),有效樣本比率不穩(wěn)定,整層高度在 60%~100%之間浮動(dòng),2 km以下00:00普遍高于12:00。綜上所述,當(dāng)風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)水平風(fēng)速差≤3 m·s-1時(shí),濟(jì)南站的有效樣本比率好于青島站。

    3 水平風(fēng)向差異性對(duì)比

    由青島站12:00和00:00水平風(fēng)向有效樣本比率隨高度變化曲線(xiàn)可看出(圖4a),總體水平風(fēng)向有效樣本比率較低,約為 20%~40%,且 00:00有效樣本比率較12:00偏大。晴天時(shí)(圖4c),12:00有效樣本比率較為穩(wěn)定,約為10%~20%;00:00有效樣本比率略高于12:00,且4.32 km以上隨著高度升高逐漸增大,由12%逐漸增大至42%以上。雨天有降水時(shí)(圖4e),有效樣本比率在0%~40%范圍內(nèi)波動(dòng)??梢?jiàn),青島站風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)水平風(fēng)向相差非常大。

    圖4 2013年9月—2015年10月青島站(a、c、e)、濟(jì)南站(b、d、f)不同天氣條件下12:00和00:00水平風(fēng)向有效樣本比率隨高度的變化Fig.4 The effective sample radios of horizonalwind direction with differentweather conditions at12:00 UTC and 00:00 UTC from September 2013 to October 2015 in Qingdao(a,c,e)and Ji'nan(b,d,f)stations of Shandong

    由濟(jì)南站12:00和00:00水平風(fēng)向有效樣本比率隨高度變化曲線(xiàn)看出(圖4b),總體水平風(fēng)向有效樣本比率基本在60%以上,且隨著高度升高比率逐漸增大,00:00略高于12:00。晴天時(shí)(圖4d),3.22 km以下,水平風(fēng)向有效樣本比率約為70%,3.22 km以上有效樣本比率增到80% ~90%,且隨著高度升高比率增大,00:00的有效樣本比率略高于12:00。雨天有降水時(shí)(圖4f),有效樣本比率約在65%~100%之間浮動(dòng),同樣00:00的有效樣本比率略高于12:00,且波動(dòng)也大。綜上所述,當(dāng)風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)水平風(fēng)向差≤20°時(shí),濟(jì)南站有效樣本比率基本在70%以上,2種雷達(dá)探測(cè)的水平風(fēng)向相差較小,資料的可用性較高;而青島站有效樣本比率僅為20%,2種雷達(dá)探測(cè)的水平風(fēng)向相差非常大。

    4 誤差原因

    通過(guò)對(duì)比濟(jì)南站和青島站風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)的水平風(fēng)速、風(fēng)向發(fā)現(xiàn),2站的誤差及有效樣本比率結(jié)果明顯不同。除個(gè)別高度層外,濟(jì)南站00:00和12:00晴天1.5 km以上及雨天0.64 km以上風(fēng)廓線(xiàn)雷達(dá)比L波段探空雷達(dá)的水平風(fēng)速偏小2 m·s-1左右,其原因主要是2種儀器測(cè)風(fēng)原理不同[1-2,11],探空氣球在上升的過(guò)程中受到氣流的影響發(fā)生漂移,而風(fēng)廓線(xiàn)雷達(dá)則依靠探測(cè)大氣中的湍流引起的大氣折射指數(shù)變化來(lái)獲取高空風(fēng)數(shù)據(jù)。

    圖5 濟(jì)南站和青島站風(fēng)廓線(xiàn)雷達(dá)和L波段探空雷達(dá)的位置Fig.5 The location of wind profile radar and L band sounding radar in Ji'nan and Qingdao stations of Shandong

    值得注意的是,當(dāng)風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)水平風(fēng)向差≤20°時(shí),青島站有效樣本比率僅為20%,誤差較大,究其原因主要有:一是青島站2種儀器沒(méi)有布設(shè)在同一位置上,直線(xiàn)距離約為28 km,而濟(jì)南站2種儀器是布設(shè)在同一位置(圖5);二是海陸風(fēng)[12-15],海陸風(fēng)是由于海陸熱力差異引起的白天由海上向內(nèi)陸吹、夜間由內(nèi)陸向海上吹的局地小尺度熱力環(huán)流。青島風(fēng)廓線(xiàn)雷達(dá)站位于膠州灣沿岸,而L波段探空雷達(dá)站位于黃海沿岸(圖5),當(dāng)有海陸風(fēng)發(fā)生時(shí),風(fēng)向明顯不同。

    5 結(jié) 論

    (1)除個(gè)別高度外,2站12:00和00:00風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)的水平風(fēng)速高度相關(guān),且濟(jì)南站的相關(guān)系數(shù)整體比青島站偏高。相關(guān)系數(shù)青島站12:00比00:00偏低,濟(jì)南站12:00比00:00偏高。

    (2)濟(jì)南站風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)的水平風(fēng)速平均絕對(duì)誤差比青島站偏小。青島站6.48 km以下大多高度層風(fēng)廓線(xiàn)雷達(dá)比L波段探空雷達(dá)水平風(fēng)速偏小2~4 m·s-1,且雨天比晴天的平均絕對(duì)誤差偏大。濟(jì)南站晴天1.5 km以上及雨天0.64 km以上大多高度層風(fēng)廓線(xiàn)雷達(dá)的水平風(fēng)速比L波段探空雷達(dá)偏小2 m·s-1左右。

    (3)當(dāng)風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)水平風(fēng)速差≤3 m·s-1時(shí),濟(jì)南站的有效樣本比率高于青島站。晴天時(shí),青島站水平風(fēng)速有效樣本比率隨高度升高先增大后減小,大多高度層比率值為65%~90%;濟(jì)南站有效樣本比率隨高度升高先減小后增大再減小,大多高度層比率值為70%~90%。雨天有降水時(shí),2站有效樣本比率均不穩(wěn)定,青島站在40%~90%之間擺動(dòng),濟(jì)南站在60%~100%之間浮動(dòng)。

    (4)當(dāng)風(fēng)廓線(xiàn)雷達(dá)與L波段探空雷達(dá)水平風(fēng)向差≤20°時(shí),濟(jì)南站有效樣本比率基本在70%以上,2種儀器探測(cè)的水平風(fēng)向相差較小,資料可用性較高;青島站有效樣本比率僅為20%,2種儀器探測(cè)的水平風(fēng)向相差非常大,除了2種儀器的地理位置有一定的距離外,海陸風(fēng)也是引起水平風(fēng)向差異性較大的主要原因。

    [1]楊馨蕊,馬舒慶,吳蕾.UHF風(fēng)廓線(xiàn)雷達(dá)降水?dāng)?shù)據(jù)判別方法的比較與評(píng)價(jià)[J].大氣科學(xué)學(xué)報(bào),2010,33(5):576-581.

    [2]王燁芳,薛云朝,李清.氣球軌跡法測(cè)風(fēng)和風(fēng)廓線(xiàn)雷達(dá)測(cè)風(fēng)的對(duì)比分析[J].氣象水文海洋儀器,2007,(3):30-34.

    [3]孫旭映,韓暉,段海霞,等.風(fēng)廓線(xiàn)儀與氣球測(cè)風(fēng)資料的對(duì)比分析[J].干旱氣象,2008,26(3):48-52.

    [4]鄧闖,阮征,魏鳴,等.風(fēng)廓線(xiàn)雷達(dá)測(cè)風(fēng)精度評(píng)估[J].應(yīng)用氣象學(xué)報(bào),2012,23(5):523-533.

    [5]萬(wàn)蓉,周志敏,崔春光,等.風(fēng)廓線(xiàn)雷達(dá)資料與探空資料的對(duì)比分析[J].暴雨災(zāi)害,2011,30(2):130-136.

    [6]陳添宇,陳躍,陳乾.風(fēng)廓線(xiàn)雷達(dá)資料的誤差及對(duì)祁連山地形云風(fēng)場(chǎng)檢測(cè)的初步分析[J].干旱氣象,2011,29(4):416-422.

    [7]楊梅,李玉林,單九生,等.探空風(fēng)廓線(xiàn)與雷達(dá)風(fēng)廓線(xiàn)資料相關(guān)分析[J].氣象,2006,32(6):20-24.

    [8]董保舉,張曄,徐安倫.高原地區(qū)風(fēng)廓線(xiàn)雷達(dá)資料評(píng)估[J].氣象科技,2009,37(5):580-583.

    [9]董保舉,付志嘉,李明,等.風(fēng)廓線(xiàn)雷達(dá)資料在暴雨天氣過(guò)程特征分析中的應(yīng)用[J].氣象科技,2012,40(1):74-86.

    [10]胡明寶,鄭國(guó)光,肖文建.風(fēng)廓線(xiàn)雷達(dá)數(shù)據(jù)獲取率的統(tǒng)計(jì)分析[J].現(xiàn)代雷達(dá),2008,30(10):14-16.

    [11]宋呈文,張晨飛,劉書(shū)友,等.風(fēng)廓線(xiàn)雷達(dá)探測(cè)資料分析[J].氣象水文海洋儀器,2011,(1):1-10.

    [12]吳增茂,龍寶森.青島局地風(fēng)特征的分析[J].海洋湖沼通報(bào),1993,(1):16-22.

    [13]盛春巖,王建林,刁秀廣.2006年8月青島國(guó)際帆船賽期間海陸風(fēng)特征及三維結(jié)構(gòu)分析[J].中國(guó)海洋大學(xué)學(xué)報(bào),2007,37 (4):609-614.

    [14]薛德強(qiáng),鄭全嶺,錢(qián)喜鎮(zhèn),等.山東半島的海陸風(fēng)環(huán)流及其影響[J].南京氣象學(xué)院學(xué)報(bào),1995,18(2):293-299.

    [15]邱曉暖,范紹佳.海陸風(fēng)研究進(jìn)展與我國(guó)沿海三地海陸風(fēng)主要特征[J].氣象,2013,39(2):186-193.

    Com parison of HorizontalW ind Between W ind Profile Radar and L Band Sounding Radar

    QU Qiaona1,SHENG Chunyan1,SUN Qing2,CHE Junhui1

    (1.Shandong Provincial Meteorological Science Institute,Ji'nan 250031,China;2.Shandong Provincial Meteorological Bureau,Ji'nan 250031,China)

    In order to explore the data availability of wind profile radar,the correlation,error and effective sample ratio of horizontal wind speed and direction between wind profile radar and L band sounding radar Ji'nan and Qingdao stations of Shandong Province from September 2013 to October 2015 were analyzed.The results are as follows:(1)The positive correlations of horizontal wind speed at most altitudes between wind profile radar and L band sounding radar Qingdao and Ji'nan stations at00:00 UTC and 12:00 UTC from September 2013 to October 2015 were highly significant at the 0.05 and above significance level.(2)The horizontal wind speeds of wind profile radar atmost heights above 1.5 km on sunny weather and 0.64 km on rainy weatherwere about2m·s-1less than thatof L band souding radar at00:00 UTC and 12:00 UTC from September 2013 to October2015 Ji'nan station,and when the differences of horizontalwind direction between wind profile radar and L band sounding radarwere less than or equal to 20°,the effective sample ratioswere above 70%,which indicated that the data quality of wind profile radar Ji'nan station was very well.(3)The horizonal wind speeds ofwind profile radarwere 2-4m·s-1smaller than that of L band souding radar atmostaltitudes below 6.48 km Qingdao station at00:00 UTC and 12:00 UTC from September 2013 to October 2015,the data of wind profile radar were useful.However when the differences of horizontalwind direction between wind profile radar and L band sounding radarwere less than or equal to20°,the effective sample ratioswere only 20%due to the sea-land breeze and the longer distance between wind profile radar and L band sounding radar Qingdao,the availability of wind direction from wind profile radar wasworse.

    wind profile radar;L band sounding radar;error analysis;effective sample ratio

    1006-7639(2016)-06-1078-09

    10.11755/j.issn.1006-7639(2016)-06-1078

    P412.16

    A

    2016-01-07;改回日期:2016-04-06

    公益性行業(yè)(氣象)科研專(zhuān)項(xiàng)(GYHY201106006)及山東省氣象局青年科研基金項(xiàng)目“風(fēng)廓線(xiàn)雷達(dá)資料可用性分析”(2014SDQN03)共同資助

    曲巧娜(1984-),女,碩士,工程師,主要從事天氣及模式檢驗(yàn)工作.E-mail:quqiaona77@163.com

    猜你喜歡
    高度層風(fēng)廓探空
    探空氣球?yàn)楦呖諝庀蟆鞍衙}”
    高郵邊界層風(fēng)廓線(xiàn)雷達(dá)數(shù)據(jù)獲取率分析
    一種綜合的風(fēng)廓線(xiàn)雷達(dá)數(shù)據(jù)質(zhì)量控制方法
    用L波段探空測(cè)風(fēng)雷達(dá)評(píng)估風(fēng)廓線(xiàn)雷達(dá)測(cè)風(fēng)準(zhǔn)確性
    四川盆地風(fēng)廓線(xiàn)雷達(dá)大氣折射率結(jié)構(gòu)常數(shù)特征分析
    基于高度層的航路短時(shí)利用率模型研究
    TK-2GPS人影火箭探空數(shù)據(jù)與L波段探空數(shù)據(jù)對(duì)比分析
    飛機(jī)最佳航路爬升時(shí)機(jī)研究
    淺談凈舉力對(duì)探空氣球升速及施放高度的影響
    河南科技(2014年14期)2014-02-27 14:12:08
    鄭州探空數(shù)據(jù)庫(kù)的建設(shè)簡(jiǎn)介
    河南科技(2014年4期)2014-02-27 14:06:59
    18禁裸乳无遮挡免费网站照片| 成人国产麻豆网| 亚洲国产日韩欧美精品在线观看| 久久久久精品国产欧美久久久| 中文字幕精品亚洲无线码一区| 国产一区二区激情短视频| 99久久中文字幕三级久久日本| 不卡一级毛片| 精品久久久久久久人妻蜜臀av| 韩国av在线不卡| 久久6这里有精品| 99久久中文字幕三级久久日本| 国产午夜精品久久久久久一区二区三区 | 婷婷色综合大香蕉| 日本成人三级电影网站| 免费av观看视频| 久久久久精品国产欧美久久久| 国产中年淑女户外野战色| 99riav亚洲国产免费| 一夜夜www| 十八禁网站免费在线| 少妇人妻一区二区三区视频| 日韩中文字幕欧美一区二区| 国产精品亚洲美女久久久| 淫妇啪啪啪对白视频| 亚洲精品乱码久久久v下载方式| 好男人在线观看高清免费视频| 一个人免费在线观看电影| 亚洲国产精品sss在线观看| 又黄又爽又免费观看的视频| 国产在视频线在精品| 人妻少妇偷人精品九色| 乱人视频在线观看| 国产高清视频在线观看网站| 99国产极品粉嫩在线观看| 亚洲电影在线观看av| 国产一级毛片七仙女欲春2| 搞女人的毛片| 亚洲三级黄色毛片| 精品国内亚洲2022精品成人| 啦啦啦韩国在线观看视频| 蜜桃久久精品国产亚洲av| 国产久久久一区二区三区| 精品99又大又爽又粗少妇毛片 | 色综合色国产| 内射极品少妇av片p| 亚洲国产精品久久男人天堂| 国产亚洲欧美98| 日本撒尿小便嘘嘘汇集6| 男女视频在线观看网站免费| 舔av片在线| 香蕉av资源在线| 久久精品国产亚洲av香蕉五月| 一级黄色大片毛片| 日韩欧美国产在线观看| 一级黄片播放器| 日韩欧美精品免费久久| 观看美女的网站| 99在线人妻在线中文字幕| 美女大奶头视频| 1000部很黄的大片| 亚洲无线观看免费| 国产色爽女视频免费观看| 一进一出抽搐动态| а√天堂www在线а√下载| 亚洲av成人av| 深夜a级毛片| 亚洲电影在线观看av| 国产精品亚洲美女久久久| 亚洲人成网站高清观看| 免费看av在线观看网站| 国产视频内射| 亚洲欧美清纯卡通| 免费观看人在逋| 国产黄a三级三级三级人| 男女下面进入的视频免费午夜| 亚洲精品亚洲一区二区| 麻豆成人av在线观看| 国产视频一区二区在线看| 中文资源天堂在线| 婷婷丁香在线五月| 亚洲人与动物交配视频| 色精品久久人妻99蜜桃| 制服丝袜大香蕉在线| 免费高清视频大片| www日本黄色视频网| 免费在线观看影片大全网站| 黄色一级大片看看| 色5月婷婷丁香| 十八禁国产超污无遮挡网站| 亚洲天堂国产精品一区在线| 欧美日韩亚洲国产一区二区在线观看| 亚洲成a人片在线一区二区| 国产日本99.免费观看| 女人被狂操c到高潮| 日日啪夜夜撸| 国产精品一区二区性色av| 小说图片视频综合网站| 亚洲av中文字字幕乱码综合| 啦啦啦韩国在线观看视频| 欧美最黄视频在线播放免费| 久久久久久久亚洲中文字幕| 又黄又爽又刺激的免费视频.| 亚洲图色成人| 国产aⅴ精品一区二区三区波| 国产精品久久久久久久久免| 国产麻豆成人av免费视频| 国产欧美日韩一区二区精品| 久久久色成人| 午夜精品久久久久久毛片777| 一本一本综合久久| 久久久色成人| 久久午夜福利片| 久久6这里有精品| 男女之事视频高清在线观看| 亚洲精品国产成人久久av| 亚洲精品粉嫩美女一区| 国产精品免费一区二区三区在线| 日韩精品有码人妻一区| 久久人人精品亚洲av| 久久久久久九九精品二区国产| 男女啪啪激烈高潮av片| 午夜福利18| 国产三级在线视频| 国产精品无大码| 精品久久久久久,| 3wmmmm亚洲av在线观看| 日韩欧美精品v在线| 一级a爱片免费观看的视频| 久久国产乱子免费精品| 亚洲成人中文字幕在线播放| 最新中文字幕久久久久| 日韩欧美精品v在线| 男女做爰动态图高潮gif福利片| 久久中文看片网| 国产欧美日韩精品亚洲av| 亚洲无线在线观看| 欧美一区二区亚洲| 欧美zozozo另类| 少妇裸体淫交视频免费看高清| 亚洲狠狠婷婷综合久久图片| 麻豆国产97在线/欧美| 亚洲熟妇熟女久久| 亚洲国产欧洲综合997久久,| 日日撸夜夜添| 久久久久久久久大av| 国内揄拍国产精品人妻在线| eeuss影院久久| 国产午夜精品久久久久久一区二区三区 | 男人舔女人下体高潮全视频| 国产 一区精品| 日韩欧美精品v在线| 人人妻人人澡欧美一区二区| 亚洲中文字幕一区二区三区有码在线看| 白带黄色成豆腐渣| 亚洲最大成人中文| 久久久久久久久久成人| 亚洲中文日韩欧美视频| 欧美潮喷喷水| av在线老鸭窝| 久久精品国产亚洲av天美| 最新在线观看一区二区三区| 三级毛片av免费| 一进一出抽搐动态| 最近在线观看免费完整版| 亚洲国产色片| 三级毛片av免费| 老司机福利观看| 日日摸夜夜添夜夜添小说| 精品久久久久久成人av| 蜜桃久久精品国产亚洲av| 国产真实乱freesex| 国产亚洲欧美98| 亚洲avbb在线观看| 日日夜夜操网爽| 亚洲欧美日韩高清在线视频| 亚洲在线观看片| 免费一级毛片在线播放高清视频| 99久久精品一区二区三区| 国产亚洲精品综合一区在线观看| 一个人看视频在线观看www免费| 中文亚洲av片在线观看爽| 成人毛片a级毛片在线播放| 欧美一区二区精品小视频在线| 日本免费一区二区三区高清不卡| 日韩人妻高清精品专区| 国产精品久久久久久亚洲av鲁大| 婷婷色综合大香蕉| 免费高清视频大片| 最近中文字幕高清免费大全6 | 一边摸一边抽搐一进一小说| 久久国内精品自在自线图片| 成人鲁丝片一二三区免费| 久久精品国产亚洲网站| 啦啦啦啦在线视频资源| 亚洲人成网站在线播| 精品99又大又爽又粗少妇毛片 | 色在线成人网| 久久婷婷人人爽人人干人人爱| 日韩欧美三级三区| 精品久久久久久久久av| 一个人免费在线观看电影| 久久亚洲真实| 日韩高清综合在线| 久久国内精品自在自线图片| 最近在线观看免费完整版| 热99re8久久精品国产| 97超级碰碰碰精品色视频在线观看| 干丝袜人妻中文字幕| 久久人妻av系列| 国产伦精品一区二区三区四那| 亚洲无线观看免费| 又紧又爽又黄一区二区| 久久欧美精品欧美久久欧美| 成人三级黄色视频| 久久精品国产清高在天天线| 国产精品日韩av在线免费观看| АⅤ资源中文在线天堂| 综合色av麻豆| 欧美zozozo另类| 人妻少妇偷人精品九色| 国产精品伦人一区二区| 亚洲天堂国产精品一区在线| 国产在线男女| 久久香蕉精品热| 他把我摸到了高潮在线观看| 少妇的逼好多水| 少妇猛男粗大的猛烈进出视频 | 99热这里只有是精品在线观看| 九九爱精品视频在线观看| 久久久久免费精品人妻一区二区| 97人妻精品一区二区三区麻豆| 深爱激情五月婷婷| 别揉我奶头~嗯~啊~动态视频| 国产精品永久免费网站| a级一级毛片免费在线观看| 日本成人三级电影网站| 国产精品无大码| 少妇裸体淫交视频免费看高清| 亚洲成人久久爱视频| 欧洲精品卡2卡3卡4卡5卡区| 99精品在免费线老司机午夜| 婷婷精品国产亚洲av| 欧美日本亚洲视频在线播放| 亚洲18禁久久av| 成人无遮挡网站| 亚洲avbb在线观看| 嫩草影院新地址| 欧美性猛交╳xxx乱大交人| 国产美女午夜福利| 一个人免费在线观看电影| a级毛片a级免费在线| 亚洲黑人精品在线| 在线观看午夜福利视频| 亚洲七黄色美女视频| 他把我摸到了高潮在线观看| 91久久精品国产一区二区三区| 欧美激情在线99| 亚洲欧美日韩高清专用| 成年免费大片在线观看| 女的被弄到高潮叫床怎么办 | 韩国av一区二区三区四区| 看免费成人av毛片| 国产精品亚洲一级av第二区| av.在线天堂| 欧美日韩精品成人综合77777| 深夜精品福利| 亚洲aⅴ乱码一区二区在线播放| 成人性生交大片免费视频hd| 精品日产1卡2卡| 小蜜桃在线观看免费完整版高清| 国产高清不卡午夜福利| 午夜福利在线在线| 国产主播在线观看一区二区| 亚洲国产欧洲综合997久久,| 国产一区二区亚洲精品在线观看| 亚洲精品国产成人久久av| 超碰av人人做人人爽久久| 欧美在线一区亚洲| 欧美成人一区二区免费高清观看| 国产精华一区二区三区| 亚洲人与动物交配视频| 在现免费观看毛片| 亚洲av不卡在线观看| 国产高清三级在线| 中文字幕高清在线视频| 日韩,欧美,国产一区二区三区 | 免费观看精品视频网站| 日韩中文字幕欧美一区二区| 久久久成人免费电影| 日本精品一区二区三区蜜桃| 久久久久久大精品| 精品免费久久久久久久清纯| 啪啪无遮挡十八禁网站| 国产亚洲精品久久久com| 可以在线观看的亚洲视频| 男人狂女人下面高潮的视频| 国产精品永久免费网站| 免费人成在线观看视频色| netflix在线观看网站| 成年版毛片免费区| 国内精品一区二区在线观看| 搡老熟女国产l中国老女人| 国内精品一区二区在线观看| 欧美日韩综合久久久久久 | av在线天堂中文字幕| 亚洲av第一区精品v没综合| 久久久久国内视频| 一进一出抽搐动态| 91麻豆av在线| av黄色大香蕉| 成人亚洲精品av一区二区| 乱码一卡2卡4卡精品| 赤兔流量卡办理| 国产精品99久久久久久久久| 国产精品不卡视频一区二区| 亚洲欧美日韩东京热| 亚洲国产高清在线一区二区三| 国内少妇人妻偷人精品xxx网站| 国产男靠女视频免费网站| ponron亚洲| 日本色播在线视频| 国产蜜桃级精品一区二区三区| 床上黄色一级片| 两个人视频免费观看高清| 美女cb高潮喷水在线观看| 亚洲人成伊人成综合网2020| 亚洲图色成人| 看片在线看免费视频| 国语自产精品视频在线第100页| 悠悠久久av| 亚洲精品一区av在线观看| 一进一出好大好爽视频| 久久久久久国产a免费观看| 日本免费a在线| 久久久久九九精品影院| 成人性生交大片免费视频hd| 欧美丝袜亚洲另类 | 亚洲av中文字字幕乱码综合| 国产探花极品一区二区| 亚洲中文字幕日韩| 亚洲真实伦在线观看| 精品一区二区免费观看| 国产精品日韩av在线免费观看| 日日干狠狠操夜夜爽| 嫁个100分男人电影在线观看| 中国美女看黄片| 国产人妻一区二区三区在| 免费黄网站久久成人精品| a级一级毛片免费在线观看| 婷婷亚洲欧美| 永久网站在线| 免费看a级黄色片| 亚洲av.av天堂| 制服丝袜大香蕉在线| 在现免费观看毛片| 成年免费大片在线观看| 国产久久久一区二区三区| 美女被艹到高潮喷水动态| 欧美绝顶高潮抽搐喷水| 免费在线观看影片大全网站| 亚洲精品日韩av片在线观看| 亚洲一级一片aⅴ在线观看| 欧美性猛交╳xxx乱大交人| 久久久久九九精品影院| 久久久国产成人精品二区| 中文在线观看免费www的网站| 免费av观看视频| 精品不卡国产一区二区三区| 九九久久精品国产亚洲av麻豆| 日本a在线网址| 欧美xxxx性猛交bbbb| 啦啦啦韩国在线观看视频| 精品福利观看| 午夜福利高清视频| 美女黄网站色视频| 夜夜爽天天搞| 老司机午夜福利在线观看视频| 波多野结衣巨乳人妻| 中文亚洲av片在线观看爽| 2021天堂中文幕一二区在线观| 亚洲av免费在线观看| av视频在线观看入口| 国产精品久久久久久久久免| 午夜日韩欧美国产| 麻豆久久精品国产亚洲av| 精品一区二区免费观看| 天天躁日日操中文字幕| 高清在线国产一区| 麻豆成人av在线观看| 日韩欧美精品v在线| 欧美日韩中文字幕国产精品一区二区三区| 天天躁日日操中文字幕| 欧美色视频一区免费| 成年免费大片在线观看| 91麻豆av在线| 国产伦在线观看视频一区| 男女啪啪激烈高潮av片| 色综合站精品国产| 亚洲最大成人手机在线| 国产精品永久免费网站| 伦理电影大哥的女人| 97人妻精品一区二区三区麻豆| 听说在线观看完整版免费高清| 熟女人妻精品中文字幕| 亚洲av成人精品一区久久| 婷婷精品国产亚洲av| 久久久成人免费电影| 日韩高清综合在线| 一卡2卡三卡四卡精品乱码亚洲| 少妇熟女aⅴ在线视频| 国产乱人视频| 精品久久久久久久久av| 春色校园在线视频观看| 亚洲人与动物交配视频| 看片在线看免费视频| 身体一侧抽搐| 99热只有精品国产| 性色avwww在线观看| 亚洲久久久久久中文字幕| 永久网站在线| 成人国产麻豆网| 一个人观看的视频www高清免费观看| 又爽又黄无遮挡网站| 韩国av一区二区三区四区| 亚洲av中文av极速乱 | 校园人妻丝袜中文字幕| 欧美中文日本在线观看视频| 搞女人的毛片| 日韩中字成人| 久久6这里有精品| 男人狂女人下面高潮的视频| 免费黄网站久久成人精品| 伦理电影大哥的女人| 在线观看一区二区三区| 男插女下体视频免费在线播放| 国产探花在线观看一区二区| 日本免费一区二区三区高清不卡| 99久久精品热视频| 观看美女的网站| 老司机午夜福利在线观看视频| 看免费成人av毛片| 男人舔女人下体高潮全视频| 亚洲aⅴ乱码一区二区在线播放| 免费搜索国产男女视频| 日本一二三区视频观看| 国产精品久久电影中文字幕| 人妻制服诱惑在线中文字幕| 国产一区二区三区av在线 | 精品乱码久久久久久99久播| 亚洲精华国产精华精| 国产在视频线在精品| 一区二区三区激情视频| АⅤ资源中文在线天堂| 一区二区三区四区激情视频 | 亚洲avbb在线观看| 黄色一级大片看看| 国产精品av视频在线免费观看| 国产精品久久视频播放| 欧美高清性xxxxhd video| 国产毛片a区久久久久| 99热只有精品国产| av女优亚洲男人天堂| 国产在线男女| 美女免费视频网站| 欧美xxxx黑人xx丫x性爽| 中出人妻视频一区二区| 看免费成人av毛片| 国产高清激情床上av| 性色avwww在线观看| 国产成人a区在线观看| 中文亚洲av片在线观看爽| 五月玫瑰六月丁香| 最近在线观看免费完整版| 国产亚洲精品av在线| 热99在线观看视频| 毛片女人毛片| 波多野结衣巨乳人妻| 99国产极品粉嫩在线观看| 成人特级av手机在线观看| 国内精品久久久久久久电影| 色5月婷婷丁香| 在线观看66精品国产| 99热网站在线观看| 亚洲av免费高清在线观看| 国产精品野战在线观看| 少妇人妻一区二区三区视频| 长腿黑丝高跟| 可以在线观看的亚洲视频| 国产伦在线观看视频一区| 在线观看一区二区三区| 一边摸一边抽搐一进一小说| 又黄又爽又刺激的免费视频.| av在线蜜桃| 91久久精品国产一区二区三区| 女人被狂操c到高潮| 又黄又爽又免费观看的视频| 五月伊人婷婷丁香| 国产精品一区二区三区四区久久| 在线免费观看不下载黄p国产 | 国产男靠女视频免费网站| 91在线精品国自产拍蜜月| 欧洲精品卡2卡3卡4卡5卡区| 在线天堂最新版资源| 免费看美女性在线毛片视频| 成人鲁丝片一二三区免费| 日韩国内少妇激情av| 男女啪啪激烈高潮av片| 岛国在线免费视频观看| 久久久久久久久久久丰满 | 成年女人看的毛片在线观看| 欧美日韩中文字幕国产精品一区二区三区| x7x7x7水蜜桃| 看免费成人av毛片| 国产精品久久久久久精品电影| 简卡轻食公司| 在线观看舔阴道视频| 少妇丰满av| 男插女下体视频免费在线播放| 精品久久久噜噜| 人人妻,人人澡人人爽秒播| 天美传媒精品一区二区| 午夜久久久久精精品| 12—13女人毛片做爰片一| 极品教师在线免费播放| 亚洲欧美精品综合久久99| 欧美绝顶高潮抽搐喷水| 男女做爰动态图高潮gif福利片| 日韩亚洲欧美综合| 91在线观看av| 国产免费av片在线观看野外av| 毛片女人毛片| 真人做人爱边吃奶动态| 亚洲国产精品sss在线观看| 嫩草影院入口| 国产精品一区www在线观看 | 欧美日本视频| eeuss影院久久| 国产久久久一区二区三区| 欧美不卡视频在线免费观看| 99热只有精品国产| 国产日本99.免费观看| 网址你懂的国产日韩在线| 亚洲内射少妇av| 老司机福利观看| 麻豆久久精品国产亚洲av| 黄色视频,在线免费观看| 亚洲不卡免费看| 不卡一级毛片| 一级毛片久久久久久久久女| av黄色大香蕉| 男女啪啪激烈高潮av片| 欧美成人性av电影在线观看| 联通29元200g的流量卡| 村上凉子中文字幕在线| 久久久久国内视频| 午夜日韩欧美国产| 午夜a级毛片| 男女之事视频高清在线观看| 精品人妻一区二区三区麻豆 | 听说在线观看完整版免费高清| 色av中文字幕| 国产精品一区二区性色av| 国产伦人伦偷精品视频| 日本一本二区三区精品| 91av网一区二区| 亚洲四区av| 国产又黄又爽又无遮挡在线| 久久精品人妻少妇| 久久亚洲真实| 国产中年淑女户外野战色| 看十八女毛片水多多多| 成年女人永久免费观看视频| 免费黄网站久久成人精品| 两个人的视频大全免费| 亚洲av熟女| 啦啦啦观看免费观看视频高清| av国产免费在线观看| 亚洲国产色片| 亚洲色图av天堂| 制服丝袜大香蕉在线| 欧美潮喷喷水| 搡女人真爽免费视频火全软件 | 欧美最新免费一区二区三区| 国产精品日韩av在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 在线观看av片永久免费下载| 99久久无色码亚洲精品果冻| 99久久精品一区二区三区| 久久精品国产清高在天天线| 免费在线观看成人毛片| 国产白丝娇喘喷水9色精品| 天堂√8在线中文| 啦啦啦啦在线视频资源| 国产爱豆传媒在线观看| 亚洲av熟女| 中文资源天堂在线| 欧美日韩中文字幕国产精品一区二区三区| 日韩欧美国产在线观看| 亚洲精品日韩av片在线观看| 亚洲最大成人手机在线| 国产精品乱码一区二三区的特点| 88av欧美| 亚洲久久久久久中文字幕| 大型黄色视频在线免费观看| 日韩高清综合在线| 又黄又爽又免费观看的视频| 国产国拍精品亚洲av在线观看| 中文字幕免费在线视频6| 小说图片视频综合网站| 黄色视频,在线免费观看| 天堂网av新在线| 美女黄网站色视频| 日韩大尺度精品在线看网址| 亚洲熟妇熟女久久|