基于LUCC視角山水林田湖草沙一體化保護和修復(fù)工程對生態(tài)系統(tǒng)服務(wù)價值的影響

摘要：山水林田湖草沙一體化保護和修復(fù)工程對于提升生態(tài)環(huán)境質(zhì)量、實現(xiàn)新時代生態(tài)文明目標(biāo)具有重要作用，定量化測算其帶來的生態(tài)效益逐漸成為該研究的難點和熱點。因此，以臨沂市山水工程重點項目區(qū)為研究區(qū)域，基于土地利用變化，建立價值評估方法并探索山水工程實施對于生態(tài)系統(tǒng)服務(wù)價值影響。2020—2023年，研究區(qū)域生態(tài)系統(tǒng)服務(wù)價值總量分別為121.351 7億元和123.441 8億元，總體上升了

2.090 1億元，耕地面積較大，導(dǎo)致其生態(tài)系統(tǒng)服務(wù)價值在總量中的占比也較高，但是林地和水域以較小的面積，卻貢獻了約50%的生態(tài)系統(tǒng)服務(wù)價值，林地和水域面積的變動會較大影響區(qū)域總體的生態(tài)系統(tǒng)服務(wù)價值。本研究對山水工程生態(tài)績效評估，建立生態(tài)保護長效機制，持續(xù)提升生態(tài)系統(tǒng)服務(wù)具有借鑒價值。

關(guān)鍵詞：生態(tài)保護修復(fù)；土地利用變化；生態(tài)系統(tǒng)服務(wù)價值；山水工程；沂蒙山區(qū)

中圖分類號：X171.4 文獻標(biāo)識碼：A 文章編號：1008-9500（2024）10-0-07

61

The Effect of Integrated Protection and Restoration of Landscape， Forest， Lake， Grass and Sand on Ecosystem Service Value from the Perspective of Land Use and Land Cover Change

HUANG Yongxin

（Nanjing Agricultural University， Nanjing 210095， China）

Abstract： The integrated protection and restoration of landscape， forest， Lake， grass and sand plays an important role in improving the quality of the ecological environment and achieving the goal of ecological civilization in the new era.Quantitative measurement of the ecological benefits it brings has gradually become a difficult and hot topic in this research.Therefore，taking the key project area of landscape engineering in Linyi city as the research area，based on land use change，a value evaluation method is established and the impact of landscape engineering implementation on ecosystem service value is explored.From 2020 to 2023， the total value of ecosystem services in the study area will be 12.135 17 billion yuan and 12.344 18 billion yuan respectively， with an overall increase of 209.01 million yuan. Cultivated land accounts for a larger area， which results in a higher proportion of its ecosystem service value in the total amount， but woodland and watersheds， with a relatively smaller area， contribute to about 50% of the ecosystem service value. However， the changes in the area of forest land and watershed will have a greater impact on the overall ecosystem service value of the region. This study is useful for assessing the ecological performance of the Shanshui project， establishing a long-term mechanism for ecological protection， and continuously improving ecosystem services.

Keywords： ecological protection and restoration; land use and land cover change; ecosystem service value; landscape engineering; Yimeng mountain

生態(tài)系統(tǒng)服務(wù)是指人類通過生態(tài)系統(tǒng)直接或間接得到的收益，通常分為供給服務(wù)、調(diào)節(jié)服務(wù)、支持服務(wù)及文化服務(wù)[1]。不同學(xué)者借助經(jīng)濟學(xué)、生態(tài)學(xué)的方法將生態(tài)系統(tǒng)服務(wù)量化，得到生態(tài)系統(tǒng)服務(wù)價值的概念。土地作為人類生產(chǎn)生活的主要空間，經(jīng)常受到人類各種實踐活動而發(fā)生功能結(jié)構(gòu)上的改變，土地利用/覆被變化（Land Use and Land Cover Change，LUCC）在土地利用類型的基礎(chǔ)上定義這些變化[2]，進而成為決定生態(tài)系統(tǒng)服務(wù)功能整體狀況及其空間差異的重要因素，影響生態(tài)系統(tǒng)服務(wù)價值量的大小。山水林田湖草沙一體化保護和修復(fù)工程改變以往針對單一要素進行整治與修復(fù)的傳統(tǒng)方式，通過全域國土空間優(yōu)化生態(tài)系統(tǒng)整體質(zhì)量[3]。針對生態(tài)修復(fù)生態(tài)質(zhì)量的評價關(guān)注于生態(tài)修復(fù)具體過程，效應(yīng)與效果的研究多集中于微觀尺度。而基于LUCC的生態(tài)系統(tǒng)服務(wù)價值評估方法作為一種可以科學(xué)全面衡量區(qū)域生態(tài)系統(tǒng)變化的方法，易于定量測算生態(tài)修復(fù)工程對區(qū)域生態(tài)系統(tǒng)的影響[4]。

