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

    烏魯木齊市土地利用景觀格局變化及生態(tài)安全分析

    2020-04-09 01:22:28朱亞楠蒲春玲
    生態(tài)科學(xué) 2020年2期
    關(guān)鍵詞:烏魯木齊市格局土地利用

    朱亞楠, 蒲春玲

    烏魯木齊市土地利用景觀格局變化及生態(tài)安全分析

    朱亞楠1, 蒲春玲2,*

    1. 新疆農(nóng)業(yè)大學(xué)管理學(xué)院, 烏魯木齊 830052 2. 新疆農(nóng)業(yè)大學(xué)經(jīng)濟(jì)社會(huì)發(fā)展研究中心, 烏魯木齊 830052

    通過對(duì)烏魯木齊市1990—2017年間的土地利用景觀格局時(shí)空變化特征以及生態(tài)安全進(jìn)行分析, 以烏魯木齊市1990、2000、2010和2017年4期的遙感影像作為數(shù)據(jù)源, 借助ERDAS和ArcGIS軟件對(duì)4期遙感圖像進(jìn)行解譯分析, 并通過FRAGSTAT軟件分析得出烏魯木齊市景觀格局指數(shù), 結(jié)合景觀生態(tài)學(xué)的理論, 分析得出烏魯木齊市土地利用、景觀格局及生態(tài)安全變化情況。結(jié)果表明: (1)烏魯木齊市土地利用/覆被發(fā)生較為明顯的變化, 草地景觀占主導(dǎo)地位, 面積占土地利用總面積的56.22%, 但是草地景觀的優(yōu)勢(shì)度逐漸降低, 未利用地面積僅次于草地, 占土地利用總面積的28.84%。與此同時(shí), 景觀之間以耕地、林地和草地景觀轉(zhuǎn)出和建設(shè)用地和未利用地景觀的轉(zhuǎn)入為主; (2)烏魯木齊市整體土地利用變化速度逐漸加快, 由 3.13% 上升至 6.57%, 同時(shí)土地利用程度由1990年的226.53 上升至 2017 年的 235.08, 說明人類活動(dòng)對(duì)土地生態(tài)系統(tǒng)影響逐漸加深; (3)烏魯木齊市景觀格局變化明顯, 受人類干擾程度加強(qiáng), 景觀破碎化程度逐漸加深, 景觀斑塊更加不規(guī)則和復(fù)雜化, 景觀多樣性和結(jié)構(gòu)的不穩(wěn)定性逐漸增加; (4)烏魯木齊市生態(tài)安全度由1990年的0.56下降至2017年的0.48, 說明城市生態(tài)脆弱性逐漸加深。該研究結(jié)果對(duì)烏魯木齊市未來的土地合理利用, 景觀格局優(yōu)化以及城市生態(tài)建設(shè)和可持續(xù)發(fā)展提供一定的參考。

    土地利用; 景觀格局; 生態(tài)安全; 烏魯木齊

    0 前言

    土地是人類賴以生存和發(fā)展最基本的物質(zhì)保障, 為人類從事一切社會(huì)經(jīng)濟(jì)活動(dòng)提供基礎(chǔ)場(chǎng)所[1-2]。隨著城市化進(jìn)程的快速發(fā)展, 土地利用/土地覆蓋變化(Land Use/Land Cover Change, LUCC)問題引起了國際社會(huì)的普遍關(guān)注。與此同時(shí), 土地利用/土地覆蓋變化研究不僅對(duì)深入了解全球環(huán)境變化起到了重要的作用, 也影響著土地利用景觀格局的變化。將土地利用變化與景觀格局研究相結(jié)合, 可以深入了解土地利用景觀格局變化的發(fā)展現(xiàn)狀及趨勢(shì), 對(duì)促進(jìn)人類更好地理解社會(huì)與自然環(huán)境的相互關(guān)系有重要意義[3-5]。近年來, 運(yùn)用ArcGIS、Erdas和Fragstats軟件相結(jié)合的對(duì)土地利用景觀格局變化的研究已在國內(nèi)外大量展開, 并且研究的尺度由大范圍綜合研究逐步向小范圍深入研究轉(zhuǎn)變[6]。同時(shí), 對(duì)區(qū)域土地景觀生態(tài)安全進(jìn)行分析, 有利于揭示區(qū)域土地生態(tài)環(huán)境的變化趨勢(shì)及生態(tài)安全等級(jí), 對(duì)未來區(qū)域土地合理利用和開發(fā), 以及生態(tài)環(huán)境的保護(hù)具有較為深遠(yuǎn)的意義[7]。

    黨的十九大報(bào)告指出, 要加快生態(tài)文明體制改革, 要推進(jìn)綠色發(fā)展, 著力解決突出環(huán)境問題。烏魯木齊市是坐落在綠洲上的城市, 生態(tài)環(huán)境十分脆弱。烏魯木齊成為國家第三批入選“城市雙修”的試點(diǎn)城市, 通過對(duì)烏魯木齊市土地利用景觀格局的動(dòng)態(tài)變化及生態(tài)安全程度進(jìn)行分析, 揭示近30年烏魯木齊市土地利用景觀格局變化的規(guī)律、生態(tài)安全變化趨勢(shì), 對(duì)未來解決城市發(fā)展過程中所帶來的諸多景觀生態(tài)問題具有重要意義。

    1 研究區(qū)域及數(shù)據(jù)來源

    1.1 研究區(qū)概況

    烏魯木齊, 地處東經(jīng)86°46′10′′—88°59′48′′, 北緯42°54′16′′—44°58′16′′, 土地總面積為14192.78 hm2[16], 是新疆維吾爾自治區(qū)的首府, 是新疆的政治、經(jīng)濟(jì)、文化、科教和交通中心, 是亞歐大陸的“橋頭堡”。烏魯木齊市位于中國西北、新疆中部、亞歐大陸腹地, 地處北天山北麓、準(zhǔn)噶爾盆地南緣, 東西南三面環(huán)山, 北部為廣闊沖積平原, 地勢(shì)由南向西北方向逐漸降低。烏魯木齊市由于深居內(nèi)陸, 遠(yuǎn)離海洋, 屬于典型的溫帶半干旱大陸性氣候, 常年干燥少雨, 晝夜溫差較大。

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

    通過烏魯木齊市1990年、2000年、2010年和2017年云量低于0.5的遙感影像圖(見表1), 利用Erdas 9.2軟件對(duì)四期遙感影像圖進(jìn)行波段組合、圖像拼接、幾何校正等處理, 根據(jù)研究區(qū)邊界范圍對(duì)已處理完成的遙感影像圖進(jìn)行圖像裁剪, 并從中提取研究區(qū)所需的遙感影像數(shù)據(jù)。

    依據(jù)國家于2007年公布的《土地利用現(xiàn)狀分類》分類標(biāo)準(zhǔn), 結(jié)合烏魯木齊市的實(shí)際情況, 利用Erdas 9.2中的監(jiān)督分類與目視解譯相結(jié)合的方法將烏魯木齊市土地利用情況進(jìn)行分類。通過分類將烏魯木齊市土地利用景觀類型劃分為7類, 分別是耕地、園地、林地、草地、水域、建設(shè)用地和未利用地, 分類結(jié)果如圖1和圖2所示。為保證分類結(jié)果的準(zhǔn)確性, 通過對(duì)4期遙感影像分類結(jié)果進(jìn)行精度檢驗(yàn), 其Kappa系數(shù)均達(dá)到0.83以上, 說明解譯結(jié)果的可靠性, 可以滿足本次研究的需求。

    同時(shí), 利用ArcGIS 10.2軟件的空間疊加分析和統(tǒng)計(jì)分析功能對(duì)4期烏魯木齊市土地利用圖進(jìn)行空間疊加分析, 獲取研究期內(nèi)烏魯木齊市土地利用轉(zhuǎn)移矩陣和土地利用類型空間變化圖, 基于此對(duì)烏魯木齊市不同時(shí)期的不同地類的時(shí)空變化進(jìn)行分析。將解譯所得的遙感影像圖的柵格格式的烏魯木齊市4期土地利用景觀類型圖, 借助Fragstats 4.2軟件計(jì)算和分析烏魯木齊市4個(gè)不同時(shí)期景觀指數(shù), 從景觀層次分析烏魯木齊市土地利用景觀格局變化情況。

    表1 烏魯木齊市影像資料特征

    2 研究方法

    2.1 土地利用變化分析

    2.1.1 土地利用時(shí)空變化分析

    通過在ArcGIS 10.2中完成4期土地利用類型圖的疊加分析, 結(jié)合統(tǒng)計(jì)分析功能, 獲取1990—2017年烏魯木齊市土地利用景觀變化圖以及土地景觀利用轉(zhuǎn)移矩陣, 以此對(duì)烏魯木齊市不同時(shí)期的土地利用景觀類型時(shí)空變化進(jìn)行分析。

    2.1.2 土地利用變化速度分析

    通過計(jì)算烏魯木齊市各類土地景觀利用類型的單一土地利用動(dòng)態(tài)度、開發(fā)度和耗減度, 從而可以進(jìn)一步對(duì)烏魯木齊市土地景觀利用類型的動(dòng)態(tài)變化速度進(jìn)行分析。

    單一土地利用動(dòng)態(tài)度計(jì)算公式如下[8]:

    式中,U為研究初期某種土地利用類型的面積,U為研究末期某種土地利用類型的面積,為研究時(shí)段。

    綜合土地利用動(dòng)態(tài)度計(jì)算公式如下[9]:

    式中,LU表示為第類土地利用類型轉(zhuǎn)為非類土地利用類型面積的絕對(duì)值,LU表示為研究初期某一土地利用類型的總面積,為研究時(shí)段。

    土地利用開發(fā)度指在某一時(shí)段內(nèi)某一土地利用類型新開發(fā)的程度, 即揭示各土地利用類型增加的速率, 其中該指數(shù)不適用于未利用地, 計(jì)算公式如下[10]:

    式中, Dab表示由a到b時(shí)段某種土地利用類型新開發(fā)土地的面積, 即由其他土地利用類型轉(zhuǎn)移到該土地利用類型的土地總面積, Ua為a時(shí)刻某土地利用類型的總面積。

    Figure 1 Landscape changes of land use in Urumqi from 1990 to 2000

    圖2 烏魯木齊市2010—2017年土地利用景觀變化圖

    Figure 2 Landscape changes of land use in Urumqi from 2010 to 2017

    土地利用耗減度是指在某一時(shí)段內(nèi)某種土地利用類型實(shí)際消耗的程度, 即揭示各土地利用類型減少的速率, 計(jì)算公式如下[10]:

    式中,C表示由到時(shí)段某種土地利用類型轉(zhuǎn)移到其他土地利用類型的總面積,U為時(shí)刻某土地利用類型的土地總面積。

    2.1.3 土地利用程度分析

    根據(jù)劉紀(jì)遠(yuǎn)[11]所提出的綜合分析法, 結(jié)合烏魯木齊市實(shí)際情況, 分別將烏魯木齊市土地利用程度劃分為4個(gè)級(jí)別, 第1等級(jí)為園地和水域, 第2等級(jí)為耕地、林地和草地, 第3等級(jí)為未利用地, 第4等級(jí)為建設(shè)用地, 計(jì)算公式如下:

    式中,表示為土地利用程度綜合指數(shù),A表示為第級(jí)土地利用程度分級(jí)指數(shù),C表示為第級(jí)土地利用類型面積的比例。

    土地利用程度變化模型計(jì)算公式如下:

    式中,L表示為土地利用程度變化量,LL表示為研究初期和研究末期的土地利用綜合程度指數(shù),表示為土地利用程度的變化率。如果, 式中的L>0或者>0, 則表示該區(qū)域土地利用正處在發(fā)展期, 如果L<0或者<0, 則表示該區(qū)域土地利用正處于衰退期或者調(diào)整期。

    2.2 景觀格局指數(shù)選取

    景觀格局指數(shù)是能夠反映景觀結(jié)構(gòu)組成和空間配置某些方面特征的簡(jiǎn)單定量指標(biāo), 能夠高度濃縮的景觀信息[12-15]。景觀格局的分析主要是由斑塊、斑塊類型和景觀水平3個(gè)水平指數(shù)構(gòu)成。由于斑塊水平指數(shù)只是計(jì)算其他景觀的指數(shù)的基礎(chǔ), 對(duì)景觀格局變化分析無意義, 所以本文以景觀水平和斑塊類型為基礎(chǔ), 選取相關(guān)景觀格局指數(shù), 分析研究區(qū)土地利用景觀格局變化特征和趨勢(shì)。

    2.3 景觀生態(tài)安全分析指數(shù)構(gòu)建

    本文在對(duì)烏魯木齊市景觀格局變化分析的基礎(chǔ)上, 通過了解和掌握景觀格局指數(shù)所具有的生態(tài)學(xué)涵義, 擬采用景觀破碎度、景觀分離度、景觀優(yōu)勢(shì)度等基本判別指數(shù)的基礎(chǔ)上, 構(gòu)建景觀干擾度和景觀脆弱度指數(shù), 與景觀生態(tài)安全度之間建立定量化表大, 以期研究烏魯木齊市土地生態(tài)安全的時(shí)空分異特征。

    2.3.1 景觀干擾度指數(shù)

    在自然以及人為活動(dòng)的干擾下, 使得景觀由最初的單一、均質(zhì)和連續(xù)逐漸趨向復(fù)雜、異質(zhì)以及不連續(xù)的斑塊鑲嵌體。本文通過將景觀破碎度、分離度以及優(yōu)勢(shì)度指數(shù)進(jìn)行疊加構(gòu)建景觀干擾度指數(shù), 計(jì)算公式如下[17]:

    式中,U為景觀干擾度,C為景觀破碎度,D為景觀優(yōu)勢(shì)度,分別為景觀破碎度、分離度和優(yōu)勢(shì)度的權(quán)重,L為景觀類型相對(duì)密度,P為景觀類型相對(duì)蓋度,分別為景觀類型相對(duì)密度和相對(duì)蓋度的權(quán)重,S為景觀類型距離指數(shù),N為景觀斑塊數(shù)量,為景觀斑塊總數(shù),A為景觀斑塊面積,為景觀總面積。

    權(quán)重反映了各個(gè)景觀指數(shù)對(duì)景觀生態(tài)環(huán)境不同的影響程度, 根據(jù)前人的研究成果同時(shí)與研究區(qū)實(shí)際情況相結(jié)合, 綜合分析不同的景觀指數(shù)對(duì)生態(tài)環(huán)境的貢獻(xiàn)程度, 分析得出在景觀干擾度的計(jì)算中認(rèn)為景觀的優(yōu)勢(shì)度最為重要, 其次是分離度和破碎度, 其權(quán)重分別賦值為0.5、0.3和0.2, 在景觀優(yōu)勢(shì)度的計(jì)算中認(rèn)為相對(duì)蓋度的權(quán)重為0.6, 相對(duì)密度的權(quán)重為0.4[18]。

    2.3.2 景觀脆弱度指數(shù)

    景觀脆弱度是指不同景觀類型所代表的生態(tài)系統(tǒng)當(dāng)受到外界干擾后的易損性, 景觀脆弱度值的大小是由景觀在自然演替中所處的階段決定的[18–20]。土地利用程度變化不僅受到土地本身自然屬性的影響, 同時(shí)反映了人為因素和自然因素的綜合干擾的結(jié)果。研究中將土地利用類型與景觀脆弱性相聯(lián)系, 結(jié)合前人的研究成果[19]與研究區(qū)現(xiàn)狀, 分別將土地利用類型賦值如下: 草地-7、耕地-6、未利用地-5、林地-4、水域-3、園地-2、建設(shè)用地-1, 其中草地最為脆弱, 建設(shè)用地最為穩(wěn)定, 通過對(duì)其進(jìn)行分等定級(jí)并標(biāo)準(zhǔn)化處理后作為景觀脆弱度指數(shù)的賦值。

    生態(tài)安全程度的值是由景觀干擾度和景觀脆弱度指數(shù)來表示, 計(jì)算公式如下[17]:

    (14)

    式中,ES為景觀類型生態(tài)安全度,U為景觀干擾度,Q為景觀脆弱度,為景觀安全度,為景觀類型個(gè)數(shù)。

    3 結(jié)果與分析

    3.1 土地利用變化分析

    3.1.1 土地利用結(jié)構(gòu)變化分析

    基于對(duì)烏魯木齊市四期遙感影像進(jìn)行分析, 獲得烏魯木齊市各類土地利用變化情況。通過分析得出, 烏魯木齊市土地利用占比較大的是草地和未利用地, 占比分別為56.22%和28.84%。1990—2017年耕地、林地和草地的面積有所減少, 分別減少了30870.9公頃、9669.24公頃和58065.84公頃, 其中草地面積減少的比例尤為明顯。同時(shí), 園地、水域、建設(shè)用地和未利用地的面積有所增加, 分別增加了455.13公頃、1790.46公頃、27151.47公頃和69208.92公頃。在西部大開發(fā)以及城市化不斷發(fā)展的背景下, 烏魯木齊市建設(shè)用地面積有了明顯的增加, 米東區(qū)荒漠面積較大, 受土地荒漠化的影響, 導(dǎo)致未利用地面積有所增加。

    3.1.2 土地利用時(shí)空變化分析

    對(duì)研究區(qū)土地景觀利用變化分析不僅體現(xiàn)在不同時(shí)段土地景觀類型面積的變化上, 還體現(xiàn)在不同土地景觀利用類型面積的相互轉(zhuǎn)化中。通過ArcGIS軟件將不同時(shí)期土地景觀利用類型圖進(jìn)行疊加分析, 獲得烏魯木齊市土地景觀利用類型轉(zhuǎn)移矩陣, 以此分析烏魯木齊市土地景觀利用時(shí)空變化情況(見表3)。

    通過分析得出烏魯木齊市整體景觀類型面積以耕地、林地、草地景觀類型轉(zhuǎn)出和建設(shè)用地、未利用地景觀類型的轉(zhuǎn)入為主。1990—2000年, 烏魯木齊市耕地、林地和牧草地的轉(zhuǎn)出面積分別為11140.2公頃、1402.56公頃和12339公頃; 建設(shè)用地和未利用地轉(zhuǎn)入面積分別為8945.46公頃和13027.59公頃。2000—2010年, 烏魯木齊市耕地、林地和牧草地的轉(zhuǎn)出面積分別為19160.82公頃、7573.32公頃和52363.26公頃; 建設(shè)用地和未利用地轉(zhuǎn)入面積分別為6576.84公頃和11613.78公頃。2010—2017年, 烏魯木齊市耕地、林地和牧草地的轉(zhuǎn)出面積分別為14167.44公頃、7824.96公頃和36567.81公頃; 建設(shè)用地和未利用地轉(zhuǎn)入面積分別為2191.86公頃和6692.31公頃。

    表2 1990—2017年烏魯木齊市土地利用結(jié)構(gòu)變化

    表3 1990—2017年烏魯木齊市景觀類型轉(zhuǎn)移矩陣(ha)

    通過對(duì)烏魯木齊市土地景觀利用類型面積轉(zhuǎn)移矩陣進(jìn)行分析, 烏魯木齊市土地景觀利用類型轉(zhuǎn)化主要發(fā)生在耕地與草地、耕地與建設(shè)用地、耕地與未利用地、林地與草地、草地與建設(shè)用地、草地與未利用地、建設(shè)用地與未利用地之間相互轉(zhuǎn)化, 1990—2017年, 耕地、林地與草地轉(zhuǎn)出面積最大, 轉(zhuǎn)出面積分別為40060.8公頃、13227.57公頃和81921.24公頃。建設(shè)用地與未利用地轉(zhuǎn)入面積最多, 轉(zhuǎn)入面積分別為30786.84公頃和75324.69公頃。

    在西部大開發(fā)政策影響下, 烏魯木齊市城市化進(jìn)程的速度不斷加快, 建設(shè)用地面積也不斷, 主要是以耕地、草地和未利用地的轉(zhuǎn)入為主。其中, 未利用地轉(zhuǎn)入面積相對(duì)較少, 主要是由于烏魯木齊市是綠洲上的城市, 未利用地大多是沙漠, 可利用面積相對(duì)較少。同時(shí), 未利用地面積有明顯的增加, 絕大部分是由耕地和草地轉(zhuǎn)為未利用地, 在城市快速發(fā)展的背景下, 人口大量向城市流動(dòng), 烏魯木齊市對(duì)耕地和草地過度開發(fā)利用, 城市生態(tài)環(huán)境較為脆弱, 耕地撂荒嚴(yán)重, 過度放牧加劇草地荒漠化, 從而致使未利用地面積增加。總體來說, 烏魯木齊土地景觀利用是以建設(shè)用地和未利用地面積增加, 耕地、林地和草地面積減少為特征。

    3.1.3 土地利用速度變化分析

    不同時(shí)期內(nèi)土地景觀利用類型變化速度存在明顯差異, 由表4分析可知, 烏魯木齊市1990—2017年綜合土地動(dòng)態(tài)度呈上升趨勢(shì), 由3.13%上升至6.57%, 說明烏魯木齊市土地景觀利用年變化率處于一個(gè)較高的水平。其中, 建設(shè)用地的動(dòng)態(tài)度和開發(fā)度均最高, 開發(fā)度由1990—2000年的4.18%上升至2000—2010年的5.24%, 雖然2010—2017年建設(shè)用地的開發(fā)度有所下降, 但與其他土地景觀類型相比仍處于較高水平, 說明隨著城市不斷發(fā)展, 建設(shè)用地增加速度明顯上升。1990—2017年烏魯木齊市耕地、林地和草地的變化速度呈下降趨勢(shì), 且土地利用耗減度大于土地利用開發(fā)度, 說明烏魯木齊市耕地、林地和草地減少的速度大于增加的速度。園地、水域的開發(fā)度都呈上升的趨勢(shì), 說明園地和水域面積的增加速度有了明顯的提升。烏魯木齊市未利用地1990—2010年變化速度和耗減度呈上升趨勢(shì), 2010—2017年變化速度和耗減度呈下降趨勢(shì), 但整體變化速度較慢, 說明雖然未利用地面積得以利用, 但是由于烏魯木齊市未利用地多為沙漠、荒灘, 所以可利用面積較少。

    綜合分析, 烏魯木齊市土地景觀利用類型在研究后期變化幅度趨于平緩, 說明研究后期烏魯木齊市土地景觀利用變化逐漸平穩(wěn)。

    表4 1990—2017年烏魯木齊市土地景觀利用類型速度變化指數(shù)(%)

    注: 表格中代表單一土地動(dòng)態(tài)度;代表綜合土地利用動(dòng)態(tài)度;代表土地利用開發(fā)度;代表土地利用耗減度。

    3.1.4 土地利用程度分析

    由表5可知, 烏魯木齊市1990—2017年土地利用程度由226.53上升至235.08, 說明隨著城市的發(fā)展需求, 對(duì)城市土地開發(fā)和利用程度逐漸增加, 人類活動(dòng)對(duì)土地生態(tài)系統(tǒng)影響也逐漸加深。同時(shí)1990—2017年烏魯木齊市土地利用程度變化量和土地利用程度變化率均大于0, 說明烏魯木齊市土地利用正處于發(fā)展期。

    3.2 景觀格局時(shí)空變化分析

    3.2.1 景觀斑塊類型格局變化分析

    通過對(duì)景觀斑塊類型相關(guān)指數(shù)變化進(jìn)行分析, 有利于了解各景觀類型在城市變化過程中的演變趨勢(shì)和城市景觀整體變化規(guī)律。由圖3可見, 烏魯木齊市27年間城市景觀結(jié)構(gòu)發(fā)生了較大程度上的變化。耕地的斑塊總數(shù)由1990年的2605增加到2017年的11626, 平均斑塊大小由1990年的54.5 ha減少為2017年的9.4 ha, 說明耕地景觀破碎化程度逐漸增大。最大斑塊指數(shù)由1990年的7.67%減少為2017年的1.62%, 說明耕地在城市景觀的優(yōu)勢(shì)度在逐漸降低, 同時(shí)景觀形狀指數(shù)、邊緣密度由1990年的58.33和6.1 ha增加為2017年的114.92和10.66 ha, 說明耕地景觀斑塊形狀更加不規(guī)則以及邊界形狀趨于復(fù)雜。平均形狀指數(shù)、平均斑塊分維度和景觀連通度指數(shù)都存在不同程度的減少, 說明耕地景觀斑塊由聚集趨于分散, 受人類的干擾逐漸增大。

    園地景觀由于所占面積較少, 景觀斑塊個(gè)數(shù)無明顯變化, 但是其平均斑塊大小、景觀連通度由1990年的54 ha和96.24%增加為2017年的82.1 ha和97%, 說明園地景觀破碎化程度逐漸降低, 景觀結(jié)合更為緊密, 同時(shí)平均斑塊分維度更趨近于1, 說明景觀形狀更簡(jiǎn)單, 景觀分布更為集中。

    林地景觀的斑塊總數(shù)由1990年的7630減少為2017年的6531, 同時(shí)林地最大斑塊指數(shù)也有一定程度的增加, 說明林地景觀破碎化程度降低, 景觀優(yōu)勢(shì)度增強(qiáng)。與此同時(shí), 景觀形狀指數(shù)和邊緣密度分別由99.4、6.44 ha減少為93.1和5.8 ha, 說明景觀的自然特征和內(nèi)部生態(tài)過程增強(qiáng), 斑塊復(fù)雜程度降低。景觀斑塊分維指數(shù)和連通度由1.16和97%增加為1.2和98%, 說明原本分離的斑塊之前出現(xiàn)了連接的廊道, 林地趨向于均衡分布的空間結(jié)構(gòu)。

    草地景觀是烏魯木齊市主要景觀之一, 通過分析該景觀斑塊總數(shù)由1990年6632增加為2017年的24490, 同時(shí)平均斑塊大小也由129.07 ha減少32.54 ha, 斑塊總數(shù)和平均斑塊大小都有大幅度的增減和減少, 說明隨著城市的發(fā)展草地景觀破碎化程度呈大幅度地增大。景觀形狀指數(shù)、邊緣密度、平均形狀指數(shù)以及斑塊分維指數(shù)均有所增加, 說明草地景觀斑塊更加不規(guī)則, 景觀形狀趨于復(fù)雜, 景觀受人類干擾水平增加。1990—2017年間草地景觀最大斑塊指數(shù)和景觀連通度分別由56.44%和99.97%減少為51.8%和99.95%, 說明景觀優(yōu)勢(shì)度和聚集度降低, 但是相較于其他景觀草地景觀的優(yōu)勢(shì)度和聚集程度仍處于較高水平。

    水域景觀1990—2017年斑塊總數(shù)由1990年的165增加至2017年的267, 平均斑塊大小和最大斑塊指數(shù)由34.6和0.21%減少至28和0.19%, 說明水域景觀破碎化程度有所增加, 同時(shí)景觀優(yōu)勢(shì)度逐漸降低。邊緣密度、景觀形狀指數(shù)和平均形狀指數(shù)由有所增加, 說明水域景觀斑塊受人類影響更加不規(guī)則和復(fù)雜化, 例如水庫、人造湖等景觀的修建。同時(shí), 景觀斑塊分維指數(shù)和連通度都有所增加, 斑塊分維指數(shù)趨于1, 水域景觀分布較為集中。

    建設(shè)用地景觀斑塊總數(shù)由1990年的4150減少至2017年的2778, 同時(shí)連通度由99.38%增加至99.75%, 說明建設(shè)用地景觀更加聚集。由于受人類活動(dòng)和城市化水平的不斷提升, 導(dǎo)致景觀邊緣密度和平均斑塊大小的值有所增加, 說明景觀破碎化程度隨著城市活動(dòng)的加強(qiáng)逐漸增加。與此同時(shí), 1990—2017年建設(shè)用地景觀形狀指數(shù)、平均形狀指數(shù)以及斑塊分維指數(shù)都有所增加, 說明受人類活動(dòng)干擾, 景觀斑塊更加不規(guī)則和復(fù)雜化。最大斑塊指數(shù)由1990年的1.04%增加至2017年的2.8%, 說明隨著城市的不斷擴(kuò)張和發(fā)展, 建設(shè)用地的景觀優(yōu)勢(shì)度逐漸增加。

    表5 烏魯木齊市1990—2017年土地利用程度變化

    注: NP—斑塊總數(shù), LPI—最大斑塊指數(shù), LSI—景觀形狀指數(shù), ED—邊緣密度, MPS—平均斑塊面積, AWMSI—平均形狀指數(shù), AWMPDF—平均斑塊分維指數(shù), COHESION—景觀連通度。

    Figure 3 Landscape pattern index change of various landscape types in Urumqi from 1990 to 2017

    未利用地景觀斑塊總數(shù)由1990年的4218增加為2017年的13320, 平均斑塊大小由80.65減小為30.72, 說明未利用地景觀破碎化程度不斷加劇。同時(shí)景觀形狀指數(shù), 平均形狀指數(shù)以及邊緣密度值均有所增加, 說明景觀斑塊更加不規(guī)則, 復(fù)雜化程度逐漸增加, 同1990—2017年間景觀斑塊分維指數(shù)由1.2增加至1.3, 說明未利用地受人類的干預(yù)程度較大。同時(shí)最大斑塊指數(shù)由13.62%增加至16%, 說明受人類活動(dòng)的影響, 加之城市生態(tài)環(huán)境脆弱, 導(dǎo)致土地荒漠化程度加劇, 景觀優(yōu)勢(shì)度也有所增加。由于地域的特殊性, 以及受到氣候的影響, 烏魯木齊市可利用土地面積較少, 未利用地多為戈壁和沙漠, 可利用的土地面積較少, 主要集中在烏魯木齊市北部和中部地區(qū)。

    3.2.2 景觀水平格局變化分析

    通過分析表6可知, 烏魯木齊市1990—2017年景觀水平尺度格局發(fā)生明顯的變化。景觀斑塊總數(shù)、斑塊密度、景觀形狀指數(shù)、周長(zhǎng)面積分維數(shù)、景觀分離指數(shù)、Shannon多樣性指數(shù)和均勻度指數(shù)呈上升趨勢(shì), 而最大斑塊指數(shù)、優(yōu)勢(shì)度指數(shù)以及連通度指數(shù)逐漸呈下降趨勢(shì)。這表明了烏魯木齊市景觀破碎化程度逐漸增高, 景觀優(yōu)勢(shì)度不斷降低, 尤其是以草地為主導(dǎo)的景觀格局正在發(fā)生較大變化, 建設(shè)用地不斷擴(kuò)張逐步侵占草地、耕地及其他用地, 景觀異質(zhì)性逐漸增加, 斑塊形狀趨于復(fù)雜化。通過分析得出, 人類活動(dòng)對(duì)烏魯木齊市景觀格局的干擾度不斷增加, 景觀格局逐漸復(fù)雜化, 景觀逐漸趨向離散分布, 結(jié)構(gòu)的不穩(wěn)定性逐漸加大, 景觀自我調(diào)節(jié)的能力不斷降低。目前, 烏魯木齊市的城市化過程逐步加快, 從而對(duì)于城市的可持續(xù)發(fā)展、生態(tài)環(huán)境以及生物多樣性的保護(hù)帶來較大的挑戰(zhàn)。

    3.3 景觀生態(tài)安全分析

    由圖4可知, 烏魯木齊市27年間耕地、草地和未利用地的生態(tài)安全度有所下降, 由1990年的0.746、0.371和0.747下降為2017年的0.449、0.268和0.663, 主要由于城市化過程不斷加快, 建設(shè)用地不斷擴(kuò)張導(dǎo)致耕地、草地面積不斷減少, 土地利用程度不斷提升, 導(dǎo)致耕地、草地等自然景觀受干擾程度增加, 生態(tài)安全度降低。同時(shí), 烏魯木齊市受地域的影響, 生態(tài)環(huán)境較為脆弱, 人為干擾導(dǎo)致城市土地荒漠化程度加劇, 耕地撂荒、草地退化, 使得未利用地面積增加, 生態(tài)安全度降低。林地和水域呈倒“V”字變化, 說明城市自然景觀逐漸被人工景觀所代替, 如人造水庫、園林綠化等工程, 這種人為干擾對(duì)人工景觀的影響是有利的, 但隨著城市化水平的不斷提升, 高效益經(jīng)濟(jì)作物土地種植面積的增加, 以及建設(shè)用地需求的增長(zhǎng), 導(dǎo)致建設(shè)用地和園地的生態(tài)安全度有所增加。

    根據(jù)參考文獻(xiàn)[17-19]并結(jié)合烏魯木齊市實(shí)際研究結(jié)果, 將研究區(qū)景觀格局生態(tài)安全劃為4個(gè)等級(jí): 重警區(qū)[0—0.3]、預(yù)警區(qū)[0.3—0.5]、較安全區(qū)[0.5—0.8]、安全區(qū)[0.8—1.0]。

    通過對(duì)表7、8分析, 27年間, 烏魯木齊市景觀類型安全等級(jí)發(fā)生了較的變化, 重警區(qū)和預(yù)警區(qū)的面積有所增加, 較安全區(qū)和安全區(qū)的面積有所減少。其中重警區(qū)和預(yù)警區(qū)增加的總面積為42671.07 ha。其中, 草地景觀的安全等級(jí)由預(yù)警區(qū)變?yōu)橹鼐瘏^(qū), 耕地景觀的安全等級(jí)由較安全區(qū)變?yōu)轭A(yù)警區(qū), 建設(shè)用地景觀由較安全區(qū)變?yōu)榘踩珔^(qū)。由于人類活動(dòng)的干擾, 加之烏魯木齊市生態(tài)的脆弱性, 導(dǎo)致草地、耕地和建設(shè)用地的安全等級(jí)發(fā)生了明顯變化。

    表6 1990—2017年烏魯木齊市景觀水平格局變化指數(shù)

    注:—斑塊總數(shù),—斑塊密度,—景觀形狀指數(shù),—周長(zhǎng)面積分維數(shù),—最大斑塊指數(shù),—景觀分離度,—景觀連通度,—香濃多樣性指數(shù),—香濃均勻度指數(shù),—景觀優(yōu)勢(shì)度指數(shù)。

    圖4 烏魯木齊市土地景觀類型生態(tài)安全度

    Figure 4 Ecological security of land landscape types in Urumqi

    將各景觀類型生態(tài)安全度與該景觀類型脆弱度的乘積合計(jì), 可得出烏魯木齊市景觀格局生態(tài)安全的整體水平。如圖5所示, 烏魯木齊市景觀生態(tài)安全度呈波動(dòng)下降的趨勢(shì), 由1990年處于較安全區(qū)的0.56下降為2017年處于預(yù)警區(qū)的0.48, 說明烏魯木齊市生態(tài)脆弱性逐漸加深, 生態(tài)安全度有明顯的下降趨勢(shì)。

    4 討論

    1990—2017年間烏魯木齊市土地利用發(fā)生較明顯變化, 園地、水域、建設(shè)用地和未利用地面積呈上升趨勢(shì), 耕地、林地和草地面積呈下降趨勢(shì), 尤其是草地面積下降幅度較大, 所占的比例下降4.09%。同時(shí), 土地景觀類型面積以耕地、林地和草地轉(zhuǎn)出和建設(shè)用地和未利用地轉(zhuǎn)入為主, 烏魯木齊市土地利用正處于發(fā)展期且土地景觀利用年變化率處于一個(gè)較高的水平。隨著城市化進(jìn)程以及經(jīng)濟(jì)快速發(fā)展, 人口不斷增加對(duì)建設(shè)用地的需求也相應(yīng)增加, 經(jīng)濟(jì)的發(fā)展和人口的增加導(dǎo)致耕地、林地和草地面積有所減少。引起景觀相關(guān)指數(shù)變化的主要原因是人類活動(dòng)的干擾, 受干擾程度不同景觀指數(shù)也有不同程度的變化。在景觀斑塊類型格局變化上, 耕地、草地景觀的優(yōu)勢(shì)度逐漸降低, 邊界和形狀逐漸趨于復(fù)雜化, 破碎化程度增加, 連通度下降, 受人類干擾較大, 但草地景觀仍是優(yōu)勢(shì)景觀。與此同時(shí), 受“城市園林綠化工程”等政府政策措施的影響, 林地景觀破碎化程度降低, 景觀優(yōu)勢(shì)度和連通度增加。隨著經(jīng)濟(jì)社會(huì)的不斷發(fā)展, 建設(shè)用地優(yōu)勢(shì)度和聚集度逐漸增加, 受人類活動(dòng)干擾, 景觀斑塊更加不規(guī)則和復(fù)雜化。在城市土地開發(fā)與城鎮(zhèn)發(fā)展的加快的背景下, 人類活動(dòng)加劇了城市土地荒漠化程度, 從而導(dǎo)致城市未利用地景觀優(yōu)勢(shì)度上升。在景觀水平格局變化上, 分析得出人類活動(dòng)對(duì)烏魯木齊市景觀格局的干擾度不斷增加, 景觀格局逐漸復(fù)雜化, 景觀逐漸趨向離散分布, 結(jié)構(gòu)的不穩(wěn)定性逐漸加大, 景觀自我調(diào)節(jié)的能力不斷降低。烏魯木齊市27年間耕地、草地和未利用地的生態(tài)安全度有所下降, 同時(shí)研究區(qū)內(nèi)處于重警區(qū)和預(yù)警區(qū)的面積逐漸增加, 處于較安全區(qū)和安全區(qū)的面積逐漸減少。草地景觀的安全等級(jí)由預(yù)警區(qū)變?yōu)橹鼐瘏^(qū), 耕地景觀的安全等級(jí)由較安全區(qū)變?yōu)轭A(yù)警區(qū), 建設(shè)用地景觀由較安全區(qū)變?yōu)榘踩珔^(qū)。烏魯木齊市總體景觀生態(tài)安全度呈波動(dòng)下降的趨勢(shì), 由0.56下降為0.48, 說明烏魯木齊市生態(tài)脆弱性逐漸加深, 生態(tài)安全度有明顯的下降趨勢(shì)。

    表7 1990—2000年烏魯木齊市土地景觀格局生態(tài)安全預(yù)警

    表8 2010—2017年烏魯木齊市土地景觀格局生態(tài)安全預(yù)警

    圖5 烏魯木齊市整體景觀格局生態(tài)安全度

    Figure 5 Ecological security of the overall landscape pattern in Urumqi

    本文是通過土地利用、景觀格局變化和生態(tài)安全對(duì)烏魯木齊市土地景觀類型變化進(jìn)行梳理和分析, 研究結(jié)果對(duì)烏魯木齊市未來的土地合理利用, 景觀格局優(yōu)化以及城市生態(tài)建設(shè)和可持續(xù)發(fā)展提供一定的參考。但考慮到數(shù)據(jù)的可獲取性本文選擇四個(gè)時(shí)段進(jìn)行分析, 并非持續(xù)性時(shí)間數(shù)據(jù), 獲得的結(jié)果也具有一定的局限性。與此同時(shí), 區(qū)域生態(tài)安全主要是由自然、經(jīng)濟(jì)和社會(huì)等因素共同構(gòu)成的, 景觀格局是其重要影響因素之一, 本文構(gòu)建的景觀生態(tài)安全度分析不具有絕對(duì)性, 未來可以利用數(shù)學(xué)模型或生態(tài)模型等方法構(gòu)建評(píng)價(jià)指標(biāo)體系對(duì)烏魯木齊市生態(tài)安全從自然、經(jīng)濟(jì)和社會(huì)等影因素完善對(duì)烏魯木齊市生態(tài)安全的分析。但是人類的社會(huì)經(jīng)濟(jì)活動(dòng)對(duì)景觀格局的變化會(huì)產(chǎn)生重要的影響, 因此, 本文選用景觀空間結(jié)構(gòu)來研究烏魯木齊市的生態(tài)安全也有重要的意義。

    [1] 張利利, 佘濟(jì)云, 李銳, 等. 1998—2010年五指山市土地利用景觀格局變化分析[J]. 西北林學(xué)院學(xué)報(bào), 2016, 31(1): 221–225.

    [2] 胡曉明. 2011-2016年鞍山市土地利用及景觀格局演變特征[J]. 環(huán)??萍? 2019, 25(2): 44–48.

    [3] 王琦, 王曉娜, 牟惟勇, 等. 黃石市土地利用景觀格局變化研究[J]. 中國農(nóng)學(xué)通報(bào), 2019, 35(9): 81–85.

    [4] 王鵬, 王亞娟, 劉小鵬, 等. 干旱區(qū)生態(tài)移民土地利用景觀格局變化分析—以寧夏紅寺堡區(qū)為例[J]. 干旱區(qū)資源與環(huán)境, 2018, 32(12): 69–74.

    [5] 張旭. 蘭州市土地利用與景觀格局時(shí)空演變動(dòng)態(tài)分析[J]. 甘肅科技, 2018, 34(19): 41–44.

    [6] 王鵬, 王亞娟, 劉小鵬, 等. 寧夏沙坡頭區(qū)土地利用景觀格局變化及其驅(qū)動(dòng)力分析[J]. 西北林學(xué)院學(xué)報(bào), 2018, 33(6): 197–203.

    [7] 王文玲, 高永, 楊光, 等. 基于RS和GIS的武川縣景觀格局變化研究[J]. 北方園藝, 2014(9): 97–100.

    [8] 范麗娟, 田廣星. 1995—2015年上海市土地利用及其景觀格局變化[J]. 水土保持通報(bào), 2018, 38(1): 287–292.

    [9] 韋海航, 鄭芊卉, 莊家堯, 等. 竹山縣土地利用景觀格局動(dòng)態(tài)變化分析[J]. 林業(yè)資源管理, 2018(6): 76–83.

    [10] 郭力宇, 李怡, 郭昭. 寶雞市渭北旱塬區(qū)土地利用及景觀格局動(dòng)態(tài)變化分析[J]. 西北林學(xué)院學(xué)報(bào), 2018, 33(5): 263–268.

    [11] 劉紀(jì)遠(yuǎn), 張?jiān)鱿? 徐新良, 匡文慧, 周萬村, 張樹文, 李仁東, 顏長(zhǎng)珍, 于東升, 吳世新, 江南. 21世紀(jì)初中國土地利用變化的空間格局與驅(qū)動(dòng)力分析. 地理學(xué)報(bào), 2009, 64(12): 1411–1420.

    [12] 佟光臣, 林杰, 陳杭, 等. 1986–2013年南京市土地利用/覆被景觀格局時(shí)空變化及驅(qū)動(dòng)力因素分析[J]. 水土保持研究, 2017, 24(2): 240–245.

    [13] 徐允, 張培, 周寶同. 基于Shannon多樣性t檢驗(yàn)法的土地利用景觀分區(qū)研究—以重慶市永川區(qū)為例[J]. 中國農(nóng)業(yè)資源與區(qū)劃, 2019, 40(1): 134–141+151.

    [14] 馬明, 顧康康, 李詠. 基于生態(tài)安全格局的城鄉(xiāng)生態(tài)空間布局與優(yōu)化—以宣城市為例[J]. 中國農(nóng)業(yè)資源與區(qū)劃, 2019, 40(4): 93–102.

    [15] 冉晨, 白曉永, 譚秋, 等. 典型喀斯特地區(qū)石漠化景觀格局對(duì)土地利用變化的響應(yīng)[J]. 生態(tài)學(xué)報(bào), 2018, 38(24): 8901–8910.

    [16] 哈孜亞·包浪提將, 毋兆鵬, 陳學(xué)剛, 等. 烏魯木齊市景觀格局變化及驅(qū)動(dòng)力分析[J]. 生態(tài)科學(xué), 2018, 37(1): 62–70.

    [17] 梁二敏, 張軍民. 新疆瑪納斯河流域景觀格局變化的生態(tài)安全分析[J]. 水土保持研究, 2016, 23(3): 170–175.

    [18] 郝彩蓮, 王勇, 肖偉華, 尹軍. 基于景觀格局的承德市武烈河流域生態(tài)安全分析. 南水北調(diào)與水利科技, 2012, 10(5): 67–71.

    [19] 李曉燕, 張樹文. 基于景觀結(jié)構(gòu)的吉林西部生態(tài)安全動(dòng)態(tài)分析. 干旱區(qū)研究, 2005(1): 57–62.

    [20] 王瑋, 馮琦勝, 郭鈮, 等. 基于長(zhǎng)時(shí)間序列NDVI資料的我國西北干旱區(qū)植被覆蓋動(dòng)態(tài)監(jiān)測(cè)[J]. 草業(yè)科學(xué), 2015, 32(12): 1969–1979.

    [21] 張濤, 張穎, 楊力鵬, 等. 內(nèi)蒙古自治區(qū)呼和浩特市1990—2010年景觀格局變化及其驅(qū)動(dòng)力[J]. 水土保持通報(bào), 2018, 38(2): 217–222.

    Analysis on landscape pattern change and ecological security of land use in Urumqi

    ZHU Yanan1, PU Chunling2,*

    1. School of management of Xinjiang Agricultural University, Urumqi 830052, China 2. Research Center for Economic and Social Development, Xinjiang Agricultural University, Urumqi 830052, China

    Based on the analysis of the spatial and temporal changes of land use landscape pattern and ecological security in Urumqi from 1990 to 2017, using the remote sensing images of Urumqi in 1990, 2000, 2010 and 2017 as the data source, the remote sensing images of Urumqi were interpreted and analyzed by ERDAS and ArcGIS software, and the landscape pattern index of Urumqi was obtained by FRAGSTAT software. Combined with the theory of landscape ecology, the changes of land use, landscape pattern and ecological security in Urumqi were analyzed. The results showed that: (1) The land use/cover change was obvious in Urumqi; the grassland landscape was dominant, accounting for 56.22% of the total land use area, but the dominance of grassland landscape decreased gradually; the unused land area was second only to grassland, accounting for 28.84% of the total land use area. At the same time, the transfer of cultivated land, woodland and grassland landscape and the transfer of construction land and unused land landscape were the main landscape types. (2) The change rate of land use in Urumqi was gradually accelerated from 3.13% to 6.57%, and the degree of land use was increased from 226.53 in 1990 to 235.08 in 2017, which indicated that the impact of human activities on land ecosystem was gradually deepened. (3) The change of landscape pattern in Urumqi was obvious, the degree of landscape fragmentation was gradually deepened, the landscape patches were more irregular and complicated, and the diversity and instability of landscape structure were gradually increased. (4) The ecological security degree of Urumqi decreased from 0.56 in 1990 to 0.48 in 2017, which indicated that the ecological vulnerability of Urumqi was gradually deepened. The results of this study will provide some reference for the future rational land use, landscape pattern optimization, urban ecological construction and sustainable development in Urumqi.

    land use; landscape pattern; ecological security; Urumqi

    10.14108/j.cnki.1008-8873.2020.02.017

    S731.2

    A

    1008-8873(2020)02-133-12

    2019-06-18;

    2019-07-31

    國家社會(huì)科學(xué)基金項(xiàng)目(14XGL005); 國家社會(huì)科學(xué)基金項(xiàng)目(08BJY061)

    朱亞楠(1995—), 女, 新疆烏魯木齊人, 碩士研究生, 主要從事土地資源利用研究, E-mail:727526126@qq.com

    蒲春玲(1961—), 女, 新疆阿克蘇人, 教授, 主要從事土地資源可持續(xù)利用和區(qū)域經(jīng)濟(jì)研究, E-mail:puchunling@163.com

    朱亞楠, 蒲春玲. 烏魯木齊市土地利用景觀格局變化及生態(tài)安全分析[J]. 生態(tài)科學(xué), 2020, 39(2): 133–144.

    ZHU Yanan, PU Chunling. Analysis on landscape pattern change and ecological security of land use in Urumqi[J]. Ecological Science, 2020, 39(2): 133–144.

    猜你喜歡
    烏魯木齊市格局土地利用
    格局
    聯(lián)手共建 努力打造大調(diào)解工作格局
    烏魯木齊市園林綠化養(yǎng)護(hù)管理存在問題及對(duì)策
    土地利用生態(tài)系統(tǒng)服務(wù)研究進(jìn)展及啟示
    近5年烏魯木齊市PM2.5變化分析
    烏魯木齊市土地征收與融資問題探討
    濱??h土地利用挖潛方向在哪里
    烏魯木齊市將建立報(bào)廢汽車管理長(zhǎng)效機(jī)制
    小人物的大格局
    小說月刊(2015年12期)2015-04-23 08:51:10
    論低碳經(jīng)濟(jì)與轉(zhuǎn)變土地利用方式
    亚洲欧美精品自产自拍| 欧美成人午夜免费资源| 好男人视频免费观看在线| 国产免费一区二区三区四区乱码| 美女福利国产在线| 亚洲精品乱久久久久久| 国产一区亚洲一区在线观看| 国产精品秋霞免费鲁丝片| 国产精品无大码| 香蕉国产在线看| 亚洲精品视频女| 天天操日日干夜夜撸| 午夜福利网站1000一区二区三区| 久久久国产精品麻豆| 成人免费观看视频高清| 蜜桃在线观看..| 伦精品一区二区三区| 少妇高潮的动态图| 色网站视频免费| 一级毛片黄色毛片免费观看视频| videosex国产| 精品人妻熟女毛片av久久网站| 国产成人午夜福利电影在线观看| av视频免费观看在线观看| 日韩欧美一区视频在线观看| 精品视频人人做人人爽| 精品一区二区三卡| 亚洲少妇的诱惑av| 亚洲国产精品国产精品| 内地一区二区视频在线| 午夜福利视频在线观看免费| 十八禁高潮呻吟视频| 国产亚洲精品第一综合不卡 | 成人毛片a级毛片在线播放| 美女大奶头黄色视频| 国产精品久久久久久久电影| 女性生殖器流出的白浆| 色94色欧美一区二区| 99国产综合亚洲精品| 免费女性裸体啪啪无遮挡网站| 精品国产露脸久久av麻豆| 在线精品无人区一区二区三| 午夜福利,免费看| 久久亚洲国产成人精品v| 国产成人一区二区在线| 人人妻人人爽人人添夜夜欢视频| 午夜福利,免费看| 另类亚洲欧美激情| 成人综合一区亚洲| 两个人免费观看高清视频| 两个人免费观看高清视频| 国产成人一区二区在线| 大香蕉久久网| 九九在线视频观看精品| 欧美日韩一区二区视频在线观看视频在线| 香蕉丝袜av| 国产有黄有色有爽视频| 人成视频在线观看免费观看| freevideosex欧美| 成人黄色视频免费在线看| 国产精品国产三级专区第一集| 少妇被粗大猛烈的视频| 我的女老师完整版在线观看| 亚洲欧美成人综合另类久久久| 国产精品国产三级国产av玫瑰| 免费久久久久久久精品成人欧美视频 | 久久久欧美国产精品| 久久久欧美国产精品| 岛国毛片在线播放| 全区人妻精品视频| av.在线天堂| 视频在线观看一区二区三区| 国产一区二区在线观看日韩| 汤姆久久久久久久影院中文字幕| 男女啪啪激烈高潮av片| 成人手机av| 极品人妻少妇av视频| 日韩中文字幕视频在线看片| 丰满乱子伦码专区| 成人漫画全彩无遮挡| 日本av手机在线免费观看| 9191精品国产免费久久| 老司机亚洲免费影院| 三级国产精品片| 内地一区二区视频在线| 国产色婷婷99| 久久av网站| 国产男女超爽视频在线观看| 亚洲一码二码三码区别大吗| 永久免费av网站大全| 香蕉国产在线看| 一区二区av电影网| 亚洲内射少妇av| 伦理电影免费视频| 国产1区2区3区精品| 亚洲精品日本国产第一区| 日韩欧美精品免费久久| 我的女老师完整版在线观看| 中国三级夫妇交换| 久久久久人妻精品一区果冻| xxxhd国产人妻xxx| 在线观看免费高清a一片| 国语对白做爰xxxⅹ性视频网站| 下体分泌物呈黄色| 中文字幕亚洲精品专区| 男人添女人高潮全过程视频| 亚洲欧美成人精品一区二区| 欧美激情极品国产一区二区三区 | 男女边吃奶边做爰视频| 亚洲国产精品一区三区| 中国国产av一级| 91精品伊人久久大香线蕉| www日本在线高清视频| 黄色视频在线播放观看不卡| 久久久精品区二区三区| 精品亚洲成a人片在线观看| 成年女人在线观看亚洲视频| 97超碰精品成人国产| 99视频精品全部免费 在线| 欧美精品人与动牲交sv欧美| 欧美bdsm另类| 国产亚洲欧美精品永久| 丰满乱子伦码专区| 亚洲高清免费不卡视频| 欧美国产精品va在线观看不卡| 久久久精品区二区三区| 乱人伦中国视频| 久热这里只有精品99| 国产探花极品一区二区| 高清黄色对白视频在线免费看| 2022亚洲国产成人精品| 欧美人与性动交α欧美精品济南到 | 一区二区三区乱码不卡18| 日韩 亚洲 欧美在线| 国产欧美日韩一区二区三区在线| 天天影视国产精品| 深夜精品福利| 免费观看在线日韩| 国产精品偷伦视频观看了| 精品一区二区三区四区五区乱码 | 91午夜精品亚洲一区二区三区| 两性夫妻黄色片 | 日本av手机在线免费观看| 亚洲四区av| 国产精品久久久久久久久免| 在线观看人妻少妇| 久久亚洲国产成人精品v| 精品熟女少妇av免费看| 亚洲精品日韩在线中文字幕| 少妇人妻 视频| 久久久国产一区二区| 久久韩国三级中文字幕| 性色av一级| 日本欧美视频一区| 捣出白浆h1v1| 在线看a的网站| 亚洲欧洲国产日韩| 精品一区二区三区四区五区乱码 | 成人毛片60女人毛片免费| 欧美日韩视频精品一区| 七月丁香在线播放| 亚洲国产精品一区三区| 观看美女的网站| 欧美精品亚洲一区二区| 精品99又大又爽又粗少妇毛片| 国产色爽女视频免费观看| 久久久国产精品麻豆| 亚洲av国产av综合av卡| 少妇被粗大猛烈的视频| 成人影院久久| 午夜老司机福利剧场| 亚洲精品第二区| 国产亚洲精品第一综合不卡 | 久久人人爽人人爽人人片va| 91久久精品国产一区二区三区| 少妇人妻 视频| 又大又黄又爽视频免费| 免费观看a级毛片全部| 激情五月婷婷亚洲| 精品久久久精品久久久| 欧美亚洲日本最大视频资源| 日韩精品免费视频一区二区三区 | 男女下面插进去视频免费观看 | 成人二区视频| h视频一区二区三区| 99精国产麻豆久久婷婷| 99香蕉大伊视频| 黑人猛操日本美女一级片| 一区二区三区乱码不卡18| 纯流量卡能插随身wifi吗| 中国美白少妇内射xxxbb| 免费观看在线日韩| 最新中文字幕久久久久| 欧美性感艳星| 色婷婷久久久亚洲欧美| 日日啪夜夜爽| 日韩,欧美,国产一区二区三区| 天天躁夜夜躁狠狠久久av| 久久av网站| 成年美女黄网站色视频大全免费| 五月伊人婷婷丁香| 亚洲内射少妇av| 久久午夜福利片| 久久精品久久久久久久性| 一区二区日韩欧美中文字幕 | 日本与韩国留学比较| 日韩三级伦理在线观看| 精品卡一卡二卡四卡免费| 国产日韩欧美视频二区| 美女中出高潮动态图| 欧美xxⅹ黑人| 精品国产乱码久久久久久小说| av天堂久久9| 中文乱码字字幕精品一区二区三区| 少妇精品久久久久久久| 久久精品国产亚洲av涩爱| 精品福利永久在线观看| 成人国产av品久久久| 国产视频首页在线观看| 女性被躁到高潮视频| 男人舔女人的私密视频| 久久久久精品人妻al黑| 在线观看免费日韩欧美大片| 九色成人免费人妻av| 国产精品熟女久久久久浪| 国产亚洲av片在线观看秒播厂| 免费大片黄手机在线观看| 视频在线观看一区二区三区| 日韩中文字幕视频在线看片| 亚洲婷婷狠狠爱综合网| 国产成人欧美| 人妻一区二区av| 1024视频免费在线观看| 丰满乱子伦码专区| 亚洲欧美中文字幕日韩二区| 国产日韩欧美在线精品| 国产男女超爽视频在线观看| 日本vs欧美在线观看视频| 国产精品女同一区二区软件| 国产黄色视频一区二区在线观看| 一级爰片在线观看| 插逼视频在线观看| 国产精品秋霞免费鲁丝片| 卡戴珊不雅视频在线播放| 如何舔出高潮| 欧美激情国产日韩精品一区| 欧美日韩成人在线一区二区| 青青草视频在线视频观看| 啦啦啦啦在线视频资源| 人妻一区二区av| 久久久亚洲精品成人影院| 国产男人的电影天堂91| 人妻人人澡人人爽人人| 日韩制服丝袜自拍偷拍| 亚洲经典国产精华液单| 国产精品偷伦视频观看了| 国产高清国产精品国产三级| 下体分泌物呈黄色| 精品人妻在线不人妻| 人成视频在线观看免费观看| 高清av免费在线| 午夜福利影视在线免费观看| 老女人水多毛片| 亚洲av日韩在线播放| 精品久久久久久电影网| 国产1区2区3区精品| 纯流量卡能插随身wifi吗| 人人妻人人爽人人添夜夜欢视频| av在线观看视频网站免费| 青春草亚洲视频在线观看| 黄网站色视频无遮挡免费观看| 蜜桃在线观看..| 国产精品成人在线| 亚洲欧洲精品一区二区精品久久久 | 美女国产高潮福利片在线看| 国产一区有黄有色的免费视频| 日本-黄色视频高清免费观看| 国产精品免费大片| 国产乱来视频区| 久久av网站| 一区在线观看完整版| 看免费av毛片| 晚上一个人看的免费电影| 国产一区二区激情短视频 | 日本91视频免费播放| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av福利一区| 国语对白做爰xxxⅹ性视频网站| 极品少妇高潮喷水抽搐| 日本-黄色视频高清免费观看| 大码成人一级视频| 在现免费观看毛片| 久久精品熟女亚洲av麻豆精品| 男的添女的下面高潮视频| av天堂久久9| 国产亚洲一区二区精品| 精品少妇黑人巨大在线播放| 免费在线观看黄色视频的| 18禁在线无遮挡免费观看视频| 欧美日韩亚洲高清精品| 女人久久www免费人成看片| 极品人妻少妇av视频| 国产成人a∨麻豆精品| 高清在线视频一区二区三区| 日本免费在线观看一区| 在线观看三级黄色| 久久久久精品性色| 国产又爽黄色视频| 欧美精品一区二区大全| 国产片内射在线| 亚洲精品中文字幕在线视频| 不卡视频在线观看欧美| av福利片在线| 性色avwww在线观看| 精品一区二区三卡| 这个男人来自地球电影免费观看 | 午夜免费鲁丝| 妹子高潮喷水视频| 亚洲中文av在线| 一级毛片 在线播放| 久久久久网色| 免费在线观看完整版高清| 国产一区二区在线观看日韩| 精品卡一卡二卡四卡免费| 国产一区有黄有色的免费视频| 亚洲欧美精品自产自拍| 国产精品免费大片| 亚洲成人一二三区av| 成人午夜精彩视频在线观看| 不卡视频在线观看欧美| 男女免费视频国产| 伦理电影免费视频| 国产欧美日韩一区二区三区在线| 亚洲精品美女久久久久99蜜臀 | 欧美国产精品va在线观看不卡| 国产深夜福利视频在线观看| 国产深夜福利视频在线观看| 久久精品久久精品一区二区三区| 欧美97在线视频| 日本猛色少妇xxxxx猛交久久| 黑人猛操日本美女一级片| 亚洲精品成人av观看孕妇| 久久久久国产网址| 免费播放大片免费观看视频在线观看| 亚洲婷婷狠狠爱综合网| 国产又爽黄色视频| 视频在线观看一区二区三区| 人人妻人人添人人爽欧美一区卜| 日韩制服丝袜自拍偷拍| 老司机影院毛片| 丰满乱子伦码专区| 香蕉精品网在线| 插逼视频在线观看| 成人影院久久| 亚洲国产成人一精品久久久| 中文字幕人妻丝袜制服| 国产片内射在线| 国产精品久久久久久久久免| 老司机影院毛片| 97在线视频观看| 成人二区视频| 国产免费现黄频在线看| 大片电影免费在线观看免费| 成人二区视频| 国产有黄有色有爽视频| 夫妻午夜视频| 成年av动漫网址| 大话2 男鬼变身卡| 黑人欧美特级aaaaaa片| 久久久久久人妻| 满18在线观看网站| 青青草视频在线视频观看| 精品亚洲成a人片在线观看| 国产免费一级a男人的天堂| 免费大片18禁| 老女人水多毛片| 91国产中文字幕| 国产成人a∨麻豆精品| 国产免费又黄又爽又色| 男女边摸边吃奶| 亚洲,一卡二卡三卡| 免费人成在线观看视频色| 一区在线观看完整版| 热99国产精品久久久久久7| 久久精品国产亚洲av涩爱| 亚洲精品乱久久久久久| 一二三四在线观看免费中文在 | 久久99精品国语久久久| 久久韩国三级中文字幕| 考比视频在线观看| av播播在线观看一区| 亚洲四区av| 国产在线视频一区二区| 黑人高潮一二区| 亚洲av综合色区一区| 国产成人精品久久久久久| 日韩成人av中文字幕在线观看| 老女人水多毛片| 国产精品偷伦视频观看了| 亚洲,一卡二卡三卡| 精品人妻在线不人妻| 免费在线观看黄色视频的| av女优亚洲男人天堂| 亚洲国产日韩一区二区| 亚洲精品乱久久久久久| 欧美 日韩 精品 国产| 制服丝袜香蕉在线| 久久精品久久精品一区二区三区| 成人综合一区亚洲| 日韩伦理黄色片| 侵犯人妻中文字幕一二三四区| 黑丝袜美女国产一区| 欧美国产精品一级二级三级| 日韩制服丝袜自拍偷拍| 日产精品乱码卡一卡2卡三| 高清欧美精品videossex| 欧美激情国产日韩精品一区| av免费在线看不卡| 丰满饥渴人妻一区二区三| 成年美女黄网站色视频大全免费| 精品少妇久久久久久888优播| 又黄又爽又刺激的免费视频.| 街头女战士在线观看网站| 热99国产精品久久久久久7| www.av在线官网国产| 制服人妻中文乱码| 欧美精品高潮呻吟av久久| 久久久亚洲精品成人影院| 18禁在线无遮挡免费观看视频| 91精品三级在线观看| 国产亚洲欧美精品永久| 欧美亚洲日本最大视频资源| 精品福利永久在线观看| 精品国产国语对白av| 欧美人与善性xxx| 中文精品一卡2卡3卡4更新| 久久精品国产鲁丝片午夜精品| 久久99蜜桃精品久久| 国产老妇伦熟女老妇高清| 成人亚洲精品一区在线观看| 久久久久人妻精品一区果冻| 午夜影院在线不卡| 免费播放大片免费观看视频在线观看| 韩国av在线不卡| 亚洲精品aⅴ在线观看| 性色avwww在线观看| 免费av不卡在线播放| 精品少妇久久久久久888优播| 成年女人在线观看亚洲视频| 99热网站在线观看| 亚洲精品一区蜜桃| 成人黄色视频免费在线看| 一级爰片在线观看| 亚洲国产最新在线播放| 久久久久久久久久成人| 自拍欧美九色日韩亚洲蝌蚪91| 如何舔出高潮| 18在线观看网站| 我要看黄色一级片免费的| 欧美激情极品国产一区二区三区 | 久久久久久人妻| 母亲3免费完整高清在线观看 | 午夜激情久久久久久久| 婷婷成人精品国产| 2022亚洲国产成人精品| 在线 av 中文字幕| 99视频精品全部免费 在线| av国产精品久久久久影院| 亚洲精华国产精华液的使用体验| 精品午夜福利在线看| av黄色大香蕉| 亚洲精品久久久久久婷婷小说| 曰老女人黄片| 亚洲第一区二区三区不卡| 9191精品国产免费久久| 黄色配什么色好看| 天美传媒精品一区二区| 久久精品熟女亚洲av麻豆精品| 国产精品国产三级国产av玫瑰| 18禁观看日本| 免费观看a级毛片全部| 久久99热6这里只有精品| 国产精品 国内视频| 中文乱码字字幕精品一区二区三区| 欧美日韩综合久久久久久| videossex国产| 五月开心婷婷网| 一区二区三区四区激情视频| 亚洲精华国产精华液的使用体验| 少妇的逼好多水| 日本91视频免费播放| 亚洲av男天堂| 国产乱人偷精品视频| 国产欧美日韩综合在线一区二区| 欧美丝袜亚洲另类| 国产精品成人在线| 免费黄色在线免费观看| 亚洲精品自拍成人| 黄色视频在线播放观看不卡| 久久人妻熟女aⅴ| 一区二区三区乱码不卡18| 国产男女超爽视频在线观看| 毛片一级片免费看久久久久| www.av在线官网国产| 亚洲,一卡二卡三卡| 人人妻人人爽人人添夜夜欢视频| 亚洲精品国产色婷婷电影| 免费少妇av软件| av福利片在线| 欧美日韩亚洲高清精品| 日产精品乱码卡一卡2卡三| 最近的中文字幕免费完整| 日韩一本色道免费dvd| 成年人午夜在线观看视频| 永久网站在线| 一本大道久久a久久精品| 免费av不卡在线播放| 亚洲成人手机| 女人久久www免费人成看片| av天堂久久9| 久久久精品免费免费高清| 少妇的逼好多水| 欧美精品高潮呻吟av久久| 最近最新中文字幕免费大全7| av黄色大香蕉| 亚洲国产av影院在线观看| 成年av动漫网址| 久久精品久久久久久噜噜老黄| 熟妇人妻不卡中文字幕| 亚洲av成人精品一二三区| 男人舔女人的私密视频| 草草在线视频免费看| 欧美成人午夜免费资源| 亚洲欧洲国产日韩| 飞空精品影院首页| 日本色播在线视频| 国产成人欧美| 看非洲黑人一级黄片| 国产精品久久久久久久久免| 国产一区二区三区av在线| 国产无遮挡羞羞视频在线观看| 免费高清在线观看视频在线观看| 三级国产精品片| 亚洲激情五月婷婷啪啪| 久久精品夜色国产| 亚洲av国产av综合av卡| 97精品久久久久久久久久精品| 男女啪啪激烈高潮av片| 国产日韩欧美视频二区| 久久99精品国语久久久| 99热全是精品| 亚洲精品,欧美精品| 侵犯人妻中文字幕一二三四区| 欧美亚洲 丝袜 人妻 在线| 激情视频va一区二区三区| 亚洲精品自拍成人| 少妇 在线观看| 国产精品偷伦视频观看了| 亚洲精品中文字幕在线视频| 国产精品一国产av| 国产欧美日韩一区二区三区在线| 婷婷色综合大香蕉| 亚洲国产欧美日韩在线播放| 91aial.com中文字幕在线观看| 99热这里只有是精品在线观看| 国产精品一国产av| 18禁观看日本| 亚洲av男天堂| 国产xxxxx性猛交| 欧美国产精品一级二级三级| 少妇被粗大的猛进出69影院 | 中国国产av一级| 99re6热这里在线精品视频| 亚洲精品久久午夜乱码| 岛国毛片在线播放| 91精品伊人久久大香线蕉| 99精国产麻豆久久婷婷| 国产亚洲精品久久久com| 人成视频在线观看免费观看| 久久精品夜色国产| 日本与韩国留学比较| 黄色 视频免费看| 亚洲av成人精品一二三区| 亚洲少妇的诱惑av| 各种免费的搞黄视频| 在线观看一区二区三区激情| 亚洲内射少妇av| 亚洲人成网站在线观看播放| 黑人欧美特级aaaaaa片| 黄色配什么色好看| 少妇猛男粗大的猛烈进出视频| 亚洲人成77777在线视频| 美国免费a级毛片| 满18在线观看网站| 黄色视频在线播放观看不卡| 亚洲一区二区三区欧美精品| 三级国产精品片| av免费观看日本| 精品人妻一区二区三区麻豆| 又粗又硬又长又爽又黄的视频| 亚洲人成网站在线观看播放| 午夜免费鲁丝| 老司机影院毛片| www日本在线高清视频| 国产精品 国内视频| 啦啦啦啦在线视频资源| 26uuu在线亚洲综合色| av在线app专区| 黑人高潮一二区| 国产不卡av网站在线观看| 久久人人爽av亚洲精品天堂| 少妇高潮的动态图| 最近手机中文字幕大全| 男女下面插进去视频免费观看 | 边亲边吃奶的免费视频|