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

    基于擾動(dòng)指數(shù)的天水南部植被擾動(dòng)時(shí)空序列研究

    2019-04-02 02:02:10周妍妍郭曉娟鄒明亮陳冠光郭建軍岳東霞
    生態(tài)科學(xué) 2019年1期
    關(guān)鍵詞:天水擾動(dòng)植被

    周妍妍, 郭曉娟, 鄒明亮, 陳冠光, 郭建軍, 岳東霞*

    ?

    基于擾動(dòng)指數(shù)的天水南部植被擾動(dòng)時(shí)空序列研究

    周妍妍1, 郭曉娟1, 鄒明亮1, 陳冠光1, 郭建軍2, 岳東霞1*

    1. 蘭州大學(xué)資源環(huán)境學(xué)院, 蘭州 730000 2. 中國(guó)科學(xué)院西北生態(tài)環(huán)境資源研究院沙漠與沙漠化重點(diǎn)實(shí)驗(yàn)室, 蘭州 730000

    利用全球擾動(dòng)指數(shù)(The MODIS Global Disturbance Index, MGDI)模型并以長(zhǎng)時(shí)間序列的MODIS的增強(qiáng)型植被指數(shù)(Enhanced Vegetation Index, EVI)和地表溫度(Land Surface Temperature, LST)數(shù)據(jù)作為數(shù)據(jù)源對(duì)天水南部地區(qū)2001—2016年植被擾動(dòng)指數(shù)時(shí)空分布變化特征和驅(qū)動(dòng)因子進(jìn)行研究。結(jié)果表明: (1)16年來(lái)植被擾動(dòng)指數(shù)的均值空間分布不均, 其值分布于0.66—1.09之間, 天水市南部地區(qū)的西北和西部部分地區(qū)植被擾動(dòng)指數(shù)明顯低于其它地區(qū), 且主要為積極擾動(dòng)。(2)植被擾動(dòng)指數(shù)的年際變化呈降低趨勢(shì), 占總面積的89.58%, 其中研究區(qū)中農(nóng)耕區(qū)減小的趨勢(shì)最為明顯。植被擾動(dòng)指數(shù)的變異系數(shù)分析結(jié)果表明, 2001—2016年以來(lái)天水南部大部分地區(qū)植被擾動(dòng)指數(shù)變化較穩(wěn)定。其中, 西北部的農(nóng)耕區(qū)擾動(dòng)指數(shù)波動(dòng)幅度較大, 而林區(qū)波動(dòng)相對(duì)較小。(3)植被擾動(dòng)情況受海拔和坡度、河流和道路的影響??臻g分布上, 擾動(dòng)主要發(fā)生在海拔相對(duì)較高、坡度較緩及近河流和道路地區(qū)。

    MODIS; 植被; 擾動(dòng)指數(shù); 時(shí)空變化; 天水南部

    0 前言

    植被作為陸地生態(tài)系統(tǒng)必不可少的組成部分, 能夠根據(jù)陸地生態(tài)系統(tǒng)中水、氣、熱等的狀況, 調(diào)控其內(nèi)部與外部的全球能量物質(zhì)平衡、生物化學(xué)循環(huán)和水循環(huán)[1]。它對(duì)區(qū)域生態(tài)系統(tǒng)環(huán)境變化有著重要指示作用, 氣候變化、人類活動(dòng), 以及大氣中CO2的施肥效應(yīng)都會(huì)對(duì)植被產(chǎn)生擾動(dòng)[2]。生態(tài)擾動(dòng)是一個(gè)導(dǎo)致生態(tài)系統(tǒng)結(jié)構(gòu)和功能持續(xù)中斷的事件, 其中植被擾動(dòng)是指生態(tài)系統(tǒng)葉面積指數(shù)(Leaf Area Index, LAI)在一年以上發(fā)生重大變化的事件, 擾動(dòng)可分為積極擾動(dòng)(如農(nóng)業(yè)灌溉)和消極擾動(dòng)(如火災(zāi))[3–6]。陸地表層植被擾動(dòng)的研究有助于跟蹤生物圈對(duì)氣候變化的響應(yīng), 其時(shí)空特征更是區(qū)域植被生物量和碳儲(chǔ)量研究的關(guān)鍵。獲取高精度植被擾動(dòng)信息不僅有利于區(qū)域乃至全球生物量的精確估算[7], 而且可以提高區(qū)域或者全球碳匯的估算精度, 將在國(guó)家政策的制定、碳預(yù)算的評(píng)估以及改善資源管理中發(fā)揮重要的作用。隨著全球生態(tài)環(huán)境的不斷惡化, 量化擾動(dòng)對(duì)于了解生態(tài)系統(tǒng)的動(dòng)態(tài)至關(guān)重要[8], 因此量化擾動(dòng)已成為全球生態(tài)系統(tǒng)環(huán)境的研究熱點(diǎn)。近半個(gè)世紀(jì)以來(lái), 不斷發(fā)展的遙感(Remote Sensing, RS)和地理信息系統(tǒng)(Geographic Information System, GIS)技術(shù)為大面積長(zhǎng)時(shí)間序列植被擾動(dòng)的研究提供可能[9]。MODIS等遙感影像可提供數(shù)據(jù)來(lái)提取關(guān)鍵物候參數(shù), 并監(jiān)測(cè)植被動(dòng)態(tài)的趨勢(shì)[10]。2007年, Mildrexler等[11]根據(jù)植被指數(shù)和地表溫度的負(fù)相關(guān)關(guān)系提出了擾動(dòng)指數(shù), 用于監(jiān)測(cè)基于像元尺度LST/EVI的長(zhǎng)期變化; 2009年, Mildrexler等[12]改進(jìn)了先前提出的擾動(dòng)指數(shù), 將其更名為MODIS全球擾動(dòng)指數(shù)(The MODIS Global Disturbance Index, MGDI), 監(jiān)測(cè)北美森林火災(zāi)發(fā)生的位置、嚴(yán)重程度以及颶風(fēng)和砍伐森林對(duì)植被造成的破壞。由于MGDI使用的是年最大合成數(shù)據(jù), 這就避免了地表溫度可能存在的自然變化, 因此MGDI在國(guó)內(nèi)外植被擾動(dòng)研究中得到廣泛的應(yīng)用。然而, 國(guó)內(nèi)外學(xué)者主要用MGDI模型監(jiān)測(cè)短時(shí)間內(nèi)的植被變化。如COOPS等[13]研究表明, MGDI監(jiān)測(cè)得到的擾動(dòng)區(qū)域與使用其它衛(wèi)星數(shù)據(jù)獲取的林火區(qū)域具有很高的相關(guān)性, 并且指出MGDI還可用于區(qū)域尺度植被病蟲害的監(jiān)測(cè); Sun等[14]應(yīng)用MGDI監(jiān)測(cè)汶川地震森林?jǐn)_動(dòng), 研究表明MGDI能夠較好的應(yīng)用于地震擾動(dòng)的監(jiān)測(cè); Dutta等[15]對(duì)比分析印度和孟加拉國(guó)紅樹林旋風(fēng)前后的擾動(dòng), 計(jì)算出了擾動(dòng)的面積; 尤慧等[16]基于MODIS數(shù)據(jù)利用MGDI監(jiān)測(cè)了加拿大地區(qū)森林火燒跡地的面積; 此外, Phompila等[17]應(yīng)用擾動(dòng)指數(shù)檢測(cè)老撾熱帶森林土地覆被空間變化, 研究表明, 植被擾動(dòng)指數(shù)能夠高精度的監(jiān)測(cè)7年內(nèi)的植被變化, 但缺少分析說(shuō)明其變化的驅(qū)動(dòng)因素。天水南部地區(qū)是一個(gè)典型的農(nóng)林牧交錯(cuò)復(fù)合區(qū), 屬于我國(guó)氣候變化的敏感區(qū)和生態(tài)環(huán)境的脆弱區(qū), 易受人類活動(dòng)和滑坡、泥石流等自然災(zāi)害的擾動(dòng)[18]。本文基于MGDI模型, 采用MODIS EVI數(shù)據(jù)和LST數(shù)據(jù)監(jiān)測(cè)天水南部的植被擾動(dòng)狀況, 分析了研究區(qū)長(zhǎng)時(shí)間序列擾動(dòng)指數(shù)的時(shí)空動(dòng)態(tài)分布特征和驅(qū)動(dòng)因素, 全面客觀地認(rèn)識(shí)了植被擾動(dòng)變化的內(nèi)在聯(lián)系和作用機(jī)理, 對(duì)減緩研究區(qū)植被破壞和植被保護(hù)提供了科學(xué)依據(jù)[19], 建立了區(qū)域生態(tài)環(huán)境與社會(huì)經(jīng)濟(jì)協(xié)調(diào)的發(fā)展模式。

    1 研究區(qū)概況

    天水市南部地區(qū)位于黃土高原西南部, 介于東經(jīng)105°42′—106°42′, 北緯34°4′—34°37′之間, 地處陜、甘兩省交界, 包括黨川鄉(xiāng)、利橋鄉(xiāng)、東岔鎮(zhèn)、李子鄉(xiāng)、娘娘壩等12個(gè)鄉(xiāng)鎮(zhèn)。天水南部地區(qū)總面積約為3262 km2, 海拔主要分布在1000—2000 m之間, 占天水市總面積的22.77%。天水市南部地區(qū)林草占總面積的85.72%, 南部小隴山林區(qū), 樹種類型多樣, 森林資源十分豐富。農(nóng)耕區(qū)主要分布于研究區(qū)的西北部地區(qū), 占總面積的12.70%。其屬隴中南部冷溫帶半濕潤(rùn)氣候, 年平均降水量574 mm, 時(shí)空分布不均勻, 自西北向東南逐漸增加, 中東部山區(qū)雨量在600 mm以上, 且降雨多集中在4月下旬至10月上旬, 冬半年降水較少[20]。這一地區(qū)土壤蓄水性能很差, 每當(dāng)遇有強(qiáng)降水, 雨水就會(huì)迅速在溝壑處匯集成洪流, 沖刷表土, 沖毀坡地, 引發(fā)山體滑坡、泥石流等災(zāi)害, 造成交通和通訊中段、房屋倒塌、作物受災(zāi), 甚至人員傷亡等[21]。2013年7月25日天水市發(fā)生群發(fā)性山洪地質(zhì)災(zāi)害, 其中本研究區(qū)所在的天水市麥積山區(qū)發(fā)生200多個(gè)地質(zhì)災(zāi)害[22]。

    圖1 研究區(qū)概況

    Figure 1 Sketch map of the southern part of Tianshui

    2 研究方法

    2.1 HANTS

    時(shí)間序列諧波分析法(Harmonic Analysis of Time Series, HANTS)是改進(jìn)傅立葉分析的一種算法[23]。長(zhǎng)時(shí)間序列數(shù)據(jù)被解為幾個(gè)最為顯著的頻率, 利用最小二乘法擬合各頻率的正弦或余弦分量, 其中通過不停的迭代得到各頻率的振幅和相位。該方法是平滑和濾波兩種方法的結(jié)合, 它通過利用遙感影像的時(shí)空特點(diǎn), 將其時(shí)間上的變化規(guī)律和空間上的分布規(guī)律聯(lián)系起來(lái), 其主要目的是遙感圖像的去云處理[24], 具體內(nèi)容可參考文獻(xiàn)[25–26]。本文針對(duì)研究區(qū)特殊的氣候地形環(huán)境, 在ENVI軟件下利用HANTS方法對(duì)MODIS EVI和MODIS LST數(shù)據(jù)進(jìn)行無(wú)云處理, 以此提高數(shù)據(jù)質(zhì)量, 獲得整個(gè)時(shí)序數(shù)據(jù)長(zhǎng)期變化趨勢(shì), 以為擾動(dòng)指數(shù)的計(jì)算奠定基礎(chǔ)。

    2.2 擾動(dòng)指數(shù)

    擾動(dòng)指數(shù)為每年合成的最大地表溫度和同年合成的最大植被指數(shù)的比率與目標(biāo)監(jiān)測(cè)年之前多年的擾動(dòng)比率均值的比值[27]。MODIS全球擾動(dòng)指數(shù)(MGDI)是使用植被指數(shù)和地表溫度數(shù)據(jù)來(lái)檢測(cè)生態(tài)擾動(dòng)的一種方法(Mildrexler et al., 2009), MGDI的提出是基于通過潛熱能量傳遞, 植被密度降低而表面溫度增加的事實(shí)[28]。此外, 兩個(gè)輸入變量, 即地表溫度(LST)和增強(qiáng)植被指數(shù)(EVI), 反應(yīng)不同的生物物理過程, 從而豐富指標(biāo)的內(nèi)容[29]。該指數(shù)在長(zhǎng)時(shí)間序列數(shù)據(jù)集的基礎(chǔ)上, 計(jì)算像元尺度長(zhǎng)期變化, 并且檢測(cè)表面能量的顯著年際變化。生態(tài)系統(tǒng)在未發(fā)生擾動(dòng)時(shí), 在正常年份的地表能量平衡都會(huì)在自然變異范圍內(nèi)波動(dòng)(擾動(dòng)指數(shù)一般接近1.0)。當(dāng)生態(tài)系統(tǒng)受人類灌溉活動(dòng)擾動(dòng)時(shí), 擾動(dòng)指數(shù)就會(huì)小于1, 且擾動(dòng)強(qiáng)度越大, 其值越小; 當(dāng)生態(tài)系統(tǒng)受到干旱, 火災(zāi)等自然災(zāi)害擾動(dòng)時(shí), 擾動(dòng)指數(shù)大于1, 且擾動(dòng)程度越高, 其值越大。擾動(dòng)指數(shù)根據(jù)擾動(dòng)方式可分為瞬時(shí)擾動(dòng)指數(shù)(如火災(zāi))和非瞬時(shí)擾動(dòng)指數(shù)(如蟲災(zāi), 過度放牧, 氣候變化和干旱)[30]。根據(jù)擾動(dòng)的影響可分為積極擾動(dòng)(灌溉)和消極擾動(dòng)(砍伐、火災(zāi))。在本研究中, 使用MGDI方法分析了2001—2016年天水南部地區(qū)的因人類活動(dòng)和滑坡泥石流等自然災(zāi)害造成的非瞬時(shí)擾動(dòng)。非瞬時(shí)MGDI表達(dá)式如下[31]:

    公式中, DI為擾動(dòng)指數(shù); EVImax為每個(gè)像元的EVI一年內(nèi)最大值; LSTmax為每個(gè)像元的LST一年內(nèi)最大值; 分子表示目標(biāo)年份的(LSTmax/EVImax), 分母表示除目標(biāo)年份以外之前多年(LSTmax/EVImax)的均值。

    2.3 線性趨勢(shì)分析

    為研究2001—2016年擾動(dòng)指數(shù)變化的空間分布特征采用一元線性回歸方程方法, 分析每個(gè)像元的DI變化趨勢(shì), 定量計(jì)算研究區(qū)范圍內(nèi)中所有像元的DI與時(shí)間的回歸斜率。斜率K即為DI線性趨勢(shì), 其值為正說(shuō)明擾動(dòng)指數(shù)有增加趨勢(shì), 為負(fù)說(shuō)明擾動(dòng)指數(shù)呈減少趨勢(shì), K絕對(duì)值越大, 則說(shuō)明擾動(dòng)指數(shù)變化越明顯[32]。趨勢(shì)斜率用最小二乘計(jì)算, 公式如下[33]:

    公式中, K為擾動(dòng)指數(shù)變化的線性趨勢(shì);為研究年數(shù); DI為每個(gè)像元在第年第個(gè)像元的DI值。

    2.4 穩(wěn)定性分析

    為了解天水南部地區(qū)DI在空間上的波動(dòng)程度, 計(jì)算2001—2016年每個(gè)柵格DI的變異系數(shù)(The coefficient of variation, C)。變異系數(shù)是反映數(shù)據(jù)變異程度的一個(gè)統(tǒng)計(jì)量, 可消除因單位和平均值不同對(duì)變異程度比較的影響[34], 用其來(lái)評(píng)估DI數(shù)據(jù)在時(shí)間與空間上的相對(duì)變化程度, 變異系數(shù)越大, 變異程度越大; 反之, 變異程度越小[35]。計(jì)算式為[36]:

    公式中, C為變異系數(shù),`為2001—2016年DI平均值,s為2001—2016年DI標(biāo)準(zhǔn)差。C值越大, 表示數(shù)據(jù)分布越離散, 數(shù)據(jù)時(shí)間序列波動(dòng)較大; 反之則表示數(shù)據(jù)分布較聚集, 數(shù)據(jù)時(shí)間序列波動(dòng)較小, 穩(wěn)定性較好。

    3 數(shù)據(jù)源與數(shù)據(jù)處理

    (1)MODIS EVI數(shù)據(jù)和LST數(shù)據(jù): 搭載在Terra和Aqua兩顆衛(wèi)星上的中分辨率成像光譜儀(Moderate-Resolution Imaging Spectroradiometer, MODIS)是美國(guó)地球觀測(cè)系統(tǒng)(EOS)計(jì)劃中用于觀測(cè)全球生物和物理過程的重要儀器。它獲取的影像具有36個(gè)中等分辨率水平(0.25 um—1 um)的光譜波段, 每1—2 d對(duì)地球表面觀測(cè)一次。獲取陸地和海洋溫度、初級(jí)生產(chǎn)率、陸地表面覆蓋、云、汽溶膠、水汽和火情等目標(biāo)的圖像。NASA提供的MODIS全球數(shù)據(jù)產(chǎn)品, 共有44種標(biāo)準(zhǔn)產(chǎn)品, 具有不同的時(shí)間和空間分辨率。其中MOD13Q1產(chǎn)品的增強(qiáng)型植被指數(shù)EVI數(shù)據(jù), 是空間分辨率為250 m的16 d合成產(chǎn)品; MOD11A2產(chǎn)品的地表溫度LST數(shù)據(jù), 是空間分辨率為1000 m的8 d合成產(chǎn)品。利用MRT軟件把所有MODIS遙感影像投影為相同的坐標(biāo)系, 通過最鄰近法重采樣為250 m空間分辨率。本文選取MODIS數(shù)據(jù)產(chǎn)品行列號(hào)為h26v05的2000—2016年所有的EVI和LST數(shù)據(jù)產(chǎn)品用于擾動(dòng)指數(shù)的計(jì)算。利用研究區(qū)每年的EVI和LST數(shù)據(jù), 在ArcGIS軟件的像元統(tǒng)計(jì)數(shù)據(jù)工具中, 計(jì)算出每個(gè)像元一年內(nèi)的最大值, 取該值為目標(biāo)年的EVImax, 每年合成一幅當(dāng)年EVImax影像。根據(jù)同樣的算法, 每年合成一幅當(dāng)年LSTmax。

    (2)數(shù)字高程模型(DEM): 來(lái)源于中國(guó)地理空間數(shù)據(jù)云的數(shù)字高程數(shù)據(jù)產(chǎn)品, 其為90 m空間分辨率。根據(jù)DEM數(shù)據(jù)利用ArcGIS軟件提取研究區(qū)的河流和坡度。

    (3)土地利用類型數(shù)據(jù): 2015年的土地利用數(shù)據(jù)來(lái)源于中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http:// www.resdc.cn), 其空間分辨率為30 m[37]。

    4 結(jié)果與分析

    4.1 天水南部地區(qū)2001—2016年擾動(dòng)指數(shù)空間分布特征

    根據(jù)公式(1), 在ArcGIS軟件支持下, 將空間化的各參數(shù)帶入MGDI模型進(jìn)行柵格地圖代數(shù)運(yùn)算, 得到2001—2016年DI時(shí)間序列數(shù)據(jù)集, 制作2001— 2016年平均值的天水南部擾動(dòng)指數(shù)圖(圖3a)。根據(jù)構(gòu)建的DI時(shí)間序列數(shù)據(jù)集, 計(jì)算每個(gè)像元點(diǎn)多年數(shù)據(jù)集的平均值和標(biāo)準(zhǔn)差, 確定每個(gè)像元點(diǎn)的平均狀態(tài)和自然變異范圍, 從而制作2016年擾動(dòng)情況圖[38](圖3b)。由圖3可知,擾動(dòng)指數(shù)的空間分布不均一, 其值分布于0.66—1.09之間, 絕大部分的值小于1, 且為擾動(dòng)較小區(qū)域。從區(qū)域分布看, 西北及北部邊緣區(qū)域、西部娘娘壩鎮(zhèn)的一部分及北部邊緣地區(qū)擾動(dòng)指數(shù)明顯低于其它地區(qū), 主要是積極擾動(dòng), 占總面積的31.36%, 擾動(dòng)面積最大。由于這些地區(qū)植被類型主要為耕地, 耕地由于受人類灌溉, 耕種等農(nóng)業(yè)活動(dòng)影響較大, 導(dǎo)致地表溫度較低, 植被覆蓋較好, 擾動(dòng)指數(shù)較低[39–40]。河谷地區(qū), 地勢(shì)平坦, 土壤肥沃, 水資源豐富, 土地利用類型主要為耕地, 擾動(dòng)指數(shù)也普遍低于1。除河谷地區(qū)其它大部分地區(qū)主要為林地, 受自然和人為擾動(dòng)較少, 擾動(dòng)指數(shù)接近于1, 主要為無(wú)擾動(dòng)區(qū)。消極擾動(dòng)的區(qū)域主要分布于西北部伯陽(yáng)鄉(xiāng)及西部的娘娘壩鎮(zhèn)和李子鄉(xiāng), 擾動(dòng)面積較小, 主要是位于滑坡泥石流災(zāi)害的多發(fā)區(qū)。

    圖2 2015年天水南部土地覆被狀況

    Figure 2 Land cover of the southern part of Tianshui in 2015

    圖3 2001—2016年天水南部擾動(dòng)指數(shù)和2016年擾動(dòng)情況空間分布圖

    Figure 3 Spatial distribution map of average DI from 2001 to 2016 and DI in 2016 over the southern part of Tianshui

    4.2 天水南部地區(qū)擾動(dòng)指數(shù)的時(shí)空變化特征

    4.2.1 擾動(dòng)指數(shù)年際變化趨勢(shì)

    根據(jù)公式(2), 用16年逐一像元的DI值與年份計(jì)算線性趨勢(shì), DI回歸斜率能夠反映植被DI變化的方向和速率。圖4為2001—2016年天水南部擾動(dòng)指數(shù)的變化趨勢(shì), 擾動(dòng)指數(shù)變化較小, 平均趨勢(shì)系數(shù)為-0.0038, 呈緩慢減小的趨勢(shì)。從面積比例上看, 擾動(dòng)指數(shù)減小的區(qū)域占總面積的89.58%, 明顯降低的區(qū)域主要分布于研究區(qū)西北部的伯陽(yáng)鄉(xiāng)和元龍鄉(xiāng)及娘娘壩鎮(zhèn)的西北部, 這主要由于人類的灌溉耕作等農(nóng)業(yè)活動(dòng)使植被覆蓋相對(duì)升高, 地表溫度降低[41]。中部大面積的林區(qū)擾動(dòng)指數(shù)呈增加趨勢(shì), 占總面積10.42%, 且DI逐漸接近于1。說(shuō)明該區(qū)域受自然和人為擾動(dòng)減少, 植被覆蓋度逐漸增加, 植被狀況不斷改善。這主要由于自1999年以來(lái), 該區(qū)域積極實(shí)施天然林保護(hù)、封山育林、退耕還林等工程的實(shí)施, 該區(qū)域森林覆蓋率呈顯著上升趨勢(shì), 且林地面積不斷擴(kuò)大[42]。DI多年基本不變的區(qū)域主要分布于現(xiàn)代冰川上, 該區(qū)海拔高且遠(yuǎn)離城區(qū), 受自然和人為擾動(dòng)較少[43]。總之, 研究區(qū)DI增加的區(qū)域遠(yuǎn)高于減少區(qū)域的面積, 同時(shí)增加的強(qiáng)度也高于減少的強(qiáng)度。

    圖4 2001—2016年天水南部擾動(dòng)指數(shù)變化趨勢(shì)空間分布圖

    Figure 4 Changing trends of DI in the southern part of Tianshui from 2001 to 2016

    4.2.2 擾動(dòng)指數(shù)變化的穩(wěn)定性分析

    變異系數(shù)能夠反映植被DI的時(shí)空變動(dòng)程度。圖5顯示的是天水南部地區(qū)16年DI逐像元計(jì)算出的變異系數(shù)。研究區(qū)的變異系數(shù)分布0.03—0.38之間, 平均值為0.08, 整體DI年際變異幅度較小, 變化較為穩(wěn)定。DI的空間分布與其年際變異系數(shù)的空間分布存在高度的負(fù)相關(guān)性, 區(qū)域不同地區(qū)的DI穩(wěn)定性存在明顯的空間差異。DI高波動(dòng)區(qū)多分布在西北部伯陽(yáng)鄉(xiāng), 元龍鄉(xiāng)等鄉(xiāng)鎮(zhèn)及白家河和紅崖河等河流的兩側(cè), 且是擾動(dòng)指數(shù)低值區(qū), 該區(qū)域受到人類生產(chǎn)生活的影響[44], DI的穩(wěn)定程度低。草地的生長(zhǎng)受水熱條件以及人為干預(yù)的影響較大, 故近16年草地DI表現(xiàn)出較大的波動(dòng)性。林地主要分布在丘陵和山區(qū), 自然環(huán)境好且對(duì)氣候和人為影響具有較強(qiáng)的適應(yīng)性, 生長(zhǎng)穩(wěn)定, DI波動(dòng)較小[45]。

    4.3 天水南部地區(qū)植被擾動(dòng)驅(qū)動(dòng)因子分析

    4.3.1 海拔和坡度

    地形是影響植被分布及擾動(dòng)的重要因素, 本文按照400 m等間隔, 將海拔由低到高分為五個(gè)等級(jí), 按照10°等間隔, 將坡度由低到高劃分為五個(gè)等級(jí)[46], 利用ArcGIS10.2軟件中面積制表工具統(tǒng)計(jì)在不同等級(jí)海拔、坡度中擾動(dòng)區(qū)域的分布情況[47]。由表1可知, 研究區(qū)植被擾動(dòng)區(qū)域的海拔分別主要分布在等級(jí)三和等級(jí)四, 分別占擾動(dòng)面積的37.76%、50.68%, 由此可知, 較高海拔即1200—2000 m范圍內(nèi)的植被易受擾動(dòng)。由于研究區(qū)中高海拔區(qū)占研究區(qū)面積最大, 且是人類主要活動(dòng)的區(qū)域, 隨著海拔升高, 植被增加, 人類活動(dòng)較少, 擾動(dòng)面積少[48]。植被擾動(dòng)區(qū)域的坡度分別主要分布于等級(jí)二和等級(jí)三, 分別占擾動(dòng)面積的41.19%、30.08%, 由此可知坡度較緩的植被更加易受到擾動(dòng)影響, 并且在天水地區(qū)滑坡災(zāi)害多發(fā)生在10°—40°之間, 增加該區(qū)域的擾動(dòng)[49], 結(jié)果說(shuō)明地形是影響植被擾動(dòng)的重要因素。

    4.3.2 距河流和道路的距離

    近河流和道路的區(qū)域是工農(nóng)業(yè)生產(chǎn)主要區(qū)域, 其影響著植被的蓋度和擾動(dòng)狀況。利用ArcGIS平臺(tái)的緩沖區(qū)分析工具, 按1 km間隔(緩沖區(qū)半徑)建立水體、道路邊界緩沖區(qū), 進(jìn)行區(qū)域分析, 得到各個(gè)緩沖區(qū)等級(jí)擾動(dòng)制表。由表2可知, 在靠近道路和水系的緩沖區(qū)一和緩沖區(qū)二區(qū)域, 植被擾動(dòng)發(fā)生的所占比例較高[50]。緩沖區(qū)之間植被擾動(dòng)的比較顯示, 無(wú)論是河流還是道路附近, 植被擾動(dòng)區(qū)域都由內(nèi)向外顯著遞減。河流和道路附近植被擾動(dòng)都較大, 而不同范圍內(nèi)的變化有所差別, 在一定程度上反映出人為活動(dòng)與植被擾動(dòng)存在重要的本質(zhì)聯(lián)系[51]??傊? 在空間分布上, 天水南部植被擾動(dòng)主要發(fā)生在海拔相對(duì)較高、坡度較緩及近河流和道路區(qū)域。

    圖5 2001—2016年天水南部擾動(dòng)指數(shù)變化穩(wěn)定程度

    Figure 5 Standard deviation of DI changes in the southern part of Tianshui from 2001 to 2016

    表1 天水南部地區(qū)植被擾動(dòng)面積在不同海拔和坡度的分布

    表2 天水南部植被擾動(dòng)面積在道路和水系不同緩沖區(qū)的分布

    5 結(jié)論

    本文利用MODIS EVI和MODIS LST遙感數(shù)據(jù), 基于MGDI模型監(jiān)測(cè)了天水南部地區(qū)2001—2016年植被擾動(dòng)狀況, 分析擾動(dòng)指數(shù)的時(shí)空變化特征及其對(duì)自然和人為因素的響應(yīng), 得出以下結(jié)論:

    1)16年來(lái)擾動(dòng)指數(shù)的均值空間分布不均, 其值分布于0.66—1.09之間, 研究區(qū)西北和西部部分地區(qū)擾動(dòng)指數(shù)明顯低于其它地區(qū); 2)擾動(dòng)指數(shù)的年際變化呈降低趨勢(shì), 占總面積的89.58%, 其中研究區(qū)中農(nóng)耕區(qū)減小的趨勢(shì)最為明顯。擾動(dòng)指數(shù)的變異系數(shù)分析結(jié)果表明, 2001—2016年以來(lái)天水南部大部分地區(qū)擾動(dòng)指數(shù)變化較穩(wěn)定。其中, 西北農(nóng)耕區(qū)擾動(dòng)指數(shù)波動(dòng)較大, 而林區(qū)波動(dòng)較小; 3)植被擾動(dòng)受海拔、坡度與河流和道路距離因素的影響, 空間上, 擾動(dòng)主要發(fā)生在海拔相對(duì)較高、坡度較緩和近水體和道路地區(qū)。植被擾動(dòng)是自然和人為因素綜合作用的結(jié)果, 其中合理的安排農(nóng)業(yè)生產(chǎn)活動(dòng)以及退耕還林、還草等政策, 有力促進(jìn)了綠化工作的開展, 大大改善了生態(tài)與環(huán)境。此外, 植被的擾動(dòng)變化與氣候變化、土地利用方式、城市化和人口密度等因素密切相關(guān), 這將是下一步研究的重點(diǎn)[52]。

    [1] SCHIMEL D S. Terrestrial biogeochemical cycles: Global estimates with remote sensing[J]. Remote Sensing of Environment, 1995, 51(1): 49–56.

    [2] PIAO S, MOHAMMAT A, FANG J, et al. NDVI-based increase in growth of temperate grasslands and its responses to climate changes in China[J]. Global Environmental Change, 2006, 16(4): 340–348.

    [3] PICKETT S T A,WHITE P S. The ecology of natural disturbance and patch dynamics[M]. Academic Press, 1985.

    [4] TILMAN D. The Resource- Ratio hypothesis of plant succession[J]. American Naturalist, 1985, 125(6): 827–852.

    [5] WARING R H, RUNNING S W. Forest ecosystems: analysis at multiple scales[J]. Forest Ecosystems Analysis at Multiple Scales, 2007.

    [6] DALE V H, BROWN S, HAEUBER R A, et al. Ecological principles and guidelines for managing the use of land[J]. Ecological Applications, 2014, 10(3): 639–670.

    [7] YANG C, SHEN R. Progress in the study of forest disturbance by remote sensing[J]. Remote Sensing for Land & Resources, 2015, 27(1): 1–8.

    [8] BADRELDIN N, GOOSSENS R. A satellite–based disturbance index algorithm for monitoring mitigation strategies effects on desertification change in an arid environment[J]. Mitigation and Adaptation Strategies for Global Change, 2015, 20(2): 263–276.

    [9] 鄧朝平, 郭鈮, 王介民, 等. 近20余年來(lái)西北地區(qū)植被變化特征分析[J]. 冰川凍土, 2006, 28(5): 686–693.

    [10] PéREZHOYOS A, MARTíNEZ B, GILABERT M A, et al. A multi-temporal analysis of vegetation dynamics in the Iberian peninsula using MODIS-NDVI data[J]. Earsel Eproceedings, 2010(1).

    [11] MILDREXLER D J. A new satellite-based methodology for continental-scale disturbance detection[J]. Ecological Applications A Publication of the Ecological Society of America, 2007, 17(1): 235–50.

    [12] MILDREXLER D J, ZHAO Maosheng, et al. Testing a MODIS global disturbance index across North America[J]. Remote sensing of environment, 2009, 113(10): 2103– 3566.

    [13] COOPS N C, WULDER M A, IWANICKA D. Large area monitoring with a MODIS–based disturbance index(DI) sensitive to annual and seasonal variations[J]. Remote Sensing of Environment, 2009, 113(6): 1250–1261.

    [14] SUN L, CHEN J, LI T. A MODIS- based method for detecting large-scale vegetation disturbance due to natural hazards: a case study of Wenchuan earthquake stricken regions in China[J]. Stochastic Environmental Research & Risk Assessment, 2015, 30(8): 1–12.

    [15] DUTTA D, DAS P K, PAUl S, et al. Assessment of ecological disturbance in the mangrove forest of Sundar-bans caused by cyclones using MODIS time–series data(2001—2011)[J]. Natural Hazards, 2015, 79(2): 775– 790.

    [16] 尤慧, 劉榮高 , 祝善友, 等. 加拿大北方森林火燒跡地遙感分析[J]. 地球信息科學(xué), 2013, (4): 597–603.

    [17] PHOMPILA C. Applying the Global disturbance index for detecting vegetation changes in Lao Tropical Forests[J]. Advances in Remote Sensing, 2015, 05(1): 1–9.

    [18] 胡利平, 姚延鋒, 裴古娥, 等. 天水地區(qū)近50年氣溫與降水變化特征[J]. 地理科學(xué)進(jìn)展, 2009, 28(4): 651–656.

    [19] 李月臣, 宮鵬, 劉春霞, 等. 北方13省1982年~1999年植被變化及其與氣候因子的關(guān)系[J]. 資源科學(xué), 2006, 28(2): 109–117.

    [20] 郭慧, 黃濤, 鄧茂芝, 等. 甘肅天水地區(qū)45 a來(lái)強(qiáng)降水與洪澇災(zāi)害特征分析[J]. 冰川凍土, 2007, 29(5): 808–812.

    [21] 趙洪濤. 甘肅中南部滑坡、泥石流災(zāi)害機(jī)制及防災(zāi)可持續(xù)發(fā)展研究[D]. 蘭州: 蘭州大學(xué), 2010.

    [22] 于國(guó)強(qiáng), 張茂省, 胡煒. 天水市“7.25”群發(fā)性山洪地質(zhì)災(zāi)害發(fā)育特征及成因分析[J]. 西北地質(zhì), 2014(3): 185– 191.

    [23] 趙偉, 李召良. 利用MODIS/EVI時(shí)間序列數(shù)據(jù)分析干旱對(duì)植被的影響[J]. 地理科學(xué)進(jìn)展, 2007, 26(6): 40–47.

    [24] 劉惠峰. 基于時(shí)序NDVI的疏勒河流域植被覆蓋分類研究[D]. 蘭州: 蘭州大學(xué), 2014.

    [25] VERHOEF W, MENENTI M, AZZAH S. A colour composite of NOAA—AVHRR—NDVI based on time series analysis 1981— 11992. International Journal of Remote Sensing, 1996, 17: 231–235.

    [26] ROERINK G J, MENENTI M, VERHOEF W. Reconstructing cloudfree NDVI composites using Fourier analysis of time series[J]. International Journal of Remote Sensing, 2000, 21(9): 1911–1917.

    [27] 尤慧. 近十年?yáng)|南亞森林?jǐn)_動(dòng)的遙感分析[D]. 南京: 南京信息工程大學(xué), 2013.

    [28] RAMAKRISHNA R N.Estimation of Regional Surface Resistance to Evapotranspiration from NDVI and Therma–IR AVHRR Data[J]. Journal of Applied Meteo-rology, 1989, 28(4): 276–284.

    [29] LAMBIN E F, EHRLICH D. The surface temperature- vegetation index space for land cover and land–cover change analysis[J]. International Journal of Remote Sensing, 1996, 17(3): 463–487.

    [30] GOETZ S J. Multi–sensor analysis of NDVI, surface temperature and biophysical variables at a mixed grassland site[J]. International Journal of Remote Sensing, 1997, 18(1): 71–94.

    [31] 曹云. 基于MODIS數(shù)據(jù)產(chǎn)品與Landsat遙感影像的云南省近十年森林?jǐn)_動(dòng)監(jiān)測(cè)方法研究[D]. 南京: 南京信息工程大學(xué), 2015.

    [32] 于泉洲, 梁春玲, 劉煜杰, 等. 基于MODIS的山東省植被覆蓋時(shí)空變化及其原因分析[J]. 生態(tài)環(huán)境學(xué)報(bào), 2015, 24(11): 1799–1807.

    [33] MA M, FRANK V. Interannual variability of vegetation cover in the Chinese Heihe River Basin and its relation to meterological parameters[J]. International Journal of Remote Sensing, 2006, 27(16): 3473–3486.

    [34] 張仁平. 2000—2012年中國(guó)北方草地NDVI和氣候因子時(shí)空變化[J]. 中國(guó)沙漠, 2015, 35(5): 1403–1412.

    [35] 楊慧. 空間分析與建模[M]. 北京: 清華大學(xué)出版社, 2013.

    [36] 潘竟虎. 2001—2010年疏勒河流域植被凈初級(jí)生產(chǎn)力時(shí)空變化及其與氣候因子的關(guān)系[J]. 生態(tài)學(xué)報(bào), 2017, 37(6): 1888–1899.

    [37] 徐新良. 土地利用/覆被變化時(shí)空信息分析方法及應(yīng)用[M]. 北京: 科學(xué)技術(shù)文獻(xiàn)出版社, 2014.

    [38] 吳立葉. 區(qū)域森林?jǐn)_動(dòng)遙感監(jiān)測(cè)研究[D]. 南京: 南京信息工程大學(xué), 2014.

    [39] FRIEDL M A, DAVIS F W. Sources of variation in radiometric surface temperature over a tallgrass prairie[M]// Handbook of Africa's international relations/. Routledge, 1994: 1–17.

    [40] BICHERON P, DEFOURNY P, BROCKMANN C. GLOBCOVER: products description and validation report[J]. Foro Mundial De La Salud, 2011, 17(3): 285–287.

    [41] 趙晶, 周耀治, 鄧興耀. 西北干旱區(qū)植被覆蓋時(shí)空動(dòng)態(tài)特征[J]. 林業(yè)資源管理, 2017, 2(1): 118–126.

    [42] BADRELDIN N, GOOSSENS R. A satellite–based disturbance index algorithm for monitoring mitigation strategies effects on desertification change in an arid environment[J]. Mitigation and Adaptation Strategies for Global Change, 2015, 20(2): 263–276.

    [43] 任沖, 鞠洪波, 張懷清, 等. 天水市近30年林地動(dòng)態(tài)變化遙感監(jiān)測(cè)研究[J]. 林業(yè)科學(xué)研究, 2017, 30(1): 25–33.

    [44] CHANG Yaping, DING Yongjian, ZHAO Qiudong, et al. Remote estimation of terrestrial evapotranspiration by Landsat 5 TM and the SEBAL model in cold and high–altitude regions: a case study of the upper reach of the Shule River Basin, China[J]. Hydrological Processes, 2016: n/a-n/a.

    [45] 張東海, 任志遠(yuǎn), 王曉峰, 等. 基于MODIS的陜西黃土高原植被覆蓋度變化特征及其驅(qū)動(dòng)分析[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 2013, 29(1): 29–35.

    [46] 奚硯濤, 徐勇, 劉欣婷. 2000-2013年江蘇省不同植被類型NDVI時(shí)空變化特征[J]. 水土保持研究, 2016, 23(1): 86-91.

    [47] 楊辰. 區(qū)域森林植被擾動(dòng)遙感監(jiān)測(cè)研究[D]. 南京: 南京信息工程大學(xué), 2013.

    [48] 張宇欣. 基于擾動(dòng)指數(shù)的岷江上游森林?jǐn)_動(dòng)時(shí)間序列研究——以理縣為例[J]. 浙江林業(yè)科技, 2017, 37(3): 48–53.

    [49] 齊敬輝, 牛叔文, 馬利邦, 等. 2000-2014年疏勒河流域植被覆蓋時(shí)空變化[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 2016, 32(5): 757–766.

    [50] 趙芳. 三江源區(qū)草地MODIS植被指數(shù)時(shí)空變異及驅(qū)動(dòng)因子分析[D]. 西寧: 青海大學(xué), 2012.

    [51] 李小亞, 張勃, 靳自寶. 基于MODIS-NDVI的甘肅河?xùn)|地區(qū)植被覆蓋度動(dòng)態(tài)監(jiān)測(cè)[J]. 水土保持研究, 2013, 20(1): 112–115.

    [52] 王文輝. 福建長(zhǎng)汀植被覆蓋度變化的主要驅(qū)動(dòng)影響因子及影響力分析[J]. 福建農(nóng)林大學(xué)學(xué)報(bào)(自然版), 2017, 46(3): 277–283.

    Study on spatio-temporal variation of vegetation disturbance in the southern part of Tianshui by disturbance index

    ZHOU Yanyan1, GUO Xiaojuan1, ZOU Mingliang1, CHEN Guanguang1, GUO Jianjun2, YUE Dongxia1*

    1. College of Earth and Environmental Sciences, Lanzhou University, Lanzhou 730000, China 2. Key Laboratory of Desert and Desertification, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China

    The study investigated the temporal-spatial distribution characteristics and driving factors of vegetation disturbance index in the southern part of Tianshui from 2001 to 2016 based on MGDI model, using the Enhanced Vegetation Index(EVI) and Land Surface Temperature(LST) in long time series as data sources. The results are as follows. (1)The spatial variation range of mean vegetation disturbance index value was 0.66-1.09mm on the grid during 2001-2016, which had a significant difference in spatial patterns. And the vegetation disturbance index of the northwest and western parts of the study area was obviously lower than that of other areas, and mainly for positive disturbance. (2)The interannual variability of vegetation disturbance index was decreasing, accounting for 89.58% of the total area. Among all the study area, the reducing trend was most obvious in the cultivated area. The analysis of variance coefficient of vegetation disturbance index showed that the change of vegetation disturbance index in most parts of southern Tianshui since 2001-2016 was relatively stable while the disturbance index of the farming area in the northwest was larger, and the forest area fluctuated less. (3)Vegetation disturbances were affected mainly by altitude, slope, river, and road. Judging from spatial distribution, the disturbance occurred mainly at relatively low elevations, with slower slopes and near rivers and roads.

    MODIS; vegetation; disturbance index; spatial and temporal variation; south of Tianshui

    10.14108/j.cnki.1008-8873.2019.01.014

    P9

    A

    1008-8873(2019)01-102-09

    2017-12-07;

    2018-3-20

    甘肅省民生科技計(jì)劃(1503FCME006)

    周妍妍(1992—), 女, 河北邯鄲人, 在讀碩士研究生, 主要從事生態(tài)安全評(píng)價(jià)研究, E-mail:yyzhou16@lzu.edu.cn

    岳東霞, 女, 博士, 教授, 主要從事生態(tài)安全評(píng)價(jià)研究, E-mail: dxyue@lzu.edu.cn

    周妍妍, 郭曉娟, 鄒明亮, 等. 基于擾動(dòng)指數(shù)的天水南部植被擾動(dòng)時(shí)空序列研究[J]. 生態(tài)科學(xué), 2019, 38(1): 102-110.

    ZHOU Yanyan, GUO Xiaojuan, ZOU Mingliang, et al. Study on spatio-temporal variation of vegetation disturbance in the southern part of Tianshui by disturbance index[J]. Ecological Science, 2019, 38(1): 102-110.

    猜你喜歡
    天水擾動(dòng)植被
    基于植被復(fù)綠技術(shù)的孔植試驗(yàn)及應(yīng)用
    Bernoulli泛函上典則酉對(duì)合的擾動(dòng)
    天水嬸與兩岸商貿(mào)
    天水地區(qū)的『秦與戎』
    (h)性質(zhì)及其擾動(dòng)
    重返絲綢之路—從天水到青海湖
    美食(2018年10期)2018-10-18 08:10:58
    綠色植被在溯溪旅游中的應(yīng)用
    《天水之鏡像》
    小噪聲擾動(dòng)的二維擴(kuò)散的極大似然估計(jì)
    基于原生植被的長(zhǎng)山群島植被退化分析
    成人特级黄色片久久久久久久| 老司机在亚洲福利影院| 国产精品久久视频播放| 超碰av人人做人人爽久久 | 一区二区三区国产精品乱码| 黄色丝袜av网址大全| tocl精华| 国产精品,欧美在线| 婷婷精品国产亚洲av| 人妻丰满熟妇av一区二区三区| 午夜福利在线观看吧| 最近最新中文字幕大全免费视频| 国产黄a三级三级三级人| 国产精品,欧美在线| 亚洲美女视频黄频| 亚洲精品日韩av片在线观看 | 19禁男女啪啪无遮挡网站| 成人特级av手机在线观看| 女生性感内裤真人,穿戴方法视频| 哪里可以看免费的av片| 国产一区二区亚洲精品在线观看| 亚洲18禁久久av| 亚洲内射少妇av| 亚洲精品美女久久久久99蜜臀| 国产精品香港三级国产av潘金莲| 久久人人精品亚洲av| 亚洲精品一区av在线观看| 国产三级中文精品| 午夜免费成人在线视频| ponron亚洲| 精品国产超薄肉色丝袜足j| 亚洲国产高清在线一区二区三| 一级黄色大片毛片| 亚洲专区国产一区二区| 久久香蕉精品热| 亚洲 欧美 日韩 在线 免费| 三级国产精品欧美在线观看| 日本黄色视频三级网站网址| 亚洲中文字幕日韩| 人妻夜夜爽99麻豆av| 在线观看av片永久免费下载| 成人国产一区最新在线观看| 亚洲精品在线观看二区| 亚洲精品亚洲一区二区| 国产av在哪里看| 1000部很黄的大片| 人人妻,人人澡人人爽秒播| 久久久久久人人人人人| 精品国产超薄肉色丝袜足j| 久久久久久人人人人人| 欧美又色又爽又黄视频| 亚洲真实伦在线观看| 亚洲成人中文字幕在线播放| 日韩成人在线观看一区二区三区| 国产野战对白在线观看| 国产99白浆流出| 动漫黄色视频在线观看| 成人特级黄色片久久久久久久| 我的老师免费观看完整版| 嫩草影视91久久| 一区二区三区高清视频在线| 精品国产超薄肉色丝袜足j| 国产高潮美女av| 国产精品综合久久久久久久免费| 国产乱人视频| 午夜福利18| 色在线成人网| 啦啦啦免费观看视频1| 精华霜和精华液先用哪个| 丰满的人妻完整版| 欧美黄色片欧美黄色片| 免费高清视频大片| 国产av不卡久久| 亚洲狠狠婷婷综合久久图片| 在线观看美女被高潮喷水网站 | 成人三级黄色视频| 精品国产美女av久久久久小说| 亚洲av美国av| 欧美性猛交黑人性爽| 欧美性猛交╳xxx乱大交人| 亚洲av一区综合| 亚洲av熟女| 51午夜福利影视在线观看| 制服丝袜大香蕉在线| 久久香蕉国产精品| 久久精品国产清高在天天线| 精品乱码久久久久久99久播| 亚洲av免费在线观看| 精品乱码久久久久久99久播| 19禁男女啪啪无遮挡网站| 欧美一区二区国产精品久久精品| 听说在线观看完整版免费高清| 亚洲七黄色美女视频| 一区二区三区免费毛片| 男人和女人高潮做爰伦理| 成人性生交大片免费视频hd| 日本撒尿小便嘘嘘汇集6| 老司机午夜福利在线观看视频| 内地一区二区视频在线| 一本精品99久久精品77| 欧美一级a爱片免费观看看| 亚洲欧美日韩高清在线视频| 91在线观看av| 在线播放无遮挡| 少妇熟女aⅴ在线视频| 国产亚洲精品久久久久久毛片| 神马国产精品三级电影在线观看| 精品不卡国产一区二区三区| 91字幕亚洲| 男女下面进入的视频免费午夜| 欧美乱码精品一区二区三区| 麻豆一二三区av精品| 9191精品国产免费久久| 国产黄片美女视频| 欧美三级亚洲精品| 精品电影一区二区在线| 国产真实伦视频高清在线观看 | 亚洲精品乱码久久久v下载方式 | 久久中文看片网| 女生性感内裤真人,穿戴方法视频| 亚洲av一区综合| 国产精华一区二区三区| 精品无人区乱码1区二区| 国产在视频线在精品| 人人妻人人看人人澡| 岛国在线观看网站| 成年女人看的毛片在线观看| 999久久久精品免费观看国产| 国产三级中文精品| 国产不卡一卡二| 美女高潮喷水抽搐中文字幕| 国产精品爽爽va在线观看网站| 亚洲aⅴ乱码一区二区在线播放| 啪啪无遮挡十八禁网站| 欧美乱妇无乱码| 免费观看的影片在线观看| 久久久色成人| x7x7x7水蜜桃| 两个人看的免费小视频| 成人欧美大片| 精品福利观看| 久久香蕉精品热| 亚洲精品一区av在线观看| 草草在线视频免费看| 午夜福利免费观看在线| 欧美成人a在线观看| 操出白浆在线播放| 老司机福利观看| 久久久久久久久中文| 亚洲色图av天堂| 波多野结衣高清无吗| 午夜福利欧美成人| 国产av不卡久久| 午夜老司机福利剧场| 少妇的丰满在线观看| 夜夜躁狠狠躁天天躁| 无限看片的www在线观看| 身体一侧抽搐| 在线观看一区二区三区| 中文字幕熟女人妻在线| www.www免费av| 亚洲欧美激情综合另类| 欧美中文综合在线视频| 精品久久久久久久人妻蜜臀av| 亚洲无线观看免费| 香蕉久久夜色| 欧美一区二区精品小视频在线| 午夜免费成人在线视频| 国产精品一区二区三区四区久久| 欧美成人a在线观看| 丰满乱子伦码专区| 在线观看午夜福利视频| 熟女人妻精品中文字幕| 亚洲国产中文字幕在线视频| 国产激情欧美一区二区| 白带黄色成豆腐渣| 午夜影院日韩av| 老汉色av国产亚洲站长工具| 久久6这里有精品| 非洲黑人性xxxx精品又粗又长| 精品日产1卡2卡| 国产成人av激情在线播放| 久久精品国产综合久久久| 999久久久精品免费观看国产| 中文字幕人成人乱码亚洲影| or卡值多少钱| 午夜老司机福利剧场| 亚洲欧美精品综合久久99| 免费观看精品视频网站| 搞女人的毛片| 搞女人的毛片| 一区二区三区高清视频在线| 日本黄大片高清| 午夜a级毛片| 19禁男女啪啪无遮挡网站| 午夜福利在线观看吧| 91麻豆av在线| 9191精品国产免费久久| 午夜福利在线观看免费完整高清在 | h日本视频在线播放| 男人舔奶头视频| 亚洲精品色激情综合| 日本 欧美在线| 狂野欧美激情性xxxx| 91在线观看av| 国产免费一级a男人的天堂| av欧美777| 舔av片在线| 一本综合久久免费| 中文字幕人妻熟人妻熟丝袜美 | 90打野战视频偷拍视频| 欧美乱妇无乱码| 最新在线观看一区二区三区| 一级毛片高清免费大全| 日本成人三级电影网站| 国产精品久久久人人做人人爽| 岛国视频午夜一区免费看| 亚洲欧美日韩卡通动漫| 久久久国产成人精品二区| 18+在线观看网站| 日韩免费av在线播放| 在线免费观看的www视频| 国产精品 国内视频| 日本撒尿小便嘘嘘汇集6| 少妇的逼水好多| 1000部很黄的大片| 一进一出好大好爽视频| 啪啪无遮挡十八禁网站| 欧美中文日本在线观看视频| 色综合婷婷激情| 精品人妻偷拍中文字幕| 亚洲精品美女久久久久99蜜臀| 18禁裸乳无遮挡免费网站照片| 欧洲精品卡2卡3卡4卡5卡区| 国产精品亚洲美女久久久| 中文字幕av成人在线电影| 麻豆国产97在线/欧美| 久久午夜亚洲精品久久| 亚洲天堂国产精品一区在线| 亚洲色图av天堂| 一a级毛片在线观看| 亚洲aⅴ乱码一区二区在线播放| 首页视频小说图片口味搜索| 在线国产一区二区在线| 日本免费一区二区三区高清不卡| 欧美黑人欧美精品刺激| 午夜福利在线观看吧| 精品久久久久久成人av| 手机成人av网站| 搡老岳熟女国产| 日本免费a在线| 欧美3d第一页| 国产伦一二天堂av在线观看| 免费看a级黄色片| 国产一区二区三区在线臀色熟女| av女优亚洲男人天堂| 欧美最黄视频在线播放免费| 亚洲成人免费电影在线观看| 久久99热这里只有精品18| 69人妻影院| 老司机午夜福利在线观看视频| 无限看片的www在线观看| 女人高潮潮喷娇喘18禁视频| 欧美日韩中文字幕国产精品一区二区三区| 国产欧美日韩一区二区精品| 国产色爽女视频免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲欧美日韩无卡精品| 亚洲狠狠婷婷综合久久图片| 国产成人av教育| 国产av麻豆久久久久久久| 久久欧美精品欧美久久欧美| 久久精品影院6| 毛片女人毛片| 两个人看的免费小视频| 亚洲成人免费电影在线观看| 久久中文看片网| 国产精品嫩草影院av在线观看 | 精品不卡国产一区二区三区| 欧美av亚洲av综合av国产av| 男女之事视频高清在线观看| 亚洲一区高清亚洲精品| 国产精品99久久久久久久久| 日韩成人在线观看一区二区三区| 精品免费久久久久久久清纯| 最新在线观看一区二区三区| 国产成+人综合+亚洲专区| 成人三级黄色视频| 日韩欧美免费精品| 亚洲欧美精品综合久久99| 手机成人av网站| 他把我摸到了高潮在线观看| 男插女下体视频免费在线播放| 久久人妻av系列| 国产一区二区三区视频了| 九色成人免费人妻av| 日韩有码中文字幕| 亚洲成av人片在线播放无| 亚洲美女黄片视频| 51午夜福利影视在线观看| 久久草成人影院| 老司机在亚洲福利影院| 国产蜜桃级精品一区二区三区| 久久天躁狠狠躁夜夜2o2o| 亚洲国产精品久久男人天堂| 精品一区二区三区视频在线观看免费| 国产国拍精品亚洲av在线观看 | 精品国内亚洲2022精品成人| 中文字幕av成人在线电影| 精品久久久久久久毛片微露脸| 国产69精品久久久久777片| 欧美日本视频| 免费看a级黄色片| 亚洲av第一区精品v没综合| 成熟少妇高潮喷水视频| 国产精品98久久久久久宅男小说| 国产精品久久久久久精品电影| 亚洲欧美一区二区三区黑人| 亚洲不卡免费看| 一个人看视频在线观看www免费 | 国产v大片淫在线免费观看| 欧美日韩精品网址| 亚洲av成人av| av片东京热男人的天堂| 色老头精品视频在线观看| 香蕉av资源在线| 真人一进一出gif抽搐免费| 在线国产一区二区在线| 五月伊人婷婷丁香| 一区二区三区免费毛片| 亚洲av成人精品一区久久| 午夜福利在线在线| 在线观看美女被高潮喷水网站 | 亚洲五月天丁香| 国产高清三级在线| 精品人妻偷拍中文字幕| 日韩欧美免费精品| 国产精品永久免费网站| 悠悠久久av| 亚洲美女视频黄频| 青草久久国产| 国产伦精品一区二区三区四那| 亚洲精品一区av在线观看| 日本a在线网址| 每晚都被弄得嗷嗷叫到高潮| 亚洲av美国av| 女警被强在线播放| 久久精品亚洲精品国产色婷小说| 亚洲成人久久性| 99国产综合亚洲精品| 亚洲人成网站在线播| 九九在线视频观看精品| 国产黄a三级三级三级人| x7x7x7水蜜桃| 国产精品野战在线观看| av国产免费在线观看| 亚洲精品影视一区二区三区av| 国产高清三级在线| 国产精品综合久久久久久久免费| 日本免费a在线| 国产高清videossex| 精品人妻偷拍中文字幕| 一个人免费在线观看电影| 男人舔奶头视频| 国产亚洲欧美在线一区二区| 非洲黑人性xxxx精品又粗又长| 国产av在哪里看| 久久精品综合一区二区三区| 国产成人aa在线观看| 99久久成人亚洲精品观看| 又紧又爽又黄一区二区| 国产欧美日韩精品亚洲av| 国产精品亚洲av一区麻豆| 又紧又爽又黄一区二区| 亚洲无线在线观看| 婷婷丁香在线五月| 国产精品亚洲一级av第二区| 亚洲国产欧美人成| 欧美一级毛片孕妇| 国产亚洲精品久久久久久毛片| 亚洲精品国产精品久久久不卡| 久久精品国产99精品国产亚洲性色| 成年女人看的毛片在线观看| or卡值多少钱| 国产av一区在线观看免费| 亚洲专区中文字幕在线| 亚洲七黄色美女视频| 国产精品一区二区免费欧美| 欧美日本亚洲视频在线播放| 亚洲精品影视一区二区三区av| 成人一区二区视频在线观看| 91av网一区二区| 好男人在线观看高清免费视频| 亚洲国产精品合色在线| 亚洲熟妇熟女久久| 91麻豆av在线| 每晚都被弄得嗷嗷叫到高潮| 色播亚洲综合网| 在线国产一区二区在线| 热99在线观看视频| 国产美女午夜福利| 成人高潮视频无遮挡免费网站| 波多野结衣高清作品| 欧美成人免费av一区二区三区| 精品一区二区三区人妻视频| 国产高清有码在线观看视频| 免费人成在线观看视频色| 五月伊人婷婷丁香| 热99re8久久精品国产| 在线a可以看的网站| 亚洲性夜色夜夜综合| 高清在线国产一区| ponron亚洲| 久久久久久久精品吃奶| 色尼玛亚洲综合影院| 老熟妇乱子伦视频在线观看| 麻豆国产97在线/欧美| 在线观看日韩欧美| 脱女人内裤的视频| 久久99热这里只有精品18| 国产美女午夜福利| 久久国产精品影院| 18禁在线播放成人免费| 免费观看精品视频网站| 国产一区二区三区视频了| 最好的美女福利视频网| 亚洲国产高清在线一区二区三| 国产精品亚洲美女久久久| 观看免费一级毛片| 亚洲欧美精品综合久久99| 日韩亚洲欧美综合| 一夜夜www| 首页视频小说图片口味搜索| 小说图片视频综合网站| 欧美性感艳星| 精品久久久久久久久久免费视频| 亚洲美女黄片视频| 99在线视频只有这里精品首页| 老司机在亚洲福利影院| 日韩欧美免费精品| 99久久精品一区二区三区| www.色视频.com| 国产v大片淫在线免费观看| 亚洲精品在线观看二区| 可以在线观看的亚洲视频| 欧美三级亚洲精品| 亚洲欧美一区二区三区黑人| 精品久久久久久久久久免费视频| 久久精品国产自在天天线| 少妇的丰满在线观看| 日本免费一区二区三区高清不卡| 国产精品久久视频播放| 国产探花极品一区二区| 亚洲成人免费电影在线观看| 欧美激情在线99| 国产又黄又爽又无遮挡在线| 成年版毛片免费区| 国产一区二区激情短视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 免费看十八禁软件| 好看av亚洲va欧美ⅴa在| 岛国视频午夜一区免费看| 法律面前人人平等表现在哪些方面| 人人妻,人人澡人人爽秒播| 午夜精品在线福利| 伊人久久精品亚洲午夜| 乱人视频在线观看| 欧美最黄视频在线播放免费| 91麻豆精品激情在线观看国产| 亚洲av不卡在线观看| 此物有八面人人有两片| 欧美成狂野欧美在线观看| 在线a可以看的网站| 最近最新中文字幕大全免费视频| 免费在线观看日本一区| 18禁国产床啪视频网站| 性欧美人与动物交配| 嫩草影视91久久| 久久九九热精品免费| 国产视频一区二区在线看| 欧美在线一区亚洲| 亚洲午夜理论影院| 国产av不卡久久| 欧美大码av| 久久婷婷人人爽人人干人人爱| 久久草成人影院| 亚洲乱码一区二区免费版| 国产成人av教育| 国产精品综合久久久久久久免费| 毛片女人毛片| 国产精品1区2区在线观看.| 两个人视频免费观看高清| 久久精品国产综合久久久| 美女黄网站色视频| 欧美中文综合在线视频| 美女免费视频网站| 欧美色欧美亚洲另类二区| 中文字幕av成人在线电影| 欧美成狂野欧美在线观看| 全区人妻精品视频| 亚洲精华国产精华精| 国产乱人伦免费视频| 日日干狠狠操夜夜爽| 91九色精品人成在线观看| 亚洲精品美女久久久久99蜜臀| 国产在线精品亚洲第一网站| 热99在线观看视频| 色av中文字幕| 91av网一区二区| 桃色一区二区三区在线观看| 99久久久亚洲精品蜜臀av| 日本一二三区视频观看| 久久国产精品影院| 天天躁日日操中文字幕| 成年版毛片免费区| 两个人视频免费观看高清| 久久精品国产亚洲av涩爱 | 成年女人看的毛片在线观看| 久久久久国内视频| 岛国在线观看网站| 91在线精品国自产拍蜜月 | 听说在线观看完整版免费高清| 亚洲av中文字字幕乱码综合| 亚洲熟妇熟女久久| 国产色婷婷99| 国产精品,欧美在线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品1区2区在线观看.| 99国产精品一区二区蜜桃av| 99久久九九国产精品国产免费| 最近最新中文字幕大全免费视频| 免费看光身美女| 欧美一区二区亚洲| 久久久久免费精品人妻一区二区| 国产精品爽爽va在线观看网站| 成人精品一区二区免费| 国产三级黄色录像| 欧美最黄视频在线播放免费| 天堂√8在线中文| 亚洲在线观看片| 99精品欧美一区二区三区四区| 国产成人av激情在线播放| 国产成年人精品一区二区| 国产色爽女视频免费观看| 国产精品,欧美在线| 国产在线精品亚洲第一网站| 99热只有精品国产| 国产色爽女视频免费观看| 亚洲精品色激情综合| 观看免费一级毛片| 在线免费观看的www视频| 亚洲人成网站高清观看| 99精品欧美一区二区三区四区| 免费观看精品视频网站| 国产精品三级大全| 久久久久久久久中文| 丰满乱子伦码专区| 好看av亚洲va欧美ⅴa在| 午夜亚洲福利在线播放| 成人永久免费在线观看视频| 亚洲欧美一区二区三区黑人| 日韩精品中文字幕看吧| 天天躁日日操中文字幕| 熟女电影av网| 久久久国产成人免费| 波多野结衣高清无吗| 老司机午夜十八禁免费视频| 亚洲乱码一区二区免费版| 欧美3d第一页| 国产一区在线观看成人免费| 亚洲国产精品久久男人天堂| 香蕉av资源在线| 国产探花在线观看一区二区| 亚洲av日韩精品久久久久久密| 国产亚洲精品综合一区在线观看| 亚洲av一区综合| 老熟妇乱子伦视频在线观看| www.www免费av| 999久久久精品免费观看国产| 久久久国产成人免费| 又粗又爽又猛毛片免费看| 国产日本99.免费观看| 精品久久久久久,| 人妻久久中文字幕网| 色综合亚洲欧美另类图片| 国产在线精品亚洲第一网站| 麻豆国产97在线/欧美| 搡老熟女国产l中国老女人| а√天堂www在线а√下载| 网址你懂的国产日韩在线| 久久久久九九精品影院| 12—13女人毛片做爰片一| bbb黄色大片| 国产野战对白在线观看| 欧美另类亚洲清纯唯美| 中文字幕人成人乱码亚洲影| 亚洲精品乱码久久久v下载方式 | 国产精品久久久久久久久免 | 麻豆久久精品国产亚洲av| 亚洲三级黄色毛片| 亚洲欧洲日产国产| 国产亚洲精品av在线| 成人午夜高清在线视频| av天堂中文字幕网| 在线免费十八禁| 中文欧美无线码| 亚洲av在线观看美女高潮| 五月玫瑰六月丁香| 人人妻人人看人人澡| 欧美一区二区亚洲| 九草在线视频观看| 性色avwww在线观看| 国内少妇人妻偷人精品xxx网站| 干丝袜人妻中文字幕| 菩萨蛮人人尽说江南好唐韦庄|