山東省臨沂市山水林田湖草沙一體化保護和修復(fù)工程涉及臨沂市10個區(qū)縣，其中沂水縣、費縣、平邑縣、蒙陰縣為國家重點生態(tài)功能區(qū)，以“水源涵養(yǎng)”為生態(tài)主體類型，是山水工程的主要項目區(qū)。2021—2023年，當(dāng)?shù)卣凑铡吧剿痔锖萆呈且粋€生命共同體”理念，通過小流域綜合治理、廢棄礦山修復(fù)等多項具體修復(fù)工程，大幅提高了區(qū)域生態(tài)環(huán)境質(zhì)量。基于此，以沂水縣、費縣、平邑縣、蒙陰縣為例，構(gòu)建山水工程對區(qū)域生態(tài)系統(tǒng)服務(wù)價值影響的定量化評估方法，運用當(dāng)量因子法科學(xué)評估研究區(qū)域的生態(tài)系統(tǒng)服務(wù)總價值和各單項生態(tài)服務(wù)功能價值。

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

1.1 研究區(qū)域概況

研究區(qū)域為沂水縣、蒙陰縣、費縣、平邑縣（地理坐標(biāo)為東經(jīng)117°24′～119°11′，北緯34°22′～36°13′）位于山東省中部偏東南、臨沂市東北部，地近黃海，東連日照，西接棗莊市、濟寧市、泰安市，北靠淄博市、濰坊市。研究區(qū)地處魯中南低山丘陵區(qū)，海拔一般在400 m以上，控制著沂河上游及其主要支流流向。研究區(qū)域概況如圖1所示，區(qū)域內(nèi)地貌主要是中低山丘陵地貌，山高谷深，河流縱橫交錯。研究區(qū)地處中緯度區(qū)，屬于暖溫帶半濕潤大陸性季風(fēng)氣候，具有氣候適宜、四季分明、光照充足等特點。多年平均氣溫為13.4 ℃，多年平均降水量為819.7 mm。研究區(qū)內(nèi)具備森林-水源-農(nóng)田生態(tài)系統(tǒng)類型，主要服務(wù)于水源涵養(yǎng)、洪水調(diào)蓄、土壤保持等生態(tài)功能。

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

研究數(shù)據(jù)包括土地利用遙感影像解譯數(shù)據(jù)、行政邊界數(shù)據(jù)等。研究區(qū)域采用的空間分辨率為30 m的土地利用數(shù)據(jù)，分別來自中國科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心2020年和2023年土地利用遙感監(jiān)測數(shù)據(jù)集，表征山水工程實施前后土地利用狀況，并按照行政邊界，提取研究區(qū)域的土地利用類型數(shù)據(jù)，同時參照其2017年土地利用現(xiàn)狀分類體系，將研究區(qū)域土地利用類型劃分為耕地、林地、草地、水域、建設(shè)用地與未利用地6個一級類。但由于研究區(qū)域的未利用地面積較小，本次研究不選入。此外，研究所涉及的矢量地圖來源于標(biāo)準(zhǔn)地圖服務(wù)系統(tǒng)，研究區(qū)相關(guān)社會經(jīng)濟數(shù)據(jù)來自2020年和2023年《臨沂市統(tǒng)計年鑒》《全國農(nóng)產(chǎn)品成本收益資料匯編》等。

