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

    近30年福建省土地利用時空演變的海陸梯度分異規(guī)律

    2023-08-08 09:01:16吳燕芳陳耀亮陸燈盛
    熱帶地理 2023年7期
    關(guān)鍵詞:分異林地梯度

    吳燕芳,陳耀亮,陸燈盛

    (福建師范大學(xué) a.地理研究所;b.濕潤亞熱帶山地生態(tài)國家重點實驗室培育基地,福州 350007)

    土地利用變化是在人類活動和自然因素驅(qū)動下地表類型及其結(jié)構(gòu)和功能發(fā)生或恢復(fù)性、或退化性、或徹底的改變(史培軍 等,2000;羅格平 等,2005)。近現(xiàn)代以來,地球表面的變化主要是由人類社會經(jīng)濟活動所導(dǎo)致的,是人類活動與生態(tài)環(huán)境相互作用的重要體現(xiàn)(劉紀遠 等,2018)。隨著中國工業(yè)化和城鎮(zhèn)化進程的加速推進,土地利用類型發(fā)生明顯的時空變化,使得城鎮(zhèn)、農(nóng)業(yè)、生態(tài)空間矛盾愈演愈烈,在一定程度上限制了國土空間的可持續(xù)發(fā)展(范業(yè)婷 等,2019)。因此,探究土地利用時空演變特征,可為土地資源的可持續(xù)發(fā)展提供有效參考。

    目前國內(nèi)外學(xué)者在不同的數(shù)據(jù)源、研究尺度的基礎(chǔ)上對土地利用變化進行了大量研究,已形成一套比較成熟的研究方法與理論體系(劉紀遠 等,2003;Lu et al., 2013)。大多研究以遙感影像數(shù)據(jù)、長時間序列土地利用分類數(shù)據(jù)集等為數(shù)據(jù)源(張宇碩 等,2015;Pan et al., 2017; Jiang et al., 2018),以全球、國家、城市等為研究尺度(杜國明 等,2015;Feng et al., 2017; Liu et al., 2021),采用轉(zhuǎn)移矩陣、動態(tài)度、空間分析等動態(tài)變化分析方法(張麗 等,2014;付建新 等,2020),主成分分析法、地理探測器等驅(qū)動分析方法(李進濤 等,2018),從時空變化特征、驅(qū)動因素、生境效應(yīng)等方面揭示了土地利用變化的演變規(guī)律及影響作用。如基于中國土地利用數(shù)據(jù)集的分析,發(fā)現(xiàn)中國不同分區(qū)之間土地利用變化存在顯著的差異性,國家政策是引起變化的一個重要因素(劉紀遠 等,2018;Kuang,2020);中國土地利用變化的生態(tài)環(huán)境效應(yīng)存在顯著的空間分異性(陳萬旭 等,2019)。

    緩沖區(qū)分析法可將研究區(qū)域以某一中心點或邊界等距劃分為一定數(shù)量的子區(qū)域,是目前地理學(xué)領(lǐng)域中常被廣泛應(yīng)用于土地利用變化研究的一種梯度分析法(王德智 等,2014;于濤 等,2019)。特別是沿海地區(qū),相關(guān)研究從不透水地表(ISA)擴張、土地利用景觀格局變化、海岸帶土地利用變化等方面探究沿海區(qū)域土地利用變化的海陸梯度分異規(guī)律(Kuang et al., 2014; Lu et al., 2018;孟紫琪 等,2018;Yan et al., 2020;賈艷艷 等,2020)。部分ISA 擴張研究結(jié)合研究區(qū)地理特征,以河流岸線、海岸線為基準進行緩沖區(qū)分析,如陳明發(fā)等(2021)研究得出粵港澳大灣區(qū)的城市建設(shè)以珠江支流與海岸線為軸心,在縱向和橫向上具有明顯的擴張?zhí)荻忍卣鳎灰灿袑W(xué)者采用景觀格局指數(shù)與緩沖帶分析等方法對福建東山島海陸梯度下的景觀破碎化變化特征進行分析,發(fā)現(xiàn)距離海岸線500和1 000 m 處是東山島沿海人為干擾的強烈地帶(巫麗蕓等,2020)。大多研究表明不同行政區(qū)劃單元下的海岸帶土地利用變化在空間上存在顯著的海陸梯度分異規(guī)律(吳莉 等,2015;Ding et al., 2019),如中國海岸帶不同土地利用類型之間的海陸梯度特征不同,其變化存在省域差異(侯西勇 等,2011);江蘇省海岸帶土地利用多樣性和土地利用程度沿陸海方向分別呈現(xiàn)“低—高—低”“高—中—低”的不同規(guī)律,且南北梯度規(guī)律差異顯著(許艷 等,2012);煙臺市海岸帶土地利用形態(tài)在空間上呈現(xiàn)類似圈層形狀的海陸梯度分異特征(鄒敏 等,2016);青島市海岸帶港口區(qū)域的土地利用格局在空間上沿海陸方向按一定順序位置分布,且不同時期土地開發(fā)強度梯度特征會發(fā)生變化(Yan et al.,2021);長三角地區(qū)海岸帶的土地利用動態(tài)度存在明顯的建設(shè)用地向海、海洋用地向陸變化規(guī)律(王曼曼 等,2020)。綜上所述,大多數(shù)對于土地利用變化海陸梯度分異規(guī)律的研究多限于海岸帶、島域、單一城市等中小尺度,僅能分析局部梯度變化規(guī)律,對于整個省域視角的土地利用變化海陸梯度分異規(guī)律的研究還比較少。

    福建省地處中國東南沿海,自改革開放以來,全省社會經(jīng)濟建設(shè)從沿海到內(nèi)陸地區(qū)取得不同程度的成就,城鎮(zhèn)化程度不斷提升,對福建省的土地利用空間結(jié)構(gòu)產(chǎn)生很大的影響。近年來,國內(nèi)學(xué)者對福建省土地利用變化展開了相關(guān)研究,但現(xiàn)有研究大多以行政區(qū)劃為研究單元分析福建省沿海到內(nèi)陸地區(qū)的土地利用變化時空演變規(guī)律、驅(qū)動機制及生境效應(yīng)(曲勝秋 等,2020;李勝鵬 等,2020),較少以緩沖區(qū)梯度帶為研究單元探索福建省沿海陸方向上土地利用變化趨勢及其空間分異規(guī)律。以緩沖帶為研究單元分析整個省域不同梯度帶下土地利用變化空間分異規(guī)律,能準確直觀地反映土地利用變化、分布特征,可為政府部門和國土規(guī)劃部門制定科學(xué)合理的土地管理政策提供依據(jù);作為實例研究,還能豐富福建省土地利用實證分析研究,為推進福建省土地利用類型結(jié)構(gòu)調(diào)整及布局優(yōu)化提供參考。

    因此,本文以福建省為研究區(qū),基于1990、2000、2010、2020年4期的土地利用數(shù)據(jù)集,以海岸線為基準線,采用緩沖區(qū)分析法與土地利用動態(tài)變化指標,分析沿海陸方向不同梯度帶的土地利用類型的時空演變特征,并結(jié)合格網(wǎng)單元探究海陸梯度下福建省土地利用變化的空間分異規(guī)律。本文擬回答以下問題:1)沿海陸方向,近30年研究區(qū)各土地利用類型是否存在梯度變化規(guī)律?2)不同梯度帶呈現(xiàn)的土地利用時空格局是如何變化的?3)海陸梯度下土地利用變化是否具有空間分異特征?

    1 研究區(qū)概況與數(shù)據(jù)來源

    1.1 研究區(qū)概況

    福建省位于中國東南沿海,陸地總面積12.14萬km2。地勢西北高,東南低;境內(nèi)山地、丘陵面積約占全省總面積的80%以上,西部和中部形成兩大山帶,其間為互不貫通的河谷、盆地,東部沿海為丘陵、臺地和濱海平原。全省屬亞熱帶海洋性季風氣候,但氣候區(qū)域差異較大,閩東南沿海地區(qū)屬南亞熱帶氣候,其他地區(qū)屬中亞熱帶氣候;全省70%區(qū)域≥10℃的積溫為5 000~7 600℃,年均溫17~21℃,平均降雨量1 400~2 000 mm,是中國雨量最豐富的省份之一。全省水系密布,主要河流有閩江、九龍江、晉江、交溪、汀江;森林覆蓋率連續(xù)42 年為全國第一。2020 年末全省常住人口為4 161萬人,其中,城鎮(zhèn)常住人口2 861萬人,常住人口城鎮(zhèn)化率為68.8%;同年,福建省GDP達到4.39萬億元;到2020 年,全省基本形成鐵路網(wǎng)和公路網(wǎng),2020年末全省鐵路營業(yè)里程和公路通車里程分別為3 774 和11 萬km(福建省統(tǒng)計局,2021)。研究區(qū)范圍包括福州市、廈門市、泉州市、漳州市、莆田市、寧德市、南平市、三明市和龍巖市9個設(shè)區(qū)市,共有29 個市轄區(qū),12 個縣級市,43 個縣(不包含金門縣)。

    1.2 數(shù)據(jù)收集和整理

    1.2.1 福建省土地利用數(shù)據(jù)的獲取 1990、2000、2010、2020年4期的福建省土地利用矢量數(shù)據(jù)集來源于中國科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心①https://www.resdc.cn。該數(shù)據(jù)集是劉紀遠團隊基于Landsat TM/ETM+數(shù)據(jù),采用人機交互解譯的方法研發(fā)的全國土地利用現(xiàn)狀1∶10萬矢量數(shù)據(jù)庫,經(jīng)大量野外調(diào)查驗證,其產(chǎn)品總體精度達到90%以上(劉紀遠 等,2003,2018;Kuang et al., 2016;匡文慧 等,2022)。根據(jù)研究區(qū)的實際情況,采用的土地利用分類體系參考中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心的中國土地利用遙感監(jiān)測分類體系②https://www.resdc.cn/DOI/DOI.aspx?DOIID=54,將福建省土地利用類型劃分為6個一級類,即耕地、林地、草地、水域、城鄉(xiāng)建設(shè)用地和未利用地。

    1.2.2 政策事件收集整理 主要通過中央人民政府門戶網(wǎng)站③http://www.gov.cn或福建省人民政府門戶網(wǎng)站④http://www.fujian.gov.cn收集16個主要區(qū)域發(fā)展戰(zhàn)略及政策,并根據(jù)研究時段進行匯總整理。

    2 研究方法

    2.1 梯度分析法

    基于福建省海岸線,采用緩沖區(qū)梯度帶法(侯西勇 等,2011;陳明發(fā) 等,2021;Yan et al.,2021),分析該省海陸梯度下土地利用變化的時空演變規(guī)律及空間分異規(guī)律。首先,根據(jù)2020 年Landsat8 OLI 影像進行目視解譯提取該省海岸線,其中,海岸線直接以水陸分界線為解譯原則,不將平潭島、東山島等島嶼包含在內(nèi);其次,借鑒已有研究(陳明發(fā) 等,2021),并結(jié)合福建省實際情況,以該海岸線為基準按照20 km 間隔向陸地一側(cè)進行緩沖分析,最終沿海陸方向生成15個梯度帶,在此基礎(chǔ)上劃分三大梯度帶,分別為沿海梯度帶(距海岸線0~60 km)、過渡梯度帶(距海岸線60~180 km)和內(nèi)陸梯度帶(距海岸線180~300 km),如圖1所示;第三,結(jié)合4期土地利用數(shù)據(jù)集,疊加分析不同梯度的土地利用的數(shù)量變化規(guī)律;第四,利用轉(zhuǎn)移矩陣,分析不同時期海陸梯度下土地利用類型的空間變化分布特征,以及通過動態(tài)度分析不同土地利用類型的梯度變化幅度及土地利用變化的劇烈程度;最后,以2 km × 2 km 的矢量格網(wǎng)作為研究單元、綜合土地利用動態(tài)度為觀測指標,利用空間自相關(guān)、熱點分析等方法,分析1990—2020年該省海陸梯度下土地利用變化的冷點和熱點在空間分布上的差異性。

    圖1 福建省梯度帶劃分Fig.1 The buffer zone division in Fujian Province

    2.2 土地利用動態(tài)分析法

    2.2.1 土地利用類型的面積變化分析 面積變化量是反映研究期內(nèi)某一種土地利用類型的面積變化幅度(劉鳳蓮 等,2021)。表達式為:

    式中:ΔC為單一土地利用類型的面積變化量,ΔC> 0 表示面積增加,ΔC< 0 表示面積減少;Ua和Ub分別表示研究期初和期末土地利用類型的面積。

    2.2.2 轉(zhuǎn)移矩陣分析 轉(zhuǎn)移矩陣不僅可以體現(xiàn)研究區(qū)各土地利用類型互相轉(zhuǎn)移的方向和面積變化情況,還可以從空間角度直觀反映不同土地利用類型在不同時期的格局分布規(guī)律(孫智斌 等,2018)?;?期的土地利用數(shù)據(jù)集,利用ArcGIS的空間分析和統(tǒng)計分析功能,分別對臨近的2期數(shù)據(jù)進行空間疊加分析,得到1990—2020年不同時期的土地利用轉(zhuǎn)移矩陣。

    2.2.3 動態(tài)度分析 結(jié)合轉(zhuǎn)移矩陣分析,采用土地利用動態(tài)度對1990—2020年福建省土地利用類型動態(tài)變化進行分析。土地利用動態(tài)度分為單一土地利用動態(tài)度和綜合土地利用動態(tài)度,二者可用于反映在一定時期內(nèi)研究區(qū)土地利用類型在數(shù)量上的動態(tài)變化特征(王宏亮 等,2017;付建新 等,2020)。

    單一土地利用動態(tài)度(K)反映一定時期內(nèi)研究區(qū)某一種土地利用類型的變化速率。表達式為:

    式中:T表示研究時長(年)。

    綜合土地利用動態(tài)度反映在一定時期內(nèi)研究區(qū)內(nèi)整體土地利用類型變化的劇烈程度。表達式為:

    式中:LC 為區(qū)域內(nèi)的綜合土地利用動態(tài)度;ΔLA(i,t0)為研究期初第i類土地利用類型的總面積(km2);ΔLA(i,t1-t0)為研究期初至期末第i類土地利用類型轉(zhuǎn)為其他類型面積的總和。

    2.3 海陸梯度土地利用變化空間分異性分析

    為探究1990—2020年福建省海陸梯度下土地利用變化的空間分異性,選取綜合土地利用動態(tài)度指數(shù)作為觀測指標,但由于福建省各縣市區(qū)行政區(qū)劃之間面積大小存在明顯差異,為了消除面積差異帶來的影響,結(jié)合研究區(qū)實際情況,以2 km×2 km格網(wǎng)作為研究單元對其進行分析。

    2.3.1 空間自相關(guān) 空間自相關(guān)是對整個研究區(qū)某一觀測值在空間上是否存在聚集性的顯著性檢驗分析(郭椿陽 等,2019)。采用全局Moran'sI指數(shù)對綜合土地利用動態(tài)度指數(shù)進行全局空間自相關(guān)分析,基于其值大小從全局角度判斷沿海陸方向土地利用變化是否具有空間聚集性。其中,I>0,表示存在正的空間自相關(guān),在空間上表現(xiàn)為聚集特征且值越接近1,聚集現(xiàn)象越明顯;I<0,表示存在負的空間自相關(guān),在空間上表現(xiàn)為分散分布;I=0,表示不存在空間相關(guān)性。

    2.3.2 土地利用變化熱點分析 為進一步探究局部相鄰區(qū)域之間土地利用變化的空間分異性,采用局部自相關(guān)分析中的熱點分析(Gi*指數(shù))測度其空間差異特征,分析沿海陸方向的土地利用變化是否具有統(tǒng)計顯著意義上的冷熱點聚集現(xiàn)象(劉光旭 等,2021)。其中,冷熱點聚集主要根據(jù)表示概率的P值、標準差的倍數(shù)的z得分和置信度進行顯著性聚類,結(jié)合綜合土地利用動態(tài)度指標,將研究區(qū)劃分為5 個區(qū)域(陳萬旭 等,2019):冷點區(qū)、次冷點區(qū)、次熱點區(qū)、熱點區(qū)和非顯著區(qū)。

    3 結(jié)果分析

    3.1 海陸梯度土地利用變化時空格局分析

    3.1.1 海陸梯度土地利用數(shù)量變化分析 1990—2020年福建省土地利用變化沿海陸方向呈現(xiàn)明顯的梯度變化特征(圖2)。沿海梯度帶的各土地利用類型變化較為劇烈,其中,20 km 梯度帶內(nèi)的變化最為劇烈,以耕地和城鄉(xiāng)建設(shè)用地變化為主;過渡梯度帶的各土地利用類型變化程度明顯小于沿海梯度帶,以林地、草地和城鄉(xiāng)建設(shè)用地3種類型變化為主;內(nèi)陸梯度帶的各土地利用類型變化相對趨于穩(wěn)定。

    圖2 福建省海陸梯度主要土地利用類型面積變化特征Fig.2 Area change characteristics of main land use types under the sea-land gradient in Fujian Province

    近30年各土地利用類型的面積變化量隨離海岸線距離的增加具有不同的變化趨勢(圖2-a)。其中,耕地和城鄉(xiāng)建設(shè)用地的面積變化量從沿海到內(nèi)陸呈持續(xù)減少的趨勢,面積顯著變化的區(qū)域主要集中在沿海梯度帶,面積變化量遠大于其他梯度帶,特別是距海岸線20 km 范圍,耕地和城鄉(xiāng)建設(shè)用地的面積變化量均達到最大,分別為-1 317.12和1 799.8 km2;而過渡梯度帶中城鄉(xiāng)建設(shè)用地出現(xiàn)2處小峰值(120和180 km),其增加量明顯高于周邊,分別為177.85 和195.49 km2,這主要是由龍巖、三明和南平市轄區(qū)建成區(qū)擴張引起的。林地和草地的面積變化量呈現(xiàn)“減少—增加—減少”的梯度變化趨勢,林地增加量顯著增多的區(qū)域主要集中在過渡梯度帶和內(nèi)陸梯度帶的140~240 km 范圍內(nèi),而草地在距海岸線20 km 范圍的減少量最大,為362.86 km2;水域整體變化幅度不大,僅在距海岸線20 km范圍有較為明顯的減少。

    從不同時期看,各土地利用類型變化存在明顯的時序特征。林地和草地面積變化量呈持續(xù)減少的特征,2000 年以后顯著降低,尤其是1990—2000年二者的變化量均為最大,分別為1 831.21 和-1 851.76 km2,其中140~220 km 是林地和草地顯著變化的核心區(qū)域(圖2-b),該區(qū)域二者的變化量分別達到1 144.36 和-1 163.76 km2;耕地和城鄉(xiāng)建設(shè)用地的面積變化量呈先增加后減少的特征,2000 年以后二者的變化量顯著上升,且梯度分異規(guī)律更加明顯,特別是2000—2010 年變化量最大,分別為-1 082.8和1 869.35 km2(圖2-c),該時期二者的梯度變化規(guī)律與1990—2020年較為一致,其中城鄉(xiāng)建設(shè)用地在120 和180 km兩處梯度帶的增加量均在120 km2左右;水域面積變化量總體呈微弱增加,其中,2000—2010年水域以增加為主,在過渡梯度帶的120~140 km 范圍內(nèi)有明顯的增加(圖2-c),而2010—2020年水域在沿海20 km范圍內(nèi)具有顯著的減少量,達176.93 km2(圖2-d)。

    3.1.2 海陸梯度轉(zhuǎn)移矩陣分析 1990—2020 年前10個土地利用變化類型的轉(zhuǎn)移面積占總轉(zhuǎn)移面積的比例隨與海岸線距離的增加而降低(圖3-a),沿海梯度帶(44.23%)>過渡梯度帶(37.31%)>內(nèi)陸梯度帶(15.03%);土地利用變化類型隨與海岸線距離的增加,由以耕地轉(zhuǎn)城鄉(xiāng)建設(shè)用地為主逐漸向以草地轉(zhuǎn)林地、林地轉(zhuǎn)草地為主發(fā)生轉(zhuǎn)移(圖3-b)。

    圖3 1990—2020年不同梯度帶前10個土地利用變化類型轉(zhuǎn)移變化比例Fig.3 Top 10 land use change types under different gradient zones in the transfer change ratio from1990 to 2020

    從不同時期看,前10個土地利用變化類型的轉(zhuǎn)移變化比例具有顯著的時序特征和梯度分異規(guī)律(見圖3-a),表現(xiàn)為以60 km梯度帶為分界帶,沿海梯度帶內(nèi)的轉(zhuǎn)移變化比例隨著時間的推移而顯著升高,而過渡梯度帶和內(nèi)陸梯度帶則是下降,逐漸接近于0。其中,1990—2000年,前10個類型的轉(zhuǎn)移變化比例整體偏低(12%以下),呈現(xiàn)“降低—升高—降低”的梯度變化趨勢(見圖3-a);該時期的土地利用變化類型以林地增加、草地和耕地減少為主,各轉(zhuǎn)移變化類型存在不同的梯度變化特征(圖3-c)。具體而言,草地轉(zhuǎn)林地比例呈現(xiàn)隨與海岸線距離增加而升高,主要集中在140~220 km 梯度帶內(nèi),這可能與20 世紀90 年代初的“三五七”造林綠化工程等生態(tài)環(huán)境建設(shè)有關(guān);耕地轉(zhuǎn)城鄉(xiāng)建設(shè)用地僅在沿海梯度帶內(nèi)有部分占比,特別是距海岸線20 km的沿海地區(qū)。

    2000—2010 年,前10 個類型的轉(zhuǎn)移變化比例整體呈持續(xù)下降的梯度變化趨勢(見圖3-a);該時期的土地利用變化類型以城鄉(xiāng)建設(shè)用地急速增加,耕地、林地和草地減少為主(圖3-d)。具體而言,城鄉(xiāng)建設(shè)用地轉(zhuǎn)入類型主要來源于耕地(30.91%)、林地(17.54%)和草地(7.21%),這主要是由于沿海開放經(jīng)濟區(qū)以及海峽西岸經(jīng)濟區(qū)等區(qū)域發(fā)展戰(zhàn)略的實行,促使大范圍的城市擴張;但不同轉(zhuǎn)入類型存在明顯的梯度分異性,耕地轉(zhuǎn)城鄉(xiāng)建設(shè)用地比例主要表現(xiàn)為波動降低,主要集中在沿海梯度帶,林地轉(zhuǎn)城鄉(xiāng)建設(shè)用地比例呈現(xiàn)先升高后降低,主要集中在過渡梯度帶;在120~140、260~300 km 梯度帶內(nèi)出現(xiàn)小幅度的水域增加,這是由于該時期有大型水庫的建成使用,使得耕地、林地轉(zhuǎn)為水域。

    2010—2020 年,前10 個類型的轉(zhuǎn)移變化比例整體仍表現(xiàn)為持續(xù)下降的趨勢(圖3-a);該時期的土地利用變化類型以城鄉(xiāng)建設(shè)用地增加、林地小幅度增加、耕地和草地減少為主(圖3-e)。具體而言,耕地(43.63%)、林地(11.66%)、水域(9.43%)和草地(5.79%)向城鄉(xiāng)建設(shè)用地轉(zhuǎn)變,過渡及內(nèi)陸梯度帶中的耕地、林地、草地轉(zhuǎn)城鄉(xiāng)建設(shè)用地占比均顯著升高,水域的轉(zhuǎn)入主要集中在距海岸線20 km 的范圍內(nèi),由于沿海地區(qū)城市擴張引起的圍填海造地,使得一定量的水域面積減少;沿海梯度帶內(nèi)有明顯較高比例的耕地、草地轉(zhuǎn)林地,這是由于現(xiàn)代特色農(nóng)業(yè)產(chǎn)業(yè)、農(nóng)業(yè)產(chǎn)業(yè)園的建設(shè),使得草地和耕地大面積轉(zhuǎn)為果園、茶園和苗圃地等,致使林地面積增加。

    從空間分布看,不同梯度內(nèi)的土地利用變化類型存在顯著的南北差異(圖4)。距海岸線20 km范圍內(nèi)是耕地減少、城鄉(xiāng)建設(shè)用地增加的最主要區(qū)域,特別是東南部地區(qū)的沿海城市,如泉州市轄區(qū)、晉江市、石獅市等(圖4-d2);2000 年以后,過渡梯度帶部分地區(qū)的城鄉(xiāng)建設(shè)用地顯著增加,如龍巖市轄區(qū)(120 km 梯度帶),以及三明市轄區(qū)、永安市、建陽區(qū)的建成區(qū)(180 km 梯度帶)(圖4-d),這是由于主要交通干線的建成使用促使其城市邊緣發(fā)生擴張,而建陽區(qū)后期的增加(圖4-d1)也與2014年國務(wù)院批復(fù)南平市行政區(qū)劃調(diào)整,撤市設(shè)區(qū)的政策有關(guān)。過渡梯度帶、內(nèi)陸梯度帶北部地區(qū)是林地增加、草地減少的主要區(qū)域,特別是2000年以前,如位于建甌市、將樂縣等地區(qū)(圖4-b1、c1);2010—2020 年,在沿海梯度帶南部地區(qū)(漳州市的漳浦縣、云霄縣、平和縣,圖4-b2)存在顯著的林地增加。

    圖4 1990—2020年福建省海陸梯度主要土地利用類型變化的空間分布(a.耕地;b.林地;c.草地;d.城鄉(xiāng)建設(shè)用地)Fig.4 Spatial distribution of main land use types changes under land-sea gradient in Fujian Province from 1990 to 2020(a.cultivated land; b.forest land; c.grassland; d.urban and rural construction land)

    3.2 海陸梯度下土地利用空間分異特征分析

    3.2.1 海陸梯度下的動態(tài)度變化分析 1990—2020

    年不同梯度帶間的綜合土地利用動態(tài)度呈現(xiàn)“降低—升高—降低”變化規(guī)律,但沿海梯度帶>過渡梯度帶>內(nèi)陸梯度帶(表1),其中,距海岸線20 km范圍綜合動態(tài)度最高,為0.75%;單一動態(tài)度上(圖5-a),近30 年內(nèi)城鄉(xiāng)建設(shè)用地動態(tài)度具有顯著的梯度變化特征,呈現(xiàn)先升高后降低的態(tài)勢,其變化速率遠高于其他類型,這主要是由于其他3種類型的基數(shù)遠遠大于城鄉(xiāng)建設(shè)用地;耕地動態(tài)度呈現(xiàn)持續(xù)降低的趨勢,其中距海岸線20 km 范圍的耕地減少速率最大,為2.22%;草地動態(tài)度呈先降低后升高的趨勢,在220 km 梯度帶內(nèi)變化速率(-1.95%)達到最大;林地動態(tài)度最為穩(wěn)定,變化速率最低,這是由于全省林地面積基數(shù)最大。

    表1 1990—2000年福建省海陸梯度綜合土地利用動態(tài)度Table 1 Comprehensive dynamic degree of land use of sea-land gradient in Fujian Province from 1990 to 2000 %

    圖5 1990—2020年福建省海陸梯度單一土地利用類型動態(tài)度Fig.5 Dynamic degree of single land use type of sea-land gradient in Fujian Province from 1990 to 2020

    從不同時期看,1990—2000年,各梯度帶的綜合動態(tài)度平均維持在0.24%,但過渡梯度帶>內(nèi)陸梯度帶>沿海梯度帶(見表1),這是由于該時期林地、草地、耕地在過渡和內(nèi)陸梯度帶中大幅度相互轉(zhuǎn)化引起;單一動態(tài)度上,林地動態(tài)度(0.34%~0.48%)、草地動態(tài)度(1.60%左右)在140~240 km梯度帶明顯高于其他類型,城鄉(xiāng)建設(shè)用地(0.69%~1.40%)則在沿海梯度帶變化速率較高(圖5-b)。2000—2010年,綜合動態(tài)度較前期具有顯著的梯度分異規(guī)律,呈現(xiàn)持續(xù)減少的特征,但內(nèi)陸梯度帶>沿海梯度帶>過渡梯度帶(見表1),這是由于內(nèi)陸梯度帶區(qū)域面積明顯小于其他區(qū)域,但發(fā)生轉(zhuǎn)移變化面積與其他區(qū)域相差不大;單一動態(tài)度上,城鄉(xiāng)建設(shè)用地動態(tài)度呈現(xiàn)先升高后降低的特征,基本都達到6%以上的變化速率(圖5-c),是變化速率最高的時期,而耕地動態(tài)度在沿海梯度帶明顯高于林地和草地,距海岸線20 km 范圍內(nèi)達到最高變化速率,為-1.16%。2010—2020年,整體綜合動態(tài)度顯著降低,僅沿海梯度帶(0.18%)有較劇烈的土地利用變化,特別是距海岸線0~40 km 范圍,其他梯度帶的變化程度基本保持穩(wěn)定(見表1);單一動態(tài)度上,沿海梯度帶的城鄉(xiāng)建設(shè)用地變化速率明顯小于過渡、內(nèi)陸梯度帶的變化速率;其他類型的動態(tài)度與2000—2010年的趨勢相對一致(圖5-d)。

    3.2.2 空間自相關(guān)分析 福建省各個時期的全局Moran'sI值在0.46~0.52(表2),表明土地利用變化在空間上具有顯著的聚集性;3 個時期的全局Mo‐ran'sI值先升高后降低,表明不同階段土地利用變化的空間聚集效應(yīng)呈先增強后減弱的特征。從各梯度帶的全局Moran'sI值看,1990—2020 年,20、40、120、140、180 、240 km梯度帶的Moran'sI值遠大于其他梯度帶,說明這些梯度帶內(nèi)的土地利用變化空間分布具有更為顯著的聚集性;距海岸線0~40 km 范圍3 個時期的Moran'sI值在0.45~0.50,說明這2個梯度帶的土地利用變化始終保持著較強的空間聚集效應(yīng),而其他梯度帶則隨時間表現(xiàn)出不同程度的空間聚集性,如120 km梯度帶呈現(xiàn)先增強后減弱的空間聚集效應(yīng),140和180 km梯度帶的空間聚集效應(yīng)則是持續(xù)減弱。

    3.2.3 海陸梯度下冷熱點分析 根據(jù)海陸梯度下冷熱點的空間聚集結(jié)果(圖6)可知,空間上,各個時期的福建省土地利用變化沿海陸方向均有顯著的空間分異特征;時間上,以2000年為時間分界點,熱點區(qū)域大致呈現(xiàn)“由內(nèi)陸地區(qū)向沿海地區(qū)轉(zhuǎn)移”、由“北多南少”向“南多北少”轉(zhuǎn)變的空間聚集特征,2000年以后,熱點區(qū)域逐漸向沿海梯度帶聚集分布,過渡、內(nèi)陸梯度帶的熱點區(qū)域則是少數(shù)且分散分布。總體上,1990—2020年,全省土地利用變化具有明顯的冷點與熱點分區(qū)(圖6-a),其中,土地利用變化劇烈的熱點區(qū)域主要集中在沿海梯度帶,呈現(xiàn)“南多北少”的特征,過渡、內(nèi)陸梯度帶則相反,分散分布在主要縣市區(qū),如三明市轄區(qū)、建陽區(qū)等;冷點區(qū)域主要分布在過渡梯度帶的中部地區(qū),如永泰縣、德化縣等,而內(nèi)陸梯度帶主要分布在中、南部的邊界地區(qū),如武平縣、長汀縣等,其變化程度均較為平緩或穩(wěn)定。

    圖6 1990—2020年福建省海陸梯度下土地利用變化冷熱點空間分布Fig.6 Spatial distribution of cold and hot spots of land use change under the land-sea gradient in Fujian Province from 1990 to 2020

    從不同時期看,1990—2000年,熱點區(qū)域呈現(xiàn)“北多南少”的空間分異特征,其中熱點區(qū)所占比例較大、次熱點區(qū)分布在其四周,熱點區(qū)域主要分布在福建省偏北部地區(qū),特別是過渡梯度帶區(qū)域,如柘榮縣、閩清縣、古田縣、建甌市、將樂縣等,這是由于這些區(qū)域的類型轉(zhuǎn)移多且面積變化較大(見圖4),土地利用變化程度劇烈,從而綜合土地利用動態(tài)度值較高;冷點區(qū)域在空間上呈現(xiàn)“四周多中部少”的特征,其中次冷點區(qū)所占比例大,主要聚集在沿海梯度帶的中部、過渡梯度帶的南部邊界地區(qū)(圖6-b)。

    2000—2010年,全省土地利用變化只呈現(xiàn)熱點區(qū)域的空間分布特征,這很可能是與該時期城鄉(xiāng)建設(shè)用地和耕地大面積變化,而其他類型變化相對不明顯有關(guān)(見圖2、5)。該時期的熱點區(qū)域呈現(xiàn)顯著的“南多北少”的空間分異特征(圖6-c)。相較1990—2000年,該期間的核心熱點區(qū)域由西北部的過渡梯度帶轉(zhuǎn)移到東南部的沿海梯度帶,呈現(xiàn)“一核心多分散”的分布特征。其中,“一個核心”的熱點區(qū)域主要集中分布在距海岸線20 km 范圍的廈漳泉地區(qū),包括晉江市、石獅市、惠安縣、廈門市轄區(qū)、漳州市轄區(qū)等;“多個分散”的熱點區(qū)域則位于沿海梯度帶北部的福州市轄區(qū),過渡梯度帶中、南部的德化縣、永安市、沙縣、新羅區(qū)等地區(qū)。

    2010—2020 年,熱點區(qū)域的空間分異特征與2000—2010 年相似,但海陸差異更為顯著(圖6-d)。熱點區(qū)域更加聚集分布在沿海梯度帶,呈大范圍、均勻分布;而過渡梯度帶呈現(xiàn)“北多南少”的分散分布特征,主要分布在過渡梯度帶的主要縣市區(qū)建成區(qū),包括南平市轄區(qū)、建甌市、三明市轄區(qū)等。該時期整體的熱點區(qū)域分布明顯發(fā)生“由南往北”轉(zhuǎn)移,其中沿海梯度帶表現(xiàn)為福州市轄區(qū)、莆田市轄區(qū)的熱點區(qū)域明顯擴大,而廈漳泉地區(qū)明顯縮??;過渡梯度帶表現(xiàn)為武夷山市、建陽區(qū)、建甌市等出現(xiàn)小范圍熱點聚集區(qū),而南部地區(qū)無顯著的聚集區(qū)。

    4 土地利用變化空間分異原因分析

    沿海梯度帶的土地利用變化隨時間呈現(xiàn)“小幅度變化—急劇變化—減緩變化”的特征,具有顯著的空間分異規(guī)律,這與福建省沿海地區(qū)城市發(fā)展在政策方面受到很大支持有關(guān)(圖7)。自1985 年,國家支持閩南廈漳泉地區(qū)建設(shè)發(fā)展,并列入沿海經(jīng)濟開放區(qū),閩南金三角地區(qū)率先發(fā)展(李永實,2001),使得耕地等類型被占用,城鄉(xiāng)建設(shè)用地開始小幅度增加;隨后2004年國家提出了“海峽西岸經(jīng)濟區(qū)”戰(zhàn)略設(shè)想,并在2009年由區(qū)域發(fā)展戰(zhàn)略上升為國家戰(zhàn)略,這10年間城鎮(zhèn)化加速推進,開啟了以東南沿海的廈漳泉區(qū)域為核心、主要交通干線為紐帶,輻射帶動過渡梯度帶南部交通沿線附近的建成區(qū)(如龍巖市的新羅區(qū))加速發(fā)展的局面(白如山 等,2009);2010—2020 年,“中國(福建)自由貿(mào)易試驗區(qū)”“21世紀海上絲綢之路核心區(qū)”“福州新區(qū)設(shè)立”“福廈泉國家自主創(chuàng)新示范區(qū)”等一些系列國家政策支持,以及“三縱六橫”鐵路網(wǎng)、“六縱十橫”高速公路網(wǎng)已基本形成使用,全省區(qū)域發(fā)展重心逐漸從廈漳泉地區(qū)轉(zhuǎn)移向福州地區(qū),形成以福州地區(qū)為核心,通過交通路網(wǎng)連通廈漳泉地區(qū)協(xié)同發(fā)展,帶動福州市周圍及過渡梯度帶北部主要縣市區(qū)的建成區(qū)(如建陽區(qū)等)建設(shè)發(fā)展(李勝鵬 等,2020);而自十八大以來,中國經(jīng)濟不斷推進供給側(cè)結(jié)構(gòu)性改革,著重生態(tài)建設(shè)發(fā)展,城鎮(zhèn)化進程進入發(fā)展轉(zhuǎn)型期,全省的城鄉(xiāng)建設(shè)用地增速明顯較前期有所減緩,特別是廈漳泉地區(qū)。

    圖7 20世紀80年代以來福建省主要區(qū)域發(fā)展戰(zhàn)略及政策匯總Fig.7 Summary of major regional development strategies and policies in Fujian Province since the 1980s

    過渡梯度帶的土地利用變化隨時間變化具有明顯的區(qū)域差異性。早期過渡梯度帶具有較為劇烈的土地利用變化,這是由于1981年“大念山海經(jīng)”到1989年“三五七造林綠化工程”等生態(tài)環(huán)境建設(shè)政策的實施(見圖7),極大地推動了草地和耕地等宜林地、荒山大面積轉(zhuǎn)為林地,使得1990—2020年植被變化相當活躍,林地轉(zhuǎn)入面積大幅增加;自2000年之后,過渡梯度帶的土地利用變化格局由植被變化轉(zhuǎn)變?yōu)槌青l(xiāng)擴張變化,但其南北地區(qū)的城市擴張程度不同,這很可能與沿海梯度帶城市發(fā)展的輻射影響有關(guān)。而內(nèi)陸梯度帶由于距離海岸線最遠,且丘陵山地分布居多,以林地、草地變化為主,城市擴張最不顯著,僅在2010—2020年有小幅度的城鄉(xiāng)建設(shè)用地增加。

    5 結(jié)論與啟示

    基于1990、2000、2010、2020 年土地利用數(shù)據(jù)集,采用緩沖區(qū)梯度帶法、動態(tài)度指標、轉(zhuǎn)移矩陣等方法分析了近30年福建省海陸梯度下的土地利用變化時空演變特征,通過空間自相關(guān)、熱點分析等方法結(jié)合格網(wǎng)單元探究了1990—2020年福建省土地利用變化的海陸梯度分異規(guī)律。得到如下主要結(jié)論:

    1)從面積變化上看,總體上近30 年內(nèi)主要土地利用類型的面積變化隨與海岸線距離的增加具有不同的變化趨勢,其中,耕地和城鄉(xiāng)建設(shè)用地面積變化量呈持續(xù)減少的梯度變化趨勢,距海岸線20 km 范圍二者的面積變化量最大,分別為-1 317.12和1 799.8 km2;而林地和草地面積變化量呈現(xiàn)“減少—增加—減少”的趨勢,集中在140~240 km 梯度帶;水域整體變化幅度不大,僅在距海岸線20 km 范圍有較為明顯的減少。不同時期的各土地利用類型面積變化存在明顯的時序特征。

    2)從轉(zhuǎn)換類型上看,近30年前10個土地利用變化類型的轉(zhuǎn)移面積占總轉(zhuǎn)移變化面積的比例隨著離海岸線距離的增加而降低,沿海梯度帶(44.23%)>過渡梯度帶(37.31%)>內(nèi)陸梯度帶(15.03%);土地利用變化類型隨與海岸線距離增加,由以耕地轉(zhuǎn)城鄉(xiāng)建設(shè)用地為主逐漸向以草地轉(zhuǎn)林地、林地轉(zhuǎn)草地為主發(fā)生轉(zhuǎn)移。空間分布上,不同梯度內(nèi)的土地利用類型變化存在顯著的南北差異??臻g上,沿海梯度帶是耕地減少、城鄉(xiāng)建設(shè)用地增加的主要區(qū)域;過渡梯度帶、內(nèi)陸梯度帶北部地區(qū)是林地增加、草地減少的主要區(qū)域。

    3)從動態(tài)度上看,1990—2020 年不同梯度帶間的綜合土地利用動態(tài)度呈現(xiàn)“降低—升高—降低”變化規(guī)律,但沿海梯度帶>過渡梯度帶>內(nèi)陸梯度帶;單一動態(tài)度上,城鄉(xiāng)建設(shè)用地動態(tài)度沿海陸方向呈現(xiàn)先升高后降低的趨勢,其變化速率遠高于其他類型;耕地動態(tài)度呈現(xiàn)持續(xù)降低,其中距海岸線20 km 范圍的耕地減少速率最大,為2.22%;草地動態(tài)度呈現(xiàn)先降低后升高,在220 km梯度帶內(nèi)變化速率(-1.95%)達到最大;林地動態(tài)度最為穩(wěn)定,變化速率最低,這是由于全省林地面積基數(shù)最大。

    4)1990—2020 年全省土地利用變化在空間上具有顯著的海陸分異規(guī)律,且存在南北區(qū)域差異;時間上,以2000 年為時間分界點,熱點區(qū)域呈現(xiàn)“由內(nèi)陸地區(qū)向沿海地區(qū)轉(zhuǎn)移”、由“北多南少”向“南多北少”轉(zhuǎn)變的空間聚集特征。1990—2020 年福建省土地利用變化的海陸梯度分異特征與福建省各階段的區(qū)域發(fā)展戰(zhàn)略及生態(tài)環(huán)境建設(shè)等政策息息相關(guān)。

    目前關(guān)于土地利用變化的海陸梯度分異規(guī)律的研究大多基于特定緩沖區(qū)域或行政區(qū)劃單元進行緩沖分析,而針對沿海省域尺度的海陸梯度演變規(guī)律的研究較少。以緩沖梯度帶作為研究單元,分析省域尺度下土地利用變化的梯度分異規(guī)律是本文的一個創(chuàng)新點,具有拓寬研究思路的作用。本文結(jié)果對新時代新型城鎮(zhèn)化發(fā)展背景下福建省國土空間規(guī)劃具有重要參考意義。

    結(jié)合本研究,就新時代福建省國土空間開發(fā)保護方面提出建議:①因地制宜,保護國土資源。對于沿海梯度帶區(qū)域,要切實保護沿海地區(qū)集中連片的優(yōu)質(zhì)耕地資源,嚴格遵守耕地保護紅線,協(xié)調(diào)好耕地與城鄉(xiāng)建設(shè)用地之間的轉(zhuǎn)換關(guān)系;同時,要加強海域使用管理,兼顧沿海梯度帶區(qū)域的海洋資源的大力保護與適度開發(fā)。對于過渡、內(nèi)陸梯度帶區(qū)域,除了落實耕地保護,還要對林地資源進行嚴格管控與保護,加強對優(yōu)質(zhì)林地資源的可持續(xù)經(jīng)營管理,并積極推動造林、生林等項目工程發(fā)展。②合理規(guī)劃開發(fā),提高利用效率。結(jié)合不同梯度帶區(qū)域的城鎮(zhèn)發(fā)展規(guī)模和類型,差異化地明確城鎮(zhèn)開發(fā)邊界;同時,根據(jù)土地利用空間分異規(guī)律特征,在國土空間規(guī)劃中按照土地類型比例綜合考慮并設(shè)置混合用地,打破單一土地單一功能的局限性,促進福建省土地節(jié)約集約利用。

    本文也存在一些不足,如僅分析了海陸梯度的分異規(guī)律,對于同一梯度帶內(nèi)的南北區(qū)域的發(fā)展與特征差異未進一步深入探討;對于海陸梯度分異的原因僅從定性角度討論,后續(xù)還需要結(jié)合定量角度進一步探討具體原因。

    猜你喜歡
    分異林地梯度
    一個改進的WYL型三項共軛梯度法
    一種自適應(yīng)Dai-Liao共軛梯度法
    一類扭積形式的梯度近Ricci孤立子
    丹東市林地分類研究
    重慶市臭氧時空分異及其影響因素研究
    淺談林地保護及恢復(fù)措施
    平泉縣下營坊雜巖體分異演化及其成巖成礦
    林地流轉(zhuǎn)模式的選擇機理及其政策啟示
    小型無人飛行器用于林地監(jiān)視的嘗試
    地溫梯度判定地熱異常的探討
    河南科技(2014年3期)2014-02-27 14:05:45
    av在线老鸭窝| 国产高清不卡午夜福利| 国产色婷婷99| or卡值多少钱| 超碰av人人做人人爽久久| 插逼视频在线观看| 尾随美女入室| 国模一区二区三区四区视频| 久久精品久久久久久噜噜老黄 | 99久国产av精品国产电影| 日本欧美国产在线视频| 最后的刺客免费高清国语| 日本五十路高清| 日本黄色视频三级网站网址| 欧美成人精品欧美一级黄| 99久久成人亚洲精品观看| 在线观看午夜福利视频| 最近视频中文字幕2019在线8| 国内揄拍国产精品人妻在线| 国产精品久久久久久精品电影| 久久久精品94久久精品| 久久久精品94久久精品| 日韩成人av中文字幕在线观看 | .国产精品久久| 国产精品野战在线观看| 最近的中文字幕免费完整| 国产成人影院久久av| h日本视频在线播放| 国产亚洲精品久久久久久毛片| 成人精品一区二区免费| 一区二区三区免费毛片| 天堂av国产一区二区熟女人妻| 日本欧美国产在线视频| 香蕉av资源在线| 国产成人精品久久久久久| 可以在线观看的亚洲视频| 久久精品国产清高在天天线| 国产精品人妻久久久久久| 亚洲欧美日韩高清专用| 亚洲18禁久久av| 精品午夜福利在线看| 91在线精品国自产拍蜜月| 女人十人毛片免费观看3o分钟| 美女cb高潮喷水在线观看| 午夜福利在线观看吧| 女生性感内裤真人,穿戴方法视频| 床上黄色一级片| 日韩欧美三级三区| 国产成年人精品一区二区| 深夜精品福利| 波多野结衣高清作品| 天堂动漫精品| 亚洲熟妇熟女久久| 国产精品三级大全| 日韩强制内射视频| 国产一区二区激情短视频| а√天堂www在线а√下载| 1024手机看黄色片| 精品久久久久久久人妻蜜臀av| 国产精品综合久久久久久久免费| 激情 狠狠 欧美| 国产视频内射| 成年版毛片免费区| 亚州av有码| 国内精品宾馆在线| 真人做人爱边吃奶动态| 2021天堂中文幕一二区在线观| 嫩草影院精品99| 非洲黑人性xxxx精品又粗又长| 麻豆一二三区av精品| 国产乱人视频| 日韩成人av中文字幕在线观看 | 日韩精品有码人妻一区| 国产成人影院久久av| 最近最新中文字幕大全电影3| 欧美区成人在线视频| 亚洲av.av天堂| 最近手机中文字幕大全| 日韩一区二区视频免费看| 成人毛片a级毛片在线播放| 国产av不卡久久| 久久久国产成人免费| 成人一区二区视频在线观看| 亚洲成人久久性| 亚洲内射少妇av| 最近2019中文字幕mv第一页| 国产探花极品一区二区| 久久久成人免费电影| 欧美xxxx黑人xx丫x性爽| 精品人妻一区二区三区麻豆 | 狠狠狠狠99中文字幕| 亚洲av电影不卡..在线观看| 日韩在线高清观看一区二区三区| 大香蕉久久网| 国产精品久久久久久亚洲av鲁大| 午夜影院日韩av| 搡老岳熟女国产| 欧美又色又爽又黄视频| 色在线成人网| 蜜臀久久99精品久久宅男| www.色视频.com| 亚洲av免费在线观看| 99热只有精品国产| 久久人人爽人人爽人人片va| 国产欧美日韩精品一区二区| 欧美高清成人免费视频www| av.在线天堂| 日本色播在线视频| 国产欧美日韩精品亚洲av| 一区二区三区免费毛片| 国产精品伦人一区二区| 不卡视频在线观看欧美| 日韩三级伦理在线观看| 欧美丝袜亚洲另类| 国产精品1区2区在线观看.| 国产精品久久久久久精品电影| 春色校园在线视频观看| 97超碰精品成人国产| 91av网一区二区| 女人十人毛片免费观看3o分钟| 亚洲第一电影网av| 大型黄色视频在线免费观看| 亚洲欧美日韩卡通动漫| 精品久久久噜噜| 欧美+亚洲+日韩+国产| 在线观看美女被高潮喷水网站| 亚洲欧美精品综合久久99| 免费无遮挡裸体视频| 国产精品一及| 麻豆精品久久久久久蜜桃| 天天一区二区日本电影三级| 亚洲精品色激情综合| 熟女人妻精品中文字幕| 欧美不卡视频在线免费观看| 亚洲国产欧洲综合997久久,| 欧美成人免费av一区二区三区| 免费无遮挡裸体视频| 在线免费观看的www视频| 久久久色成人| 欧美高清成人免费视频www| 美女免费视频网站| 在现免费观看毛片| 精品少妇黑人巨大在线播放 | 亚洲激情五月婷婷啪啪| 六月丁香七月| 午夜免费男女啪啪视频观看 | av在线老鸭窝| 蜜臀久久99精品久久宅男| 99国产精品一区二区蜜桃av| 亚洲三级黄色毛片| 黄片wwwwww| 俺也久久电影网| 日韩一区二区视频免费看| 午夜福利18| 亚洲五月天丁香| 人妻丰满熟妇av一区二区三区| 97碰自拍视频| 久久午夜福利片| av专区在线播放| 亚洲美女黄片视频| 亚洲va在线va天堂va国产| 免费观看精品视频网站| 亚洲精品一区av在线观看| 女人十人毛片免费观看3o分钟| 真实男女啪啪啪动态图| 日韩欧美精品v在线| 18禁在线无遮挡免费观看视频 | 欧美日韩乱码在线| 亚洲三级黄色毛片| 欧美日韩在线观看h| 久久午夜福利片| 国产欧美日韩精品一区二区| 日韩欧美一区二区三区在线观看| av女优亚洲男人天堂| 国产成人一区二区在线| 免费高清视频大片| 夜夜夜夜夜久久久久| h日本视频在线播放| 长腿黑丝高跟| 欧美另类亚洲清纯唯美| 91久久精品电影网| av卡一久久| 青春草视频在线免费观看| 日韩精品青青久久久久久| 成人午夜高清在线视频| 看片在线看免费视频| 亚洲精品日韩av片在线观看| 可以在线观看的亚洲视频| 国产精品乱码一区二三区的特点| 久久久久久国产a免费观看| 人妻久久中文字幕网| 精品一区二区三区人妻视频| 性插视频无遮挡在线免费观看| 亚洲丝袜综合中文字幕| 亚洲av.av天堂| 韩国av在线不卡| 亚洲专区国产一区二区| 我要看日韩黄色一级片| 国内久久婷婷六月综合欲色啪| 日本免费一区二区三区高清不卡| 久久午夜福利片| 久久久欧美国产精品| 99久久精品热视频| 性欧美人与动物交配| 日本黄色视频三级网站网址| 久久6这里有精品| 国产成人a区在线观看| 国产男靠女视频免费网站| 一级毛片我不卡| 91在线观看av| 国产av不卡久久| 人妻久久中文字幕网| 亚洲成人久久性| 久久6这里有精品| 一级毛片aaaaaa免费看小| 精品人妻视频免费看| 给我免费播放毛片高清在线观看| 成年版毛片免费区| 免费观看人在逋| 国产在线精品亚洲第一网站| 国产不卡一卡二| 亚洲欧美中文字幕日韩二区| 九九在线视频观看精品| 国产精品女同一区二区软件| 淫妇啪啪啪对白视频| 亚洲欧美日韩无卡精品| 老司机影院成人| 久久久久国产网址| 99九九线精品视频在线观看视频| 亚洲欧美日韩高清在线视频| 久久99热这里只有精品18| 久久久久性生活片| 午夜福利高清视频| 男人的好看免费观看在线视频| 日韩成人伦理影院| 女生性感内裤真人,穿戴方法视频| 精品99又大又爽又粗少妇毛片| 久久热精品热| 老司机影院成人| 99久久九九国产精品国产免费| 精品久久国产蜜桃| 男女那种视频在线观看| 日韩精品有码人妻一区| 成人鲁丝片一二三区免费| 国产成人精品久久久久久| 成人午夜高清在线视频| 国产伦一二天堂av在线观看| 丰满乱子伦码专区| 成人综合一区亚洲| 久久欧美精品欧美久久欧美| av免费在线看不卡| 亚洲av五月六月丁香网| 18禁黄网站禁片免费观看直播| 嫩草影院新地址| 国产精品一区二区三区四区久久| 亚洲av第一区精品v没综合| 又爽又黄无遮挡网站| 中文资源天堂在线| 日产精品乱码卡一卡2卡三| 欧美性猛交╳xxx乱大交人| 国产亚洲av嫩草精品影院| 欧美高清性xxxxhd video| 免费看光身美女| 国产男靠女视频免费网站| 哪里可以看免费的av片| 插阴视频在线观看视频| 久久精品国产鲁丝片午夜精品| 亚洲成人av在线免费| 丰满的人妻完整版| 亚洲最大成人手机在线| 免费av不卡在线播放| 少妇裸体淫交视频免费看高清| 久久久精品大字幕| 亚州av有码| 亚洲美女搞黄在线观看 | 国产淫片久久久久久久久| .国产精品久久| 观看美女的网站| 看片在线看免费视频| 久久中文看片网| 国产69精品久久久久777片| 国产女主播在线喷水免费视频网站 | 直男gayav资源| 色在线成人网| 久久久久性生活片| 日日啪夜夜撸| 久久人人爽人人片av| 亚洲一区二区三区色噜噜| 一区二区三区免费毛片| 亚洲av二区三区四区| 我要搜黄色片| 在线国产一区二区在线| 国产一区二区在线av高清观看| 级片在线观看| 99热网站在线观看| 色尼玛亚洲综合影院| 美女大奶头视频| 午夜激情福利司机影院| 亚洲中文字幕一区二区三区有码在线看| 性插视频无遮挡在线免费观看| 欧美绝顶高潮抽搐喷水| 亚洲av一区综合| 赤兔流量卡办理| 精品午夜福利在线看| 99久久成人亚洲精品观看| 国产一区二区在线观看日韩| 欧美3d第一页| av在线观看视频网站免费| 人妻少妇偷人精品九色| 日本-黄色视频高清免费观看| 亚洲五月天丁香| 黄色配什么色好看| 国产高清有码在线观看视频| 天天一区二区日本电影三级| а√天堂www在线а√下载| 日韩成人伦理影院| 综合色av麻豆| 免费av观看视频| 亚洲精品粉嫩美女一区| 99精品在免费线老司机午夜| 国产成人福利小说| 真人做人爱边吃奶动态| 亚洲精品国产成人久久av| 啦啦啦韩国在线观看视频| 国产精品免费一区二区三区在线| 乱码一卡2卡4卡精品| 99久久精品热视频| 国产在视频线在精品| 久久精品久久久久久噜噜老黄 | 秋霞在线观看毛片| 国产精品久久久久久亚洲av鲁大| 国产美女午夜福利| 亚洲av中文av极速乱| 在线免费十八禁| 国产极品精品免费视频能看的| 99riav亚洲国产免费| 婷婷精品国产亚洲av| 丝袜美腿在线中文| 99久久九九国产精品国产免费| 最近中文字幕高清免费大全6| av在线亚洲专区| 最近视频中文字幕2019在线8| 少妇的逼水好多| 国产成人福利小说| 人人妻人人看人人澡| 亚洲国产欧美人成| 小蜜桃在线观看免费完整版高清| 日韩强制内射视频| 成人永久免费在线观看视频| 日本a在线网址| 亚洲,欧美,日韩| 中国国产av一级| 波多野结衣巨乳人妻| 男女边吃奶边做爰视频| 看黄色毛片网站| 淫秽高清视频在线观看| 亚洲国产日韩欧美精品在线观看| 深爱激情五月婷婷| 内射极品少妇av片p| 超碰av人人做人人爽久久| 日日摸夜夜添夜夜添av毛片| 日韩欧美免费精品| 亚洲欧美日韩东京热| 日韩欧美免费精品| 日韩精品中文字幕看吧| 丝袜美腿在线中文| 在线a可以看的网站| 欧美xxxx黑人xx丫x性爽| 三级男女做爰猛烈吃奶摸视频| 欧美激情在线99| 日韩精品中文字幕看吧| 美女内射精品一级片tv| 国产熟女欧美一区二区| 好男人在线观看高清免费视频| 永久网站在线| 国产人妻一区二区三区在| 一区二区三区高清视频在线| 99热6这里只有精品| 国产三级中文精品| 最新在线观看一区二区三区| 婷婷色综合大香蕉| 日本在线视频免费播放| 91麻豆精品激情在线观看国产| 国产大屁股一区二区在线视频| 免费av毛片视频| 老女人水多毛片| 黑人高潮一二区| 欧美一级a爱片免费观看看| 国产一区亚洲一区在线观看| 国产探花在线观看一区二区| 91久久精品国产一区二区三区| 97超碰精品成人国产| av在线天堂中文字幕| 国产成人一区二区在线| 丝袜美腿在线中文| 国产亚洲精品久久久久久毛片| 免费观看的影片在线观看| 亚洲电影在线观看av| 日本欧美国产在线视频| 欧美+亚洲+日韩+国产| 国产日本99.免费观看| 午夜福利18| 晚上一个人看的免费电影| 男女下面进入的视频免费午夜| 久久精品国产亚洲av天美| 日日撸夜夜添| 精品久久久久久久久久免费视频| 色综合站精品国产| 精品少妇黑人巨大在线播放 | 美女xxoo啪啪120秒动态图| 国产伦在线观看视频一区| 校园人妻丝袜中文字幕| av福利片在线观看| 久久九九热精品免费| 成人av一区二区三区在线看| 国产中年淑女户外野战色| 日本五十路高清| 亚洲精品456在线播放app| 欧美bdsm另类| 丝袜喷水一区| a级毛片免费高清观看在线播放| 久久亚洲国产成人精品v| 嫩草影院新地址| 国语自产精品视频在线第100页| 中国美女看黄片| 亚洲电影在线观看av| av女优亚洲男人天堂| 久久久午夜欧美精品| 欧美中文日本在线观看视频| 男女做爰动态图高潮gif福利片| 国产大屁股一区二区在线视频| 天堂动漫精品| 亚洲国产日韩欧美精品在线观看| 日韩av在线大香蕉| 国产精品久久久久久久电影| 亚洲人与动物交配视频| 美女免费视频网站| 最近最新中文字幕大全电影3| 免费搜索国产男女视频| 国产极品精品免费视频能看的| 亚洲成人久久爱视频| 免费看a级黄色片| 欧美极品一区二区三区四区| 97超视频在线观看视频| 赤兔流量卡办理| 天美传媒精品一区二区| 变态另类成人亚洲欧美熟女| 亚洲欧美中文字幕日韩二区| 国产一区二区三区av在线 | 在线观看免费视频日本深夜| 成人av在线播放网站| 欧美日本亚洲视频在线播放| 别揉我奶头 嗯啊视频| 天堂动漫精品| 亚洲性久久影院| 天天躁日日操中文字幕| 国产视频一区二区在线看| 国产淫片久久久久久久久| 自拍偷自拍亚洲精品老妇| 亚洲欧美精品综合久久99| 香蕉av资源在线| 欧美潮喷喷水| 少妇裸体淫交视频免费看高清| 婷婷精品国产亚洲av| 精品欧美国产一区二区三| 97碰自拍视频| 国产精品久久电影中文字幕| 精品午夜福利在线看| 亚洲精华国产精华液的使用体验 | 内地一区二区视频在线| 五月玫瑰六月丁香| 在线播放无遮挡| 少妇猛男粗大的猛烈进出视频 | 欧美性猛交黑人性爽| 久久精品国产亚洲av天美| 国产成人精品久久久久久| aaaaa片日本免费| 自拍偷自拍亚洲精品老妇| 亚洲精品乱码久久久v下载方式| 波多野结衣高清作品| 中国美白少妇内射xxxbb| 欧美成人a在线观看| 免费在线观看影片大全网站| 男人狂女人下面高潮的视频| 美女黄网站色视频| 91在线观看av| 中国美女看黄片| 22中文网久久字幕| 欧美成人精品欧美一级黄| 亚洲一级一片aⅴ在线观看| 看黄色毛片网站| 69人妻影院| 成年女人看的毛片在线观看| 少妇人妻精品综合一区二区 | 看黄色毛片网站| 亚洲aⅴ乱码一区二区在线播放| 欧美日韩一区二区视频在线观看视频在线 | 中国美女看黄片| 天天躁夜夜躁狠狠久久av| 老司机福利观看| 丰满人妻一区二区三区视频av| 日韩欧美在线乱码| 午夜福利在线观看吧| 日韩人妻高清精品专区| 蜜桃久久精品国产亚洲av| 国产亚洲91精品色在线| 国产精品乱码一区二三区的特点| 午夜精品在线福利| 亚洲成人中文字幕在线播放| 精品久久久久久久久久免费视频| 在线观看美女被高潮喷水网站| 亚洲人成网站在线播放欧美日韩| 在线看三级毛片| 一进一出好大好爽视频| 免费av观看视频| 欧美成人一区二区免费高清观看| 久久久午夜欧美精品| 国内久久婷婷六月综合欲色啪| 一级毛片电影观看 | 国产在视频线在精品| 淫秽高清视频在线观看| 成人亚洲欧美一区二区av| 精品无人区乱码1区二区| 搡女人真爽免费视频火全软件 | 亚洲18禁久久av| 国产精品99久久久久久久久| 亚洲最大成人手机在线| 成人三级黄色视频| 亚洲欧美日韩高清在线视频| 亚洲七黄色美女视频| 不卡一级毛片| 精品日产1卡2卡| 亚洲内射少妇av| 成年版毛片免费区| 中文字幕精品亚洲无线码一区| 日韩制服骚丝袜av| 久久亚洲国产成人精品v| 午夜免费激情av| 亚洲乱码一区二区免费版| 国产色爽女视频免费观看| 我要搜黄色片| 人妻夜夜爽99麻豆av| 午夜福利成人在线免费观看| 国产精品人妻久久久影院| 国产午夜精品论理片| 国产亚洲av嫩草精品影院| 97人妻精品一区二区三区麻豆| 亚洲欧美清纯卡通| 22中文网久久字幕| 午夜精品国产一区二区电影 | 草草在线视频免费看| 可以在线观看毛片的网站| av视频在线观看入口| ponron亚洲| 精品99又大又爽又粗少妇毛片| 色哟哟·www| 人人妻人人澡人人爽人人夜夜 | 寂寞人妻少妇视频99o| 国产精品久久久久久久久免| 99在线人妻在线中文字幕| 色综合色国产| 非洲黑人性xxxx精品又粗又长| 一进一出好大好爽视频| 国产伦一二天堂av在线观看| 精品欧美国产一区二区三| 一级a爱片免费观看的视频| 免费观看精品视频网站| 男女视频在线观看网站免费| 人人妻,人人澡人人爽秒播| 真实男女啪啪啪动态图| 狠狠狠狠99中文字幕| 国内少妇人妻偷人精品xxx网站| 久久精品国产亚洲av天美| 激情 狠狠 欧美| 99视频精品全部免费 在线| 国产亚洲欧美98| 亚洲av成人av| 精品福利观看| 99热这里只有是精品在线观看| 国产一区亚洲一区在线观看| 日韩,欧美,国产一区二区三区 | 久久午夜亚洲精品久久| 久久精品久久久久久噜噜老黄 | 欧美潮喷喷水| 午夜福利在线观看吧| 亚洲成人精品中文字幕电影| 日韩大尺度精品在线看网址| 国产成人a∨麻豆精品| 亚洲精品一区av在线观看| 欧美一级a爱片免费观看看| 亚洲国产精品成人综合色| 99国产极品粉嫩在线观看| 成人美女网站在线观看视频| 免费搜索国产男女视频| 人人妻人人澡欧美一区二区| 99热这里只有是精品50| 99国产精品一区二区蜜桃av| 日日啪夜夜撸| 国内精品久久久久精免费| 免费搜索国产男女视频| av在线观看视频网站免费| 中文字幕久久专区| videossex国产| 美女免费视频网站| 久久久久久久久中文| 午夜精品在线福利| 国产伦精品一区二区三区视频9| 国产精品不卡视频一区二区| 一本一本综合久久| 99久久精品热视频| 国产免费男女视频| 三级毛片av免费| 国产高清视频在线观看网站| 欧美不卡视频在线免费观看| 亚洲成a人片在线一区二区|