• <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年來全球背景下青藏高原氣候變化特征分析
    基于小波變換的三江平原旬降水量主周期識別
    国产片内射在线| 麻豆一二三区av精品| 国产97色在线日韩免费| 一本久久中文字幕| 天堂√8在线中文| 纯流量卡能插随身wifi吗| 又黄又粗又硬又大视频| www.精华液| а√天堂www在线а√下载| 亚洲精品美女久久久久99蜜臀| 在线免费观看的www视频| 亚洲av电影在线进入| 日韩视频一区二区在线观看| 免费搜索国产男女视频| 可以在线观看的亚洲视频| 国产精品亚洲av一区麻豆| av片东京热男人的天堂| 免费在线观看黄色视频的| 亚洲,欧美精品.| 亚洲精品国产区一区二| 好看av亚洲va欧美ⅴa在| 级片在线观看| 免费不卡黄色视频| 禁无遮挡网站| 免费观看人在逋| 亚洲激情在线av| 涩涩av久久男人的天堂| 国产乱人伦免费视频| 女人爽到高潮嗷嗷叫在线视频| 一级毛片精品| 精品国产乱子伦一区二区三区| 国产麻豆69| 男人的好看免费观看在线视频 | 午夜成年电影在线免费观看| 97超级碰碰碰精品色视频在线观看| 18禁裸乳无遮挡免费网站照片 | 12—13女人毛片做爰片一| av免费在线观看网站| 咕卡用的链子| 日本三级黄在线观看| 人人妻人人澡人人看| 久久午夜亚洲精品久久| 十八禁人妻一区二区| 日本a在线网址| 免费观看人在逋| av有码第一页| 亚洲男人天堂网一区| 1024香蕉在线观看| 国产免费av片在线观看野外av| 国产真人三级小视频在线观看| 伊人久久大香线蕉亚洲五| 纯流量卡能插随身wifi吗| 精品国产亚洲在线| 亚洲一区二区三区色噜噜| 久久久久九九精品影院| 亚洲精品美女久久久久99蜜臀| 国产精品二区激情视频| 精品国产一区二区三区四区第35| 人人妻人人澡欧美一区二区 | 亚洲久久久国产精品| 国产真人三级小视频在线观看| 午夜免费成人在线视频| 国产欧美日韩一区二区三区在线| 妹子高潮喷水视频| 中文亚洲av片在线观看爽| 国产精品98久久久久久宅男小说| 身体一侧抽搐| 日韩精品中文字幕看吧| 99国产精品一区二区三区| 成人18禁在线播放| av福利片在线| 国产一卡二卡三卡精品| 成人国产综合亚洲| 丰满人妻熟妇乱又伦精品不卡| 久久久久九九精品影院| 老汉色av国产亚洲站长工具| 欧美一区二区精品小视频在线| 中文字幕高清在线视频| 18禁黄网站禁片午夜丰满| 日韩视频一区二区在线观看| 成人三级黄色视频| 国产精品98久久久久久宅男小说| 级片在线观看| 亚洲欧美激情综合另类| 欧美日韩亚洲综合一区二区三区_| 午夜激情av网站| 操美女的视频在线观看| 久久性视频一级片| 亚洲国产精品成人综合色| 国产精品av久久久久免费| 又大又爽又粗| 成人手机av| 日韩有码中文字幕| 丰满人妻熟妇乱又伦精品不卡| 99精品在免费线老司机午夜| 欧美绝顶高潮抽搐喷水| 97人妻天天添夜夜摸| 嫩草影视91久久| 亚洲片人在线观看| or卡值多少钱| 免费看美女性在线毛片视频| 久久久水蜜桃国产精品网| 亚洲av美国av| 国产精品久久久久久精品电影 | 免费在线观看影片大全网站| 侵犯人妻中文字幕一二三四区| 亚洲成人久久性| 精品电影一区二区在线| 久久人人97超碰香蕉20202| 成人三级黄色视频| 久久婷婷人人爽人人干人人爱 | 视频在线观看一区二区三区| 国产三级黄色录像| avwww免费| 999久久久国产精品视频| 激情在线观看视频在线高清| 又紧又爽又黄一区二区| 色哟哟哟哟哟哟| 少妇粗大呻吟视频| avwww免费| 国产三级黄色录像| 午夜亚洲福利在线播放| 国产精华一区二区三区| 黄片大片在线免费观看| 99国产精品一区二区蜜桃av| 久久人妻福利社区极品人妻图片| 深夜精品福利| 悠悠久久av| 黄色a级毛片大全视频| 午夜精品久久久久久毛片777| 91av网站免费观看| 99re在线观看精品视频| 日韩欧美国产一区二区入口| 女人被躁到高潮嗷嗷叫费观| 精品不卡国产一区二区三区| 午夜免费鲁丝| 黄频高清免费视频| 欧美日韩黄片免| 精品国内亚洲2022精品成人| 大码成人一级视频| 成人18禁在线播放| 黄色a级毛片大全视频| 久久国产亚洲av麻豆专区| 黄色毛片三级朝国网站| 狂野欧美激情性xxxx| 狠狠狠狠99中文字幕| 成人国语在线视频| 18禁黄网站禁片午夜丰满| 一区在线观看完整版| 欧美人与性动交α欧美精品济南到| 久久久久久亚洲精品国产蜜桃av| 日日爽夜夜爽网站| 成年版毛片免费区| 搡老岳熟女国产| 国产精品日韩av在线免费观看 | 国产片内射在线| 在线观看舔阴道视频| 欧美日韩黄片免| 日韩欧美三级三区| 久久精品91无色码中文字幕| 中出人妻视频一区二区| 一级毛片女人18水好多| 亚洲国产精品合色在线| 天天躁夜夜躁狠狠躁躁| 人人妻人人爽人人添夜夜欢视频| 国产高清有码在线观看视频 | 成年版毛片免费区| 国产男靠女视频免费网站| 女人被躁到高潮嗷嗷叫费观| 国产精品免费视频内射| 亚洲专区字幕在线| 黄色毛片三级朝国网站| 午夜久久久久精精品| 啦啦啦观看免费观看视频高清 | 一个人免费在线观看的高清视频| 欧美在线一区亚洲| 国产成人av教育| 黄色a级毛片大全视频| 精品少妇一区二区三区视频日本电影| 久久婷婷人人爽人人干人人爱 | 色播在线永久视频| 999精品在线视频| 国产不卡一卡二| 变态另类丝袜制服| 欧美日韩中文字幕国产精品一区二区三区 | 久久国产精品人妻蜜桃| 可以在线观看毛片的网站| 一级黄色大片毛片| 男女做爰动态图高潮gif福利片 | 国产av又大| 少妇的丰满在线观看| 中文字幕av电影在线播放| 一边摸一边做爽爽视频免费| 亚洲精品久久成人aⅴ小说| 日本免费a在线| 又大又爽又粗| 亚洲avbb在线观看| 亚洲午夜理论影院| 色老头精品视频在线观看| 国产av在哪里看| 欧美黑人精品巨大| av有码第一页| 亚洲av成人av| 国产高清激情床上av| 久热这里只有精品99| 宅男免费午夜| 日韩三级视频一区二区三区| 久久中文看片网| 国产xxxxx性猛交| 国产av在哪里看| 欧美+亚洲+日韩+国产| 大型黄色视频在线免费观看| 大型黄色视频在线免费观看| 久久欧美精品欧美久久欧美| 一进一出抽搐gif免费好疼| 久久婷婷成人综合色麻豆| 精品国产超薄肉色丝袜足j| 黄网站色视频无遮挡免费观看| 中文亚洲av片在线观看爽| 神马国产精品三级电影在线观看 | 欧美黑人精品巨大| 人人妻人人爽人人添夜夜欢视频| videosex国产| 如日韩欧美国产精品一区二区三区| 麻豆国产av国片精品| 国产97色在线日韩免费| 怎么达到女性高潮| 亚洲av五月六月丁香网| 少妇粗大呻吟视频| 露出奶头的视频| 三级毛片av免费| 最近最新中文字幕大全电影3 | 国产精品九九99| 一夜夜www| 国产成人影院久久av| 日本vs欧美在线观看视频| 级片在线观看| 一级毛片高清免费大全| 老汉色∧v一级毛片| 婷婷六月久久综合丁香| 高清黄色对白视频在线免费看| 亚洲人成电影观看| 国产极品粉嫩免费观看在线| 国内精品久久久久精免费| 亚洲午夜精品一区,二区,三区| 国产亚洲精品第一综合不卡| 亚洲精品粉嫩美女一区| 黑人巨大精品欧美一区二区蜜桃| 国产高清视频在线播放一区| 午夜免费观看网址| 国产成人精品久久二区二区91| 一进一出抽搐动态| 国产人伦9x9x在线观看| 老司机深夜福利视频在线观看| 久久婷婷成人综合色麻豆| 精品一品国产午夜福利视频| 少妇裸体淫交视频免费看高清 | 国产精品久久电影中文字幕| 黄色毛片三级朝国网站| 国产亚洲av高清不卡| 国产一区在线观看成人免费| 变态另类丝袜制服| 亚洲国产欧美一区二区综合| 操美女的视频在线观看| 黄片播放在线免费| 亚洲色图av天堂| 婷婷精品国产亚洲av在线| 中文字幕人妻丝袜一区二区| 男女之事视频高清在线观看| 人人妻人人澡人人看| 成人18禁高潮啪啪吃奶动态图| 曰老女人黄片| 国产男靠女视频免费网站| 午夜福利免费观看在线| 18禁裸乳无遮挡免费网站照片 | 久久精品国产清高在天天线| 久久久国产精品麻豆| 99riav亚洲国产免费| 法律面前人人平等表现在哪些方面| 中文亚洲av片在线观看爽| 国产亚洲欧美在线一区二区| 久久久精品国产亚洲av高清涩受| 可以在线观看毛片的网站| 亚洲成人免费电影在线观看| 久久精品aⅴ一区二区三区四区| 国内毛片毛片毛片毛片毛片| 禁无遮挡网站| 色综合欧美亚洲国产小说| 无限看片的www在线观看| 亚洲av成人不卡在线观看播放网| 黄色视频,在线免费观看| 一级作爱视频免费观看| 黑人欧美特级aaaaaa片| 一边摸一边做爽爽视频免费| 丰满的人妻完整版| 嫩草影院精品99| 中文字幕人妻熟女乱码| 最新美女视频免费是黄的| www日本在线高清视频| 18禁观看日本| 欧美日本亚洲视频在线播放| 亚洲av日韩精品久久久久久密| 99在线视频只有这里精品首页| 日日爽夜夜爽网站| 99久久久亚洲精品蜜臀av| 亚洲片人在线观看| 国产成人精品久久二区二区91| 桃色一区二区三区在线观看| 精品久久久精品久久久| 欧美日本亚洲视频在线播放| 88av欧美| 国产色视频综合| 国产精品一区二区在线不卡| 欧美日韩精品网址| 国产一区在线观看成人免费| 亚洲精品中文字幕一二三四区| 中文字幕色久视频| 亚洲专区字幕在线| 97人妻精品一区二区三区麻豆 | 亚洲在线自拍视频| 老司机在亚洲福利影院| 午夜免费观看网址| 18禁黄网站禁片午夜丰满| 97超级碰碰碰精品色视频在线观看| 日本在线视频免费播放| 国产熟女xx| 欧美日韩黄片免| 亚洲精品国产区一区二| 国产精品久久久av美女十八| 国产主播在线观看一区二区| 亚洲av熟女| 午夜久久久在线观看| 一本大道久久a久久精品| 精品人妻1区二区| 亚洲熟妇中文字幕五十中出| 制服诱惑二区| 久久国产亚洲av麻豆专区| 亚洲伊人色综图| 琪琪午夜伦伦电影理论片6080| 日本vs欧美在线观看视频| 中文字幕av电影在线播放| 日韩精品青青久久久久久| 国产亚洲欧美在线一区二区| 久9热在线精品视频| 丰满人妻熟妇乱又伦精品不卡| 老司机深夜福利视频在线观看| 久久人人97超碰香蕉20202| 国产亚洲欧美在线一区二区| 99在线视频只有这里精品首页| 淫妇啪啪啪对白视频| 国产精品久久久久久亚洲av鲁大| 性色av乱码一区二区三区2| 精品久久久久久久久久免费视频| 日韩成人在线观看一区二区三区| 精品高清国产在线一区| 一区二区日韩欧美中文字幕| 欧美人与性动交α欧美精品济南到| 看免费av毛片| 人人妻,人人澡人人爽秒播| а√天堂www在线а√下载| 黄色丝袜av网址大全| 亚洲国产欧美日韩在线播放| 大陆偷拍与自拍| www.精华液| 亚洲一码二码三码区别大吗| 亚洲男人的天堂狠狠| 久久久久久免费高清国产稀缺| 成人亚洲精品一区在线观看| 精品午夜福利视频在线观看一区| 又紧又爽又黄一区二区| 日韩三级视频一区二区三区| 在线av久久热| 精品福利观看| 亚洲久久久国产精品| x7x7x7水蜜桃| 777久久人妻少妇嫩草av网站| 国产精品一区二区三区四区久久 | 琪琪午夜伦伦电影理论片6080| 成年女人毛片免费观看观看9| 亚洲人成网站在线播放欧美日韩| 亚洲熟妇熟女久久| 亚洲免费av在线视频| 一二三四在线观看免费中文在| 久久久水蜜桃国产精品网| 欧美丝袜亚洲另类 | 国产人伦9x9x在线观看| 亚洲成人久久性| 黄频高清免费视频| 国产日韩一区二区三区精品不卡| 97碰自拍视频| 色播亚洲综合网| 精品国内亚洲2022精品成人| 日韩精品免费视频一区二区三区| aaaaa片日本免费| 不卡一级毛片| 国产亚洲精品久久久久5区| 久久人妻av系列| 日韩精品青青久久久久久| 久久中文字幕人妻熟女| 日本在线视频免费播放| 亚洲 国产 在线| 电影成人av| 男男h啪啪无遮挡| 一区二区日韩欧美中文字幕| 欧美日韩福利视频一区二区| 国产精品乱码一区二三区的特点 | 熟女少妇亚洲综合色aaa.| 国产成人欧美| 国产成人精品久久二区二区91| 啦啦啦观看免费观看视频高清 | 两人在一起打扑克的视频| 变态另类成人亚洲欧美熟女 | 欧美久久黑人一区二区| 亚洲美女黄片视频| 一级,二级,三级黄色视频| 波多野结衣巨乳人妻| 一区福利在线观看| 可以在线观看的亚洲视频| 欧美不卡视频在线免费观看 | 国产亚洲精品一区二区www| 亚洲成人国产一区在线观看| 日韩一卡2卡3卡4卡2021年| 精品福利观看| 一个人免费在线观看的高清视频| 色尼玛亚洲综合影院| 久久欧美精品欧美久久欧美| 久久天躁狠狠躁夜夜2o2o| 免费一级毛片在线播放高清视频 | 欧美黄色淫秽网站| 国产午夜福利久久久久久| 一夜夜www| 欧美最黄视频在线播放免费| 日韩精品免费视频一区二区三区| 在线观看免费午夜福利视频| 国语自产精品视频在线第100页| 女生性感内裤真人,穿戴方法视频| 欧美中文综合在线视频| 老司机靠b影院| 人成视频在线观看免费观看| 波多野结衣巨乳人妻| 亚洲视频免费观看视频| 国产97色在线日韩免费| 99在线人妻在线中文字幕| 欧美丝袜亚洲另类 | 国产成人精品无人区| 亚洲一区高清亚洲精品| 老司机在亚洲福利影院| 亚洲精品美女久久av网站| 国产色视频综合| 亚洲国产欧美网| 午夜免费观看网址| 女人被狂操c到高潮| 亚洲国产欧美日韩在线播放| 男女下面插进去视频免费观看| 99国产精品99久久久久| 精品国产美女av久久久久小说| 免费久久久久久久精品成人欧美视频| 久久久国产欧美日韩av| 精品无人区乱码1区二区| 一本大道久久a久久精品| or卡值多少钱| 国产aⅴ精品一区二区三区波| 欧美日本亚洲视频在线播放| av福利片在线| 国产精品99久久99久久久不卡| 欧美激情极品国产一区二区三区| 午夜老司机福利片| 熟妇人妻久久中文字幕3abv| 亚洲精品久久成人aⅴ小说| 18美女黄网站色大片免费观看| 精品久久久久久久久久免费视频| 午夜福利高清视频| 日韩欧美免费精品| 欧美中文综合在线视频| 亚洲精品在线美女| 久久久久久大精品| 成人三级做爰电影| 一区二区三区激情视频| 久久久水蜜桃国产精品网| 日韩一卡2卡3卡4卡2021年| 757午夜福利合集在线观看| av视频免费观看在线观看| 波多野结衣一区麻豆| 欧美日本视频| 欧美午夜高清在线| 桃色一区二区三区在线观看| 亚洲成人免费电影在线观看| 免费看美女性在线毛片视频| 制服人妻中文乱码| 成人av一区二区三区在线看| 欧美日韩乱码在线| 免费在线观看日本一区| 午夜福利欧美成人| 操出白浆在线播放| 九色亚洲精品在线播放| 嫩草影院精品99| 正在播放国产对白刺激| 91成人精品电影| 男人操女人黄网站| 亚洲色图av天堂| or卡值多少钱| 国产精品一区二区免费欧美| 美女免费视频网站| 99久久国产精品久久久| 亚洲欧美激情在线| 精品午夜福利视频在线观看一区| 日本黄色视频三级网站网址| 最近最新中文字幕大全免费视频| 麻豆一二三区av精品| 精品乱码久久久久久99久播| 欧美成人免费av一区二区三区| 91字幕亚洲| 禁无遮挡网站| 亚洲av五月六月丁香网| 国产精品自产拍在线观看55亚洲| 国内精品久久久久精免费| 男女做爰动态图高潮gif福利片 | 国产成人免费无遮挡视频| 国产免费av片在线观看野外av| 在线av久久热| 大型黄色视频在线免费观看| 嫁个100分男人电影在线观看| 欧美日韩瑟瑟在线播放| 免费高清在线观看日韩| 久久青草综合色| 国产私拍福利视频在线观看| 午夜福利成人在线免费观看| www.自偷自拍.com| 制服诱惑二区| 欧美成人午夜精品| 欧美在线一区亚洲| 国产免费男女视频| 亚洲狠狠婷婷综合久久图片| 日日夜夜操网爽| 88av欧美| 级片在线观看| 日本 av在线| 自拍欧美九色日韩亚洲蝌蚪91| av天堂在线播放| 国产极品粉嫩免费观看在线| 欧美日韩黄片免| 亚洲人成电影观看| 一级a爱视频在线免费观看| а√天堂www在线а√下载| 亚洲国产日韩欧美精品在线观看 | 在线十欧美十亚洲十日本专区| 亚洲av电影在线进入| 一级毛片精品| 99在线人妻在线中文字幕| 两个人看的免费小视频| 丰满的人妻完整版| 欧美中文综合在线视频| 国产一区在线观看成人免费| 性欧美人与动物交配| 女人高潮潮喷娇喘18禁视频| 黄色片一级片一级黄色片| 国产精品综合久久久久久久免费 | 国产aⅴ精品一区二区三区波| 精品人妻在线不人妻| 免费久久久久久久精品成人欧美视频| 成人三级黄色视频| 亚洲一区高清亚洲精品| 人人妻人人澡欧美一区二区 | 少妇粗大呻吟视频| 久9热在线精品视频| 国产成人欧美| 99在线视频只有这里精品首页| 亚洲中文字幕一区二区三区有码在线看 | 久热这里只有精品99| 天堂影院成人在线观看| 亚洲片人在线观看| 法律面前人人平等表现在哪些方面| av在线播放免费不卡| 欧美另类亚洲清纯唯美| avwww免费| 亚洲国产精品久久男人天堂| 国产高清激情床上av| 日本在线视频免费播放| 女人高潮潮喷娇喘18禁视频| 99riav亚洲国产免费| 国产精品亚洲美女久久久| 亚洲男人天堂网一区| 久久草成人影院| 久久久久久人人人人人| 久久亚洲真实| 韩国av一区二区三区四区| 精品乱码久久久久久99久播| 狂野欧美激情性xxxx| 午夜视频精品福利| 国产精品久久视频播放| 女人爽到高潮嗷嗷叫在线视频| 免费在线观看完整版高清| 大码成人一级视频| 天天躁夜夜躁狠狠躁躁| 国产成人影院久久av| 少妇熟女aⅴ在线视频| 女性生殖器流出的白浆| 亚洲欧洲精品一区二区精品久久久| 免费看a级黄色片| 麻豆国产av国片精品| 国产精品1区2区在线观看.| 天天添夜夜摸| 人妻丰满熟妇av一区二区三区| 男人舔女人下体高潮全视频| 精品国产一区二区久久| 亚洲国产高清在线一区二区三 | 亚洲av片天天在线观看| 1024香蕉在线观看| 精品日产1卡2卡| 人人妻,人人澡人人爽秒播| 久久九九热精品免费| 777久久人妻少妇嫩草av网站| 极品教师在线免费播放| 亚洲色图综合在线观看| 国产91精品成人一区二区三区|