2 研究方法

2.1 土地利用轉(zhuǎn)移矩陣

土地利用轉(zhuǎn)移矩陣來源于系統(tǒng)分析中對系統(tǒng)狀態(tài)與狀態(tài)轉(zhuǎn)移的定量描述[5]，是量化土地利用轉(zhuǎn)移面積的分析方法，可以直觀地獲知特定時間尺度內(nèi)土地利用類型的變換方向和面積變化值，設(shè)i、j

（i， j=1，2，3，…，n）分別為研究開始時和結(jié)束時的土地利用類型，Sij為轉(zhuǎn)移之前第i種土地利用類型向第j種土地利用類型所轉(zhuǎn)移的面積。土地利用轉(zhuǎn)移矩陣表達式為

（1）

式中：n為研究區(qū)域的土地利用類型的數(shù)量。

2.2 土地利用動態(tài)度

土地利用動態(tài)度可以用來反映土地利用轉(zhuǎn)移變化的速度，單一土地利用動態(tài)度反映了某研究區(qū)在一定時間范圍內(nèi)某種土地利用類型數(shù)量變化的速度和幅度。種土地利用類型的動態(tài)度公式為

（2）

式中：K為該種土地利用類型的動態(tài)度；Ua、Ub分別為研究期初和期末該種土地利用類型的面積；T為研究時長。

2.3 生態(tài)系統(tǒng)服務(wù)價值研究

本研究借助謝高地等[6]制定的單位面積生態(tài)系統(tǒng)服務(wù)價值當(dāng)量因子表，通過計算研究區(qū)域1個標(biāo)準(zhǔn)單位生態(tài)系統(tǒng)服務(wù)價值當(dāng)量因子的經(jīng)濟價值量，并結(jié)合研究區(qū)域的實際情況進行修正，最終得到修正后價值當(dāng)量表、研究區(qū)域各土地利用類型生態(tài)系統(tǒng)服務(wù)單價和總價值。

2.3.1 研究區(qū)1個單位標(biāo)準(zhǔn)當(dāng)量因子經(jīng)濟價值計算

根據(jù)謝高地研究成果，1個標(biāo)準(zhǔn)單位生態(tài)服務(wù)價值當(dāng)量因子（以下簡稱標(biāo)準(zhǔn)當(dāng)量）是指1 hm2全國農(nóng)田平均每年單位面積自然糧食產(chǎn)量的經(jīng)濟價值，可以用量化不同類型生態(tài)系統(tǒng)對生態(tài)服務(wù)功能的潛在貢獻能力的作用。根據(jù)已有研究，確定1個標(biāo)準(zhǔn)當(dāng)量的經(jīng)濟價值量等于當(dāng)年全國平均糧食單產(chǎn)市場價值的1/7。計算公式為

（3）

式中：Ea為單位農(nóng)田生態(tài)系統(tǒng)提供生產(chǎn)服務(wù)功能的經(jīng)濟價值，元/hm2；c為作物種類；pc為第c種糧食作物的全國平均市場價格，元/t；qc為第c種糧食作物單位產(chǎn)出（t/hm2）；mc為第c種糧食作物種植面積，hm2；M為g種糧食作物總面積，hm2。

根據(jù)研究區(qū)所在地區(qū)臨沂市2020年和2023年統(tǒng)計年鑒，獲得當(dāng)?shù)刂饕Z食作物小麥和玉米種植面積和單位產(chǎn)出，且由于水稻種植面積極少，不納入計算范圍。結(jié)合2020年和2023兩年《全國農(nóng)產(chǎn)品成本收益資料匯編》中小麥和玉米的平均出售價格，最終求取平均值獲得研究區(qū)單位面積生態(tài)服務(wù)價值當(dāng)量因子

為2 224.85元/hm2。

2.3.2 生態(tài)系統(tǒng)服務(wù)價值當(dāng)量因子表修正

