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

    淮河流域表層土壤濕度時(shí)空特征及其與地面降水的關(guān)系

    2017-06-15 15:28:56黃勇邱旭敏黃國(guó)貴
    關(guān)鍵詞:淮河流域土壤濕度表層

    黃勇,邱旭敏,黃國(guó)貴

    1. 安徽省氣象科學(xué)研究所 安徽省大氣科學(xué)與衛(wèi)星遙感重點(diǎn)實(shí)驗(yàn)室,安徽 合肥 230031;2. 壽縣國(guó)家氣候觀象臺(tái),安徽 壽縣 232200;3. 天津市濱海新區(qū)氣象局,天津 300457

    淮河流域表層土壤濕度時(shí)空特征及其與地面降水的關(guān)系

    黃勇1,2,邱旭敏3,黃國(guó)貴2

    1. 安徽省氣象科學(xué)研究所 安徽省大氣科學(xué)與衛(wèi)星遙感重點(diǎn)實(shí)驗(yàn)室,安徽 合肥 230031;2. 壽縣國(guó)家氣候觀象臺(tái),安徽 壽縣 232200;3. 天津市濱海新區(qū)氣象局,天津 300457

    土壤濕度是陸面與大氣能量和物質(zhì)交換的關(guān)鍵物理量,除地面觀測(cè)外,衛(wèi)星遙感也是獲取土壤濕度的重要手段。本文針對(duì)歐洲衛(wèi)星氣象中心(EUMETSAT)基于衛(wèi)星遙感而研發(fā)的面向水文氣象衛(wèi)星應(yīng)用產(chǎn)品(H-SAF)中的土壤水分指數(shù)(SM-DAS-2)產(chǎn)品,分析產(chǎn)品在淮河流域土壤水分監(jiān)測(cè)中的適用性,為在淮河流域利用該產(chǎn)品開(kāi)展農(nóng)田生態(tài)監(jiān)測(cè)奠定基礎(chǔ)。首先,分析了該產(chǎn)品監(jiān)測(cè)得到的中國(guó)土壤表層水分空間分布特征,在此基礎(chǔ)上,總結(jié)淮河流域表層土壤濕度的時(shí)空分布特征;并結(jié)合中國(guó)氣象局發(fā)布的中國(guó)逐日降水量格點(diǎn)化產(chǎn)品,分析表層土壤濕度與降水量之間的相關(guān)關(guān)系,最后通過(guò)分析淮河流域土壤濕度與地面降水之間的對(duì)應(yīng)關(guān)系,探討土壤濕度對(duì)降水的響應(yīng)時(shí)間。從分析結(jié)果來(lái)看,對(duì)于全國(guó)而言,該產(chǎn)品能夠有效地反映出土壤濕度的分布特征,并且能夠刻畫(huà)出較細(xì)致的空間分布規(guī)律。在淮河流域,土壤濕度產(chǎn)品和人工觀測(cè)數(shù)值之間具有性較好的相關(guān)性(相關(guān)系數(shù)為0.51,通過(guò)置信度0.05的t檢驗(yàn)),但SM-DAS-2產(chǎn)品指數(shù)偏低(SMDAS-2產(chǎn)品的年均值為0.6,地面觀測(cè)站監(jiān)測(cè)的表層土壤濕度為0.69)。表層土壤濕度總體上保持“南高北低”空間分布特征,并且空間分布存在著季節(jié)性差異。山區(qū)和河流附近區(qū)域是土壤濕度的高值中心,這一特征在夏秋兩季最為顯著。表層土壤濕度與降水之間存在顯著的相關(guān)關(guān)系,而在山區(qū)和河流附近,兩者間的相關(guān)性變差。從時(shí)間上看,夏季是淮河流域表層土壤濕度與同期降水相關(guān)性最好的季節(jié),總體上降水發(fā)生3 d以后表層土壤濕度才會(huì)對(duì)降水的變化產(chǎn)生響應(yīng),土壤濕度對(duì)降水變化的響應(yīng)時(shí)間最長(zhǎng)能延長(zhǎng)到10~15 d。

    淮河流域;SM-DAS-2;土壤濕度;降水

    土壤濕度作為陸面與大氣能量和水分交換研究中的關(guān)鍵物理量,一直備受重視。Delworth et al.(1988)研究表明,在大陸上,土壤濕度和雪蓋對(duì)于氣候的“記憶效應(yīng)”的作用幾乎同等重要。Dirmeyer(2000)通過(guò)GCM模式模擬,認(rèn)為在干濕氣候的過(guò)渡地區(qū)土壤濕度對(duì)提高降水的可預(yù)報(bào)性具有較大意義。符淙斌等(1995)也指出氣候過(guò)渡帶在環(huán)境變化中比較敏感,氣候變率更加顯著。孫丞虎等(2005)通過(guò)分析中國(guó)東部1990—2000年旬平均土壤濕度的時(shí)空分布特征指出,表層土壤濕度異常的第一旋轉(zhuǎn)空間模態(tài)出現(xiàn)在淮河流域,正是淮河流域地表層對(duì)氣候變率最為敏感的反映。周娟等(2017)研究表明,土壤濕度與高原季風(fēng)之間也存在顯著相關(guān)關(guān)系。李潤(rùn)春等(2017)認(rèn)為,渭河流域前一年秋季和當(dāng)年春季土壤濕度與東亞夏季風(fēng)強(qiáng)度指數(shù)呈顯著正相關(guān),秋季和春季陸面土壤濕度狀況對(duì)東亞夏季風(fēng)的強(qiáng)弱有正反饋的作用。此外,青藏高原的土壤濕度對(duì)中國(guó)東部地區(qū)的夏季降水有著顯著影響(李登宣等,2016;王靜等,2016)。以上研究均表明,土壤濕度除對(duì)農(nóng)作物存在顯著影響(董彬,2012;沈振西等,2015;沈振西等,2016;付剛等,2017;)外,還對(duì)于天氣氣候有著重要的作用。

    歐洲中期天氣預(yù)報(bào)中心(ECMWF)和美國(guó)環(huán)境預(yù)報(bào)中心(NCEP)的再分析資料,全球土壤濕度計(jì)劃的資料(GSWP2),以及國(guó)家大氣研究中心(NCAR)利用觀測(cè)資料強(qiáng)迫“通用陸面模式”CLM所產(chǎn)生的土壤濕度資料,是幾種較為常用的土壤濕度資料。國(guó)內(nèi)外學(xué)者也對(duì)4種土壤濕度資料的應(yīng)用做了大量分析研究。張文君等(2008)對(duì)4種資料在中國(guó)的可靠性做了詳細(xì)分析。崔文瑞等(2009)對(duì)ECMWF的ERA40和NCEP的R-1中黑河流域的再分析資料進(jìn)行對(duì)比分析。胡婭敏等(2009)將GSWP2資料同化進(jìn)區(qū)域氣候模式BCC_RegCM1.0(簡(jiǎn)稱(chēng)RegCM)的土壤濕度場(chǎng)中,改進(jìn)了區(qū)域氣候模擬。郭維棟等(2007)通過(guò)將一個(gè)土壤濕度反演模型引入IAP跨季度氣候預(yù)測(cè)系統(tǒng)中,探討了土壤濕度應(yīng)用于季節(jié)降水預(yù)測(cè)的可行性及效果。王迪等(2015)對(duì)CLM4.5數(shù)據(jù)集的分析表明,CLM4.5存在高估表層土壤濕度的現(xiàn)象。

    隨著衛(wèi)星遙感技術(shù)的發(fā)展,越來(lái)越多的主被動(dòng)微波傳感器資料在土壤濕度估計(jì)中得到了深入的應(yīng)用。包括:2009年歐洲太空局(ESA)發(fā)射的土壤濕度與海水鹽度衛(wèi)星SMOS;2015年美國(guó)國(guó)家航天局(NASA)發(fā)射的SMAP(Soil Moisture Active and Passive)衛(wèi)星;以及Metop衛(wèi)星上搭載的ASCAT(Advanced Scatterometer)和Aqua衛(wèi)星上搭載的AMSR-E(Advanced Microwave Scanning Radiometer-Earth Observing System)(張秀英等,2010;蘭鑫宇等,2015;李哲等,2017)。這些衛(wèi)星和探測(cè)器都為大范圍的土壤濕度數(shù)據(jù)收集提供了基礎(chǔ)的遙感觀測(cè)數(shù)據(jù)。

    根區(qū)土壤水分指數(shù)(SM-DAS-2)產(chǎn)品,是歐洲衛(wèi)星氣象中心(EUMETSAT)組織研制的一套土壤濕度產(chǎn)品。該產(chǎn)品是面向水文的氣象衛(wèi)星應(yīng)用產(chǎn)品(H-SAF)之一,通過(guò)同化ASCAT反演的表層土壤濕度產(chǎn)品后得到。作為一種土壤濕度資料,在歐洲經(jīng)過(guò)地面觀測(cè)數(shù)據(jù)對(duì)比分析得出,同化后的土壤水分產(chǎn)品比ASCAT數(shù)據(jù)直接反演的土壤水分產(chǎn)品具有較小的偏差和均方根誤差。反演得到的地表土壤水分(表層2 cm)被向下擴(kuò)展到地表下2.89 m,提供了4層土壤水分產(chǎn)品,分別是0~7、7~28、28~100、100~289 cm。

    本文主要應(yīng)用SM-DAS-2產(chǎn)品分析土壤濕度的時(shí)空特征及其與降水的相關(guān)性,以更好地應(yīng)用這種土壤濕度數(shù)據(jù)。考慮到淮河流域?yàn)楦蓾駳夂蜻^(guò)渡帶,其表層對(duì)氣候變率最為敏感,且SMDAS-2產(chǎn)品中表層以下均由模式模擬得到,因此主要探討淮河流域表層土壤(0~7 cm)與降水的相關(guān)性。

    1 資料和方法

    所用資料包括土壤水分指數(shù)(SM-DAS-2)產(chǎn)品的表層(0~7 cm)土壤單位體積含水量,以及空間分辨率。該降水產(chǎn)品基于中國(guó)2419國(guó)家級(jí)地面氣象站日降水量觀測(cè)數(shù)據(jù),采用“基于氣候背景場(chǎng)”的最優(yōu)插值方法,實(shí)時(shí)生成格點(diǎn)化產(chǎn)品(沈艷等,2010)。

    在分析土壤濕度與降水間相關(guān)關(guān)系時(shí),分別計(jì)算時(shí)間相關(guān)系數(shù)和空間相關(guān)系數(shù),并計(jì)算時(shí)滯相關(guān)系數(shù)。由于SM-DAS-2產(chǎn)品的空間分辨率0.225°×0.225°與日降水產(chǎn)品0.25°×0.25°的空間分辨率不相匹配,在計(jì)算空間相關(guān)系數(shù)時(shí),以0.25°×0.25°降水網(wǎng)格為基準(zhǔn),采用間接重采樣的方法,將SM-DAS-2產(chǎn)品處理成0.25°×0.25°的產(chǎn)品,使之與降水產(chǎn)品的空間分辨率相一致,在此基礎(chǔ)上計(jì)算兩者的時(shí)空相關(guān)系數(shù)。

    空間相關(guān)系數(shù)計(jì)算公式:

    式中,A表示土壤濕度,B表示降水量;A 表示土壤濕度的空間平均值,B表示降水量的空間平均值;i,j為區(qū)域內(nèi)網(wǎng)格點(diǎn)的橫、縱坐標(biāo),N、M為橫、縱向網(wǎng)格點(diǎn)的總數(shù)。

    時(shí)間相關(guān)系數(shù)計(jì)算公式:

    式中,A表示土壤濕度,B表示降水量;A 表示土壤濕度的時(shí)間平均值,B表示降水量的時(shí)間平均值;i表示一年中的第i天。

    2 土壤濕度分布特征

    圖1 中國(guó)區(qū)域土壤濕度指數(shù)年平均空間分布圖Fig. 1 Spatial distrbution of soil moisture index for China and nearby

    圖1 為2012年SM-DAS-2土壤水分指數(shù)產(chǎn)品中表層土壤濕度指數(shù)全年平均值的空間分布。從圖中可以看出,表層相對(duì)濕度總體上呈現(xiàn)南高北低的空間分布,長(zhǎng)江以南的區(qū)域?yàn)橥寥栏邼穸葏^(qū)域,最大值中心區(qū)域位于中國(guó)西南中部(四川省一帶);濕度最小的區(qū)域位于西北地區(qū)新疆、甘肅到內(nèi)蒙西部一帶和青藏高原西北部。另外,中國(guó)東北地區(qū)土壤濕度也相對(duì)較高,遼寧—吉林一帶處于土壤濕度相對(duì)較高的中心。孫丞虎等(2005)利用1990—2000年?duì)顩r較好的中國(guó)114站的土壤濕度(土壤重量含水百分率)資料,得到中國(guó)土壤濕度空間分布的氣候特征:0~50 cm各層土壤濕度均呈“兩濕一干”分布,其中,東北、華東較為潮濕,華北則較為干燥,尤其河套地區(qū)是土壤干中心。與空間分布的氣候特征相比,SM-DAS-2土壤水分指數(shù)產(chǎn)品的空間分布特征總體上相似,但存在局部差異。與地面觀測(cè)站得到的氣候特征相比,SM-DAS-2產(chǎn)品的差異主要表現(xiàn)在,(1)干中心位置和濕區(qū)極值中心位置不同;SM-DAS-2產(chǎn)品的干中心主要位于西北地區(qū)和青藏高原西北部,濕區(qū)極值中心位于位于四川省一帶。(2)SM-DAS-2產(chǎn)品刻畫(huà)的空間分布更加細(xì)膩。這主要是因?yàn)榈孛嬲军c(diǎn)的數(shù)量較少(114個(gè))、空間分布不均勻(缺少西北地區(qū)的觀測(cè)站),導(dǎo)致氣候特征的空間分辨率和代表性不足,難以刻畫(huà)細(xì)致的空間分布,以及西北地區(qū)的土壤濕度分布。

    從不同季節(jié)的空間分布(圖2)來(lái)看,空間分布規(guī)律總體上與年均分布相近,呈現(xiàn)出南高北低的特征。從馬思源等(2016)對(duì)歐洲航天局(European Space Agency)和歐洲中期天氣預(yù)報(bào)中心(ECMWF)兩套再分析土壤濕度產(chǎn)品的對(duì)比分析結(jié)果來(lái)看,SM-DAS-2產(chǎn)品在春、夏、秋3個(gè)季節(jié)的土壤空間分布與這兩套再分析資料的結(jié)果相近,也就是SM-DAS-2產(chǎn)品能夠反映出中國(guó)春、夏、秋3個(gè)季節(jié)土壤濕度的分布情況。另外,就冬季而言,在40°N以北的區(qū)域以及青藏高原和西部地區(qū),SM-DAS-2產(chǎn)品所獲取的數(shù)值均低于0.1,土壤濕度小。冬季地面積雪、結(jié)冰、凍土等現(xiàn)象均會(huì)影響微波遙感探測(cè)土壤濕度的能力,而SM-DAS-2產(chǎn)品是以被動(dòng)衛(wèi)星微波遙感探測(cè)為基礎(chǔ)的土壤濕度產(chǎn)品,故該產(chǎn)品會(huì)受到冬季地面積雪、結(jié)冰、凍土等因素影響而在中國(guó)北部、西部以及青藏高原地區(qū)表現(xiàn)出大范圍的低值中心。此外,在SM-DAS-2產(chǎn)品空間分布圖中,還存著一些空白的區(qū)域。由于地表水體對(duì)土壤濕度的微波反演也存在一定影響,在計(jì)算過(guò)程中需要剔除地表水體的影響。因此,SM-DAS-2產(chǎn)品中剔除了內(nèi)陸大面積湖泊的土壤濕度反演數(shù)據(jù),在空間分布圖中形成了一些零星分布的空白區(qū)域。

    圖2 中國(guó)區(qū)域不同季節(jié)土壤濕度指數(shù)的空間分布Fig. 2 Spatial distrbution in difference season of soil moisture index for China and nearby

    圖3 淮河流域全年平均的土壤濕度指數(shù)(左)和日降水量(右,單位:mm)空間分布Fig. 3 Spatial distribution of soil moisture (Left) and precipitation (Right, Unit: mm) for Huaihe River Basin and nearby

    淮河流域作為干濕地區(qū)的過(guò)渡帶,干濕度隨季節(jié)變化明顯。從全年平均的空間分布(圖3左)來(lái)看,淮河流域的土壤濕度總體上呈現(xiàn)隨緯度升高而降低的“南高北低”空間分布規(guī)律。其中,淮河流域中部洪澤湖一帶出現(xiàn)了一個(gè)濕度大值中心。與地面觀測(cè)站監(jiān)測(cè)的表層土壤濕度相比,2012年兩者相關(guān)性較好(相關(guān)系數(shù)為0.51,通過(guò)置信度0.05的t檢驗(yàn)),但SM-DAS-2產(chǎn)品指數(shù)偏低(SM-DAS-2產(chǎn)品的年均值為0.6,地面觀測(cè)站監(jiān)測(cè)的表層土壤濕度為0.69)。

    日降水量分布(圖3右)也呈現(xiàn)“南高北低”空間分布規(guī)律,即空間分布上,土壤濕度和日降水量之間總體上有著較好的一致性:降水量大的區(qū)域,土壤濕度大;反之,降水量小的區(qū)域,土壤濕度小。

    春季,淮河流域土壤濕度具有顯著的緯向分布特征(圖4左),隨著緯度的升高,土壤濕度逐漸降低。與全年平均土壤濕度的分布(圖3左)相比,洪澤湖附近的濕中心不是很顯著。日降水量的分布(圖3右)與土壤濕度保持很好的一致性。

    夏季,淮河流域的土壤濕度呈現(xiàn)出多個(gè)高值中心,且隨緯度升高、土壤濕度降低的特征也不顯著(圖4)。淮河流域及周邊區(qū)域內(nèi),存在多個(gè)主要的濕區(qū)中心,分別位于淮河中部和東北部、大別山區(qū)、皖南山區(qū)、長(zhǎng)江中游、山東丘陵區(qū)以及伏牛山附近。從降水的空間分布(圖4右)來(lái)看,大降水中心主要位于大別山區(qū)、皖南山區(qū)、長(zhǎng)江中下游、淮河流域東部,以及淮河流域北側(cè)局部區(qū)域。從空間分布來(lái)看,除淮河?xùn)|北部、大別山區(qū)和皖南山區(qū)這3個(gè)高濕度的區(qū)域存在相對(duì)應(yīng)的大降水中心外,淮河中部、長(zhǎng)江中游、山東丘陵區(qū)以及伏牛山附近的高濕度區(qū)域,未出現(xiàn)相應(yīng)的大降水中心。也就是說(shuō),在夏季土壤濕度的空間分布與日降水量之間存在的一定的差異。長(zhǎng)江中游、淮河中部均為主河道附近,長(zhǎng)江和淮河周邊區(qū)域土壤相對(duì)較濕,即使在本地?zé)o大降水的情況下,上游的來(lái)水也能使得這些區(qū)域具有較高的土壤濕度。就山東丘陵區(qū)和伏牛山而言,山區(qū)內(nèi)的植被覆蓋降低了土壤水分的蒸發(fā),從而導(dǎo)致區(qū)域內(nèi)的土壤具有相對(duì)較高的濕度。

    秋季,土壤濕度的空間分布與夏季較為相似,只是伏牛山和山東丘陵的高值中心消失,并且在土壤濕度的數(shù)值上總體有所降低。降水的空間分布總體呈現(xiàn)“南多北少”,但是在淮河流域東部出現(xiàn)了2個(gè)降水大值中心。降水和土壤濕度的空間分布特征大體相似,但存在一些差異。主要表現(xiàn)為,(1)大別山區(qū)和大別山區(qū)的高土壤濕度區(qū),均未對(duì)應(yīng)有大降水中心;(2)沿淮區(qū)域出現(xiàn)的高濕度中心也未對(duì)應(yīng)有大降水中心?;春右阅仙絽^(qū)的植被降低了土壤水分的蒸發(fā),故雖然大別山區(qū)和皖南山區(qū)的降水量不是區(qū)域最大,但其土壤濕度卻是區(qū)域最高的?;春又械拇罅克w是造成沿淮區(qū)域土壤濕度偏高的主要因素。

    圖4 淮河流域季節(jié)平均土壤濕度(左)和降水(右,單位:mm)的空間分布Fig. 4 Spatial distribution of soil moisture index (Left) and precipitation (Right, Unit: mm) in difference season in Huaihe River Basin and nearby

    就冬季而言,土壤濕度在總體保持“南高北低”這種分布規(guī)律的基礎(chǔ)上,還在長(zhǎng)江中游、皖南山區(qū)和沿淮出現(xiàn)了3個(gè)高濕中心。受冬季固態(tài)降水的影響,多數(shù)緯度偏北的觀測(cè)都難以有效評(píng)估包含了固態(tài)降水的降水總量,因此日降水產(chǎn)品的空間分布特征僅反映降水的大致分布(維持“南高北低”的特征),難以細(xì)致地刻畫(huà)出區(qū)域內(nèi)降水的情況。因此,也就無(wú)法通過(guò)土壤濕度和降水的空間分布來(lái)評(píng)估土壤濕度和降水的相關(guān)性。

    以上分析表明,SM-DAS-2產(chǎn)品能夠有效地描繪出中國(guó)土壤表層水分含量的空間分布,就淮河流域而言,土壤濕度呈現(xiàn)出了“南高北低”的特征,并且在山區(qū)和河流附近(沿長(zhǎng)江、沿淮河),土壤表層的水分含量較高。從表層土壤濕度與降水量空間分布的對(duì)應(yīng)關(guān)系來(lái)看,總體上兩者具有較好的對(duì)應(yīng)關(guān)系,但這種對(duì)應(yīng)關(guān)系在山區(qū)和河流沿岸存在著不一致性。

    圖5為淮河流域逐日土壤濕度與降水量的時(shí)序圖,從圖中可以看出,總體上土壤濕度對(duì)降水有較好的響應(yīng),降水過(guò)程發(fā)生后,大部分土壤都出現(xiàn)了濕度增加的現(xiàn)象。但是,夏季流域內(nèi)存在著土壤濕度隨降水發(fā)生而持續(xù)降低的現(xiàn)象。如2012年7月初,雖然有間斷性的降水發(fā)生,但是土壤濕度卻出現(xiàn)了持續(xù)的降低。由于7月初淮河流域處于梅雨期,持續(xù)性的降水使得淮河流域的土壤濕度出現(xiàn)了飽和(土壤濕度指數(shù)接近1)。土壤飽和以后,降水過(guò)程無(wú)法再增加土壤濕度,而是形成地表匯流向主河道匯集。這也是造成7月初土壤濕度與降水之間相關(guān)關(guān)系較差的主要原因:土壤中水分飽和以后,降水無(wú)法再補(bǔ)充土壤水分。而在其他季節(jié),由于土壤中的水分沒(méi)有達(dá)到飽和,降水的發(fā)生能夠有效地補(bǔ)充土壤中的水分,使得土壤濕度增加,此時(shí)土壤濕度與降水量之間具有很好的對(duì)應(yīng)關(guān)系。

    3 土壤濕度與降水之間的關(guān)系

    為進(jìn)一步了解淮河流域土壤濕度與降水之間的關(guān)系,選取相關(guān)系數(shù)評(píng)判兩者之間的同步關(guān)系。圖6所示為淮河流域及周邊地區(qū)逐日土壤濕度與降水量之間相關(guān)系數(shù)的空間分布??傮w上,土壤濕度與降水呈顯著正相關(guān),且都通過(guò)了置信度0.05的t檢驗(yàn)。區(qū)域內(nèi)相關(guān)性最好的區(qū)域位于淮河流域中部和西部地區(qū)。

    圖6 全年土壤濕度與降水間相關(guān)系數(shù)的空間分布Fig. 6 Spatial distribution of the correlation coefficient between soil moisture and precipitation

    從4個(gè)季節(jié)相關(guān)系數(shù)的空間分布(圖7)來(lái)看,春季相關(guān)系數(shù)最大的區(qū)域位于長(zhǎng)江中游附近,此外在淮河流域東南側(cè)(長(zhǎng)江下游)、山東丘陵和伏牛山附近分別存在1個(gè)相關(guān)系數(shù)大于0.5的區(qū)域。沿淮流域西北部及西側(cè)山區(qū)土壤濕度和降水之間相關(guān)程度較低,相關(guān)系數(shù)在0.2以下,部分區(qū)域相關(guān)系數(shù)接近0,表示無(wú)相關(guān)性。

    夏季的相關(guān)性最好,淮河流域大部分區(qū)域的相關(guān)系數(shù)都大于0.6,北部地區(qū)相關(guān)系數(shù)超過(guò)0.9,表現(xiàn)出很好的一致性變化。但是,在長(zhǎng)江流域(包括中游和下游),兩者的相關(guān)系數(shù)較低。尤其是在長(zhǎng)江下游,相關(guān)系數(shù)小于0.2,部分區(qū)域接近0。

    秋季是一年四季中土壤濕度和降水相關(guān)性最差的季節(jié),大部分區(qū)域相關(guān)系數(shù)都在0.2~0.4之間,只有極個(gè)別小區(qū)域的相關(guān)系數(shù)超過(guò)0.4。此時(shí),長(zhǎng)江中游土壤濕度和降水的相關(guān)性最大,而大部分區(qū)域相關(guān)系數(shù)都小于0.2。

    圖5 土壤濕度指數(shù)與降水量的逐日變化Fig. 5 Time series of soil moisture index and precipitation

    圖7 不同季節(jié)土壤濕度與降水間相關(guān)系數(shù)的空間分布Fig. 7 Spatial distributions in different seasons of correlation coefficient between soil moisture and precipitation

    冬季土壤濕度和降水相關(guān)系數(shù)呈現(xiàn)“南北低、中間高”的分布特征,相關(guān)性最好的區(qū)域位于研究區(qū)域的中部,其中,淮河流域西側(cè)山區(qū)的相關(guān)系數(shù)最大,相關(guān)系數(shù)最大值超過(guò)0.6。相關(guān)系數(shù)最小的區(qū)域位于淮河流域西北部,最小值接近0,表示無(wú)明顯相關(guān)性。另外,在長(zhǎng)江下游,相關(guān)性也較差,相關(guān)系數(shù)在0.3以下,部分區(qū)域低于0.1。

    圖8 淮河流域土壤濕度和降水量相關(guān)系數(shù)時(shí)序圖Fig. 8 Time series of correlation coefficient between soil moisture and precipitation

    圖9 時(shí)滯相關(guān)系數(shù)(負(fù)X軸為降水超前,正X軸為降水滯后)Fig. 9 Lag correlation coefficient (negative X axis: precipitation ahead, positive X axis: precipitation lag)

    從相關(guān)系數(shù)的逐日變化(圖8)可知(注:流域降水為0時(shí),無(wú)法計(jì)算當(dāng)日相關(guān)系數(shù),因此對(duì)其作滑動(dòng)平均處理),土壤濕度與當(dāng)日降水間的相關(guān)關(guān)系存在著明顯起伏式變化,但總體上存在顯著的相關(guān)性。從時(shí)滯相關(guān)系數(shù)(圖9)可知,當(dāng)降水超前土壤濕度時(shí),可維持10~15 d的顯著正相關(guān);而當(dāng)降水落后土壤濕度時(shí),能維持3~55 d的顯著正相關(guān)。同時(shí),從圖9中可以看出,在降水超前3 d時(shí),相關(guān)系數(shù)最高。也就是說(shuō),降水發(fā)生后,大多數(shù)情況下3 d以后表層土壤濕度才會(huì)對(duì)降水的變化產(chǎn)生響應(yīng)。土壤濕度對(duì)降水變化的響應(yīng)時(shí)間最長(zhǎng)能延長(zhǎng)到10~15 d。

    4 結(jié)論

    利用歐洲衛(wèi)星氣象中心(EUMETSAT組織)的SM-DAS-2土壤水分指數(shù)產(chǎn)品,在對(duì)比分析了土壤濕度空間分布特征與已有研究成果之間一致性的基礎(chǔ)上,討論了淮河流域表層土壤濕度的空間分布特征,及其與降水之間的關(guān)系。結(jié)果表明:

    (1)SM-DAS-2土壤水分指數(shù)產(chǎn)品能夠較好地反映中國(guó)表層土壤濕度的空間分布規(guī)律,并且具有較高的空間分辨率,能夠刻畫(huà)出較細(xì)致的空間分布特征。

    (2)淮河流域的表層土壤濕度總體上呈現(xiàn)“南高北低”的空間分布特征,并且還存在著季節(jié)性變化。山區(qū)和河流附近土壤濕度都較高,并且這一特征在夏秋兩季最為明顯,冬季最不顯著。

    (3)淮河流域降水與表層土壤濕度有顯著的季節(jié)性特征,且二者空間分布特征較為相似;但對(duì)于山區(qū)和河流附近土壤而言,降水與土壤濕度間的關(guān)系存在著不一致性。夏季是淮河流域表層土壤濕度與同期降水相關(guān)性最好的季節(jié),春秋季節(jié)較差。

    (4)對(duì)于淮河流域而言,總體上在降水發(fā)生3 d以后表層土壤濕度才會(huì)對(duì)降水的變化產(chǎn)生響應(yīng)。土壤濕度對(duì)降水變化的響應(yīng)時(shí)間最長(zhǎng)能延長(zhǎng)到10~15 d。

    致謝:感謝歐洲衛(wèi)星氣象中心和國(guó)家氣象信息中心數(shù)據(jù)服務(wù)網(wǎng)提供的土壤濕度和網(wǎng)格化降水產(chǎn)品;感謝審稿專(zhuān)家提供的寶貴意見(jiàn)。

    DELWORTH T L,SYUKUROMANABE. 1988. The influence of potential evaporation on the variabilities of sumilated soil wetness and climate [J]. Journals of Climate, 1(5): 523-547.

    DIRMEYER P A. 2000. Using a global soil wetness dataset to improve seasonal climate simulation [J]. Journals of Climate, 13(16): 2900-2922.

    崔文瑞,高艷紅,彭雯. 2009. 兩套土壤濕度再分析資料在黑河流域的對(duì)比分析[J]. 高原氣象, 28(6): 1274-1281.

    董彬. 2012. 土壤溫度和濕度對(duì)魯東南楊樹(shù)人工林土壤呼吸的協(xié)同影響[J]. 生態(tài)環(huán)境學(xué)報(bào), 21(5): 864-869.

    符淙斌, 葉篤正. 1995. 全球變化和中國(guó)未來(lái)的生存環(huán)境[J]. 大氣科學(xué), 19(1): 110-126.

    付剛, 鐘志明. 2017. 實(shí)驗(yàn)增溫對(duì)西藏高原玉米田土壤呼吸的影響[J]. 生態(tài)環(huán)境學(xué)報(bào), 26(1): 49-54.

    郭維棟, 馬柱國(guó), 王會(huì)軍. 2007. 土壤濕度——一個(gè)跨季度降水預(yù)測(cè)中的重要因子及其應(yīng)用探討[J]. 氣候與環(huán)境研究, 12(1): 20-28.

    胡婭敏, 丁一匯, 廖菲. 2009. 土壤濕度資料同化對(duì)中國(guó)東部夏季區(qū)域氣候模擬的改進(jìn)[J]. 科學(xué)通報(bào), 54(16): 2388-2394.

    蘭鑫宇, 郭子祺, 田野, 等. 2015. 土壤濕度遙感估算同化研究綜述[J]. 地球科學(xué)進(jìn)展, 30(6): 668-679.

    李登宣, 王澄海. 2016. 青藏高原春季土壤濕度與中國(guó)東部夏季降水之間的關(guān)系[J]. 冰川凍土, 38(1): 89-99.

    李潤(rùn)春, 張秀芝, 呂心艷, 等. 2017. 渭河流域土壤濕度異常特征與我國(guó)若干氣候背景的聯(lián)系[J]. 大氣科學(xué)學(xué)報(bào), 40(1): 100-110.

    李哲, 王磊, 王林, 等. 2017. 基于AMSR-E反演青藏高原夏季表層土壤濕度[J]. 高原氣象, 36(1): 67-78.

    馬思源, 朱克云, 李明星, 等. 2016. 中國(guó)區(qū)域多源土壤濕度數(shù)據(jù)的比較研究 [J]. 氣候與環(huán)境研究, 21(2): 121-133.

    沈艷, 馮明農(nóng), 張洪政, 等. 2010. 我國(guó)逐日降水量格點(diǎn)化方法[J]. 應(yīng)用氣象學(xué)報(bào), 21(3): 279-286.

    沈振西, 付剛. 2016. 藏北高原高寒草甸水分利用效率與環(huán)境溫濕度的關(guān)系[J]. 生態(tài)環(huán)境學(xué)報(bào), 23(8): 1259-1263.

    沈振西, 孫維, 李少偉, 等. 2015. 藏北高原不同海拔高度高寒草甸植被指數(shù)與環(huán)境溫濕度的關(guān)系[J]. 生態(tài)環(huán)境學(xué)報(bào), 24(10): 1591-1598.

    孫丞虎, 李維京, 張祖強(qiáng), 等. 2005. 淮河流域土壤濕度異常的時(shí)空分布特征及其與氣候異常關(guān)系的初步研究[J]. 應(yīng)用氣象學(xué)報(bào), 16(2): 129-138.

    王迪, 沈潤(rùn)平, 師春香, 等. 2015. CLM4.5對(duì)中國(guó)區(qū)域夏季土壤濕度的模擬研究, 科學(xué)技術(shù)與工程, 15(17): 1671-1815.

    王靜, 祁莉, 何金海, 等. 2016. 青藏高原春季土壤濕度與我國(guó)長(zhǎng)江流域夏季降水的聯(lián)系及其可能機(jī)理[J]. 地球物理學(xué)報(bào), 59(11): 3985-3995.

    張文君, 周天軍, 宇如聰. 2008. 中國(guó)土壤濕度的分布與變化I.多種資料間的比較[J]. 大氣科學(xué), 32(3): 581-597.

    張秀英, 江洪, 韓英. 2010. 陸面數(shù)據(jù)同化系統(tǒng)及其在全球變化研究中的應(yīng)用[J]. 遙感信息, (4): 135-143.

    周娟, 文軍, 王欣, 等. 2017. 青藏高原季風(fēng)演變及其與土壤濕度的相關(guān)分析[J]. 高原氣象, 36(1): 45-56.

    2017. The relationship between surface soil moisture and precipitation in Huaihe River basin [J]. Ecology and Environmental Sciences, 26(4): 561-569.

    HUANG Yong, QIU Xumin, HUANG Guogui.

    The Relationship between Surface Soil Moisture and Precipitation in Huaihe River Basin

    HUANG Yong1,2, QIU Xumin3, HUANG Guogui1,2
    1. Anhui Meteorology Institute Anhui Key Lab of Atmospheric Science and Remote Sensing, Hefei 230031, China; 2. Shouxian Climatology Observatory, Shouxian 232200, China; 3. Tianjin Binhai New Area Meteorological Bureau, Tianjin, 300457, China

    Soil moisture is a key factor for the energy and matter exchange between land surface and atmosphere. In addition to the field observation, satellite remote sensing is also an important method for the soil moisture measure. In this study, the soil moisture index in the roots region product (SM-DAS-2) which is developed by EUMETSAT Hydrology Satellite Application Facility (HSAF) is applied in Huaihe River Basin, and the applicability of this produce in soil moisture monitoring is analyzed to provide better foundation for the application of this product. In this paper, the ability of monitoring the soil moisture over China is analyzed at the first step. And then, the spatial and temporal characters of soil moisture in Huaihe River basin were studied, and its relationship with the precipitation was analyzed with daily precipitation products from China Meterological Administration (CMA). At last, the relationship of SM-DAS-2 and precipitation is calculated to discuss the response time of soil moisture to precipitation downloaded from CMA. As the results shown, SM-DAS-2 data could embody the spatial characteristics of soil moisture distribution in China. As to Huaihe River Basin, the soil moisture indexes of SM-DAS-2 were similar to the observed data in Meteorology Observatory (correlation coefficient was 0.51, which was tested by T with confidence level of 0.05). But soil moisture indexes of SM-DAS-2 were lower than observation data (annual data of SM-DAS-2 was 0.6, and the observation data was 0.69).The surface soil moisture shows a “Dry North, Wet South” pattern, and it was changed with season. The wet center is located in the mountain and the river areas which are stronger in summer and autumn. The study also shows that the soil moisture was obviously increasing with the rainfall. However, the relationship was less obvious in the mountain and the river regions. In summer, this relationship is even stronger when compared with other seasons. In general, the surface soil moisture response 3 days after the rainfall. And the longest period of rain affecting surface moisture could be 10~15 days later.

    Huaihe River basin; SM-DAS-2; soil moisture; precipitation

    10.16258/j.cnki.1674-5906.2017.04.003

    P933; X16; X14

    A

    1674-5906(2017)04-0561-09

    黃勇, 邱旭敏, 黃國(guó)貴. 2017. 淮河流域表層土壤濕度時(shí)空特征及其與地面降水的關(guān)系[J]. 生態(tài)環(huán)境學(xué)報(bào), 26(4): 561-569.

    國(guó)家自然科學(xué)基金項(xiàng)目(41275030);公益性行業(yè)(氣象)科研專(zhuān)項(xiàng)(GYHY201306040);安徽省自然科學(xué)基金項(xiàng)目(1508085MD64);中國(guó)氣象局項(xiàng)目(CMAGJ2015M28)

    黃勇(1980年生),高級(jí)工程師,博士,主要從事大氣遙感應(yīng)用技術(shù)研發(fā)工作。E-mail: hy121_2000@126.com

    2017-01-06

    猜你喜歡
    淮河流域土壤濕度表層
    淮河流域省級(jí)河湖長(zhǎng)第一次聯(lián)席會(huì)議召開(kāi)
    治淮(2022年8期)2022-09-03 03:42:26
    半潛式平臺(tái)表層卡套管處理與認(rèn)識(shí)
    海洋石油(2021年3期)2021-11-05 07:43:10
    水體表層沉積物對(duì)磷的吸收及釋放研究進(jìn)展
    土壤濕度傳感器在園林綠化灌溉上的應(yīng)用初探
    基于51單片機(jī)控制花盆土壤濕度
    電子制作(2019年15期)2019-08-27 01:12:12
    基于ELECTRE III的淮河流域水資源安全評(píng)價(jià)
    淮委研究部署淮河流域推進(jìn)河長(zhǎng)制工作
    治淮(2017年2期)2017-04-17 11:54:33
    四川盆地土壤濕度時(shí)空分布及影響因子分析
    中國(guó)不同氣候區(qū)土壤濕度特征及其氣候響應(yīng)
    氬弧熔覆原位合成Ti(C,N)-WC增強(qiáng)鎳基表層復(fù)合材料的研究
    焊接(2015年6期)2015-07-18 11:02:25
    天堂影院成人在线观看| 亚洲成人免费av在线播放| 国产精品免费视频内射| 男女下面进入的视频免费午夜 | 黄色毛片三级朝国网站| 亚洲欧美精品综合一区二区三区| av在线播放免费不卡| 国产精品自产拍在线观看55亚洲| 99精品久久久久人妻精品| 伦理电影免费视频| 97人妻天天添夜夜摸| 成人三级黄色视频| 女生性感内裤真人,穿戴方法视频| 夜夜爽天天搞| 中出人妻视频一区二区| 欧美乱色亚洲激情| 亚洲自拍偷在线| 国产一区二区三区在线臀色熟女 | 美女大奶头视频| 美女高潮到喷水免费观看| 一级片免费观看大全| 久久精品国产99精品国产亚洲性色 | 欧美中文日本在线观看视频| 91九色精品人成在线观看| 亚洲专区中文字幕在线| 天堂动漫精品| tocl精华| 人妻丰满熟妇av一区二区三区| 亚洲人成网站在线播放欧美日韩| 国产欧美日韩一区二区三| 69av精品久久久久久| 中文字幕色久视频| 免费在线观看视频国产中文字幕亚洲| 免费看a级黄色片| 国产亚洲欧美在线一区二区| 欧美人与性动交α欧美精品济南到| 最近最新中文字幕大全电影3 | 国产精品 欧美亚洲| 手机成人av网站| 久久国产精品男人的天堂亚洲| 亚洲九九香蕉| 成人国语在线视频| 国产免费男女视频| 天天影视国产精品| 狠狠狠狠99中文字幕| 亚洲成国产人片在线观看| 9191精品国产免费久久| 精品人妻1区二区| 一个人免费在线观看的高清视频| 桃色一区二区三区在线观看| 国产熟女xx| 黑人巨大精品欧美一区二区蜜桃| 一个人免费在线观看的高清视频| 夜夜夜夜夜久久久久| 9191精品国产免费久久| 村上凉子中文字幕在线| 国产亚洲精品久久久久5区| av片东京热男人的天堂| 一本综合久久免费| 精品久久久久久,| 叶爱在线成人免费视频播放| 首页视频小说图片口味搜索| 超碰成人久久| 久久狼人影院| 亚洲五月天丁香| 久久伊人香网站| 国产精品1区2区在线观看.| 久久亚洲精品不卡| 中国美女看黄片| 国产人伦9x9x在线观看| 不卡一级毛片| 日韩三级视频一区二区三区| 18禁黄网站禁片午夜丰满| 精品高清国产在线一区| 免费看a级黄色片| 国产有黄有色有爽视频| 老司机午夜福利在线观看视频| 免费日韩欧美在线观看| 亚洲精品久久成人aⅴ小说| 老鸭窝网址在线观看| 久久伊人香网站| 欧美成人午夜精品| 久久久水蜜桃国产精品网| 69精品国产乱码久久久| 中文字幕高清在线视频| 村上凉子中文字幕在线| 欧美在线一区亚洲| 成人影院久久| 亚洲成人久久性| 琪琪午夜伦伦电影理论片6080| 黄色片一级片一级黄色片| netflix在线观看网站| 亚洲精品一区av在线观看| 欧美激情高清一区二区三区| 欧美成人免费av一区二区三区| 岛国在线观看网站| 精品一品国产午夜福利视频| 大香蕉久久成人网| 免费在线观看日本一区| 成人国产一区最新在线观看| 极品教师在线免费播放| 国产亚洲欧美精品永久| 亚洲精品一卡2卡三卡4卡5卡| 1024香蕉在线观看| 男人操女人黄网站| 在线国产一区二区在线| 国产野战对白在线观看| 日韩av在线大香蕉| 丰满人妻熟妇乱又伦精品不卡| 精品少妇一区二区三区视频日本电影| 亚洲男人天堂网一区| 大香蕉久久成人网| 99国产精品99久久久久| 亚洲国产精品sss在线观看 | 美国免费a级毛片| 国产精品久久久av美女十八| 电影成人av| 久久久久精品国产欧美久久久| 午夜两性在线视频| 91精品三级在线观看| 国产精品美女特级片免费视频播放器 | 欧美乱色亚洲激情| 久久狼人影院| 91国产中文字幕| 欧美最黄视频在线播放免费 | 在线观看免费日韩欧美大片| 久久国产亚洲av麻豆专区| 久久欧美精品欧美久久欧美| 亚洲国产看品久久| 大型av网站在线播放| 多毛熟女@视频| 久久精品91蜜桃| 亚洲精品国产精品久久久不卡| 国产精品美女特级片免费视频播放器 | 欧美人与性动交α欧美软件| 午夜免费观看网址| 欧美日韩视频精品一区| 黑人操中国人逼视频| 亚洲第一av免费看| 99riav亚洲国产免费| 怎么达到女性高潮| 亚洲少妇的诱惑av| av在线天堂中文字幕 | 日韩欧美国产一区二区入口| 久久久久久久精品吃奶| 午夜免费鲁丝| 亚洲国产精品999在线| 久久亚洲真实| 久久青草综合色| av天堂在线播放| 久久99一区二区三区| 欧美大码av| 日本免费一区二区三区高清不卡 | 亚洲成国产人片在线观看| 国产精品98久久久久久宅男小说| 老汉色av国产亚洲站长工具| 性欧美人与动物交配| 欧美日韩精品网址| 欧美日本中文国产一区发布| 精品欧美一区二区三区在线| 热re99久久精品国产66热6| 国产亚洲欧美在线一区二区| 精品久久久久久成人av| 999精品在线视频| 黑丝袜美女国产一区| 久久草成人影院| 国产精品亚洲一级av第二区| 淫秽高清视频在线观看| 精品久久久久久成人av| 精品电影一区二区在线| 男女下面插进去视频免费观看| 国产精品一区二区免费欧美| 久久国产乱子伦精品免费另类| 色精品久久人妻99蜜桃| 久久国产亚洲av麻豆专区| 国产真人三级小视频在线观看| 青草久久国产| av电影中文网址| 婷婷六月久久综合丁香| 精品国产一区二区久久| 国产乱人伦免费视频| 99国产精品免费福利视频| 国产区一区二久久| 黄网站色视频无遮挡免费观看| 国产精品爽爽va在线观看网站 | 午夜福利在线观看吧| 免费在线观看日本一区| 精品一区二区三卡| 别揉我奶头~嗯~啊~动态视频| 又黄又粗又硬又大视频| 91精品三级在线观看| 亚洲成国产人片在线观看| 一区二区三区激情视频| 亚洲人成77777在线视频| 欧洲精品卡2卡3卡4卡5卡区| www.自偷自拍.com| 中国美女看黄片| 看黄色毛片网站| 久久亚洲精品不卡| 精品久久久久久电影网| 亚洲中文字幕日韩| 精品国产乱码久久久久久男人| 91大片在线观看| 最好的美女福利视频网| 国产色视频综合| 无限看片的www在线观看| 国产一卡二卡三卡精品| 国产区一区二久久| 亚洲国产欧美日韩在线播放| 人成视频在线观看免费观看| 看片在线看免费视频| 国产又色又爽无遮挡免费看| 人人妻人人澡人人看| av免费在线观看网站| 视频区欧美日本亚洲| 免费久久久久久久精品成人欧美视频| 纯流量卡能插随身wifi吗| 亚洲aⅴ乱码一区二区在线播放 | 男女下面插进去视频免费观看| 欧美最黄视频在线播放免费 | 最新在线观看一区二区三区| 身体一侧抽搐| 母亲3免费完整高清在线观看| 精品高清国产在线一区| 日韩欧美三级三区| 狂野欧美激情性xxxx| 老司机靠b影院| 99国产精品一区二区蜜桃av| 热re99久久国产66热| 在线永久观看黄色视频| 欧美黄色淫秽网站| 不卡一级毛片| 老鸭窝网址在线观看| 色播在线永久视频| 免费av毛片视频| 99久久精品国产亚洲精品| 欧美成人午夜精品| 欧美成人性av电影在线观看| 99国产精品99久久久久| 免费在线观看完整版高清| 老司机午夜福利在线观看视频| 每晚都被弄得嗷嗷叫到高潮| 最近最新免费中文字幕在线| 午夜91福利影院| 美女高潮喷水抽搐中文字幕| 日韩欧美在线二视频| 亚洲精品在线观看二区| 99在线视频只有这里精品首页| 天堂俺去俺来也www色官网| 亚洲精品美女久久久久99蜜臀| 精品欧美一区二区三区在线| 乱人伦中国视频| 日韩大码丰满熟妇| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久国产一区二区| 天天躁夜夜躁狠狠躁躁| 久久久久九九精品影院| 女性被躁到高潮视频| 午夜福利免费观看在线| 黄色怎么调成土黄色| 在线观看一区二区三区| 操美女的视频在线观看| 999久久久国产精品视频| 亚洲一区二区三区欧美精品| 又黄又粗又硬又大视频| 久久久久久大精品| 我的亚洲天堂| 男女做爰动态图高潮gif福利片 | 国产精品一区二区免费欧美| 日韩欧美一区视频在线观看| 首页视频小说图片口味搜索| 国产深夜福利视频在线观看| 99在线视频只有这里精品首页| 一级a爱视频在线免费观看| 午夜精品在线福利| 视频区欧美日本亚洲| 久久午夜亚洲精品久久| 欧洲精品卡2卡3卡4卡5卡区| 国产精品国产高清国产av| 女人高潮潮喷娇喘18禁视频| 欧美日韩亚洲综合一区二区三区_| 99久久精品国产亚洲精品| 久久久久国产精品人妻aⅴ院| 18禁国产床啪视频网站| 国产熟女午夜一区二区三区| 精品日产1卡2卡| 欧美人与性动交α欧美软件| 日本五十路高清| 免费高清在线观看日韩| 精品一区二区三区四区五区乱码| 国产精品久久久久成人av| 一区福利在线观看| 大香蕉久久成人网| 99久久99久久久精品蜜桃| 男女下面进入的视频免费午夜 | 中亚洲国语对白在线视频| 两人在一起打扑克的视频| 又黄又爽又免费观看的视频| 91麻豆精品激情在线观看国产 | 一进一出抽搐gif免费好疼 | 日韩免费高清中文字幕av| 男女下面进入的视频免费午夜 | 曰老女人黄片| 亚洲av美国av| 免费在线观看亚洲国产| 国产精品国产高清国产av| 天堂动漫精品| 很黄的视频免费| 手机成人av网站| 在线观看一区二区三区| 香蕉丝袜av| 国产区一区二久久| 欧美日韩国产mv在线观看视频| 日韩免费av在线播放| 日本黄色日本黄色录像| 天堂中文最新版在线下载| 亚洲人成伊人成综合网2020| 欧美精品啪啪一区二区三区| 丰满的人妻完整版| 中文字幕人妻熟女乱码| av电影中文网址| 日韩欧美一区视频在线观看| 国产熟女午夜一区二区三区| 亚洲专区字幕在线| 一进一出好大好爽视频| 美女 人体艺术 gogo| 又黄又爽又免费观看的视频| 日本精品一区二区三区蜜桃| e午夜精品久久久久久久| 少妇 在线观看| 19禁男女啪啪无遮挡网站| 久久人人爽av亚洲精品天堂| 日本精品一区二区三区蜜桃| 人妻丰满熟妇av一区二区三区| 免费在线观看完整版高清| 久久久久精品国产欧美久久久| 亚洲精品久久成人aⅴ小说| 亚洲欧美精品综合久久99| 色在线成人网| 成年女人毛片免费观看观看9| 最好的美女福利视频网| 亚洲一区高清亚洲精品| 美女大奶头视频| 女人被狂操c到高潮| 99久久人妻综合| 日韩有码中文字幕| 天天躁夜夜躁狠狠躁躁| 日本wwww免费看| 嫩草影视91久久| bbb黄色大片| 纯流量卡能插随身wifi吗| 国产熟女午夜一区二区三区| 国产欧美日韩一区二区三区在线| 99久久国产精品久久久| 日本五十路高清| 老司机午夜十八禁免费视频| 久久人人精品亚洲av| 亚洲免费av在线视频| 午夜免费成人在线视频| 女人被狂操c到高潮| 久久人人精品亚洲av| 桃色一区二区三区在线观看| 亚洲国产精品999在线| 亚洲一区高清亚洲精品| 国产xxxxx性猛交| 中出人妻视频一区二区| 在线观看免费视频网站a站| 成人亚洲精品av一区二区 | 中文字幕最新亚洲高清| 国产精品一区二区三区四区久久 | 如日韩欧美国产精品一区二区三区| 黄网站色视频无遮挡免费观看| 久久久久九九精品影院| 亚洲中文字幕日韩| 麻豆国产av国片精品| 国产一区二区三区视频了| 后天国语完整版免费观看| 每晚都被弄得嗷嗷叫到高潮| 国产精品秋霞免费鲁丝片| 国产精品影院久久| 最新在线观看一区二区三区| 日韩 欧美 亚洲 中文字幕| 99久久综合精品五月天人人| 精品人妻在线不人妻| 欧美中文综合在线视频| 午夜免费观看网址| 免费观看精品视频网站| 亚洲中文av在线| 国产成人精品久久二区二区免费| av在线天堂中文字幕 | 欧美日韩中文字幕国产精品一区二区三区 | 性欧美人与动物交配| 国产精品久久久av美女十八| √禁漫天堂资源中文www| 男女做爰动态图高潮gif福利片 | 欧美乱码精品一区二区三区| 免费女性裸体啪啪无遮挡网站| 久久精品亚洲熟妇少妇任你| 老鸭窝网址在线观看| 亚洲五月色婷婷综合| 亚洲五月婷婷丁香| 日本免费一区二区三区高清不卡 | 好看av亚洲va欧美ⅴa在| 欧美日韩乱码在线| 水蜜桃什么品种好| 欧美乱妇无乱码| 老司机靠b影院| 黄色片一级片一级黄色片| 精品人妻1区二区| 久久久国产成人精品二区 | 欧美久久黑人一区二区| 国产三级黄色录像| 美国免费a级毛片| 最近最新中文字幕大全免费视频| 亚洲精品一卡2卡三卡4卡5卡| a在线观看视频网站| 亚洲一区二区三区色噜噜 | 一边摸一边做爽爽视频免费| 波多野结衣一区麻豆| 国产高清视频在线播放一区| 免费在线观看亚洲国产| 午夜福利欧美成人| 精品第一国产精品| 久久久国产精品麻豆| 在线观看舔阴道视频| 亚洲av第一区精品v没综合| 欧美日韩精品网址| tocl精华| 精品人妻1区二区| 中出人妻视频一区二区| 久久精品国产清高在天天线| 丁香六月欧美| 中文字幕人妻丝袜制服| 热re99久久精品国产66热6| 日韩有码中文字幕| 久久欧美精品欧美久久欧美| 国产区一区二久久| 男女做爰动态图高潮gif福利片 | 国产精品乱码一区二三区的特点 | 天堂影院成人在线观看| 国产午夜精品久久久久久| 亚洲国产欧美网| 国产亚洲精品久久久久5区| 国产乱人伦免费视频| 久久天堂一区二区三区四区| 日韩欧美在线二视频| 亚洲,欧美精品.| 脱女人内裤的视频| 国产有黄有色有爽视频| 一级a爱视频在线免费观看| 香蕉国产在线看| 黑人巨大精品欧美一区二区蜜桃| 婷婷丁香在线五月| 免费在线观看完整版高清| 国产1区2区3区精品| 中文字幕人妻丝袜制服| 在线观看66精品国产| 制服诱惑二区| 男女午夜视频在线观看| 黑人猛操日本美女一级片| 88av欧美| 咕卡用的链子| 免费搜索国产男女视频| 老汉色∧v一级毛片| 午夜老司机福利片| 久久狼人影院| 日韩大尺度精品在线看网址 | 操美女的视频在线观看| 日韩欧美一区视频在线观看| 久久青草综合色| 国产精品av久久久久免费| 亚洲av五月六月丁香网| 一区在线观看完整版| 婷婷精品国产亚洲av在线| 91麻豆av在线| 国产精品二区激情视频| 美国免费a级毛片| 色精品久久人妻99蜜桃| 999久久久国产精品视频| 亚洲第一欧美日韩一区二区三区| 精品乱码久久久久久99久播| 黑人操中国人逼视频| 久久亚洲真实| 纯流量卡能插随身wifi吗| 欧美中文日本在线观看视频| 亚洲第一青青草原| 色综合婷婷激情| 欧美在线一区亚洲| 12—13女人毛片做爰片一| 女人被狂操c到高潮| 国产欧美日韩一区二区三| 黄色毛片三级朝国网站| 99国产精品一区二区蜜桃av| 丝袜在线中文字幕| www国产在线视频色| 国产精品免费一区二区三区在线| 精品国产国语对白av| 啪啪无遮挡十八禁网站| 纯流量卡能插随身wifi吗| 中出人妻视频一区二区| 深夜精品福利| 亚洲自偷自拍图片 自拍| 嫁个100分男人电影在线观看| 久久欧美精品欧美久久欧美| 色尼玛亚洲综合影院| 国产一区二区三区综合在线观看| 欧美性长视频在线观看| 欧美乱妇无乱码| 午夜影院日韩av| 搡老乐熟女国产| 淫秽高清视频在线观看| 黄色成人免费大全| 无人区码免费观看不卡| 日韩欧美一区视频在线观看| 一区二区三区国产精品乱码| 亚洲国产欧美一区二区综合| 91字幕亚洲| 久久久国产一区二区| 长腿黑丝高跟| 日韩精品中文字幕看吧| 波多野结衣一区麻豆| 天天影视国产精品| 免费在线观看日本一区| 大型av网站在线播放| 亚洲成人久久性| 亚洲欧美精品综合久久99| 国产精品一区二区三区四区久久 | 一级片'在线观看视频| 中文字幕人妻丝袜一区二区| 亚洲自偷自拍图片 自拍| 成人手机av| 午夜影院日韩av| xxx96com| 久久人妻熟女aⅴ| 在线观看www视频免费| 国产亚洲欧美98| 国产一卡二卡三卡精品| 亚洲一区二区三区不卡视频| 中文欧美无线码| 成人三级黄色视频| 看免费av毛片| 免费搜索国产男女视频| 日韩精品青青久久久久久| 黑人猛操日本美女一级片| 国产精品爽爽va在线观看网站 | 欧美久久黑人一区二区| 国产伦人伦偷精品视频| 亚洲精品中文字幕在线视频| 老司机福利观看| 国产一区二区在线av高清观看| 久久久久国产精品人妻aⅴ院| 一级作爱视频免费观看| 一区二区日韩欧美中文字幕| 久久久久久久久中文| 中国美女看黄片| 新久久久久国产一级毛片| 9色porny在线观看| 亚洲国产精品一区二区三区在线| ponron亚洲| 多毛熟女@视频| 精品久久久久久久毛片微露脸| 一边摸一边做爽爽视频免费| 久久久国产成人精品二区 | 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品一区二区精品视频观看| 日本一区二区免费在线视频| 一区二区日韩欧美中文字幕| 久久精品国产综合久久久| 久久久久九九精品影院| 亚洲国产精品999在线| av国产精品久久久久影院| 母亲3免费完整高清在线观看| 日韩欧美国产一区二区入口| 成人18禁高潮啪啪吃奶动态图| 黄频高清免费视频| 国产欧美日韩一区二区三区在线| 在线播放国产精品三级| 天天添夜夜摸| 天堂俺去俺来也www色官网| 久久久久国产一级毛片高清牌| 欧美激情极品国产一区二区三区| 亚洲色图av天堂| 美女 人体艺术 gogo| 精品乱码久久久久久99久播| www日本在线高清视频| 欧美日韩一级在线毛片| 男女做爰动态图高潮gif福利片 | 亚洲人成网站在线播放欧美日韩| 欧美中文日本在线观看视频| 啦啦啦 在线观看视频| 在线观看66精品国产| 首页视频小说图片口味搜索| 19禁男女啪啪无遮挡网站| 侵犯人妻中文字幕一二三四区| 性欧美人与动物交配| 国产一区二区在线av高清观看| 亚洲欧美日韩另类电影网站| www国产在线视频色| 琪琪午夜伦伦电影理论片6080| 亚洲欧美激情在线| 色综合欧美亚洲国产小说| 91成年电影在线观看| 亚洲精品在线美女| 午夜久久久在线观看| 国产一区二区三区在线臀色熟女 | 欧美中文日本在线观看视频| 欧美激情极品国产一区二区三区| 一级黄色大片毛片| 中出人妻视频一区二区| 午夜福利免费观看在线| 国产无遮挡羞羞视频在线观看| 国产亚洲欧美精品永久| 亚洲性夜色夜夜综合| 女人被躁到高潮嗷嗷叫费观| 中文字幕最新亚洲高清|