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

    廣西防城金花茶國家級(jí)自然保護(hù)區(qū)建立對(duì)區(qū)域土地利用動(dòng)態(tài)的影響

    2022-06-29 01:53:32黃秋霞廖南燕楊元征俎佳星王嘉麗蔡文華蘇宏新
    熱帶地理 2022年6期
    關(guān)鍵詞:金花核心區(qū)保護(hù)區(qū)

    黃秋霞,廖南燕,楊元征,3,4,俎佳星,3,4,王嘉麗,3,4,蔡文華,3,4,楊 健,3,4,蘇宏新,3,4,楊 婷

    (1. 南寧師范大學(xué)地理科學(xué)與規(guī)劃學(xué)院,南 寧530001;2. 廣西壯族自治區(qū)防城金花茶國家級(jí)自然保護(hù)區(qū)管理處,廣西 防城港 538021;3. 北部灣環(huán)境演變與資源利用教育部重點(diǎn)實(shí)驗(yàn)室,南寧 530001;4. 廣西地表過程與智能模擬重點(diǎn)實(shí)驗(yàn)室,南寧 530001)

    自然保護(hù)區(qū)的建立是保護(hù)珍稀動(dòng)植物及其生境的重要舉措(Maiorano et al., 2008; Lung et al.,2010; Radeloff et al., 2010),其在緩解土地利用變化、保護(hù)生態(tài)多樣性方面起積極作用(Stein et al.,2008;Xu et al.,2017)。自1956年中國建立第一個(gè)自然保護(hù)區(qū)(鼎湖山國家級(jí)自然保護(hù)區(qū),面積為11.33 km2)以來,自然保護(hù)區(qū)數(shù)量和面積不斷增加,建設(shè)體系日益成熟(張鐿鋰等,2015),了解自然保護(hù)區(qū)是否有效實(shí)現(xiàn)保護(hù)目的逐漸成為熱點(diǎn)話題,但其建設(shè)成效評(píng)估體系仍不完善(祝萍等,2018)。目前,評(píng)價(jià)保護(hù)區(qū)建設(shè)成效的方法很多,至今尚未形成統(tǒng)一標(biāo)準(zhǔn)。土地利用變化能夠直觀表現(xiàn)自然保護(hù)區(qū)保護(hù)成效(郭子良等,2021),但簡(jiǎn)單分析保護(hù)區(qū)成立以來的土地利用類型面積變化并不能達(dá)到有效評(píng)價(jià)的目的。從長(zhǎng)時(shí)間序列深入了解保護(hù)區(qū)及其鄰近區(qū)土地利用類型動(dòng)態(tài)演變,可以綜合評(píng)估保護(hù)區(qū)建設(shè)成效,為土地資源規(guī)劃提供決策支撐,對(duì)維持生態(tài)系統(tǒng)平衡和資源環(huán)境的可持續(xù)發(fā)展具有重要意義(Lambert et al.,2017)。

    土地利用是鏈接人類社會(huì)經(jīng)濟(jì)活動(dòng)與自然生態(tài)過程的紐帶(Liu et al.,2014),其利用方式轉(zhuǎn)變對(duì)生物多樣性(García et al., 2007)、生態(tài)系統(tǒng)服務(wù)(田義超等,2015;劉曉娟等,2019;Chen et al.,2020; Cao et al., 2021)、資源的可持續(xù)利用(姚靜韜等,2020)以及氣候變化都有深遠(yuǎn)影響(Bu et al.,2014;Huang et al.,2015)。近年來,國內(nèi)外學(xué)者對(duì)自然保護(hù)區(qū)土地利用變化做了大量研究,主要利用地理信息系統(tǒng)技術(shù),結(jié)合數(shù)理統(tǒng)計(jì)與景觀生態(tài)學(xué)方法分析保護(hù)區(qū)建立后土地利用和景觀格局動(dòng)態(tài)變化(Zhang et al.,2017;張宏靜等,2021),探討造成該動(dòng)態(tài)變化的驅(qū)動(dòng)因素(Su et al., 2015;滿蘇爾·沙比提等,2016,2017),并預(yù)測(cè)未來土地利用變化趨勢(shì)及其對(duì)生態(tài)系統(tǒng)服務(wù)的影響(井云清等,2016;Zhang et al.,2019;Li et al.,2021)。但目前研究較多集中在保護(hù)區(qū)內(nèi)部,較少針對(duì)保護(hù)區(qū)內(nèi)部與外部土地利用動(dòng)態(tài)差異。

    廣西防城金花茶國家級(jí)自然保護(hù)區(qū)(簡(jiǎn)稱“金花茶保護(hù)區(qū)”)位于廣西壯族自治區(qū)防城港市境內(nèi),于1986年經(jīng)廣西壯族自治區(qū)人民政府批準(zhǔn),建立廣西防城上岳金花茶自治區(qū)級(jí)自然保護(hù)區(qū),1994年晉升為國家級(jí)自然保護(hù)區(qū),主要保護(hù)對(duì)象為珍稀瀕危金花茶組植物及其生境,是中國乃至世界金花茶植物的分布中心。目前對(duì)金花茶保護(hù)區(qū)的研究主要集中在植物多樣性(龐潔等,2008)、植物資源保護(hù)及開發(fā)(謝釗等,2019)、昆蟲區(qū)系分析(莫璐瑕等,2014)等方面,對(duì)其土地利用變化的研究比較少。保護(hù)區(qū)的建立使人類活動(dòng)逐漸減少,對(duì)金花茶生境質(zhì)量提升起一定的積極作用,但在當(dāng)?shù)亟?jīng)濟(jì)發(fā)展的需求下,部分人類活動(dòng)對(duì)金花茶保護(hù)區(qū)土地利用格局仍然造成影響(黃曉園等,2019)。保護(hù)區(qū)及鄰近區(qū)域土地利用變化會(huì)威脅保護(hù)區(qū)生態(tài)安全與生態(tài)系統(tǒng)服務(wù)功能。因而,有必要探索保護(hù)區(qū)的建立對(duì)區(qū)域土地利用產(chǎn)生的動(dòng)態(tài)影響。鑒于此,以1986—2019 年7 期Landsat5 TM/Landsat8 OLI 遙感影像為數(shù)據(jù)源,采用決策樹分類方法進(jìn)行土地利用分類,探討金花茶保護(hù)區(qū)非核心區(qū)和保護(hù)區(qū)外部2 km鄰近區(qū)的土地利用類型動(dòng)態(tài)演變規(guī)律。以期為保護(hù)區(qū)的建設(shè)及監(jiān)管、資源開發(fā)及生態(tài)安全保護(hù)提供一定的數(shù)據(jù)支撐。

    1 研究區(qū)概況與方法

    1.1 研究區(qū)概況

    金花茶保護(hù)區(qū)位于廣西壯族自治區(qū)防城港市防城區(qū)境內(nèi)(李廣超等,2021),總面積為9 098.62 hm2,其核心區(qū)面積為1 479.10 hm2,非核心區(qū)面積為7 619.52 hm2,非核心區(qū)包括保護(hù)區(qū)內(nèi)實(shí)驗(yàn)區(qū)和緩沖區(qū)兩部分,其中緩沖區(qū)面積為3 459.21 hm2,實(shí)驗(yàn)區(qū)面積為4 160.31 hm2,跨越那梭、華石、大菉和扶隆4個(gè)鄉(xiāng)鎮(zhèn),地理位置為21°44′-21°48′N、108°03′-108°12′E(圖1),保護(hù)區(qū)內(nèi)人口約1 900人(實(shí)地走訪調(diào)查所得),糧食作物以玉米和水稻為主,經(jīng)濟(jì)作物主要有紅薯、蔬菜、八角、肉桂。保護(hù)區(qū)地處十萬大山山脈,屬熱帶季風(fēng)氣候,全年氣候溫和濕潤,冬短夏長(zhǎng)(李廣超等,2021),年均溫約21.8℃,年均降雨量高達(dá)2 900 mm。植被類型以闊葉混交林、亞熱帶常綠闊葉林為主。為更好地評(píng)估保護(hù)區(qū)建設(shè)與發(fā)展對(duì)當(dāng)?shù)厣鷳B(tài)環(huán)境的保護(hù)力度,將研究區(qū)劃分為保護(hù)區(qū)內(nèi)的核心區(qū)(簡(jiǎn)稱“核心區(qū)”)、保護(hù)區(qū)非核心區(qū)(簡(jiǎn)稱“非核心區(qū)”),保護(hù)區(qū)外部鄰近區(qū)(簡(jiǎn)稱“鄰近區(qū)”)三部分(見圖1)。其中,鄰近區(qū)為在保護(hù)區(qū)外部建立2 km 的緩沖區(qū),外部緩沖區(qū)面積與保護(hù)區(qū)面積大致相同。通過對(duì)比保護(hù)區(qū)外部鄰近區(qū)與非核心區(qū)土地利用變化趨勢(shì),可為綜合評(píng)價(jià)保護(hù)區(qū)建設(shè)成效提供數(shù)據(jù)支撐。

    圖1 研究區(qū)示意Fig.1 Spatial distribution of the study area

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

    Landsat 系列數(shù)據(jù)具有時(shí)間跨度長(zhǎng)、重訪周期短、覆蓋范圍廣等優(yōu)點(diǎn),被廣泛用于土地利用時(shí)空格局演變研究中?;贕oogle Earth Engine(GEE)云平臺(tái)①https://earthengine.google.com/獲取金花茶保護(hù)區(qū)及其鄰近區(qū)1986、1994、1999、2004、2009、2014、2019 年7 期Landsat5 TM/Landsat8 OLI 地表反射率產(chǎn)品,影像行列號(hào)為125/45,空間分辨率為30 m。選取當(dāng)年云量最小影像,消除云雨等因素的影響。該數(shù)據(jù)產(chǎn)品已進(jìn)行大氣校正與正射校正處理。

    1.3 研究方法

    因核心區(qū)人為干擾較少,土地利用狀況穩(wěn)定,主要對(duì)比非核心區(qū)和鄰近區(qū)的土地利用變化。采用決策樹分類方法進(jìn)行土地利用遙感解譯,并基于分類結(jié)果對(duì)比非核心區(qū)和鄰近區(qū)的土地利用演變規(guī)律。為定量評(píng)價(jià)土地利用時(shí)空演變規(guī)律,選取單一土地動(dòng)態(tài)度、綜合土地利用動(dòng)態(tài)度和土地利用程度綜合指數(shù)表征土地利用變化動(dòng)態(tài)度,利用數(shù)量轉(zhuǎn)移軌跡與質(zhì)心遷移軌跡定量展示土地利用時(shí)空演變特征。

    1.3.1 土地利用類型解譯 基于金花茶保護(hù)區(qū)土地利用實(shí)際情況,并結(jié)合Landsat5 TM/Landsat8 OLI遙感數(shù)據(jù)的特點(diǎn),參考中國科學(xué)院土地利用/土地覆蓋分類系統(tǒng),利用ENVI 5.3 軟件進(jìn)行決策樹分類,將研究區(qū)土地利用劃分為林地、草地、耕地、水域、建設(shè)用地和未利用地6類(未利用地包括裸地、裸巖等),并依據(jù)野外調(diào)查數(shù)據(jù)和Google Earth高分辨率遙感影像進(jìn)行精度驗(yàn)證。

    1.3.2 土地利用類型動(dòng)態(tài)演變 不同土地利用類型的變化速度可以通過土地利用動(dòng)態(tài)度定量描述(鄔亞娟等,2020),單一土地動(dòng)態(tài)度K可用于分析某一用地類型在一定時(shí)間范圍內(nèi)的數(shù)量變化情況(張浚茂等,2019),綜合動(dòng)態(tài)度LC可以表征研究區(qū)不同時(shí)間節(jié)點(diǎn)全部用地類型的整體變化速率(張浚茂等,2019),計(jì)算公式為:

    式中:Uai、Ubi分別表示研究初期和研究末期某單一土地利用類型i的面積;T為研究時(shí)長(zhǎng);LUi為研究期初期第i類土地利用類型的面積;ΔLUi-j為研究期內(nèi)第i類土地利用類型轉(zhuǎn)為其余土地利用類型面積的絕對(duì)值(張浚茂等,2019)。

    土地利用程度可以反映人類活動(dòng)對(duì)土地利用變化的影響程度(樊嘉琦等,2018),計(jì)算公式為:

    式中:La為研究區(qū)土地利用程度綜合指數(shù),數(shù)值越大,表明該區(qū)域人類活動(dòng)越強(qiáng)烈;Ri為地類i的土地利用程度的分級(jí)指數(shù),參考自文獻(xiàn)(鄔亞娟等,2020)(表1);Ci為地類i占研究區(qū)總面積的百分比;n為土地利用類型的數(shù)量。

    表1 土地利用程度分級(jí)賦值Table 1 Grades of land use degree classification

    1.3.3 土地利用數(shù)量變化軌跡 為認(rèn)識(shí)非核心區(qū)和鄰近區(qū)不同時(shí)段內(nèi)各土地利用類型之間轉(zhuǎn)移的方向和數(shù)量,計(jì)算不同時(shí)段的土地利用變化轉(zhuǎn)移矩陣(張佰發(fā)等,2020)。轉(zhuǎn)移矩陣的公式為:

    式中:Yij表示某時(shí)期內(nèi)第i種土地利用類型轉(zhuǎn)化為第j種土地類型的面積百分比;n為土地利用類型的數(shù)量。

    通過R語言中的NetworkD3 包,采用?;鶊D(Sankey diagram)將土地利用轉(zhuǎn)移矩陣進(jìn)行可視化,構(gòu)建土地利用數(shù)量變化軌跡模型。桑基圖是一種顯示物質(zhì)或者能量流動(dòng)的可視化方法(姜婷婷等,2016),在材料分析、金融、生態(tài)學(xué)等多個(gè)領(lǐng)域被廣泛運(yùn)用。

    1.3.4 土地利用遷移軌跡 為有效反映土地利用類型的時(shí)空遷移過程,采用標(biāo)準(zhǔn)差橢圓法(The Standard Deviational Ellipse,SDE)計(jì)算各地類的質(zhì)心坐標(biāo)。其質(zhì)心坐標(biāo)計(jì)算公式為(劉濤等,2013):

    式中:SDE(x,y)表示某一地類的算數(shù)平均中心坐標(biāo),即重心;(xi,yi)表示某地類第i個(gè)圖斑的空間區(qū)位坐標(biāo);n為土地利用類型的數(shù)量。

    利用ArcGIS 10.5 疊加各時(shí)段不同地類的重心坐標(biāo),得到鄰近區(qū)和非核心區(qū)土地利用重心遷移軌跡,運(yùn)用歐式距離方法計(jì)算得到不同土地利用類型重心遷移距離。

    2 結(jié)果與分析

    采用決策樹分類法進(jìn)行遙感解譯得到金花茶保護(hù)區(qū)及鄰近區(qū)1986、1994、1999、2004、2009、2014和2019年7期土地利用數(shù)據(jù)(圖2),根據(jù)野外調(diào)查數(shù)據(jù)和Google Earth 高分辨率遙感影像對(duì)土地利用分類結(jié)果進(jìn)行精度驗(yàn)證,7 期遙感圖像的分類結(jié)果總體精度均>85%,Kappa 系數(shù)均>0.81,分類精度能滿足研究需要。

    圖2 1986—2019年研究區(qū)土地利用Fig.2 Land use pattern in the study area during 1986-2019

    2.1 土地利用總體結(jié)構(gòu)變化

    由7個(gè)時(shí)期土地利用時(shí)空分布情況可知,保護(hù)區(qū)(核心區(qū)和非核心區(qū))土地利用類型主要為林地;未利用地和草地呈片狀分布于中部、中偏西部海拔較高區(qū)域,呈現(xiàn)一定的集聚性;耕地呈點(diǎn)帶狀環(huán)保護(hù)區(qū)邊緣較平緩地勢(shì)分布(見圖2)。鄰近區(qū)土地利用類型主要為林地、耕地;林地成片狀集中分布,東部、西部分布較多,南部分布較為破碎;南部有一條明顯的耕地帶,呈西南—東北走向,建設(shè)用地成片狀分布,與耕地的分布具有一定的空間協(xié)同性。整體而言,金花茶保護(hù)區(qū)及其鄰近區(qū)土地利用類型主要為林地、耕地,二者面積之和占總面積的93%以上(圖3),呈現(xiàn)集聚性分布;林地面積分布最廣,占比為76.73%~79.46%,主要分布在研究區(qū)中部及西南部;其次是耕地,面積比例為15.60%~18.36%,呈連片狀分布,南部分布較多;水域、未利用地、建設(shè)用地、草地占總面積比例較少,面積比例均小3%,其中建設(shè)用地面積分布最少。近35 年間1986-2019 年土地利用類型發(fā)生較大變化,建設(shè)用地面積占比增長(zhǎng)幅度最大,且呈現(xiàn)逐年上升趨勢(shì),面積占比從1986年的0.19%增長(zhǎng)至2019 年的0.93%,2004—2009 年增長(zhǎng)了0.29%,較為迅速;林地呈現(xiàn)波動(dòng)性上升趨勢(shì),從1986 年的78.09%上升到了2019 年的78.82%;草地面積呈現(xiàn)略微增加趨勢(shì),面積占比從1986年的1.30%增長(zhǎng)至2019年的1.68%;未利用地面積縮減幅度最大,面積占比從1986 年的2.02%減少至2019 年的1.07%,2004-2009 年面積占比減少了1.87%;耕地和水域呈現(xiàn)先增加后減少趨勢(shì),總體上面積占比分別減少了0.76%和0.13%。

    圖3 金花茶保護(hù)區(qū)及鄰近區(qū)土地利用面積變化趨勢(shì)Fig.3 Land use changes at various stages in Golden Camellia Reserve and its adjacent areas

    2.2 土地利用動(dòng)態(tài)變化特征

    根據(jù)非核心區(qū)和鄰近區(qū)單一土地利用類型動(dòng)態(tài)度可以看出(表2),非核心區(qū)建設(shè)用地單一土地利用類型動(dòng)態(tài)度居于第一,其中,建設(shè)用地1986-2019年平均動(dòng)態(tài)度為161.60%,總體呈現(xiàn)顯著擴(kuò)張趨勢(shì),1994—1999年動(dòng)態(tài)度最高為56%,2009年之后動(dòng)態(tài)度為負(fù)值且呈現(xiàn)下降趨勢(shì);草地的動(dòng)態(tài)度相對(duì)于其余用地類型一直處于較高水平,年平均動(dòng)態(tài)度為1.84%,總體位居第二且呈現(xiàn)擴(kuò)張趨勢(shì);林地呈現(xiàn)微增趨勢(shì);未利用地1986—2019 年平均動(dòng)態(tài)度為-1.58%,呈現(xiàn)波動(dòng)減少趨勢(shì);水域和耕地的動(dòng)態(tài)度波動(dòng)較平穩(wěn),呈微小降低趨勢(shì)。鄰近區(qū)建設(shè)用地的動(dòng)態(tài)度一直處于較大值(見表2),1986-2019年平均動(dòng)態(tài)度為11.42%,呈現(xiàn)擴(kuò)張趨勢(shì);其余5 種土地利用類型的動(dòng)態(tài)度波動(dòng)較為平緩,林地呈現(xiàn)微小增加趨勢(shì),草地、未利用地、水域和耕地呈微小縮小趨勢(shì)。綜上,非核心區(qū)和鄰近區(qū)建設(shè)用地的動(dòng)態(tài)度最高,但非核心區(qū)的建設(shè)用地動(dòng)態(tài)度更大,2個(gè)地區(qū)近35年的單一土地利用類型動(dòng)態(tài)變化除草地外其他類型趨于一致,非核心區(qū)草地為正值且排名第二,鄰近區(qū)草地則為負(fù)值。

    表2 鄰近區(qū)和非核心區(qū)單一土地利用動(dòng)態(tài)度Table 2 The rate of land use change at various stages in adjacent areas and non-core areas%

    1986—2019年,鄰近區(qū)綜合土地利用動(dòng)態(tài)度曲線波動(dòng)較大(圖4-a),呈現(xiàn)先急劇減少再上升到最高值,之后波動(dòng)下降的趨勢(shì)。1994—1999年處于較低值,1999—2004 年上升到最高值0.63%,表明該時(shí)期人類活動(dòng)對(duì)鄰近區(qū)土地利用變化造成較大影響;2004年以后呈波動(dòng)下降趨勢(shì),土地利用變化趨于平緩。非核心區(qū)1986—2019年綜合土地利用動(dòng)態(tài)度呈現(xiàn)先增加后遞減的趨勢(shì),2004—2009年處于最高值0.55%,該時(shí)期土地利用變化劇烈。整體而言,鄰近區(qū)綜合土地利用動(dòng)態(tài)度比非核心區(qū)變化更劇烈,表明保護(hù)區(qū)內(nèi)部非核心區(qū)土地利用格局更為穩(wěn)定,受人為活動(dòng)影響較小,鄰近區(qū)人地矛盾更為劇烈,2 個(gè)區(qū)域2009—2019 年綜合動(dòng)態(tài)度都趨于下降,表示土地利用變化都趨于平緩。

    圖4 鄰近區(qū)和非核心區(qū)1986—2019年綜合土地利用動(dòng)態(tài)度(a)和土地利用程度綜合指數(shù)(b)Fig.4 Comprehensive land use dynamic degree(a)and comprehensive index of land use degree(b)in the adjacent and non-core areas from 1986 to 2019

    鄰近區(qū)和非核心區(qū)土地利用程度綜合指數(shù)變化趨勢(shì)相似,均呈現(xiàn)先上升再下降,之后較為平緩波動(dòng)趨勢(shì)(圖4-b)。鄰近區(qū)在1999 年達(dá)到最大值229.58后略微下降,2004年以后平緩上升,這與該時(shí)期土地利用方式轉(zhuǎn)換相一致,表示土地利用程度高的建設(shè)用地在擴(kuò)張,未利用地等土地利用程度較低的地類面積占比減少,土地利用強(qiáng)度顯著提高。非核心區(qū)土地利用程度綜合指數(shù)在1999年達(dá)到最大值,隨后呈波動(dòng)下降趨勢(shì),該區(qū)域土地利用程度綜合指數(shù)始終比鄰近區(qū)低,表明非核心的土地利用受人為干擾較少,建設(shè)用地、耕地等高土地利用強(qiáng)度所占比例得到控制并逐漸減少。

    2.3 土地利用數(shù)量轉(zhuǎn)移軌跡

    從金花茶保護(hù)區(qū)土地利用數(shù)量轉(zhuǎn)換軌跡看(圖5) (空間上的土地利用類型轉(zhuǎn)移情況見文末附圖1),鄰近區(qū)1986—2019 年土地利用類型轉(zhuǎn)換頻繁(圖5-a)。其中,林地轉(zhuǎn)換量最為顯著,主要流向耕地、草地,1986—1994 年流向耕地面積高達(dá)1 423.33 hm2,1994—1999 年變化相對(duì)較少,其他時(shí)期轉(zhuǎn)移面積波動(dòng)減緩;耕地以流出為主,主要流向?yàn)榱值亍⒔ㄔO(shè)用地、草地,1986—1994年流向林地面積最大,達(dá)到1 149.14 hm2;建設(shè)用地流入行為明顯,主要由耕地、未利用地、林地流入,由耕地流入面積逐年增加。就土地利用類型轉(zhuǎn)換的軌跡線而言,1986—1994 年不同地類之間轉(zhuǎn)化最為活躍,轉(zhuǎn)換軌跡多樣化;1994—1999年不同地類之間轉(zhuǎn)換相對(duì)較少;1999—2019年4個(gè)時(shí)段不同土地利用類型間的轉(zhuǎn)變幅度較為一致。整體上,1986—2019年鄰近區(qū)地類轉(zhuǎn)換軌跡線比較稠密,各地類相互轉(zhuǎn)換比較多樣,不同類型間變化較活躍。

    附圖1 金花茶國家級(jí)自然保護(hù)區(qū)鄰近區(qū)、非核心區(qū)1986—2019年土地利用變化情況Accompanying Fig.1 Land use change in adjacent and non-core areas of the Fangcheng Camellia National Nature Reserve from 1986 to 2019

    非核心區(qū)1986—2019年土地利用類型之間轉(zhuǎn)換較少(圖5-b)。林地流入行為最為明顯,主要由耕地、草地、未利用地流入,耕地轉(zhuǎn)入面積最大,1986—1994 年流入面積達(dá)最大值212.75 hm2,其余年間轉(zhuǎn)換面積相對(duì)較少;草地以流入為主,主要由林地、未利用地、耕地轉(zhuǎn)入,2004—2009年由未利用地轉(zhuǎn)入52.73 hm2;未利用地流出明顯,主要轉(zhuǎn)化為草地、林地,2004—2009 年轉(zhuǎn)換為草地104.22 hm2;耕地、水域和建設(shè)用地與其他土地利用類型之間的轉(zhuǎn)移變化相對(duì)較少。與鄰近區(qū)對(duì)比,非核心區(qū)軌跡線較為稀疏,說明該區(qū)域土地利用類型變化較少,且轉(zhuǎn)變類型更簡(jiǎn)單。綜上,保護(hù)區(qū)的建設(shè)與管理顯著效果,當(dāng)?shù)厣鷳B(tài)環(huán)境得到有效保護(hù)。

    圖5 1986—2019年土地利用數(shù)量變化軌跡Fig.5 Land use change trajectory from 1986 to 2019

    2.4 土地利用時(shí)空遷移分析

    土地利用重心遷移軌跡分析結(jié)果表明(圖6),鄰近區(qū)林地重心呈現(xiàn)三角形狀遷移特征,總體向東北方向遷移,1999—2004年遷移幅度最大;耕地重心呈現(xiàn)“V”字型遷移特征,總體向西南方向遷移,1999—2004年遷移幅度最大;建設(shè)用地重心向西北方向遷移,1986—1994年遷移幅度最大;草地重心向東北方向遷移,2004—2009年遷移幅度最大;未利用地重心向東北方向遷移,1986—1994年遷移幅度最大;水域重心遷移幅度較小,向東北方向遷移。耕地與建設(shè)耕地總體遷移方向都為向西遷移。

    圖6 鄰近區(qū)、非核心區(qū)土地利用重心遷移軌跡Fig.6 Migration trajectory of land use gravity center in adjacent areas and non-core areas

    非核心區(qū)6種土地利用類型重心整體上都是向東北方向遷移,建設(shè)用地重心遷移幅度最大,在1986—1994年遷移距離最遠(yuǎn);林地在1999—2004年重心遷移幅度最大;耕地重心呈現(xiàn)“W”字型遷移特征,2004—2009年遷移幅度最大;草地在1986—1994年重心遷移幅度最大;未利用地在2014—2019年重心遷移幅度最大;水域重心呈現(xiàn)“Z”字型遷移特征。

    3 討論

    綜上,保護(hù)區(qū)周邊鄰近區(qū)土地利用變化劇烈,土地利用類型轉(zhuǎn)換頻繁,非核心區(qū)土地利用變化較平緩,土地利用類型相互轉(zhuǎn)換較少,表明保護(hù)區(qū)建立后管護(hù)政策的貫徹實(shí)施,管護(hù)力度不斷加強(qiáng),減少了人為活動(dòng)的干擾,保護(hù)區(qū)的生態(tài)環(huán)境質(zhì)量得到改善。1986—2019年非核心區(qū)和鄰近區(qū)建設(shè)用地動(dòng)態(tài)度變化最大,都為擴(kuò)張趨勢(shì),主要由耕地流入,保護(hù)區(qū)為加強(qiáng)開發(fā)管理對(duì)自然保護(hù)區(qū)內(nèi)的公路進(jìn)行擴(kuò)建,這是非核心區(qū)用地動(dòng)態(tài)度變化的主要原因。鄰近區(qū)土地利用變化主要受那梭農(nóng)場(chǎng)和沖皇溝生態(tài)旅游區(qū)的發(fā)展建設(shè)影響,西南方向原有居民地向西北方向不斷擴(kuò)張,說明西部人類活動(dòng)逐漸增多,社會(huì)經(jīng)濟(jì)活動(dòng)活躍,需對(duì)保護(hù)區(qū)西部加強(qiáng)管理;退耕還林、沼氣池建設(shè)工程等政策的實(shí)施,使非核心區(qū)和鄰近區(qū)林地面積都有所增加,非核心區(qū)的林地因得到保護(hù)增幅更為明顯,2 個(gè)區(qū)域林地都往東北方向遷移,與森林資源消耗的減少、東北方向耕地退耕有一定關(guān)聯(lián);兩個(gè)區(qū)域的耕地均呈減少趨勢(shì),與建設(shè)用地、林地的動(dòng)態(tài)演變具有一定協(xié)同性;非核心區(qū)未利用地主要分布在海拔較高山頂區(qū)域;隨著保護(hù)區(qū)生態(tài)環(huán)境的改善,西南方向較低海拔的未利用地轉(zhuǎn)換為林地和草地,鄰近區(qū)的未利用地主要轉(zhuǎn)換為林地與耕地。

    從時(shí)間梯度看,鄰近區(qū)綜合動(dòng)態(tài)度在1994—1999 年處于較低值,土地類型轉(zhuǎn)換較少,這與90年代以后外出務(wù)工人員大量增加,土地利用受人為干擾下降有一定關(guān)系,當(dāng)?shù)鼐用裢ㄟ^務(wù)工收入水平顯著提高,擴(kuò)大住宅面積,土地利用程度高的建設(shè)用地在擴(kuò)張,使得土地利用程度上升;1999—2004年,那灣生態(tài)園的建設(shè)、那梭農(nóng)場(chǎng)特色產(chǎn)業(yè)的開發(fā),建設(shè)用地在西部顯著擴(kuò)張,人類活動(dòng)劇烈,土地利用類型之間轉(zhuǎn)換劇烈,以致土地綜合動(dòng)態(tài)度較高,2002年退耕還林政策的實(shí)施,使得耕地面積減少,土地利用程度下降;2003年國家重點(diǎn)公益林補(bǔ)助試行之后,當(dāng)?shù)鼐用窠ㄔO(shè)沼氣池的積極性提高,減少了森林資源消耗,退耕還林政策的繼續(xù)實(shí)施,有效控制了耕地的面積,但建設(shè)用地不斷擴(kuò)張,土地利用動(dòng)態(tài)度趨于波動(dòng)增長(zhǎng),土地利用程度平緩提升。非核心區(qū)受道路建設(shè)、退耕還林等政策的影響,土地利用變化較強(qiáng)烈,土地利用類型之間轉(zhuǎn)換較頻繁,1986—2009年土地利用動(dòng)態(tài)度一直處于增長(zhǎng)趨勢(shì),2009年之后呈現(xiàn)下降趨勢(shì),表明人類活動(dòng)逐漸較少,土地利用格局趨于穩(wěn)定。

    整體而言,保護(hù)區(qū)的設(shè)立,退耕還林還草等工程的實(shí)施,管護(hù)工作力度的不斷加強(qiáng),使得保護(hù)區(qū)內(nèi)部土地利用受人為干擾下降,土地利用格局更為穩(wěn)定,建設(shè)用地、耕地等高土地利用強(qiáng)度所占比例得到控制并逐漸降低。同時(shí),保護(hù)區(qū)的建立也促進(jìn)了周邊生態(tài)旅游、礦產(chǎn)資源、水資源的開發(fā)以及經(jīng)濟(jì)作物產(chǎn)業(yè)的發(fā)展,促使周邊人類活動(dòng)增加,土地利用變化劇烈,土地利用動(dòng)態(tài)度一直處于較高水平??傮w上,保護(hù)區(qū)的管理取得較為顯著的成效,內(nèi)部生態(tài)環(huán)境得到改善,但內(nèi)部建設(shè)用地的緩慢增加與資源的不合理生產(chǎn)開發(fā),仍然對(duì)保護(hù)區(qū)的生境質(zhì)量產(chǎn)生負(fù)面影響。保護(hù)區(qū)在未來的管理建設(shè)中,應(yīng)合理控制區(qū)內(nèi)人類活動(dòng),強(qiáng)化土地集約化發(fā)展,有效實(shí)施生態(tài)移民搬遷、異地安置等政策。

    金花茶保護(hù)區(qū)位于海拔較高區(qū)域,其地形特點(diǎn)在一定程度上也驅(qū)動(dòng)著土地利用類型及其動(dòng)態(tài)演變(崔步禮等,2011)。為更精準(zhǔn)地評(píng)估保護(hù)區(qū)的建設(shè)成效,提取1986—2019年非核心區(qū)坡度<6°的土地利用數(shù)據(jù),并進(jìn)行土地利用動(dòng)態(tài)度的計(jì)算,1986-2019 年7 個(gè)年份土地利用強(qiáng)度分別為224.41、221.54、223.96、223.02、224.15、222.01、221.20,均小于鄰近區(qū)的226.99、229.58、229.59、226.51、228.72、227.68、228.86,證明保護(hù)區(qū)的建設(shè)對(duì)生態(tài)環(huán)境的提高有促進(jìn)作用。

    4 結(jié)論

    利用1986—2019年遙感影像對(duì)金花茶保護(hù)區(qū)及其鄰近區(qū)進(jìn)行土地利用遙感解譯,并基于GIS空間分析與R語言統(tǒng)計(jì)分析方法,探討了研究區(qū)近35年來土地利用動(dòng)態(tài)演變規(guī)律,得到的主要結(jié)論有:

    1)研究區(qū)土地利用類型整體分布具有一定的集聚性,保護(hù)區(qū)土地利用類型以林地為主,外部鄰近區(qū)以林地和耕地為主,近35年間土地利用類型發(fā)生較大變化,建設(shè)用地面積占比增長(zhǎng)幅度最大,未利用地面積縮減幅度最大,其余地類呈現(xiàn)波動(dòng)增減趨勢(shì)。

    2)自保護(hù)區(qū)建立以來,鄰近區(qū)土地利用變化強(qiáng)烈,非核心區(qū)土地利用格局更為穩(wěn)定;非核心區(qū)和鄰近區(qū)建設(shè)用地的動(dòng)態(tài)度一直處于較大值,其余5 種土地利用類型的動(dòng)態(tài)度波動(dòng)較為平緩;鄰近區(qū)綜合土地利用動(dòng)態(tài)度比非核心區(qū)變化更劇烈;鄰近區(qū)和非核心區(qū)土地利用程度綜合指數(shù)變化趨勢(shì)相似,均呈現(xiàn)先上升再下降,之后較為平緩波動(dòng)趨勢(shì)。

    3)鄰近區(qū)近35年間土地利用類型之間轉(zhuǎn)換頻繁,相對(duì)而言,非核心區(qū)土地利用類型相互轉(zhuǎn)換較少,2 個(gè)區(qū)域均是耕地與林地之間的轉(zhuǎn)移交換最為劇烈。鄰近區(qū)轉(zhuǎn)換軌跡線比較稠密,各地類相互轉(zhuǎn)換呈多樣化,非核心區(qū)軌跡線較為稀疏。

    4)土地利用重心遷移具有一定的方向性,非核心區(qū)6 種土地利用類型總體上都是向東北遷移,鄰近區(qū)建設(shè)用地與耕地重心向西遷移,其余地類重心向東北遷移,表明建設(shè)用地、耕地在已有基礎(chǔ)上不斷擴(kuò)張。

    5)保護(hù)區(qū)的管理取得較為顯著的成效,保護(hù)區(qū)內(nèi)部土地利用受人為干擾下降,土地利用格局更為穩(wěn)定,生態(tài)環(huán)境得到改善,但內(nèi)部建設(shè)用地的緩慢增加與資源的不合理生產(chǎn)開發(fā),仍然對(duì)保護(hù)區(qū)的生境質(zhì)量產(chǎn)生負(fù)面影響。

    本研究定量對(duì)比分析了金花茶保護(hù)區(qū)非核心區(qū)和外部鄰近區(qū)在保護(hù)區(qū)建設(shè)與發(fā)展歷程中土地利用動(dòng)態(tài)演變,探究保護(hù)區(qū)的建立對(duì)區(qū)域土地利用動(dòng)態(tài)的影響,可為保護(hù)區(qū)的建設(shè)及監(jiān)管、資源開發(fā)及生態(tài)安全保護(hù)提供一定的數(shù)據(jù)支撐。

    猜你喜歡
    金花核心區(qū)保護(hù)區(qū)
    故鄉(xiāng)的金花茶
    黃河之聲(2021年18期)2021-12-14 04:10:40
    河南董寨國家級(jí)自然保護(hù)區(qū)鳥類名錄與分布
    品牌研究(2021年7期)2021-11-28 03:40:50
    黃河三角洲保護(hù)區(qū)自然資源的開發(fā)與保護(hù)
    金花
    打造永久“大糧倉” 助推糧食核心區(qū)建設(shè)
    打造永久"大糧倉"助推糧食核心區(qū)建設(shè)
    打造永久"大糧倉"助推糧食核心區(qū)建設(shè)
    五朵金花
    如何保護(hù)好興隆山保護(hù)區(qū)內(nèi)的馬麝
    涉及各類保護(hù)區(qū)的項(xiàng)目建設(shè)的探討
    河南科技(2014年7期)2014-02-27 14:11:33
    99九九在线精品视频| 麻豆乱淫一区二区| 一区二区日韩欧美中文字幕| 亚洲欧美一区二区三区黑人| 亚洲国产av新网站| 777米奇影视久久| av国产久精品久网站免费入址| 飞空精品影院首页| 美女高潮到喷水免费观看| 青草久久国产| 天堂8中文在线网| 成年人午夜在线观看视频| 国产精品蜜桃在线观看| 在线观看一区二区三区激情| 日日啪夜夜爽| 国产日韩欧美视频二区| 嫩草影视91久久| 久久久精品免费免费高清| 在线精品无人区一区二区三| 波野结衣二区三区在线| 999久久久国产精品视频| 9热在线视频观看99| 中文字幕av电影在线播放| 欧美乱码精品一区二区三区| 老司机亚洲免费影院| 麻豆精品久久久久久蜜桃| 男人舔女人的私密视频| 国产精品麻豆人妻色哟哟久久| 秋霞伦理黄片| 在线精品无人区一区二区三| 久久精品人人爽人人爽视色| 欧美激情高清一区二区三区 | 精品久久蜜臀av无| 国产一区二区 视频在线| 热99国产精品久久久久久7| 亚洲国产最新在线播放| 精品亚洲成国产av| 黑丝袜美女国产一区| 久久这里只有精品19| 国产成人a∨麻豆精品| 悠悠久久av| 欧美黑人欧美精品刺激| 操出白浆在线播放| 精品一区在线观看国产| av又黄又爽大尺度在线免费看| 熟女av电影| 欧美黄色片欧美黄色片| 久久久久精品人妻al黑| 亚洲成人av在线免费| 亚洲国产欧美一区二区综合| 哪个播放器可以免费观看大片| 国产成人免费观看mmmm| 丰满乱子伦码专区| 丁香六月欧美| 99久国产av精品国产电影| 亚洲三区欧美一区| av在线app专区| 久久精品aⅴ一区二区三区四区| 老汉色av国产亚洲站长工具| 亚洲一码二码三码区别大吗| 亚洲欧洲国产日韩| 男女国产视频网站| 亚洲欧美一区二区三区久久| 亚洲精品国产区一区二| xxx大片免费视频| xxx大片免费视频| 天天躁日日躁夜夜躁夜夜| 精品一区在线观看国产| 亚洲国产精品一区三区| 中文字幕精品免费在线观看视频| 97人妻天天添夜夜摸| 中文字幕制服av| 亚洲,欧美,日韩| 午夜av观看不卡| 久久青草综合色| 欧美成人精品欧美一级黄| 精品人妻熟女毛片av久久网站| av免费观看日本| av天堂久久9| 高清欧美精品videossex| 亚洲精品在线美女| av线在线观看网站| 熟女av电影| 久久精品国产亚洲av高清一级| 人人妻人人澡人人看| 在线观看国产h片| 国产成人a∨麻豆精品| 亚洲精品第二区| 欧美精品高潮呻吟av久久| 欧美黄色片欧美黄色片| 大话2 男鬼变身卡| 日本一区二区免费在线视频| 国产精品 欧美亚洲| 国产成人a∨麻豆精品| 亚洲欧美日韩另类电影网站| 91成人精品电影| 国产精品二区激情视频| 午夜福利网站1000一区二区三区| 久久性视频一级片| 成年女人毛片免费观看观看9 | 欧美黑人精品巨大| 99精国产麻豆久久婷婷| 国产欧美日韩综合在线一区二区| xxxhd国产人妻xxx| 免费黄色在线免费观看| 欧美精品高潮呻吟av久久| 国产不卡av网站在线观看| 宅男免费午夜| 天堂中文最新版在线下载| 国产熟女午夜一区二区三区| 亚洲国产欧美在线一区| 亚洲精品av麻豆狂野| 久久人人爽人人片av| 在线观看人妻少妇| 少妇被粗大的猛进出69影院| 久久精品亚洲熟妇少妇任你| 免费看av在线观看网站| 日韩精品免费视频一区二区三区| 丝袜美腿诱惑在线| 又大又黄又爽视频免费| 亚洲欧洲国产日韩| 久久热在线av| 9热在线视频观看99| av在线app专区| 成年人免费黄色播放视频| 亚洲精品日韩在线中文字幕| e午夜精品久久久久久久| 人人澡人人妻人| 水蜜桃什么品种好| av在线app专区| 国产欧美亚洲国产| 超色免费av| 国产精品三级大全| 一二三四中文在线观看免费高清| 国产亚洲av片在线观看秒播厂| 亚洲精品aⅴ在线观看| 久久久久精品国产欧美久久久 | 99热网站在线观看| 国产亚洲av高清不卡| 男女床上黄色一级片免费看| 日韩大片免费观看网站| 久久久久国产精品人妻一区二区| 天天添夜夜摸| 国产精品 欧美亚洲| 亚洲精品国产区一区二| 久久精品国产亚洲av高清一级| 男女国产视频网站| 97在线人人人人妻| 国产精品国产三级专区第一集| 国产99久久九九免费精品| av线在线观看网站| 国产精品久久久久久人妻精品电影 | svipshipincom国产片| 国产亚洲一区二区精品| 成人亚洲欧美一区二区av| 99精品久久久久人妻精品| 精品午夜福利在线看| 日日摸夜夜添夜夜爱| 中文字幕精品免费在线观看视频| 男人舔女人的私密视频| 亚洲精品美女久久av网站| 亚洲精品视频女| 免费久久久久久久精品成人欧美视频| 亚洲欧美一区二区三区国产| 亚洲欧美色中文字幕在线| 热99久久久久精品小说推荐| 久久久国产精品麻豆| 男人操女人黄网站| √禁漫天堂资源中文www| 亚洲一区中文字幕在线| www.自偷自拍.com| 另类亚洲欧美激情| 在线天堂最新版资源| 人人妻人人澡人人看| 亚洲 欧美一区二区三区| 亚洲成人一二三区av| 欧美黑人精品巨大| 人成视频在线观看免费观看| 日日撸夜夜添| 青青草视频在线视频观看| 一区二区三区乱码不卡18| 三上悠亚av全集在线观看| 亚洲精品美女久久久久99蜜臀 | 久久久精品国产亚洲av高清涩受| 极品人妻少妇av视频| 天堂8中文在线网| 一边摸一边做爽爽视频免费| 日韩电影二区| 制服丝袜香蕉在线| 狂野欧美激情性xxxx| 美国免费a级毛片| tube8黄色片| 免费高清在线观看视频在线观看| 又黄又粗又硬又大视频| 在线精品无人区一区二区三| 日本猛色少妇xxxxx猛交久久| 黄色 视频免费看| 51午夜福利影视在线观看| 成人18禁高潮啪啪吃奶动态图| 自线自在国产av| 国产成人精品久久二区二区91 | 一区二区日韩欧美中文字幕| 亚洲免费av在线视频| 国产亚洲av高清不卡| 久久国产精品大桥未久av| 大码成人一级视频| www.自偷自拍.com| 中文欧美无线码| 婷婷色综合www| 国产成人一区二区在线| 大码成人一级视频| 亚洲伊人色综图| 久久天躁狠狠躁夜夜2o2o | 男人舔女人的私密视频| 大香蕉久久成人网| 2021少妇久久久久久久久久久| 免费在线观看完整版高清| 一区在线观看完整版| av国产精品久久久久影院| 男男h啪啪无遮挡| 韩国av在线不卡| 在线观看免费视频网站a站| 亚洲av日韩在线播放| 成年女人毛片免费观看观看9 | 免费久久久久久久精品成人欧美视频| 新久久久久国产一级毛片| 这个男人来自地球电影免费观看 | 国产精品久久久久久人妻精品电影 | 欧美在线一区亚洲| 亚洲综合色网址| 亚洲自偷自拍图片 自拍| 精品亚洲成a人片在线观看| 欧美激情高清一区二区三区 | 波野结衣二区三区在线| av一本久久久久| 韩国精品一区二区三区| 制服诱惑二区| 亚洲精品日本国产第一区| 成年人免费黄色播放视频| 午夜久久久在线观看| 国产深夜福利视频在线观看| 国产熟女午夜一区二区三区| 国产乱来视频区| 精品国产乱码久久久久久男人| 综合色丁香网| 另类亚洲欧美激情| 日韩中文字幕视频在线看片| 亚洲av在线观看美女高潮| 18在线观看网站| 一本久久精品| 在线天堂中文资源库| www.精华液| 女人高潮潮喷娇喘18禁视频| 91aial.com中文字幕在线观看| 国产老妇伦熟女老妇高清| 亚洲av日韩精品久久久久久密 | 日日摸夜夜添夜夜爱| 成人午夜精彩视频在线观看| 亚洲综合色网址| 中文字幕制服av| 免费观看人在逋| 女人精品久久久久毛片| xxxhd国产人妻xxx| 久久久久久人妻| 国产av码专区亚洲av| xxx大片免费视频| 国产精品99久久99久久久不卡 | 伦理电影免费视频| 亚洲av成人不卡在线观看播放网 | 亚洲av综合色区一区| 捣出白浆h1v1| 国产精品秋霞免费鲁丝片| 一本久久精品| 少妇精品久久久久久久| 色网站视频免费| 色综合欧美亚洲国产小说| 国产 精品1| 日韩视频在线欧美| 男的添女的下面高潮视频| 国产精品 欧美亚洲| 青春草国产在线视频| 国产探花极品一区二区| 99热网站在线观看| 黄网站色视频无遮挡免费观看| 久久精品国产亚洲av高清一级| av卡一久久| 欧美人与善性xxx| 一边亲一边摸免费视频| 午夜影院在线不卡| 国产成人欧美在线观看 | 性高湖久久久久久久久免费观看| 色视频在线一区二区三区| 成人午夜精彩视频在线观看| 免费久久久久久久精品成人欧美视频| 国产色婷婷99| xxx大片免费视频| 高清在线视频一区二区三区| 亚洲国产av影院在线观看| 婷婷成人精品国产| 91国产中文字幕| 久久久久精品性色| 精品福利永久在线观看| 在线 av 中文字幕| 国产免费福利视频在线观看| 亚洲自偷自拍图片 自拍| 久久久国产精品麻豆| 我的亚洲天堂| 尾随美女入室| 9色porny在线观看| 欧美黑人欧美精品刺激| 久久热在线av| 人人妻人人爽人人添夜夜欢视频| 久久午夜综合久久蜜桃| 欧美老熟妇乱子伦牲交| 成人18禁高潮啪啪吃奶动态图| 中文字幕亚洲精品专区| 亚洲av日韩精品久久久久久密 | 午夜福利网站1000一区二区三区| 久久影院123| 黑人巨大精品欧美一区二区蜜桃| 在线免费观看不下载黄p国产| 男女床上黄色一级片免费看| 精品国产一区二区久久| 国产精品熟女久久久久浪| 赤兔流量卡办理| 操出白浆在线播放| 五月天丁香电影| 深夜精品福利| 桃花免费在线播放| 一二三四在线观看免费中文在| 精品人妻熟女毛片av久久网站| 亚洲av男天堂| 国产亚洲av高清不卡| 性高湖久久久久久久久免费观看| 精品午夜福利在线看| 久久精品久久久久久久性| 亚洲人成77777在线视频| 2018国产大陆天天弄谢| 婷婷成人精品国产| 精品免费久久久久久久清纯 | 只有这里有精品99| 亚洲色图综合在线观看| 欧美 日韩 精品 国产| 精品酒店卫生间| 色婷婷av一区二区三区视频| 国产精品嫩草影院av在线观看| 国产精品香港三级国产av潘金莲 | 欧美日韩综合久久久久久| 桃花免费在线播放| 丝袜脚勾引网站| 午夜福利,免费看| 欧美日韩一级在线毛片| 亚洲婷婷狠狠爱综合网| 国产成人系列免费观看| 母亲3免费完整高清在线观看| 午夜激情久久久久久久| 韩国av在线不卡| www.av在线官网国产| 菩萨蛮人人尽说江南好唐韦庄| 热99国产精品久久久久久7| 两个人免费观看高清视频| 中文字幕人妻熟女乱码| av电影中文网址| 考比视频在线观看| 欧美精品一区二区免费开放| 最近中文字幕高清免费大全6| 99久久人妻综合| 久久久国产精品麻豆| 各种免费的搞黄视频| 狂野欧美激情性bbbbbb| 午夜福利,免费看| 99精品久久久久人妻精品| 婷婷色综合www| 热99久久久久精品小说推荐| 婷婷色综合大香蕉| 欧美 亚洲 国产 日韩一| 1024香蕉在线观看| 亚洲精品国产色婷婷电影| 日韩电影二区| av网站在线播放免费| 久久久国产欧美日韩av| 亚洲免费av在线视频| 黑人巨大精品欧美一区二区蜜桃| 日本一区二区免费在线视频| 午夜精品国产一区二区电影| 操出白浆在线播放| 狠狠精品人妻久久久久久综合| 国产精品无大码| 亚洲熟女毛片儿| 亚洲激情五月婷婷啪啪| 欧美日韩精品网址| 19禁男女啪啪无遮挡网站| 国产成人91sexporn| 一区二区三区乱码不卡18| 国产成人系列免费观看| 在线观看免费日韩欧美大片| 亚洲国产欧美网| 国产精品香港三级国产av潘金莲 | 电影成人av| 老司机影院成人| 肉色欧美久久久久久久蜜桃| 午夜日韩欧美国产| 国产高清不卡午夜福利| 欧美乱码精品一区二区三区| 国产精品av久久久久免费| 午夜日韩欧美国产| 黑丝袜美女国产一区| 午夜日韩欧美国产| kizo精华| 中文欧美无线码| 免费看不卡的av| 亚洲av日韩精品久久久久久密 | 国产亚洲午夜精品一区二区久久| 色婷婷久久久亚洲欧美| 国产日韩欧美在线精品| 国产女主播在线喷水免费视频网站| 综合色丁香网| 久久精品熟女亚洲av麻豆精品| 日本色播在线视频| 精品少妇内射三级| 丁香六月天网| 国产精品成人在线| 午夜日韩欧美国产| 日本欧美视频一区| 国产成人精品无人区| 午夜免费鲁丝| 日韩电影二区| 少妇人妻久久综合中文| 丝袜在线中文字幕| 亚洲精品自拍成人| 日本欧美国产在线视频| 丝袜美腿诱惑在线| 欧美乱码精品一区二区三区| 黄色一级大片看看| 两个人看的免费小视频| 赤兔流量卡办理| 国产xxxxx性猛交| 只有这里有精品99| 中文乱码字字幕精品一区二区三区| 国产在线免费精品| 欧美精品av麻豆av| 十八禁高潮呻吟视频| 精品国产乱码久久久久久男人| 亚洲 欧美一区二区三区| 亚洲国产精品999| av一本久久久久| 精品一区二区三卡| 欧美日韩福利视频一区二区| 一级黄片播放器| 精品国产一区二区三区久久久樱花| 久久99热这里只频精品6学生| 国产精品久久久久成人av| 极品少妇高潮喷水抽搐| 欧美日韩视频精品一区| 国产精品免费视频内射| 亚洲第一区二区三区不卡| 中文字幕最新亚洲高清| 国产日韩欧美视频二区| 母亲3免费完整高清在线观看| 婷婷色综合www| 久久久国产精品麻豆| 校园人妻丝袜中文字幕| 亚洲成人国产一区在线观看 | 99热网站在线观看| 18禁观看日本| 午夜av观看不卡| 亚洲三区欧美一区| 色精品久久人妻99蜜桃| 国产片特级美女逼逼视频| 欧美日韩一区二区视频在线观看视频在线| 卡戴珊不雅视频在线播放| 亚洲久久久国产精品| 精品第一国产精品| 国产精品久久久人人做人人爽| 又粗又硬又长又爽又黄的视频| 午夜福利视频在线观看免费| 一区二区av电影网| 丝袜在线中文字幕| 国产 一区精品| 欧美 亚洲 国产 日韩一| 少妇精品久久久久久久| 80岁老熟妇乱子伦牲交| 欧美亚洲 丝袜 人妻 在线| 午夜91福利影院| 国产精品.久久久| 久久国产亚洲av麻豆专区| 在线天堂中文资源库| 日本爱情动作片www.在线观看| 午夜福利一区二区在线看| 国产成人啪精品午夜网站| 一区福利在线观看| 成年女人毛片免费观看观看9 | 国产欧美日韩综合在线一区二区| 国产爽快片一区二区三区| 日本wwww免费看| 97在线人人人人妻| 性少妇av在线| 国产在线一区二区三区精| 七月丁香在线播放| 欧美精品av麻豆av| 欧美97在线视频| 国产成人免费观看mmmm| 男女午夜视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 日韩一区二区视频免费看| 国产激情久久老熟女| 国产不卡av网站在线观看| 久久久久人妻精品一区果冻| 精品人妻在线不人妻| 亚洲欧洲日产国产| av视频免费观看在线观看| 丰满迷人的少妇在线观看| 十八禁网站网址无遮挡| 日日摸夜夜添夜夜爱| 欧美日韩综合久久久久久| 欧美乱码精品一区二区三区| 欧美日韩综合久久久久久| 日韩欧美精品免费久久| 精品国产一区二区三区久久久樱花| 午夜91福利影院| 欧美激情高清一区二区三区 | 久久av网站| 青青草视频在线视频观看| 99re6热这里在线精品视频| 91老司机精品| 国产在线免费精品| 男男h啪啪无遮挡| 男女床上黄色一级片免费看| 日韩av在线免费看完整版不卡| 亚洲欧美清纯卡通| 男的添女的下面高潮视频| 黄网站色视频无遮挡免费观看| 伦理电影大哥的女人| 中文字幕高清在线视频| 男女之事视频高清在线观看 | 成人国产av品久久久| 成人黄色视频免费在线看| 国产无遮挡羞羞视频在线观看| 午夜av观看不卡| 在线亚洲精品国产二区图片欧美| 操出白浆在线播放| 国产深夜福利视频在线观看| 日日撸夜夜添| 精品国产一区二区三区四区第35| 欧美国产精品一级二级三级| 男女午夜视频在线观看| 亚洲成av片中文字幕在线观看| 两个人看的免费小视频| 国产欧美亚洲国产| 精品少妇内射三级| 精品国产露脸久久av麻豆| 不卡av一区二区三区| 久久久久久久久久久久大奶| 亚洲欧美激情在线| 国产熟女午夜一区二区三区| 亚洲专区中文字幕在线 | 久久人人爽人人片av| 国产精品成人在线| 黄色 视频免费看| 精品国产乱码久久久久久小说| 日韩制服丝袜自拍偷拍| 国产一区二区在线观看av| 无遮挡黄片免费观看| 国产精品人妻久久久影院| 日本欧美视频一区| 久久久精品94久久精品| 欧美人与善性xxx| 久久韩国三级中文字幕| 一级毛片黄色毛片免费观看视频| av女优亚洲男人天堂| 又大又爽又粗| 18禁动态无遮挡网站| 午夜老司机福利片| 新久久久久国产一级毛片| 中文字幕色久视频| 国产精品国产av在线观看| 久久综合国产亚洲精品| 精品人妻熟女毛片av久久网站| 毛片一级片免费看久久久久| 国产伦人伦偷精品视频| 亚洲欧美中文字幕日韩二区| 午夜日本视频在线| 中文字幕亚洲精品专区| 最新在线观看一区二区三区 | 亚洲专区中文字幕在线 | 在线观看一区二区三区激情| √禁漫天堂资源中文www| 欧美黑人欧美精品刺激| 热re99久久国产66热| 亚洲欧美一区二区三区黑人| 国产黄色视频一区二区在线观看| 免费高清在线观看视频在线观看| 999精品在线视频| 久久这里只有精品19| 国产日韩欧美在线精品| 伦理电影免费视频| 岛国毛片在线播放| 老司机深夜福利视频在线观看 | 国产日韩欧美在线精品| 久久久久精品国产欧美久久久 | 免费黄频网站在线观看国产| 又大又黄又爽视频免费| 亚洲专区中文字幕在线 | 一区二区三区乱码不卡18| 大陆偷拍与自拍| a级毛片黄视频| 久久性视频一级片| xxx大片免费视频| 欧美97在线视频| 亚洲欧洲精品一区二区精品久久久 | 亚洲熟女毛片儿| 久久久久精品久久久久真实原创| 在线观看人妻少妇| 欧美另类一区| 天堂中文最新版在线下载|