本研究基于謝高地等最新修訂的單位面積生態(tài)系統(tǒng)服務(wù)價值當(dāng)量因子表，包括供給服務(wù)、調(diào)節(jié)服務(wù)、支持服務(wù)、文化服務(wù)4個一級類和11個二級類的生態(tài)服務(wù)分類體系，結(jié)合研究區(qū)域?qū)嶋H情況，修正各生態(tài)系統(tǒng)的服務(wù)價值系數(shù)。

研究區(qū)耕地類型全部為旱地；針葉林、闊葉林、灌叢面積比例為19.7∶66.4∶13.9，根據(jù)加權(quán)方法算的林地對應(yīng)價值當(dāng)量系數(shù)[7]。研究區(qū)草地以灌叢為主，故草地當(dāng)量系數(shù)選擇當(dāng)量表中灌叢當(dāng)量系數(shù)；水域主要為天然池塘、人工水庫和灘地，面積分布為1∶1∶1，前者按照水域系數(shù)，后二者按照濕地系數(shù)加權(quán)計算；建設(shè)用地主要包括城鄉(xiāng)居民點、工礦用地等，其生態(tài)服務(wù)價值影響范圍和程度很小，且現(xiàn)有缺少將建設(shè)用地生態(tài)服務(wù)價值完全量化的方法，暫不納入研究范圍。得到研究區(qū)域生態(tài)服務(wù)價值當(dāng)量表如表1所示。

2.3.3 生態(tài)系統(tǒng)服務(wù)價值測算

根據(jù)以往研究成果，得到研究區(qū)生態(tài)系統(tǒng)服務(wù)總價值和各項服務(wù)價值的計算公式為

（4）

（5）

式中：E為生態(tài)系統(tǒng)服務(wù)總價值；Ak為第k種土地利用類型面積；Vk為第k種土地利用類型的生態(tài)系統(tǒng)服務(wù)價值系數(shù)；Ef為生態(tài)系統(tǒng)第f項服務(wù)功能價值；Vfk為第k種土地利用類型的第f項服務(wù)功能價值系數(shù)；x為土地利用類型總數(shù)量。

為了進一步研究和表征研究區(qū)生態(tài)系統(tǒng)服務(wù)價值空間分布特征，借助ArcGIS 10.8格網(wǎng)尺度法[8]，利用漁網(wǎng)工具，根據(jù)研究區(qū)域的地形、面積、高程等情況，構(gòu)建格網(wǎng)，每個格網(wǎng)內(nèi)的生態(tài)系統(tǒng)服務(wù)價值計算公式為

（6）

式中：Ed為第d個格網(wǎng)內(nèi)的生態(tài)系統(tǒng)服務(wù)價值；As為第s類土地利用類型的面積；Vst為第s類土地利用類型第t種生態(tài)系統(tǒng)服務(wù)功能的單位服務(wù)價值；y為生態(tài)系統(tǒng)服務(wù)功能的總類別數(shù)量。

3 研究結(jié)果和分析

3.1 土地利用變化分析

2020—2023年研究區(qū)土地利用類型如圖2所示。2020年項目實行期初，研究區(qū)域面積占比最多的土地利用類型為耕地，占比為69.82%，2023年項目實施后，面積減少6 015.6 hm2，是面積減少最多的土地利用類型，占比下降了0.8%，動態(tài)度為-0.38%，但仍然是研究區(qū)域內(nèi)占比最高的土地利用類型，這是因為研究區(qū)域整體較為平緩，具有溫暖濕潤的氣候、肥沃的土壤、充足的水源等優(yōu)越的自然條件。雖然區(qū)域內(nèi)部分布較多低山丘陵，但是所在區(qū)域人口數(shù)量眾多，對于糧食和農(nóng)產(chǎn)品需求旺盛，土地開發(fā)利用強度大，存在較多坡耕地和沿坡修建的梯田。

建設(shè)用地為占比第二大的土地利用類型，占比為13.45%，區(qū)域內(nèi)人口多，導(dǎo)致居住用地等建設(shè)用地面積較大，同時從2020—2023年，為面積增加最多的土地利用類型，占比上升0.32%。其次是林地和草地，占比分別為10.22%和4.62%，研究區(qū)內(nèi)中部為沂蒙山系，森林覆蓋較多，并分布有低矮灌叢草地。水域是占比最小的土地利用類型，為1.89%，雖然境內(nèi)有多條河流，并且為水源發(fā)源地，但是位于上游地區(qū)，河道較窄，沒有自然積蓄形成大型湖泊，散布較多零星池塘。總體看，除了耕地處在下降趨勢，其他4個土地利用類型面積均有所上升，動態(tài)度較高的為草地和水域，分別為2.3%和1.95%，受到了重視和保護。

