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

    CLDAS融合土壤濕度產(chǎn)品在東北地區(qū)的適用性評估及訂正*

    2019-10-18 09:03:32崔園園敬文琪
    中國農(nóng)業(yè)氣象 2019年10期
    關(guān)鍵詞:陸面土壤濕度適用性

    崔園園,張 強,覃 軍,敬文琪

    CLDAS融合土壤濕度產(chǎn)品在東北地區(qū)的適用性評估及訂正*

    崔園園1,張 強2**,覃 軍3,敬文琪4

    (1. 河北省氣象臺,石家莊 050000; 2. 中國氣象局國家氣候中心氣候研究開放實驗室,北京 100081;3. 中國地質(zhì)大學(xué)(武漢)環(huán)境學(xué)院,武漢 430074; 4. 94758部隊氣象臺,福建寧德 355103)

    通過計算CLDAS的0?20cm土壤相對濕度融合產(chǎn)品與農(nóng)業(yè)氣象觀測站逐旬10cm土壤相對濕度觀測值的相關(guān)性和偏差,綜合分析CLDAS在東北地區(qū)的適用性。結(jié)果表明:CLDAS融合土壤相對濕度能較合理地反映東北地區(qū)10cm深度土壤相對濕度的空間分布,但是在對極值的描述上有所欠缺,其中CLDAS對遼寧西部和黑龍江西南部的土壤相對濕度模擬偏高較明顯,而在黑龍江的東北部地區(qū)對實際土壤相對濕度的模擬偏低;總體上,CLDAS融合土壤相對濕度資料在東北地區(qū)的適用性由西南向東北方向遞減,在吉林和遼寧省的適用性更好?;貧w訂正和7旬滑動平均訂正法對CLDAS產(chǎn)品的誤差訂正表明,訂正結(jié)果大大提高了CLDAS產(chǎn)品的精度,其訂正后產(chǎn)品可以應(yīng)用于土壤濕度干旱監(jiān)測業(yè)務(wù)。

    東北地區(qū);CLDAS融合土壤相對濕度產(chǎn)品;回歸訂正;7旬滑動平均訂正;適用性評估

    在全球氣候變暖的背景影響下,地表升溫,水資源的時空配置也相應(yīng)發(fā)生變化,大范圍干旱事件頻發(fā),與洪澇災(zāi)害相比,干旱不僅發(fā)展緩慢而且具有較長的時間尺度,更難被察覺和監(jiān)控[1?3]。

    目前,主要采取土壤水分自動觀測站和衛(wèi)星遙感兩種手段來對土壤干旱進(jìn)行監(jiān)測[4]。由于觀測站點東西部地區(qū)分布不均勻、觀測儀器本身的安裝校正以及季節(jié)等客觀因素的影響,觀測數(shù)據(jù)的質(zhì)量還存在各種問題,無法單純依靠站點數(shù)據(jù)對全國范圍的土壤墑情變化開展精細(xì)、準(zhǔn)確、有效的監(jiān)測。光學(xué)遙感也受到晝夜變化、大氣云雨遮擋的影響無法保證長時間序列、大面積實時監(jiān)測[5];而且微波遙感信號也局限于反演0?5cm深度表層土壤濕度情況,得到的土壤濕度空間準(zhǔn)確率和分辨率亟需提高。

    隨著資料同化技術(shù)的發(fā)展,產(chǎn)生了一系列陸面融合資料,用來描述陸面物理量的時空演變。目前國際上相對成熟的陸面數(shù)據(jù)同化系統(tǒng)包括:北美陸面數(shù)據(jù)同化系統(tǒng)NLDAS(North American Land Data Assimilation System)[6],歐洲陸面數(shù)據(jù)同化系統(tǒng)ELDAS(European Land Data Assimilation System)[7]和全球陸面數(shù)據(jù)同化系統(tǒng)GLDAS(Global Land Data Assimilation System)[8]。國內(nèi)李新等[9]建立了西部陸面數(shù)據(jù)同化系統(tǒng)WCLDAS(West China Land DataAssimilation System),師春香[10]發(fā)展了中國氣象局陸面數(shù)據(jù)同化系統(tǒng)CLDAS(CMA Land Data Assimilation System)。近年一些學(xué)者也對不同的陸面模式資料的適用性進(jìn)行了評估。Chen等[11]利用青藏高原中部土壤濕度觀測網(wǎng)資料對GLDAS中4個陸面模式土壤水分產(chǎn)品進(jìn)行了評估,研究表明GLDAS-CLM(Global Land Data Assimilation System- Common Land Model)和GLDAS-NOAH(Global Land Data Assimilation System-NCEP,OSU,Air Force and Office of Hydrology)模式資料則對10?40cm土壤層濕度具有較好的代表性。韓帥等[12]通過對CLDAS陸面融合產(chǎn)品的評估發(fā)現(xiàn)CLDAS土壤濕度產(chǎn)品精度較高,對中國區(qū)域干旱監(jiān)測具有重要的實用價值。崔園園等[13?14]在分析GLDAS和CLDAS在青藏高原的適用性時指出,CLDAS在對四川東南地區(qū)和新疆地區(qū)大部分站點的土壤濕度描述能力上有了較大的改善。

    自21世紀(jì)初,東北已經(jīng)成為干旱發(fā)生頻率較高的地區(qū)之一,損失嚴(yán)重、影響范圍廣、持續(xù)時間長的干旱災(zāi)情頻繁發(fā)生,形成了跨月、跨季節(jié)甚至跨年度的不同時間尺度的干旱災(zāi)害,因此,準(zhǔn)確預(yù)測和評估東北地區(qū)干旱的發(fā)生發(fā)展規(guī)律具有重要的現(xiàn)實意義[15?18]。但經(jīng)過同化技術(shù)得到的資料存在系統(tǒng)誤差,在不同地區(qū)的適用性即精度無法保證,需要對其合理性和適用性進(jìn)行評價及訂正[19]。本研究利用中國氣象局國家氣象信息中心提供的2008?2017年CLDAS融合土壤相對濕度產(chǎn)品,并結(jié)合農(nóng)業(yè)氣象觀測站點逐旬土壤相對濕度定量評估CLDAS在東北地區(qū)的適用性,再利用多種回歸方法對CLDAS融合土壤相對濕度產(chǎn)品進(jìn)行訂正,以期得到一套訂正后的CLDAS融合土壤相對濕度產(chǎn)品,能夠應(yīng)用于農(nóng)業(yè)干旱監(jiān)測,提高干旱的監(jiān)測精度,為農(nóng)業(yè)干旱防御提供服務(wù)。

    1 資料與方法

    1.1 資料

    1.1.1 觀測資料

    經(jīng)初步檢查,國家氣象信息中心整理的2008?2017年東北地區(qū)各層土壤相對濕度逐旬觀測數(shù)據(jù)中,10cm土層的資料具有較好的完整性和可靠性,且對農(nóng)作物的生長影響較大,農(nóng)業(yè)干旱發(fā)生的指示意義明顯[19],因此,選取該層土壤相對濕度作為研究對象。在對站點觀測進(jìn)行質(zhì)量控制后,選取東北地區(qū)共230個站點進(jìn)行分析。

    1.1.2 CMA陸面數(shù)據(jù)同化數(shù)據(jù)(CLDAS)

    選取國家氣象信息中心發(fā)布的中國氣象局陸面同化系統(tǒng)第二版本(CLDAS-V2.0)土壤相對濕度產(chǎn)品,CLDAS-V2.0產(chǎn)品為覆蓋亞洲區(qū)域(0?65°N,60?160°E),0.0625°×0.0625°(空間分辨率為)與1h(時間分辨率)的等經(jīng)緯度網(wǎng)格融合分析產(chǎn)品,該數(shù)據(jù)集利用多種來源地面、衛(wèi)星等觀測資料,采用多重網(wǎng)格變分同化(STMAS)、最優(yōu)插值(OI)、概率密度函數(shù)匹配(CDF)、物理反演、地形校正等技術(shù)研制而成,在中國區(qū)域質(zhì)量優(yōu)于國際同類產(chǎn)品,且時空分辨率更高。實時更新發(fā)布滯后1h。

    1.2 數(shù)據(jù)處理

    1.2.1 觀測資料質(zhì)量控制

    由于各站觀測資料質(zhì)量參差不齊,在使用前需要對土壤相對濕度資料進(jìn)行質(zhì)量控制,主要包括以下兩步:(1)考慮到觀測場地不會出現(xiàn)土壤中沒有水分或長時間不變的情況,以及在水中浸泡時,傳感器測量不準(zhǔn),所以剔除觀測資料中土壤相對濕度等于0、長時間維持同一個值,及大于100%的數(shù)據(jù);(2)觀測站點中存在灌溉的情況,無降水發(fā)生土壤濕度突然增大明顯,從觀測資料和CLDAS融合產(chǎn)品的相對誤差的概率密度分布圖上看出,90%以上站點的相對誤差絕對值<70%,故剔除相對誤差絕對值>70%的匹配樣本。

    1.2.2 CLDAS融合數(shù)據(jù)與實測數(shù)據(jù)的匹配

    (1)時間尺度

    在分析CLDAS融合土壤相對濕度產(chǎn)品與站點資料之間的誤差時,考慮到站點觀測資料以旬為時間尺度,故將CLDAS產(chǎn)品按旬作平均,后續(xù)與站點觀測資料進(jìn)行分析對比。

    (2)空間尺度

    在分析誤差的空間分布時,采用雙線性插值法將CLDAS融合相對濕度產(chǎn)品插值到東北地區(qū)230個站點上。

    (3)層次匹配

    由于CLDAS融合產(chǎn)品中的土壤相對濕度是一定厚度土壤層的平均值,在與觀測值進(jìn)行對比時,將CLDAS融合產(chǎn)品的0?20cm土壤相對濕度產(chǎn)品與10cm層次的觀測值進(jìn)行對比。

    1.3 適用性評估方法

    采用的統(tǒng)計量和評估指標(biāo)包括CLDAS產(chǎn)品與觀測資料的平均偏差BIAS、相關(guān)系數(shù)R,以及基于空間技巧評分復(fù)合統(tǒng)計量?時間技巧評分(ST)[20]。

    (1)平均偏差BIAS

    (2)相關(guān)系數(shù)R

    (3)時間技巧評分

    設(shè)Yt和Xt分別為CLDAS融合土壤相對濕度資料和站點觀測土壤濕度資料的時間序列,則CLDAS和觀測序列的平方誤差定義為

    式中,n為樣本時間長度,對式(3)進(jìn)行無量綱化,則無量綱時間技巧評分為

    1.4 訂正方法介紹

    1.4.1 回歸訂正法

    利用回歸訂正方法對CLDAS融合土壤相對濕度產(chǎn)品進(jìn)行訂正。針對東北地區(qū)的第i個站點,記第i個站點2008?2017年的觀測時間序列為MSRobs,相對應(yīng)的CLDAS融合土壤相對濕度資料的時間序列記為MSRCLDAS,建立回歸方程,即

    記MSRCLDAS_update為訂正后的CLDAS融合土壤相對濕度產(chǎn)品序列,則第i個站點訂正的融合土壤相對濕度序列為

    同理可得東北地區(qū)230個站點經(jīng)過回歸訂正后2008?2017年CLDAS融合相對濕度產(chǎn)品。

    1.4.2 7旬滑動平均訂正法

    將當(dāng)前時刻前7旬的誤差平均值(即前7旬CLDAS融合土壤相對濕度值與站點觀測值之差的平均值)作為訂正值,再用當(dāng)前時刻CLDAS融合土壤相對濕度值減去訂正值,即得該時刻訂正后的CLDAS產(chǎn)品,即

    式中,CLDASupdate(t)為第t個時刻經(jīng)過7旬滑動平均訂正法訂正后的CLDAS產(chǎn)品,CLDAS(i)、sta(i)分別為第t個時刻CLDAS融合土壤相對濕度值和站點觀測值。

    2 結(jié)果與分析

    2.1 CLDAS融合土壤相對濕度產(chǎn)品質(zhì)量評估

    2.1.1 土壤相對濕度的空間分布

    圖1給出了2008?2017年10cm深度各站點觀測值及CLDAS融合土壤相對濕度產(chǎn)品在東北地區(qū)的空間分布,由圖可知,CLDAS融合土壤相對濕度能較合理地反映土壤濕度在東北地區(qū)的空間分布,其中,10cm深度站點觀測值在東北地區(qū)的分布(圖1a)顯示土壤相對濕度基本在50%~95%,相對濕度的高值區(qū)主要位于遼寧中東部、吉林中東部以及黑龍江東北部地區(qū),土壤相對濕度的低值區(qū)主要位于遼寧西部、吉林西部以及黑龍江省的西南部地區(qū);CLDAS融合土壤相對濕度在0?20cm層次能較好地反映東北地區(qū)10cm土層土壤相對濕度的分布,但是在極值的描述上有所欠缺,主要表現(xiàn)為CLDAS對遼寧西部和黑龍江西南部的土壤相對濕度偏低區(qū)域模擬明顯偏高,以及在黑龍江的東北部地區(qū)對實際土壤相對濕度又有一個模擬偏低的表現(xiàn)。

    2.1.2 誤差的空間分布

    為了評估CLDAS融合土壤相對濕度在東北地區(qū)的適用性,根據(jù)式(1)?式(4)計算得出2008?2017年東北地區(qū)時間技巧評分(ST)、相關(guān)性(R)和平均偏差(BIAS)的空間分布(圖2、圖3)。

    一個站點的ST指標(biāo)反映的是該站點的觀測值與CLDAS兩種時間序列演變之間的相關(guān)性和偏差信息,ST值越接近于1,資料的精度越高,與站點的變化越一致。從圖2 可知,ST值由東北地區(qū)的西南向東北方向遞減,其中,黑龍江的西部及東北部地區(qū)是ST值的低值區(qū),其ST

    ST值反映的是相關(guān)性和偏差的綜合信息,為了進(jìn)一步分析CLDAS資料在東北地區(qū)ST值偏低時相關(guān)性和偏差的表現(xiàn),圖3給出了2008?2017年CLDAS融合土壤相對濕度與站點觀測之間的相關(guān)性和偏差分布圖,由圖可見,除了在黑龍江省的東北部地區(qū)CLDAS與站點觀測之間有較大偏差外(達(dá)到25%),大部分地區(qū)偏差均較小,基本位于±15%之間。

    總體而言,CLDAS土壤相對濕度資料在東北地區(qū)大部(主要是吉林和遼寧二省)適用性較好,黑龍江省適用性相對較差。

    圖1 2008?2017年東北地區(qū)各站點10cm土層土壤相對濕度觀測值(a)和0?20cm土層CLDAS融合平均土壤相對濕度(b)的空間分布

    圖2 2008?2017年CLDAS融合土壤相對濕度與站點觀測值的時間技巧評分(ST)的空間分布

    2.2 訂正結(jié)果評估

    利用回歸訂正法和7旬滑動平均訂正對CLDAS內(nèi)插到東北地區(qū)站點上的土壤相對濕度產(chǎn)品進(jìn)行訂正,分別得到訂正后的CLDAS土壤相對濕度旬資料,回歸訂正法部分站點的回歸系數(shù)見表1,CLDAS融合土壤相對濕度產(chǎn)品訂正后在東北地區(qū)的適用性分布見圖4?圖6,兩種訂正方法分別在東北三省的誤差統(tǒng)計見表2。

    由2008?2017年CLDAS土壤相對濕度經(jīng)過訂正后在東北地區(qū)的ST分布圖(圖4)可以看出,相對未訂正的CLDAS融合土壤濕度產(chǎn)品(圖4a),經(jīng)過7旬滑動平均訂正和回歸訂正后的ST評分由未訂正的?0.26分別增加至0.23和0.4,訂正后CLDAS土壤相對濕度ST評分在東北地區(qū)均得到了顯著增加,尤其是在黑龍江省,回歸訂正后的ST評分基本位于0~0.6之間,相比于7旬滑動平均法有顯著優(yōu)勢。結(jié)合表2分析可知,兩種訂正方法在吉林和遼寧省效果比較一致,ST值分別由訂正前的?0.21和?0.41提高到了0.3和0.4以上,其中7旬滑動平均訂正法稍顯優(yōu)勢,而對于黑龍江省,回歸訂正法的ST值明顯高于7旬滑動平均法,表明回歸訂正法效果更好。

    圖5和圖6分別為訂正后CLDAS融合土壤相對濕度與站點觀測的BIAS和R的空間分布圖,結(jié)合表2可知,訂正后CLDAS融合土壤濕度產(chǎn)品在東北地區(qū)的BIAS得到了顯著降低,7旬滑動平均訂正和回歸訂正后的BIAS由未訂正的6.3%分別下降到0.3%和1.3%,回歸訂正方法在消除系統(tǒng)偏差中優(yōu)勢明顯;而對于相關(guān)性,7旬滑動平均法比回歸訂正法效果更好。

    綜合來看,兩種訂正方法在東北地區(qū)均有較好的效果,在吉林和遼寧可以選擇7旬滑動平均訂正法,而黑龍江省選擇回歸訂正法可達(dá)到最優(yōu)的訂正效果;通過對訂正后的效果評估可以較好地提高CLDAS融合產(chǎn)品對東北地區(qū)土壤相對濕度的模擬精度,平均偏差(BIAS)在?0.5%~0.5%,誤差在允許范圍內(nèi),說明訂正后的資料能夠應(yīng)用于干旱監(jiān)測業(yè)務(wù),能有效提高土壤干旱監(jiān)測業(yè)務(wù)能力。

    圖3 2008?2017年CLDAS融合土壤相對濕度與站點觀測值的相關(guān)系數(shù)(R)和平均偏差(BIAS)的空間分布

    注:相關(guān)系數(shù)均通過0.05水平的顯著性檢驗。下同。

    Note: The R all passed the significance test at the level of 0.05. The same as below.

    表1 部分站點信息及其回歸訂正系數(shù)

    Table 1 The regression coefficient of partial sites

    表2 兩種訂正方法的誤差統(tǒng)計

    Table 2 Error statistics of the two revised methods

    注Note:J:吉林Jilin, L:遼寧Liaoning, H:黑龍江Heilongjiang, NE:東北地區(qū)Northeast of China.

    下同。The same as below.

    圖4 2008?2017年CLDAS融合土壤相對濕度與站點觀測值的時間技巧評分(ST)分布

    Fig. 4 Distribution of ST between soil relative moisture observation and CLDAS merged soil relative moisture during 2008?2017

    圖5 2008?2017年CLDAS融合土壤相對濕度與站點觀測之間的平均偏差(BIAS)分布 Fig. 5 Distribution of BIAS between soil relative moisture observation and CLDAS merged soil relative moisture during 2008?2017

    圖6 2008?2017年CLDAS融合土壤相對濕度與站點觀測之間的相關(guān)系數(shù)(R)分布 Fig. 6 Distribution of R between soil relative moisture observation and CLDAS merged soil relative moisture during 2008?2017

    3 結(jié)論與討論

    隨著各陸面數(shù)據(jù)同化系統(tǒng)的不斷深入,通過陸面過程模擬得到中國區(qū)域的土壤濕度成為可能[21?22]。國家氣象信息中心在CLDAS-1.0的基礎(chǔ)上,通過開發(fā)多陸面模式模擬、衛(wèi)星亮溫資料同化等,可進(jìn)一步提升土壤濕度的模擬精度[4]。但是相關(guān)學(xué)者指出同化產(chǎn)品的精度在很大程度上也取決于背景場的精度及分辨率,對于站點觀測資料稀少的地區(qū)適用性仍然較低[23],保鴻燕等[24]指出,土壤濕度的誤差在非凍結(jié)期主要由模式模擬的降水不準(zhǔn)確及降水對土壤濕度的強迫存在一定問題而引起,在凍結(jié)期主要由模式的凍融參數(shù)化方案導(dǎo)致對土壤的凍融過程模擬不準(zhǔn)確而引起,所以在后期應(yīng)用這些同化產(chǎn)品時,進(jìn)行誤差訂正是很有必要的。

    (1)CLDAS融合土壤相對濕度能較合理地反映土壤濕度在東北地區(qū)的實際空間分布,但是在對極值的描述上有所欠缺,表現(xiàn)在CLDAS對遼寧西部和黑龍江西南部的土壤相對濕度偏低區(qū)域模擬明顯偏高,以及在黑龍江的東北部地區(qū)對實際土壤相對濕度又偏低,說明CLDAS融合土壤相對濕度的空間誤差存在較大差異性。

    (2)由時間技巧評分的分布可知,CLDAS融合土壤相對濕度資料在東北地區(qū)的適用性由西南向東北方向遞減,在東北地區(qū)大部(主要是吉林和遼寧二?。┻m用性較好,黑龍江省適用性相對較差,說明融合資料與觀測的時間變化偏差各地區(qū)也存在差異,是否是凍融參數(shù)化方案導(dǎo)致,還需深入研究。

    (3)兩種訂正方法在東北地區(qū)均有較好的效果,在吉林和遼寧可以選擇7旬滑動平均訂正法,而黑龍江省選擇回歸訂正法可達(dá)到最優(yōu)的訂正效果。

    參考文獻(xiàn) References

    [1] 費振宇,孫宏巍,金菊良,等.近50年中國氣象干旱危險性的時空格局探討[J].水電能源科學(xué),2014,32(15):5-10.

    Fei Z Y,Sun H W,Jin J L,et alTemporal and spatial patterns of meteorological drought hazard in China for recent 50years[J]. Hydropower Energy Science,2014,32(15):5-10.(in Chinese)

    [2] 李憶平,李耀輝.氣象干旱指數(shù)在中國的適應(yīng)性研究進(jìn)展[J].干旱氣象,2017,35(5):709-723.

    Li Y P,Li Y H.Advances in adaptability of meteorological drought indices in China[J].Journal of Arid Meteorology, 2017,35(5):709-723.(in Chinese)

    [3] 徐一丹.基于四種干旱指數(shù)的東北地區(qū)多時間尺度干旱特征對比分析[D].沈陽:沈陽農(nóng)業(yè)大學(xué),2017.

    Xu Y D.Comparative analysis of multi-time scales drought characteristics based on four drought indices in Northeast China[D].Shenyang:Shenyang Agricultural University,2017. (in Chinese)

    [4] 韓帥,師春香,林泓錦,等.CLDAS土壤濕度業(yè)務(wù)產(chǎn)品的干旱監(jiān)測應(yīng)用[J].冰川凍土,2015,37(2):446-453.

    Han S,Shi C X,Lin J H,et al.The CLDAS soil moisture operation products applied to minitor soil drought[J].Journal of Glaciology and Geocryology,2015,37(2):446-453.(in Chinese)

    [5] 朱嫻韻,蘇布達(dá),黃金龍,等.云南氣候變化高分辨率模擬與RCP4.5情景預(yù)估[J].長江流域資源與環(huán)境,2015,24(3): 476-481.

    Zhu X Y,Su B D,Huang J L,et al.Simulation of climatic change in Yunnan Province and RCP4.5 scenario projected trend by CCLM[J].Resources and Environment in the Yangtze Basin,2015,24(3):476-481.(in Chinese)

    [6] Mitchell K E,Lohmann D,Houser P R,et al.The multi- institution North American Land Data Assimilation System (NLDAS):utilizing multiple GCIP products and partners in a continental distributed hydrological modeling system[J]. Journal of Geophysical Research,2004,109(D07):101-130.

    [7] Van den Hurk B,et al.Overview of the European Land Data Assimilation(ELDAS) project[J].EOS Trans,2002,83(47):H62D- 0888.

    [8] Rodell M,Houser P,Jambor U,et al.Status and availability of results from NASAs Global Land Data Assimilation System [J].EOS Trans,2002,83(19),B41A-04.

    [9] 李新,南卓銅,吳立宗,等.中國西部環(huán)境與生態(tài)科學(xué)數(shù)據(jù)中心:面向西部環(huán)境與生態(tài)科學(xué)的數(shù)據(jù)集成與共享[J]. 地球科學(xué)進(jìn)展,2008,23(6):628-637.

    Li X,Nan Z T,Wu L Z,et al.Environmental and Ecological Science Data Center for West China:integration and sharing of environmental and ecological data[J].Advances in Earth Science,2008,23(6):628-637.(in Chinese)

    [10] 師春香.基于ENKF算法的衛(wèi)星遙感土壤濕度同化研究[D].北京:中國科學(xué)院研究生院,2008.

    Shi C X.China land soil moisture ENKF data assimilation based on satellite remote sensing data[D].Beijing:Graduate School of Chinese Academy of Sciences,2008.(in Chinese)

    [11] Chen Y Y,Yang K,Qin J,et al.Evaluation of AMSR-Eretrievals and GLDAS simulations against observations of a soil moisture net work on the central Qinghai-Xizang Plateau [J].Journal of Geophysical Research:Atmospheres,2013,118 (10):4466-4475.

    [12] 韓帥,師春香,姜立鵬,等.CLDAS土壤濕度模擬結(jié)果及評估[J].應(yīng)用氣象學(xué)報,2017,28(3):369-378.

    Han S,Shi C X,Jiang L P,et al.The simulation and evaluation of soil moisture based on CLDAS[J].J Appl Meteor Sci,2017, 28(3): 369-378. (in Chinese)

    [13] 崔園園,覃軍,敬文琪,等.GLDAS和CLDAS融合土壤水分產(chǎn)品在青藏高原地區(qū)的適用性評估[J].高原氣象,2018,37 (1):123-136.

    Cui Y Y,Qin J,Jing W Q,et al.Applicability evaluation of merged soil moisture in GLDAS and CLDAS products over Qinghai-Tibetan[J].Plateau Meteorology,2018,37(1):123-136. (in Chinese)

    [14] 崔園園,敬文琪,覃軍.基于TIPEX Ⅲ資料對CLDAS-V2.0和GLDAS-NOAH陸面模式產(chǎn)品在青藏高原地區(qū)的適用性評估[J].高原氣象,2018,37(5):1143-1160.

    Cui Y Y,Jing W Q,Qin J.Applicability evaluation of merged soil moisture in GLDAS-NOAH and CLDAS-V2.0 products over Qinghai-Xizang Plateau of 2015 based on TIPEX Ⅲ observations[J].Plateau Meteorology,2018,37(5):1143-1160. (in Chinese)

    [15] 錢正安,宋敏紅,李萬源,等.全球、中蒙干旱區(qū)及其部分地區(qū)降水分布細(xì)節(jié)[J].高原氣象,2011,30(1):1-12.

    Qian Z A,Song M H,Li W Y,et al.Global,China-Mongolia Arid-and Semiarid-Areas(CMASA) and the details of precipitation distribution in part regions of CMASA[J]. Plateau Meteorology,2011,30(1):1-12.(in Chinese)

    [16] 張強,張良,崔顯成,等.干旱監(jiān)測與評價技術(shù)的發(fā)展及其科學(xué)挑戰(zhàn)[J].地球科學(xué)進(jìn)展,2011,26(7):763-778.

    Zhang Q,Zhang L,Cui X C,et al.Progresses and challenges in drought assessment and monitoring[J].Advances in Earth Science,2011,26(7):763-778.(in Chinese)

    [17] 張強,韓蘭英,張立陽,等.論氣候變暖背景下干旱和干旱災(zāi)害風(fēng)險特征與管理策略[J].地球科學(xué)進(jìn)展,2014,29(1): 80-91.

    Zhang Q,Han L Y,Zhang L Y,et al.Analysis on the character and management strategy of drought disaster and risk under the climatic warming[J].Advances in Earth Science,2014,29 (1):80-91.(in Chinese)

    [18] 鄒旭愷,張強.近半個世紀(jì)我國干旱變化的初步研究[J].應(yīng)用氣象學(xué)報,2008,19(6):679-686.

    Zou X K,Zhang Q.Preliminary studies on variations in droughts over China during past 50 years[J].Journal of Applied Meteorological Science,2008,19(6):679-686.(in Chinese)

    [19] 王飛,朱仲元,郝祥云,等.雪深再分析資料在錫林浩特地區(qū)的適用性評價[J].中國農(nóng)業(yè)氣象,2019,40(2):76-85.

    Wang F,Zhu Z Y,Hao X Y,et al.Applicability evaluation of snow depth reanalysis data in Xilinhot and its surrounding areas[J].Chinese Journal of Agrometeorology,2019,40(2): 76-85.(in Chinese)

    [20] Pierce D W,Barnett T P,Santer B D,et al.Selecting global climate models for regional climate change studies[J]. Proceedings of the National Academy of Sciences,2009, 106(21):8441-8446.

    [21] 羅立輝,張耀南,周劍,等.基于WRF驅(qū)動的CLM模型對青藏高原地區(qū)陸面過程模擬研究[J].冰川凍土,2013,35(3): 553-564.

    Luo L H,Zhang Y N,Zhou J,et al.Simulation and application of the land surface model CLM driven by WRF in the Tibetan Plateau[J].Journal of Glaciology and Geocryology, 2013,35(3):553-564.(in Chinese)

    [22] 李新,黃春林,車濤,等.中國陸面數(shù)據(jù)同化系統(tǒng)研究的進(jìn)展與前瞻[J].自然科學(xué)進(jìn)展,2007,17(2):163-173.

    Li X,Huang C L,Che T,et al.Progress and foresight in the study of Chinese land data assimilation system[J].Natural Science Progress,2007,17(2):163-173.(in Chinese)

    [23] 龔建東.同化技術(shù):數(shù)值天氣預(yù)報突破的關(guān)鍵:以歐洲中期天氣預(yù)報中心同化技術(shù)演進(jìn)為例[J].氣象科技進(jìn)展,2013,3 (3):6-13.

    Gong J D.Data assimilation:a key technology for NWP:technical review of data assimilation in ECMWF[J]. Advances in Meteorological Science and Technology,2013, 3(3):6-13.(in Chinese)

    [24] 保鴻燕,楊凱,李一玲,等. GLDAS土壤濕度在青藏高原地區(qū)的誤差及其原因[A].S8數(shù)值模式產(chǎn)品應(yīng)用與評估[C].西安:中國氣象學(xué)會年會,2016.

    Bao H Y,Yang K,Li Y L,et al.Error of GLDAS soil moistrue product in Tibetan Plateau region and its cause[A].S8 application and evaluation of numerical model products [A].Xi’an:Chinese Meteorology Society,2016.(in Chinese)

    Applicability Evaluation of CLDAS Merged Soil Relative Moisture in Northeast of China and Its Correction Research

    CUI Yuan-yuan1,ZHANG Qiang2, QIN Jun3,JING Wen-qi4

    (1. Hebei Meteorological Observatory, Shijiazhuang 050000, China;2. Laboratory for Climate Studies, Chinese Meteorological Administration, National Climate Center, Beijing 100081;3. School of Environmental Studies, China University of Geosciences(Wuhan), Wuhan 430074;4. Meteorological Observatory of 94758 Troops, Ningde, Fujian 355103)

    Abstract: Using observed soil relative moisture data at 10cm layer, the applicability of CLDAS(CMA Land Data Assimilation System)merged soil relative moisture at 0?20cm layer was tested and verified by calculating the correlation coefficient and bias between the CLDAS product and in-situ soil moisture observation in Northeast of China during 2008?2017. The results showed that CLDAS merged soil relative moisture data could reasonably reflect the actual spatial distribution of soil moisture at a depth of 10cm in Northeast of China, however, it was not very excellent at the extreme values’ performance. CLDAS merged soil relative moisture data was obviously higher than actual site observations in western of Liaoning province and southwestern of Heilongjiang province, while the simulation was lower than the actual soil relative moisture in northeastern of Heilongjiang province. The applicability of CLDAS product in Northeast of China decreases from southwest to northeast and was better in Jilin and Liaoning provinces. Furthermore, further analysis found that the poor applicability was mainly caused by the low correlation between CLDAS and site observation in Heilongjiang province. Finally, the error of CLDAS product was corrected by regression correction method and 7 ten-day sliding average correction method. The result greatly improves the accuracy of CLDAS product, and the revised products could be applied to soil moisture drought monitoring business.

    Key words: Northeast of China; CLDAS merged soil relative moisture product; CLDAS merged product; Regress correction method; 7 ten-day sliding average correction method; Applicability evaluate

    doi:10.3969/j.issn.1000-6362.2019.10.006

    崔園園,張強,覃軍,等.CLDAS融合土壤濕度產(chǎn)品在東北地區(qū)的適用性評估及訂正[J].中國農(nóng)業(yè)氣象,2019,40(10):660-668

    *收稿日期:2019?02?16

    **通訊作者。E-mail:zhq62@cma.gov.cn

    基金項目:國家重點研發(fā)計劃重大自然災(zāi)害監(jiān)測預(yù)警與防范重點專項(2017YFC1502402)

    作者簡介:崔園園(1991?),碩士生,助理工程師,主要從事氣象資料質(zhì)量控制研究。E-mail:doubleyuancui@foxmail.com

    猜你喜歡
    陸面土壤濕度適用性
    強調(diào)簡潔和適用性 MICHI by Rotel X5/X3合并功放
    土壤濕度傳感器在園林綠化灌溉上的應(yīng)用初探
    戴永久: 砥礪前行,永遠(yuǎn)奮斗
    相似廓線法在星載雷達(dá)陸面衰減訂正中的應(yīng)用
    基于51單片機控制花盆土壤濕度
    電子制作(2019年15期)2019-08-27 01:12:12
    環(huán)保技術(shù)在土木工程領(lǐng)域的適用性探討
    小議階級分析理論的適用性
    基于區(qū)域氣候模式RegCM4的長江流域降水預(yù)報效果分析
    新一代車上轉(zhuǎn)換裝置在廠礦企業(yè)鐵路運輸中的適用性研究
    四川盆地土壤濕度時空分布及影響因子分析
    精品不卡国产一区二区三区| 欧美日本中文国产一区发布| 久久亚洲精品不卡| 午夜精品国产一区二区电影| 亚洲专区中文字幕在线| 97人妻天天添夜夜摸| 一级黄色大片毛片| 欧美日韩精品网址| 亚洲av五月六月丁香网| 成人欧美大片| av天堂久久9| 欧美大码av| 午夜久久久在线观看| 日本一区二区免费在线视频| 最近最新中文字幕大全免费视频| 精品久久久久久久毛片微露脸| 久久精品人人爽人人爽视色| 久久久水蜜桃国产精品网| 成人亚洲精品一区在线观看| 18禁国产床啪视频网站| 一级a爱片免费观看的视频| 韩国精品一区二区三区| 18禁观看日本| 18禁观看日本| 久久人人精品亚洲av| 精品久久久久久,| 亚洲成人精品中文字幕电影| 一级a爱片免费观看的视频| 中文字幕人成人乱码亚洲影| 久久狼人影院| av在线播放免费不卡| 日韩 欧美 亚洲 中文字幕| 一进一出好大好爽视频| 免费在线观看完整版高清| www.999成人在线观看| 亚洲无线在线观看| 亚洲av日韩精品久久久久久密| 国产精品二区激情视频| 宅男免费午夜| 一进一出好大好爽视频| 日韩欧美一区二区三区在线观看| 亚洲成人国产一区在线观看| 亚洲男人天堂网一区| 久久久久久人人人人人| 午夜激情av网站| 国产精品亚洲美女久久久| 久久精品亚洲精品国产色婷小说| 窝窝影院91人妻| 久久青草综合色| 搞女人的毛片| 国产乱人伦免费视频| 激情在线观看视频在线高清| 黄网站色视频无遮挡免费观看| 妹子高潮喷水视频| 亚洲国产高清在线一区二区三 | 中文字幕色久视频| 国产成人欧美在线观看| 免费不卡黄色视频| 成年女人毛片免费观看观看9| 精品日产1卡2卡| 级片在线观看| 9色porny在线观看| 亚洲精品av麻豆狂野| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩精品中文字幕看吧| 看片在线看免费视频| 一区二区三区国产精品乱码| 19禁男女啪啪无遮挡网站| 国产男靠女视频免费网站| 成在线人永久免费视频| 久久香蕉国产精品| 日本a在线网址| 麻豆久久精品国产亚洲av| 日韩免费av在线播放| 中文字幕人妻丝袜一区二区| 中国美女看黄片| 亚洲avbb在线观看| 黄色毛片三级朝国网站| 少妇熟女aⅴ在线视频| 一级黄色大片毛片| 天天躁夜夜躁狠狠躁躁| 亚洲成人精品中文字幕电影| 欧美日本视频| 18禁观看日本| 午夜日韩欧美国产| 一区二区三区国产精品乱码| 女人爽到高潮嗷嗷叫在线视频| 99国产精品一区二区三区| 国产精品一区二区精品视频观看| 每晚都被弄得嗷嗷叫到高潮| 女人被躁到高潮嗷嗷叫费观| 好男人电影高清在线观看| 一区福利在线观看| 一二三四社区在线视频社区8| 免费高清在线观看日韩| 少妇裸体淫交视频免费看高清 | 在线观看免费午夜福利视频| aaaaa片日本免费| 色av中文字幕| 国产精品电影一区二区三区| 好看av亚洲va欧美ⅴa在| 一级黄色大片毛片| 18禁观看日本| 天天躁狠狠躁夜夜躁狠狠躁| 无遮挡黄片免费观看| 男女下面进入的视频免费午夜 | 亚洲一区二区三区色噜噜| 欧美av亚洲av综合av国产av| 亚洲成人国产一区在线观看| 亚洲人成77777在线视频| 69av精品久久久久久| 欧美成人一区二区免费高清观看 | 高清毛片免费观看视频网站| 波多野结衣一区麻豆| 亚洲人成电影免费在线| 亚洲av电影在线进入| 老汉色av国产亚洲站长工具| 日韩三级视频一区二区三区| 欧美激情高清一区二区三区| 亚洲精华国产精华精| 一进一出抽搐动态| 女人被狂操c到高潮| 欧美另类亚洲清纯唯美| 国产人伦9x9x在线观看| 亚洲精品av麻豆狂野| 成人永久免费在线观看视频| 少妇 在线观看| 亚洲精品在线美女| 一区二区三区高清视频在线| 免费看a级黄色片| 久久久水蜜桃国产精品网| 高清毛片免费观看视频网站| av片东京热男人的天堂| 成人18禁在线播放| 国产精品久久电影中文字幕| 亚洲欧美日韩高清在线视频| 国产熟女xx| 黄网站色视频无遮挡免费观看| 涩涩av久久男人的天堂| av视频免费观看在线观看| 亚洲久久久国产精品| 免费高清在线观看日韩| 亚洲人成电影观看| 亚洲成a人片在线一区二区| 欧洲精品卡2卡3卡4卡5卡区| 色尼玛亚洲综合影院| 成人欧美大片| 两性午夜刺激爽爽歪歪视频在线观看 | 波多野结衣av一区二区av| 两个人免费观看高清视频| 视频区欧美日本亚洲| videosex国产| 亚洲精华国产精华精| 国产麻豆成人av免费视频| 国内精品久久久久久久电影| 日本三级黄在线观看| 免费在线观看日本一区| 欧美国产精品va在线观看不卡| 亚洲欧美激情在线| 日韩欧美国产一区二区入口| 亚洲片人在线观看| 久久香蕉国产精品| 欧美日韩乱码在线| 亚洲色图 男人天堂 中文字幕| 桃色一区二区三区在线观看| www.精华液| 啦啦啦韩国在线观看视频| 一本大道久久a久久精品| 99久久久亚洲精品蜜臀av| 久久久国产欧美日韩av| 日韩欧美三级三区| 国产亚洲欧美在线一区二区| 国产精品 国内视频| 亚洲人成电影免费在线| 成人亚洲精品av一区二区| 午夜免费鲁丝| 亚洲午夜精品一区,二区,三区| 亚洲电影在线观看av| 国产成人系列免费观看| 女生性感内裤真人,穿戴方法视频| 午夜福利,免费看| 视频区欧美日本亚洲| 欧美成人性av电影在线观看| 久久香蕉激情| 国产激情久久老熟女| 亚洲av片天天在线观看| 亚洲 欧美 日韩 在线 免费| 此物有八面人人有两片| 少妇被粗大的猛进出69影院| 久久人妻福利社区极品人妻图片| 午夜免费激情av| 亚洲av电影不卡..在线观看| 精品国产国语对白av| 亚洲专区中文字幕在线| 18禁国产床啪视频网站| 国产高清有码在线观看视频 | 精品高清国产在线一区| 激情视频va一区二区三区| 日韩 欧美 亚洲 中文字幕| 美女高潮喷水抽搐中文字幕| 国产高清有码在线观看视频 | 亚洲欧美精品综合一区二区三区| 在线视频色国产色| 长腿黑丝高跟| 青草久久国产| 天堂影院成人在线观看| 搞女人的毛片| 亚洲 国产 在线| 1024香蕉在线观看| av视频在线观看入口| 在线av久久热| 一级毛片精品| 在线观看www视频免费| 男女下面进入的视频免费午夜 | 日本五十路高清| АⅤ资源中文在线天堂| 午夜福利视频1000在线观看 | 最好的美女福利视频网| 在线天堂中文资源库| 又大又爽又粗| 色av中文字幕| 不卡av一区二区三区| av中文乱码字幕在线| 91精品三级在线观看| 亚洲五月天丁香| 亚洲第一欧美日韩一区二区三区| 亚洲中文日韩欧美视频| 亚洲 欧美 日韩 在线 免费| 日本 av在线| 两个人视频免费观看高清| 很黄的视频免费| 亚洲 欧美一区二区三区| 午夜两性在线视频| 久久精品aⅴ一区二区三区四区| 一级片免费观看大全| 成人亚洲精品av一区二区| 精品久久久精品久久久| ponron亚洲| 黄色视频不卡| 欧美日韩福利视频一区二区| 免费在线观看日本一区| 日日夜夜操网爽| 精品不卡国产一区二区三区| 男女下面进入的视频免费午夜 | 日韩精品免费视频一区二区三区| 嫩草影院精品99| 黄色片一级片一级黄色片| 99精品在免费线老司机午夜| 国产一区二区三区视频了| 欧美在线黄色| 午夜亚洲福利在线播放| 精品国产一区二区久久| 精品国内亚洲2022精品成人| 亚洲电影在线观看av| 久久精品国产清高在天天线| 国产精品 国内视频| 成人亚洲精品av一区二区| 一本大道久久a久久精品| 国产亚洲精品第一综合不卡| 国产精品爽爽va在线观看网站 | 99国产精品免费福利视频| av超薄肉色丝袜交足视频| e午夜精品久久久久久久| 真人做人爱边吃奶动态| 精品电影一区二区在线| 黄色视频不卡| 一进一出抽搐动态| 怎么达到女性高潮| 日韩中文字幕欧美一区二区| 男女做爰动态图高潮gif福利片 | 婷婷丁香在线五月| 欧美性长视频在线观看| 男人舔女人的私密视频| 日韩三级视频一区二区三区| 久9热在线精品视频| 国产高清videossex| 亚洲av五月六月丁香网| 日本欧美视频一区| 黑丝袜美女国产一区| 亚洲人成伊人成综合网2020| 看黄色毛片网站| 亚洲专区中文字幕在线| 国产一区二区三区综合在线观看| 99riav亚洲国产免费| 人妻丰满熟妇av一区二区三区| 亚洲欧美激情在线| 岛国视频午夜一区免费看| 久久午夜亚洲精品久久| 日韩欧美在线二视频| 欧美激情久久久久久爽电影 | 精品欧美一区二区三区在线| 757午夜福利合集在线观看| 91在线观看av| 日本黄色视频三级网站网址| 国产私拍福利视频在线观看| 18美女黄网站色大片免费观看| www.精华液| 美女国产高潮福利片在线看| 国产精品久久久av美女十八| 免费高清在线观看日韩| 性欧美人与动物交配| 老司机深夜福利视频在线观看| 亚洲伊人色综图| 色综合站精品国产| 国产国语露脸激情在线看| 精品高清国产在线一区| 乱人伦中国视频| 国产伦人伦偷精品视频| 亚洲一区高清亚洲精品| 国产成+人综合+亚洲专区| 国产在线观看jvid| 一夜夜www| 成人亚洲精品av一区二区| 亚洲欧美精品综合一区二区三区| 少妇被粗大的猛进出69影院| a在线观看视频网站| 九色国产91popny在线| 亚洲欧美日韩另类电影网站| 1024香蕉在线观看| 久久久久久久久久久久大奶| 多毛熟女@视频| 午夜两性在线视频| 久久久水蜜桃国产精品网| 国产极品粉嫩免费观看在线| 99国产精品99久久久久| 自线自在国产av| 免费高清在线观看日韩| 中文字幕色久视频| 成人国产综合亚洲| 美女免费视频网站| 国产av精品麻豆| 性欧美人与动物交配| 午夜激情av网站| 国产精品,欧美在线| 丝袜人妻中文字幕| 国产精品一区二区免费欧美| 不卡一级毛片| 国产欧美日韩精品亚洲av| 操出白浆在线播放| 91精品国产国语对白视频| av中文乱码字幕在线| 国内毛片毛片毛片毛片毛片| 国产97色在线日韩免费| 欧美国产精品va在线观看不卡| 亚洲片人在线观看| 99riav亚洲国产免费| 9色porny在线观看| 黄色 视频免费看| 午夜福利免费观看在线| 日本撒尿小便嘘嘘汇集6| 搡老熟女国产l中国老女人| 国产真人三级小视频在线观看| 亚洲一码二码三码区别大吗| 丰满人妻熟妇乱又伦精品不卡| 999久久久精品免费观看国产| 色综合站精品国产| 狂野欧美激情性xxxx| 久久香蕉精品热| 国产午夜精品久久久久久| 国产区一区二久久| 熟妇人妻久久中文字幕3abv| 丁香欧美五月| 亚洲国产欧美一区二区综合| 脱女人内裤的视频| 日韩欧美国产一区二区入口| 久久中文字幕一级| 中文字幕色久视频| 老熟妇乱子伦视频在线观看| 久久精品国产亚洲av香蕉五月| 真人一进一出gif抽搐免费| 天堂影院成人在线观看| 中文亚洲av片在线观看爽| 在线观看免费视频网站a站| 亚洲成av人片免费观看| 久久青草综合色| 亚洲国产毛片av蜜桃av| 满18在线观看网站| 亚洲天堂国产精品一区在线| 男人舔女人的私密视频| 国产成人精品无人区| 制服人妻中文乱码| 国产精品久久久久久亚洲av鲁大| 精品国产乱码久久久久久男人| 欧美国产日韩亚洲一区| 麻豆久久精品国产亚洲av| 国产精品自产拍在线观看55亚洲| 亚洲免费av在线视频| 老熟妇乱子伦视频在线观看| 国产精品,欧美在线| 国产麻豆69| 欧美国产精品va在线观看不卡| 国产主播在线观看一区二区| 黄片大片在线免费观看| 精品卡一卡二卡四卡免费| 国产精品久久久久久精品电影 | 高清毛片免费观看视频网站| 怎么达到女性高潮| 妹子高潮喷水视频| 久久久久久大精品| 精品日产1卡2卡| 久久久久久人人人人人| 搡老妇女老女人老熟妇| 成人18禁高潮啪啪吃奶动态图| 91字幕亚洲| 亚洲黑人精品在线| 黄色a级毛片大全视频| 涩涩av久久男人的天堂| 亚洲五月天丁香| 黄色毛片三级朝国网站| 久久婷婷人人爽人人干人人爱 | 国产亚洲精品久久久久久毛片| 成在线人永久免费视频| 正在播放国产对白刺激| 女人高潮潮喷娇喘18禁视频| a在线观看视频网站| 黑人操中国人逼视频| 欧美黑人精品巨大| 中文字幕人妻熟女乱码| 亚洲av五月六月丁香网| 久久香蕉国产精品| svipshipincom国产片| 嫩草影视91久久| 日本 欧美在线| 亚洲av成人一区二区三| 国产私拍福利视频在线观看| 国内精品久久久久久久电影| 一级a爱视频在线免费观看| 亚洲人成77777在线视频| 国产精品久久电影中文字幕| 久久久久久人人人人人| av超薄肉色丝袜交足视频| 国产精品野战在线观看| 一区二区日韩欧美中文字幕| 国产麻豆成人av免费视频| 首页视频小说图片口味搜索| 一级a爱视频在线免费观看| av在线天堂中文字幕| 精品久久久久久久久久免费视频| 精品久久久久久,| av网站免费在线观看视频| 一区二区三区国产精品乱码| 亚洲成人国产一区在线观看| 亚洲免费av在线视频| 亚洲中文字幕一区二区三区有码在线看 | 久久久久亚洲av毛片大全| 一边摸一边抽搐一进一小说| 女人被狂操c到高潮| 非洲黑人性xxxx精品又粗又长| 国产精品1区2区在线观看.| 国产精品av久久久久免费| 最近最新中文字幕大全电影3 | 99精品在免费线老司机午夜| 午夜免费观看网址| 一卡2卡三卡四卡精品乱码亚洲| 日韩精品青青久久久久久| 久久久久九九精品影院| 在线十欧美十亚洲十日本专区| 欧美日韩中文字幕国产精品一区二区三区 | 大码成人一级视频| 大型黄色视频在线免费观看| 在线观看午夜福利视频| 无限看片的www在线观看| 国产一区二区在线av高清观看| 国产高清videossex| 欧美在线黄色| 无人区码免费观看不卡| 亚洲熟女毛片儿| 欧美成人性av电影在线观看| 精品国产国语对白av| 欧美日韩中文字幕国产精品一区二区三区 | 国产片内射在线| 国产成人系列免费观看| 18禁黄网站禁片午夜丰满| 12—13女人毛片做爰片一| 黄片大片在线免费观看| 国产伦一二天堂av在线观看| 国产一区二区激情短视频| 50天的宝宝边吃奶边哭怎么回事| av视频免费观看在线观看| 老司机午夜福利在线观看视频| 日本免费一区二区三区高清不卡 | 国产精品99久久99久久久不卡| 国产精品av久久久久免费| 成人欧美大片| 亚洲,欧美精品.| 国内久久婷婷六月综合欲色啪| 美女大奶头视频| 国产亚洲欧美98| 神马国产精品三级电影在线观看 | 美女免费视频网站| 国产麻豆69| 天天躁夜夜躁狠狠躁躁| tocl精华| 动漫黄色视频在线观看| 十八禁人妻一区二区| 欧美精品亚洲一区二区| 久久中文字幕人妻熟女| 亚洲 欧美一区二区三区| 午夜福利,免费看| 欧洲精品卡2卡3卡4卡5卡区| 成人三级黄色视频| 日韩一卡2卡3卡4卡2021年| 日韩免费av在线播放| 欧美人与性动交α欧美精品济南到| 亚洲 欧美 日韩 在线 免费| 一级,二级,三级黄色视频| 超碰成人久久| 91成年电影在线观看| 中文字幕av电影在线播放| 日本精品一区二区三区蜜桃| 国产片内射在线| 久久精品aⅴ一区二区三区四区| 老司机午夜福利在线观看视频| 久久久久亚洲av毛片大全| 欧美人与性动交α欧美精品济南到| 久久香蕉精品热| 国产成人影院久久av| 亚洲午夜精品一区,二区,三区| а√天堂www在线а√下载| 一本大道久久a久久精品| 99国产综合亚洲精品| 9191精品国产免费久久| 一个人免费在线观看的高清视频| 免费观看精品视频网站| 亚洲熟妇中文字幕五十中出| 国产aⅴ精品一区二区三区波| 日韩欧美一区二区三区在线观看| 美女扒开内裤让男人捅视频| 色在线成人网| 桃色一区二区三区在线观看| 国产欧美日韩一区二区三区在线| 国产色视频综合| 亚洲成人免费电影在线观看| 在线观看免费日韩欧美大片| 亚洲电影在线观看av| 给我免费播放毛片高清在线观看| 成人特级黄色片久久久久久久| 亚洲精品久久国产高清桃花| 黑人巨大精品欧美一区二区mp4| 最近最新免费中文字幕在线| 久久精品国产综合久久久| 亚洲,欧美精品.| 十分钟在线观看高清视频www| 亚洲国产高清在线一区二区三 | 午夜精品在线福利| 国产日韩一区二区三区精品不卡| 亚洲国产高清在线一区二区三 | 一本综合久久免费| 免费高清在线观看日韩| 一级毛片精品| 国产精品,欧美在线| 国产精品一区二区在线不卡| 国产主播在线观看一区二区| 法律面前人人平等表现在哪些方面| 一区福利在线观看| 高清在线国产一区| 1024香蕉在线观看| 国产精品一区二区免费欧美| 国产精品一区二区三区四区久久 | 一区二区三区精品91| 亚洲三区欧美一区| 桃红色精品国产亚洲av| 黑人巨大精品欧美一区二区mp4| 这个男人来自地球电影免费观看| 免费无遮挡裸体视频| 99国产综合亚洲精品| 91九色精品人成在线观看| 精品国产乱码久久久久久男人| 国产一区二区在线av高清观看| 制服人妻中文乱码| АⅤ资源中文在线天堂| 中文字幕av电影在线播放| 久久精品人人爽人人爽视色| 亚洲国产欧美网| 国产精品自产拍在线观看55亚洲| 精品久久久精品久久久| 久久人妻福利社区极品人妻图片| 搞女人的毛片| 99国产精品一区二区蜜桃av| 久热爱精品视频在线9| 高清毛片免费观看视频网站| 妹子高潮喷水视频| 99香蕉大伊视频| 亚洲精品一区av在线观看| 女人被狂操c到高潮| 国产又色又爽无遮挡免费看| 国产精品永久免费网站| 他把我摸到了高潮在线观看| 制服人妻中文乱码| 亚洲精品美女久久久久99蜜臀| 怎么达到女性高潮| 国产成人精品在线电影| 国产亚洲精品一区二区www| 最近最新中文字幕大全免费视频| 国产精品98久久久久久宅男小说| 久久性视频一级片| 在线观看66精品国产| 免费看十八禁软件| 国产午夜福利久久久久久| 欧美人与性动交α欧美精品济南到| 97碰自拍视频| 女人精品久久久久毛片| 久久狼人影院| 午夜久久久久精精品| 变态另类丝袜制服| 色播在线永久视频| 亚洲av第一区精品v没综合| 亚洲专区国产一区二区| 久久午夜亚洲精品久久| 纯流量卡能插随身wifi吗| 欧美性长视频在线观看| 国产欧美日韩一区二区三区在线| 亚洲欧美激情综合另类| 色在线成人网| 久久伊人香网站|