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

    1998-2013年釜溪河流域植被覆蓋時(shí)空演變及其驅(qū)動(dòng)因子研究

    2017-11-02 08:26:28王二麗劉小芳候亞妮陳光建萬(wàn)新
    關(guān)鍵詞:氣候因子溪河氣溫

    王二麗, 劉小芳, 候亞妮,陳光建,萬(wàn)新

    (1.四川理工學(xué)院計(jì)算機(jī)學(xué)院,四川自貢643000;2.四川理工學(xué)院江河流域生態(tài)環(huán)境的集成感知與應(yīng)用院士(專(zhuān)家)工作站,四川自貢643000)

    1998-2013年釜溪河流域植被覆蓋時(shí)空演變及其驅(qū)動(dòng)因子研究

    王二麗1,2, 劉小芳1,2, 候亞妮1,2,陳光建1,萬(wàn)新1

    (1.四川理工學(xué)院計(jì)算機(jī)學(xué)院,四川自貢643000;2.四川理工學(xué)院江河流域生態(tài)環(huán)境的集成感知與應(yīng)用院士(專(zhuān)家)工作站,四川自貢643000)

    植被覆蓋是評(píng)價(jià)生態(tài)環(huán)境的重要指標(biāo)之一?;?998-2013年SPOT VEGETATION數(shù)據(jù)集,采用趨勢(shì)分析、偏相關(guān)分析及計(jì)算Hurst指數(shù)等方法,反演了釜溪河流域NDVI分布的時(shí)空格局及其變化趨勢(shì);并結(jié)合該地區(qū)同期降水量和溫度數(shù)據(jù),利用相關(guān)系數(shù)法對(duì)釜溪河流域植被覆蓋變化的驅(qū)動(dòng)因子進(jìn)行分析。研究表明1998-2013年間釜溪河流域植被覆蓋表現(xiàn)出上升趨勢(shì),增速為0.24%/16a;Hurst指數(shù)計(jì)算結(jié)果顯示,流域內(nèi)植被反向變化特征強(qiáng)于同向特征,結(jié)合NDVI趨勢(shì)圖,未來(lái)植被改善區(qū)域面積所占比例為60.98%,持續(xù)退化面積所占比重為5.97%,主要分布于自貢市南部及其他環(huán)繞城鎮(zhèn)地區(qū);植被覆蓋變化主要以非氣候因子驅(qū)動(dòng)型為主,其面積所占比例為88.48%,而影響植被生長(zhǎng)的氣候因子中氣溫的作用強(qiáng)于降水的作用。

    植被覆蓋變化;Hurst指數(shù);偏相關(guān)性分析;驅(qū)動(dòng)因子;釜溪河流域

    引言

    植被不僅是陸地生態(tài)系統(tǒng)的重要組成部分,也是水循環(huán)、生物化學(xué)等能量循環(huán)的重要媒介,在減少溫室氣體排放、調(diào)節(jié)碳平衡、維持氣候穩(wěn)定方面具有十分重要的作用[1-2],由于其對(duì)氣候變化的敏感性,植被覆蓋動(dòng)態(tài)變化已成為全球陸地生態(tài)系統(tǒng)監(jiān)測(cè)的關(guān)鍵問(wèn)題之一[3]。目前國(guó)內(nèi)外學(xué)者針對(duì)不同地域,利用植被指數(shù)對(duì)植被相關(guān)變化及其與氣候因子的關(guān)系進(jìn)行了深入的探討[4-7],研究結(jié)果均表明氣候?qū)χ脖簧L(zhǎng)起到了重要作用,但部分地區(qū)植被變化的決定性因子與人類(lèi)活動(dòng)密切相關(guān)。毛德華[8]等研究表明溫度是影響東北地區(qū)不同植被類(lèi)型NDVI最重要的因子,而嚴(yán)恩萍[9]等研究成果表明對(duì)于京津風(fēng)沙源重點(diǎn)治理區(qū)植被覆蓋上升最主要的驅(qū)動(dòng)因素是人類(lèi)活動(dòng)。不同區(qū)域植被變化驅(qū)動(dòng)力通常有所不同,大多數(shù)研究在選擇研究區(qū)域時(shí)主要集中于省市[10]、各類(lèi)生態(tài)區(qū)[11-12]、國(guó)家[13]及全球[14]等更為宏觀(guān)的層面,但對(duì)以流域?yàn)閱卧M(jìn)行研究的并不多見(jiàn)[15]。

    釜溪河流經(jīng)四川省威遠(yuǎn)縣、自貢市,是沿途區(qū)域農(nóng)業(yè)灌溉、工業(yè)及生活用水的主要來(lái)源,而近些年隨著工業(yè)建設(shè),流域生態(tài)環(huán)境狀況愈下,部分區(qū)域植被退化嚴(yán)重,亟需快速監(jiān)測(cè)手段進(jìn)行應(yīng)對(duì)管理。為探究釜溪河流域植被覆蓋情況,本文采用SPOT VEGETATION數(shù)據(jù)集,研究流域1998-2013年植被覆蓋變化,同時(shí)結(jié)合氣溫、降水?dāng)?shù)據(jù),分析流域植被氣候變化的響應(yīng)特征,并嘗試對(duì)植被變化的驅(qū)動(dòng)力進(jìn)行分析,旨在為該地區(qū)生態(tài)系統(tǒng)的恢復(fù)、穩(wěn)定及可持續(xù)發(fā)展提供科學(xué)的依據(jù)。

    1 研究區(qū)域概況

    釜溪河在沱江下游右岸,由旭水河和威遠(yuǎn)河在自貢市境內(nèi)鳳凰壩雙河口處匯成干流,河長(zhǎng)190公里,天然落差360余米,于富順縣釜溪口匯入沱江。流域面積3490平方公里,地勢(shì)自西北向東南傾斜。釜溪河流域?qū)儆趤啛釒貪駳夂?,具有冬暖、春早、夏熱、秋涼四季分明的特點(diǎn),年平均氣溫約為18 ℃,年平均降水量為1000 mm~1080 mm,集中在6-9月份。研究區(qū)域如圖1所示。

    圖1 研究區(qū)域

    2 研究方法

    2.1數(shù)據(jù)來(lái)源

    研究所選NDVI數(shù)據(jù)集來(lái)源于SPOT VGT數(shù)據(jù)分發(fā)中心,時(shí)間跨度為2000年到2013年。根據(jù)公式NDVI=DN×0.004-0.1將像元值轉(zhuǎn)換為對(duì)應(yīng)的NDVI值,并采用最大值合成法對(duì)產(chǎn)品進(jìn)行處理,獲取逐月最大NDVI。降水和氣溫?cái)?shù)據(jù)來(lái)源于中國(guó)氣象數(shù)據(jù)共享網(wǎng)月平均溫度和月降水資料,進(jìn)行Kriging空間插值,獲取與NDVI數(shù)據(jù)集像元大小一致、投影相同的氣象柵格圖像,并通過(guò)數(shù)據(jù)裁剪獲得流域月平均溫度和月降水量數(shù)據(jù)。

    2.2趨勢(shì)分析

    一元線(xiàn)性回歸分析針對(duì)每個(gè)柵格點(diǎn)逐像元擬合NDVI的斜率,反映植被覆蓋時(shí)空演變特征[16]。計(jì)算公式為:

    (1)

    式中:n為監(jiān)測(cè)時(shí)間序列的長(zhǎng)度,NDVIi第i年NDVI的均值;slope是趨勢(shì)線(xiàn)增加的斜率。

    2.3偏相關(guān)分析

    地理要素之間的相關(guān)關(guān)系密切程度的測(cè)定,主要是通過(guò)對(duì)相關(guān)系數(shù)的計(jì)算與檢驗(yàn)完成。用以度量除卻其他因素影響,而研究單獨(dú)兩個(gè)要素之間的相互關(guān)系的密切程度,稱(chēng)為偏相關(guān)系數(shù)[17]。計(jì)算公式為

    (2)

    式中:rxy,z為變量固定后變量x和y的偏相關(guān)系數(shù),rxy、rxz、ryz分別表示變量x和y、變量x和變量z、變量y和變量z的相關(guān)系數(shù)。

    2.4復(fù)相關(guān)分析

    多要素間相關(guān)程度的測(cè)定通過(guò)復(fù)相關(guān)系數(shù)的計(jì)算與檢驗(yàn)實(shí)現(xiàn),復(fù)相關(guān)系數(shù)計(jì)算公式為:

    (3)

    式中:Rx,yz表示因變量x和自變量y,z的復(fù)相關(guān)系數(shù),rxz,y表示固定變量后y變量x和z的偏相關(guān)系數(shù),rxy表示變量x和變量y的相關(guān)系數(shù)。

    2.5 Hurst指數(shù)及R/S分析

    Hurst指數(shù)是描述自然界自相似性和長(zhǎng)期依賴(lài)性最有效的方法[18],由Hurst提出的重標(biāo)極差法(R/S)是目前估測(cè)Hurst指數(shù)最常用的方法[19],其計(jì)算主要過(guò)程[20]:

    (1)將時(shí)間序列{ξ(τ)}(τ=1,2,3,...,n)分解為τ個(gè)子時(shí)間序列X(t),其中t=1,...,τ。

    (2)計(jì)算每個(gè)子時(shí)間序列的均值:

    (4)

    其中,τ=1,2,3,...,n。

    (3)計(jì)算累計(jì)離差:

    (5)

    其中,t=1,...,τ。

    (4)計(jì)算極差:

    R(τ)=maxX(t,τ)-minX(t,τ)

    (6)

    其中,t=1,...,τ,τ=1,2,3,...,n。

    (5)計(jì)算時(shí)間序列的標(biāo)準(zhǔn)差:

    (7)

    其中,τ=1,2,3,...,n。

    (6)求取H:

    (8)

    根據(jù)Hurst及Mandelbrot等人的研究[20-21],Hurst指數(shù)H的值域范圍在0~1之間,當(dāng)H值等于0.5時(shí),表明所研究的時(shí)間序列沒(méi)有持續(xù)性,為隨機(jī)序列;當(dāng)H大于0.5時(shí),該時(shí)間序列表現(xiàn)出一致性,而H小于0.5時(shí),則呈現(xiàn)相反的跡象。

    3 結(jié)果與分析

    3.1植被覆蓋時(shí)間變化特征

    1998-2013年釜溪河NDVI變化趨勢(shì)如圖2所示,研究區(qū)域內(nèi)植被覆蓋表現(xiàn)出上升趨勢(shì),但并不顯著(R2=0.1001)。2013年較1998年相比NDVI值增加0.038,平均每年以0.0024速度增加。1998年到2004年植被變化較為劇烈,呈波浪式推進(jìn),2004年達(dá)到最高值0.752,之后幾年植被覆蓋略有起伏,2011年下降較為明顯,整體的植被覆蓋呈現(xiàn)良好狀態(tài)。

    圖2 1998-2013年釜溪河流域NDVI值變化趨勢(shì)

    3.2植被覆蓋空間變化特征

    3.2.1 植被空間分布特征

    以每年最大NDVI值計(jì)算得到釜溪河流域NDVI均值,空間分布圖如圖3(a)所示,16年均值為0.66,最大值為0.73。釜溪河流域植被覆蓋格局整體呈現(xiàn)片狀聚集的分布特征,高值區(qū)主要分布于流域南部、威遠(yuǎn)河上游及自貢市北部地區(qū),環(huán)繞旭水河、威遠(yuǎn)河中間夾集區(qū)也有較高的NDVI值分布,由于這些區(qū)域主要以耕地、林地為主,植被長(zhǎng)勢(shì)較好,NDVI值較大;而低值區(qū)則主要分布于自貢市區(qū)、威遠(yuǎn)河中游(威遠(yuǎn)縣城)及旭水河中上游(榮縣縣城)部分地區(qū),該區(qū)域土地利用類(lèi)型均以城鎮(zhèn)用地為主,植被NDVI值較小。

    圖3 1998-2013釜溪河流域NDVI空間分布及變化趨勢(shì)圖

    為監(jiān)測(cè)釜溪河流域植被覆蓋變化特征,計(jì)算1998-2013年植被NDVI趨勢(shì),并將結(jié)果進(jìn)行標(biāo)準(zhǔn)差分類(lèi),NDVI變化趨勢(shì)空間分布圖如圖3(b)所示。明顯改善區(qū)域主要分布于流域東北部,該地區(qū)地形較高,受人類(lèi)活動(dòng)干預(yù)較少,特別是近年退耕還林等政策的實(shí)施使得區(qū)域植被覆蓋度得到明顯的提高。而嚴(yán)重退化區(qū)域與植被覆蓋度較低區(qū)域相一致,主要以環(huán)繞城鎮(zhèn)地區(qū)為主。中度退化及輕微退化區(qū)域則集中于旭水河中下游、威遠(yuǎn)河中下游及釜溪河中下游地區(qū)。釜溪河流域NDVI變化趨勢(shì)統(tǒng)計(jì)結(jié)果見(jiàn)表1,結(jié)果表明流域植被覆蓋所改善的面積大于退化面積,整體改善呈上升趨勢(shì),這與文中3.1所述植被時(shí)間變化特征相一致。流域植被覆蓋所改善的面積(27.44%)大于退化面積(25.74%),流域內(nèi)保持不變的面積所占比重46.82%。

    表1 釜溪河流域NDVI變化趨勢(shì)

    3.2.2流域未來(lái)趨勢(shì)特征

    上述分析主要針對(duì)近16年來(lái)釜溪河流域植被空間和時(shí)間上的分布特征,但未來(lái)趨勢(shì)變化特征并不明朗,為此進(jìn)一步計(jì)算植被NDVI的Hurst指數(shù),以指示其未來(lái)變化趨勢(shì),Hurst空間分布圖如圖4(a)所示。釜溪河流域植被NDVI的Hurst指數(shù)值域區(qū)間為[0.22,0.75],平均值為0.46,其中Hurst指數(shù)大于0.5的像元數(shù)所占百分比為33.24%,表明研究區(qū)域植被反向變化特征強(qiáng)于同向特征。將釜溪河流域Hurst指數(shù)與趨勢(shì)slope圖進(jìn)行疊加,NDVI變化特征如圖4(b)所示,統(tǒng)計(jì)顯示流域內(nèi)由退化區(qū)域轉(zhuǎn)為改善區(qū)域所占面積為46.23%,主要分布于旭水河、威遠(yuǎn)河、釜溪河中下游地區(qū)。而由改善轉(zhuǎn)為退化區(qū)域所占面積為32.99%,持續(xù)退化面積所占比重為5.97%,主要分布于自貢市南部及其它環(huán)繞城鎮(zhèn)地區(qū)。持續(xù)增加面積所占比重為14.75%,以流域邊緣地區(qū)分布居多。

    圖4 Hurst指數(shù)空間分布圖及NDVI變化特征

    3.3氣溫、降水變化特征

    1998-2013年釜溪河流域平均年氣溫為18.02 ℃,平均年降水量為9608.23 mm,流域平均氣溫變化及降雨量變化如圖5所示。流域氣溫和降水總體呈現(xiàn)遞增趨勢(shì),但趨勢(shì)并不明顯,而年平均氣溫增長(zhǎng)顯著程度(R2=0.2671)遠(yuǎn)大于年降水量增長(zhǎng)程度(R2=0.0021)。平均氣溫總體呈現(xiàn)上升趨勢(shì),從2000年最低氣溫(17.11 ℃)至2002年持續(xù)增長(zhǎng),2003年略有下降后流域氣溫一直遞增到2006年出現(xiàn)峰值(18.94 ℃),2009年出現(xiàn)次高溫后氣溫略有下降。年降雨量變化較為反復(fù),2013年降雨量為頂峰達(dá)到13104.78 mm,2006年降雨量為低谷值7080.64 mm,其余時(shí)段降雨量波動(dòng)較為劇烈。

    圖5 1998-2013年釜溪河流域平均氣溫變化及降雨量變化

    3.4 NDVI與氣候間的關(guān)系

    根據(jù)研究區(qū)域氣溫、降水?dāng)?shù)據(jù),計(jì)算NDVI與氣溫和降水的偏相關(guān)系數(shù),偏相關(guān)系數(shù)空間分布如圖6所示。結(jié)果顯示,NDVI與氣溫的相關(guān)性為0.532(-0.58~0.88),與降水的相關(guān)性為0.28(-0.60~0.72)。盡管流域內(nèi)NDVI與氣溫和降水的相關(guān)性并不高,但NDVI與氣溫的偏相關(guān)系數(shù)大于NDVI與降水的偏相關(guān)系數(shù),這在一定程度表明釜溪河流域植被生長(zhǎng)受氣溫的影響要強(qiáng)與降水的影響。

    圖6 NDVI與氣溫、降水偏相關(guān)系數(shù)空間分布

    流域內(nèi)NDVI最高值出現(xiàn)在2004年,最低值出現(xiàn)在2000年,而同期的氣溫?cái)?shù)據(jù)最低值與NDVI最低值相符,均出現(xiàn)在2000年,但氣溫最高值與NDVI最高值年份有所差異,出現(xiàn)在2006年。NDVI在2006年后呈現(xiàn)緩慢波動(dòng)上升趨勢(shì),但在2011年陡然下降,該年份降水量出現(xiàn)一個(gè)較小值,但氣溫與同期相比并沒(méi)有太大差異。這在一定程度上表明流域內(nèi)植被覆蓋變化并不僅僅是氣候因素主導(dǎo)作用的結(jié)果,同時(shí)NDVI與氣溫和降水的偏相關(guān)系數(shù)較低也能從側(cè)面闡述氣候因子與NDVI變化趨勢(shì)并不一致的現(xiàn)象。

    為進(jìn)一步闡明氣候因子與NDVI的關(guān)系,計(jì)算復(fù)相關(guān)系數(shù)以反映NDVI與降水、溫度間的復(fù)相關(guān)程度,復(fù)相關(guān)系數(shù)空間分布結(jié)果如圖7(a)所示。釜溪河流域平均復(fù)相關(guān)系數(shù)為0.34,復(fù)相關(guān)系數(shù)較大區(qū)域主要分布于威遠(yuǎn)河中上游。同時(shí)參照陳云浩、王永財(cái)[22-23]等研究植被覆蓋變化驅(qū)動(dòng)分區(qū)原則及標(biāo)準(zhǔn),對(duì)釜溪河流域植被覆蓋變化進(jìn)行驅(qū)動(dòng)分區(qū)。流域植被覆蓋驅(qū)動(dòng)分區(qū)結(jié)果如圖7(b)所示,10.1%區(qū)域?qū)儆跉鉁仳?qū)動(dòng)型,主要呈片狀集中于流域東北部、威遠(yuǎn)河上游地帶。降水驅(qū)動(dòng)型面積僅占流域面積的0.38%,零星狀分布于旭水河中下游地帶。降水和氣溫驅(qū)動(dòng)型所占流域面積的1.04%,主要位于威遠(yuǎn)河中部?jī)缮綂A集處。而非氣候因子驅(qū)動(dòng)型所占比例最大為88.48%,基本覆蓋整個(gè)流域。以上分析表明,釜溪河流域植被覆蓋動(dòng)態(tài)變化在一定程度、一定區(qū)域范圍內(nèi)受到氣候因素的變異,但主要驅(qū)動(dòng)因素還是非氣候因子主導(dǎo)作用。已有研究表明退耕還林等生態(tài)工程的實(shí)施,使得流域植被覆蓋率上升,而城市擴(kuò)展將導(dǎo)致城市周?chē)脖辉獾狡茐闹脖桓采w率降低[24],這與釜溪河流域城鎮(zhèn)NDVI值較小且有持續(xù)惡化趨勢(shì)相一致。

    圖7 復(fù)相關(guān)系數(shù)分布及驅(qū)動(dòng)力分區(qū)圖

    4 結(jié)論與討論

    本文以SPOT VEGETATION數(shù)據(jù)集,輔以趨勢(shì)分析、Hurst指數(shù)及偏相關(guān)分析等方法,反演釜溪河流域1998-2013年NDVI分布時(shí)間、空間格局及其變化趨勢(shì),并結(jié)合該地區(qū)同期降水量和溫度數(shù)據(jù),采用相關(guān)系數(shù)法對(duì)釜溪河流域植被覆蓋變化的驅(qū)動(dòng)因子進(jìn)行分析。研究表明釜溪河流域植被覆蓋表現(xiàn)出上升、改善趨勢(shì)。流域植被反向變化特征強(qiáng)于同向特征,植被生長(zhǎng)受氣溫的影響強(qiáng)于降水影響,而植被覆蓋驅(qū)動(dòng)因素分區(qū)中非氣候因子驅(qū)動(dòng)型所占比例為88.48%。

    2001年,自貢市委市政府提出“東進(jìn)南移”城市發(fā)展戰(zhàn)略,隨著城市建設(shè)的不斷推進(jìn),使得自貢市南部和東部地區(qū)呈現(xiàn)由改善區(qū)域轉(zhuǎn)為退化區(qū)域及持續(xù)退化區(qū)域的趨勢(shì)。城市發(fā)展早期由于缺少對(duì)環(huán)境的保護(hù)易導(dǎo)致植被遭到較嚴(yán)重的破壞,相關(guān)部門(mén)應(yīng)該做好城市可持續(xù)發(fā)展的規(guī)劃,增加園林綠化面積,加強(qiáng)對(duì)環(huán)境的治理,提升植被覆蓋度。城鎮(zhèn)化進(jìn)程的加快促進(jìn)人口的聚集,而人類(lèi)對(duì)于生態(tài)環(huán)境的干預(yù)程度遠(yuǎn)高于氣候因子的作用,更應(yīng)在保護(hù)環(huán)境中發(fā)揮正向主導(dǎo)作用。針對(duì)植被覆蓋度持續(xù)下降區(qū)域,政府部門(mén)應(yīng)積極實(shí)施生態(tài)工程,引導(dǎo)農(nóng)民合理種植農(nóng)作物,退耕還林,擴(kuò)大植被覆蓋率,防止水土流失及生態(tài)環(huán)境的進(jìn)一步惡化。

    植被覆蓋變化監(jiān)測(cè)是一個(gè)長(zhǎng)期過(guò)程,文中僅針對(duì)釜溪河流域近16年植被空間格局變化進(jìn)行了探討,如有更長(zhǎng)時(shí)間序列數(shù)據(jù)則可進(jìn)一步揭示流域植被時(shí)空變化及其歸因。植被覆蓋變化必然受到除氣候之外的因素影響,本文中僅僅對(duì)氣候因素進(jìn)行了重點(diǎn)闡述和說(shuō)明,并未明確其他因素在流域植被動(dòng)態(tài)變化中所產(chǎn)生的影響,如何定量區(qū)分其他因素如人類(lèi)活動(dòng)對(duì)植被空間格局的變化及其影響仍是進(jìn)一步需要探討的內(nèi)容。

    [1] GODINEZ-ALVAREZ H,HERRICK J,MATTOCKS M,et al.Comparison of Three Vegetation Monitoring Methods:Their Relative Utility for Ecological Assessment and Monitoring[J].Ecological Indicators,2009,9(5):1001-1008.

    [2] PENG J,LIU Z,LIU Y,et al.Trend Analysis of Vegetation Dynamics in Qinghai-Tibet Plateau Using Hurst Exponent[J].Ecological Indicators,2012,14(1):28-39.

    [3] KELLY M,TUXEN K A,STRALBERG D.Mapping Changes to Vegetation Pattern in a Restoring Wetland:Finding Pattern Metrics that are Consistent Across Spatial Scale and Time[J].Ecological Indicators,2011,11(2):263-273.

    [4] PARK H S,SOHN B.Recent Rrends in Changes of Vegetation Over East Asia Coupled with Temperature and Rainfall Variations[J].Journal of Geophysical Research:Atmospheres(1984-2012),2010,115(D14).

    [5] ZHAO J,ZHANG H,ZHANG Z,et al.Spatial and Temporal Changes in Vegetation Phenology at Middle and High Latitudes of the Northern Hemisphere over the Past Three Decades[J].Remote Sensing,2015,7(8):10973-10995.

    [6] 孫慶齡,李寶林,許麗麗,等.2000-2013年三江源植被NDVI變化趨勢(shì)及影響因素分析[J].地球信息科學(xué)學(xué)報(bào),2016,18(12):1707-1716.

    [7] 朱明明,侯西勇,毋亭.基于MODIS NDVI的中國(guó)東部植被覆蓋時(shí)空變化特征[J].遙感技術(shù)與應(yīng)用,2013,28(6):1027-1032.

    [8] 毛德華,王宗明,羅玲,等.基于 MODIS 和 AVHRR 數(shù)據(jù)源的東北地區(qū)植被NDVI變化及其與氣溫和降水間的相關(guān)分析[J].遙感技術(shù)與應(yīng)用,2012,27(1):77-85.

    [9] 嚴(yán)恩萍,林輝,黨永峰,等.2000~2012 年京津風(fēng)沙源治理區(qū)植被覆蓋時(shí)空演變特征[J].生態(tài)學(xué)報(bào),2014,34(17):5007-5020.

    [10] 何月,樊高峰,張小偉,等.浙江省植被NDVI動(dòng)態(tài)及其對(duì)氣候的響應(yīng)[J].生態(tài)學(xué)報(bào),2012,32(14):4352-4362.

    [11] 宋富強(qiáng),邢開(kāi)雄,劉陽(yáng),等.基于 MODIS/NDVI 的陜北地區(qū)植被動(dòng)態(tài)監(jiān)測(cè)與評(píng)價(jià)[J].生態(tài)學(xué)報(bào),2011,31(2):354-363.

    [12] 王園香,唐世浩,鄭照軍.1982~2006年中國(guó)5-9月的NDVI變化與人類(lèi)活動(dòng)影響分析[J].地球信息科學(xué)學(xué)報(bào),2015,17(11):1333-1340.

    [13] 劉憲鋒,朱秀芳,潘耀忠,等.1982-2012年中國(guó)植被覆蓋時(shí)空變化特征[J].生態(tài)學(xué)報(bào),2015,35(16):5331-5342.

    [14] SOUTHWORTH J,ZHU L,BUNTING E,et al.Changes in Vegetation Persistence Across Global Savanna Landscapes 1982~2010[J].Journal of Land Use Science,2016,11(1):7-32.

    [15] 符靜,秦建新,張猛,等.2001~2013年洞庭湖流域植被覆蓋度時(shí)空變化特征[J].地球信息科學(xué)學(xué)報(bào),2016,18(9):1209-1216.

    [16] STOW D,DAESCHNER S,HOPE A,et al.Variability of the Seasonally Integrated Normalized Difference Vegetation Index Across the North Slope of Alaska in the 1990s[J].International Journal of Remote Sensing,2003,24(5):1111-1117.

    [17] 徐建華.計(jì)量地理學(xué)[M].北京:高等教育出版社,2006.

    [18] 謝平,陳廣才,雷紅富.基于Hurst系數(shù)的水文變異分析方法[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào),2009,17(1):32-39.

    [19] 江田漢,鄧蓮堂.Hurst 指數(shù)估計(jì)中存在的若干問(wèn)題[J].地理科學(xué),2004,24(2):177-182.

    [20] HURST H E.Long-term Storage Capacity of Reservoirs[J].American Society of Civil Engineering,1951,116:770-808.

    [21] MANDELBROT B B,WALLIS J R.Robustness of the Rescaled Range R/S in the Measurement of Noncyclic Long Run Statistical Dependence[J].Water Resources Research,1969,5(5):967-988.

    [22] 陳云浩,李曉兵,史培軍.1983~1992年中國(guó)陸地NDVI變化的氣候因子驅(qū)動(dòng)分析[J].植物生態(tài)學(xué)報(bào),2001,25(6):716-720.

    [23] 王永財(cái),孫艷玲,王中良.1998-2011 年海河流域植被覆蓋變化及氣候因子驅(qū)動(dòng)分析[J].資源科學(xué),2014,36(3):594-620.

    [24] 劉憲峰,潘耀忠,朱秀芳,等.2000-2014年秦巴山區(qū)植被覆蓋時(shí)空變化特征及其歸因[J].地理學(xué)報(bào),2015,70(5):705-716.

    TheSpatiotemporalChangesofVegetationCoverandClimateFactorDriversofVariationinFuxiRiverBasinFrom1998-2013

    WANGErli1,2,LIUXiaofang1,2*,HOUYani1,2,CHENGuangjian1,WANXin1

    (1.School of Computer Science, Sichuan University of Science & Engineering, Zigong 643000, China; 2.Sichuan Provincial Academician(Expert) Workstation-Integrated Perception and Application of Eco-Environment in River Basins, Zigong 643000, China)

    Vegetation coverage is an important index to evaluate the environment. On the basis of SPOT VEGETATION from 1998 to 2013, distribution, spatial pattern and trends of NDVI were analyzed by the methods of trend analysis, Hurst exponent and partial correlation analysis in Fuxi River Basin. Combined with precipitation and temperature, vegetation cover change driving factors were analyzed using correlation coefficient. The result indicates that vegetation cover shows a rising trend and the growth rate is 0.24%/16a. Hurst exponent shows that the reverse characteristics of vegetation coverage change are stronger than the same characteristic in Fuxi River Basin, and 60.98% of the total area will follow a continuous increasing trend with the slope of NDVI. The proportion of persistent degradation area is 5.97%, mainly distributed in the southern part of Zigong and surrounding towns. Vegetation cover change is mainly driven by non-climatic factors, which accounting for 88.48% of the entire area. The growth of vegetation is more effected by temperature than precipitation in the climatic factors.

    vegetation coverage; Hurst index; partial correlation analysis; driven factor; Fuxi River Basin

    TP79

    A

    2017-07-21

    四川省科技廳項(xiàng)目(2017GZ0303);自貢市科技局項(xiàng)目(2016DZ11);江河流域生態(tài)環(huán)境的集成感知與應(yīng)用院士(專(zhuān)家)工作站項(xiàng)目(2016YSGZZ02);企業(yè)信息化與物聯(lián)網(wǎng)測(cè)控技術(shù)四川省重點(diǎn)實(shí)驗(yàn)室項(xiàng)目(2016WYY01);四川省智慧旅游研究基地項(xiàng)目(ZHY17-01)

    王二麗(1991-),女,河北石家莊人,碩士,主要從事環(huán)境遙感方面的研究,(E-mail)welsuse@163.com

    1673-1549(2017)05-0083-07

    10.11863/j.suse.2017.05.15

    猜你喜歡
    氣候因子溪河氣溫
    基于氣候變量的苧麻產(chǎn)量SSA-BP預(yù)測(cè)模型
    基于FY-3D和FY-4A的氣溫時(shí)空融合
    深冬氣溫多變 蔬菜管理要隨機(jī)應(yīng)變
    安徽裕溪河流域發(fā)現(xiàn)的石制品研究
    東方考古(2020年0期)2020-11-06 05:34:02
    美麗的泗溪河
    美麗的泗溪河
    四川水稻產(chǎn)量及其構(gòu)成要素對(duì)不同生育期氣候因子的響應(yīng)分析
    與氣溫成反比的東西
    天山北坡NDVI對(duì)氣候因子響應(yīng)的敏感性分析
    沙漠地區(qū)微波地表發(fā)射率年內(nèi)變化規(guī)律與氣候因子的關(guān)系分析
    欧美激情 高清一区二区三区| 久久精品国产a三级三级三级| 在线观看免费高清a一片| 国产成人av教育| 90打野战视频偷拍视频| 老汉色av国产亚洲站长工具| 两个人免费观看高清视频| 变态另类成人亚洲欧美熟女 | 久久99一区二区三区| 亚洲精品久久成人aⅴ小说| 久久久精品区二区三区| av国产精品久久久久影院| 可以免费在线观看a视频的电影网站| 中亚洲国语对白在线视频| 色94色欧美一区二区| 久久久欧美国产精品| 欧美黑人精品巨大| av欧美777| 热99久久久久精品小说推荐| 欧美黑人精品巨大| 久久精品aⅴ一区二区三区四区| 久久亚洲精品不卡| 国产亚洲一区二区精品| 啦啦啦视频在线资源免费观看| 国产99久久九九免费精品| 人人澡人人妻人| 高潮久久久久久久久久久不卡| 亚洲 国产 在线| 久久人人爽av亚洲精品天堂| 91麻豆av在线| 欧美日韩视频精品一区| 国产亚洲精品久久久久5区| 久久中文字幕一级| 叶爱在线成人免费视频播放| 露出奶头的视频| 无遮挡黄片免费观看| 国产精品一区二区精品视频观看| 啦啦啦中文免费视频观看日本| 18禁裸乳无遮挡动漫免费视频| 中国美女看黄片| 亚洲av日韩精品久久久久久密| 精品一品国产午夜福利视频| 欧美日韩黄片免| 一区二区三区乱码不卡18| 法律面前人人平等表现在哪些方面| 欧美人与性动交α欧美精品济南到| 欧美日韩亚洲综合一区二区三区_| 国产欧美日韩一区二区三| 成人国产一区最新在线观看| 成年版毛片免费区| 欧美另类亚洲清纯唯美| 无人区码免费观看不卡 | 美国免费a级毛片| 欧美乱码精品一区二区三区| 丰满少妇做爰视频| 日韩免费高清中文字幕av| 超碰成人久久| 亚洲午夜理论影院| 一本色道久久久久久精品综合| 超色免费av| 国产精品美女特级片免费视频播放器 | 亚洲一区二区三区欧美精品| 黄色视频,在线免费观看| 99精品欧美一区二区三区四区| 宅男免费午夜| 亚洲精品粉嫩美女一区| 两个人免费观看高清视频| 大型黄色视频在线免费观看| 国产精品国产高清国产av | 多毛熟女@视频| 久久久精品国产亚洲av高清涩受| av超薄肉色丝袜交足视频| 精品久久久久久久毛片微露脸| 亚洲精品一二三| 在线观看一区二区三区激情| 久久久久久免费高清国产稀缺| 在线十欧美十亚洲十日本专区| 我的亚洲天堂| 午夜久久久在线观看| 视频区图区小说| a在线观看视频网站| 男女免费视频国产| 午夜福利免费观看在线| a级片在线免费高清观看视频| 国产一区有黄有色的免费视频| 在线观看www视频免费| 中文亚洲av片在线观看爽 | 久久国产精品影院| 国产成人欧美| 女性被躁到高潮视频| 欧美激情久久久久久爽电影 | 亚洲精品国产区一区二| 18禁观看日本| 国产一区有黄有色的免费视频| 亚洲国产欧美日韩在线播放| 飞空精品影院首页| 大型黄色视频在线免费观看| 十八禁高潮呻吟视频| 亚洲国产毛片av蜜桃av| 久久精品国产a三级三级三级| 亚洲全国av大片| 亚洲av成人不卡在线观看播放网| 热re99久久精品国产66热6| 天堂俺去俺来也www色官网| 国产欧美日韩一区二区三区在线| 精品亚洲乱码少妇综合久久| 亚洲欧美色中文字幕在线| 蜜桃国产av成人99| 下体分泌物呈黄色| 久久精品国产综合久久久| 国产精品电影一区二区三区 | 国产精品免费大片| 99九九在线精品视频| 精品高清国产在线一区| 男男h啪啪无遮挡| 无限看片的www在线观看| 水蜜桃什么品种好| 成人av一区二区三区在线看| videos熟女内射| 欧美黄色片欧美黄色片| 91av网站免费观看| 天堂中文最新版在线下载| 亚洲一码二码三码区别大吗| 桃红色精品国产亚洲av| 久9热在线精品视频| 在线av久久热| 国产精品欧美亚洲77777| 高清av免费在线| 黑人猛操日本美女一级片| 丰满迷人的少妇在线观看| 久久影院123| 人人妻人人添人人爽欧美一区卜| 我要看黄色一级片免费的| 九色亚洲精品在线播放| 欧美日韩亚洲高清精品| 精品国产亚洲在线| 亚洲欧美日韩另类电影网站| 老熟妇乱子伦视频在线观看| 成人永久免费在线观看视频 | 国产成人系列免费观看| 无遮挡黄片免费观看| 国产男靠女视频免费网站| 伊人久久大香线蕉亚洲五| 动漫黄色视频在线观看| 成人特级黄色片久久久久久久 | 亚洲美女黄片视频| 女人高潮潮喷娇喘18禁视频| 乱人伦中国视频| 91字幕亚洲| 男人舔女人的私密视频| 国产高清视频在线播放一区| 动漫黄色视频在线观看| 午夜精品国产一区二区电影| 国产真人三级小视频在线观看| 午夜福利,免费看| 国产精品一区二区在线观看99| av网站免费在线观看视频| 操出白浆在线播放| 国产精品 欧美亚洲| 亚洲免费av在线视频| 免费观看a级毛片全部| 久久久国产一区二区| 免费在线观看日本一区| 激情在线观看视频在线高清 | 曰老女人黄片| 超碰97精品在线观看| 亚洲欧美色中文字幕在线| 国产高清激情床上av| 深夜精品福利| 性色av乱码一区二区三区2| 久9热在线精品视频| 亚洲欧洲精品一区二区精品久久久| 久久精品国产综合久久久| 亚洲一区中文字幕在线| 9热在线视频观看99| 黄色视频在线播放观看不卡| 国产精品98久久久久久宅男小说| 午夜激情av网站| 欧美 亚洲 国产 日韩一| 欧美国产精品va在线观看不卡| 老司机靠b影院| 丁香欧美五月| 亚洲精品国产精品久久久不卡| 国产精品 国内视频| 久久久精品国产亚洲av高清涩受| 午夜成年电影在线免费观看| 男男h啪啪无遮挡| 免费在线观看日本一区| 老司机亚洲免费影院| 又黄又粗又硬又大视频| 免费日韩欧美在线观看| 在线亚洲精品国产二区图片欧美| 日韩欧美一区视频在线观看| av有码第一页| 国产亚洲精品久久久久5区| a在线观看视频网站| 国产成人免费观看mmmm| 国产成人啪精品午夜网站| 啦啦啦在线免费观看视频4| 水蜜桃什么品种好| 精品亚洲成国产av| 久久久久精品人妻al黑| 丝袜人妻中文字幕| 一区二区av电影网| 国产精品秋霞免费鲁丝片| 久久这里只有精品19| 精品久久久久久电影网| 久久久久国产一级毛片高清牌| 亚洲精品国产一区二区精华液| 黄色丝袜av网址大全| 精品国产乱子伦一区二区三区| 窝窝影院91人妻| 一区二区三区乱码不卡18| 下体分泌物呈黄色| 精品国产超薄肉色丝袜足j| 国产欧美日韩一区二区三区在线| 高清欧美精品videossex| 成人精品一区二区免费| 国产精品久久久久成人av| 欧美日韩av久久| 日本黄色视频三级网站网址 | 欧美日韩av久久| 欧美亚洲日本最大视频资源| 午夜免费成人在线视频| 亚洲精品久久成人aⅴ小说| 黑丝袜美女国产一区| 亚洲黑人精品在线| 操出白浆在线播放| 午夜免费鲁丝| 久久久久网色| 天天躁狠狠躁夜夜躁狠狠躁| 成年人免费黄色播放视频| av超薄肉色丝袜交足视频| 黄色片一级片一级黄色片| 国产人伦9x9x在线观看| 欧美性长视频在线观看| 国产真人三级小视频在线观看| 日韩视频一区二区在线观看| 国产极品粉嫩免费观看在线| 免费一级毛片在线播放高清视频 | 一二三四在线观看免费中文在| 啦啦啦免费观看视频1| 亚洲熟妇熟女久久| 久久精品成人免费网站| 国产激情久久老熟女| 欧美精品啪啪一区二区三区| 午夜福利视频在线观看免费| 国产精品.久久久| 91国产中文字幕| 老司机午夜十八禁免费视频| av一本久久久久| 大片免费播放器 马上看| 激情在线观看视频在线高清 | 黄色视频在线播放观看不卡| 少妇粗大呻吟视频| 亚洲精品在线观看二区| 成人特级黄色片久久久久久久 | 亚洲成人国产一区在线观看| 一进一出好大好爽视频| av欧美777| 久久毛片免费看一区二区三区| 国产真人三级小视频在线观看| 色94色欧美一区二区| 国产福利在线免费观看视频| 后天国语完整版免费观看| www.熟女人妻精品国产| 亚洲五月色婷婷综合| 午夜福利视频在线观看免费| 国产av精品麻豆| e午夜精品久久久久久久| 天天躁日日躁夜夜躁夜夜| 国产在线一区二区三区精| 五月开心婷婷网| 欧美人与性动交α欧美精品济南到| 亚洲精品乱久久久久久| 久久99热这里只频精品6学生| 欧美激情久久久久久爽电影 | 色在线成人网| 亚洲精品一二三| 亚洲一码二码三码区别大吗| 中文字幕色久视频| 看免费av毛片| 欧美成人午夜精品| 80岁老熟妇乱子伦牲交| 亚洲一码二码三码区别大吗| 成年女人毛片免费观看观看9 | 亚洲黑人精品在线| 免费看十八禁软件| 香蕉国产在线看| 久久ye,这里只有精品| 色播在线永久视频| 成人特级黄色片久久久久久久 | 亚洲精品美女久久久久99蜜臀| 久久99一区二区三区| 日本一区二区免费在线视频| 一级毛片电影观看| 女人爽到高潮嗷嗷叫在线视频| 999久久久精品免费观看国产| 午夜福利欧美成人| 亚洲av电影在线进入| 超色免费av| 午夜激情久久久久久久| 亚洲成人手机| 一区二区三区激情视频| 欧美在线黄色| 亚洲av日韩精品久久久久久密| 这个男人来自地球电影免费观看| 成人国产av品久久久| 亚洲精品国产区一区二| 日韩精品免费视频一区二区三区| 午夜福利免费观看在线| bbb黄色大片| 亚洲色图综合在线观看| 午夜两性在线视频| 久久精品aⅴ一区二区三区四区| 汤姆久久久久久久影院中文字幕| 天堂8中文在线网| 国产亚洲午夜精品一区二区久久| 午夜91福利影院| 国产一区有黄有色的免费视频| 国产伦人伦偷精品视频| 大片电影免费在线观看免费| 黄色视频在线播放观看不卡| 91麻豆精品激情在线观看国产 | 最新美女视频免费是黄的| 九色亚洲精品在线播放| 日韩欧美一区二区三区在线观看 | 亚洲 国产 在线| 亚洲中文av在线| 亚洲国产av新网站| 一区二区三区乱码不卡18| 成人av一区二区三区在线看| 啦啦啦 在线观看视频| 69精品国产乱码久久久| 啦啦啦在线免费观看视频4| 少妇粗大呻吟视频| 国产一区二区三区在线臀色熟女 | 国产精品亚洲一级av第二区| 日韩 欧美 亚洲 中文字幕| 日韩视频一区二区在线观看| 午夜精品久久久久久毛片777| av天堂久久9| 亚洲精品一卡2卡三卡4卡5卡| 国产日韩欧美在线精品| 99精国产麻豆久久婷婷| 精品国产乱码久久久久久小说| 欧美久久黑人一区二区| 侵犯人妻中文字幕一二三四区| 午夜福利欧美成人| 18禁美女被吸乳视频| 国产成人一区二区三区免费视频网站| 50天的宝宝边吃奶边哭怎么回事| 在线观看免费视频日本深夜| 一区二区日韩欧美中文字幕| 男女免费视频国产| 可以免费在线观看a视频的电影网站| netflix在线观看网站| 亚洲专区中文字幕在线| 久久久久精品人妻al黑| 成人手机av| 久久热在线av| 国产亚洲av高清不卡| 人妻 亚洲 视频| 国产精品久久久久久人妻精品电影 | 最新美女视频免费是黄的| 两性午夜刺激爽爽歪歪视频在线观看 | 99国产精品一区二区三区| 一区福利在线观看| 最新的欧美精品一区二区| 日日摸夜夜添夜夜添小说| 美女午夜性视频免费| 青草久久国产| 国产亚洲精品一区二区www | 国产精品1区2区在线观看. | 91九色精品人成在线观看| 在线观看一区二区三区激情| 亚洲av国产av综合av卡| 黄色怎么调成土黄色| 女人被躁到高潮嗷嗷叫费观| 日韩三级视频一区二区三区| 成年女人毛片免费观看观看9 | 欧美日韩亚洲国产一区二区在线观看 | 国产三级黄色录像| 日韩中文字幕视频在线看片| 亚洲熟女毛片儿| 人人妻人人澡人人爽人人夜夜| 亚洲伊人久久精品综合| 久久精品国产综合久久久| videos熟女内射| 久久人妻熟女aⅴ| 精品午夜福利视频在线观看一区 | 91字幕亚洲| 午夜福利在线免费观看网站| e午夜精品久久久久久久| 欧美激情久久久久久爽电影 | 亚洲国产毛片av蜜桃av| 免费高清在线观看日韩| 91麻豆av在线| 国产无遮挡羞羞视频在线观看| 极品教师在线免费播放| 男女高潮啪啪啪动态图| 18禁观看日本| 久久亚洲真实| 窝窝影院91人妻| 久久狼人影院| 在线观看免费午夜福利视频| 免费日韩欧美在线观看| 91精品三级在线观看| 99riav亚洲国产免费| 亚洲欧洲日产国产| 一区二区三区激情视频| 国产精品自产拍在线观看55亚洲 | 日韩免费av在线播放| 国产不卡av网站在线观看| 亚洲五月婷婷丁香| 精品少妇内射三级| 巨乳人妻的诱惑在线观看| 成人手机av| 国产亚洲精品一区二区www | 精品亚洲成a人片在线观看| 香蕉丝袜av| 大片电影免费在线观看免费| 黄色成人免费大全| 亚洲视频免费观看视频| 中文字幕人妻熟女乱码| 久久久水蜜桃国产精品网| 一边摸一边抽搐一进一小说 | 亚洲精品av麻豆狂野| 国产又色又爽无遮挡免费看| 精品国产乱码久久久久久小说| 老熟妇乱子伦视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲成人国产一区在线观看| 视频区图区小说| 国产免费视频播放在线视频| av视频免费观看在线观看| 电影成人av| 久久午夜综合久久蜜桃| 欧美 日韩 精品 国产| 久久中文字幕人妻熟女| 国产精品免费大片| 另类精品久久| 90打野战视频偷拍视频| 自拍欧美九色日韩亚洲蝌蚪91| 一级a爱视频在线免费观看| 99re6热这里在线精品视频| 精品一区二区三区四区五区乱码| 免费日韩欧美在线观看| 亚洲精品一二三| 人人妻人人澡人人看| 在线 av 中文字幕| 日韩欧美三级三区| 精品久久久久久电影网| 黄色视频在线播放观看不卡| 亚洲精品美女久久久久99蜜臀| 国产成人免费观看mmmm| 久热爱精品视频在线9| 国产在线免费精品| 欧美成人免费av一区二区三区 | 老司机午夜福利在线观看视频 | 国产又色又爽无遮挡免费看| 久久性视频一级片| 久久精品国产综合久久久| 精品人妻在线不人妻| 女人精品久久久久毛片| 999久久久精品免费观看国产| 99国产综合亚洲精品| 正在播放国产对白刺激| 免费看a级黄色片| 香蕉国产在线看| 国产精品久久久久久精品电影小说| 91成年电影在线观看| 精品第一国产精品| 男女无遮挡免费网站观看| 久久这里只有精品19| 天堂俺去俺来也www色官网| 最近最新中文字幕大全免费视频| 亚洲中文日韩欧美视频| 99re6热这里在线精品视频| 可以免费在线观看a视频的电影网站| 麻豆成人av在线观看| 久久人妻福利社区极品人妻图片| av国产精品久久久久影院| 精品人妻在线不人妻| 桃花免费在线播放| 老熟妇仑乱视频hdxx| 久久精品成人免费网站| 成人永久免费在线观看视频 | 精品午夜福利视频在线观看一区 | 精品人妻熟女毛片av久久网站| 成人国语在线视频| 999精品在线视频| cao死你这个sao货| 欧美乱码精品一区二区三区| 国产欧美日韩精品亚洲av| 法律面前人人平等表现在哪些方面| 女人高潮潮喷娇喘18禁视频| 嫁个100分男人电影在线观看| 国产成+人综合+亚洲专区| 人妻一区二区av| 日韩中文字幕视频在线看片| 丝袜在线中文字幕| 亚洲黑人精品在线| 精品卡一卡二卡四卡免费| 久久天堂一区二区三区四区| aaaaa片日本免费| 另类精品久久| 天天操日日干夜夜撸| 精品高清国产在线一区| a级片在线免费高清观看视频| 十分钟在线观看高清视频www| 国产无遮挡羞羞视频在线观看| 国产成人精品无人区| 精品少妇内射三级| 黄片小视频在线播放| 老司机福利观看| 欧美在线黄色| 熟女少妇亚洲综合色aaa.| 精品国产一区二区三区四区第35| av又黄又爽大尺度在线免费看| 人人妻,人人澡人人爽秒播| 桃红色精品国产亚洲av| 久久精品亚洲熟妇少妇任你| 五月开心婷婷网| 9色porny在线观看| 日韩 欧美 亚洲 中文字幕| 亚洲国产av新网站| 老汉色∧v一级毛片| 中国美女看黄片| 五月天丁香电影| 露出奶头的视频| 欧美久久黑人一区二区| av视频免费观看在线观看| 美女扒开内裤让男人捅视频| 啦啦啦免费观看视频1| 777米奇影视久久| 日韩成人在线观看一区二区三区| 日韩一卡2卡3卡4卡2021年| www.自偷自拍.com| 精品一区二区三区四区五区乱码| 久久久久久人人人人人| 最黄视频免费看| 纯流量卡能插随身wifi吗| 欧美日韩福利视频一区二区| 欧美 日韩 精品 国产| 美女主播在线视频| 亚洲国产成人一精品久久久| 制服诱惑二区| 欧美成人午夜精品| 最新的欧美精品一区二区| 水蜜桃什么品种好| 人成视频在线观看免费观看| 99国产精品一区二区三区| 俄罗斯特黄特色一大片| 在线观看舔阴道视频| 国产精品自产拍在线观看55亚洲 | 国产男女内射视频| 国产精品久久久av美女十八| 在线av久久热| 午夜福利视频在线观看免费| 久久久精品国产亚洲av高清涩受| 黄色 视频免费看| 国产黄频视频在线观看| 国产日韩欧美视频二区| 女同久久另类99精品国产91| 男女下面插进去视频免费观看| 国产精品熟女久久久久浪| 国产免费福利视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 久久人妻福利社区极品人妻图片| 亚洲av片天天在线观看| 五月开心婷婷网| 99精品在免费线老司机午夜| 99精品欧美一区二区三区四区| 久久午夜综合久久蜜桃| 国产欧美日韩一区二区三| √禁漫天堂资源中文www| 国产淫语在线视频| av网站在线播放免费| 丁香欧美五月| 大码成人一级视频| 天天影视国产精品| 大码成人一级视频| 啦啦啦免费观看视频1| 亚洲第一欧美日韩一区二区三区 | 99re6热这里在线精品视频| 国产一区二区 视频在线| 欧美激情 高清一区二区三区| 国产精品一区二区免费欧美| 欧美老熟妇乱子伦牲交| 两个人看的免费小视频| 一进一出好大好爽视频| 欧美久久黑人一区二区| 岛国毛片在线播放| 色综合婷婷激情| 国产有黄有色有爽视频| 热99国产精品久久久久久7| 啪啪无遮挡十八禁网站| 亚洲五月色婷婷综合| 久久九九热精品免费| 三上悠亚av全集在线观看| 亚洲人成电影免费在线| 久久精品国产a三级三级三级| 最新在线观看一区二区三区| 国产精品久久久av美女十八| svipshipincom国产片| 亚洲av美国av| 亚洲精品av麻豆狂野| videosex国产| 香蕉久久夜色| 久久人妻熟女aⅴ| av免费在线观看网站| 免费人妻精品一区二区三区视频| 黄色 视频免费看| 亚洲人成伊人成综合网2020|