山水工程對不同土地利用類型的變化和轉(zhuǎn)移方向具體作用如表3所示。耕地受坡耕地還草還林、水土保持、流域治理等工程措施影響，向草地、林地、水域轉(zhuǎn)移10 859.58 hm2，占比為78.13%，縣鄉(xiāng)等建設(shè)用地持續(xù)擴張，導(dǎo)致有3 039.66 hm2土地被侵占。草地主要轉(zhuǎn)向耕地和林地；建設(shè)用地由于社會經(jīng)濟持續(xù)發(fā)展，轉(zhuǎn)移面積不大，主要為耕地和水域；林地受到保育等工程，也未發(fā)生較大轉(zhuǎn)移，主要轉(zhuǎn)向耕地和草地。水域由于難以開發(fā)利用，并且受到工程高度保護和重視，受小流域綜合治理和水環(huán)境治理等工程，轉(zhuǎn)移面積最小。

3.2 生態(tài)系統(tǒng)服務(wù)價值分析

3.2.1 生態(tài)系統(tǒng)服務(wù)價值總體變化

2020—2023年研究區(qū)不同土地利用類型的生態(tài)系統(tǒng)服務(wù)價值如表4所示，2020年研究區(qū)域的生態(tài)系統(tǒng)服務(wù)價值為121.351 7億，生態(tài)修復(fù)工程實施之后，生態(tài)系統(tǒng)服務(wù)價值增加了2.090 1億元，至123.441 8億元，變化率為1.72%。除了耕地由于面積減少導(dǎo)致價值其生態(tài)系統(tǒng)服務(wù)價值下降0.536 7億元外，其他4個土地利用類型均增加，草地和水域增加值最多為1.049 1億元和1.409 7億元。但總體來看，耕地在總體價值的占比最高，為37.42%；其次是林地，貢獻28.71%；雖然水域的面積最小，但是其單位面積生態(tài)系統(tǒng)服務(wù)價值較大，貢獻度占比為20.71%，說明水資源型生態(tài)系統(tǒng)對于生態(tài)服務(wù)功能的重要性程度，需要持續(xù)貫徹執(zhí)行流域治理、設(shè)置水源自然保護區(qū)等保護和修復(fù)工程。林地和水域生態(tài)系統(tǒng)應(yīng)是改善區(qū)域生態(tài)環(huán)境、保障區(qū)域生態(tài)安全的關(guān)鍵因素。

3.2.2 生態(tài)系統(tǒng)服務(wù)價值單項變化

2020—2023年研究區(qū)單項生態(tài)系統(tǒng)服務(wù)價值變化如表5所示?？梢园l(fā)現(xiàn)，修復(fù)工程實施之后，生態(tài)服務(wù)系統(tǒng)的單項服務(wù)功能價值除了食物生產(chǎn)、原料生產(chǎn)、維持養(yǎng)分循環(huán)3項由于耕地面積下降較多而發(fā)生相應(yīng)降低之后，其他各項服務(wù)價值均有不同程度上的增加，其中水文調(diào)節(jié)服務(wù)功能價值增加最多為1.127 2億元，

同時水文調(diào)節(jié)服務(wù)功能在總體生態(tài)系統(tǒng)服務(wù)總價值中的占比最高，為24.6%。這是由于林地、草地、水域等的單位面積水文調(diào)節(jié)功能均較高，3類平均值為43 254.72元/hm2，尤其是水域為111 761.63元/hm2。

各項功能中，占比最少的為水資源供給和維持養(yǎng)分循環(huán)，美學(xué)景觀、原料生產(chǎn)、凈化環(huán)境適中，食物生產(chǎn)、氣體調(diào)節(jié)、氣候調(diào)節(jié)、土壤保持以及生物多樣性對生態(tài)系統(tǒng)總體服務(wù)功能價值的貢獻較高。

