沿海資源富集區(qū)“三生空間”演變與生態(tài)環(huán)境響應(yīng)

陳寶強(qiáng)， 楊永崇， 楊梅煥

(西安科技大學(xué)測繪科學(xué)與技術(shù)學(xué)院， 西安 710054)

為構(gòu)建生產(chǎn)空間-生活空間-生態(tài)空間相互協(xié)調(diào)的國土空間布局，十八大強(qiáng)調(diào)“調(diào)整空間結(jié)構(gòu)，推進(jìn)生產(chǎn)集約高效、生活宜居適度、生態(tài)山清水秀”，十九大指出要合理確定“三區(qū)三線”，是參照土地利用功能對土地所涵蓋功能分類成生產(chǎn)空間、生活空間和生態(tài)空間(“三生空間”)。“三生空間”區(qū)域性研究是協(xié)調(diào)自然資源保護(hù)與土地合理利用的基礎(chǔ)性工作，可以為國土空間規(guī)劃與生態(tài)紅線提供借鑒[1-2]。隨著城鎮(zhèn)化外延式、粗放式[3-4]的模式發(fā)展，會推動土地利用功能發(fā)生轉(zhuǎn)型[5]，隨之會影響區(qū)域的生態(tài)環(huán)境。土地利用轉(zhuǎn)型[6-7](land-use and land-cover change, LUCC)始終是中外關(guān)注的重點內(nèi)容之一，當(dāng)前諸多研究的是土地利用變化與生態(tài)系統(tǒng)服務(wù)價值變化[8-10]，同時針對“三生空間”土地利用功能演變與生態(tài)環(huán)境關(guān)系大多數(shù)以省市、區(qū)域性展開研究，沿海低丘陵資源富集區(qū)研究尺度綜合性較少。

山東省招遠(yuǎn)市隸屬于膠東半島沿海低丘陵資源富集區(qū)西北部，“黃金之都”，礦產(chǎn)資源最豐富，儲量和產(chǎn)量占全國的1/10，其中花崗巖礦是第二大優(yōu)勢礦產(chǎn)資源。21世紀(jì)以來，膠東半島沿海低丘陵區(qū)域步入快速發(fā)展階段，區(qū)域性土地利用結(jié)構(gòu)隨礦產(chǎn)資源開發(fā)、城鎮(zhèn)擴(kuò)張等發(fā)生明顯變化，導(dǎo)致區(qū)域性土地利用功能也發(fā)生變化，對沿海低丘陵資源富集區(qū)生態(tài)環(huán)境質(zhì)量的研究是在自然與人類活動的長期相互作用過程演變形成的，對沿海低丘陵資源富集區(qū)生態(tài)環(huán)境質(zhì)量探討是為揭示土地利用功能轉(zhuǎn)型、生態(tài)環(huán)境質(zhì)量演變規(guī)律，明確土地利用轉(zhuǎn)型的生態(tài)環(huán)境響應(yīng)，為資源富集區(qū)土地資源整治和生態(tài)環(huán)境保護(hù)提供借鑒。因此，探討研究膠東半島沿海低丘陵資源富集區(qū)招遠(yuǎn)市土地利用轉(zhuǎn)型的生態(tài)環(huán)境響應(yīng)，基于“三生空間”視角，選取沿海低丘陵資源富集區(qū)招遠(yuǎn)市為對象，定量分析研究區(qū)域2008年、2013年、2018年土地利用轉(zhuǎn)型及生態(tài)環(huán)境響應(yīng)。

1 數(shù)據(jù)來源與研究方法

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

研究區(qū)以2008年、2013年和2018年招遠(yuǎn)市土地利用變更調(diào)查矢量數(shù)據(jù)為基礎(chǔ)數(shù)據(jù)，嚴(yán)格按照《土地利用現(xiàn)狀分類標(biāo)準(zhǔn)(GB/T 21010—2007)》分類，以“三生空間”土地利用主導(dǎo)功能進(jìn)行整理(圖1)，參考諸多學(xué)者的研究成果[11]，采用生態(tài)環(huán)境質(zhì)量指數(shù)面積加權(quán)賦值的方法,如表1所示。

