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

    瑪納斯河流域近14年多尺度粒度變化對土地利用景觀格局的影響

    2017-09-12 13:22:53唐湘玲
    水土保持研究 2017年2期
    關(guān)鍵詞:瑪納斯格局粒度

    韓 蓉, 唐湘玲

    (石河子大學(xué) 理學(xué)院 地理系, 新疆 石河子 832003)

    瑪納斯河流域近14年多尺度粒度變化對土地利用景觀格局的影響

    韓 蓉, 唐湘玲

    (石河子大學(xué) 理學(xué)院 地理系, 新疆 石河子 832003)

    以2000年、2014年瑪納斯河流域30 m空間分辨率的遙感影像為基礎(chǔ)數(shù)據(jù)源,借助GIS和ENVI軟件,制作瑪納斯河流域土地利用覆被圖,選取景觀指數(shù),在類型和整體景觀的水平上,探討了瑪納斯河流域景觀格局指數(shù)隨粒度變化的基本規(guī)律,并分析了各景觀格局指數(shù)間的相關(guān)性。結(jié)果表明:100~3 000 m粒度范圍內(nèi),流域景觀格局指數(shù)的“臨界粒度”“粗?;爆F(xiàn)象明顯,總體而言,在400~800 m,2 000~3 000 m分辨率下瑪納斯河流域的景觀格局指數(shù)的“臨界粒度”為400 m和2 000 m,適宜計(jì)算的粒度范圍為400~2 000 m,所以在利用遙感影像對瑪納斯河流域景觀進(jìn)行分析研究時(shí),需進(jìn)行一定的粒度轉(zhuǎn)換及注意粒度的影響?,敿{斯河流域景觀形態(tài)具有分形特征,各類景觀斑塊的分維數(shù)對粒度變化的響應(yīng)不同,分維數(shù)隨粒度的增大呈逐漸下降趨勢,表明景觀類型邊界趨于簡單化。相關(guān)性分析定量反映了所選景觀指數(shù)受粒度變化影響的相關(guān)性程度,可為瑪河流域景觀格局變化研究提供參考。2000年、2014年間瑪河流域景觀多樣性由復(fù)雜逐漸變?yōu)楹唵?,多樣性指?shù)逐步減小。

    土地利用; 景觀格局; 多尺度; 粒度變化; 相關(guān)分析; 瑪納斯河流域

    在景觀生態(tài)學(xué)中,把尺度定義為對某一研究對象或現(xiàn)象在空間或時(shí)間上的量度,往往以粒度大小和幅度變化來表示[1]。在地學(xué)和生態(tài)學(xué)研究中,研究尺度效應(yīng)需要解答的問題主要在于:(1) 在哪種尺度上,可以正確地表達(dá)特定的地理現(xiàn)象?(2) 如何有效地將數(shù)據(jù)和信息從一種尺度轉(zhuǎn)換為另一種尺度[2]?(3) 原始數(shù)據(jù)和信息經(jīng)過尺度轉(zhuǎn)換后,會出現(xiàn)何種信息的損失或效應(yīng)?即不同尺度的數(shù)據(jù)反映相同的地物和現(xiàn)象時(shí)的差異如何[3]?多尺度空間分析是進(jìn)行尺度效應(yīng)分析的基礎(chǔ),是發(fā)現(xiàn)和識別景觀等級結(jié)構(gòu)和特征尺度的主要方法。從不同的空間尺度研究土地利用空間數(shù)據(jù)對土地利用景觀格局的表達(dá),有利于從不同的層次,認(rèn)識土地利用與覆蓋變化的機(jī)理與過程[4]。目前,格局與過程的時(shí)空尺度變化是景觀生態(tài)學(xué)的研究熱點(diǎn),尺度效應(yīng)分析受到格外重視和發(fā)展[5],其文獻(xiàn)在生態(tài)學(xué)的相關(guān)研究中占相當(dāng)重大的比重[6-7]。景觀指數(shù)在景觀生態(tài)學(xué)中應(yīng)用很廣泛,但景觀指數(shù)對不同景觀特征和分析尺度的響應(yīng)及其生態(tài)學(xué)意義尚缺乏明確的闡述[8]。景觀格局作為景觀異質(zhì)性在空間上的綜合表現(xiàn),具有很強(qiáng)的尺度依賴性,分析同一景觀格局在不同尺度下的特征及變化規(guī)律,明確區(qū)域景觀格局特征的尺度效應(yīng),可以更加準(zhǔn)確地反映土地景觀格局的真實(shí)狀況[9]。只有在連續(xù)的尺度序列上對景觀格局進(jìn)行考察和研究,才能把握其內(nèi)在的演變規(guī)律[10-11],所以選擇適宜的研究尺度顯得極為重要。尺度選擇不當(dāng),就不能科學(xué)地揭示研究對象的本質(zhì)。尺度選擇過大,往往導(dǎo)致大量細(xì)節(jié)被忽略;尺度選擇過小,就會陷入局部,容易忽略總體規(guī)律[10-12]。許多學(xué)者對基于1∶50 000以下小比例尺的景觀指數(shù)粒度效應(yīng)進(jìn)行了研究[13-20],對不同比例尺的矢量數(shù)據(jù)如何選擇適宜粒度的研究也有涉及[21-23]。土地調(diào)查數(shù)據(jù)作為開展土地資源合理管理和優(yōu)化配置的基礎(chǔ),也是開展景觀生態(tài)分析數(shù)據(jù)的重要來源,為景觀格局的分析與研究提供了有利條件[24-26],但當(dāng)前涉及基于此數(shù)據(jù)開展的景觀格局對粒度變化的影響及其機(jī)理研究較少。而且目前大多數(shù)都是針對單一尺度上的研究,多尺度空間粒度研究較少。

    本文以瑪納斯河流域30 m空間分辨率的遙感影像為基礎(chǔ)數(shù)據(jù)源,通過GIS和ENVI軟件技術(shù)支持,運(yùn)用景觀生態(tài)學(xué)相關(guān)方法,從類型與景觀2個(gè)水平上選取多個(gè)經(jīng)典景觀格局指數(shù)選擇重要的景觀指數(shù)[27],主要分析土地利用景觀指數(shù)的差異和各土地利用景觀類型格局變化趨勢特征機(jī)理,并對所選景觀指數(shù)在粒度上的差異進(jìn)行相關(guān)性分析??偨Y(jié)相同景觀在不同粒度下的景觀特征差異,使得景觀格局指數(shù)在測度過程中更具說服力,對不同空間粒度下景觀格局的差異做出正確解釋,為土地利用現(xiàn)狀數(shù)據(jù)開展景觀格局分析,指導(dǎo)土地資源利用規(guī)劃和保護(hù)提供理論依據(jù)。

    1 研究區(qū)概況

    瑪納斯和流域位于新疆準(zhǔn)噶爾盆地南部,西南部分別相鄰于準(zhǔn)噶爾盆地邊緣和天山山脈,地理位置為85°01′—86°32′E,43°27′—45°21′N,最高海拔為5 242.5 m,平原區(qū)海拔為300~500 m。東西198.7 km,南北260.8 km,流域面積3.1萬km2。地勢由東南向西北傾斜,海拔最高5 259 m,最低168 m。由南向北依次分為山地、山前平原和沙漠三大地貌類型區(qū)(比例約為2.83∶1∶1.07)[28],流域?qū)賰?nèi)陸干旱區(qū),夏季炎熱干燥,冬季寒冷多風(fēng),年均氣溫6.8℃,年降水量110~200 mm[29],年蒸發(fā)量約為1 500~2 000 mm,瑪納斯河流域人口約為1.3×105人。景觀格局十分復(fù)雜,即以草地、林地、耕地、荒漠為主,是流域自然狀況的體現(xiàn),瑪納斯河流域?qū)俅箨懶愿珊禋夂颍瑲夂蚨湎臒?,光、熱、風(fēng)、水資源豐富。

    2 資料及方法

    2.1 數(shù)據(jù)來源與處理

    以瑪納斯河小流域全色波段遙感2000年的Landsat ETM+和2014年的Landsat OLI數(shù)據(jù);分辨率為30 m的ASTER GDEM V2版的DEM數(shù)據(jù)(數(shù)據(jù)來源于中國科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國際科學(xué)數(shù)據(jù)鏡像網(wǎng)站(http:∥www.gscloud.cn))。選取無云或少云的遙感影像和DEM在ENVI 5.1的支持下,進(jìn)行鑲嵌和裁剪,并將坐標(biāo)系統(tǒng)統(tǒng)一Albers_Conic_Equal_Area投影,基準(zhǔn)面為WGS_1984,中央經(jīng)線為105°E。

    本文遙感影像分類在ENVI 5.1中完成。首先數(shù)據(jù)預(yù)處理,分別進(jìn)行2000年和2014年瑪納斯河流域遙感影像圖的NDVI和ISODATA非監(jiān)督分類計(jì)算,在進(jìn)行ISODATA分類時(shí),最大分類數(shù)量設(shè)置為10,迭代次數(shù)設(shè)置為10。將獲得的瑪納斯河流域的DEM,NDVI,ISODATA和近紅外、紅、綠波段進(jìn)行波段組合,得到一個(gè)6波段的文件。然后,選取訓(xùn)練樣本。根據(jù)瑪納斯河流域的實(shí)際土地利用情況和遙感影像的紋理信息和光譜特征,本文將研究區(qū)土地分類為:(1) 草地包括草甸、綠地;(2) 林地包括高覆蓋灌木林、常綠灌木林;(3)水體包括湖泊、河流;(4) 建設(shè)用地;(5)未利用地包括裸地和荒漠、鹽堿地等;(6) 耕地。在兩期影像中,分別選取一定數(shù)量的6類土地利用的訓(xùn)練樣本,并保證訓(xùn)練樣本的可分離度在1.8以上。最后,獲取規(guī)則和應(yīng)用規(guī)則分類。根據(jù)選擇的訓(xùn)練樣本,基于CART算法獲取決策樹規(guī)則,進(jìn)行決策樹分類,最終提取出瑪納斯河流域土地利用信息。本文中,使用混淆矩陣方法,檢驗(yàn)遙感影像分類的精度。經(jīng)過精度評價(jià),研究區(qū)2000年、2014年分類圖總體分類精度為91.6%,88.5%,可以滿足本研究的精度要求。最后,在ArcGIS 10.1支持下,通過格式轉(zhuǎn)換、空間疊加等處理,進(jìn)行數(shù)據(jù)統(tǒng)計(jì)分析。

    2.2 指標(biāo)計(jì)算與技術(shù)路線

    基于ArcGIS 10.1軟件平臺,利用Spatial Analyst模塊將瑪納斯河流域土地利用類型圖(2000年和2014年)的矢量數(shù)據(jù)轉(zhuǎn)為柵格數(shù)據(jù)。采用粒度推譯法,粒度推譯包括上推和下推。其中最常用的為粒度上推,即把精微尺度上的觀察、試驗(yàn)及模擬結(jié)果外推到較大尺度的過程。它是研究成果的“粗粒化”,即不斷降低景觀格局的空間分辨率[30]。由于研究區(qū)面積較小,所以采用的粒度大小為柵格大小,從100 m到3 000 m,即100,200,400,600,800,1000,2 000,3 000 m為像元大小,得到8幅景觀類型柵格圖(Grid),用Fragstats 3.3計(jì)算景觀格局指數(shù),景觀指數(shù)依據(jù)前人的研究[31]和研究區(qū)特點(diǎn),景觀指數(shù)包括5個(gè)面積/周長/密度指數(shù),7個(gè)形狀指數(shù),8個(gè)蔓延度指數(shù);景觀水平上選取的景觀指數(shù)包括4個(gè)面積/周長/密度指數(shù),4個(gè)形狀指數(shù),3個(gè)蔓延度指數(shù)和4個(gè)多樣性指數(shù)。各指數(shù)計(jì)算公式及生態(tài)學(xué)意義見有關(guān)參考文獻(xiàn)[32]。

    2.3 相關(guān)性分析

    在SPSS 11.5環(huán)境下,主要采用二元定距變量相關(guān)性分析。二元定距變量的相關(guān)分析是指通過計(jì)算定距變量間兩兩相關(guān)的相關(guān)系數(shù),對2個(gè)或2個(gè)以上定距變量之間兩兩相關(guān)的程度進(jìn)行分析。在二元定距變量的相關(guān)分析過程中常用的幾個(gè)相關(guān)系數(shù)是Person 簡單相關(guān)系數(shù)、Spearman 和Kendall′s tua-b等級相關(guān)系數(shù)[33]。本文運(yùn)用的是person簡單相關(guān)系數(shù)來進(jìn)行相關(guān)性分析。

    3 結(jié)果與分析

    3.1 類型水平景觀指數(shù)的空間粒度變化行為測度

    空間粒度“粗?;睂σ恍┚坝^格局指數(shù)有明顯的影響,即同一景觀不同粒度時(shí),用同一種景觀格局指數(shù)表示的結(jié)果不同。類型水平景觀指數(shù)的相關(guān)性越強(qiáng),指數(shù)值越相近[34]。根據(jù)景觀指數(shù)對空間粒度的響應(yīng)關(guān)系和粒度效應(yīng)關(guān)系,選取20個(gè)景觀指數(shù)中的14個(gè)重要景觀指數(shù),分為3類,從地類景觀單個(gè)指標(biāo)特征柱狀圖(圖1)來看,其粒度效應(yīng)表現(xiàn)出斑塊密度(圖1A)、景觀形狀指數(shù)(圖1D)、平均周長面積比(圖1M)、聚集度指數(shù)(圖1K)、相似鄰近比指數(shù)(圖1J),聚合指數(shù)(圖1G)都表現(xiàn)出逐漸減少的趨勢。而從14個(gè)類型水平景觀指數(shù)變化柱狀圖變化特征來看,其粒度效應(yīng)有以下特征:

    (1) 第1類景觀指數(shù)隨著粒度的增加而減小,且具有比較明顯的粒度效應(yīng)關(guān)系,可預(yù)測性較強(qiáng)。這類景觀指數(shù)主要有斑塊密度、景觀形狀指數(shù)、平均周長面積比、聚集度指數(shù)、相似鄰近比指數(shù),聚合指數(shù)(圖1G)。以聚合指數(shù)為例,聚合指數(shù)表示景觀類型內(nèi)部的團(tuán)聚程度[29],如果一個(gè)景觀由許多離散的小斑塊組成,其聚合指數(shù)較?。划?dāng)景觀中以少數(shù)大斑塊為主或同一類型斑塊高度相連時(shí),其聚合指數(shù)較大。從聚合指數(shù)的變化趨勢來看(圖1G),其地類景觀指數(shù)下降較為明顯,其中下降最明顯的是建設(shè)用地和水域,其次是未利用地、林地,可見隨著粒度的增加,類型景觀的形狀愈加簡單,并且在1 000 m粒度時(shí),建設(shè)用地、耕地和水域都出現(xiàn)了拐點(diǎn)。

    (2) 第2類景觀指數(shù)隨著粒度的增加呈現(xiàn)增加的趨勢,同樣具有比較明顯的粒度效應(yīng)關(guān)系,代表性指數(shù)有歸一化形狀指數(shù)(圖1F)。歸一化形狀指數(shù)隨著粒度的增加,其柱狀圖呈逐漸上升的趨勢,可見隨著粒度的增加景觀類型的形狀便越來越簡單,由圖1F可以看出,所有地類的歸一化形狀指數(shù)都大于0.1,表明研究區(qū)所有地類景觀形狀比較簡單,并且隨著粒度的增加邊界趨勢于平滑、規(guī)則。

    (3) 第3類景觀格局指數(shù)隨粒度的變化柱狀圖差異較大,無明顯規(guī)律,代表性指數(shù)如周長面積變異系數(shù)、斑塊密度、景觀形狀變異系數(shù)、最大斑塊指數(shù)、有效粒度面積和平均周長面積比(圖1)。有效粒度面積用于比較斑塊的平均面積大小。景觀總面積不變時(shí),有效粒度面積變大,反映該類型景觀面積增加,其在景觀中的比例加大。通過分析,粒度變化對耕地的有效粒度面積指數(shù)影響最大,尤其在800~2 000 m粒度其值變化較為劇烈,而其他地類的有效粒度面積指數(shù)隨粒度的增加變化不明顯。其他幾類景觀指數(shù)隨粒度同樣表現(xiàn)出不規(guī)律的變化。

    3.2 整體水平景觀指數(shù)的空間粒度變化行為測度

    利用Fragstats 3.3軟件分析的結(jié)果,進(jìn)行整體水平上的景觀特征對比分析?,敿{斯河流域2014年整體景觀特征的粒度效應(yīng)見圖2。

    圖1 類型水平景觀特征指數(shù)的粒度變化曲線

    在整體水平上,選取特征較明顯的10個(gè)景觀指數(shù)進(jìn)行分析,其粒度變化柱狀圖見圖2。由圖2可以看出,2014年瑪納斯河流域的類型景觀指數(shù)大部分隨著粒度的增加呈現(xiàn)減小的趨勢。隨著粒度的增加,蔓延度指數(shù)維變化差值為2~3,景觀形狀變異系數(shù)變化差值維持在5以內(nèi)、凝聚度指數(shù)變化差值維持為2~3,聚合指數(shù)變化差值維持在5以內(nèi)、景觀形狀指數(shù)變化差值維持在2以內(nèi),說明這些景觀格局指數(shù)粒度效應(yīng)受空間格局特征的影響不大,主要表現(xiàn)為逐漸下降勢降,也有少數(shù)情況下指數(shù)的粒度效應(yīng)曲線表現(xiàn)為基本趨于穩(wěn)定,如平均分維數(shù)維持在1.2左右、平均形狀指數(shù)維持在2.5左右、香農(nóng)多樣性指數(shù)維持2.5左右、香農(nóng)均勻度指數(shù)維持在0.8左右、景觀豐度維持在6左右,可見景觀中斑塊形狀逐漸規(guī)則且斑塊之間相互鄰接。

    圖2 整體水平景觀特征指數(shù)的粒度變化曲線

    為了更好地研究各個(gè)時(shí)期景觀指數(shù)隨粒度的變化特征,此次研究將8種粒度在2個(gè)時(shí)期的景觀特征指數(shù)進(jìn)行對比,其結(jié)果見表1。從表1看出,隨著粒度的增加,景觀形狀變異系數(shù)呈現(xiàn)減小趨勢,2014年隨著粒度增加景觀形狀變異系數(shù)從66.28減少至53.01,相較2014年,2000年景觀形狀變異系數(shù)最大,最大可達(dá)到71.626,可見2000年景觀的形狀指數(shù)最為不規(guī)則。平均分維數(shù)、蔓延度指數(shù)、凝聚度和景觀豐度、聚合度指數(shù)幾乎沒有變化。香農(nóng)多樣性指數(shù)由2000年的1.38增加到2014年的1.45,說明2014年拼塊鄰接程度達(dá)到最高。隨著粒度的增加,景觀形狀指數(shù)逐漸降低,2000年最大達(dá)到56.285,2014年最小值達(dá)到9.03,其特征指數(shù)變化與景觀形狀變異系數(shù)變化趨勢相一致。聚合度指數(shù)與香農(nóng)均勻度和香農(nóng)多樣性表現(xiàn)出同樣的變化規(guī)律,都在2014年達(dá)到最大值,可見2014年景觀類型面積增加,并趨于均勻。凝聚度、景觀豐度在2000年和2014年處于穩(wěn)定,幾乎沒有變化波動。蔓延度在2000年達(dá)到最大,2014年達(dá)到最小。總的來看:(1)粒度發(fā)生變化時(shí),多數(shù)指標(biāo)年際間的變化規(guī)律不一樣,使得指標(biāo)的協(xié)同變化也不一樣;(2)聚合度指數(shù)與香農(nóng)均勻度和香農(nóng)多性、景觀形狀指數(shù)與景觀形狀變異系數(shù)、聚合指數(shù)和景觀豐度的變化在粒度間變化趨于一致。

    3.3 景觀形態(tài)的自相似性及其對粒度效應(yīng)的響應(yīng)

    在各種粒度下,瑪納斯河流域2000—2014年各類景觀斑塊的平均分維數(shù)都大于1(表2),這說明各類景觀形態(tài)都具有分形特征,呈現(xiàn)出一定的自相似性和復(fù)雜性。

    表1 各時(shí)期整體水平景觀特征指數(shù)的差異

    表2 各類景觀斑塊的分維數(shù)對粒度變化的響應(yīng)

    從平均水平來看,草地、耕地和未利用地分維數(shù)平均值達(dá)到1.03,說明景觀相對比較不規(guī)則和復(fù)雜,林地、建設(shè)用地、水域景觀平均分維數(shù)達(dá)到1.02,說明景觀相對比較規(guī)則簡單,這也符合研究區(qū)的實(shí)際情況,因?yàn)?(1) 瑪納斯河流域流域地處干旱區(qū),有“沙漠綠洲”之稱,環(huán)境惡劣,該區(qū)的未利用地所占面積比重比較大,而且草地的管理也較少,從而導(dǎo)致草地景觀比較復(fù)雜的格局;(2) 2000年退耕還林(草)政策的實(shí)施,使得耕地大面積減少,林草地大面積增加,并且得到了有序的種植和管理;(3) 瑪納斯河流域內(nèi)由于居住地的轉(zhuǎn)移和工礦用地的建設(shè),使得建設(shè)用地景觀的形態(tài)格局復(fù)雜化程度接近于荒草地。表2還表明,各類景觀斑塊的分維數(shù)對粒度變化的響應(yīng)不同,它們的分維數(shù)隨著粒度的變化基本保持不變;在粒度較小時(shí),各類景觀斑塊之間的分維數(shù)差異較大,而隨著粒度增大,各類景觀斑塊之間的分維數(shù)差異減小。其中,草地和耕地、水域、建設(shè)用地分維數(shù)在400~2 000 m 粒度內(nèi),具有明顯的隨粒度增加先減小后增加的趨勢。2000年和2014年其他類型景觀的分維數(shù)在400~800 m,800~1 000 m 粒度范圍內(nèi)變化較小,說明這些景觀的形態(tài)結(jié)構(gòu)在這一粒度范圍內(nèi)具有自相似性,在1 000~2 000 m,2 000~3 000 m粒度時(shí)各類型景觀分維數(shù)略有波動,但整體呈下降趨勢。

    3.4 類型景觀指數(shù)與粒度的相關(guān)分析

    采用相關(guān)性分析,獲得了2000—2014年瑪納斯河流域各地類景觀指數(shù)與粒度的相關(guān)性(表4)。從表3,4可以看出,聚集度指數(shù)、斑塊密度和相似鄰近比、聚合指數(shù)、平均面積周長比、散步與并列指數(shù)與粒度的相關(guān)性在類型水平上呈現(xiàn)極顯著負(fù)相關(guān),歸一化形狀指數(shù)與粒度呈極顯著正相關(guān),散步與并列指數(shù)呈兩級分化,即呈正相關(guān)又有負(fù)相關(guān)。水域、草地、耕地和林地的平均面積周長比與粒度的相關(guān)性不高,而且平均面積周長比與粒度的相關(guān)性呈負(fù)相關(guān)。對于瑪納斯河流域而言,建設(shè)用地和林草地之間的轉(zhuǎn)化頻繁。所以重點(diǎn)研究三者在斑塊密度、歸一化形狀指數(shù)、聚集度指數(shù)與粒度之間的相關(guān)性。從草地和林地類型景觀指數(shù)與粒度的相關(guān)性來看,斑塊密度與粒度為極顯著負(fù)相關(guān),歸一化形狀指數(shù)與粒度極顯著正相關(guān),聚集度指數(shù)與粒度則呈現(xiàn)極顯著負(fù)相關(guān)。從建設(shè)用地類型景觀指數(shù)與粒度的相關(guān)性來看,斑塊密度與粒度呈極顯著負(fù)相關(guān),而歸一化形狀指數(shù)和聚集度指數(shù)與粒度的相關(guān)性則與草地、林地類似。

    表3 2000年各地類景觀指數(shù)與粒度的相關(guān)性

    注:*p<0.05時(shí),相關(guān)性是顯著的,**p<為0.01時(shí),相關(guān)性是極顯著的。下表同。

    表4 2014年各地類景觀指數(shù)與粒度的相關(guān)性

    3.5 整體景觀指數(shù)與粒度的相關(guān)性

    從表5可以看出,香農(nóng)均勻度指數(shù)、香農(nóng)多樣性指數(shù)與粒度相關(guān)性在2000年、2014年呈兩極分化,在2000年呈負(fù)相關(guān),相關(guān)性系數(shù)都為-0.227,在2014年呈正相關(guān),相關(guān)性系數(shù)為0.066,景觀形狀指數(shù)、平均分維數(shù)、蔓延度指數(shù)與粒度相關(guān)性在2000年和2014年均呈負(fù)相關(guān),說明周邊景觀多樣性降低,景觀格局由復(fù)雜變?yōu)楹唵巍?/p>

    表5 各時(shí)期景觀指數(shù)與粒度的相關(guān)性

    4 結(jié)論與討論

    尺度效應(yīng)在土地利用數(shù)據(jù)上表現(xiàn)明顯,土地利用景觀特性的表現(xiàn)在不同的空間尺度上具有較大的差異。土地利用景觀隨著尺度的變化而發(fā)生改變,呈現(xiàn)出明顯的相關(guān)性和復(fù)雜性。通過對瑪河流域2000年和2014年景觀格局指數(shù)的粒度“粗?;毖芯浚治隽爽敿{斯河流域景觀格局的多粒度尺度效應(yīng)。從景觀的類型水平和整體水平分別進(jìn)行分析,實(shí)現(xiàn)了瑪納斯河流域景觀格局研究從部分到整體的結(jié)合。研究發(fā)現(xiàn)“臨界粒度”及“粗粒化”現(xiàn)象在瑪納斯河流域景觀格局普遍存在?,敽恿饔蚓坝^格局指數(shù)的“臨界粒度”“粗?;爆F(xiàn)象較為明顯。在100~3 000 m 粒度范圍內(nèi),400~800 m,2 000~3 000分辨率下瑪河流域的景觀格局指數(shù)的“臨界粒度”為400 m和2 000 m。景觀形態(tài)自相似性的研究可以很好地反映景觀格局的邊界特征,研究中,瑪納斯河流域景觀形態(tài)具有分形特征,并且呈現(xiàn)出一定的自相似性和復(fù)雜性。各類景觀斑塊的分維數(shù)對粒度變化的響應(yīng)不同,分維數(shù)隨粒度的增大呈現(xiàn)逐漸下降趨勢,表明景觀類型邊界趨于簡單化。在粒度較小時(shí),各類景觀斑塊之間的分維數(shù)差異較大,而隨著粒度增大,各類景觀斑塊之間的分維數(shù)差異減小。從平均水平看,草地、林地和水域景觀相對比較規(guī)則簡單,耕地和建設(shè)用地景觀比較復(fù)雜,而未利用地介于它們之間。本研究表明,瑪納斯河流域的景觀分維數(shù)隨粒度的增加而逐漸減小。通過相關(guān)性分析,一方面定量反映了所選景觀格局指數(shù)受粒度變化影像的相關(guān)性程度,另一方面可為后續(xù)瑪河流域因“退耕還林(草)”工程引起的景觀格局變化研究提供參考。

    本文在對景觀指數(shù)的粒度效應(yīng)研究中,針對從測量尺度向本征尺度的過渡進(jìn)行了說明,但文章只從空間上分析了各地類景觀指數(shù)的粒度效應(yīng),并沒有從時(shí)間變化角度分析景觀指數(shù)的粒度效應(yīng),這一點(diǎn)還需深入研究。因?yàn)榫坝^格局是不同景觀單元和生態(tài)過程在一定時(shí)間和一定空間內(nèi)相互作用的表現(xiàn),它在空間或時(shí)間單方面的特征都不能代表其本質(zhì)的規(guī)律性。判別景觀格局的特征尺度是進(jìn)行格局分析的前提,而這個(gè)特征尺度不僅包括空間尺度,還應(yīng)包含時(shí)間尺度,二者缺一不可。目前在我國開展的大量景觀格局分析中,更多的側(cè)重于景觀格局指數(shù)的計(jì)算與分析,而將景觀格局分析與實(shí)際問題相結(jié)合的研究相對較少。

    [1] 鄔建國.景觀生態(tài)學(xué):概念和理論[J].生態(tài)學(xué)雜志,2000,19(1):42-52.

    [2] Marceau D I, Hay G J. Remote sensing contributions to the scale issue[J]. Canadian Journal of Remote Sensing, 1999,25(4):357-366.

    [3] 劉明亮,唐先明,劉紀(jì)遠(yuǎn),等.基于1 km格網(wǎng)的空間數(shù)據(jù)尺度效應(yīng)研究[J].遙感學(xué)報(bào),2001,5(3):183-190.

    [4] 趙磊.土地利用空間數(shù)據(jù)轉(zhuǎn)換的尺度效應(yīng)研究[D].北京:中國農(nóng)業(yè)大學(xué),2005.

    [5] 肖篤寧.論現(xiàn)代景觀科學(xué)的形成與發(fā)展[J].地理科學(xué),1999,19(4):379-384.

    [6] 肖篤寧,李秀珍.當(dāng)代景觀生態(tài)學(xué)的進(jìn)展和展望[J].地理科學(xué),1997,17(4):356-363.

    [7] 李秀珍.從第十五屆美國景觀生態(tài)學(xué)年會看當(dāng)前景觀生態(tài)學(xué)發(fā)展的熱點(diǎn)和前沿[J].生態(tài)學(xué)報(bào),2000,20(6):1113-1115.

    [8] 鄔建國.景觀生態(tài)學(xué)中的十大研究論題[J].生態(tài)學(xué)報(bào),2004,24(9) ;2074-2076.

    [8] Peternson G D. Scaling ecological dynamics: Self-organization hierarchical structure and ecological resilience[J]. Climate Change, 2000,44:291-309.

    [9] Cai B F, Yu R. Comparison on spatial scale analysis methods in landscape ecology[J]. Acta Ecological Sinica, 2008,25(5):2279-2287.[10] Levin S A. The problem of pattern and scale in ecology[J]. Ecology, 1992,73:1943-1967.

    [11] Wu J G. Landscape ecology: Concept and theory[J]. Chinese Journal of Ecology, 2000,19(1):42-52.

    [12] Zhao W W, Fu B, Chen L D. Some fundamental issues in scaling[J]. Advance in Earth Sciences, 2002,17(6):905-911.

    [13] Lam N, Quattrochi D A. On the issues of scale, resolution, and fractural analysis in the mapping sciences[J]. The Professional Geographer, 1992,(44):88-987.

    [14] Qi Y, Wu J. G. Effects of changing spatial resolution on the results of landscape pattern analysis using spatial auto-correlation indice[J]. Landscape Ecology, 1996,11(1):39-49.

    [15] Wu J, Shen W, Sun W, et al. Empirical patterns of the effects of changing scale on landscape metrics[J]. Landscape Ecology, 2002,(17):761-782.

    [12] Wu J. Effects of change scale on landscape pattern analysis: Scaling relations[J]. Landscape Ecology, 2004,19:125-138.

    [17] 呂一河,傅伯杰.生態(tài)學(xué)中的尺度及尺度轉(zhuǎn)換方法[J].生態(tài)學(xué)報(bào),2001,21(12):2097-2104.

    [18] 蔡博峰,于嶸.景觀生態(tài)學(xué)中的尺度分析方法[J].生態(tài)學(xué)報(bào),2008,28(5):2279-2287.

    [19] 申衛(wèi)軍,鄔建國,林永標(biāo),等.空間粒度變化對景觀格局分析的影響[J].生態(tài)學(xué)報(bào),2003,23(12):2506-2519.

    [20] 趙文武,傅伯杰,陳利頂.景觀指數(shù)的粒度變化效應(yīng)[J].第四紀(jì)研究,2003,23(3):326-333.

    [21] 楊存建,劉紀(jì)元,張?jiān)鱿?等.土地利用數(shù)據(jù)尺度轉(zhuǎn)換的精度損失分析[J].山地學(xué)報(bào),2001,19(3):258-264.

    [22] 楊麗,甄霖,謝高地,等.涇河流域景觀指數(shù)的粒度效應(yīng)分析[J].資源科學(xué),2007,29(2):183-187.

    [23] 孟陳,李俊祥,朱穎,等.粒度變化對上海市景觀格局分析的影響[J].生態(tài)學(xué)雜志,2007,26(3):1138-1142.

    [24] 何丙輝,徐霞,辜世賢.潼南縣土地利用景觀格局特征分析[J].水土保持研究,2005,12(5):130-133.

    [25] 焦峰,溫仲明,王飛,等.黃土丘陵縣域尺度整體景觀格局分析[J].水土保持學(xué)報(bào),2005,19(2):167-170.

    [26] 呂一河,陳利頂,傅伯杰.縣域人類活動與景觀格局分析[J].生態(tài)學(xué)報(bào),2004,24(9):1833-1838.

    [27] 張金屯,邱揚(yáng),鄭鳳英.景觀格局的數(shù)量研究方法[J].山地學(xué)報(bào),2000,18(4):346-352.

    [28] 鄔建國.景觀生態(tài)學(xué):格局、過程、尺度與等級 [M].2版.北京:高等教育出版社,2007:117-119.

    [29] 賴先齊.綠洲鹽漬化棄耕地生態(tài)重建研究[M].北京:中國農(nóng)業(yè)出版社,2007.

    [30] 游麗平,林廣發(fā),楊陳照,等.景觀指數(shù)的空間尺度效應(yīng)分析:以廈門島土地利用格局為例[J].地球信息科學(xué),2008,10(1):74-79.

    [31] 周偉,鐘星,袁春.1∶10000比例尺土地利用景觀指數(shù)的粒度效應(yīng)分析[J].中國土地科學(xué),2010,24(11):20-26.

    [32] 盧玲.景觀格局指數(shù)及各指數(shù)意義[EB/OL]http:∥wenku. baidu. com/view/869 ea62 ccfc789 eb172 dc8 a1. html.2011年6月7號上傳,2015年1月4日下載

    [33] 余建英,何旭宏.數(shù)據(jù)統(tǒng)計(jì)分析與SPSS應(yīng)用[M].北京:人民郵電出版社,2003.

    [34] 趙文武,傅伯杰,陳利頂.景觀指數(shù)的粒度變化效應(yīng)[J].第四紀(jì)研究,2003,23(3):326-333.

    InfluenceofGrainSizeChangeonLandUseLandscapePatterninManasiRiverBasininRecent14Years

    HAN Rong, TANG Xiangling

    (DepartmentofGeography,CollegeofScience,ShiheziUniversity,Shihezi,Xinjiang832003,China)

    Based on landscape heterogeneity and spatial differences in land use, landscape pattern indices of measure in time and space ‘coarse’, using RS and GIS software and landscape ecology methods, based on Manas River Basin of the remote sensing image in the period from 2000 to 2014, landscape type in the Manas River Basin, landscape type map grid, the important landscape index, from the overall level of landscape types and the landscape, land using effect and mechanism of the changes of the landscape pattern due to changing grain size in the Manas River Basin were analyzed, and the correlation of the landscape pattern indices with the particle sizes was analyzed. The results showed that in 100~3 000 m particle size range, Manas River Basin landscape pattern indices critical granularity ‘coarse’ phenomenon is obvious, overall, in the 400~800 m resolution critical size of Manas River Basin landscape pattern index was 400 m, suitable for calculating the size range for 200~800 m. So in the use of remote sensing image prediction, comparison and evaluation of the Manas River Valley landscape, the attention should be paid to the effects of grain size and of certain size conversion. Manas River Basin landscape form has the characteristic of fractal, fractal dimension of various types of landscape patches of size change in response to different fractal dimensions with particle size increase gradually decline trend, indicating that the types of landscape boundary tended to be simple. Quantitative correlation analysis reflected the degree of correlation of landscape index influenced by the change of the particle size, which can provide a reference for the subsequent Manas River Basin due to the project of returning farmland to forest (grass) caused by the study on the change of landscape pattern, in nearly 14 years from 2000 to 2014, landscape diversity in Manas River Basin became from the complexity to simplicity, diversity index decreased gradually.

    land use; landscape pattern; multi-scale; grain size change; correlation analysis; Manasi River Basin

    2016-05-24

    :2016-06-25

    石河子大學(xué)高層次人才項(xiàng)目(RCZX201417)

    韓蓉(1990—),女,新疆庫爾勒人,碩士研究生,研究方向:景觀生態(tài)學(xué)。E-mail:1227677398@qq.com

    唐湘玲(1976—),女,湖南人,副教授,博士, 研究方向:全球變化與生態(tài)響應(yīng)。E-mail:Txl@sina.com

    P901

    :A

    :1005-3409(2017)02-0194-08

    猜你喜歡
    瑪納斯格局粒度
    粉末粒度對純Re坯顯微組織與力學(xué)性能的影響
    基于矩陣的多粒度粗糙集粒度約簡方法
    格局
    聯(lián)手共建 努力打造大調(diào)解工作格局
    民族史詩入夢來——歌劇《瑪納斯》觀感
    歌劇(2017年11期)2018-01-23 03:41:04
    基于粒度矩陣的程度多粒度粗糙集粒度約簡
    近60年瑪納斯河徑流變化特征分析
    瑪納斯河流域水文特征值及演變規(guī)律——以瑪納斯河為例
    談?wù)劇艾敿{斯”和“瑪納斯奇”這兩個(gè)詞
    語言與翻譯(2015年1期)2015-07-18 11:10:08
    小人物的大格局
    小說月刊(2015年12期)2015-04-23 08:51:10
    春色校园在线视频观看| 中国三级夫妇交换| 国产精品欧美亚洲77777| 卡戴珊不雅视频在线播放| 少妇 在线观看| 久久精品国产亚洲av天美| 国产不卡av网站在线观看| 伦理电影大哥的女人| 午夜日本视频在线| 免费久久久久久久精品成人欧美视频 | 欧美日韩视频精品一区| 日韩中字成人| 丰满迷人的少妇在线观看| 日韩欧美精品免费久久| 97在线人人人人妻| 免费av不卡在线播放| 午夜av观看不卡| 国产片内射在线| 老司机亚洲免费影院| 亚洲精品一区蜜桃| 性高湖久久久久久久久免费观看| 亚洲国产av影院在线观看| 亚洲欧美日韩卡通动漫| 在线亚洲精品国产二区图片欧美| 黄色配什么色好看| 一区在线观看完整版| 大片电影免费在线观看免费| 亚洲欧美中文字幕日韩二区| 国产又色又爽无遮挡免| av播播在线观看一区| 狂野欧美激情性xxxx在线观看| 亚洲国产精品成人久久小说| 午夜影院在线不卡| 日本黄大片高清| 十分钟在线观看高清视频www| 美女福利国产在线| 国产日韩欧美亚洲二区| 制服丝袜香蕉在线| √禁漫天堂资源中文www| 丰满少妇做爰视频| 一个人免费看片子| 国产永久视频网站| 亚洲情色 制服丝袜| 99热这里只有是精品在线观看| 国产成人午夜福利电影在线观看| 视频在线观看一区二区三区| 青青草视频在线视频观看| 国产精品一二三区在线看| 午夜免费男女啪啪视频观看| 国产一区二区三区av在线| 精品第一国产精品| 天天躁夜夜躁狠狠躁躁| 国产有黄有色有爽视频| 国产激情久久老熟女| 精品一区二区三区视频在线| 国产不卡av网站在线观看| freevideosex欧美| 日日爽夜夜爽网站| 国产亚洲精品第一综合不卡 | 热99国产精品久久久久久7| 乱码一卡2卡4卡精品| 欧美日韩综合久久久久久| 亚洲综合色惰| 热99国产精品久久久久久7| 日本欧美视频一区| 视频在线观看一区二区三区| av在线播放精品| 欧美日韩亚洲高清精品| 插逼视频在线观看| 日韩av免费高清视频| 久久久久精品性色| 午夜福利,免费看| 只有这里有精品99| 午夜免费鲁丝| 久久久久久人人人人人| 一级毛片黄色毛片免费观看视频| 爱豆传媒免费全集在线观看| 亚洲人与动物交配视频| 精品少妇黑人巨大在线播放| 色婷婷av一区二区三区视频| 中文乱码字字幕精品一区二区三区| 成人综合一区亚洲| 天天躁夜夜躁狠狠久久av| 黑人猛操日本美女一级片| 国产乱来视频区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一区在线观看完整版| 亚洲欧美日韩卡通动漫| 免费少妇av软件| 国产免费一级a男人的天堂| 亚洲成国产人片在线观看| 国产精品久久久久成人av| 成人综合一区亚洲| 中文字幕另类日韩欧美亚洲嫩草| 如日韩欧美国产精品一区二区三区| av片东京热男人的天堂| av又黄又爽大尺度在线免费看| 男女边吃奶边做爰视频| 免费av中文字幕在线| 欧美另类一区| 国产综合精华液| 久久久久网色| 黑人巨大精品欧美一区二区蜜桃 | 一区二区三区四区激情视频| 热re99久久精品国产66热6| 国产一区有黄有色的免费视频| 看十八女毛片水多多多| 美女主播在线视频| 在线观看美女被高潮喷水网站| 99热6这里只有精品| 精品一品国产午夜福利视频| 日本黄大片高清| av福利片在线| 在线观看免费高清a一片| 亚洲av综合色区一区| 巨乳人妻的诱惑在线观看| 波野结衣二区三区在线| 久久精品久久久久久噜噜老黄| 国产日韩一区二区三区精品不卡| 国产老妇伦熟女老妇高清| av线在线观看网站| 中国美白少妇内射xxxbb| 18+在线观看网站| 免费看光身美女| 老司机影院毛片| 亚洲av日韩在线播放| 这个男人来自地球电影免费观看 | 免费在线观看完整版高清| 亚洲精品一二三| 香蕉丝袜av| 777米奇影视久久| 国产老妇伦熟女老妇高清| 两个人免费观看高清视频| 国产精品成人在线| 人人澡人人妻人| 成人二区视频| 成人免费观看视频高清| 成人午夜精彩视频在线观看| 这个男人来自地球电影免费观看 | av又黄又爽大尺度在线免费看| 国产黄色视频一区二区在线观看| 免费看不卡的av| 哪个播放器可以免费观看大片| av卡一久久| 亚洲欧美色中文字幕在线| 亚洲一码二码三码区别大吗| 黄色毛片三级朝国网站| 免费人成在线观看视频色| 街头女战士在线观看网站| 男男h啪啪无遮挡| 一级爰片在线观看| 欧美日韩精品成人综合77777| 欧美日韩国产mv在线观看视频| 日韩视频在线欧美| 成人国产麻豆网| 热99国产精品久久久久久7| 蜜桃在线观看..| 久久精品国产自在天天线| 亚洲国产最新在线播放| 高清视频免费观看一区二区| 亚洲精品色激情综合| 一二三四在线观看免费中文在 | 国产精品三级大全| 免费黄网站久久成人精品| 你懂的网址亚洲精品在线观看| 有码 亚洲区| 妹子高潮喷水视频| 亚洲av日韩在线播放| 久久午夜综合久久蜜桃| 国产乱人偷精品视频| 国产女主播在线喷水免费视频网站| 男女边吃奶边做爰视频| 九九爱精品视频在线观看| 国产免费福利视频在线观看| 啦啦啦视频在线资源免费观看| 在线观看人妻少妇| 五月伊人婷婷丁香| 热99国产精品久久久久久7| 国产毛片在线视频| 大码成人一级视频| 三上悠亚av全集在线观看| 免费女性裸体啪啪无遮挡网站| 狂野欧美激情性xxxx在线观看| 韩国精品一区二区三区 | 免费看av在线观看网站| 最近中文字幕2019免费版| 黑人猛操日本美女一级片| 狠狠精品人妻久久久久久综合| 国产免费福利视频在线观看| 午夜福利视频精品| 午夜免费观看性视频| 成人黄色视频免费在线看| 又黄又粗又硬又大视频| 久久久久精品人妻al黑| 国产激情久久老熟女| 欧美日韩亚洲高清精品| 亚洲精品乱码久久久久久按摩| 国产av国产精品国产| 天天影视国产精品| 国产欧美日韩综合在线一区二区| 97人妻天天添夜夜摸| 亚洲精品美女久久av网站| av又黄又爽大尺度在线免费看| 久久99热6这里只有精品| 国产免费现黄频在线看| 成人国语在线视频| 久久久精品区二区三区| 成人毛片60女人毛片免费| 黄色怎么调成土黄色| 一个人免费看片子| 视频区图区小说| 国产精品久久久久久精品电影小说| 在线精品无人区一区二区三| 中文精品一卡2卡3卡4更新| 亚洲精品456在线播放app| 一级毛片黄色毛片免费观看视频| 久久国产亚洲av麻豆专区| 亚洲丝袜综合中文字幕| 国产成人免费观看mmmm| a级毛色黄片| 中文字幕另类日韩欧美亚洲嫩草| 美女福利国产在线| 人妻一区二区av| 精品一区二区免费观看| 国产精品国产三级国产专区5o| 免费高清在线观看日韩| 国产片内射在线| 国产无遮挡羞羞视频在线观看| 欧美老熟妇乱子伦牲交| 国产高清不卡午夜福利| 欧美日韩精品成人综合77777| 9色porny在线观看| 在线观看www视频免费| 国产xxxxx性猛交| 免费观看av网站的网址| 视频区图区小说| 久久人人爽av亚洲精品天堂| 国产精品一区www在线观看| 午夜影院在线不卡| 深夜精品福利| 两性夫妻黄色片 | 热re99久久精品国产66热6| 少妇的逼水好多| 久久人妻熟女aⅴ| 精品久久蜜臀av无| 狠狠精品人妻久久久久久综合| 男女边摸边吃奶| 国产免费一区二区三区四区乱码| 精品人妻熟女毛片av久久网站| 又粗又硬又长又爽又黄的视频| 中文字幕免费在线视频6| 男女啪啪激烈高潮av片| 国产精品熟女久久久久浪| 亚洲中文av在线| 久久久精品94久久精品| 国产色婷婷99| 国产成人免费观看mmmm| 丝瓜视频免费看黄片| 视频在线观看一区二区三区| 我要看黄色一级片免费的| 十八禁网站网址无遮挡| 国产高清国产精品国产三级| 成年动漫av网址| 国产麻豆69| 免费看av在线观看网站| 只有这里有精品99| 少妇人妻精品综合一区二区| 午夜福利网站1000一区二区三区| 丰满少妇做爰视频| 夫妻性生交免费视频一级片| 久久韩国三级中文字幕| 黄色一级大片看看| 色视频在线一区二区三区| 日韩一区二区视频免费看| www.色视频.com| 国产成人欧美| 久久免费观看电影| 美女国产高潮福利片在线看| 免费观看av网站的网址| 咕卡用的链子| 女人久久www免费人成看片| 久久午夜综合久久蜜桃| 精品一区二区三卡| 中文欧美无线码| 久久热在线av| 一边摸一边做爽爽视频免费| av国产久精品久网站免费入址| 麻豆精品久久久久久蜜桃| 久久av网站| 91国产中文字幕| 在线观看国产h片| 欧美日韩av久久| 国产日韩欧美在线精品| av.在线天堂| 看非洲黑人一级黄片| 校园人妻丝袜中文字幕| 少妇熟女欧美另类| 久久av网站| 边亲边吃奶的免费视频| 777米奇影视久久| 午夜激情av网站| 在线亚洲精品国产二区图片欧美| 熟女av电影| 亚洲欧美成人综合另类久久久| 在线 av 中文字幕| 一级毛片黄色毛片免费观看视频| 国产免费一区二区三区四区乱码| 成人免费观看视频高清| 亚洲成人手机| 五月开心婷婷网| 久久综合国产亚洲精品| 人妻一区二区av| √禁漫天堂资源中文www| 免费黄频网站在线观看国产| 我的女老师完整版在线观看| 丝瓜视频免费看黄片| 国产一区二区三区av在线| 少妇熟女欧美另类| 这个男人来自地球电影免费观看 | 9191精品国产免费久久| 国产一区二区激情短视频 | 欧美成人精品欧美一级黄| 久久久亚洲精品成人影院| 欧美 亚洲 国产 日韩一| 90打野战视频偷拍视频| 免费观看无遮挡的男女| 看十八女毛片水多多多| 少妇猛男粗大的猛烈进出视频| 91精品三级在线观看| √禁漫天堂资源中文www| 青春草视频在线免费观看| 亚洲国产精品专区欧美| 黑人巨大精品欧美一区二区蜜桃 | 亚洲欧美色中文字幕在线| 亚洲欧美成人综合另类久久久| 久久久久久人人人人人| 最黄视频免费看| 国产精品久久久久久久久免| av免费观看日本| 久久99热6这里只有精品| 日韩一本色道免费dvd| 9热在线视频观看99| 亚洲欧美精品自产自拍| 中文字幕精品免费在线观看视频 | 国产黄色视频一区二区在线观看| 男的添女的下面高潮视频| 国产片特级美女逼逼视频| 免费观看在线日韩| 春色校园在线视频观看| 亚洲国产精品一区二区三区在线| av国产久精品久网站免费入址| 精品酒店卫生间| 免费观看性生交大片5| kizo精华| 9色porny在线观看| 99re6热这里在线精品视频| 在线观看美女被高潮喷水网站| 国产又色又爽无遮挡免| 久久影院123| 国产精品人妻久久久久久| 日韩av免费高清视频| 国产免费又黄又爽又色| 2018国产大陆天天弄谢| 日韩电影二区| 日韩av免费高清视频| 晚上一个人看的免费电影| 久久久久精品人妻al黑| 777米奇影视久久| 久久久国产欧美日韩av| 亚洲av.av天堂| 欧美成人精品欧美一级黄| 亚洲色图 男人天堂 中文字幕 | 国产高清不卡午夜福利| 精品国产国语对白av| 欧美日韩综合久久久久久| 在线观看免费高清a一片| 久久精品国产亚洲av涩爱| www.熟女人妻精品国产 | 在线观看国产h片| 国产亚洲精品第一综合不卡 | 亚洲欧美中文字幕日韩二区| 不卡视频在线观看欧美| 日韩大片免费观看网站| 女性被躁到高潮视频| 国产不卡av网站在线观看| 亚洲av国产av综合av卡| 高清av免费在线| 国产精品蜜桃在线观看| 久久精品人人爽人人爽视色| 我的女老师完整版在线观看| 一区二区三区乱码不卡18| 久热这里只有精品99| 国产免费一级a男人的天堂| 综合色丁香网| 熟妇人妻不卡中文字幕| 日本-黄色视频高清免费观看| 精品久久国产蜜桃| 在现免费观看毛片| 一区二区日韩欧美中文字幕 | 精品国产乱码久久久久久小说| 韩国av在线不卡| 欧美精品一区二区大全| 狂野欧美激情性bbbbbb| 久久99一区二区三区| 99久久综合免费| 久久久久久久久久久免费av| 午夜福利视频精品| 黄色配什么色好看| 人人妻人人添人人爽欧美一区卜| 少妇被粗大的猛进出69影院 | 中文乱码字字幕精品一区二区三区| 插逼视频在线观看| 最后的刺客免费高清国语| 国产男女超爽视频在线观看| 高清不卡的av网站| 宅男免费午夜| 中文字幕最新亚洲高清| 赤兔流量卡办理| 一级片'在线观看视频| 国产精品 国内视频| 久久 成人 亚洲| 伦理电影大哥的女人| 午夜免费鲁丝| 亚洲性久久影院| 久久影院123| 亚洲,欧美精品.| 在线观看三级黄色| 夜夜爽夜夜爽视频| 2018国产大陆天天弄谢| 午夜影院在线不卡| 亚洲av.av天堂| 亚洲少妇的诱惑av| 精品一区二区免费观看| 九九爱精品视频在线观看| 99热国产这里只有精品6| 麻豆精品久久久久久蜜桃| 亚洲精品aⅴ在线观看| 五月天丁香电影| 日韩av在线免费看完整版不卡| www日本在线高清视频| 国产永久视频网站| 人妻系列 视频| 欧美人与性动交α欧美精品济南到 | 久久青草综合色| 黄色配什么色好看| 国产在视频线精品| 免费观看在线日韩| 麻豆乱淫一区二区| 搡女人真爽免费视频火全软件| 少妇精品久久久久久久| 久久ye,这里只有精品| 久久鲁丝午夜福利片| 王馨瑶露胸无遮挡在线观看| 在线看a的网站| 人妻少妇偷人精品九色| 欧美日韩成人在线一区二区| 晚上一个人看的免费电影| 观看av在线不卡| 最近的中文字幕免费完整| 伊人亚洲综合成人网| 一边摸一边做爽爽视频免费| 日产精品乱码卡一卡2卡三| av天堂久久9| 午夜福利乱码中文字幕| 国产亚洲最大av| 亚洲国产毛片av蜜桃av| 国产精品三级大全| 亚洲性久久影院| xxxhd国产人妻xxx| 免费av不卡在线播放| 国产精品久久久久久av不卡| 国产免费现黄频在线看| 久久国内精品自在自线图片| 伊人亚洲综合成人网| 九色亚洲精品在线播放| 超色免费av| 天堂俺去俺来也www色官网| 最新的欧美精品一区二区| 建设人人有责人人尽责人人享有的| 丝袜人妻中文字幕| 韩国高清视频一区二区三区| 观看av在线不卡| 亚洲国产欧美日韩在线播放| 少妇被粗大猛烈的视频| 青春草国产在线视频| 99热国产这里只有精品6| 久久久久久久国产电影| 中文字幕人妻熟女乱码| 国产69精品久久久久777片| 成人毛片60女人毛片免费| 国产1区2区3区精品| 另类亚洲欧美激情| 日韩精品有码人妻一区| 国产成人aa在线观看| 久久人妻熟女aⅴ| 国产成人a∨麻豆精品| 免费看不卡的av| 国产精品国产三级国产av玫瑰| 国产日韩欧美视频二区| 伦理电影大哥的女人| 91久久精品国产一区二区三区| 成人二区视频| 全区人妻精品视频| 国产成人a∨麻豆精品| 精品人妻在线不人妻| 亚洲伊人久久精品综合| 久久久精品区二区三区| 边亲边吃奶的免费视频| 热re99久久国产66热| 中国美白少妇内射xxxbb| 国产成人精品婷婷| 成年av动漫网址| 国产精品国产三级专区第一集| 午夜福利影视在线免费观看| 午夜福利,免费看| 黑人猛操日本美女一级片| 天天影视国产精品| 五月玫瑰六月丁香| 欧美+日韩+精品| 久久97久久精品| 免费高清在线观看日韩| 男人爽女人下面视频在线观看| 男女下面插进去视频免费观看 | 久久这里只有精品19| 国产成人精品一,二区| 两个人看的免费小视频| 国产熟女欧美一区二区| 日韩一区二区三区影片| 99久久中文字幕三级久久日本| 亚洲三级黄色毛片| 欧美精品国产亚洲| 久久精品国产亚洲av涩爱| 99久国产av精品国产电影| 欧美精品一区二区大全| 亚洲国产成人一精品久久久| 九色成人免费人妻av| 不卡视频在线观看欧美| 久久久久网色| 99香蕉大伊视频| av线在线观看网站| 日本vs欧美在线观看视频| 男人舔女人的私密视频| 日本色播在线视频| 建设人人有责人人尽责人人享有的| 天天操日日干夜夜撸| 亚洲国产色片| 22中文网久久字幕| 九色亚洲精品在线播放| 成人毛片a级毛片在线播放| 精品一区二区三区视频在线| 日韩一本色道免费dvd| 亚洲欧美中文字幕日韩二区| 午夜福利网站1000一区二区三区| 国产老妇伦熟女老妇高清| 交换朋友夫妻互换小说| 人妻少妇偷人精品九色| 亚洲精品视频女| 免费大片18禁| 久久亚洲国产成人精品v| 制服人妻中文乱码| 宅男免费午夜| 久久午夜综合久久蜜桃| 国产午夜精品一二区理论片| 亚洲精品国产色婷婷电影| 亚洲性久久影院| 国产精品女同一区二区软件| 亚洲,欧美精品.| 久久人人爽av亚洲精品天堂| 亚洲成av片中文字幕在线观看 | 一级黄片播放器| 日韩精品有码人妻一区| 考比视频在线观看| 美国免费a级毛片| 最新的欧美精品一区二区| 国产高清三级在线| 在线观看国产h片| 亚洲丝袜综合中文字幕| 丰满乱子伦码专区| 久久久久国产精品人妻一区二区| 久久久久久久亚洲中文字幕| 少妇人妻 视频| 中文字幕免费在线视频6| 国产一区二区在线观看av| 国产深夜福利视频在线观看| 欧美人与性动交α欧美精品济南到 | 欧美最新免费一区二区三区| 最近手机中文字幕大全| 国产福利在线免费观看视频| 99热6这里只有精品| 成人黄色视频免费在线看| 日韩制服丝袜自拍偷拍| 各种免费的搞黄视频| 久久狼人影院| 国产精品女同一区二区软件| 在线天堂中文资源库| 久久午夜福利片| 大话2 男鬼变身卡| 亚洲欧美一区二区三区国产| 中文欧美无线码| 丝瓜视频免费看黄片| 2018国产大陆天天弄谢| 蜜桃在线观看..| 黄色毛片三级朝国网站| 日本黄大片高清| 麻豆乱淫一区二区| 国产精品秋霞免费鲁丝片| 亚洲国产毛片av蜜桃av| 国产精品欧美亚洲77777| 久久久久精品人妻al黑| 丝袜脚勾引网站| 精品一区二区免费观看| 精品少妇黑人巨大在线播放| 狂野欧美激情性bbbbbb| 99久久中文字幕三级久久日本| 大片电影免费在线观看免费|