總之，該地區(qū)生態(tài)系統(tǒng)服務(wù)功能主要以氣體調(diào)節(jié)、氣候調(diào)節(jié)、水文調(diào)節(jié)以及土壤保持為主，因此應(yīng)當(dāng)繼續(xù)堅守水域、林地面積不被轉(zhuǎn)換，劃定并嚴(yán)格落實林地生態(tài)保護紅線和生態(tài)保護自然區(qū)，對進一步提升區(qū)域生態(tài)系統(tǒng)服務(wù)價值有重要意義。

3.2.3 生態(tài)系統(tǒng)服務(wù)價值空間分異

為更好地展現(xiàn)研究區(qū)域生態(tài)系統(tǒng)服務(wù)價值的空間分異規(guī)律[9]，本研究通過ArcGIS 10.8軟件測算格網(wǎng)內(nèi)的生態(tài)系統(tǒng)服務(wù)價值，來表示不同地區(qū)的單位生態(tài)系統(tǒng)服務(wù)價值并繪制，并通過自然斷面法將生態(tài)系統(tǒng)服務(wù)價值劃分為低、較低、中、較高、高5個級別，如圖3所示。

結(jié)合研究區(qū)域土地利用現(xiàn)狀可發(fā)現(xiàn)，區(qū)域生態(tài)系統(tǒng)服務(wù)價值呈現(xiàn)生態(tài)服務(wù)價值低值區(qū)域集中連片分布、高值區(qū)域零星分布的總體布局，零散水體區(qū)域的價值高，中部林地覆蓋沂蒙山區(qū)價值較高，南部和北部草地集中分布區(qū)價值適中，耕地與建設(shè)用地價值低的空間分異規(guī)律。

4 討論與分析

本研究首先探討了各項修復(fù)工程具體措施對于土地利用類型變化的作用效果，并且利用了謝高地等[6]研究的生態(tài)系統(tǒng)服務(wù)價值當(dāng)量因子法進行服務(wù)價值測算，主要用于評估多尺度區(qū)域生態(tài)系統(tǒng)服務(wù)價值變化。該方法的關(guān)鍵是確定區(qū)域地類面積及單位面積生態(tài)系統(tǒng)服務(wù)價值，因此先根據(jù)研究區(qū)域?qū)嶋H情況修正并建立單位面積生態(tài)系統(tǒng)服務(wù)價值當(dāng)量表，后借助ArcGIS 10.8軟件工具，在600 m×600 m格網(wǎng)尺度下解譯分析工程實施前后2020年和2023年兩年的土地利用類型面積和變化狀況，評估出臨沂市山水林田湖草沙一體化保護和修復(fù)工程實施前后的總體生態(tài)系統(tǒng)服務(wù)價值、各項生態(tài)功能服務(wù)價值和變化量，進而展現(xiàn)和說明該項生態(tài)修復(fù)工程對于當(dāng)?shù)厣鷳B(tài)系統(tǒng)服務(wù)價值的影響。項目實施前后（2020—2023年），研究區(qū)域（沂水縣、蒙陰縣、平邑縣、費縣）各類土地利用類型占比排序均為耕地＞建設(shè)用地＞林地＞草地＞水域，土地利用空間分布格局基本保持穩(wěn)定，各土地利用類型相互之間產(chǎn)生了不同程度轉(zhuǎn)移。從面積變化的絕對值來看，耕地失去了最多的土地，而草地和建設(shè)用地則增加了最多，水域和林地也獲得提升，這與生態(tài)修復(fù)工程的退耕還林還草、水土保持、小流域綜合治理密不可分，但由于經(jīng)濟發(fā)展和城市化的必然性，建設(shè)用地依然獲得面積增加。從面積變化的動態(tài)度來看，草地和水域動態(tài)度較高，說明面積較原來獲得增加幅度大，并且有繼續(xù)擴大趨勢[10]。未來，為保持工程實施效果的長期性，要繼續(xù)落實林地草地育種和水土保持等工程，防止草地等地類轉(zhuǎn)化回建設(shè)用地和耕地，同時減緩建設(shè)用地向其他地類侵占的速度。