圖1 沿海低丘陵區(qū)“三生空間”土地利用現(xiàn)狀空間分布Fig.1 The land use status of Ecologicl-Production-Living Spaces in coastal low hilly areas distributed spatially

表1 沿海低丘陵區(qū)域“三生空間”土地利用主導(dǎo)功能分類與生態(tài)環(huán)境質(zhì)量指數(shù)表Table 1 Classification of the leading functions of land use in Ecologicl-Production-Living Spaces in low hilly coastal areas and the index table of ecological environment quality

1.2 土地利用空間擴(kuò)展模型

土地利用空間擴(kuò)展模型[12-13]主要分析“三生空間”用地的空間擴(kuò)展與收縮。

(1)

式(1)中：UII為土地利用空間擴(kuò)展；ULAi,t+m為t+m年空間單元i土地利用面積；ULAi,t分別為t年空間單元i土地利用面積；TLAi為研究區(qū)域土地利用總面積。

1.3 土地利用轉(zhuǎn)移模型

土地利用轉(zhuǎn)移模型[14]主要用來分析不同時期土地利用主導(dǎo)功能之間轉(zhuǎn)移情況。

(2)

式(2)中：Sij為由i轉(zhuǎn)換成j類型土地利用主導(dǎo)功能的面積；n為某“三生空間”類型數(shù)量。

1.4 土地利用功能重心模型

土地利用功能重心模型[15]反映“三生空間”功能用地在某時間內(nèi)在空間重心的變化。

(3)

式(3)中：Xt、Yt分別為t時間段某土地利用功能用地地理坐標(biāo)；Xh、Yh分別為某k個單元的重心幾何坐標(biāo)；Shi為hi區(qū)域“三生空間”土地利用主導(dǎo)功能用地的面積；D為相鄰年份之間土地利用功能用地偏移距離。

1.5 “三生空間”時空轉(zhuǎn)型生態(tài)環(huán)境響應(yīng)模型

1.5.1 生態(tài)單元綜合環(huán)境質(zhì)量指數(shù)

生態(tài)單元綜合環(huán)境質(zhì)量指數(shù)[16-18]是“三生空間”用地生態(tài)環(huán)境質(zhì)量與面積的比值。

(4)

式(4)中：Eh為h區(qū)域所代表的生態(tài)環(huán)境質(zhì)量指數(shù)；Sh為h研究區(qū)域“三生空間”土地利用主導(dǎo)功能的面積之和；eh為第h類“三生空間”類型用地的生態(tài)環(huán)境指數(shù)。

1.5.2 半變異函數(shù)分析模型

研究尺度不同，結(jié)論不同。現(xiàn)采用等間距200 m×200 m格網(wǎng)作為空間單元，將每個生態(tài)環(huán)境質(zhì)量指數(shù)賦值給空間單元，利用半變異函數(shù)分析模型、克里金空間，形成研究區(qū)整個生態(tài)環(huán)境質(zhì)量等級空間分布，半變異函數(shù)[19]公式為

(5)

