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

    河南省近30年植被指數(shù)與氣象因子的相關(guān)性分析

    2018-11-07 06:48:18郭磊王世東
    生態(tài)科學(xué) 2018年5期
    關(guān)鍵詞:年際植被指數(shù)降雨量

    郭磊, 王世東

    ?

    河南省近30年植被指數(shù)與氣象因子的相關(guān)性分析

    郭磊, 王世東*

    河南理工大學(xué)測(cè)繪與國(guó)土信息工程學(xué)院, 焦作 454000

    GIMMS NDVI3g數(shù)據(jù)因其具有半月合成及長(zhǎng)時(shí)間序列的優(yōu)勢(shì), 被國(guó)內(nèi)外學(xué)者廣泛應(yīng)用于植被指數(shù)研究。因此利用河南省1982—2013年該數(shù)據(jù)集和省內(nèi)及周邊共32個(gè)氣象站點(diǎn)的逐月氣象數(shù)據(jù), 采用均值法、克里金插值和相關(guān)性分析等方法, 分析了河南省近30年來的植被指數(shù)與主要?dú)庀笠蜃?溫度、降雨量)的相關(guān)性。結(jié)果表明: (1)近30年來, 河南省年平均NDVI呈波動(dòng)性上升趨勢(shì), 其年增長(zhǎng)率為0.002, 總體植被覆蓋增加; (2)在全球變暖的大環(huán)境下, 省內(nèi)年均氣溫顯著增加, 其年平均增長(zhǎng)率為0.038, 而降雨量則逐年減少, 其年下降率為2.151; (3)年平均NDVI與年平均溫度呈一般正相關(guān), 兩者增減基本同步, 而與降雨量的年際變化則存在一定的滯后性, 體現(xiàn)在當(dāng)年降雨增多作用于來年植被NDVI上升明顯。(4)在植被覆蓋較差的中部地區(qū), 植被NDVI與氣溫、降雨的相關(guān)性較高, 在植被覆蓋較好的東部和西部地區(qū), 其NDVI與氣溫、降雨的相關(guān)性相對(duì)較低, 該結(jié)果為土地遭破壞區(qū)域的生態(tài)恢復(fù)提供了理論參考。

    歸一化植被指數(shù)(NDVI); 氣象因子; 相關(guān)性分析; 生態(tài)恢復(fù)

    1 前言

    近些年來, 我國(guó)受氣候變暖的影響日益劇增[1]。作為中原大地的河南省, 其植被生長(zhǎng)變化與氣候的關(guān)系也越來越緊密。因此, 探討植被覆蓋變化與氣象因子之間的相互作用關(guān)系對(duì)植被保護(hù)、氣候變化和生態(tài)恢復(fù)等研究都有著重要的指示意義。但是, 受制于研究區(qū)面積和自然經(jīng)濟(jì)條件, 傳統(tǒng)的常規(guī)觀測(cè)手段很難實(shí)現(xiàn)對(duì)大面積區(qū)域內(nèi)的植被生長(zhǎng)狀況進(jìn)行觀測(cè)[2]。

    遙感技術(shù)的發(fā)展實(shí)現(xiàn)了對(duì)某一地區(qū)的植被覆蓋變化進(jìn)行動(dòng)態(tài)觀測(cè)[3]。其中, GIMMS NDVI3g數(shù)據(jù)因其時(shí)間跨度長(zhǎng)、覆蓋范圍廣, 被越來越多的應(yīng)用于植被與氣候變化研究。

    NDVI(歸一化植被指數(shù))是表示植被生長(zhǎng)狀況的系數(shù)[4], 是反映植被動(dòng)態(tài)變化的有效指標(biāo)[5]。學(xué)者們對(duì)NDVI與氣象因子的研究已取得顯著成果。例如, 在全國(guó)區(qū)域, 王茜等[6]利用GIMMS NDVI3g數(shù)據(jù)發(fā)現(xiàn)國(guó)內(nèi)NDVI的年內(nèi)變化顯著, 劉可等[7]發(fā)現(xiàn)近30年來我國(guó)東部季風(fēng)區(qū)的農(nóng)田和森林生態(tài)系統(tǒng)的植被活動(dòng)顯著增強(qiáng), 而西北的荒漠生態(tài)系統(tǒng)和東北的森林生態(tài)系統(tǒng)的植被活動(dòng)顯著下降; 在西北地區(qū), 劉振元等[8]對(duì)青藏高原植被進(jìn)行分區(qū), 發(fā)現(xiàn)不同區(qū)影響NDVI下降的氣候因子不同, 而趙紫薇[9]發(fā)現(xiàn)青藏高原植被有改善趨勢(shì), 同時(shí)陸晴等[10]也發(fā)現(xiàn)高原草地生長(zhǎng)季NDVI整體有增加趨勢(shì); 在東北地區(qū), 王彥穎[11]發(fā)現(xiàn)植被在過去32年里NDVI趨勢(shì)以負(fù)向變化為主, 郭金停等[12]卻發(fā)現(xiàn)在凍土區(qū)植被生長(zhǎng)季平均NDVI呈顯著增加趨勢(shì); 在西南地區(qū), 張笑鶴[13]發(fā)現(xiàn)絕大部分地區(qū)NDVI與溫度呈正相關(guān), 與降水以負(fù)相關(guān)為主, 而鄭朝菊等[14]卻發(fā)現(xiàn)年最大植被覆蓋度(FVC)與夏季均溫和秋季降雨的相關(guān)性最好。從以上的研究中可以發(fā)現(xiàn)影響NDVI的氣候因子具有空間異質(zhì)性, 在不同地區(qū), 氣候因子對(duì)NDVI的響應(yīng)程度不同, 然而針對(duì)我國(guó)腹地, 特別是以河南省為主的中原地區(qū)的NDVI與氣候因子的相關(guān)性研究相對(duì)較少。

    因此, 為進(jìn)一步了解我國(guó)中原地區(qū)的植被覆蓋變化特征, 本文利用1982—2013年GIMMS NDVI3g數(shù)據(jù)和省內(nèi)及周邊共32個(gè)氣象站點(diǎn)的逐月氣象資料(溫度、降雨量), 分析了河南省近30年來的NDVI生長(zhǎng)季(4—9月)和年際(1—12月)變化規(guī)律, 同時(shí)進(jìn)一步對(duì)溫度和降雨與NDVI的相關(guān)性進(jìn)行了分析。

    2 數(shù)據(jù)與研究方法

    2.1 數(shù)據(jù)來源

    GIMMS NDVI3g數(shù)據(jù)來自美國(guó)航天局全球監(jiān)測(cè)和模型研究組(Global Inventor Modeling and Mapping Studies, GIMMS)發(fā)布的半月合成數(shù)據(jù), 經(jīng)最大合成法合成所得的月數(shù)據(jù)能最大限度的減少云層和氣溶膠的影響[15], 其空間分辨率為8km, 時(shí)間跨度為1982—2013年, 且該數(shù)據(jù)已經(jīng)過投影轉(zhuǎn)換、輻射定標(biāo)、幾何校正等一系列處理[16]。

    氣象數(shù)據(jù)下載自中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http: //cdc.gov.cn/home.do)。根據(jù)河南省內(nèi)及周邊區(qū)域篩選出32個(gè)氣象站點(diǎn)得到研究區(qū)1982—2013年逐月氣象數(shù)據(jù), 包括月平均溫度、月降水量2個(gè)主要?dú)庀笠蜃拥臄?shù)據(jù)資料。

    2.2 數(shù)據(jù)分析

    (1) NDVI

    NDVI(歸一化植被指數(shù))是一種表示植被生長(zhǎng)狀況的系數(shù)。其多光譜遙感數(shù)據(jù)主要是近紅外和紅光波段, 計(jì)算原理是根據(jù)植被的光譜反射特征來評(píng)價(jià)地表植被覆蓋度[17]。NDVI的計(jì)算公式為:

    式中:是近紅外波段的數(shù)值,是紅光波段的數(shù)值。計(jì)算得到的NDVI值在[-1, 1]之間, 當(dāng)?shù)乇碛兄脖桓采w時(shí)其值為正, 且覆蓋度越高值越大; 當(dāng)?shù)乇肀辉茖踊蜓┱趽鯐r(shí)其值為負(fù); 當(dāng)?shù)乇頌槁懵稁r石或土壤時(shí)其值為零。

    (2) GIMMS NDVI3g

    GIMMS NDVI3g數(shù)據(jù)集是相對(duì)標(biāo)準(zhǔn)的數(shù)據(jù), 其通過轉(zhuǎn)化公式可以將柵格單元的灰度值轉(zhuǎn)化成NDVI真實(shí)值[18]。對(duì)圖像進(jìn)行旋轉(zhuǎn)、添加投影信息(WGS-84)后, 利用河南省的邊界矢量文件裁剪得到研究區(qū)圖像, 然后創(chuàng)建掩膜文件進(jìn)行掩膜處理并進(jìn)行歸一化, 歸一化完成之后可以進(jìn)行信息提取, 最后對(duì)圖像進(jìn)行最大合成法合成, 以上過程皆由IDL編程完成。河南省1982—2013年半月數(shù)據(jù)共768景, 合成月數(shù)據(jù)384景, 再由月數(shù)據(jù)合成季度數(shù)據(jù)128景, 最終由季度數(shù)據(jù)合成年數(shù)據(jù)32景。GIMMS NDVI3g數(shù)據(jù)合成方法如下:

    ①半月合成:

    式中:表示上半月,表示下半月,表示半月NDVI值。

    ②月合成:

    式中:表示月最大NDVI值,、表示月份。

    ③季度合成:

    式中:表示季度最大NDVI值,表示季度,表示月份。

    ④年合成:

    式中:表示年最大NDVI值,表示年份,表示月份。

    2.3 研究方法

    2.3.1 克里金插值法

    為使氣象數(shù)據(jù)和植被指數(shù)能進(jìn)行空間計(jì)算, 本文對(duì)溫度和降雨數(shù)據(jù)進(jìn)行克里金插值, 使離散的氣象站點(diǎn)數(shù)據(jù)轉(zhuǎn)化為連續(xù)的柵格數(shù)據(jù)[19], 得到河南省1982—2013年的年平均溫度柵格數(shù)據(jù)和年降雨量柵格數(shù)據(jù)。

    2.3.2 均值法

    為反映河南省近30年來的總體植被覆蓋變化特征, 本文采用年平均植被指數(shù)和生長(zhǎng)季平均植被指數(shù)來對(duì)河南省植被NDVI的變化進(jìn)行趨勢(shì)計(jì)算, 即計(jì)算相應(yīng)時(shí)間區(qū)間內(nèi)所有象元NDVI的平均值。同時(shí), 對(duì)于氣象站點(diǎn)數(shù)據(jù), 本文所采用的生長(zhǎng)季和年際劃分標(biāo)準(zhǔn)分別為4—9月份和1—12月份, 將由氣象站點(diǎn)所得的逐月氣象數(shù)據(jù), 對(duì)溫度累加求平均得生長(zhǎng)季平均溫度和年平均溫度, 對(duì)降雨量累加除氣象站點(diǎn)個(gè)數(shù)得生長(zhǎng)季降雨量和年降雨量。

    2.3.3 相關(guān)性分析法

    相關(guān)性分析是用來研究要素與要素之間相互作用強(qiáng)弱的關(guān)系量[20], 本文主要采用偏相關(guān)分析和復(fù)相關(guān)分析的方法。偏相關(guān)分析是指忽略其他因素影響, 僅考慮某一種要素對(duì)另一種要素作用時(shí)產(chǎn)生的影響程度, 其結(jié)果以偏相關(guān)系數(shù)的形式展現(xiàn), 偏相關(guān)系數(shù)越大, 表明這種要素對(duì)另一種要素的作用效果越強(qiáng), 兩者間的相關(guān)性就越大。當(dāng)要素類型多于兩個(gè)時(shí), 會(huì)遇到考慮其中多種要素同時(shí)對(duì)某一種要素作用的情況, 這時(shí)就會(huì)用到復(fù)相關(guān)分析, 同理, 復(fù)相關(guān)系數(shù)越大, 表明它們之間的相關(guān)性越高。

    本文在討論溫度和降雨與植被指數(shù)之間的響應(yīng)關(guān)系時(shí), 采用偏相關(guān)分析和復(fù)相關(guān)分析的方法, 通過討論單一氣象因子對(duì)植被指數(shù)的影響和兩個(gè)氣象因子同時(shí)對(duì)植被指數(shù)影響的方法來分析氣象因子與植被指數(shù)的相關(guān)性。其計(jì)算公式如下:

    (1)偏相關(guān)分析

    式中:k、k表示兩個(gè)變量,r表示兩個(gè)變量之間的相關(guān)系數(shù)。假如有三個(gè)變量、、,固定求變量與變量的偏相關(guān)系數(shù), 其計(jì)算公式為:

    式中:r、r、r表示變量?jī)蓛砷g的相關(guān)系數(shù),r表示變量、間的偏相關(guān)系數(shù)。

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

    式中:為因變量,、為自變量,r為k與、的復(fù)相關(guān)系數(shù)。

    3 結(jié)果與分析

    3.1 NDVI的時(shí)間變化特點(diǎn)

    為研究河南省近30年來整體植被指數(shù)年際間的變化趨勢(shì), 本文主要從植被生長(zhǎng)季平均NDVI和年平均NDVI的變化特征進(jìn)行分析。生長(zhǎng)季是植被在全年里生長(zhǎng)變化特征最明顯的時(shí)段, 分析生長(zhǎng)季植被的NDVI可以預(yù)測(cè)該年植被的長(zhǎng)勢(shì), 而植被NDVI的年際變化可以反映出全年植被覆蓋的增減。如圖2所示, 近30年來河南省植被NDVI逐生長(zhǎng)季變化顯著, 變化范圍在0.656—0.780之間, 其增長(zhǎng)率為0.002, 在2011年達(dá)到峰值0.780。NDVI的逐年變化特性也呈現(xiàn)出波動(dòng)性上升趨勢(shì), 其年增長(zhǎng)率為0.002, 變化范圍在0.717—0.820之間, 總體上與生長(zhǎng)季NDVI的變化特性一致。由生長(zhǎng)季平均NDVI和年平均NDVI都呈波動(dòng)性上升趨勢(shì), 表明河南省近30年來總體植被覆蓋度有所增加, 植被生長(zhǎng)狀況逐漸得到改善。

    圖1 河南省省內(nèi)及周邊氣象站點(diǎn)分布圖

    圖2 河南省1982—2013年生長(zhǎng)季平均NDVI和年平均NDVI變化趨勢(shì)圖

    3.2 NDVI的空間變化特點(diǎn)

    圖3反映了河南省1982到2013年多年平均NDVI的空間分布狀況。從圖中可以看出, 植被覆蓋最高的地區(qū)在河南省洛陽(yáng)、南陽(yáng)和三門峽三市交界處的伏牛山地區(qū), 其NDVI為0.859—0.952; 其次是河南省東部的大部分地區(qū), 其NDVI為0.742—0.859; 植被覆蓋較差的地區(qū)主要分布在河南省中部的鄭州市、平頂山市, 北部的濟(jì)源市、焦作市, 以及洛陽(yáng)北部和南陽(yáng)南部等地區(qū), 其NDVI在0.551—0.742之間。

    圖3 河南省1982—2013年平均NDVI空間分布圖

    3.3 河南省氣候因子的時(shí)空變化特點(diǎn)

    河南省1982到2013年的2種主要?dú)夂蛞蜃?溫度和降雨量)在植被生長(zhǎng)季和年際的變化狀況分別如圖4、圖5所示。在生長(zhǎng)季(圖4), 平均溫度上升趨勢(shì)明顯, 增長(zhǎng)率為0.035, 降雨量逐生長(zhǎng)季也有輕微上升, 但遠(yuǎn)沒溫度上升明顯, 大部分年份的生長(zhǎng)季降雨量都保持在400 mm到800 mm之間; 而在年際(圖5), 平均溫度逐年明顯上升, 年增長(zhǎng)率為0.038, 與生長(zhǎng)季的溫度增長(zhǎng)率十分相近, 而年降雨量則呈下降趨勢(shì), 年減少率為2.151, 大部分年份的年降雨量保持在600 mm到1000 mm之間。

    在空間上, 對(duì)32年溫度和降雨量逐象元求平均得圖6和圖7。由圖6可以看出, 高溫地區(qū)主要集中在河南省的南部, 年平均氣溫可達(dá)15 ℃以上, 低溫地區(qū)主要集中在西部和北部的部分地區(qū), 年平均溫度在13 ℃左右, 而中部地區(qū)的年平均溫度則介于兩者之間; 由圖7可得, 河南省年降雨量由南向北逐漸減少, 南部地區(qū)年降雨量較高, 在1000 mm左右, 中部處于過渡區(qū), 在600—1000 mm之間, 而北部地區(qū)則在600 mm以下。

    圖4 1982—2013年生長(zhǎng)季溫度和降雨變化曲線

    Fig. 4 1982-2013 curve of growing season temperature and rainfall

    圖5 1982—2013年年際溫度和降雨變化曲線

    Fig. 5 1982-2013 curve of interannual temperature and rainfall

    圖6 河南省1982—2013年平均溫度空間分布圖

    3.4 NDVI對(duì)氣候因子的響應(yīng)

    3.4.1 NDVI與溫度和降雨的年際關(guān)系

    圖8和圖9是河南省1982—2013年植被NDVI分別與溫度和降水的生長(zhǎng)季和年際變化曲線。如圖8所示, 河南省近30年來植被NDVI在生長(zhǎng)季和年際的變化趨勢(shì)都和相應(yīng)時(shí)段溫度的變化趨勢(shì)相似, 兩者在大部分年份具有相同的增減性, 但也有少部分年份不符合這一特性。這一方面說明溫度的升高確實(shí)促進(jìn)了植被生長(zhǎng), 但另一方面也表明溫度并不是影響植被NDVI變化的唯一限制因子。由圖9可以看出, NDVI在生長(zhǎng)季和年際間的變化趨勢(shì)與相應(yīng)時(shí)段的降雨量變化存在一定的滯后性, 大部分年份NDVI隨降雨量的增加呈現(xiàn)減小趨勢(shì), 但作用于來年植被NDVI上升明顯, 這表明植被對(duì)降水更加敏感。

    3.4.2 NDVI與溫度和降雨的空間關(guān)系

    一般而言, 氣候變化和人類活動(dòng)是影響植被覆蓋變化的主要因子, 其中, 氣候變化更是作為驅(qū)動(dòng)力, 通過改變其自身的條件從而對(duì)植被產(chǎn)生影響[21–22], 然而僅在時(shí)序上討論其對(duì)NDVI的作用并不能反映出植被覆蓋變化的空間差異, 在空間上對(duì)其與NDVI的相互作用進(jìn)行分析能更直接的反映出它們之間的相關(guān)性。因此, 本文對(duì)1982—2013年河南省植被NDVI與年平均溫度、年降雨量進(jìn)行基于象元的偏相關(guān)性分析和復(fù)相關(guān)分析。由圖10可以看出, 河南省絕大部分區(qū)域的NDVI都與年平均溫度在空間上呈現(xiàn)正相關(guān), 且在中部地區(qū)相關(guān)性最高, 達(dá)到0.8以上, 僅在西部的三門峽、洛陽(yáng)和南陽(yáng)三市交界處以及北部的安陽(yáng)和南部的信陽(yáng)的少量區(qū)域內(nèi)呈現(xiàn)出負(fù)相關(guān); 由圖11可以看出, NDVI與年降雨量呈正相關(guān)的區(qū)域也占了河南省的大部分面積, 但要少于溫度與NDVI呈正相關(guān)的區(qū)域面積, 降雨對(duì)NDVI影響最顯著的區(qū)域主要在河南省東北部的安陽(yáng)市、濮陽(yáng)市、新鄉(xiāng)市以及中部的鄭州、平頂山與洛陽(yáng)三市交匯之地和南陽(yáng)南部的部分地區(qū), 其相關(guān)系數(shù)達(dá)0.6以上, 而在信陽(yáng)市北部和駐馬店南部以及周口西北部的NDVI顯示出與降雨量成明顯負(fù)相關(guān); 圖12是NDVI與年平均溫度和年降雨量進(jìn)行復(fù)相關(guān)分析的結(jié)果圖, 從圖中可以看出在, 在河南省中部地區(qū), NDVI與溫度和降雨的復(fù)相關(guān)系數(shù)較高, 其值在0.8左右, 而在東部和西部地區(qū)的復(fù)相關(guān)系數(shù)相對(duì)較低, 其值在0.5以下。

    圖7 河南省1982—2013年年均降雨量空間分布圖

    圖8 生長(zhǎng)季和年際NDVI與溫度的變化曲線

    圖9 生長(zhǎng)季和年際NDVI與降雨量的變化曲線

    綜上所述, 溫度和降雨與NDVI的相關(guān)性存在空間差異。在河南省內(nèi)的絕大部分地區(qū), 兩種氣象因子與NDVI都呈現(xiàn)出正相關(guān), 且NDVI與溫度相關(guān)性要高于NDVI與降雨的相關(guān)性。同時(shí), 通過對(duì)溫度和降雨與NDVI進(jìn)行偏相關(guān)分析與復(fù)相關(guān)分析, 結(jié)果表明在河南省植被覆蓋較好的東部和西部地區(qū), 其偏相關(guān)系數(shù)和復(fù)相關(guān)系數(shù)相對(duì)較小, 而在植被覆蓋較差的中部地區(qū), 其偏相關(guān)系數(shù)和復(fù)相關(guān)系數(shù)相對(duì)較大。

    4 結(jié)論與討論

    本文基于GIMMS NDVI3g數(shù)據(jù), 結(jié)合1982—2013年氣象數(shù)據(jù), 利用多種分析和計(jì)算方法, 研究了河南省近30年來的NDVI時(shí)空變化規(guī)律, 并探討了溫度和降雨與NDVI的相關(guān)性, 所得結(jié)論如下:

    圖10 河南省平均NDVI與平均溫度偏相關(guān)系數(shù)空間分布圖

    圖11 河南省平均NDVI與年均降雨量偏相關(guān)系數(shù)空間分布圖

    (1) 河南省近30年來, 在植被生長(zhǎng)季和年際, 其NDVI都以0.002的增長(zhǎng)率不斷上升, 表明總體植被覆蓋有所增加。NDVI的空間分布存在較明顯的區(qū)域差異, 總體上呈現(xiàn)出東、西高, 中部低的分布特點(diǎn)。

    (2) 近30年來的氣象數(shù)據(jù)表明, 省內(nèi)年均溫度在以0.038的增長(zhǎng)率逐漸升高, 而年降雨量則以每年2.151的減少率在逐漸下降。在空間上表現(xiàn)為, 河南省東南地區(qū)的氣候條件要遠(yuǎn)遠(yuǎn)優(yōu)于中部以及西北地區(qū)。

    (3) 通過對(duì)NDVI與氣象因子的相關(guān)性分析, 發(fā)現(xiàn)河南省植被NDVI變化是受這兩種氣象因子共同作用的結(jié)果。河南省大部分面積的NDVI與年平均溫度和年降雨量均呈較明顯的正相關(guān)關(guān)系, 與年平均溫度的正相關(guān)系數(shù)最高達(dá)0.811, 與年降雨量的正相關(guān)系數(shù)最高達(dá)0.673, 溫度和降雨共同對(duì)NDVI作用的復(fù)相關(guān)系數(shù)最高達(dá)0.859??梢? 河南省NDVI年際間0.002的增長(zhǎng)率主要?dú)w功于年平均溫度的升高, 同時(shí)降雨也促進(jìn)了年平均NDVI的提升, 但其促進(jìn)作用沒有溫度的促進(jìn)效果大。

    圖12 河南省年均NDVI與年均溫度和年均降雨量復(fù)相關(guān)系數(shù)空間分布圖

    (4) 相關(guān)性分析結(jié)果表明, 在植被覆蓋較好的東部和西部地區(qū), 溫度和降雨與NDVI的復(fù)相關(guān)系數(shù)相對(duì)較低, 在0.002—0.489之間, 而在植被覆蓋較差的中部地區(qū), 其與NDVI的復(fù)相關(guān)系數(shù)相對(duì)較高, 在0.489—0.859之間。這表明, 在植被覆蓋較差區(qū)域, 溫度和降雨對(duì)植被生長(zhǎng)起著決定性作用, 該結(jié)果為生態(tài)環(huán)境遭破壞區(qū)域的植被恢復(fù)提供了理論依據(jù)。

    綜上所述, 本研究主要說明了河南省近30年來植被NDVI的時(shí)空變化特點(diǎn), 闡述了其與溫度和降雨這兩種主要?dú)庀笠蜃拥南嚓P(guān)關(guān)系。但由于缺乏人類活動(dòng)和自然災(zāi)害數(shù)據(jù), 也未考慮其它氣象因子影響, 因而存在部分年份的變化規(guī)律不符現(xiàn)象。這些有待后續(xù)進(jìn)一步研究, 為我國(guó)中原地區(qū)的植被恢復(fù)及生態(tài)重建提出合理建議。

    [1] 杜衛(wèi), 張?jiān)鲂? 郝雨杉, 等. 江蘇省植被覆蓋動(dòng)態(tài)變化及影響因素分析[J]. 科技通報(bào), 2017, 33(5): 24–29.

    [2] 阿多, 趙文吉, 宮兆寧, 等. 1981–2013華北平原氣候時(shí)空變化及其對(duì)植被覆蓋的影響[J]. 生態(tài)學(xué)報(bào), 2017, 37(2): 576–592.

    [3] 郭繼凱, 吳秀芹, 董貴華, 等. 基于MODIS/NDVI的塔里木河流域植被覆蓋變化驅(qū)動(dòng)因素相對(duì)作用分析[J]. 干旱區(qū)研究, 2017, 34(3): 621–629.

    [4] ECKERT S, HUSLER F, LINIGER H, et al. Trend analysis of MODIS NDVI time series for detecting land degradation and regeneration in Mongolia[J]. Journal of Arid Environments, 2015, 113(2): 16–28.

    [5] 徐麗萍, 郭鵬, 劉林, 等. 天山北坡NDVI對(duì)氣候因子響應(yīng)的敏感性分析[J]. 湖北農(nóng)業(yè)科學(xué), 2014, 53(21): 5116–5125.

    [6] 王茜, 陳瑩, 阮璽睿, 等. 1982–2012年中國(guó)NDVI變化及其與氣候因子的關(guān)系[J]. 草地學(xué)報(bào), 2017, 25(4): 691–700.

    [7] 劉可, 杜靈通, 候靜, 等. 近30年中國(guó)陸地生態(tài)系統(tǒng)NDVI時(shí)空變化特征[J]. 生態(tài)學(xué)報(bào), 2017, 38(6): 1–12.

    [8] 劉振元, 張杰, 陳立. 青藏高原植被指數(shù)最新變化特征及其與氣候因子的關(guān)系[J]. 氣候與環(huán)境研究, 2017, 22(3): 289–300.

    [9] 趙紫薇. 1982–2013年青藏高原植被動(dòng)態(tài)變化時(shí)序分析[J]. 測(cè)繪科學(xué), 2017, 42(6): 62–70.

    [10] 陸晴, 吳紹洪, 趙東升. 1982–2013年青藏高原高寒草地覆蓋變化及與氣候之間的關(guān)系[J]. 地理科學(xué), 2017, 37(2): 292–300.

    [11] 王彥穎. 中國(guó)東北植被時(shí)空動(dòng)態(tài)變化及其對(duì)氣候響應(yīng)研究[D]. 長(zhǎng)春: 東北師范大學(xué), 2016.

    [12] 郭金停, 胡遠(yuǎn)滿, 熊在平, 等. 中國(guó)東北多年凍土區(qū)土地植被生長(zhǎng)季NDVI時(shí)空變化及其對(duì)氣候變化的響應(yīng)[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2017, 28(8): 2413–2422.

    [13] 張笑鶴. 西南地區(qū)NDVI和NPP時(shí)空動(dòng)態(tài)及其與氣候因子相關(guān)性分析[D]. 北京: 中國(guó)林業(yè)科學(xué)研究院, 2011.

    [14] 鄭朝菊, 曾源, 趙玉金, 等. 近15年中國(guó)西南地區(qū)植被覆蓋度動(dòng)態(tài)變化[J]. 國(guó)土資源遙感, 2017, 29(3): 128–136.

    [15] BECK H E, MCVICAR T R, DIJK A I J M V, et al. Global evaluation of four AVHRR–NDVI data sets: Intercomparison and assessment against Landsat imagery[J]. Remote Sensing of Environment, 2011, 115(10): 2547–2563.

    [16] 時(shí)忠杰, 高吉喜, 徐麗宏. 內(nèi)蒙古地區(qū)近25年植被對(duì)氣溫和降水變化的影響[J]. 生態(tài)環(huán)境學(xué)報(bào), 2011, 20(11): 1594–1601.

    [17] 郭鈮. 植被指數(shù)及其研究進(jìn)展[J]. 干旱氣象, 2003, 21(4): 71–75.

    [18] 單楠. 京津風(fēng)沙源區(qū)植被指數(shù)(NDVI)對(duì)氣候變化響應(yīng)研究[D]. 北京: 中國(guó)林業(yè)科學(xué)研究院, 2013.

    [19] 蔣育昊, 劉鵬舉, 夏智武, 等. 站點(diǎn)密度對(duì)復(fù)雜地形PRISM月降雨空間插值精度的影響[J].南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2017, 41(4): 115–120.

    [20] 潘秋玲. 那曲地區(qū)植被指數(shù)時(shí)空變化及其影響因素分析[D]. 北京: 中國(guó)地質(zhì)大學(xué)(北京), 2017.

    [21] 華東師范大學(xué). 植物學(xué)[M]. 北京: 高等教育出版社, 1982.

    [22] 李娜. 1999–2006石羊河流域植被對(duì)氣候變化的響應(yīng)研究[D]. 蘭州: 蘭州大學(xué), 2010.

    Correlation analysis of vegetation index and meteorological factors in Henan Province in the last 30 years

    GUO Lei, WANG Shidong*

    School of Surveying and mapping and land information engineering, Henan Polytechnic University, Jiaozuo 454000, China

    GIMMS NDVI3gdata has been widely used by scholars both at home and abroad for the study of vegetation index because of its advantages of semilunar synthesis and long time series. This paper used the 1982-2013 year GIMMS NDVI3gdata set in Henan and the data of 32 meteorological stations in and around the province, andadopted many methodsabout average, Kriging interpolation, correlation analysis and so on toanalyze the correlation between the vegetation index and the main meteorological factors (temperature, rainfall) in Henan Province in the last 30 years. The results are as follows. (1) In the past 30 years, the average annual NDVI in Henan province was fluctuating upward trend, the annual growth rate of 0.002; the overall vegetation coverage increased. (2) In the environment of global warming, the annual average temperature increased significantly, with the average annual growth rate of 0.038; while the rainfall was decreasing year by year and it fell at an annual rate of 2.151. (3) The annual average NDVI was generally positively correlated with annual mean temperature change, and the annual variation of precipitation was lagging behind. (4) In the areas of poor vegetation cover, there was a high correlation between NDVI and temperature or rainfall.In the area of vegetation cover was good, the correlation between NDVI and temperature or rainfall was relatively low. The result can provide a theoretical reference for ecological restoration.

    normalized vegetation index (NDVI); meteorological factors; correlation analysis; ecological restoration

    10.14108/j.cnki.1008-8873.2018.05.014

    S46

    A

    1008-8873(2018)05-102-11

    2018-03-02;

    2018-07-27

    國(guó)家自然科學(xué)基金項(xiàng)目(41301617); 中國(guó)博士后科學(xué)基金(2016M590679); 河南省高等學(xué)校重點(diǎn)科研項(xiàng)目(17A420001); 河南省高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金(NSFRF1630); 河南理工大學(xué)創(chuàng)新性科研團(tuán)隊(duì)項(xiàng)目(B2017-16); 河南省高校科技創(chuàng)新團(tuán)隊(duì)支持計(jì)劃(18IRTSTHN008)資助

    郭磊(1993—), 碩士生, 主要研究方向?yàn)閲?guó)土資源遙感, E-mail: 1174283475@qq.com

    通信作者:王世東(1978—), 副教授, 博士, 碩士生導(dǎo)師, 主要研究方向?yàn)閲?guó)土資源信息化與遙感, 土地復(fù)墾, E-mail: wsd0908@163.com

    郭磊, 王世東. 河南省近30年植被指數(shù)與氣象因子的相關(guān)性分析[J]. 生態(tài)科學(xué), 2018, 37(5): 102-112.

    GUO Lei, WANG Shidong. Correlation analysis of vegetation index and meteorological factors in Henan Province in the last 30 years[J]. Ecological Science, 2018, 37(5): 102-112.

    猜你喜歡
    年際植被指數(shù)降雨量
    北緯30°中層頂區(qū)域鈉與鐵原子層的結(jié)構(gòu)和年際變化
    降雨量與面積的關(guān)系
    AMSR_2微波植被指數(shù)在黃河流域的適用性對(duì)比與分析
    河南省冬小麥產(chǎn)量遙感監(jiān)測(cè)精度比較研究
    亞洲夏季風(fēng)的年際和年代際變化及其未來預(yù)測(cè)
    與北大西洋接壤的北極海冰和年際氣候變化
    洞庭湖區(qū)降雨特性分析
    主要植被指數(shù)在生態(tài)環(huán)評(píng)中的作用
    西藏科技(2015年1期)2015-09-26 12:09:29
    基于MODIS數(shù)據(jù)的植被指數(shù)與植被覆蓋度關(guān)系研究
    基于M-K法對(duì)圖們江干流含沙量年際變化的分析
    精品久久久精品久久久| 亚洲五月天丁香| 国产精品综合久久久久久久免费 | bbb黄色大片| 国产成年人精品一区二区| 色在线成人网| 国产精品一区二区精品视频观看| 伦理电影免费视频| 亚洲一区高清亚洲精品| 国产成人影院久久av| 欧美日韩乱码在线| 啦啦啦 在线观看视频| 一a级毛片在线观看| 大型黄色视频在线免费观看| 亚洲人成网站在线播放欧美日韩| 国产精品野战在线观看| 亚洲第一欧美日韩一区二区三区| 制服诱惑二区| 日本 欧美在线| 亚洲av熟女| 免费高清在线观看日韩| 久久精品91蜜桃| 一级作爱视频免费观看| 一进一出抽搐gif免费好疼| 亚洲欧美激情综合另类| 国产免费av片在线观看野外av| 久久热在线av| 母亲3免费完整高清在线观看| 国产野战对白在线观看| 欧美日韩亚洲国产一区二区在线观看| 午夜a级毛片| 50天的宝宝边吃奶边哭怎么回事| 久久久久久大精品| 欧美精品亚洲一区二区| 女人精品久久久久毛片| 亚洲成人国产一区在线观看| 久久国产精品影院| 亚洲色图综合在线观看| 国产精品亚洲av一区麻豆| 最近最新免费中文字幕在线| 免费观看人在逋| 黄色视频不卡| 国产精品免费一区二区三区在线| 村上凉子中文字幕在线| 黑人巨大精品欧美一区二区mp4| 日韩精品中文字幕看吧| 成人精品一区二区免费| 人人澡人人妻人| 亚洲国产中文字幕在线视频| 黄色片一级片一级黄色片| 亚洲欧美一区二区三区黑人| 亚洲精品粉嫩美女一区| 欧美黄色片欧美黄色片| 一级a爱片免费观看的视频| 精品少妇一区二区三区视频日本电影| 日本精品一区二区三区蜜桃| 亚洲欧美精品综合一区二区三区| 一级毛片女人18水好多| 国产精品 欧美亚洲| 国产国语露脸激情在线看| 午夜免费激情av| 性色av乱码一区二区三区2| 1024香蕉在线观看| 精品一品国产午夜福利视频| 日韩 欧美 亚洲 中文字幕| 久久狼人影院| 男女下面进入的视频免费午夜 | 国产欧美日韩综合在线一区二区| 久久久久久久精品吃奶| 亚洲av日韩精品久久久久久密| 深夜精品福利| 午夜福利一区二区在线看| 国产免费av片在线观看野外av| 国产精品久久久av美女十八| 亚洲精品一区av在线观看| 乱人伦中国视频| 老熟妇仑乱视频hdxx| 日本黄色视频三级网站网址| 在线观看66精品国产| 少妇 在线观看| 熟妇人妻久久中文字幕3abv| 不卡av一区二区三区| 中文字幕av电影在线播放| 99国产精品99久久久久| 丝袜美足系列| 国内精品久久久久精免费| 黄片播放在线免费| 最近最新中文字幕大全电影3 | 亚洲成人国产一区在线观看| 午夜福利一区二区在线看| 91麻豆精品激情在线观看国产| 男女下面进入的视频免费午夜 | 国产高清videossex| 午夜福利成人在线免费观看| 99久久精品国产亚洲精品| 国产亚洲欧美98| 亚洲第一欧美日韩一区二区三区| 国产97色在线日韩免费| 男女之事视频高清在线观看| 成人三级黄色视频| 欧美日韩瑟瑟在线播放| 多毛熟女@视频| 免费搜索国产男女视频| 国产成年人精品一区二区| 久久久久久免费高清国产稀缺| 中文字幕最新亚洲高清| 久久性视频一级片| 精品乱码久久久久久99久播| 一二三四在线观看免费中文在| 夜夜躁狠狠躁天天躁| 淫妇啪啪啪对白视频| 亚洲精品av麻豆狂野| 黄色毛片三级朝国网站| 一级毛片精品| 又黄又粗又硬又大视频| 国产av一区在线观看免费| av天堂久久9| 精品久久久久久久久久免费视频| 91九色精品人成在线观看| 亚洲专区字幕在线| 亚洲一区二区三区不卡视频| 欧美不卡视频在线免费观看 | 久久久国产精品麻豆| 黄色视频不卡| 黑丝袜美女国产一区| 亚洲性夜色夜夜综合| 男女下面插进去视频免费观看| 男女下面插进去视频免费观看| 国产成人精品在线电影| 长腿黑丝高跟| 天堂动漫精品| 亚洲熟妇熟女久久| 欧美日韩中文字幕国产精品一区二区三区 | 黄片大片在线免费观看| 一个人免费在线观看的高清视频| 少妇被粗大的猛进出69影院| 欧美在线黄色| 精品国产一区二区三区四区第35| 亚洲精品久久成人aⅴ小说| 欧美午夜高清在线| 精品乱码久久久久久99久播| 亚洲熟妇中文字幕五十中出| 99在线人妻在线中文字幕| 美女高潮喷水抽搐中文字幕| 日韩欧美一区二区三区在线观看| 欧美午夜高清在线| 99精品欧美一区二区三区四区| 日韩欧美一区二区三区在线观看| 国产主播在线观看一区二区| 免费av毛片视频| 久久狼人影院| 色综合亚洲欧美另类图片| av视频免费观看在线观看| 国产精品久久久人人做人人爽| 亚洲av美国av| 免费在线观看完整版高清| 中文字幕人妻丝袜一区二区| 男女床上黄色一级片免费看| 免费高清视频大片| 一二三四社区在线视频社区8| 日本精品一区二区三区蜜桃| 欧美日本中文国产一区发布| 国产一卡二卡三卡精品| 黄色女人牲交| 亚洲精品美女久久av网站| 黄色视频不卡| 99精品欧美一区二区三区四区| 国产一区二区激情短视频| 国产精品九九99| 亚洲一码二码三码区别大吗| 亚洲七黄色美女视频| 这个男人来自地球电影免费观看| 午夜精品久久久久久毛片777| 18禁黄网站禁片午夜丰满| 欧美黄色淫秽网站| 两人在一起打扑克的视频| 老司机福利观看| 怎么达到女性高潮| 超碰成人久久| 女同久久另类99精品国产91| 女同久久另类99精品国产91| 国产成人欧美| 男人舔女人的私密视频| 国产亚洲精品第一综合不卡| 久久国产乱子伦精品免费另类| 亚洲国产精品sss在线观看| 精品电影一区二区在线| 村上凉子中文字幕在线| 欧美在线黄色| 欧美老熟妇乱子伦牲交| 琪琪午夜伦伦电影理论片6080| 久久久久久国产a免费观看| 啪啪无遮挡十八禁网站| 精品熟女少妇八av免费久了| 人妻久久中文字幕网| 亚洲在线自拍视频| 中文字幕高清在线视频| 一级a爱视频在线免费观看| 一区二区三区激情视频| 午夜福利高清视频| 亚洲av成人不卡在线观看播放网| 欧美av亚洲av综合av国产av| 亚洲国产看品久久| 深夜精品福利| 黄频高清免费视频| 欧美成人性av电影在线观看| 最近最新中文字幕大全免费视频| 麻豆av在线久日| 亚洲人成伊人成综合网2020| 多毛熟女@视频| 伦理电影免费视频| 亚洲人成电影免费在线| 如日韩欧美国产精品一区二区三区| 人人妻人人澡人人看| 亚洲国产中文字幕在线视频| 国产av一区在线观看免费| 正在播放国产对白刺激| 大香蕉久久成人网| 亚洲色图 男人天堂 中文字幕| 99久久精品国产亚洲精品| 在线播放国产精品三级| 成人18禁高潮啪啪吃奶动态图| 国产精品一区二区免费欧美| 一区二区三区国产精品乱码| 91在线观看av| 在线观看免费视频日本深夜| 女人高潮潮喷娇喘18禁视频| 日韩 欧美 亚洲 中文字幕| АⅤ资源中文在线天堂| 嫩草影院精品99| 亚洲国产精品成人综合色| 高清毛片免费观看视频网站| 操出白浆在线播放| 精品久久久久久成人av| 日日夜夜操网爽| 色老头精品视频在线观看| 日韩欧美国产一区二区入口| 在线观看免费视频网站a站| 亚洲av成人av| 中文字幕另类日韩欧美亚洲嫩草| 99riav亚洲国产免费| 老司机午夜十八禁免费视频| 国产激情久久老熟女| 成人精品一区二区免费| 一级a爱视频在线免费观看| 亚洲专区中文字幕在线| 成人三级黄色视频| 国产欧美日韩精品亚洲av| cao死你这个sao货| 两性夫妻黄色片| 狂野欧美激情性xxxx| 高清在线国产一区| 美国免费a级毛片| 色哟哟哟哟哟哟| 国产精品久久久人人做人人爽| 免费观看人在逋| 别揉我奶头~嗯~啊~动态视频| 婷婷丁香在线五月| 久久青草综合色| 国内久久婷婷六月综合欲色啪| 色播在线永久视频| 777久久人妻少妇嫩草av网站| 欧美精品啪啪一区二区三区| 777久久人妻少妇嫩草av网站| 亚洲专区字幕在线| 国产精品1区2区在线观看.| 一进一出抽搐gif免费好疼| 精品久久久久久,| 少妇被粗大的猛进出69影院| 久久精品国产亚洲av高清一级| 在线观看免费午夜福利视频| 国产精品久久电影中文字幕| 久久国产亚洲av麻豆专区| 777久久人妻少妇嫩草av网站| 天天躁夜夜躁狠狠躁躁| 国产在线观看jvid| 淫秽高清视频在线观看| 亚洲人成电影观看| 亚洲欧美精品综合久久99| 欧美最黄视频在线播放免费| 欧美日韩瑟瑟在线播放| ponron亚洲| 美女免费视频网站| 亚洲久久久国产精品| 久久久精品欧美日韩精品| 国产91精品成人一区二区三区| av免费在线观看网站| 淫秽高清视频在线观看| 日日摸夜夜添夜夜添小说| 一边摸一边抽搐一进一出视频| 麻豆国产av国片精品| av网站免费在线观看视频| av片东京热男人的天堂| 国产高清videossex| 老司机在亚洲福利影院| 久久狼人影院| 久久中文字幕一级| 99在线视频只有这里精品首页| 国产精品 欧美亚洲| 女人爽到高潮嗷嗷叫在线视频| 亚洲成人免费电影在线观看| 久久婷婷人人爽人人干人人爱 | 99久久精品国产亚洲精品| 亚洲国产中文字幕在线视频| 韩国av一区二区三区四区| 老熟妇仑乱视频hdxx| 精品乱码久久久久久99久播| 国产精品 国内视频| 老司机午夜福利在线观看视频| 国产99久久九九免费精品| 一级黄色大片毛片| 亚洲自拍偷在线| 久久久精品欧美日韩精品| 变态另类丝袜制服| 国产精品自产拍在线观看55亚洲| 操出白浆在线播放| 国产av一区在线观看免费| 亚洲色图综合在线观看| 午夜福利高清视频| 十八禁网站免费在线| 成年人黄色毛片网站| 久久草成人影院| 欧美一级a爱片免费观看看 | 国产精品久久久人人做人人爽| 男女之事视频高清在线观看| 中文字幕人妻熟女乱码| 久久久久精品国产欧美久久久| 欧美国产精品va在线观看不卡| 国产一区在线观看成人免费| 国内毛片毛片毛片毛片毛片| bbb黄色大片| 成人欧美大片| 国产一级毛片七仙女欲春2 | 每晚都被弄得嗷嗷叫到高潮| av视频免费观看在线观看| 一级片免费观看大全| 国产成人啪精品午夜网站| 精品卡一卡二卡四卡免费| 伊人久久大香线蕉亚洲五| 国产亚洲欧美在线一区二区| 69精品国产乱码久久久| 久久 成人 亚洲| 久久人人97超碰香蕉20202| 国产精品自产拍在线观看55亚洲| 两性夫妻黄色片| 在线播放国产精品三级| 欧美日韩瑟瑟在线播放| 美女午夜性视频免费| 久久久久久久久久久久大奶| 正在播放国产对白刺激| 少妇 在线观看| 男女做爰动态图高潮gif福利片 | 日韩有码中文字幕| 精品久久久久久久人妻蜜臀av | 日韩欧美一区二区三区在线观看| 久久久国产精品麻豆| 99久久精品国产亚洲精品| 亚洲午夜理论影院| 最好的美女福利视频网| 国产精品自产拍在线观看55亚洲| 国产成人av教育| 看黄色毛片网站| 黄色 视频免费看| 亚洲av成人av| 如日韩欧美国产精品一区二区三区| 最新美女视频免费是黄的| 国产亚洲av高清不卡| 国产精品免费一区二区三区在线| 电影成人av| 在线播放国产精品三级| 精品福利观看| 国产不卡一卡二| 两个人视频免费观看高清| 欧美成人免费av一区二区三区| 天天一区二区日本电影三级 | 两个人视频免费观看高清| 大码成人一级视频| 动漫黄色视频在线观看| 欧美日韩精品网址| 亚洲美女黄片视频| 亚洲成av片中文字幕在线观看| 香蕉国产在线看| 久久热在线av| 亚洲 欧美 日韩 在线 免费| 亚洲熟妇熟女久久| 操出白浆在线播放| 国产在线精品亚洲第一网站| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美国产精品va在线观看不卡| 女人高潮潮喷娇喘18禁视频| 国产精品,欧美在线| 桃红色精品国产亚洲av| 国产人伦9x9x在线观看| 国产乱人伦免费视频| 国产91精品成人一区二区三区| 亚洲国产欧美一区二区综合| 久久久国产成人精品二区| 精品久久久精品久久久| 欧美av亚洲av综合av国产av| 亚洲一码二码三码区别大吗| 精品久久蜜臀av无| 亚洲色图av天堂| 精品人妻在线不人妻| 欧美乱妇无乱码| 热re99久久国产66热| 波多野结衣高清无吗| 亚洲欧美一区二区三区黑人| 亚洲国产中文字幕在线视频| 国产精品98久久久久久宅男小说| 一a级毛片在线观看| 男女之事视频高清在线观看| 久久伊人香网站| 在线观看www视频免费| 一区二区日韩欧美中文字幕| 一级片免费观看大全| 精品第一国产精品| 18禁观看日本| 免费看美女性在线毛片视频| 大码成人一级视频| 老熟妇乱子伦视频在线观看| 国产免费av片在线观看野外av| 久久精品91无色码中文字幕| 国产精品一区二区在线不卡| 国产视频一区二区在线看| 一边摸一边抽搐一进一出视频| 久久精品aⅴ一区二区三区四区| 欧美色视频一区免费| 亚洲va日本ⅴa欧美va伊人久久| 国产私拍福利视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 国产区一区二久久| 美女国产高潮福利片在线看| 天天躁夜夜躁狠狠躁躁| 一个人免费在线观看的高清视频| 欧美中文综合在线视频| 一区二区三区高清视频在线| 制服丝袜大香蕉在线| 国产片内射在线| 电影成人av| 性欧美人与动物交配| 国产成人精品久久二区二区免费| 搡老岳熟女国产| 亚洲五月婷婷丁香| 亚洲午夜理论影院| 久久天堂一区二区三区四区| 国产极品粉嫩免费观看在线| 欧美 亚洲 国产 日韩一| 丰满的人妻完整版| 国产麻豆成人av免费视频| 九色亚洲精品在线播放| 欧美日韩中文字幕国产精品一区二区三区 | 嫩草影视91久久| av视频在线观看入口| 国产激情欧美一区二区| 美女午夜性视频免费| 给我免费播放毛片高清在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 这个男人来自地球电影免费观看| 99精品在免费线老司机午夜| 美女午夜性视频免费| 波多野结衣一区麻豆| 一区二区三区国产精品乱码| 在线视频色国产色| 一个人观看的视频www高清免费观看 | 亚洲伊人色综图| 国产亚洲精品久久久久5区| 精品第一国产精品| 久久人人精品亚洲av| av视频在线观看入口| 两个人看的免费小视频| 91字幕亚洲| 精品免费久久久久久久清纯| 九色国产91popny在线| 国产av在哪里看| 欧美成人一区二区免费高清观看 | 两个人视频免费观看高清| 日日爽夜夜爽网站| 91字幕亚洲| 久久人妻福利社区极品人妻图片| 少妇的丰满在线观看| 亚洲国产精品合色在线| 在线观看一区二区三区| 日本一区二区免费在线视频| 黄色成人免费大全| 制服丝袜大香蕉在线| 亚洲色图综合在线观看| 久久精品国产清高在天天线| 精品午夜福利视频在线观看一区| 亚洲 欧美 日韩 在线 免费| 欧美激情高清一区二区三区| 正在播放国产对白刺激| 丁香欧美五月| 欧美日韩一级在线毛片| 少妇 在线观看| 亚洲av第一区精品v没综合| 国产99久久九九免费精品| 久久婷婷人人爽人人干人人爱 | 国产一级毛片七仙女欲春2 | 日本免费a在线| 一级毛片精品| 男女做爰动态图高潮gif福利片 | 国产激情久久老熟女| 人人妻人人澡欧美一区二区 | 久久香蕉精品热| 国产伦人伦偷精品视频| av天堂久久9| 99re在线观看精品视频| 久久香蕉激情| 欧美成人午夜精品| 免费女性裸体啪啪无遮挡网站| 国产精品野战在线观看| 他把我摸到了高潮在线观看| 狂野欧美激情性xxxx| 日本五十路高清| 欧美色视频一区免费| 国产精品1区2区在线观看.| а√天堂www在线а√下载| 色播在线永久视频| 老汉色av国产亚洲站长工具| 午夜久久久久精精品| 亚洲第一av免费看| 老司机午夜十八禁免费视频| 免费在线观看影片大全网站| 久久久久国内视频| 日韩精品中文字幕看吧| 国产精品自产拍在线观看55亚洲| 久久青草综合色| 午夜福利,免费看| √禁漫天堂资源中文www| 18美女黄网站色大片免费观看| 91九色精品人成在线观看| 日本精品一区二区三区蜜桃| 在线播放国产精品三级| 中文字幕人成人乱码亚洲影| av片东京热男人的天堂| 亚洲 欧美 日韩 在线 免费| 免费搜索国产男女视频| 午夜福利影视在线免费观看| 国产精品日韩av在线免费观看 | 久久狼人影院| 国产成人欧美在线观看| 国产精品影院久久| 午夜老司机福利片| 亚洲国产欧美日韩在线播放| 高清毛片免费观看视频网站| 变态另类成人亚洲欧美熟女 | 亚洲av片天天在线观看| 黄色视频不卡| 可以在线观看毛片的网站| 午夜福利,免费看| 午夜久久久久精精品| av视频免费观看在线观看| 法律面前人人平等表现在哪些方面| 日本欧美视频一区| 国产区一区二久久| 老汉色∧v一级毛片| 香蕉国产在线看| 亚洲精品在线美女| 嫁个100分男人电影在线观看| 亚洲欧美精品综合一区二区三区| а√天堂www在线а√下载| 熟妇人妻久久中文字幕3abv| 黑丝袜美女国产一区| 黄片播放在线免费| 色综合亚洲欧美另类图片| 国产成人啪精品午夜网站| 亚洲国产看品久久| 性欧美人与动物交配| 亚洲av成人不卡在线观看播放网| 美女 人体艺术 gogo| 国产精品影院久久| 免费在线观看亚洲国产| 午夜免费成人在线视频| 国产熟女午夜一区二区三区| 黄网站色视频无遮挡免费观看| 久久精品91无色码中文字幕| 免费无遮挡裸体视频| 久久亚洲精品不卡| 久久久久亚洲av毛片大全| 热99re8久久精品国产| 免费在线观看影片大全网站| 久久久久久久久久久久大奶| 久久影院123| 午夜福利一区二区在线看| 精品无人区乱码1区二区| 黑人巨大精品欧美一区二区蜜桃| 午夜福利视频1000在线观看 | 69精品国产乱码久久久| 国产一区二区三区视频了| 大香蕉久久成人网| 日韩精品中文字幕看吧| 国产高清有码在线观看视频 | 精品日产1卡2卡| 精品久久久久久,| 欧美黄色片欧美黄色片| 欧洲精品卡2卡3卡4卡5卡区| 国产一区在线观看成人免费| 欧美日韩福利视频一区二区| ponron亚洲| 老鸭窝网址在线观看| 久久草成人影院| 国产精品综合久久久久久久免费 | 亚洲专区字幕在线| 欧美午夜高清在线| 中文字幕最新亚洲高清| 亚洲成av片中文字幕在线观看| 欧美大码av| 中文字幕人成人乱码亚洲影| 香蕉国产在线看| 一区二区三区高清视频在线| 一级,二级,三级黄色视频| 一二三四在线观看免费中文在| 大陆偷拍与自拍| 久久精品亚洲熟妇少妇任你| 99在线人妻在线中文字幕|