5 結(jié)論

文章以臨沂市山水工程重點項目區(qū)為研究區(qū)域，探究基于LUCC視角山水林田湖草沙一體化保護和修復(fù)工程對生態(tài)系統(tǒng)服務(wù)價值的影響。得出結(jié)論：工程項目實施后，研究區(qū)土地利用總體格局保持穩(wěn)定，耕地面積減少最多，林地、草地、水域面積因為相關(guān)工程措施獲得不同程度增加；山水工程實施前，地區(qū)生態(tài)系統(tǒng)服務(wù)價值為121.351 7億元，實施后增加至123.441 8億元，漲幅1.02%，林地、水域以較少的面積貢獻了較高的生態(tài)系統(tǒng)服務(wù)價值，占49.42%，以林地和水域作為生態(tài)修復(fù)和保護工程的重點單元是保障生態(tài)安全、提升生態(tài)質(zhì)量的關(guān)鍵；生態(tài)服務(wù)價值高值區(qū)集中在中部沂蒙山區(qū)，以及河流和湖泊等水環(huán)境用地，應(yīng)以此優(yōu)化生態(tài)修復(fù)布局和形成長期專項保護治理工作。臨沂市山水工程要繼續(xù)堅持以河流、濕地等水生生態(tài)系統(tǒng)和沂蒙山區(qū)林地生態(tài)系統(tǒng)作為生態(tài)修復(fù)和空間綜合治理的核心單元，杜絕水域和林地被城鎮(zhèn)居民點等建設(shè)用地侵占。

參考文獻

1 COSTANZA R，D’ARGE R，DEGROOT R.et al.

The value of the world’s ecosystem services and natural capital[J].Nature，1997（6630）：253-260.

2 傅伯杰，張立偉.土地利用變化與生態(tài)系統(tǒng)服務(wù)：概念、方法與進展[J].地理科學(xué)進展，2014（4）：441-446.

3 唐秀美，劉 敏，吳 悠，等.國土綜合整治與生態(tài)修復(fù)對區(qū)域生態(tài)系統(tǒng)服務(wù)價值的影響評估[J].生態(tài)學(xué)報，2024（14）：5974-5984.

4 汪燕林，岳文澤，張水軍，等.促進人與自然和諧共生的實踐與啟示：以浙江甌江源頭區(qū)域山水林田湖草沙一體化保護修復(fù)工程為例[J].中國國土資源經(jīng)濟，2023（41nBtN2YDViedb1MBxPUtVZ2Y9SfQyP+PDGepfTX4nug=）：37-43.

5 劉 瑞，朱道林.基于轉(zhuǎn)移矩陣的土地利用變化信息挖掘方法探討[J].資源科學(xué)，2010（8）：1544-1550.

6 謝高地，張彩霞，張雷明，等.基于單位面積價值當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價值化方法改進[J].自然資源學(xué)報，2015（8）：1243-1254.

7 姜 晗，吳 群.基于LUCC的江蘇省生態(tài)系統(tǒng)服務(wù)價值評估及時空演變特征研究[J].長江流域資源與環(huán)境，2021（11）：2712-2725.

8 危小建，辛思怡，張穎藝，等.不同格網(wǎng)尺度下生態(tài)系統(tǒng)服務(wù)價值空間分異及其影響因素差異：以大南昌都市圈為例[J].生態(tài)學(xué)報，2023

（18）：7585-7597.

9 孟雪源，陳 剛，鄭志元.國家重點生態(tài)功能區(qū)生態(tài)系統(tǒng)服務(wù)價值對土地利用變化的時空響應(yīng)[J].浙江農(nóng)林大學(xué)學(xué)報，2024（4）：830-840.

10 王 晶，楊亮彥，孔 輝，等.治溝造地背景下延安市土地利用與生態(tài)系統(tǒng)服務(wù)價值變化[J].生態(tài)科學(xué)，2023（5）：57-65.