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

    CLDAS氣溫和降水產(chǎn)品在內(nèi)蒙古地區(qū)適用性分析

    2023-11-14 01:57:52董祝雷趙艷麗馮曉晶劉詩夢
    干旱氣象 2023年5期
    關鍵詞:內(nèi)蒙古地區(qū)平均氣溫降水量

    董祝雷,趙艷麗,馮曉晶,劉詩夢

    (內(nèi)蒙古自治區(qū)氣候中心,內(nèi)蒙古 呼和浩特 010000)

    引言

    近年來公眾對氣象服務精細化程度和準確率的要求不斷提高,各類氣象業(yè)務及服務對格點化實況資料需求越來越強烈。隨著氣象觀測手段不斷增加,氣象數(shù)據(jù)種類大大增多,利用數(shù)據(jù)融合和同化技術,綜合多源觀測資料及多模式模擬數(shù)據(jù),能夠獲取高精度、高質(zhì)量、時空連續(xù)的多源融合格點產(chǎn)品,其在天氣預報、氣候研究、防災減災等工作中發(fā)揮重要作用(師春香等,2019;張君霞等,2022;郭陽等,2023;楊柳等,2023)。

    我國多源氣象數(shù)據(jù)融合研究起步相對較晚,現(xiàn)階段業(yè)務服務主要依賴于國外一些再分析產(chǎn)品數(shù)據(jù)集,如美國國家環(huán)境預報中心(National Centers for Environmental Prediction,NCEP)、國家大氣研究中心(National Center for Atmospheric Research,NCAR)與能源部(Department Of Energy,DOE)聯(lián)合發(fā)布的NCAR(NCEP1)、NCEP/DOE(NCEP2)數(shù)據(jù)(Kalnary et al.,1996)以及英國東英吉利大學(University of East Anglia)氣候研究中心的CRU(聞新宇等,2006)、歐洲中期天氣預報中心的ERA-Interim(Dee et al.,2011)、日本氣象廳JRA55(Kobayashi and Iwasaki,2016)等再分析數(shù)據(jù),這些數(shù)據(jù)一般具有時間序列較長、適用性較好等特征。大氣再分析資料側(cè)重于同化大氣參數(shù),用以提高大氣參數(shù)模擬質(zhì)量,但隨著精細化服務需求的不斷提升,這些數(shù)據(jù)無論在時空分辨率上還是精確度上已無法很好地滿足區(qū)域分析的需求。

    中國氣象局2014年啟動了國家氣象科技創(chuàng)新工程“氣象資料質(zhì)量控制及多源數(shù)據(jù)融合與再分析”,建成了業(yè)務化的亞洲區(qū)域中國氣象局陸面數(shù)據(jù)同化系統(tǒng)(CMA Land Data Assimilation System,CLDAS)(師春香等,2019)。相較于再分析資料,CLDAS產(chǎn)品側(cè)重于同化地面觀測要素,進而獲得較為精確的陸面過程相關變量。由于CLDAS數(shù)據(jù)時空分辨率較高,推出后得到了全國各地科研和業(yè)務工作者的本地化釋用。研究表明,CLDAS的數(shù)據(jù)質(zhì)量優(yōu)于數(shù)值模式資料,能夠較好地反映觀測要素的時空分布特征,在我國東部地區(qū)的適用性高于西部地區(qū)(龔偉偉,2014;孫帥等,2017),且土壤溫度、濕度和氣溫產(chǎn)品優(yōu)于全球陸面數(shù)據(jù)同化系統(tǒng)產(chǎn)品(韓帥等,2017;崔園園等,2018;崔園園等,2019;Chen and Yuan,2020;陳艷麗等,2020;劉佩佩等,2021),但是也不可避免地存在一定的誤差,尤其高海拔地區(qū)誤差較大,因此在地形復雜、氣象站稀疏區(qū)域應用該產(chǎn)品需謹慎(Han et al.,2020)。如CLDAS氣溫產(chǎn)品在貴州、山西等測站分布稀疏和地形起伏較大地區(qū)誤差較大(董春卿等,2021;楊富燕等,2023);CLDAS的10 m風產(chǎn)品在江蘇地區(qū)誤差較大,與業(yè)務服務需求有一定差距(俞劍蔚等,2019)。

    國務院氣象高質(zhì)量發(fā)展綱要(2022—2035年)強調(diào)“要提高氣象服務經(jīng)濟高質(zhì)量發(fā)展水平”,這對精細化氣象服務提出了更高要求,而高精度數(shù)據(jù)集是提升精細化服務的基礎。內(nèi)蒙古幅員遼闊、地形復雜,既有山丘林地,又有高原草地和戈壁沙漠,氣象站點較為稀疏,無法滿足精細化服務要求。因此,挖掘適用于內(nèi)蒙古地區(qū)的高分辨率、高精度數(shù)據(jù)產(chǎn)品對提高內(nèi)蒙古氣象精細化服務水平有重要意義。

    陸面數(shù)據(jù)同化產(chǎn)品和大氣再分析產(chǎn)品彌補了觀測資料時空分布不均的缺陷,為大氣科學領域研究提供了重要的數(shù)據(jù)支撐。然而,不同數(shù)據(jù)集在運算、同化的觀測資料及方法等方面存在差異,使得不同數(shù)據(jù)集的可信度也有差異(葉夢姝,2018),因此有必要比較多種再分析數(shù)據(jù)之間的差異(Simmons et al.,2004;趙天保和符淙斌,2009)。鑒于此,本文對2008—2019年CLDAS的2 m氣溫和降水產(chǎn)品在內(nèi)蒙古地區(qū)的適用性進行評估,同時與CRU TS和ERA5兩套再分析資料進行對比,評估CLDAS產(chǎn)品在內(nèi)蒙古地區(qū)的適用性,以期為開展精細化氣象服務和氣候資源利用等提供科技支撐。

    1 資料與方法

    1.1 資 料

    內(nèi)蒙古位于我國北部邊疆,地勢由東北向西南斜伸,呈狹長形(圖1),東西長約2 400 km,南北最大跨度超過1 700 km,總面積118.3萬km2,橫跨東北、華北、西北地區(qū)。內(nèi)蒙古地貌主要以高原為主,大部分地區(qū)海拔在1 000 m以上,東部的大興安嶺,經(jīng)燕山、陰山與西部的賀蘭山蜿蜒相連,呈拉長的“S”型中脊橫貫全區(qū),成為分隔南部平原和北部高原的屏障。

    圖1 內(nèi)蒙古地形高度和國家級氣象觀測站分布Fig.1 The distribution of terrain height and national meteorological stations in Inner Mongolia

    所用資料包括:(1)2008—2019年內(nèi)蒙古119個國家級氣象觀測站逐日氣溫和降水數(shù)據(jù),站點分布如圖1。(2)國家氣象信息中心下發(fā)的2008—2019年月尺度的CN05.1格點化氣溫和降水觀測數(shù)據(jù)集。該數(shù)據(jù)集采用距平逼近法,由氣候場和距平場分別插值后疊加得到(吳佳和高學杰,2013),使用了中國境內(nèi)2 400個臺站觀測資料,空間分辨率為0.25°×0.25°。(3)2008—2019年逐小時2 m氣溫和降水格點數(shù)據(jù),空間分辨率為0.062 5°×0.062 5°,來源于CLDAS產(chǎn)品數(shù)據(jù)集。該產(chǎn)品數(shù)據(jù)集采用數(shù)據(jù)融合與同化技術,對地面觀測數(shù)據(jù)、衛(wèi)星遙感觀測資料和數(shù)值模式產(chǎn)品等多源數(shù)據(jù)進行融合同化,可以提供逐小時降水、氣溫等6個要素共計15種網(wǎng)格實況分析產(chǎn)品(孫小龍等,2015;單帥等,2021)。(4)2008—2019年歐洲中期天氣預報中心ERA5逐月氣溫和降水數(shù)據(jù)集(https://www.ecmwf.int/en/forecasts/datasets/browse-reanalysis-datasets),空間分辨率0.1°×0.1°。ERA5再分析資料是歐洲中期天氣預報中心2016年采用新一代四維變分同化技術制作的第五代大氣再分析產(chǎn)品,采用實時更新方式提供1950年至今的高分辨率全球大氣再分析數(shù)據(jù),其時空分辨率和精度相比上一代產(chǎn)品有明顯提升(呂潤清和李響,2021;張俊兵等,2021)。(5)CRU TS4.05數(shù)據(jù)集中2008—2019年月尺度氣溫和降水數(shù)據(jù)(https://crudata.uea.ac.uk/cru/data/hrg/cru_ts_4.05/),空間分辨率為0.5°×0.5°。CRU TS數(shù)據(jù)集由英國東英吉利大學氣候研究小組制作,是世界公認的質(zhì)量較高的數(shù)據(jù),具有很高的可靠性(聞新宇等,2006;季飛,2015)。

    文中地圖涉及的內(nèi)蒙古自治區(qū)行政邊界基于內(nèi)蒙古自然資源廳標準地圖網(wǎng)站下載的審圖號為蒙S(2019)33號的標準地圖制作,底圖無修改。

    1.2 方法

    為便于CLDAS、CRU TS和ERA5三種數(shù)據(jù)集與CN05.1格點觀測數(shù)據(jù)的空間誤差比較,采用局部區(qū)域平均分別將CLDAS和ERA5數(shù)據(jù)插值為與CN05.1數(shù)據(jù)相同的網(wǎng)格(0.25°×0.25°),將CN05.1觀測數(shù)據(jù)插值為與CRU TS數(shù)據(jù)相同的網(wǎng)格(0.5°×0.5°)。此外,基于內(nèi)蒙古119站坐標信息,分別提取對應坐標位置CLDAS、CRU TS和ERA5數(shù)據(jù)產(chǎn)品的氣溫和降水序列,并與站點觀測資料對比,計算平均絕對誤差(Mean Absolute Error,MAE)、均方根誤差(Root Mean Square Error,RMSE)、相關系數(shù),計算公式參考《氣象統(tǒng)計分析與預報方法》(黃嘉佑,2016)。

    2 降水和氣溫產(chǎn)品適用性評估

    2.1 空間分布對比

    圖2 是2008—2019年內(nèi)蒙古地區(qū)CN05.1觀測數(shù)據(jù)集與CLDAS、CRU TS和ERA5產(chǎn)品的平均年降水量空間分布及其差值。可以看到,CLDAS、CRU TS和ERA5降水產(chǎn)品數(shù)據(jù)均表現(xiàn)出內(nèi)蒙古年降水量從東向西逐漸遞減的分布特征[圖2(b)、(c)、(d)]。其中,CLDAS數(shù)據(jù)呼倫貝爾市和興安盟東部為降水大值區(qū),鄂爾多斯市東部到呼和浩特市南部為次大值區(qū),阿拉善盟為降水低值區(qū)。與CN05.1觀測數(shù)據(jù)相比,CLDAS在內(nèi)蒙古中部和東部地區(qū)降水普遍偏少,在內(nèi)蒙古西部地區(qū)總體偏多,而在湖泊等低洼處降水異常偏多[圖2(e)],如呼倫湖和達里諾爾湖;CRU TS降水數(shù)據(jù)在內(nèi)蒙古北部地區(qū)描述較好,但在中東部及西南地區(qū)降水普遍偏少、東南部和西部地區(qū)降水普遍偏多[圖2(f)];ERA5降水數(shù)據(jù)低估了內(nèi)蒙古地區(qū)年降水量,大部地區(qū)年降水量偏少超過40 mm[圖2(g)]。

    圖2 2008—2019年內(nèi)蒙古地區(qū)平均年降水量觀測值(a)、三種數(shù)據(jù)集的平均年降水量(b、c、d)及其與觀測值的差值(e、f、g)空間分布(單位:mm)(a)CN05.1,(b、e)CLDAS,(c、f)CRU TS,(d、g)ERA5Fig.2 The spatial distribution of average annual precipitation from the observation (a) and three datasets (b, c, d), and difference between them and observation values (e, f, g) in Inner Mongolia during 2008-2019 (Unit: mm)(a) CN05.1, (b, e) CLDAS, (c, f) CRU TS, (d, g) ERA5

    圖3是2008—2019年內(nèi)蒙古地區(qū)CN05.1與CLDAS、CRU TS、ERA5的年平均氣溫空間分布及其差值。CN05.1氣溫數(shù)據(jù)[圖3(a)]顯示,內(nèi)蒙古地區(qū)年平均氣溫從西南向東北遞減,呼倫貝爾市北部為低值區(qū),內(nèi)蒙古西部為高值區(qū),東南部的赤峰市和通遼市一帶年平均氣溫較高。CLDAS、CRU TS和ERA5氣溫產(chǎn)品能夠較好地刻畫這一空間分布特征[圖3(b)、(c)、(d)],但在部分地區(qū)仍與CN05.1觀測數(shù)據(jù)有較大差異。CLDAS年平均氣溫在內(nèi)蒙古東北部大興安嶺及河套平原北部陰山較低,能夠表現(xiàn)出地形起伏對局地氣溫變化的影響。與CN05.1觀測數(shù)據(jù)相比,CLDAS年平均氣溫在內(nèi)蒙古大部地區(qū)偏高,尤其在阿拉善盟東部、巴彥淖爾市西部及赤峰市西部偏差較大[圖3(e)];CRU TS年平均氣溫除了呼倫貝爾市東部和阿拉善盟大部地區(qū)偏低外,其余大部地區(qū)均偏高,不能完全反映地形變化對氣溫場的影響,尤其在陰山山脈氣溫偏高明顯[圖3(f)];ERA5年平均氣溫在內(nèi)蒙古地區(qū)普遍明顯偏高,尤其在呼倫貝爾市北部偏差高達2.0 ℃左右,但呼倫貝爾市南部、興安盟大部、赤峰市北部、烏蘭察布市南部以及阿拉善盟東部等地區(qū)氣溫變化與觀測數(shù)據(jù)接近,且在陰山山脈能夠表現(xiàn)出地形對局地氣溫的影響[圖3(g)]。

    圖3 2008—2019年內(nèi)蒙古地區(qū)年平均氣溫觀測值(a)、三種數(shù)據(jù)集的年平均氣溫(b、c、d)及其與觀測值的差值(e、f、g)空間分布(單位:℃)(a)CN05.1,(b、e)CLDAS,(c、f)CRU TS,(d、g)ERA5Fig.3 The spatial distribution of annual mean temperature from the observation (a) and three datasets (b, c, d), and difference between them and observation values (e, f, g) in Inner Mongolia during 2008-2019 (Unit: ℃)(a) CN05.1, (b, e) CLDAS, (c, f) CRU TS, (d, g) ERA5

    2.2 線性趨勢對比

    進一步對比2008—2019年內(nèi)蒙古地區(qū)三種數(shù)據(jù)集降水產(chǎn)品與CN05.1格點觀測數(shù)據(jù)的線性趨勢(圖4)。CN05.1降水觀測數(shù)據(jù)[圖4(a)]顯示,內(nèi)蒙古大部地區(qū)年降水量呈增加趨勢,呼倫貝爾市東南部、興安盟至通遼市一帶以及鄂爾多斯市中東部和阿拉善盟南部增加速率較高,而呼倫貝爾市北部和巴彥淖爾市北部年降水量呈明顯減少趨勢,錫林郭勒盟北部降水量有弱的減少趨勢。CLDAS降水數(shù)據(jù)基本能夠反映2008—2019年內(nèi)蒙古年降水量線性趨勢[圖4(b)],但在呼倫貝爾市北部、錫林郭勒盟西北部和巴彥淖爾市年降水量減少地區(qū)以及阿拉善盟東部年降水量增加區(qū)域降水變率偏大,而在赤峰市、通遼市和興安盟等部分地區(qū)年降水量呈減少趨勢,與觀測數(shù)據(jù)相反。CRU TS降水數(shù)據(jù)年變化率[圖4(c)]在內(nèi)蒙古東部地區(qū)與CN05.1觀測數(shù)據(jù)相似,但呼倫貝爾市北部年降水量減少速率和呼倫貝爾市東部年降水量增加速率均小于觀測,且在內(nèi)蒙古西部地區(qū)年降水變率偏小,尤其在鄂爾多斯市東部和阿拉善盟南部明顯低于觀測。ERA5降水數(shù)據(jù)對阿拉善盟東南部及鄂爾多斯市東部到烏蘭察布市一帶年降水量增加趨勢表現(xiàn)較好,而在其他地區(qū)年降水變率偏?。蹐D4(d)]。

    圖4 2008—2019年內(nèi)蒙古地區(qū)CN05.1(a)、CLDAS(b)、CRU TS(c)和ERA5(d)年降水量線性趨勢(單位:mm·a-1)Fig.4 The linear trend of annual precipitation from CN05.1 (a), CLDAS (b), CRU TS (c) and ERA5 (d)datasets in Inner Mongolia during 2008-2019 (Unit: mm·a-1)

    CN05.1格點觀測資料[圖5(a)]顯示,2008—2019年內(nèi)蒙古大部地區(qū)年平均氣溫呈升高趨勢,增溫率為0.03~0.06 ℃·a-1,呼倫貝爾市東北部和鄂爾多斯市北部與巴彥淖爾市交界處增溫率較低。CLDAS數(shù)據(jù)[圖5(b)]顯示,內(nèi)蒙古大部地區(qū)呈增暖趨勢,且大部地區(qū)增溫率高于觀測,尤其阿拉善盟東南部增溫率與觀測偏差較大,而阿拉善盟中西部、巴彥淖爾市西部、赤峰市西部年平均氣溫呈明顯下降趨勢。CRU TS年平均氣溫線性趨勢[圖5(c)]表現(xiàn)為全區(qū)一致增暖的分布格局,大部地區(qū)增溫率超過0.07 ℃·a-1,較觀測明顯偏高。ERA5年平均氣溫線性趨勢分布[圖5(d)]與觀測較一致,除阿拉善盟西部和呼倫貝爾市東北部地區(qū)增溫率與觀測接近外,其他大部地區(qū)增溫率高于觀測。

    圖5 2008—2019年內(nèi)蒙古地區(qū)CN05.1(a)、CLDAS(b)、CRU TS(c)和ERA5(d)年平均氣溫線性趨勢(單位:℃·a-1)Fig.5 The linear trend of annual mean temperature from CN05.1 (a), CLDAS (b), CRU TS (c) and ERA5 (d) datasets in Inner Mongolia during 2008-2019 (Unit :℃·a-1)

    3 氣溫和降水產(chǎn)品誤差評估

    采用距離最近格點與站點進行匹配,提取內(nèi)蒙古119個國家站上述三種數(shù)據(jù)集逐月降水量和平均氣溫數(shù)據(jù),統(tǒng)計其與站點觀測數(shù)據(jù)的相關系數(shù)、均方根誤差及平均絕對誤差(表1)??梢钥闯?,三種數(shù)據(jù)集的平均氣溫與站點觀測值相關性高,相關系數(shù)均超過0.998,CLDAS平均氣溫的相關系數(shù)最高達0.999 5,其MAE僅0.38 ℃,遠小于ERA5的0.85 ℃和CRU TS的1.13 ℃。CRU TS和ERA5的降水量與站點觀測值的相關系數(shù)低于0.9,而CLDAS降水量的相關系數(shù)達0.965 7。從誤差來看,CLDAS降水量與站點觀測值MAE為4.73 mm,而CRU TS和ERA5降水量的MAE則在10.00 mm左右。可見,無論是平均氣溫還是降水量,CLDAS產(chǎn)品與站點觀測資料的相關性最高、誤差最小。

    表1 三種數(shù)據(jù)集的月平均氣溫和降水量與內(nèi)蒙古119站觀測值相關系數(shù)和誤差統(tǒng)計Tab.1 The statistics of correlation coefficient and error of monthly mean temperature and precipitation between three datasets and observation values at 119 stations in Inner Mongolia

    進一步對三種數(shù)據(jù)集的季節(jié)平均氣溫和降水量誤差進行統(tǒng)計(表2)。與觀測數(shù)據(jù)相比,CLDAS夏季平均氣溫偏低,春、秋、冬季平均氣溫偏高,且冬、秋季MAE稍高,分別為0.45、0.33 ℃,春、夏季MAE略小,在0.30 ℃左右;CRU TS平均氣溫除冬季偏高外,其他季節(jié)均偏低,冬、春、夏季MAE均超過1.01 ℃,秋季誤差最小為0.87 ℃;ERA5春、夏季平均氣溫偏低,秋、冬季平均氣溫偏高,冬季MAE最大為1.26 ℃,其他季節(jié)誤差在0.60 ℃左右??梢?,內(nèi)蒙古地區(qū)CLDAS各季平均氣溫誤差均小于ERA5和CRU TS,ERA5平均氣溫數(shù)據(jù)質(zhì)量整體優(yōu)于CRU TS,僅冬季平均氣溫MAE略高于CRU TS。從四季降水誤差看,CLDAS對內(nèi)蒙古春季和秋季降水量低估、夏季和冬季降水量高估,而CRU TS降水產(chǎn)品夏季高估、其他季節(jié)均低估;ERA5普遍低估內(nèi)蒙古夏季降水,高估春季和冬季降水,對秋季降水也存在低估現(xiàn)象。另外,三種降水產(chǎn)品在內(nèi)蒙古地區(qū)MAE均是夏季最大,秋季次之,冬季最小。

    表2 三種數(shù)據(jù)集的季節(jié)平均氣溫和降水量與內(nèi)蒙古119站觀測值偏差和平均絕對誤差統(tǒng)計Tab.2 The statistics of bias and mean absolute error of seasonal mean temperature and precipitation between three datasets and observation values at 119 stations in Inner Mongolia

    圖6是內(nèi)蒙古119站季節(jié)平均氣溫和降水量觀測值與CLDAS、CRU TS和ERA5產(chǎn)品平均絕對誤差箱形圖。從誤差范圍、中位數(shù)、上下四分位數(shù)來看,CLDAS資料最優(yōu)。綜上可見, CLDAS氣溫和降水產(chǎn)品在內(nèi)蒙古地區(qū)準確性明顯優(yōu)于CRU TS和ERA5相關產(chǎn)品。

    圖6 內(nèi)蒙古119站季節(jié)平均氣溫(a)和降水量(b)與CLDAS、CRU TS和ERA5產(chǎn)品的平均絕對誤差分布Fig.6 The distribution of mean absolute error between seasonal mean temperature (a), precipitation (b) at 119 stations in Inner Mongolia and CLDAS, CRU TS, ERA5 products

    4 CLDAS產(chǎn)品誤差空間分布

    CLDAS氣溫和降水產(chǎn)品在內(nèi)蒙古地區(qū)的適用性整體優(yōu)于CRU TS和ERA5,那么CLDAS氣溫和降水數(shù)據(jù)誤差在內(nèi)蒙古分布如何?從CLDAS氣溫和降水產(chǎn)品與內(nèi)蒙古119站觀測值的相關系數(shù)和均方根誤差空間分布(圖7)看出,整體上CLDAS氣溫數(shù)據(jù)精度優(yōu)于降水數(shù)據(jù),CLDAS氣溫與內(nèi)蒙古站點觀測值的相關系數(shù)普遍超過0.999,僅內(nèi)蒙古中西部和河套地區(qū)局部站點相關系數(shù)為0.997~0.999,阿拉善左旗站相關系數(shù)最小為0.988 6;從氣溫均方根誤差分布看,RMSE超過1.0 ℃的站點主要位于阿拉善盟西部及赤峰市西北部與錫林郭勒盟交界地區(qū)。CLDAS降水與內(nèi)蒙古站點觀測值的相關系數(shù)大部分在0.95以上,僅4站相關系數(shù)低于0.90,其中烏海站相關系數(shù)最低為0.77,可見CLDAS降水產(chǎn)品在內(nèi)蒙古西部尤其是河套地區(qū)精度較差;從降水的均方根誤差分布看,內(nèi)蒙古北部地區(qū)誤差明顯小于南部地區(qū)。

    圖7 CLDAS氣溫(a、b)和降水(c、d)產(chǎn)品與內(nèi)蒙古119站觀測值的相關系數(shù)(a、c)和均方根誤差(b、d)空間分布Fig.7 The spatial distribution of correlation coefficient (a, c) and root mean square error (b, d) between CLDAS temperature(a, b), precipitation (c, d) products and observation values at 119 stations in Inner Mongolia

    5 結論與討論

    基于內(nèi)蒙古119個國家氣象站觀測資料和國家氣象信息中心CN05.1格點觀測數(shù)據(jù)集,對CLDAS產(chǎn)品數(shù)據(jù)集的2 m平均氣溫和降水量在內(nèi)蒙古地區(qū)的適用性進行檢驗評估,同時與CRU TS和ERA5再分析數(shù)據(jù)集進行對比,得到如下結論:

    (1)CLDAS、CRU TS和ERA5數(shù)據(jù)集均能較好表現(xiàn)內(nèi)蒙古年降水量和年平均氣溫空間分布特征。CLDAS產(chǎn)品低估了內(nèi)蒙古中東部地區(qū)年降水量、高估了西部地區(qū)年降水量,尤其對呼倫湖、達里諾爾湖等地勢低洼地區(qū)降水嚴重高估;CRU TS對內(nèi)蒙古東部和西南地區(qū)年降水量存在低估現(xiàn)象,而對東南部和西部地區(qū)年降水量存在高估現(xiàn)象。相較于CLDAS和CRU TS,ERA5嚴重低估了內(nèi)蒙古地區(qū)年降水量。三種數(shù)據(jù)集對內(nèi)蒙古大部地區(qū)年平均氣溫存在高估現(xiàn)象,其中CLDAS對阿拉善盟東部地區(qū)平均氣溫嚴重高估,CRU TS對內(nèi)蒙古中西部以及大興安嶺和陰山山脈平均氣溫嚴重高估,而ERA5嚴重高估了呼倫貝爾市北部地區(qū)平均氣溫。

    (2)CLDAS和CRU TS降水數(shù)據(jù)對內(nèi)蒙古降水變化趨勢的表現(xiàn)優(yōu)于ERA5。內(nèi)蒙古地區(qū),CRU TS和ERA5資料與CN05.1格點觀測資料增暖趨勢一致,但增溫率明顯高于觀測,而CLDAS資料還反映出內(nèi)蒙古局地變冷趨勢。

    (3)與站點觀測的誤差評估表明,CLDAS氣溫和降水產(chǎn)品質(zhì)量明顯優(yōu)于CRU TS和ERA5數(shù)據(jù)集。CLDAS平均氣溫與觀測值的相關系數(shù)為0.999 5,高于CRU TS和ERA5,平均絕對誤差(0.38 ℃)遠小于ERA5(0.85 ℃)和CRU TS(1.13 ℃);CRU TS和ERA5降水量與觀測值的相關系數(shù)低于0.88,而CLDAS降水產(chǎn)品的相關系數(shù)超過0.96,平均絕對誤差為4.7 mm,低于CRU TS和ERA5(超過9.5 mm)。季節(jié)尺度上,CLDAS氣溫產(chǎn)品在內(nèi)蒙古地區(qū)夏季偏低,其他季節(jié)偏高,而降水產(chǎn)品春季和秋季低估、夏季和冬季高估。

    (4)CLDAS氣溫數(shù)據(jù)主要在阿拉善盟東部、河套地區(qū)誤差較大,而降水數(shù)據(jù)在內(nèi)蒙古北部地區(qū)誤差明顯小于南部地區(qū),且河套平原氣溫和降水與站點觀測值的相關系數(shù)較低。

    CLDAS通過融合地面氣象觀測提高大氣驅(qū)動數(shù)據(jù)的質(zhì)量,進而提高數(shù)據(jù)產(chǎn)品的模擬精度(單帥等,2021)。CLDAS降水和氣溫產(chǎn)品空間分辨率高,其數(shù)據(jù)質(zhì)量在內(nèi)蒙古大部地區(qū)優(yōu)于CRU TS和ERA5再分析資料,具有較高的應用價值(武榮盛等,2021),但該數(shù)據(jù)產(chǎn)品在局部地區(qū)存在一定誤差,如河套地區(qū)。此外,CLDAS降水產(chǎn)品局地受地形影響較大,須加以修訂。

    致謝:感謝國家氣象信息中心數(shù)據(jù)分析科提供的CLDAS中國區(qū)域逐小時陸面融合分析產(chǎn)品及技術指導!

    猜你喜歡
    內(nèi)蒙古地區(qū)平均氣溫降水量
    繪制和閱讀降水量柱狀圖
    第十三屆中國音樂金鐘獎內(nèi)蒙古地區(qū)選拔賽
    草原歌聲(2021年2期)2021-10-14 08:27:56
    降水量是怎么算出來的
    啟蒙(3-7歲)(2019年8期)2019-09-10 03:09:08
    烏蘭縣近38年氣溫變化特征分析
    從全球氣候變暖大背景看萊州市30a氣溫變化
    內(nèi)蒙古地區(qū)蒙古族常用野生植物相關傳統(tǒng)知識研究
    1988—2017年呼和浩特市降水演變特征分析
    1981—2010年拐子湖地區(qū)氣溫變化特征及趨勢分析
    近50年來全球背景下青藏高原氣候變化特征分析
    基于小波變換的三江平原旬降水量主周期識別
    一区二区三区精品91| 日韩中文字幕欧美一区二区| 国产高清视频在线播放一区| 亚洲成人精品中文字幕电影| 国产精品精品国产色婷婷| 久久久国产成人免费| 成人18禁在线播放| 国产熟女xx| 身体一侧抽搐| 亚洲精品美女久久久久99蜜臀| 中亚洲国语对白在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲午夜理论影院| 国产野战对白在线观看| 女性生殖器流出的白浆| 真人一进一出gif抽搐免费| 91大片在线观看| 国产精品99久久99久久久不卡| 丁香欧美五月| 男女下面插进去视频免费观看| e午夜精品久久久久久久| 窝窝影院91人妻| 激情视频va一区二区三区| 一级a爱片免费观看的视频| 黑丝袜美女国产一区| 日本一区二区免费在线视频| 亚洲中文日韩欧美视频| 乱人伦中国视频| 国产三级黄色录像| 国产av又大| 久久影院123| 在线观看66精品国产| 99re在线观看精品视频| a级毛片在线看网站| 熟女少妇亚洲综合色aaa.| 欧美绝顶高潮抽搐喷水| 一区福利在线观看| 亚洲人成电影观看| 亚洲 欧美一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 亚洲av美国av| 麻豆国产av国片精品| 久久人人爽av亚洲精品天堂| 亚洲成人久久性| 深夜精品福利| 一区在线观看完整版| 精品久久久久久久久久免费视频| 欧美亚洲日本最大视频资源| 亚洲中文av在线| 久久香蕉激情| 看黄色毛片网站| 亚洲专区中文字幕在线| 99国产综合亚洲精品| 久99久视频精品免费| 国产亚洲精品综合一区在线观看 | 中亚洲国语对白在线视频| av片东京热男人的天堂| 大型av网站在线播放| 日韩有码中文字幕| 亚洲欧美精品综合久久99| 色综合亚洲欧美另类图片| 午夜福利高清视频| 国产视频一区二区在线看| 一边摸一边抽搐一进一出视频| 亚洲情色 制服丝袜| 欧洲精品卡2卡3卡4卡5卡区| 国产97色在线日韩免费| 久久精品亚洲熟妇少妇任你| 亚洲成a人片在线一区二区| 99精品久久久久人妻精品| 51午夜福利影视在线观看| 美女高潮喷水抽搐中文字幕| 国产一区在线观看成人免费| 久久久国产欧美日韩av| 久久香蕉精品热| 亚洲一区中文字幕在线| 在线十欧美十亚洲十日本专区| 美国免费a级毛片| av天堂久久9| 中出人妻视频一区二区| 国产精华一区二区三区| 欧美成人性av电影在线观看| АⅤ资源中文在线天堂| 国产成人系列免费观看| 成人国产一区最新在线观看| 午夜精品国产一区二区电影| 精品不卡国产一区二区三区| 色老头精品视频在线观看| 国产亚洲精品av在线| 91大片在线观看| 亚洲在线自拍视频| 中文字幕高清在线视频| 校园春色视频在线观看| 日本在线视频免费播放| 一边摸一边抽搐一进一小说| 欧美激情高清一区二区三区| 国产精品日韩av在线免费观看 | 久久久精品国产亚洲av高清涩受| 老鸭窝网址在线观看| 叶爱在线成人免费视频播放| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品久久电影中文字幕| 日韩欧美在线二视频| 九色国产91popny在线| 中文字幕高清在线视频| 久久天堂一区二区三区四区| 亚洲国产精品sss在线观看| 美女免费视频网站| av欧美777| 久久精品国产亚洲av香蕉五月| 亚洲第一av免费看| 丁香欧美五月| 涩涩av久久男人的天堂| 国产高清激情床上av| 欧美午夜高清在线| 美女高潮喷水抽搐中文字幕| 叶爱在线成人免费视频播放| 极品教师在线免费播放| 99国产综合亚洲精品| 国产亚洲精品久久久久5区| 久久久精品欧美日韩精品| 性色av乱码一区二区三区2| 黑丝袜美女国产一区| 精品国产乱码久久久久久男人| 国产精品 国内视频| 久久精品影院6| 韩国av一区二区三区四区| 国产亚洲欧美在线一区二区| 丰满的人妻完整版| 免费高清视频大片| 日本免费a在线| 久久久久久国产a免费观看| 91av网站免费观看| 一边摸一边抽搐一进一小说| 国产麻豆69| 欧美丝袜亚洲另类 | 欧美黄色淫秽网站| 国内久久婷婷六月综合欲色啪| 久久久国产成人精品二区| 91麻豆av在线| 黑人欧美特级aaaaaa片| 啦啦啦观看免费观看视频高清 | 少妇被粗大的猛进出69影院| 久久精品aⅴ一区二区三区四区| 欧美黑人精品巨大| 91精品三级在线观看| 免费高清在线观看日韩| 久久精品成人免费网站| 青草久久国产| 狂野欧美激情性xxxx| 一卡2卡三卡四卡精品乱码亚洲| 国内精品久久久久精免费| 91九色精品人成在线观看| www日本在线高清视频| 日日摸夜夜添夜夜添小说| 中文字幕高清在线视频| 欧美激情久久久久久爽电影 | 久久天躁狠狠躁夜夜2o2o| 国产精品秋霞免费鲁丝片| 国产亚洲欧美精品永久| 正在播放国产对白刺激| 国产精品野战在线观看| 午夜福利影视在线免费观看| 亚洲av五月六月丁香网| 欧美一区二区精品小视频在线| 午夜a级毛片| 国产亚洲欧美在线一区二区| 高清黄色对白视频在线免费看| 视频在线观看一区二区三区| 午夜精品久久久久久毛片777| bbb黄色大片| √禁漫天堂资源中文www| 亚洲色图综合在线观看| 免费在线观看视频国产中文字幕亚洲| 国产精品1区2区在线观看.| 亚洲精品国产区一区二| 久久国产乱子伦精品免费另类| 午夜a级毛片| 长腿黑丝高跟| 免费人成视频x8x8入口观看| 国产成人影院久久av| 一区二区三区精品91| 国产伦人伦偷精品视频| 欧美成狂野欧美在线观看| 97超级碰碰碰精品色视频在线观看| 国产亚洲欧美精品永久| 麻豆一二三区av精品| 女性被躁到高潮视频| 动漫黄色视频在线观看| 此物有八面人人有两片| 亚洲,欧美精品.| 色综合站精品国产| 一级,二级,三级黄色视频| 99香蕉大伊视频| 亚洲国产中文字幕在线视频| 高潮久久久久久久久久久不卡| 老司机深夜福利视频在线观看| 欧美+亚洲+日韩+国产| 亚洲中文字幕一区二区三区有码在线看 | 淫秽高清视频在线观看| 国产精品九九99| 欧美激情高清一区二区三区| 嫁个100分男人电影在线观看| 欧美日韩一级在线毛片| 国产精品永久免费网站| 精品高清国产在线一区| 老司机深夜福利视频在线观看| 精品国产超薄肉色丝袜足j| 国产又爽黄色视频| 亚洲熟妇熟女久久| 亚洲精品中文字幕一二三四区| 国产av又大| www.精华液| 又黄又爽又免费观看的视频| 韩国精品一区二区三区| 精品乱码久久久久久99久播| 女生性感内裤真人,穿戴方法视频| 久久天躁狠狠躁夜夜2o2o| 久久久久国产一级毛片高清牌| 国产成人av激情在线播放| 曰老女人黄片| 亚洲av成人不卡在线观看播放网| 久久精品国产亚洲av香蕉五月| 一个人免费在线观看的高清视频| 国产精华一区二区三区| 欧美中文日本在线观看视频| 久久久久精品国产欧美久久久| 精品国产国语对白av| 少妇粗大呻吟视频| 亚洲激情在线av| 国产高清有码在线观看视频 | √禁漫天堂资源中文www| 真人做人爱边吃奶动态| 亚洲熟女毛片儿| 免费少妇av软件| 国产一区二区激情短视频| 这个男人来自地球电影免费观看| 久久久国产成人精品二区| 免费观看精品视频网站| 午夜两性在线视频| 国产又爽黄色视频| 国产一区二区激情短视频| 又大又爽又粗| av片东京热男人的天堂| 91老司机精品| 超碰成人久久| 黄色 视频免费看| 国产av一区二区精品久久| 亚洲九九香蕉| 日韩一卡2卡3卡4卡2021年| 精品高清国产在线一区| 手机成人av网站| 亚洲专区中文字幕在线| 亚洲国产精品sss在线观看| 长腿黑丝高跟| а√天堂www在线а√下载| 两人在一起打扑克的视频| 久久人人爽av亚洲精品天堂| 亚洲欧美日韩另类电影网站| 国产不卡一卡二| 国产一区二区三区在线臀色熟女| 中亚洲国语对白在线视频| 在线观看免费午夜福利视频| www.熟女人妻精品国产| 日韩精品中文字幕看吧| 91成人精品电影| 俄罗斯特黄特色一大片| 69精品国产乱码久久久| 中文字幕最新亚洲高清| 欧美不卡视频在线免费观看 | 欧美另类亚洲清纯唯美| 97人妻天天添夜夜摸| 亚洲精品久久成人aⅴ小说| 制服丝袜大香蕉在线| 精品无人区乱码1区二区| 亚洲色图综合在线观看| 悠悠久久av| 美女午夜性视频免费| 俄罗斯特黄特色一大片| 女人被狂操c到高潮| 熟妇人妻久久中文字幕3abv| 神马国产精品三级电影在线观看 | 久久这里只有精品19| 久久影院123| 在线十欧美十亚洲十日本专区| av中文乱码字幕在线| 嫩草影院精品99| 中国美女看黄片| 亚洲精品国产色婷婷电影| 九色亚洲精品在线播放| 天天躁狠狠躁夜夜躁狠狠躁| www.999成人在线观看| 午夜精品在线福利| 精品第一国产精品| 日韩精品中文字幕看吧| 一级a爱片免费观看的视频| 欧美日韩乱码在线| 老司机午夜十八禁免费视频| 成人欧美大片| 亚洲欧美精品综合一区二区三区| 村上凉子中文字幕在线| 成人亚洲精品一区在线观看| 国产av又大| 黄色毛片三级朝国网站| 亚洲精品中文字幕在线视频| 精品电影一区二区在线| 婷婷精品国产亚洲av在线| 色综合亚洲欧美另类图片| 亚洲熟妇熟女久久| 日韩高清综合在线| 麻豆国产av国片精品| 久久久久久久久免费视频了| 日本黄色视频三级网站网址| 国产一区二区三区综合在线观看| 在线观看免费日韩欧美大片| 国产aⅴ精品一区二区三区波| 日韩中文字幕欧美一区二区| 在线天堂中文资源库| 精品国产一区二区久久| 一级毛片高清免费大全| 国产精品久久久久久亚洲av鲁大| 久久中文字幕人妻熟女| 可以免费在线观看a视频的电影网站| 午夜成年电影在线免费观看| 久久精品国产亚洲av高清一级| 欧美最黄视频在线播放免费| 亚洲成人国产一区在线观看| 一级,二级,三级黄色视频| 麻豆成人av在线观看| 一级a爱视频在线免费观看| 亚洲美女黄片视频| 一夜夜www| 在线av久久热| 国产成+人综合+亚洲专区| 女人高潮潮喷娇喘18禁视频| 极品教师在线免费播放| 最新在线观看一区二区三区| 亚洲一区中文字幕在线| 在线永久观看黄色视频| 又大又爽又粗| 淫妇啪啪啪对白视频| 免费在线观看影片大全网站| 免费高清在线观看日韩| 狂野欧美激情性xxxx| 欧美精品亚洲一区二区| 精品人妻在线不人妻| 在线观看免费日韩欧美大片| 午夜免费鲁丝| 日韩精品青青久久久久久| 一夜夜www| www国产在线视频色| av电影中文网址| 久久国产乱子伦精品免费另类| 91av网站免费观看| 可以免费在线观看a视频的电影网站| 国产精品国产高清国产av| 在线观看免费视频网站a站| 这个男人来自地球电影免费观看| 精品一区二区三区av网在线观看| 色哟哟哟哟哟哟| 国产三级黄色录像| 首页视频小说图片口味搜索| 欧美绝顶高潮抽搐喷水| 日韩欧美一区视频在线观看| 亚洲精华国产精华精| 精品午夜福利视频在线观看一区| 女性被躁到高潮视频| 国产真人三级小视频在线观看| 国产精品久久久久久人妻精品电影| 久久影院123| 妹子高潮喷水视频| 久久久水蜜桃国产精品网| 久久久久久亚洲精品国产蜜桃av| 9热在线视频观看99| 国产主播在线观看一区二区| 男人的好看免费观看在线视频 | 搡老岳熟女国产| 十八禁人妻一区二区| 啦啦啦韩国在线观看视频| 亚洲美女黄片视频| 免费看十八禁软件| 国产99久久九九免费精品| 一进一出抽搐gif免费好疼| 欧美日本中文国产一区发布| 国产麻豆69| www日本在线高清视频| 亚洲成av片中文字幕在线观看| 一区二区三区精品91| 在线观看舔阴道视频| 久久精品国产亚洲av香蕉五月| 国产欧美日韩一区二区精品| 成人永久免费在线观看视频| 国产av一区二区精品久久| 真人一进一出gif抽搐免费| 91在线观看av| 国产一卡二卡三卡精品| 国产熟女午夜一区二区三区| 久久精品亚洲精品国产色婷小说| 免费av毛片视频| 非洲黑人性xxxx精品又粗又长| 亚洲色图av天堂| 欧美激情极品国产一区二区三区| 国产精品亚洲av一区麻豆| 一边摸一边做爽爽视频免费| 久久久久久国产a免费观看| 午夜久久久在线观看| ponron亚洲| 俄罗斯特黄特色一大片| 成人18禁高潮啪啪吃奶动态图| 国产精品久久久久久亚洲av鲁大| 亚洲欧美一区二区三区黑人| 精品欧美一区二区三区在线| 国产午夜福利久久久久久| 日韩欧美免费精品| 国产日韩一区二区三区精品不卡| ponron亚洲| av超薄肉色丝袜交足视频| 一区二区三区国产精品乱码| 精品一品国产午夜福利视频| 别揉我奶头~嗯~啊~动态视频| 国产成人av激情在线播放| 村上凉子中文字幕在线| 制服诱惑二区| 麻豆一二三区av精品| 久久精品国产亚洲av香蕉五月| 美女 人体艺术 gogo| 国产乱人伦免费视频| or卡值多少钱| av天堂久久9| 亚洲少妇的诱惑av| cao死你这个sao货| 看片在线看免费视频| 免费在线观看完整版高清| 成人精品一区二区免费| 丝袜人妻中文字幕| 国产亚洲精品一区二区www| 久久精品国产亚洲av香蕉五月| 免费在线观看亚洲国产| 国产精品1区2区在线观看.| 成人免费观看视频高清| 在线播放国产精品三级| 999久久久国产精品视频| 69av精品久久久久久| 国产激情欧美一区二区| 亚洲精品久久成人aⅴ小说| 午夜老司机福利片| 国产精品99久久99久久久不卡| 亚洲av日韩精品久久久久久密| www.熟女人妻精品国产| 国产精品影院久久| tocl精华| 女性被躁到高潮视频| 久久精品人人爽人人爽视色| 亚洲avbb在线观看| 精品电影一区二区在线| 午夜福利在线观看吧| 91九色精品人成在线观看| 深夜精品福利| 成年女人毛片免费观看观看9| 少妇粗大呻吟视频| 久久天躁狠狠躁夜夜2o2o| 香蕉久久夜色| 岛国视频午夜一区免费看| 久久久久国产一级毛片高清牌| 国产欧美日韩一区二区三区在线| 午夜福利在线观看吧| 一边摸一边抽搐一进一出视频| 黄色视频,在线免费观看| 国产1区2区3区精品| 亚洲av美国av| 午夜福利成人在线免费观看| 99国产极品粉嫩在线观看| 国产成人影院久久av| 国产精品乱码一区二三区的特点 | 午夜久久久在线观看| 亚洲av成人不卡在线观看播放网| 男人的好看免费观看在线视频 | 一区二区三区精品91| 国产精品国产高清国产av| 久久精品人人爽人人爽视色| 久久人人爽av亚洲精品天堂| 啦啦啦免费观看视频1| 日本一区二区免费在线视频| 成人特级黄色片久久久久久久| 电影成人av| 99国产精品一区二区三区| 日韩三级视频一区二区三区| 国产精品永久免费网站| 搡老岳熟女国产| 极品人妻少妇av视频| 妹子高潮喷水视频| 国产精品久久久久久人妻精品电影| 女人爽到高潮嗷嗷叫在线视频| 黄色视频不卡| 久久国产精品人妻蜜桃| 一a级毛片在线观看| 大码成人一级视频| 97碰自拍视频| 如日韩欧美国产精品一区二区三区| 欧美日韩精品网址| 性色av乱码一区二区三区2| 搡老熟女国产l中国老女人| 国产三级黄色录像| 久久青草综合色| 日本精品一区二区三区蜜桃| 久久伊人香网站| 精品国产一区二区三区四区第35| 精品无人区乱码1区二区| 亚洲全国av大片| 欧美日韩一级在线毛片| 色哟哟哟哟哟哟| 亚洲人成电影观看| 激情在线观看视频在线高清| 成在线人永久免费视频| 成人精品一区二区免费| 欧美日本中文国产一区发布| 国产精品一区二区免费欧美| 国产精品久久久久久精品电影 | 国产成人一区二区三区免费视频网站| 精品国产亚洲在线| av超薄肉色丝袜交足视频| 亚洲av五月六月丁香网| 亚洲aⅴ乱码一区二区在线播放 | 午夜福利视频1000在线观看 | 黑丝袜美女国产一区| 亚洲午夜理论影院| 久久久久久国产a免费观看| 日日摸夜夜添夜夜添小说| 亚洲成人国产一区在线观看| 亚洲精品国产色婷婷电影| 在线永久观看黄色视频| 日韩欧美在线二视频| 国产xxxxx性猛交| 成人三级做爰电影| 午夜福利欧美成人| 自拍欧美九色日韩亚洲蝌蚪91| 国产在线观看jvid| 精品日产1卡2卡| 亚洲第一av免费看| 国产精品亚洲av一区麻豆| 久久久国产欧美日韩av| 亚洲免费av在线视频| 乱人伦中国视频| 国产欧美日韩综合在线一区二区| 999精品在线视频| 久久九九热精品免费| 精品国产亚洲在线| 禁无遮挡网站| 亚洲精品在线美女| 禁无遮挡网站| 亚洲中文字幕一区二区三区有码在线看 | 激情视频va一区二区三区| 亚洲国产欧美一区二区综合| 岛国在线观看网站| 欧美成人性av电影在线观看| 欧美日韩瑟瑟在线播放| 大型av网站在线播放| 亚洲va日本ⅴa欧美va伊人久久| 精品国产乱码久久久久久男人| 欧美绝顶高潮抽搐喷水| 国产高清有码在线观看视频 | 18禁国产床啪视频网站| 亚洲第一电影网av| 欧美午夜高清在线| svipshipincom国产片| 日韩精品免费视频一区二区三区| 精品日产1卡2卡| 亚洲国产看品久久| 色av中文字幕| 午夜福利欧美成人| 色尼玛亚洲综合影院| 两个人免费观看高清视频| 天堂影院成人在线观看| 国产精品免费视频内射| 露出奶头的视频| 99精品在免费线老司机午夜| 中文字幕人妻丝袜一区二区| 夜夜躁狠狠躁天天躁| 长腿黑丝高跟| 午夜免费鲁丝| 99在线人妻在线中文字幕| 一级毛片精品| videosex国产| 亚洲午夜理论影院| 日本a在线网址| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩亚洲综合一区二区三区_| av电影中文网址| 日韩视频一区二区在线观看| 91九色精品人成在线观看| 久久精品国产亚洲av高清一级| 又黄又粗又硬又大视频| 欧美在线一区亚洲| 变态另类成人亚洲欧美熟女 | 老司机深夜福利视频在线观看| 亚洲九九香蕉| 一本大道久久a久久精品| 国产成人精品久久二区二区91| 777久久人妻少妇嫩草av网站| 69av精品久久久久久| 国产精品亚洲一级av第二区| 夜夜爽天天搞| 九色亚洲精品在线播放| 在线观看午夜福利视频| 日本黄色视频三级网站网址| 欧美色视频一区免费| 日韩免费av在线播放| 精品熟女少妇八av免费久了| 伊人久久大香线蕉亚洲五| 桃色一区二区三区在线观看| 久久狼人影院| 国内久久婷婷六月综合欲色啪| 欧美黑人欧美精品刺激| 亚洲精品av麻豆狂野| 国产一区二区三区在线臀色熟女|