式(5)中：r(h)為半方差；h為樣本的間距；M(xi)為xi處的生態(tài)環(huán)境質(zhì)量指數(shù)；M(xi+h)為xi+h處的生態(tài)環(huán)境質(zhì)量指數(shù)。采用自然間斷法分為5級；低質(zhì)量區(qū)(E≤0.2)、較低質(zhì)量區(qū)(0.2

1.5.3 生態(tài)貢獻(xiàn)率模型

生態(tài)貢獻(xiàn)率[20-25]是指研究區(qū)域土地利用功能轉(zhuǎn)型對區(qū)域環(huán)境生態(tài)質(zhì)量改變的影響。

(6)

式(6)中：L為生態(tài)貢獻(xiàn)度；LEi為初期某“三生空間”土地主導(dǎo)功能用地的生態(tài)環(huán)境質(zhì)量指數(shù)；LEj為末期某“三生空間”土地主導(dǎo)功能用地的生態(tài)環(huán)境質(zhì)量指數(shù)；ΔSh為“三生空間”土地利用主導(dǎo)功能變化初期面積與末期面積差。

2 分析與討論

2.1 土地利用轉(zhuǎn)型空間總體分析

依據(jù)式(1)計算2008—2018年土地利用主導(dǎo)功能用地的空間擴(kuò)展。農(nóng)業(yè)生產(chǎn)用地、其他生態(tài)用地、林地、水域生態(tài)用地擴(kuò)展強(qiáng)度表現(xiàn)為負(fù)值(表2)，表現(xiàn)為空間收縮；城鎮(zhèn)用地、農(nóng)村用地和工礦生產(chǎn)用地表現(xiàn)為空間擴(kuò)張，服務(wù)及其他生產(chǎn)用地變化較小。沿海低丘陵資源富集區(qū)空間擴(kuò)展與空間收縮強(qiáng)度不高，處于一種低速變化的狀態(tài)。其中工礦生產(chǎn)用地的空間擴(kuò)展強(qiáng)度最高，2008—2013年為0.01%，2013—2018年為0.06%。

表2 沿海低丘陵區(qū)域土地利用空間擴(kuò)展強(qiáng)度Table 2 Spatial expansion intensity of land use in coastal low hilly areas

2.2 土地利用轉(zhuǎn)型時空演變特征分析

為了研究分析沿海低丘陵資源富集區(qū)的土地利用功能演變過程與特征，根據(jù)式(2)，計算出土地利用轉(zhuǎn)移矩陣(表3、表4)。21世紀(jì)以來，膠東半島沿海低丘陵步入快速發(fā)展軌道，區(qū)域性土地利用結(jié)構(gòu)隨礦產(chǎn)資源開發(fā)、城鎮(zhèn)擴(kuò)張等發(fā)生明顯變化，城鎮(zhèn)用地面積增多295.92 hm2，變化10.2%，工礦生產(chǎn)用地面積增加495.68 hm2，變化17.4%，農(nóng)村用地面積增多1 309.07 hm2，變化15.31%，服務(wù)及其他生產(chǎn)用地變化比較平緩。

2008—2013年，林地生態(tài)用地、農(nóng)業(yè)生產(chǎn)用地的轉(zhuǎn)入是城鎮(zhèn)用地增多的主導(dǎo)因素，分別轉(zhuǎn)入33.64 hm2、170.07 hm2，分別變化1.08%、5.5%；隨工礦生產(chǎn)用地修復(fù)與復(fù)墾，大量工礦用地轉(zhuǎn)型為農(nóng)業(yè)生產(chǎn)用地，轉(zhuǎn)移面積144.56 hm2，轉(zhuǎn)移率為4.9%；由于經(jīng)濟(jì)高速發(fā)展，人口越來越密集，大量的林地、耕地等生態(tài)用地被建設(shè)用地占用，主要表現(xiàn)為林地生態(tài)用地、農(nóng)業(yè)生產(chǎn)用地、其他生態(tài)用地占主導(dǎo)因素轉(zhuǎn)型為農(nóng)村用地，轉(zhuǎn)入面積分別為165.09 hm2、636.63 hm2、41.96 hm2，轉(zhuǎn)入率分別為1.8%、6.8%、0.4%。2013—2018年，城鎮(zhèn)用地轉(zhuǎn)入主要類型為農(nóng)業(yè)生產(chǎn)用地，較上一時段相比，轉(zhuǎn)入率由5.5%降低至2.13%。工礦生產(chǎn)用地的增加與林地生態(tài)用地、農(nóng)業(yè)生產(chǎn)用地的增加存在互逆情況，分別轉(zhuǎn)入35.62 hm2、333.15 hm2，分別轉(zhuǎn)入率為1.07%、10.0%。較上一時間段相比，農(nóng)村用地增多主要是林地生態(tài)用地、農(nóng)業(yè)生產(chǎn)用地、其他生態(tài)用地轉(zhuǎn)入，轉(zhuǎn)入幅度較小，分別為72.97 hm2，390.74 hm2，21.23 hm2。

2008—2018年，沿海低丘陵資源富集區(qū)林地、耕地、未利用地隨時間的推移呈降低趨勢，農(nóng)業(yè)生產(chǎn)用地最明顯，減少1 325.92 hm2。2008—2013年，林地生態(tài)用地降低主要是農(nóng)村用地、農(nóng)業(yè)生產(chǎn)用地增多，分別轉(zhuǎn)出165.09 hm2、93.02 hm2，轉(zhuǎn)出率分別為0.60%、0.34%。農(nóng)業(yè)生產(chǎn)用地降低與城鎮(zhèn)用地、工礦生產(chǎn)用地、農(nóng)村用地的增多有關(guān)，轉(zhuǎn)出0.19%、0.16%、0.73%。因裸土地、沙地等生態(tài)用地受人類活動干擾影響，主要表現(xiàn)為其他生態(tài)用地隨時間推移轉(zhuǎn)型為農(nóng)村用地、農(nóng)村生產(chǎn)用地，轉(zhuǎn)移面積分別為41.96 hm2、181.51 hm2，轉(zhuǎn)出率分別為0.64%、2.77%。水域生態(tài)用地減少與工礦生態(tài)用地、農(nóng)村用地增多有關(guān)，分別轉(zhuǎn)出16.15 hm2、17.90 hm2。2013—2018年，農(nóng)業(yè)生產(chǎn)用地轉(zhuǎn)出主要為工礦生產(chǎn)用地、農(nóng)村用地，轉(zhuǎn)出率為0.38%、0.45%。較上一時間段相比，林地生態(tài)用地轉(zhuǎn)出主要為農(nóng)村用地、農(nóng)業(yè)生產(chǎn)用地，分別為轉(zhuǎn)移面積72.97 hm2、79.35 hm2，其他生態(tài)用地轉(zhuǎn)化為農(nóng)村用地、農(nóng)業(yè)生產(chǎn)用地，整體幅度降低，轉(zhuǎn)出面積分別為21.23 hm2、73.65 hm2；水域生態(tài)用地減少與工礦生產(chǎn)用地、農(nóng)業(yè)生產(chǎn)用地的增多有關(guān)，轉(zhuǎn)移面積14.22 hm2、14.19 hm2。整體來說，生活空間與生產(chǎn)空間功能增強(qiáng)，生態(tài)空間減弱是研究“三生空間”的轉(zhuǎn)型的主基調(diào)，其轉(zhuǎn)型由礦產(chǎn)資源的開發(fā)、土地管理政策等多因素造成。

表3 沿海低丘陵區(qū)域2008—2013年土地利用轉(zhuǎn)移矩陣表Table 3 Land use transfer matrix from 2008 to 2013 in coastal low hilly areas

表4 沿海低丘陵區(qū)域2013—2018年土地利用轉(zhuǎn)移矩陣表Table 4 Land use transfer matrix of coastal low hilly region from 2013 to 2018

2.3 土地利用空間轉(zhuǎn)型重心遷移分析

利用式(3)整理出2008—2018“三生空間”土地利用主導(dǎo)功能用地的重心位置與轉(zhuǎn)移方向，如圖2～圖9所示，表示不同時期“三生空間”各功能用地重心遷移情況。

從各功能用地的重心位置來分析，2008—2018年，農(nóng)業(yè)生產(chǎn)用地、水域生態(tài)用地、林地生態(tài)用地重心均位于羅峰街道內(nèi)，主要表現(xiàn)為羅峰街道內(nèi)具有豐富的水資源，龍王水庫與界河橫穿整個研究區(qū)域，西高東低，南北平坦，適合農(nóng)作物種植。

工礦生產(chǎn)用地、其他生態(tài)用地重心均位于夢芝街道內(nèi)，主要表現(xiàn)為夢芝街道是沿海開放地區(qū)之一，水陸空四通八達(dá)，城市規(guī)劃合理，具有完善的基礎(chǔ)建設(shè)。

圖2 農(nóng)業(yè)生產(chǎn)用地重心遷移圖Fig.2 Gravity center transfer map of agricultural production land

圖3 水域生態(tài)用地重心遷移圖Fig.3 Gravity center transfer map of water ecological land

服務(wù)及其他用地重心由夏甸鎮(zhèn)轉(zhuǎn)向齊山鎮(zhèn)內(nèi)，由于齊山鎮(zhèn)地處南北水系大沽河與界河水系的分水嶺地帶，北高南低，丘陵起伏，有著名“十字道”歷史意義。

城鎮(zhèn)用地、農(nóng)村用地重心位于羅峰街道幾何中心的東北方向，朝溫泉街道經(jīng)濟(jì)新區(qū)與老城區(qū)靠攏，區(qū)位條件優(yōu)越、交通服務(wù)便利。

圖4 林地生態(tài)用地重心遷移圖Fig.4 Gravity center transfer map of woodland ecological land

圖5 工礦生產(chǎn)用地重心遷移圖Fig.5 Gravity center transfer map of industrial and mining production land

圖6 其他生態(tài)用地重心遷移圖Fig.6 Gravity center transfer map of other ecological land

圖7 服務(wù)及其他生產(chǎn)用地重心遷移圖Fig.7 Gravity center transfer map of service and other production land

圖8 城鎮(zhèn)用地重心遷移圖Fig.8 Gravity center transfer map of urban land

圖9 農(nóng)村用地重心遷移Fig.9 Gravity center transfer map of rural land

從各功能“三生空間”土地利用主導(dǎo)功能的轉(zhuǎn)移方向與偏移距離來分析，服務(wù)及其他生產(chǎn)用地、工礦生產(chǎn)用地、農(nóng)村用地、城鎮(zhèn)用地偏移最顯著，主要表現(xiàn)為城鎮(zhèn)粗放式與外延式擴(kuò)張[16-17]、資源開發(fā)等，使土地利用變化活躍，而水域生態(tài)用地、其他生態(tài)用地等受到外界環(huán)境干擾不顯著，短期內(nèi)發(fā)生劇烈變化的可能性較小。其中，2008—2013年，水域生態(tài)用地的重心偏移距離較小，整體呈現(xiàn)向西偏南方向轉(zhuǎn)移趨勢；其他生態(tài)用地重心先朝北偏西方向，偏移距離106 m，然后再朝西偏南方向轉(zhuǎn)移，偏移距離39 m；農(nóng)業(yè)生產(chǎn)用地重心先朝南偏東轉(zhuǎn)移，偏移距離112 m，在朝西偏南方向偏移，偏移距離為60 m；朝北偏西方向遷移是農(nóng)村用地、城鎮(zhèn)用地主要趨勢，分別偏移距離790 m、130 m；林地生態(tài)用地總體偏移距離66 m,整體呈現(xiàn)出先向東再向北偏東移動。工礦生產(chǎn)用地朝北偏西方向，偏移距離314 m，再朝南偏東方向，偏移距離958 m；服務(wù)及其他生產(chǎn)用地整體呈現(xiàn)往北偏東方向移動的趨勢，偏移距離2 630 m。

2.4 生態(tài)環(huán)境綜合質(zhì)量時序變化

依據(jù)式(4)、式(5)，計算沿海低丘陵2008年、2013年、2018年生態(tài)環(huán)境綜合質(zhì)量指數(shù)表(表5)與區(qū)域生態(tài)質(zhì)量等級表(表6)。

2008—2018年，沿海低丘陵生態(tài)質(zhì)量綜合指數(shù)整體水平平穩(wěn)，分別為0.426 6、0.423 8、0.420 7。由于城鎮(zhèn)外延外擴(kuò)張、經(jīng)濟(jì)新區(qū)的建設(shè)、工礦生產(chǎn)用地等影響，低質(zhì)量區(qū)、較低質(zhì)量區(qū)面積在增加，分別由729.99 hm2增加至956.61 hm2、48 324.6 hm2增加至50 891.22 hm2，中質(zhì)量區(qū)、較高質(zhì)量、高質(zhì)量區(qū)面積隨著時間的變化在縮減,分別降低2.9%、2.8%、3.8%。從空間上來看，較高質(zhì)量區(qū)、高質(zhì)量區(qū)主要分布在阜山鎮(zhèn)、玲瓏鎮(zhèn)，主要表現(xiàn)為阜山鎮(zhèn)、玲瓏鎮(zhèn)有覆蓋率高達(dá)90%以上的森林，對生態(tài)環(huán)境具有保護(hù)作用；低質(zhì)量區(qū)、較低質(zhì)量區(qū)主要分布在人口集聚區(qū)羅峰街道、溫泉街道、夢芝街道，中質(zhì)量區(qū)呈零星分布。中質(zhì)量區(qū)和較低質(zhì)量區(qū)構(gòu)成了沿海低丘陵土地利用生態(tài)環(huán)境質(zhì)量主體，高質(zhì)量區(qū)、較高質(zhì)量區(qū)轉(zhuǎn)型為低質(zhì)量區(qū)、較低質(zhì)量區(qū)，促進(jìn)了膠東半島沿海低丘陵資源富集區(qū)生態(tài)環(huán)境質(zhì)量平穩(wěn)發(fā)展。

表5 研究區(qū)域各年份生態(tài)環(huán)境綜合質(zhì)量指數(shù)Table 5 Ecological environment comprehensive quality index of each year in the study area

2.5 影響生態(tài)環(huán)境質(zhì)量的因素

區(qū)域生態(tài)環(huán)境質(zhì)量整體水平保持平穩(wěn)不等同于區(qū)域生態(tài)環(huán)境質(zhì)量沒有發(fā)生任何變化，如表5、表6所示，膠東半島低丘陵沿海區(qū)域生態(tài)環(huán)境質(zhì)量隨著時間的變化呈下降趨勢。如圖10所示，可能的原因如下。

(1)研究區(qū)域礦產(chǎn)資源富集區(qū)集中分布在蠶莊鎮(zhèn)、玲瓏鎮(zhèn)、張星鎮(zhèn)，其他呈零星分布，礦產(chǎn)資源屬于開發(fā)階段，短時間內(nèi)不會對環(huán)境產(chǎn)生影響，但是隨時間推移，會改變周圍土地利用主導(dǎo)功能。

表6 研究區(qū)域各年份生態(tài)環(huán)境質(zhì)量等級Table 6 Ecological environment quality grade of each year in the study area

(2)生態(tài)環(huán)境質(zhì)量降低占主導(dǎo)地位，夢芝街道、羅峰街道與溫泉街道，水陸空交通極為便利，同時是沿海開放地區(qū)之一，招商引資，發(fā)展輕工業(yè)，推動城鄉(xiāng)一體化快速發(fā)展，使生態(tài)環(huán)境質(zhì)量指數(shù)降低；張星鎮(zhèn)的粉絲、石材、機(jī)械發(fā)展迅速，導(dǎo)致城鎮(zhèn)用地不斷向外擴(kuò)張，使得區(qū)域生態(tài)環(huán)境質(zhì)量指數(shù)降低。因此，沿海低丘陵資源富集區(qū)隨時間變化區(qū)域環(huán)境質(zhì)量指數(shù)呈降低的趨勢。

2008—2018年，生態(tài)環(huán)境改善主導(dǎo)因素是農(nóng)業(yè)生產(chǎn)用地轉(zhuǎn)型為林地生態(tài)用地和水域生態(tài)用地、工礦生產(chǎn)用地和農(nóng)村用地轉(zhuǎn)為農(nóng)業(yè)生產(chǎn)用地，其他生態(tài)用地轉(zhuǎn)為林地生態(tài)用地(表7)。2008—2013年，工礦生產(chǎn)用地轉(zhuǎn)型為農(nóng)業(yè)生產(chǎn)用地、城鎮(zhèn)用地、農(nóng)村用地，貢獻(xiàn)度分別為39.30%、0.88%、0.29%，主要表現(xiàn)為礦山復(fù)墾對生態(tài)環(huán)境的保護(hù)的貢獻(xiàn)。農(nóng)業(yè)生產(chǎn)用地轉(zhuǎn)為林地生產(chǎn)用地，轉(zhuǎn)移面積56.66 hm2,單一類型生態(tài)貢獻(xiàn)指數(shù)為0.000 157，貢獻(xiàn)率為46.04%，該時期受退耕還林政策的影響，使得大量的耕地轉(zhuǎn)為林地，推動著生態(tài)環(huán)境質(zhì)量的改善。2013—2018年，首次農(nóng)村用地向農(nóng)業(yè)生產(chǎn)用地轉(zhuǎn)型，貢獻(xiàn)度10.40%，表現(xiàn)為農(nóng)村建設(shè)用地轉(zhuǎn)為耕地，這有可能是近年來易地扶貧搬遷、農(nóng)村宅基地整治綜合作用對耕地的保護(hù)作出的貢獻(xiàn)；積極響應(yīng)國家退耕還林政策，農(nóng)業(yè)生產(chǎn)用地轉(zhuǎn)型為林地生態(tài)用地，貢獻(xiàn)度為68.40%；農(nóng)業(yè)生產(chǎn)用地轉(zhuǎn)為其他生態(tài)用地主要表現(xiàn)為存在撂荒的耕地轉(zhuǎn)為其他生態(tài)用地，貢獻(xiàn)度為0.80%；生態(tài)貢獻(xiàn)率為7.6%的其他生態(tài)用地植樹造林。2008—2018年農(nóng)業(yè)生產(chǎn)用地轉(zhuǎn)水域生態(tài)用地，貢獻(xiàn)度分別為2.64%、10%，主要是該時期將耕地開挖為坑塘水面，發(fā)展水產(chǎn)養(yǎng)殖。

林地生態(tài)用地轉(zhuǎn)型為農(nóng)業(yè)生產(chǎn)用地、農(nóng)業(yè)生產(chǎn)用地轉(zhuǎn)型為農(nóng)村用地、農(nóng)業(yè)生產(chǎn)用地轉(zhuǎn)型為工礦生產(chǎn)用地是導(dǎo)致生態(tài)環(huán)境質(zhì)量降低的主要因素。2008—2013年，在社會經(jīng)濟(jì)快速發(fā)展背景下，城鎮(zhèn)擴(kuò)張占周圍的林地、耕地等，林地生態(tài)用地、農(nóng)業(yè)生產(chǎn)用地轉(zhuǎn)型為農(nóng)村用地，生態(tài)貢獻(xiàn)度分別為27.10%、29.40%，與上一時間段相比，整體增幅度變化不大，林地生態(tài)用地轉(zhuǎn)型為城鎮(zhèn)用地生態(tài)貢獻(xiàn)度由5.52%降低至2.62%，農(nóng)業(yè)生產(chǎn)用地轉(zhuǎn)型為農(nóng)村用地由29.40%降低至25.33%；農(nóng)業(yè)生產(chǎn)用地轉(zhuǎn)型為工礦生產(chǎn)用地由8.84%增加至28.55%，主要原因是為推動山東半島城市群體與膠東半島制造業(yè)崛起，連通渤海與黃海區(qū)域，縱貫整個膠東半島，建設(shè)龍青高速。

圖10 “三生空間”生態(tài)環(huán)境質(zhì)量等級空間分布Fig.10 Spatial distribution of ecological environment quality levels in Ecologicl-Production-Living Spaces

表7 影響生態(tài)環(huán)境質(zhì)量的貢獻(xiàn)度Table 7 Contribution to ecological environment quality

3 結(jié)論

基于“三生空間”研究視角，通過對膠東半島沿海低丘陵資源富集區(qū)招遠(yuǎn)市2008年、2013年、2018年土地利用轉(zhuǎn)型及生態(tài)環(huán)境響應(yīng)的研究分析，得到以下結(jié)論：

(1)2008—2018年，沿海低丘陵資源富集區(qū)土地空間擴(kuò)展與收縮處于低速變化。其中農(nóng)村用地擴(kuò)展明顯，擴(kuò)張1 309.07 hm2，主要表現(xiàn)為招遠(yuǎn)經(jīng)濟(jì)新區(qū)建設(shè)影響周圍建設(shè)用地，增長15.31%。整體來說，生活空間與生產(chǎn)空間功能增強(qiáng)，生態(tài)空間減弱是研究“三生空間”的轉(zhuǎn)型的主基調(diào)，其轉(zhuǎn)型由礦產(chǎn)資源的開發(fā)、土地管理政策等多因素造成。

(2)從各功能用地的轉(zhuǎn)移方向與偏移距離來看，服務(wù)及其他生產(chǎn)用地、工礦生產(chǎn)用地、農(nóng)村用地、城鎮(zhèn)用地偏移最顯著，主要表現(xiàn)為城鎮(zhèn)化擴(kuò)張、礦產(chǎn)資源開發(fā)等，使土地利用變化活躍，而水域、耕地、未利用地等受到外界環(huán)境干擾較小，短時期內(nèi)不會發(fā)生劇烈變化。

(3)膠東半島低丘陵沿海區(qū)域生態(tài)環(huán)境質(zhì)量隨時間變化呈降低趨勢，主要表現(xiàn)為礦產(chǎn)資源富集區(qū)集中分布在蠶莊鎮(zhèn)、玲瓏鎮(zhèn)、張星鎮(zhèn)，其他礦產(chǎn)資源零星分布，礦產(chǎn)資源屬于開發(fā)階段，短時間內(nèi)不會對生態(tài)環(huán)境質(zhì)量指數(shù)造成影響;生態(tài)環(huán)境質(zhì)量降低占主導(dǎo)地位，夢芝街道、羅峰街道與溫泉街道作為政治中心、經(jīng)濟(jì)中心、文化中心，水陸空交通極為便利，同時是沿海開放地區(qū)之一，招商引資，發(fā)展輕工業(yè)，推動經(jīng)濟(jì)新區(qū)快速發(fā)展，使生態(tài)環(huán)境質(zhì)量指數(shù)降低；張星鎮(zhèn)的粉絲、石材、機(jī)械三大經(jīng)濟(jì)體發(fā)展迅速，導(dǎo)致城鎮(zhèn)用地不斷向外擴(kuò)張，使得區(qū)域生態(tài)環(huán)境質(zhì)量指數(shù)降低。

(4)2008—2018年，沿海低丘陵資源富集區(qū)生態(tài)環(huán)境質(zhì)量降低的土地利用功能轉(zhuǎn)移類型隨時間變化不斷增多，生態(tài)空間、生產(chǎn)空間轉(zhuǎn)生活空間、生產(chǎn)空間轉(zhuǎn)型是導(dǎo)致生態(tài)環(huán)境質(zhì)量降低的主要因素。