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

    兩種高分辨率氣溫實(shí)況融合產(chǎn)品在四川省的質(zhì)量評估

    2022-05-18 07:03:12蔣雨荷黃曉龍
    高原山地氣象研究 2022年1期
    關(guān)鍵詞:融合產(chǎn)品

    蔣雨荷 , 黃曉龍 , 杜 冰

    (1. 高原與盆地暴雨旱澇災(zāi)害四川省重點(diǎn)實(shí)驗(yàn)室, 成都 610072;2. 四川省氣象探測數(shù)據(jù)中心, 成都 610072)

    引言

    多年以來,地面氣象要素及土壤濕度等數(shù)據(jù)的獲取主要依靠地面人工站和自動站的觀測儀器進(jìn)行定時觀測,但限于站點(diǎn)離散且分布不均勻,難以覆蓋整個中國區(qū)域。利用數(shù)學(xué)方法將站點(diǎn)觀測值插值成格點(diǎn)數(shù)據(jù),在站點(diǎn)密集且空間變異少的區(qū)域,效果良好。但在地形復(fù)雜、區(qū)域廣闊、站點(diǎn)稀少的區(qū)域,尤其是在中國西部的川西高原等地區(qū),單純的數(shù)學(xué)方法難以獲得高質(zhì)量的格點(diǎn)實(shí)況數(shù)據(jù)[1]。近年來,國內(nèi)外學(xué)者將陸面數(shù)值模式和衛(wèi)星、雷達(dá)、站點(diǎn)觀測等信息進(jìn)行融合,以獲取更為可靠的地面氣象要素、陸表土壤溫濕度和地表熱通量等數(shù)據(jù)[2-3]。如美國大氣研究中心(NCAR)的高分辨率陸面數(shù)據(jù)同化系統(tǒng)(HRLDAS)[4-5]已經(jīng)具備提供土壤濕度的預(yù)報能力;韓國氣象局(KMA)發(fā)展的韓國陸面數(shù)據(jù)同化系統(tǒng)(KLDAS)[6]已用于提供中尺度天氣預(yù)報模式(WRF)模式的陸面初始場。

    我國在引進(jìn)國際先進(jìn)融合技術(shù)的基礎(chǔ)上,消化吸收并自主創(chuàng)新,建成了業(yè)務(wù)化的亞洲區(qū)域CMA陸面數(shù)據(jù)同化業(yè)務(wù)系統(tǒng)(CLDAS)和中國區(qū)域融合降水分析系統(tǒng)(CMPS-Hourly V2.1),以及實(shí)時運(yùn)行的全球海表溫度產(chǎn)品系統(tǒng)(CODAS-SST)和中國區(qū)域三維云融合系統(tǒng)(3DCloudAS)。其中CLDAS-V2.0實(shí)現(xiàn)了亞洲區(qū)域時間分辨率為1h、空間分辨率為0.0625°×0.0625°(5 km×5 km)的大氣驅(qū)動場(2 m氣溫、地面氣壓、2 m比濕、風(fēng)速、小時降水、短波輻射)和陸面要素集合分析(土壤濕度、地表溫度、土壤相對濕度和土壤體積含水量)等陸面產(chǎn)品的實(shí)時生產(chǎn)和發(fā)布[7-8]。2017年3月19日,為進(jìn)一步滿足氣象業(yè)務(wù)的需求,國家氣象信息中心研發(fā)的HRCLDAS-V1.0(高分辨率CMA陸面數(shù)據(jù)同化系統(tǒng)第一版)進(jìn)入業(yè)務(wù)化試運(yùn)行,該系統(tǒng)采用多種來源資料,包括地面站點(diǎn)觀測數(shù)據(jù)、ECMWF數(shù)值分析和預(yù)報產(chǎn)品以及DEM數(shù)字高程模型數(shù)據(jù),采用多重網(wǎng)格變分同化(STMAS)技術(shù)研制了2 m氣溫、2 m比濕、10 m風(fēng)、地面氣壓、降水等融合分析產(chǎn)品,時間分辨率為1 h,空間分辨率為0.01°×0.01°(1 km×1 km),該產(chǎn)品與國內(nèi)外同類產(chǎn)品比較,時空分辨率更高[9]。

    近幾年,國內(nèi)氣象學(xué)者在上述5 km和1 km模式融合產(chǎn)品的評估釋用方面做了一系列工作[9-12]。例如,龍柯吉等[9]對中國區(qū)域1 km氣溫實(shí)況融合產(chǎn)品進(jìn)行評估分析,結(jié)果表明:逐小時實(shí)況融合產(chǎn)品與站點(diǎn)實(shí)況基本一致,有較高的參考性,但在海拔較高,地形復(fù)雜的區(qū)域需要細(xì)致評估后使用。陳玉瑋等[10]通過檢驗(yàn)評估廣東省智能網(wǎng)格實(shí)況分析產(chǎn)品得出結(jié)論:氣溫網(wǎng)格實(shí)況的趨勢與站點(diǎn)實(shí)況相近,最低氣溫誤差大于最高氣溫。俞劍蔚等[11]分析評估了5 km國家級格點(diǎn)實(shí)況分析產(chǎn)品在江蘇省的適用性,表明日最高氣溫的誤差值較小,在江蘇的丘陵地帶誤差較大。王利亞[12]利用處于3種不同地形的站點(diǎn)數(shù)據(jù)對內(nèi)蒙古地區(qū)5 km氣溫實(shí)況融合產(chǎn)品進(jìn)行評估分析,發(fā)現(xiàn)氣溫融合產(chǎn)品的質(zhì)量從高到低依次為平原、丘陵、山地。

    四川省地處我國西南部,西部為高原山地,東部多為盆地、丘陵,下墊面復(fù)雜多樣,且高原臺站分布相對稀疏,對實(shí)況網(wǎng)格融合產(chǎn)品的適用性評估與服務(wù)應(yīng)用有著迫切的需求。2 m氣溫實(shí)況融合產(chǎn)品是CLDAS與HRCLDAS的重要組成部分,其可信度是評估實(shí)況融合產(chǎn)品質(zhì)量的重要參考依據(jù)。本文基于四川省156個國家級考核站、1282個區(qū)域考核站、1411個區(qū)域非考核站的氣溫觀測數(shù)據(jù), 通過相關(guān)系數(shù)、平均值誤差、均方根誤差等指標(biāo)對CLDAS(5 km×5 km)、HRCLDAS(1 km×1 km)2 m氣溫實(shí)況融合格點(diǎn)分析產(chǎn)品進(jìn)行質(zhì)量評估,以了解高、低分辨率實(shí)況融合產(chǎn)品對四川省氣溫的刻畫能力,以期為融合產(chǎn)品在區(qū)域氣候研究、防災(zāi)減災(zāi)、決策服務(wù)等方面的科學(xué)應(yīng)用提供技術(shù)支撐。

    1 資料與方法

    1.1 資料來源

    (1)本文所用的2 m氣溫實(shí)況融合格點(diǎn)分析產(chǎn)品(以下簡稱“融合產(chǎn)品”)數(shù)據(jù)為國家氣象信息中心提供的四川省2019年8月1日00時~2020年7月31日23時(北京時)1 km×1 km和5 km×5 km高分辨率逐小時氣溫數(shù)據(jù)。

    (2)參與檢驗(yàn)的站點(diǎn)數(shù)據(jù)為四川省氣象探測數(shù)據(jù)中心業(yè)務(wù)質(zhì)控流程質(zhì)控過的四川省地面逐小時氣溫觀測資料,并選取質(zhì)量控制碼為0(數(shù)據(jù)正確)的觀測數(shù)據(jù)進(jìn)行研究,包括國家站156個、區(qū)域考核站1282個、區(qū)域非考核站1411個。

    1.2 方法

    將四川省156個國家氣象觀測站、1282個區(qū)域考核站以及1411個區(qū)域非考核站逐小時氣溫觀測資料作為實(shí)況真值,把融合產(chǎn)品按雙線性插值的方法插值到站點(diǎn),通過統(tǒng)計比較一段時間內(nèi)兩者的評估指標(biāo)對數(shù)據(jù)進(jìn)行評估。評估檢驗(yàn)期間所有參與檢驗(yàn)的樣本均以小時觀測數(shù)據(jù)結(jié)果累積計算得到。

    雙線性插值方法即先在緯(經(jīng))向后在經(jīng)(緯)向進(jìn)行一元一次線性插值。

    評估指標(biāo)主要有:相關(guān)系數(shù)(COR)、平均值誤差(ME)、均方根誤差(RMSE)。具體公式如下:

    式中:Oi為站點(diǎn)觀測值,Gi為實(shí)況網(wǎng)格分析產(chǎn)品插值到檢驗(yàn)站點(diǎn)的數(shù)值,N為參與檢驗(yàn)的樣本總樣本數(shù)(站次數(shù))。

    2 兩種融合產(chǎn)品的質(zhì)量對比評估

    從四川省2019年8月1日00時~2020年7月31日23時的整體評估情況(表1)來看,1 km融合產(chǎn)品與實(shí)況資料的相關(guān)程度更高,其質(zhì)量較5 km有明顯的改進(jìn),兩種產(chǎn)品的平均值誤差分別為-0.069℃和-0.083℃,均方根誤差分別為1.419℃和1.87℃。綜合全年評估情況來看,兩種融合產(chǎn)品與實(shí)況站點(diǎn)資料的相關(guān)程度高,誤差小,具有很高的參考性。

    表1 兩種氣溫融合產(chǎn)品評估結(jié)果

    2.1 誤差的空間評估

    參考地區(qū)地貌類型和市州管轄范圍,將四川分為攀西地區(qū)(攀枝花和涼山州)、川西高原(甘孜州和阿壩州)、盆地中部(資陽、遂寧和內(nèi)江)、盆地南部(自貢、瀘州和宜賓)、盆地西南部(雅安、眉山和樂山)、盆地東北部(南充、巴中、達(dá)州和廣安)、盆地西北部(成都、德陽、綿陽和廣元)7個地區(qū)(圖1)。

    圖1 四川不同地區(qū)劃分

    2.1.1 分站點(diǎn)評估

    2.1.1.1 國家站

    如圖2a、b所示,國家站的相關(guān)系數(shù)基本都在0.98以上,表明兩種融合產(chǎn)品與實(shí)況資料的相關(guān)程度較高。其中,相關(guān)系數(shù)較小的臺站主要分布在甘孜州的九龍和德格、阿壩州的松潘和茂縣、涼山州的美姑(以下簡稱“三州”)以及盆地內(nèi)的峨眉山,大都位于攀西平原與青藏高原的過渡帶和高山河谷地區(qū),海拔高且地勢起伏大。

    如圖2c、d所示,平均值誤差主要集中在-1~1℃,位于三州西部邊緣、三州-盆地過渡帶的臺站平均值誤差大都在-1~0℃,而位于盆地中部、三州中部地區(qū)的臺站平均值誤差大都在0~1℃。造成這一現(xiàn)象的原因是: 位于三州-盆地過渡帶、三州西部邊緣的臺站,如寶興、雅江、九寨溝等,大都位于高山河谷區(qū)域,地勢較低,氣溫的觀測值也會略高于四周,使融合產(chǎn)品的插值結(jié)果小于觀測值,造成平均值誤差為負(fù),出現(xiàn)低估值;反之,位于盆地中部的臺站為避免受周圍建筑物的遮擋,選址相對較高,如漢源、寧南、峨眉山等,較高的地勢,使得氣溫略低于周邊區(qū)域,融合產(chǎn)品的插值結(jié)果大于觀測值,造成平均值誤差為正,出現(xiàn)高估值。融合產(chǎn)品對地形因子的影響比較敏感,而1 km融合產(chǎn)品由于分辨率更高,在這些地形復(fù)雜區(qū)域的評估效果要優(yōu)于5 km融合產(chǎn)品。

    如圖2e、f所示,均方根誤差介于0~1.8℃的臺站數(shù)較多,主要分布在盆地內(nèi)部、三州中部;三州-盆地過渡帶、甘孜州西部、阿壩州北部地區(qū)的平均絕對誤差值較大,在此區(qū)域內(nèi)5 km融合產(chǎn)品的誤差值多≥0.6℃,而1 km融合產(chǎn)品有了明顯的優(yōu)化,誤差值主要集中在0~0.6℃。

    圖2 1 km(左)和5 km(右)融合產(chǎn)品與國家站的對比評估結(jié)果(a、b. COR,c、d. ME,e、f. RMSE)

    2.1.1.2 區(qū)域考核站

    如圖3a、b所示,區(qū)域考核站的相關(guān)性在盆地內(nèi)部最好,相關(guān)系數(shù)基本≥0.98;位于盆地與高原過渡帶、甘孜州西南部、阿壩州北部臺站的相關(guān)性較差,相關(guān)系數(shù)<0.92的臺站數(shù)較多。如圖3c、d所示,盆地中部、西南部、南部的平均值誤差介于-1~1℃的臺站最多,高估值略多于低估值;盆地西北部、東北部以及盆地高原過渡帶上的平均值誤差<-1℃的臺站數(shù)較多,此區(qū)域地勢起伏大,地形較復(fù)雜,位于高山河谷的臺站數(shù)較多,多出現(xiàn)低估值。如圖3e、f所示,盆地東北部、中部的均方根誤差主要集中在0~1.8℃,三州中部、攀枝花地區(qū)主要集中在0.6~2.7℃,三州與盆地過渡帶、甘孜州西部、阿壩州南部主要集中在1.8~3.6℃,且個別臺站的均方根誤差≥3.6℃。1 km融合產(chǎn)品在三州、三州-盆地過渡帶等地形較復(fù)雜區(qū)域的誤差值較5 km融合產(chǎn)品有明顯的減小。

    圖3 同圖2,但為區(qū)域考核站

    2.1.1.3 區(qū)域非考核站

    區(qū)域非考核站在盆地西南部、東北部分布較考核站稠密,填補(bǔ)了國家站與區(qū)域考核站在此區(qū)域的評估空白。區(qū)域非考核站在盆地中部、南部的相關(guān)性與考核站較為相近(圖4a、b),主要集中在0.98~1,三州地區(qū)的相關(guān)性主要集中在0.92~0.98,位于三州盆地過渡帶臺站的相關(guān)性主要集中在0.92~0.96,且個別臺站的相關(guān)系數(shù)<0.92。非考核站的平均值誤差在盆地西南部、東北部多為低估值,誤差大都≤-1℃(圖4c、d),均方根誤差在此區(qū)域主要集中在0.6~2.7℃(圖4e、f)。

    2.1.2 分市州評估

    通過對各市州的評估可以看出:兩種融合產(chǎn)品與實(shí)況值的相關(guān)系數(shù)都比較高,在0.92以上,且盆地內(nèi)的相關(guān)性高于三州地區(qū),1 km融合產(chǎn)品明顯優(yōu)于5 km(圖5a)。盆地內(nèi)大部分市為正平均值誤差,其原因可能是:這些臺站為了減小周圍建筑物對測站環(huán)境的影響多選址在地勢相對較高的山坡上,具有海拔較低、地勢平坦開闊、地形復(fù)雜度低等特點(diǎn),從而出現(xiàn)小的高估值。而盆地內(nèi)的巴中、雅安、眉山、廣元以及三州、攀枝花等地則為負(fù)平均值誤差,這些市(州)的臺站多位于高原山區(qū)的峽谷地區(qū)、大巴山與四川盆地的交界處以及低山丘陵等地形較為復(fù)雜的區(qū)域,具有地勢起伏大、海拔較高、地形復(fù)雜度較大等特點(diǎn)(圖5b)。均方根誤差(圖5c)則呈現(xiàn)出盆地 <攀枝花<三州的趨勢,盆地內(nèi)的均方根誤差介于0.5~2.5℃,攀枝花地區(qū)介于1.5~2.8℃,三州地區(qū)介于2.0~3.5℃。相對而言,1 km融合產(chǎn)品在三州與攀枝花地區(qū)的誤差值明顯低于5 km產(chǎn)品。

    2.1.3 分海拔評估

    為研究海拔高度與評估結(jié)果的關(guān)系,按照分級海拔高度對溫度進(jìn)行評估,將海拔高度分為五級:≤1000 m(低海拔)、1000~2000 m(中海拔)、2000~4000 m(亞高海拔)、4000~6000 m(高海拔)和≥6000 m(極高海拔),由于參與評估站點(diǎn)的最高海拔為4304.1 m,因此評估海拔高度只有前四級(表2)。經(jīng)統(tǒng)計,1000 m以下的站點(diǎn)有2247個,1000~2000 m的站點(diǎn)有424個,2000~4000 m的站點(diǎn)有296個,4000 m以上的站點(diǎn)有9個(圖6)。可以看出,兩種融合產(chǎn)品隨著海拔升高,參與評估的站點(diǎn)個數(shù)逐漸減少, 4000 m以下站點(diǎn)的相關(guān)系數(shù)隨海拔增高而逐漸減小,平均值誤差與均方根誤差逐漸增大;而4000 m以上的9個站點(diǎn)雖然海拔高,但位置較為平坦開闊,受高山、低谷等復(fù)雜地形的影響小,地形復(fù)雜度也較小,使其誤差值低于2000~4000 m的站點(diǎn);相對而言,2000~4000 m的站點(diǎn),囊括了丘陵、山地、中山峽谷等地形,海拔高且地形復(fù)雜度較大,評估效果也最差??梢姡瑢τ诟吆0吻业匦螐?fù)雜度較高的地區(qū),融合產(chǎn)品的誤差較大,效果比低海拔、地勢平坦地區(qū)差。

    圖6 四川省各站點(diǎn)海拔分布

    表2 氣溫融合產(chǎn)品在四川省分海拔評估結(jié)果

    2.2 誤差的時間評估

    2.2.1 按日評估

    如圖7所示,兩種融合產(chǎn)品各時次的相關(guān)系數(shù)都比較高, 5 km融合產(chǎn)品的相關(guān)系數(shù)都在0.965以上,1 km融合產(chǎn)品的相關(guān)系數(shù)都在0.980以上;09時達(dá)到一天中相關(guān)系數(shù)的波峰,14時達(dá)到波谷;1 km與5 km的平均值誤差都在-0.4~0.2℃,并且在09~19時誤差較大,在14時達(dá)到波谷;均方根誤差的分布與平均值誤差呈現(xiàn)負(fù)相關(guān),于14時出現(xiàn)波峰,均方根誤差值達(dá)到最大,1 km的均方根誤差值較5 km有明顯改善??傮w而言,1 km與5 km融合產(chǎn)品具有一定的日變化趨勢,但變化幅度不大,與實(shí)況值較為一致,午后融合產(chǎn)品的誤差較其余時刻略大,代表性相對較差。

    圖7 兩種融合產(chǎn)品評估結(jié)果的日變化(a. COR,b. ME,c. RMSE)

    2.2.2 按季節(jié)評估

    如圖8所示,兩種融合產(chǎn)品的相關(guān)性在春、夏、秋季都比較高,集中在0.95~0.99,冬季1 km融合產(chǎn)品的相關(guān)性較5 km融合產(chǎn)品有明顯的提升,且與秋季的相關(guān)性較為接近,分別為0.959與0.957;四季的平均值誤差較穩(wěn)定,介于-0.17~0.05℃,春、秋、冬季為負(fù)平均值誤差,夏季為正平均值誤差;均方根誤差呈現(xiàn)“U”型分布,夏季誤差值最小,春秋季較為接近,冬季的誤差值最大,介于1.2~2.2℃。整體來看,兩種融合產(chǎn)品在四季的變化趨勢較為一致,誤差值均在夏季最小,冬季最差,1 km融合產(chǎn)品較5 km在冬季的提升效果更顯著。

    圖8 兩種融合產(chǎn)品評估結(jié)果的季節(jié)變化(a. COR,b. ME,c. RMSE)

    3 結(jié)論與討論

    本文利用四川省國家級考核站、區(qū)域考核站及區(qū)域非考核站的氣溫觀測數(shù)據(jù), 通過相關(guān)系數(shù)、平均值誤差、均方根誤差等指標(biāo),對1 km和5 km的2 m氣溫融合產(chǎn)品進(jìn)行質(zhì)量評估,得到如下結(jié)論:

    (1)兩種融合產(chǎn)品的相關(guān)性都很好,1 km融合產(chǎn)品的平均值誤差與均方根誤差的評估效果更優(yōu)。地形與海拔高度對融合產(chǎn)品的影響較大,在地勢起伏較大區(qū)域,二者誤差明顯增大。地處低于周圍地形的臺站易出現(xiàn)氣溫負(fù)平均值誤差(低估值),地處高于周圍地形的臺站易出現(xiàn)正平均值誤差(高估值)。

    (2)整體評估效果在盆地內(nèi)最優(yōu),盆地內(nèi)大部臺站呈現(xiàn)正平均值誤差的高估值;盆地-三州過渡帶以及地勢起伏較大的盆地北部、西南部地區(qū)次之,呈現(xiàn)較小的負(fù)平均值誤差的低估值;三州、攀枝花地區(qū)最差,呈現(xiàn)較大的負(fù)平均值誤差的低估值。

    (3)兩種融合產(chǎn)品都具有一定的日變化趨勢,但變化幅度不大,與實(shí)況值較為一致,午后融合產(chǎn)品的誤差較其余時刻略大,代表性相對較差。

    (4)兩種融合產(chǎn)品在春、夏、秋季的評估效果優(yōu)于冬季,其中夏季為弱的高估值,春、秋、冬為低估值。1 km融合產(chǎn)品在評估效果最差的冬季較5 km融合產(chǎn)品有顯著提升。

    總體來看,1 km與5 km融合產(chǎn)品均能較好地反映四川省的氣溫變化,在評估效果較差的地形條件復(fù)雜的高海拔地區(qū),1 km融合產(chǎn)品較5 km有明顯的改善,但對該區(qū)域的評估訂正,還有待進(jìn)一步研究分析。

    猜你喜歡
    融合產(chǎn)品
    好產(chǎn)品,可持續(xù)
    從靈感出發(fā),邂逅好產(chǎn)品
    一次函數(shù)“四融合”
    村企黨建聯(lián)建融合共贏
    融合菜
    從創(chuàng)新出發(fā),與高考數(shù)列相遇、融合
    寬窄融合便攜箱IPFS500
    《融合》
    2015產(chǎn)品LOOKBOOK直擊
    Coco薇(2015年1期)2015-08-13 02:23:50
    新產(chǎn)品
    玩具(2009年10期)2009-11-04 02:33:14
    色94色欧美一区二区| 美女国产高潮福利片在线看| 99热全是精品| 精品人妻一区二区三区麻豆| 色94色欧美一区二区| 成人黄色视频免费在线看| √禁漫天堂资源中文www| 日韩免费高清中文字幕av| 天天躁夜夜躁狠狠躁躁| 久久久水蜜桃国产精品网| 俄罗斯特黄特色一大片| 日本vs欧美在线观看视频| av天堂久久9| 国产在视频线精品| 在线观看人妻少妇| 国产一区二区激情短视频 | 欧美日韩av久久| 少妇粗大呻吟视频| 精品乱码久久久久久99久播| 中亚洲国语对白在线视频| 国产区一区二久久| 国产成人av激情在线播放| www日本在线高清视频| 97在线人人人人妻| 99精国产麻豆久久婷婷| 国产欧美日韩一区二区精品| 黄色a级毛片大全视频| 亚洲欧美日韩另类电影网站| 亚洲成人国产一区在线观看| a级毛片在线看网站| 爱豆传媒免费全集在线观看| 搡老熟女国产l中国老女人| 亚洲欧美日韩另类电影网站| 1024香蕉在线观看| 欧美精品一区二区免费开放| 亚洲avbb在线观看| 菩萨蛮人人尽说江南好唐韦庄| 爱豆传媒免费全集在线观看| 热99国产精品久久久久久7| 韩国精品一区二区三区| 久久亚洲国产成人精品v| 母亲3免费完整高清在线观看| 男人添女人高潮全过程视频| 天堂中文最新版在线下载| 麻豆国产av国片精品| 欧美另类一区| 一级黄色大片毛片| 日本精品一区二区三区蜜桃| tocl精华| 久久ye,这里只有精品| 欧美性长视频在线观看| 精品国产一区二区三区四区第35| 色94色欧美一区二区| 久久精品久久久久久噜噜老黄| 美女扒开内裤让男人捅视频| 一级黄色大片毛片| 国产欧美日韩一区二区三 | 精品第一国产精品| a 毛片基地| 久久久久久免费高清国产稀缺| 亚洲精品第二区| 日本a在线网址| a在线观看视频网站| 午夜激情久久久久久久| 水蜜桃什么品种好| 他把我摸到了高潮在线观看 | 女人精品久久久久毛片| 欧美 亚洲 国产 日韩一| 无遮挡黄片免费观看| 成人免费观看视频高清| 亚洲国产av影院在线观看| 黄色怎么调成土黄色| 丰满饥渴人妻一区二区三| 国产伦理片在线播放av一区| 成人手机av| 91大片在线观看| 日韩一区二区三区影片| 亚洲久久久国产精品| 成人手机av| 看免费av毛片| av在线app专区| 俄罗斯特黄特色一大片| 中文字幕色久视频| 老司机在亚洲福利影院| 久久毛片免费看一区二区三区| 99精国产麻豆久久婷婷| 中文字幕色久视频| 精品福利观看| 99精品欧美一区二区三区四区| 国产有黄有色有爽视频| 国产激情久久老熟女| 精品福利观看| 19禁男女啪啪无遮挡网站| 午夜老司机福利片| 久久 成人 亚洲| 韩国精品一区二区三区| 亚洲七黄色美女视频| 两人在一起打扑克的视频| 日韩免费高清中文字幕av| tube8黄色片| 高清黄色对白视频在线免费看| 美女主播在线视频| 精品久久久精品久久久| 亚洲成人免费av在线播放| 91av网站免费观看| 午夜成年电影在线免费观看| 婷婷色av中文字幕| 国产av精品麻豆| 一本—道久久a久久精品蜜桃钙片| 精品少妇内射三级| 久久亚洲精品不卡| 如日韩欧美国产精品一区二区三区| 人人妻人人澡人人看| 王馨瑶露胸无遮挡在线观看| 青春草视频在线免费观看| 中文字幕人妻丝袜制服| 欧美大码av| 国产无遮挡羞羞视频在线观看| 精品人妻在线不人妻| av超薄肉色丝袜交足视频| av一本久久久久| 超碰97精品在线观看| 在线天堂中文资源库| 色婷婷av一区二区三区视频| av一本久久久久| 国产精品熟女久久久久浪| 国产日韩欧美在线精品| 亚洲精华国产精华精| 久久影院123| 国产老妇伦熟女老妇高清| 国产主播在线观看一区二区| 中亚洲国语对白在线视频| 日本猛色少妇xxxxx猛交久久| 国产免费视频播放在线视频| 免费少妇av软件| 男女下面插进去视频免费观看| 亚洲综合色网址| 少妇 在线观看| 十八禁网站免费在线| 日本一区二区免费在线视频| 99久久99久久久精品蜜桃| av在线老鸭窝| 中文欧美无线码| 免费在线观看日本一区| 精品亚洲乱码少妇综合久久| 一级,二级,三级黄色视频| av天堂久久9| 狠狠精品人妻久久久久久综合| 亚洲久久久国产精品| 最黄视频免费看| 看免费av毛片| 亚洲第一av免费看| 亚洲熟女精品中文字幕| 亚洲av成人一区二区三| 热99久久久久精品小说推荐| bbb黄色大片| 又黄又粗又硬又大视频| 久久av网站| 老司机午夜十八禁免费视频| 国产精品久久久久久精品电影小说| 亚洲成人免费av在线播放| 电影成人av| 91老司机精品| 国产日韩欧美在线精品| 成人手机av| 在线观看舔阴道视频| 国产真人三级小视频在线观看| 日本a在线网址| 91大片在线观看| 嫩草影视91久久| 91av网站免费观看| 精品国内亚洲2022精品成人 | 亚洲精品久久午夜乱码| 国产在线观看jvid| 中文欧美无线码| 欧美日韩国产mv在线观看视频| 国产高清视频在线播放一区 | 久久久国产成人免费| 久久精品国产综合久久久| 午夜影院在线不卡| 男女边摸边吃奶| 久久 成人 亚洲| 欧美97在线视频| 亚洲熟女精品中文字幕| 色老头精品视频在线观看| 波多野结衣av一区二区av| 亚洲精品自拍成人| 亚洲av成人一区二区三| 日本精品一区二区三区蜜桃| 中文精品一卡2卡3卡4更新| 亚洲成人免费电影在线观看| 国产一区二区激情短视频 | 美女中出高潮动态图| 久久性视频一级片| 色婷婷久久久亚洲欧美| 大香蕉久久网| 久久ye,这里只有精品| av线在线观看网站| 在线观看免费日韩欧美大片| 国产精品熟女久久久久浪| 91精品三级在线观看| 菩萨蛮人人尽说江南好唐韦庄| 男女国产视频网站| 国产精品麻豆人妻色哟哟久久| 国产伦人伦偷精品视频| 国产av国产精品国产| tube8黄色片| 一级,二级,三级黄色视频| 国产成人a∨麻豆精品| 日韩电影二区| 亚洲一卡2卡3卡4卡5卡精品中文| 在线观看舔阴道视频| 91麻豆精品激情在线观看国产 | 精品人妻1区二区| av超薄肉色丝袜交足视频| 一区二区三区四区激情视频| 久久久精品免费免费高清| 国产国语露脸激情在线看| 色老头精品视频在线观看| 黄色视频不卡| 国产高清国产精品国产三级| 欧美av亚洲av综合av国产av| 69精品国产乱码久久久| 777久久人妻少妇嫩草av网站| 人妻一区二区av| 国产精品秋霞免费鲁丝片| 精品少妇一区二区三区视频日本电影| 美女高潮到喷水免费观看| 多毛熟女@视频| 自线自在国产av| 这个男人来自地球电影免费观看| 亚洲中文日韩欧美视频| 欧美黄色片欧美黄色片| 亚洲精品国产一区二区精华液| 999久久久精品免费观看国产| 成人亚洲精品一区在线观看| 人人妻人人添人人爽欧美一区卜| 精品亚洲乱码少妇综合久久| 久久香蕉激情| 国产精品1区2区在线观看. | 成年av动漫网址| 免费看十八禁软件| 老司机亚洲免费影院| 久久人妻熟女aⅴ| 亚洲欧美日韩高清在线视频 | 成在线人永久免费视频| 一级毛片精品| 亚洲av美国av| 免费久久久久久久精品成人欧美视频| 12—13女人毛片做爰片一| 在线看a的网站| 黄色a级毛片大全视频| 久久狼人影院| 国产成人免费观看mmmm| 自线自在国产av| 丝袜人妻中文字幕| 日韩 亚洲 欧美在线| 国产精品一区二区在线观看99| www.精华液| 韩国精品一区二区三区| av福利片在线| 精品福利观看| 中文字幕另类日韩欧美亚洲嫩草| 搡老乐熟女国产| 极品少妇高潮喷水抽搐| 在线观看免费午夜福利视频| 欧美亚洲日本最大视频资源| av超薄肉色丝袜交足视频| 女人爽到高潮嗷嗷叫在线视频| av视频免费观看在线观看| 日韩中文字幕视频在线看片| 每晚都被弄得嗷嗷叫到高潮| 黄片大片在线免费观看| 精品少妇一区二区三区视频日本电影| 黄频高清免费视频| 亚洲国产欧美网| 免费在线观看完整版高清| 大片电影免费在线观看免费| 男人舔女人的私密视频| 欧美日韩视频精品一区| 啦啦啦啦在线视频资源| 国产黄频视频在线观看| 97精品久久久久久久久久精品| 90打野战视频偷拍视频| 看免费av毛片| 一本一本久久a久久精品综合妖精| 久久久久视频综合| 嫩草影视91久久| 日本av免费视频播放| 亚洲欧美精品综合一区二区三区| 免费女性裸体啪啪无遮挡网站| 欧美久久黑人一区二区| 午夜免费观看性视频| 极品人妻少妇av视频| 不卡一级毛片| 嫁个100分男人电影在线观看| 搡老熟女国产l中国老女人| 成年人免费黄色播放视频| 搡老乐熟女国产| 美女中出高潮动态图| 十分钟在线观看高清视频www| 激情视频va一区二区三区| 母亲3免费完整高清在线观看| 亚洲成人手机| 国产精品一区二区在线观看99| 91大片在线观看| av欧美777| 久久久久视频综合| 欧美精品高潮呻吟av久久| 日本wwww免费看| 亚洲国产欧美日韩在线播放| a在线观看视频网站| 欧美日韩精品网址| 99国产精品一区二区蜜桃av | 午夜免费鲁丝| 热99久久久久精品小说推荐| 亚洲精品国产一区二区精华液| 成年av动漫网址| 后天国语完整版免费观看| 一级a爱视频在线免费观看| 免费观看av网站的网址| 99久久99久久久精品蜜桃| 国产在线免费精品| 亚洲 欧美一区二区三区| 蜜桃在线观看..| 国产亚洲午夜精品一区二区久久| 十八禁高潮呻吟视频| 一进一出抽搐动态| 三上悠亚av全集在线观看| 国产精品偷伦视频观看了| 欧美日韩av久久| 男女免费视频国产| 免费日韩欧美在线观看| 天堂中文最新版在线下载| 一本综合久久免费| 欧美日韩亚洲综合一区二区三区_| 久久热在线av| 黑人巨大精品欧美一区二区蜜桃| 最近最新免费中文字幕在线| a 毛片基地| 亚洲精品粉嫩美女一区| 欧美另类一区| 欧美变态另类bdsm刘玥| 一级毛片电影观看| 欧美老熟妇乱子伦牲交| 91老司机精品| 在线天堂中文资源库| 一级毛片女人18水好多| 十八禁网站网址无遮挡| 亚洲精品久久午夜乱码| 十八禁网站网址无遮挡| 免费一级毛片在线播放高清视频 | tocl精华| 亚洲国产成人一精品久久久| 99精品久久久久人妻精品| 久久国产精品影院| 视频在线观看一区二区三区| 建设人人有责人人尽责人人享有的| av天堂久久9| 亚洲国产毛片av蜜桃av| 精品欧美一区二区三区在线| 国产xxxxx性猛交| 欧美一级毛片孕妇| 精品国产一区二区三区四区第35| 亚洲第一欧美日韩一区二区三区 | 美女午夜性视频免费| 亚洲免费av在线视频| 另类精品久久| 亚洲中文av在线| 精品熟女少妇八av免费久了| www.999成人在线观看| 久久久久久久大尺度免费视频| tube8黄色片| 91国产中文字幕| 日日夜夜操网爽| 久久国产精品男人的天堂亚洲| av网站免费在线观看视频| 久久久久久久久免费视频了| 婷婷成人精品国产| 人人妻人人澡人人爽人人夜夜| 免费不卡黄色视频| 色婷婷久久久亚洲欧美| 久久久精品免费免费高清| 欧美日韩视频精品一区| 国产成+人综合+亚洲专区| 国产xxxxx性猛交| 国产色视频综合| 亚洲成人手机| 亚洲伊人色综图| 久久精品成人免费网站| 一本大道久久a久久精品| 视频区欧美日本亚洲| 午夜日韩欧美国产| 51午夜福利影视在线观看| 岛国在线观看网站| 国产黄色免费在线视频| 在线观看www视频免费| 久久性视频一级片| 亚洲精品国产av蜜桃| 免费在线观看黄色视频的| 伊人久久大香线蕉亚洲五| 热99re8久久精品国产| 精品人妻一区二区三区麻豆| 成人三级做爰电影| 国产淫语在线视频| 久久久久精品国产欧美久久久 | 伊人亚洲综合成人网| 97人妻天天添夜夜摸| 日韩大码丰满熟妇| 国产高清视频在线播放一区 | 涩涩av久久男人的天堂| 亚洲国产欧美在线一区| 国产欧美日韩一区二区三 | 男女无遮挡免费网站观看| 法律面前人人平等表现在哪些方面 | 国产男女超爽视频在线观看| 国产在视频线精品| 亚洲欧美一区二区三区黑人| 侵犯人妻中文字幕一二三四区| 亚洲欧美精品自产自拍| 水蜜桃什么品种好| 中文欧美无线码| 纯流量卡能插随身wifi吗| 极品人妻少妇av视频| 欧美精品高潮呻吟av久久| 午夜久久久在线观看| 两个人免费观看高清视频| 亚洲精品国产色婷婷电影| 91精品三级在线观看| 91精品伊人久久大香线蕉| 久久人妻福利社区极品人妻图片| 一本大道久久a久久精品| 天天操日日干夜夜撸| 一二三四在线观看免费中文在| 亚洲一码二码三码区别大吗| 久热爱精品视频在线9| 超色免费av| 中文字幕制服av| 国产极品粉嫩免费观看在线| 91成人精品电影| 国产成人免费观看mmmm| 久久久精品免费免费高清| 久久精品国产亚洲av高清一级| 9191精品国产免费久久| 精品福利永久在线观看| 丰满饥渴人妻一区二区三| 男人操女人黄网站| 成人国产一区最新在线观看| 大码成人一级视频| 国产91精品成人一区二区三区 | 麻豆乱淫一区二区| 久久香蕉激情| 国产一区有黄有色的免费视频| e午夜精品久久久久久久| 国产又爽黄色视频| 亚洲精品国产一区二区精华液| 欧美黑人精品巨大| 91老司机精品| 国产精品av久久久久免费| 欧美精品啪啪一区二区三区 | 中文字幕高清在线视频| 免费女性裸体啪啪无遮挡网站| 在线十欧美十亚洲十日本专区| 精品福利永久在线观看| 老汉色∧v一级毛片| 国产成人av教育| 免费黄频网站在线观看国产| 亚洲欧洲日产国产| 国产欧美亚洲国产| 国产1区2区3区精品| 久久人妻熟女aⅴ| 成年人黄色毛片网站| 中文字幕高清在线视频| 欧美激情极品国产一区二区三区| 国产男人的电影天堂91| av在线老鸭窝| 国产精品秋霞免费鲁丝片| 美女视频免费永久观看网站| 亚洲欧美成人综合另类久久久| 女人高潮潮喷娇喘18禁视频| 欧美另类一区| 亚洲精品国产精品久久久不卡| 男男h啪啪无遮挡| 中文字幕另类日韩欧美亚洲嫩草| a级毛片黄视频| 国产福利在线免费观看视频| 一本综合久久免费| 欧美+亚洲+日韩+国产| 国产在视频线精品| videosex国产| 亚洲精品av麻豆狂野| 成人亚洲精品一区在线观看| 亚洲精品国产av成人精品| 欧美精品一区二区大全| 狂野欧美激情性bbbbbb| 一级,二级,三级黄色视频| 国产精品国产三级国产专区5o| 伊人亚洲综合成人网| 国产精品一区二区精品视频观看| 天堂中文最新版在线下载| 欧美国产精品va在线观看不卡| 一个人免费在线观看的高清视频 | 老司机深夜福利视频在线观看 | 日韩免费高清中文字幕av| 中文字幕人妻丝袜一区二区| 脱女人内裤的视频| 91精品三级在线观看| 久久热在线av| 亚洲成av片中文字幕在线观看| 久久久久精品国产欧美久久久 | 久久精品成人免费网站| 国产精品国产三级国产专区5o| 亚洲av欧美aⅴ国产| 99热全是精品| 一区在线观看完整版| 每晚都被弄得嗷嗷叫到高潮| 日韩大码丰满熟妇| 999久久久精品免费观看国产| 欧美成狂野欧美在线观看| 性色av一级| 汤姆久久久久久久影院中文字幕| 久久久久国产精品人妻一区二区| 男女边摸边吃奶| 天堂俺去俺来也www色官网| 国产av又大| 亚洲全国av大片| 久久精品久久久久久噜噜老黄| 丰满饥渴人妻一区二区三| 欧美成狂野欧美在线观看| 老熟妇仑乱视频hdxx| 汤姆久久久久久久影院中文字幕| 国产视频一区二区在线看| 一区在线观看完整版| 久久久久久人人人人人| 麻豆国产av国片精品| 国产欧美日韩一区二区三区在线| tocl精华| 国产真人三级小视频在线观看| 久久国产精品人妻蜜桃| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品第二区| 国产麻豆69| 亚洲视频免费观看视频| 人妻人人澡人人爽人人| 亚洲美女黄色视频免费看| 国产亚洲精品第一综合不卡| 久久午夜综合久久蜜桃| 91av网站免费观看| 婷婷成人精品国产| 欧美日韩亚洲国产一区二区在线观看 | 女人精品久久久久毛片| 亚洲视频免费观看视频| 黄色a级毛片大全视频| 日韩,欧美,国产一区二区三区| 午夜福利一区二区在线看| 国产97色在线日韩免费| 啦啦啦免费观看视频1| 亚洲精品国产色婷婷电影| 久久天堂一区二区三区四区| 亚洲欧洲精品一区二区精品久久久| 精品乱码久久久久久99久播| 男男h啪啪无遮挡| 亚洲一码二码三码区别大吗| 一边摸一边做爽爽视频免费| 黄色片一级片一级黄色片| 免费少妇av软件| 中文字幕最新亚洲高清| 国产免费av片在线观看野外av| 一级片'在线观看视频| 久久亚洲精品不卡| 久久影院123| 国产在线视频一区二区| 色播在线永久视频| 国内毛片毛片毛片毛片毛片| 我要看黄色一级片免费的| 黄频高清免费视频| 欧美精品啪啪一区二区三区 | 制服诱惑二区| 嫩草影视91久久| 18禁裸乳无遮挡动漫免费视频| 日本91视频免费播放| 成年女人毛片免费观看观看9 | 黄片大片在线免费观看| 欧美 日韩 精品 国产| av片东京热男人的天堂| 99国产精品免费福利视频| 老熟妇乱子伦视频在线观看 | 少妇 在线观看| 夜夜夜夜夜久久久久| 高清视频免费观看一区二区| 亚洲午夜精品一区,二区,三区| 操出白浆在线播放| 久久天堂一区二区三区四区| 免费在线观看视频国产中文字幕亚洲 | 国产成+人综合+亚洲专区| 欧美人与性动交α欧美精品济南到| 精品亚洲乱码少妇综合久久| 人人澡人人妻人| 久久99一区二区三区| 国产在线一区二区三区精| 亚洲精品国产色婷婷电影| 天堂中文最新版在线下载| 亚洲 国产 在线| 9色porny在线观看| 啪啪无遮挡十八禁网站| 国产野战对白在线观看| 一区二区三区四区激情视频| 这个男人来自地球电影免费观看| 涩涩av久久男人的天堂| 脱女人内裤的视频| 久久精品成人免费网站| 国产又爽黄色视频| 亚洲av成人一区二区三| 男人操女人黄网站| 欧美变态另类bdsm刘玥|