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

    南充市轄區(qū)近15年植被覆蓋度變化特征

    2019-04-02 02:04:32王磊王杰付林林葉彬
    生態(tài)科學(xué) 2019年1期
    關(guān)鍵詞:市轄區(qū)南充坡向

    王磊, 王杰, 付林, 林葉彬,3*

    ?

    南充市轄區(qū)近15年植被覆蓋度變化特征

    王磊1, 王杰1, 付林2, 林葉彬1,3*

    1. 西華師范大學(xué),南充 637009 2. 西南交通大學(xué),成都 611756 3. 區(qū)域環(huán)境演變與保護(hù)研究中心, 西華師范大學(xué),南充 637009

    植被覆蓋度是自然條件和社會(huì)經(jīng)濟(jì)活動(dòng)的綜合反映, 研究一個(gè)區(qū)域植被覆蓋度的時(shí)空變化對(duì)于認(rèn)識(shí)該地區(qū)的社會(huì)發(fā)展?fàn)顩r、制定更加合理的生態(tài)環(huán)境政策、實(shí)現(xiàn)經(jīng)濟(jì)可持續(xù)發(fā)展具有重要意義。基于2002年Landsat 7 ETM、2014年Landsat 8影像和DEM數(shù)據(jù), 利用ENVI和ArcGIS軟件進(jìn)行處理, 運(yùn)用歸一化植被指數(shù)(NDVI)方法對(duì)南充市轄區(qū)植被覆蓋度進(jìn)行估算, 并劃分為5個(gè)不同蓋度等級(jí)。根據(jù)DEM數(shù)據(jù)提取海拔、坡度和坡向因子, 然后將植被覆蓋度圖與海拔、坡度、坡向圖進(jìn)行疊加, 定量分析了南充市轄區(qū)植被的空間分布特征及其變化情況。結(jié)果表明, 南充市轄區(qū)植被覆蓋度在2002—2014年總體變化不大, 其中以Ⅴ級(jí)(fc ≥ 0.8)植被覆蓋度為主, 其覆蓋面積達(dá)42.49%。2002—2014年, Ⅳ級(jí)和Ⅴ級(jí)增加了60 km2和191 km2, 增幅率分別為7.61%和21.39%, 而Ⅰ、Ⅱ和Ⅲ級(jí)分別減少了84 km2、124 km2和42 km2, 降幅為31.7%、56.62%和10.91%。其中植被集中分布在海拔300 m—400 m和坡度2°—6°的地帶, 植被覆蓋度坡向差異呈現(xiàn)半陽坡>半陰坡>陽坡>陰坡的特征。

    植被覆蓋度; 歸一化植被指數(shù); 空間分布; 南充市轄區(qū)

    0 前言

    植被覆蓋度是指植被(包括葉、莖、枝)在單位面積內(nèi)垂直投影面積所占比例[1]。植被作為生態(tài)系統(tǒng)的重要組成部分, 是生態(tài)系統(tǒng)中物質(zhì)循環(huán)與能量流動(dòng)的中樞, 是隨人類社會(huì)經(jīng)濟(jì)活動(dòng)有重要貢獻(xiàn)的資源, 是生態(tài)系統(tǒng)存在的基礎(chǔ), 也是連接土壤、大氣和水分的自然“紐帶”[2]。在生態(tài)環(huán)境系統(tǒng)中, 植被覆蓋狀況在很大程度上影響甚至直接決定著區(qū)域生態(tài)環(huán)境中的第一性生產(chǎn)力、環(huán)境承載力、環(huán)境潔凈與美化、水土流失強(qiáng)度等生態(tài)環(huán)境系統(tǒng)的狀態(tài)與功能[3]。因此, 植被覆蓋度作為衡量地表植被覆蓋的一個(gè)最重要的指標(biāo), 對(duì)水文、生態(tài)、全球變化等都具有重要的意義[4]。

    目前利用遙感影像數(shù)據(jù)來提取植被覆蓋度的方法主要有經(jīng)驗(yàn)?zāi)P头ā⒅脖恢笖?shù)法和像元法[5-6]。經(jīng)驗(yàn)?zāi)P头ㄖ饕峭ㄟ^建立實(shí)測(cè)植被覆蓋率數(shù)據(jù)與植被指數(shù)的經(jīng)驗(yàn)?zāi)P蛠砬笕〈竺娣e植被覆蓋率[7], 又可分為線性回歸模型法與非線性回歸模型法。線性回歸模型法, 如: Graetz與Pech對(duì)Landsat Miss第5波段與植被覆蓋度的實(shí)測(cè)數(shù)據(jù)進(jìn)行線性回歸分析, 并利用得到的回歸模型測(cè)量了稀疏草地的植被覆蓋度[8]。Dymond等為研究新西蘭退化草地的植被覆蓋度, 使用NDVI計(jì)算出植被指數(shù), 并用該植被指數(shù)與植被覆蓋度進(jìn)行擬合, 得到非線性回歸模型[9]。經(jīng)驗(yàn)?zāi)P头ㄖ贿m用于特定的區(qū)域與特定的植被類型, 一般不易推廣, 不具有普遍意義[10]。植被指數(shù)法是指直接利用植被指數(shù)近似估算植被覆蓋度, 所使用的植被指數(shù)一般都經(jīng)過驗(yàn)證與植被覆蓋度具有良好的相關(guān)關(guān)系, 不需要建立回歸模型。楊勝天等使用NDVI將植被覆蓋度分為四種類型: 高覆蓋度類型(植被覆蓋度>75%), 中高覆蓋度類型(植被覆蓋度在60%與75%之間)、中覆蓋度類型(植被覆蓋度在45%與60%之間)、低覆蓋度類型(植被覆蓋度≤45%)[11]; Muriithi等使用NDVI來研究肯尼亞中部高地植被狀況的變化[12]。Worku Zewdie等通過MODIS影像計(jì)算NDVI, 通過NDVI監(jiān)測(cè)埃塞俄比亞西北部的生態(tài)系統(tǒng)動(dòng)態(tài)和氣候變量來評(píng)估旱地植被的長期趨勢(shì)變化性[13]; Chaithanya等用Landsat 7 ETM圖像研究NDVI和LST的關(guān)系和空間變化[14]。像元分解模型法目前已開發(fā)出的模型主要有種, 分別為線性模型、概率模型、幾何光學(xué)模型、隨機(jī)幾何模型和模糊分析模型[15]。其中最常用的就是線性模型中的像元二分模型。李苗苗在像元二分模型的基礎(chǔ)上, 對(duì)已有的模型進(jìn)行改進(jìn), 利用NDVI定量估算植被覆蓋度的模型[4]。田義超等利用像元二分模型對(duì)北部灣沿海地區(qū)植被覆蓋計(jì)算, 及植被覆蓋對(duì)氣溫和降水的旬響應(yīng)特征[16]。Li等用像元二分模型計(jì)算河南省的植被覆蓋度, 再結(jié)合地形因子計(jì)算各個(gè)坡度、坡向和海拔上的植被覆蓋度的差異[17]; Wang等利用2002—2013年的Landsat影像, 根據(jù)植被覆蓋度變化對(duì)上海的地表城市熱島進(jìn)行分析[18]。

    大量研究表明, NDVI在使用遙感圖像進(jìn)行植被研究以及植物物候研究中得到廣泛應(yīng)用, 它是植物生長狀態(tài)以及植被空間分布密度的最佳指示因子[19]。很多學(xué)者通過NDVI來描述一個(gè)地區(qū)的植被覆蓋度, 并分析植被覆蓋對(duì)氣象因子的響應(yīng)特征, 以及對(duì)植被變化的驅(qū)動(dòng)力進(jìn)行分析, 像元二分模型也廣泛應(yīng)用于植被覆蓋度的提取。

    南充市轄區(qū)位于四川盆地, 全年云量較大, 影像大多難以利用, 少有學(xué)者對(duì)該區(qū)域植被覆蓋度開展研究工作。本文應(yīng)用像元二分模型、遙感與GIS技術(shù)對(duì)南充市轄區(qū)的植被覆蓋度進(jìn)行估算, 結(jié)合地形(海拔、坡度、坡向)量化了植被覆蓋度與地形因子的關(guān)系, 對(duì)研究區(qū)域過去15年來的植被覆蓋變化進(jìn)行了初步探究。

    1 研究區(qū)概況

    南充市轄區(qū)包括順慶區(qū)、嘉陵區(qū)和高坪區(qū), 位于南充市境內(nèi)(30°35′—31°51′N, 105°27′—106°58′E), 為亞熱帶溫暖濕潤季風(fēng)氣候區(qū), 四季分明, 氣候溫和, 雨量充沛。年均氣溫17℃, 年均降水量1100 mm, 地形以低山、丘陵和平壩為主。植被以亞熱帶常綠闊葉林和亞熱帶落葉闊葉林為主, 境內(nèi)多為馬尾松、柏樹和杉木等構(gòu)成的林地。林地遭到破壞后, 形成黃荊、馬桑、小果薔薇、白櫟等組成的灌叢。

    南充森林覆蓋率在經(jīng)歷了大煉鋼鐵和農(nóng)村辦公共食堂后, 樹木大量減少, 森林覆蓋率大幅度降低。在“文化大革命”期間, 再次出現(xiàn)大面積的亂砍濫伐, 南充全區(qū)林業(yè)生產(chǎn)急轉(zhuǎn)直下, 全市森林覆蓋度下降到歷史最低點(diǎn), 出現(xiàn)大量的水土流失, 生態(tài)失調(diào), 環(huán)境惡化。在中國共產(chǎn)黨第十一屆三中全會(huì)后, 南充市的森林覆蓋率不斷上升, 20世紀(jì)80年代后, 南充市興林業(yè)經(jīng)濟(jì)體制改革, 廣大群眾造林積極性空前高漲, 速生豐產(chǎn)林得到快速發(fā)展。同時(shí)該時(shí)期南充市加大了生態(tài)環(huán)境保護(hù)的投入, 實(shí)施了包括長江防護(hù)林、環(huán)城綠化、退耕還林和天然林保護(hù)等系列生態(tài)工程, 全市植被得到很好的恢復(fù)和發(fā)育, 到2003年, 順慶、高坪、嘉陵三區(qū)森林覆蓋率達(dá)到30.7%。

    2 數(shù)據(jù)來源和研究方法

    2.1 數(shù)據(jù)選擇與預(yù)處理

    本研究選擇2014年8月6日南充市轄區(qū)Landsat 8影像和2002年8月29日南充市轄區(qū)Landsat ETM+影像的多光譜影像為數(shù)據(jù)源, 分辨率為30 m, 軌道號(hào)為128/39, 云量為2.92%和0.03%, 選用DEM數(shù)字高程數(shù)據(jù), 分辨率為30 m, 數(shù)據(jù)均來自地理空間數(shù)據(jù)云(http://www.gscloud.cn)。對(duì)Landsat 8的多光譜影像在ENVI 5.1中進(jìn)行輻射定標(biāo)和大氣校正, 再根據(jù)研究區(qū)的矢量圖層對(duì)影像圖進(jìn)行裁剪。具體流程如下:

    (1)輻射定標(biāo)。在提取NDVI之前要對(duì)影像數(shù)據(jù)進(jìn)行輻射定標(biāo)。即將傳感器記錄的電壓或數(shù)字化值(DN)轉(zhuǎn)換成絕對(duì)輻射亮度值的過程或者轉(zhuǎn)換與地表反射率、表觀溫度等有關(guān)的相對(duì)值的處理過程[20]。本文采用的ENVI 5.1 Radiometric Correction模塊進(jìn)行輻射定標(biāo)。

    (2)大氣校正。由于電磁波在大氣中的傳輸和遙感器觀測(cè)過程中光照條件以及大氣作用的影響, 遙感器的測(cè)量值與地物實(shí)際的光譜輻射率是不一樣的, 使光譜分布發(fā)生變化, 測(cè)量值發(fā)生輻射失真[21]。所以, 大氣校正的目的是消除大氣和光照等因素對(duì)地物反射的影響。ENVI大氣校正模塊提供兩種大氣校正工具: FLAASH和Quick Atmospheric Correction兩種校正工具。本文采用FLAASH 模塊進(jìn)行大氣校正。得到大氣校正后的波普曲線圖(圖1)。

    (3)將shapefile格式面狀文件加載在ENVI 5.1中, 使用南充市轄區(qū)行政界線矢量數(shù)據(jù)對(duì)大氣校正后的影像進(jìn)行裁剪, 得到南充市轄區(qū)的影像(圖2)。

    2.2 NDVI計(jì)算

    NDVI根據(jù)植被反射波段的特性, 由遙感傳感器獲取的多光譜數(shù)據(jù)經(jīng)線性和非線性組合計(jì)算出來的各種數(shù)值, 對(duì)植被覆蓋具有一定的指示意義, 一般選用強(qiáng)吸收的可見光紅波段和強(qiáng)反射的近紅外波段組合計(jì)算[22]。

    即NDVI=(NIR-R)/(NIR+R) (1)

    2.3 植被覆蓋度的提取

    運(yùn)用歸一化植被指數(shù)(NDVI)估算一定區(qū)域的植被覆蓋度的研究有很多, 參考已有的研究成果, 本文采用公式(2)進(jìn)行計(jì)算, 即

    圖1 大氣校正前后光譜曲線對(duì)比

    Figure 1 Spectral curve comparison before and after atmospheric correction

    圖2 南充市轄區(qū)影像圖

    Figure 2 Image of Nanchong Jurisdiction

    FC=(NDVI-NDVImin)/(NDVImax-NDVImin)[6](2)

    式中: FC是植被覆蓋度, NDVImin和NDVImax分別為最小和最大NDVI值, 考慮到遙感影像中不可避免的存在噪聲, 它可能產(chǎn)生過低或過高的NDVI值。因此, 并不能直接取NDVI得最大最小值, 而是取一定置信度區(qū)間的最大最小值, 通過對(duì)影像的NDVI植被指數(shù)進(jìn)行直方圖分布圖分析, 在NDVI頻率累積表上取一定頻率的NDVI為NDVImin和NDVImax[19]。本文根據(jù)研究區(qū)的實(shí)際情況, 選擇5%和95%的NDVI值作為NDVImin和NDVImax(表1)。

    2.4 植被覆蓋度等級(jí)劃分

    NDVI的值在-1—1之間變動(dòng), NDVI<0表示地表覆蓋為水、雪等高反射的事物, NDVI=0表示巖石或者裸土等, NDVI>0表示植被覆蓋度的稀疏與濃密情況, NDVI的值越高, 說明植被覆蓋度越好。將得到的FC圖在ARCGIS中用重采樣的方法對(duì)研究區(qū)的植被覆蓋度進(jìn)行分類。本文結(jié)合研究區(qū)的實(shí)際情況, 并參考其他相關(guān)研究成果[23-24], 將南充市轄區(qū)的植被覆蓋度分為5級(jí)(表2)。

    2.5 植被覆蓋度分級(jí)結(jié)果

    按照公式(1)到(2)計(jì)算研究區(qū)的植被覆蓋度, 按表2進(jìn)行分級(jí), 并對(duì)各個(gè)分類進(jìn)行統(tǒng)計(jì), 得到2002年和2014年植被覆蓋度分級(jí)圖(圖3 a、b)。

    表1 研究區(qū)NDVI最小值和最大值

    2.6 地形因子類型劃分

    運(yùn)用ArcGIS 10.2軟件對(duì)DEM數(shù)據(jù)進(jìn)行裁剪, 對(duì)裁剪出的南充市轄區(qū)DEM進(jìn)行坡度、坡向提取, 然后對(duì)坡度、坡向和海拔進(jìn)行重分類。分別得到研究區(qū)的海拔、坡度、坡向圖(圖4), 研究區(qū)的坡度在0°—53.4612°之間, 根據(jù)研究區(qū)的實(shí)際情況, 將坡度劃分為5個(gè)等級(jí)[25]: 0°—2°, 2°—5°, 5°—15°, 15°—25°和>25°。海拔等級(jí)劃分, 研究區(qū)的海拔為191 m—805 m, 考慮到南充市轄區(qū)的具體的地勢(shì)情況, 將海拔劃分為5個(gè)等級(jí): <300 m, 300 m—400 m, 400 m—500 m, 500 m—600 m, >600 m。坡向劃分: 把坡向分為9中類型, 即: 平面(—1°), 北坡(0°—22.5°, 337.5°—360°)、東北坡(22.5°—67.5°)、東坡(67.5°—112.5°)、東南坡(112.5°—157.5°)、南坡(157.5°—202.5°)、西南坡(202.5°—247.5°)、西坡(247.5°—292.5°)和西北坡(292.5°—337.5°)。

    表2 植被覆蓋度分級(jí)

    3 結(jié)果與分析

    3.1 植被覆蓋度時(shí)空特征及面積統(tǒng)計(jì)

    通過2002年和2014年植被覆蓋度分級(jí)圖(圖3 a、b), 統(tǒng)計(jì)各級(jí)相應(yīng)的植被覆蓋度面積及其變化面積(表3)。由表3可知, 研究區(qū)2002和2014年植被覆蓋度都是Ⅴ級(jí)所占的比重最大, 分別為42.49%和35.02%, 說明研究區(qū)植被覆蓋狀況很好。2002—2014年研究區(qū)植被覆蓋度總體變化不大, 但其中Ⅳ級(jí)和Ⅴ級(jí)增加了60 km2和191 km2, 增幅率分別為7.61%和21.39%, 而Ⅰ、Ⅱ和Ⅲ級(jí)分別減少了84 km2、124 km2和42 km2, 降幅分別為31.7%、56.62%和10.91%, Ⅱ級(jí)降幅比較大, 說明研究區(qū)注重生態(tài)環(huán)境的保護(hù), 加強(qiáng)了生態(tài)環(huán)境建設(shè), 植被得到了有效恢復(fù), 植被覆蓋度由低植被覆蓋度向高植被覆蓋度轉(zhuǎn)化的特征。

    圖3 南充市轄區(qū)植被覆蓋度空間分布圖

    Figure 3 Spatial distribution of vegetation coverage in Nanchong Jurisdiction

    圖4 南充市轄區(qū)地形分級(jí)圖

    Figure 4 Classification map of Nanchong Jurisdiction terrain

    表3 各等級(jí)植被覆蓋度占比統(tǒng)計(jì)

    3.2 不同海拔帶植被覆蓋度分布及面積統(tǒng)計(jì)

    將植被覆蓋度圖(圖3 a、b)與海拔分級(jí)圖(圖4a)通過柵格計(jì)算器進(jìn)行計(jì)算, 得到各等級(jí)植被覆蓋度在不同海拔下的分布面積(表4)。從表4中可以看出, 在大于300 m的同一海拔范圍內(nèi), 植被覆蓋度面積隨等級(jí)的升高而增加, 在海拔<300 m時(shí), 植被呈先減少后增加的趨勢(shì)。各個(gè)等級(jí)的植被覆蓋度隨著海拔的升高而持續(xù)減少, 在海拔600 m以上幾乎減少至0。在2002年, 海拔300 m以下以Ⅰ級(jí)植被覆蓋為主, 海拔300 m以上都以Ⅴ級(jí)植被覆蓋為主; 2014年, 各個(gè)海拔帶都以Ⅴ級(jí)植被覆蓋為主。2002—2014年Ⅰ、Ⅱ和Ⅲ級(jí)植被覆蓋度區(qū)域在不同海拔帶上都呈下降趨勢(shì), 其中Ⅳ和Ⅴ級(jí)呈增加趨勢(shì), 在海拔<300 m和300—400 m時(shí), Ⅱ級(jí)下降的比重較大, 分別達(dá)到53%和62%, Ⅴ級(jí)增加40%和18%; 在海拔400—500 m時(shí), Ⅰ級(jí)降幅最大, 達(dá)到57%, Ⅳ和Ⅴ級(jí)的增幅都是16%。500—600 m地帶和>600 m地帶的區(qū)域植被覆蓋度面積變化相對(duì)較小, 說明海拔相對(duì)較低的地區(qū)人類活動(dòng)影響較大, 海拔相對(duì)較高的地方受人類活動(dòng)的影響相對(duì)較小。

    表4 不同海拔帶植被覆蓋度的變化特征

    表5 不同坡度帶植被覆蓋度的變化特征

    3.3 不同坡度帶植被覆蓋度分布及面積統(tǒng)計(jì)

    將植被覆蓋度圖(圖3 a、b)與坡度分級(jí)圖(圖4 b)通過柵格計(jì)算器進(jìn)行計(jì)算, 得到各等級(jí)植被覆蓋度在不同坡度下的分布面積(表5)。根據(jù)表5統(tǒng)計(jì)的數(shù)據(jù)顯示, 研究區(qū)植被覆蓋度具有明顯的坡度分異特征。隨著坡度的增加, 植被覆蓋度呈現(xiàn)先增加后降低的特征。2002年0°—2°的區(qū)域以Ⅰ級(jí)植被覆蓋為主, 其面積比重達(dá)到26.88%; 坡度在2°—6°的區(qū)域以Ⅳ級(jí)植被覆蓋為主, 其面積比重達(dá)到27.89%; 6°以上區(qū)域則以Ⅴ級(jí)植被覆蓋為主, 面積所占比重隨坡度的增加而增加; 在2014年, 各坡度帶中均以Ⅴ級(jí)植被覆蓋為主, 最高比重達(dá)到37.31%。2002—2014年, 研究區(qū)Ⅰ、Ⅱ、Ⅲ級(jí)的植被覆蓋區(qū)域面積總體呈現(xiàn)下降的特征, Ⅳ和Ⅴ級(jí)植被覆蓋呈現(xiàn)增加的特征, 但各個(gè)坡度帶上總的植被覆蓋并沒有發(fā)生較大變化。

    3.4 不同坡向植被覆蓋度分布及面積統(tǒng)計(jì)

    將植被覆蓋度圖(圖2 a、b)與海拔分級(jí)圖(圖3 c)通過柵格計(jì)算器進(jìn)行計(jì)算, 并將坡向按照陰坡(0± 22.5°)、半陰坡(45±22.5°、315°±22.5°、90°±22.5°)、半陽坡(270°±22.5°、135°±22.5°、225°±22.5°)、陽坡(180°±22.5°)分為 4 個(gè)坡度, 統(tǒng)計(jì)各等級(jí)植被覆蓋在不同坡向下的分布面積(表6)。根據(jù)表6 的數(shù)據(jù)顯示, 在各個(gè)坡向上, 隨著植被覆蓋度等級(jí)的升高, 植被覆蓋度呈現(xiàn)先減小后增加的趨勢(shì)。2002年和2014年研究區(qū)域植被覆蓋度在不同坡向上所占的比例總體上都呈現(xiàn)為半陽坡(39%)>半陰坡(38%)>陽坡(11%)>陰坡(12%), 且半陽坡和半陰坡、陽坡和陰坡都相差無幾。各坡向上不同植被覆蓋度等級(jí)所占面積均有不同程度的升降, Ⅰ、Ⅱ和Ⅲ級(jí)植被覆蓋減少, Ⅳ和Ⅴ級(jí)比重增加。

    4 結(jié)論與討論

    4.1 結(jié)論

    本文利用ENVI和ArcGIS軟件, 對(duì)遙感影像進(jìn)行處理, 將影響植被分布的因子進(jìn)行疊加, 得到了研究區(qū)植被覆蓋的空間分布。

    表6 不同坡向植被覆蓋度的變化特征

    (1)研究區(qū)植被覆蓋度總體較好, 具有明顯的空間差異。低植被覆蓋區(qū)域主要集中在轄區(qū)市區(qū), 高植被覆蓋區(qū)域主要集中在嘉陵區(qū)和高坪區(qū)的山區(qū), 其中Ⅴ級(jí)植被覆蓋所占的區(qū)域面積最大, 高達(dá)42.49%以上。2002—2014年, Ⅳ級(jí)和Ⅴ級(jí)增加了60 km2和191 km2, 增幅率分別為7.61%和21.39%, 而Ⅰ、Ⅱ和Ⅲ級(jí)分別減少了84 km2、124 km2和42 km2, 降幅分別為31.7%、56.62%和10.91%。

    (2)研究區(qū)植被覆蓋度因海拔、坡度和坡向等地形因子的影響呈現(xiàn)出分異的特征。其中海拔300 m—400 m和坡度2°—6°的地帶植被覆蓋度相對(duì)較高。植被覆蓋度坡向差異呈現(xiàn)半陽坡>半陰坡>陽坡>陰坡的特征。

    4.2 討論

    研究區(qū)植被覆蓋度總體變化不大, 極低植被覆蓋、低植被覆蓋和中植被覆蓋區(qū)域面積減少, 中高植被覆蓋和高植被覆蓋所占面積增加。1993年, 三區(qū)城市建成區(qū)面積僅20.3 km2, 2005年, 建成區(qū)城區(qū)面積53.74 km2, 2010年, 市轄三區(qū)建成區(qū)面積78 km2, 2013年三區(qū)建成區(qū)面積擴(kuò)張到109 km2, 極低植被覆蓋、低植被覆蓋和中植被覆蓋區(qū)域面積減少, 可能是該時(shí)期城市建設(shè)用地?cái)U(kuò)張的結(jié)果。也可能是極低植被覆蓋、低植被覆蓋和中植被覆蓋中向高植被覆蓋和高植被覆蓋演替的結(jié)果。20世紀(jì)80年代后, 南充市興林業(yè)經(jīng)濟(jì)體制改革, 廣大群眾造林積極性空前高漲, 速豐林得到快速發(fā)展。同時(shí)該時(shí)期南充市加大了生態(tài)環(huán)境保護(hù)的投入, 具體表現(xiàn)為長江防護(hù)林、環(huán)城綠化、退耕還林和天保工程, 研究區(qū)植被得到很好的恢復(fù)和發(fā)育也是中高植被覆蓋和高植被覆蓋所占面積增加的原因。另外, 本文尚有諸多不足之處:

    (1)本研究只對(duì)2002—2014年兩個(gè)節(jié)點(diǎn)上的植被景觀進(jìn)行了分析, 無法掌握研究時(shí)段內(nèi)植被覆蓋度的具體變化過程。未來如有條件, 合理增加影像數(shù)量, 以便更加精確、真實(shí)的分析南充市轄區(qū)的生態(tài)變化情況。

    (2 缺少實(shí)地測(cè)量數(shù)據(jù)NDVIveg與NDVIsoil及其對(duì)應(yīng)的NDVImax與NDVImin值。

    (3)本研究對(duì)南充市轄區(qū)近15a的植被覆蓋時(shí)空變化特征和內(nèi)部轉(zhuǎn)移變化特征進(jìn)行了分析, 但對(duì)其變化的主要影響因素未做深入討論, 未來應(yīng)結(jié)合社會(huì)經(jīng)濟(jì)因子和氣溫、降水等自然因素來分析影像植被覆蓋度變化的驅(qū)動(dòng)因子, 為南充市轄區(qū)生態(tài)系統(tǒng)監(jiān)測(cè)提供依據(jù)。

    [1] 王春梅, 顧行發(fā), 余濤, 等. 農(nóng)田春小麥葉面積指數(shù)和覆蓋度時(shí)空變異性研究[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2014, 45(08): 254–261, 235.

    [2] 袁靜文. 宜昌地區(qū)植被信息提取與變化檢測(cè)研究[J]. 交通科技與經(jīng)濟(jì), 2014, 16(04): 117–120.

    [3] 左偉著. 基于RS、GIS的區(qū)域生態(tài)安全綜合評(píng)價(jià)研究[M]. 北京: 測(cè)繪出版社, 2004..

    [4] 李苗苗, 吳炳方, 顏長珍, 等. 密云水庫上游植被覆蓋度的遙感估算[J]. 資源科學(xué), 2004(04): 153–159.

    [5] 李成范, 蘇迎春, 周廷剛, 等. 基于遙感的重慶市植被覆蓋格局變化研究[J]. 西南師范大學(xué)學(xué)報(bào): 自然科學(xué)版, 2009, 34(01): 123–128.

    [6] 易勝. 基于RS和GIS巖溶地區(qū)植被覆蓋度分析[D]. 廣西師范大學(xué), 2008.

    [7] 陳晉, 陳云浩, 何春陽, 等. 基于土地覆蓋分類的植被覆蓋率估算亞像元模型與應(yīng)用[J]. 遙感學(xué)報(bào), 2001(06): 416– 422+481.

    [8] GRAETZ R.D, PECH R.P, DAVIS A.W. The assessment and monitoring of sparsely vegetated rangelands using calibrated Landsat data[J]. International Journal of Remote Sensing, 1988, 9(7): 1201–1222.

    [9] DYMOND J.R, STEPHENS P.R, NEWSOME P.F, et al. Percent vegetation cover of a degrading rangeland from SPOT[J]. International Journal of Remote Sensing, 1992, 13(11): 1999–2007.

    [10] 程紅芳, 章文波, 陳鋒. 植被覆蓋度的遙感估算方法研究進(jìn)展[J]. 國土資源遙感, 2008(01): 13–18.

    [11] 楊勝天, 劉昌明, 楊志峰, 等. 南水北調(diào)西線調(diào)水工程區(qū)的自然生態(tài)環(huán)境評(píng)價(jià)[J]. 地理學(xué)報(bào), 2002(01): 11–18.

    [12] MURIITHI F.K, YU D, ROBIL A S. Vegetation response to intensive commercial horticulture and environmental changes within watersheds in central highlands, Kenya, using AVHRR NDVI data[J]. Remote Sensing, 2016, 53(1): 1–22.

    [13] ZEWDIE W, CSAPLOVICS E, INOSTROZA L. Monitoring ecosystem dynamics in northwestern Ethiopia using NDVI and climate variables to assess long term trends in dryland vegetation variability[J]. Applied Geography, 2017, 79: 167–178.

    [14] CHAITHANVA V.V, BINOY B.V, VINOD T.R. Estimation of the Relationship between Urban Vegetation and Land Surface Temperature of Calicut City and Suburbs, Kerala, India using GIS and Remote Sensing data[J]. International Journal of Advanced Remote Sensing and GIS, 2017, 6(1): 2088–2096.

    [15] 呂長春, 王忠武, 錢少猛. 混合像元分解模型綜述[J]. 遙感信息, 2003(03): 55–58, 60.

    [16] 田義超, 梁銘忠. 北部灣沿海地區(qū)植被覆蓋對(duì)氣溫和降水的旬響應(yīng)特征[J]. 自然資源學(xué)報(bào), 2016, 31(03): 488– 502.

    [17] LI Jingzhong, LIU Yongmei, CAO Mingming, et al. Space-Time Characteristics of Vegetation Cover and Distribution: Case of the Henan Province in China[J]. Sustainability, 2015, 7(9): 11967–11979.

    [18] WANG Haiting, ZHANG Yuanzhi, TSOU Jin Yeu, et al. Surface Urban Heat Island Analysis of Shanghai (China) Based on the Change of Land Use and Land Cover [J]. Sustainability, 2017, 9(09): 1538.

    [19] 李苗苗. 植被覆蓋度的遙感估算方法研究[D]. 北京:中國科學(xué)院研究生院(遙感應(yīng)用研究所), 2003.

    [20] 鄧書斌. ENVI遙感圖像處理方法[M]. 北京: 科學(xué)出版社, 2017.

    [21] 陳君穎, 田慶久. 高分辨率遙感植被分類研究[J]. 遙感學(xué)報(bào), 2007(02): 221–227.

    [22] 岳瑋, 劉慧明, 孫國鈞. 基于遙感和GIS技術(shù)的祖厲河流域植被覆蓋動(dòng)態(tài)變化監(jiān)測(cè)[J]. 蘭州大學(xué)學(xué)報(bào)(自然科學(xué)版), 2009, 45(S1): 6–11, 18.

    [23] 趙翠娥, 丁文榮. 基于ENVI和GIS技術(shù)的龍川江流域植被覆蓋度動(dòng)態(tài)監(jiān)測(cè)[J]. 林業(yè)調(diào)查規(guī)劃, 2013, 38(05): 14-18+44.

    [24] 肖驍, 李京忠, 韓彬, 等. 東北老工業(yè)區(qū)植被覆蓋度時(shí)空特征及城市化關(guān)聯(lián)分析[J]. 生態(tài)科學(xué), 2017, 36(06): 71– 77.

    [25] 煙貫發(fā), 萬魯河, 溫智虹, 等. 基于RS和DEM的長白山天池植被分布的坡度坡向分析[J]. 測(cè)繪通報(bào), 2012(S1): 233–236, 239.

    Characteristics of vegetation coverage changes in Nanchong Jurisdiction in the past fifteen years

    WANG Lei1, WANG Jie1, FU Lin2, LIN Yebin1,3*

    1. China West Normal University, Nanchong, 637009, China 2. Southwest Jiaotong University, Chengdu, 611756, China 3. Research Center for Regional Environmental Evolution and Conservation, China West Normal University,Nanchong, 637009, China

    Vegetation coverage is a combination of natural conditions and socio-economic activities. Studying the temporal and spatial variation of vegetation coverage of a region is significant for understanding the social development in the region, formulating a more rational ecological and environmental policy and realizing sustainable economic development. In this paper, we chose Landsat 7 ETM in 2002, Landsat 8 images in 2014 and DEM as data sources, and used ENVI and ArcGIS software to process them. The NDVI method was used to estimate the vegetation coverage of Nanchong Jurisdiction and divided into five different Cover levels. According to the DEM data, altitude, slope and aspect factors were extracted. Then, the vegetation coverage maps were overlaid with the altitude, slope and aspect maps to quantitatively analyze the spatial distribution characteristics of vegetation in Nanchong Jurisdiction and the changes. Results showed that the overall vegetation coverage in Nanchong Jurisdiction was almost unchanged from 2002 to 2014, of which the vegetation coverage of Grade V (fc≥0.8) was the dominant factor, accounting for more than 42.49%. IV andⅤwere increased by 60 and 191 km2, respectively, with increases of 7.61% and 21.39%.While Levels I, II and III were decreased by 84 km2, 124 km2and 42 km2respectively , with decreases of 31.7%, 56.62% and 10.91%. The vegetation coverage was mainly distributed at elevations of 300-400 m, on slopes of 2°-6°. And the differences of vegetation coverage in slope direction show the characteristics of semi-sunny slope > semi-shady slope > sunny slope > shady slope.

    vegetation coverage; NDVI;spatial distribution;Nanchong Jurisdiction

    10.14108/j.cnki.1008-8873.2019.01.021

    K903

    A

    1008-8873(2019)01-159-09

    2018-01-18;

    2018-02-26

    四川省教育廳自然科學(xué)重點(diǎn)項(xiàng)目(17ZA0376); 西華師范大學(xué)英才基金項(xiàng)目(17YC113)

    王磊(1992—), 男, 四川宜賓人, 碩士研究生, 主要從事土地利用遙感監(jiān)測(cè)研究, E-mail:yibinwanglei@163.com

    林葉彬(1981—), 男, 博士, 副教授, 主要研究方向?yàn)橘Y源利用與環(huán)境保護(hù), E-mail:forrest@cwnu.edu.cn

    王磊, 王杰, 付林, 等. 南充市轄區(qū)近15年植被覆蓋度變化特征[J]. 生態(tài)科學(xué), 2019, 38(1): 159-167.

    WANG Lei, WANG Jie, FU Lin, et al. Characteristics of vegetation coverage changes in Nanchong Jurisdiction in the past fifteen years[J]. Ecological Science, 2019, 38(1): 159-167.

    猜你喜歡
    市轄區(qū)南充坡向
    同構(gòu)視閾下的異構(gòu)治理:市轄區(qū)體制的優(yōu)化進(jìn)路*
    四川省南充蠶具研究有限公司
    DEM地表坡向變率的向量幾何計(jì)算法
    堅(jiān)持法德并舉 助力治蜀興川
    四川南充上智農(nóng)業(yè)機(jī)械設(shè)備有限公司
    南充市委書記為什么夸玉皇觀村
    我國市轄區(qū)行政體制改革難點(diǎn)探析
    法制博覽(2016年19期)2016-02-01 05:11:07
    青藏高原東緣高寒草甸坡向梯度上植物光合生理特征研究
    土壤與作物(2015年3期)2015-12-08 00:46:55
    中國城市市轄區(qū)設(shè)置和發(fā)展評(píng)價(jià)研究
    地級(jí)市及市轄區(qū)職責(zé)劃分研究
    两性午夜刺激爽爽歪歪视频在线观看| 欧美日本视频| 欧美日韩一区二区视频在线观看视频在线 | 亚洲av熟女| 国产亚洲精品久久久com| 99热这里只有精品一区| 久久久久久久久久成人| 欧美日本亚洲视频在线播放| 国产国拍精品亚洲av在线观看| 丝袜美腿在线中文| 亚洲国产精品久久男人天堂| 亚洲专区国产一区二区| 国产精品一及| 性色avwww在线观看| 欧美成人免费av一区二区三区| 欧美日本视频| 色在线成人网| 一夜夜www| 乱系列少妇在线播放| 精品一区二区三区av网在线观看| 亚洲一区二区三区色噜噜| 男女下面进入的视频免费午夜| 一个人免费在线观看电影| 少妇熟女欧美另类| 午夜日韩欧美国产| 国产av在哪里看| 九九热线精品视视频播放| 亚洲国产精品成人久久小说 | 熟女电影av网| 深夜精品福利| 国产欧美日韩一区二区精品| 国产精品野战在线观看| 22中文网久久字幕| 亚洲熟妇熟女久久| 51国产日韩欧美| 国产欧美日韩精品一区二区| 寂寞人妻少妇视频99o| 久久中文看片网| 一区二区三区免费毛片| 精品久久久久久久久久免费视频| 亚洲国产高清在线一区二区三| 菩萨蛮人人尽说江南好唐韦庄 | 一区福利在线观看| 中国美女看黄片| 内地一区二区视频在线| 你懂的网址亚洲精品在线观看 | 欧美日韩在线观看h| 亚洲性夜色夜夜综合| 久久精品夜色国产| 一个人看视频在线观看www免费| 国产av在哪里看| 91久久精品国产一区二区三区| 午夜爱爱视频在线播放| 国产精品人妻久久久久久| 久久精品人妻少妇| 欧美一区二区国产精品久久精品| 国产精品一二三区在线看| 日韩在线高清观看一区二区三区| 97人妻精品一区二区三区麻豆| 免费看美女性在线毛片视频| 欧美另类亚洲清纯唯美| 女同久久另类99精品国产91| АⅤ资源中文在线天堂| 国产黄色视频一区二区在线观看 | 国产aⅴ精品一区二区三区波| .国产精品久久| 日韩欧美一区二区三区在线观看| 亚洲欧美清纯卡通| 国产色婷婷99| 一进一出抽搐gif免费好疼| 波多野结衣高清作品| 久久久色成人| 麻豆一二三区av精品| 91午夜精品亚洲一区二区三区| 亚洲无线在线观看| 晚上一个人看的免费电影| 成人高潮视频无遮挡免费网站| 日本色播在线视频| 一边摸一边抽搐一进一小说| 99热精品在线国产| 亚洲精品一区av在线观看| 18禁在线播放成人免费| 成熟少妇高潮喷水视频| 久久精品91蜜桃| 99热只有精品国产| 国产中年淑女户外野战色| 人人妻人人看人人澡| 日韩一本色道免费dvd| 成人特级av手机在线观看| 精品无人区乱码1区二区| 91在线观看av| 日本 av在线| 亚洲五月天丁香| 在线播放国产精品三级| 一区二区三区高清视频在线| 欧美性感艳星| 日本三级黄在线观看| 成人亚洲欧美一区二区av| 久久精品国产亚洲av香蕉五月| 岛国在线免费视频观看| 乱码一卡2卡4卡精品| 亚洲第一区二区三区不卡| 亚洲国产欧洲综合997久久,| 亚洲一区高清亚洲精品| 亚洲精品一卡2卡三卡4卡5卡| 香蕉av资源在线| 亚洲欧美日韩高清专用| 最近最新中文字幕大全电影3| 免费观看在线日韩| 哪里可以看免费的av片| 成人鲁丝片一二三区免费| 免费看光身美女| 久久亚洲精品不卡| 久久午夜亚洲精品久久| 午夜免费男女啪啪视频观看 | 九九爱精品视频在线观看| 一级黄色大片毛片| 亚洲18禁久久av| 国产美女午夜福利| 少妇熟女aⅴ在线视频| 男女做爰动态图高潮gif福利片| 免费人成视频x8x8入口观看| 国产亚洲欧美98| 麻豆一二三区av精品| 村上凉子中文字幕在线| 亚洲高清免费不卡视频| 成人欧美大片| 精品久久久久久久久亚洲| 国产成人a∨麻豆精品| 最近最新中文字幕大全电影3| 99久久久亚洲精品蜜臀av| 噜噜噜噜噜久久久久久91| 国产黄色视频一区二区在线观看 | 欧美最黄视频在线播放免费| 国产真实伦视频高清在线观看| 欧美激情在线99| 亚洲av电影不卡..在线观看| 国产免费一级a男人的天堂| 亚州av有码| 亚洲欧美中文字幕日韩二区| 成人高潮视频无遮挡免费网站| 精品一区二区免费观看| 国产美女午夜福利| 亚洲经典国产精华液单| 亚洲精品国产成人久久av| 中出人妻视频一区二区| 身体一侧抽搐| 日韩成人av中文字幕在线观看 | 国产av不卡久久| 99热6这里只有精品| 欧美激情国产日韩精品一区| 亚洲美女视频黄频| 国产中年淑女户外野战色| 亚洲国产精品合色在线| 亚洲熟妇中文字幕五十中出| 成人鲁丝片一二三区免费| .国产精品久久| 少妇的逼好多水| 国产日本99.免费观看| 日韩,欧美,国产一区二区三区 | 日本色播在线视频| 精品欧美国产一区二区三| 欧美日本视频| 亚洲自偷自拍三级| 女的被弄到高潮叫床怎么办| 国产午夜福利久久久久久| 我要看日韩黄色一级片| 99热这里只有是精品50| 婷婷精品国产亚洲av在线| 久久6这里有精品| 又黄又爽又刺激的免费视频.| 中文字幕av成人在线电影| 69人妻影院| 老女人水多毛片| 男女啪啪激烈高潮av片| 97热精品久久久久久| 大又大粗又爽又黄少妇毛片口| 自拍偷自拍亚洲精品老妇| 天天躁日日操中文字幕| 国产激情偷乱视频一区二区| 亚洲成人中文字幕在线播放| 成人性生交大片免费视频hd| 搞女人的毛片| АⅤ资源中文在线天堂| 欧美国产日韩亚洲一区| 国产精品综合久久久久久久免费| 嫩草影院入口| 欧美人与善性xxx| 丰满人妻一区二区三区视频av| 欧美xxxx性猛交bbbb| 久久鲁丝午夜福利片| 国产一区亚洲一区在线观看| 人妻久久中文字幕网| 国产极品精品免费视频能看的| 淫秽高清视频在线观看| 国产av一区在线观看免费| 久久精品人妻少妇| 99久久精品国产国产毛片| 国产亚洲精品av在线| 国产成人91sexporn| 国产美女午夜福利| h日本视频在线播放| 国产一区二区在线av高清观看| 99久国产av精品| 亚洲欧美精品自产自拍| 深爱激情五月婷婷| 国产欧美日韩精品一区二区| 99久久中文字幕三级久久日本| 麻豆乱淫一区二区| 国产免费一级a男人的天堂| 亚洲av第一区精品v没综合| 国产精品一区二区免费欧美| 成人国产麻豆网| 国产精品美女特级片免费视频播放器| 99视频精品全部免费 在线| 嫩草影院精品99| 国产精品永久免费网站| 国产一区二区三区av在线 | 精品一区二区三区av网在线观看| 男女下面进入的视频免费午夜| 熟女人妻精品中文字幕| 伊人久久精品亚洲午夜| 18禁黄网站禁片免费观看直播| 高清午夜精品一区二区三区 | 性欧美人与动物交配| 中文字幕av在线有码专区| 成人毛片a级毛片在线播放| 成人国产麻豆网| 欧美一区二区亚洲| 在线免费十八禁| 联通29元200g的流量卡| 黄色日韩在线| 亚洲乱码一区二区免费版| 欧美性猛交╳xxx乱大交人| 99热全是精品| 国产免费男女视频| www日本黄色视频网| 欧美激情在线99| 久久久久精品国产欧美久久久| 国产午夜福利久久久久久| 波多野结衣高清作品| 精品午夜福利视频在线观看一区| 日韩一区二区视频免费看| 久久精品人妻少妇| 免费观看人在逋| 欧美成人免费av一区二区三区| 亚洲最大成人中文| 欧美bdsm另类| 国产精品一二三区在线看| 国产极品精品免费视频能看的| 日日干狠狠操夜夜爽| 日韩一区二区视频免费看| www.色视频.com| 久久人人爽人人片av| 国产精品亚洲一级av第二区| 日韩av在线大香蕉| 搞女人的毛片| 偷拍熟女少妇极品色| 亚洲三级黄色毛片| 国内揄拍国产精品人妻在线| 一级av片app| 日韩欧美 国产精品| 国产成人a∨麻豆精品| 尾随美女入室| 国产一区二区激情短视频| 高清毛片免费看| 午夜影院日韩av| av女优亚洲男人天堂| 亚洲av五月六月丁香网| 国产高清有码在线观看视频| 97超碰精品成人国产| 国产成年人精品一区二区| 极品教师在线视频| 亚洲自拍偷在线| 国产人妻一区二区三区在| 亚洲av二区三区四区| 18禁在线播放成人免费| 亚洲成av人片在线播放无| 黑人高潮一二区| 成年av动漫网址| 午夜福利在线在线| 亚洲国产精品成人久久小说 | 亚洲综合色惰| 精品99又大又爽又粗少妇毛片| 亚洲精品亚洲一区二区| 亚洲最大成人手机在线| 久久久a久久爽久久v久久| eeuss影院久久| 三级毛片av免费| 日本-黄色视频高清免费观看| 国产高清有码在线观看视频| 国产探花极品一区二区| 亚洲综合色惰| 直男gayav资源| 国产精品综合久久久久久久免费| 日本爱情动作片www.在线观看 | 69av精品久久久久久| 在线播放无遮挡| 变态另类丝袜制服| 成人美女网站在线观看视频| 在线免费观看的www视频| 中文字幕免费在线视频6| 欧美又色又爽又黄视频| 精品无人区乱码1区二区| 亚洲成人久久性| 最后的刺客免费高清国语| 国产精品一区二区三区四区免费观看 | 日本三级黄在线观看| 成人三级黄色视频| 国产69精品久久久久777片| 国产爱豆传媒在线观看| 国产精品一区www在线观看| 你懂的网址亚洲精品在线观看 | 激情 狠狠 欧美| 午夜久久久久精精品| 男人舔女人下体高潮全视频| 麻豆国产av国片精品| 欧美日韩一区二区视频在线观看视频在线 | 欧美一区二区国产精品久久精品| 国产精品野战在线观看| 国产精品日韩av在线免费观看| 老女人水多毛片| 日韩制服骚丝袜av| 精品久久久久久久久久久久久| 黄片wwwwww| 欧美中文日本在线观看视频| 成人三级黄色视频| 特大巨黑吊av在线直播| 色噜噜av男人的天堂激情| 又爽又黄无遮挡网站| av中文乱码字幕在线| 国产精品一区www在线观看| 免费观看在线日韩| 午夜精品在线福利| 1000部很黄的大片| 人妻少妇偷人精品九色| 午夜福利高清视频| 能在线免费观看的黄片| 欧美+亚洲+日韩+国产| 亚洲国产色片| 日韩欧美 国产精品| 老司机影院成人| 熟女电影av网| 国内精品宾馆在线| 国产免费一级a男人的天堂| 给我免费播放毛片高清在线观看| 国产一级毛片七仙女欲春2| 99riav亚洲国产免费| 国产精品1区2区在线观看.| 亚洲国产精品成人综合色| 一个人看视频在线观看www免费| av在线天堂中文字幕| 免费av观看视频| 国产一区二区在线av高清观看| 青春草视频在线免费观看| 日韩三级伦理在线观看| 精品久久久噜噜| 97超视频在线观看视频| 日韩大尺度精品在线看网址| 国产精品三级大全| 中文字幕av在线有码专区| 国国产精品蜜臀av免费| 亚洲精华国产精华液的使用体验 | 少妇熟女欧美另类| 亚洲乱码一区二区免费版| 色视频www国产| 人人妻人人澡欧美一区二区| 老司机午夜福利在线观看视频| 国产 一区精品| 亚洲一区高清亚洲精品| 亚洲精品成人久久久久久| 性欧美人与动物交配| 亚洲自拍偷在线| 又爽又黄无遮挡网站| 深夜a级毛片| 老司机午夜福利在线观看视频| 国产真实伦视频高清在线观看| 国产单亲对白刺激| eeuss影院久久| 欧美极品一区二区三区四区| 菩萨蛮人人尽说江南好唐韦庄 | 久久久久九九精品影院| 国产极品精品免费视频能看的| 国产亚洲精品久久久com| 永久网站在线| 国产美女午夜福利| 欧美日本亚洲视频在线播放| 99久国产av精品| 少妇被粗大猛烈的视频| 久久精品国产鲁丝片午夜精品| 丰满的人妻完整版| 久久人妻av系列| 国产精品国产三级国产av玫瑰| 国产高清三级在线| 免费高清视频大片| 久久久欧美国产精品| 中国国产av一级| 欧美+日韩+精品| а√天堂www在线а√下载| 大香蕉久久网| 在线观看午夜福利视频| 国产成人精品久久久久久| 一本一本综合久久| 午夜福利高清视频| 麻豆av噜噜一区二区三区| 午夜影院日韩av| 99久国产av精品| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品乱码久久久v下载方式| 久久这里只有精品中国| 国产精品一区二区免费欧美| 亚洲国产欧美人成| 欧美+日韩+精品| 女人十人毛片免费观看3o分钟| 看十八女毛片水多多多| 午夜福利视频1000在线观看| 最近在线观看免费完整版| 国产高清视频在线播放一区| 韩国av在线不卡| 成人三级黄色视频| 国产av一区在线观看免费| 久久久久久久久久黄片| 国产伦精品一区二区三区四那| 国产精品爽爽va在线观看网站| 国产伦在线观看视频一区| av视频在线观看入口| 人人妻人人澡人人爽人人夜夜 | 国产毛片a区久久久久| 欧美成人a在线观看| 女同久久另类99精品国产91| 亚洲国产精品成人久久小说 | 国产精品女同一区二区软件| 免费无遮挡裸体视频| 午夜免费激情av| 国产免费一级a男人的天堂| 少妇人妻精品综合一区二区 | 日韩欧美一区二区三区在线观看| 日韩大尺度精品在线看网址| 乱人视频在线观看| 国产白丝娇喘喷水9色精品| 午夜老司机福利剧场| av.在线天堂| 亚洲欧美中文字幕日韩二区| 亚洲美女黄片视频| 直男gayav资源| 天堂动漫精品| 97超碰精品成人国产| 69人妻影院| 亚洲专区国产一区二区| av女优亚洲男人天堂| av天堂在线播放| 如何舔出高潮| 在线免费观看的www视频| 美女被艹到高潮喷水动态| 18禁在线播放成人免费| 欧美最新免费一区二区三区| 欧美色视频一区免费| 久久精品人妻少妇| 人妻久久中文字幕网| av免费在线看不卡| 国产欧美日韩精品一区二区| 一夜夜www| 欧美日韩综合久久久久久| 别揉我奶头~嗯~啊~动态视频| 2021天堂中文幕一二区在线观| 最后的刺客免费高清国语| 久久精品91蜜桃| 狂野欧美白嫩少妇大欣赏| 亚洲精品日韩av片在线观看| 又爽又黄a免费视频| 国产欧美日韩精品亚洲av| 91在线观看av| 亚洲精品影视一区二区三区av| 日本黄大片高清| 日本撒尿小便嘘嘘汇集6| 午夜a级毛片| 亚洲在线观看片| or卡值多少钱| 一个人看的www免费观看视频| 国产亚洲精品久久久com| 国产精品久久久久久亚洲av鲁大| 晚上一个人看的免费电影| 亚洲一级一片aⅴ在线观看| 波多野结衣巨乳人妻| 日本a在线网址| 国产精品久久久久久久电影| 99热全是精品| 精品一区二区免费观看| 午夜视频国产福利| 日日摸夜夜添夜夜爱| 亚洲欧美精品自产自拍| 欧美丝袜亚洲另类| 蜜臀久久99精品久久宅男| 成年女人毛片免费观看观看9| 男人舔女人下体高潮全视频| 真实男女啪啪啪动态图| 亚洲内射少妇av| 看片在线看免费视频| av在线天堂中文字幕| 人妻制服诱惑在线中文字幕| 久久6这里有精品| 国产精品亚洲一级av第二区| 日韩高清综合在线| 一区二区三区高清视频在线| 亚洲成a人片在线一区二区| 日韩中字成人| 亚洲国产精品久久男人天堂| 天堂动漫精品| 尾随美女入室| 成人无遮挡网站| 亚洲专区国产一区二区| 熟女人妻精品中文字幕| 国产真实伦视频高清在线观看| 亚洲久久久久久中文字幕| 亚洲av免费在线观看| 国产伦精品一区二区三区四那| 我要看日韩黄色一级片| 亚洲av五月六月丁香网| 99久国产av精品国产电影| 男女啪啪激烈高潮av片| 99久国产av精品| 91麻豆精品激情在线观看国产| 天堂av国产一区二区熟女人妻| 免费大片18禁| 一个人看的www免费观看视频| 亚洲av免费高清在线观看| 男人舔奶头视频| 晚上一个人看的免费电影| 国产成人91sexporn| 99热精品在线国产| 免费观看人在逋| 在线观看66精品国产| 老司机影院成人| 亚洲精品一区av在线观看| 欧美极品一区二区三区四区| 黑人高潮一二区| 亚洲人成网站在线观看播放| 岛国在线免费视频观看| 日本-黄色视频高清免费观看| 国产成人aa在线观看| 欧美bdsm另类| 十八禁国产超污无遮挡网站| 久久久久久久午夜电影| 欧美激情久久久久久爽电影| 精品一区二区三区视频在线| 国产午夜精品久久久久久一区二区三区 | 久久精品国产亚洲av天美| 亚洲精品粉嫩美女一区| 国产精品亚洲美女久久久| 久久午夜福利片| 精品熟女少妇av免费看| 中文字幕精品亚洲无线码一区| 国产不卡一卡二| 国产v大片淫在线免费观看| 亚洲一级一片aⅴ在线观看| 久久久欧美国产精品| 99久久无色码亚洲精品果冻| 欧美区成人在线视频| 少妇人妻精品综合一区二区 | 91午夜精品亚洲一区二区三区| 精品99又大又爽又粗少妇毛片| 亚洲18禁久久av| 免费看光身美女| 99久久精品国产国产毛片| 三级国产精品欧美在线观看| 国产单亲对白刺激| 在线免费观看的www视频| 尾随美女入室| 老司机影院成人| 日韩精品有码人妻一区| 久久久久久久久久成人| 亚洲av第一区精品v没综合| 亚洲一区二区三区色噜噜| 欧美色欧美亚洲另类二区| 久久这里只有精品中国| 国产av麻豆久久久久久久| 成人午夜高清在线视频| 国产成人精品久久久久久| 91精品国产九色| 日本精品一区二区三区蜜桃| 亚洲性久久影院| 18禁裸乳无遮挡免费网站照片| 亚洲精品亚洲一区二区| 亚洲国产精品久久男人天堂| 国产一区二区亚洲精品在线观看| 亚洲最大成人中文| 国产人妻一区二区三区在| 看非洲黑人一级黄片| 国产激情偷乱视频一区二区| 性色avwww在线观看| 99国产极品粉嫩在线观看| 午夜精品国产一区二区电影 | 寂寞人妻少妇视频99o| 精品99又大又爽又粗少妇毛片| a级毛片a级免费在线| 欧美三级亚洲精品| 校园春色视频在线观看| 成人特级av手机在线观看| 高清午夜精品一区二区三区 | 亚洲av五月六月丁香网| 在线看三级毛片| 十八禁网站免费在线| 亚洲精品一卡2卡三卡4卡5卡| 最新在线观看一区二区三区| 九色成人免费人妻av| 精品不卡国产一区二区三区| 日本黄色片子视频| 国产成人福利小说| 国产精品亚洲一级av第二区| 国产精品久久久久久久久免| 成人特级av手机在线观看| 久久精品91蜜桃| 99热这里只有精品一区| 18禁在线播放成人免费| ponron亚洲|