1980—2020年濟(jì)南市土地利用變化對生態(tài)系統(tǒng)服務(wù)價(jià)值的影響

張夢虹 姜鵬 趙麗 趙靜 劉維國 紀(jì)東 時(shí)海 牟剛 王靜

摘要：面對資源緊缺、土地利用方式不合理、生態(tài)系統(tǒng)退化等問題，探索土地利用變化對生態(tài)系統(tǒng)服務(wù)價(jià)值的影響，對于制定適宜的國土空間規(guī)劃、自然資源管理、生態(tài)文明建設(shè)等方面的政策措施具有重要意義。本文選取濟(jì)南市作為研究區(qū)域，探討了1980—2020年研究區(qū)土地利用結(jié)構(gòu)變化以及土地利用類型間轉(zhuǎn)移情況，并核算了生態(tài)系統(tǒng)服務(wù)價(jià)值，分析了不同土地利用類型對生態(tài)系統(tǒng)服務(wù)價(jià)值變化的生態(tài)貢獻(xiàn)率。研究結(jié)果表明：①在研究時(shí)段內(nèi)，研究區(qū)建設(shè)用地面積呈現(xiàn)持續(xù)增加的趨勢；在1980—2000年階段，耕地、林地、草地、水域面積呈現(xiàn)波動式變化，在2000—2020年階段處于持續(xù)下降的趨勢。②1980—2020年，生態(tài)系統(tǒng)服務(wù)總價(jià)值呈先下降后上升再下降的趨勢；從單項(xiàng)服務(wù)價(jià)值來看，供給、調(diào)節(jié)、支持和文化4種主要的服務(wù)價(jià)值均有不同程度的減少，尤其在2005—2010年，下降幅度最大。③耕地對濟(jì)南市生態(tài)系統(tǒng)服務(wù)價(jià)值變化影響最大，主要由于研究區(qū)內(nèi)受城市發(fā)展的影響，導(dǎo)致耕地面積大幅減少，其次影響較大的是林地、草地和水域。

關(guān)鍵詞：土地利用變化；生態(tài)系統(tǒng)服務(wù)價(jià)值；生態(tài)貢獻(xiàn)率；濟(jì)南市

中圖分類號：P301.2????文獻(xiàn)標(biāo)識碼：A????doi：10.12128/j.issn.1672-6979.2024.02.010

引文格式：張夢虹，姜鵬，趙麗，等.1980—2020年濟(jì)南市土地利用變化對生態(tài)系統(tǒng)服務(wù)價(jià)值的影響［J］.山東國土資源，2024，40（2）：5664. ZHANG Menghong， JIANG Peng， ZHAO Li， et al. Influence of Land Use Changes on the Value of Ecosystem Services in Ji'nan from 1980 to 2020［J］.Shandong Land and Resources，2024，40（2）：5664.

0?引言

改革開放以來，我國經(jīng)濟(jì)進(jìn)入高速發(fā)展階段，城鎮(zhèn)化進(jìn)程加快，城鎮(zhèn)化率已由17.92%提升至2021年的64.72%［12］。社會經(jīng)濟(jì)轉(zhuǎn)型發(fā)展以土地利用轉(zhuǎn)型為空間載體，對農(nóng)業(yè)空間和生態(tài)空間的占用難以避免，進(jìn)而加劇區(qū)域土地利用結(jié)構(gòu)趨于不合理，直接影響區(qū)域生態(tài)系統(tǒng)服務(wù)水平［3］。生態(tài)系統(tǒng)服務(wù)價(jià)值（Ecological Service Value，ESV）是人類直接或間接通過生態(tài)系統(tǒng)獲得的所有惠益，主要包括供給、調(diào)節(jié)、文化和支持四種服務(wù)價(jià)值［45］。生態(tài)系統(tǒng)服務(wù)對維持地球生命系統(tǒng)的正常運(yùn)轉(zhuǎn)至關(guān)重要，它們直接或間接地為人類提供福祉［6］。聯(lián)合國千年生態(tài)系統(tǒng)評估（The Millennium Ecosystem Assessment）中指出，近100多年來全球60%的生態(tài)系統(tǒng)正在或已經(jīng)退化［7］，很大程度在于對生態(tài)系統(tǒng)服務(wù)價(jià)值認(rèn)識不足而盲目開發(fā)所導(dǎo)致［8］。土地利用變化不僅會直接影響土地利用結(jié)構(gòu)、功能和景觀格局的變化［911］，而且也會間接影響生態(tài)系統(tǒng)服務(wù)功能的變化。

近年來，生態(tài)系統(tǒng)服務(wù)價(jià)值的定量評價(jià)已成為國內(nèi)學(xué)者的研究熱點(diǎn)［12］。Costanza等首先提出了生態(tài)系統(tǒng)服務(wù)價(jià)值估算的原理和方法［13］，并得到了廣泛的應(yīng)用；但由于地域性的差異，不同地區(qū)直接采用該方法進(jìn)行核算，研究結(jié)果的誤差較大，會影響生態(tài)系統(tǒng)服務(wù)價(jià)值估算結(jié)果的可靠性。鑒于此，謝高地等在該研究基礎(chǔ)上，提出了適合我國的生態(tài)系統(tǒng)服務(wù)評估體系，該方法目前應(yīng)用廣泛［1419］。國內(nèi)眾多學(xué)者采用謝高地改進(jìn)的當(dāng)量因子法，分析了土地利用/覆被變化（Land Use / Cover Change，LUCC）下的ESV變化情況及其相互間的影響［2021］；不少學(xué)者對濟(jì)南市生態(tài)系統(tǒng)服務(wù)價(jià)值開展了系列研究和討論［2225］，但研究時(shí)段主要集中在2000年以后，對更早期的研究時(shí)段較少，且研究不同土地利用類型對生態(tài)系統(tǒng)服務(wù)價(jià)值變化的生態(tài)貢獻(xiàn)率的較少，不能較好的體現(xiàn)兩者間的內(nèi)在關(guān)聯(lián)。

本文基于長時(shí)間跨度的土地利用/覆被數(shù)據(jù)，采用謝高地等提出的當(dāng)量因子法核算濟(jì)南市ESV，探討1980—2020年濟(jì)南市土地利用變化對生態(tài)系統(tǒng)服務(wù)價(jià)值的影響，研究二者間相互關(guān)系，以期為濟(jì)南市在國土空間規(guī)劃、自然資源管理、生態(tài)文明建設(shè)提供決策依據(jù)。

1?數(shù)據(jù)與方法

1.1?研究區(qū)概況

濟(jì)南市是山東省省會，地勢由南向北依次為低山丘陵、山前沖積—洪積傾斜平原和黃河沖積平原的地貌形態(tài)。濟(jì)南市呈現(xiàn)出北部臨黃帶，中部為山前平原帶，南部為丘陵山區(qū)帶的地形特征（圖1）。其四周與德州、濱州、淄博、泰安、聊城等地市相鄰，全市總面積10244.45 km2。

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

土地利用/覆被數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心，Landsat遙感影像是該LUCC的主要數(shù)據(jù)源，解譯方法以人工目視解譯為主。受數(shù)據(jù)源的影響，1980年LUCC數(shù)據(jù)主要基于Land-MSS遙感影像數(shù)據(jù)解譯得到，1990—2010年LUCC主要基于Landsat-TM/ETM遙感影像數(shù)據(jù)解譯得到，2015、2020年LUCC更新主要使用Landsat 8數(shù)據(jù)。主要糧食產(chǎn)量、市場價(jià)格、種植面積數(shù)據(jù)來源于《濟(jì)南市2021年統(tǒng)計(jì)年鑒》和《全國農(nóng)產(chǎn)品成本收益資料匯編2019》。1.3?研究方法

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

土地利用轉(zhuǎn)移矩陣用于反映在研究時(shí)段內(nèi)濟(jì)南不同地類面積間相互轉(zhuǎn)化的動態(tài)過程［2627］。公式（1）如下：

Sij=S11?S12?…?S1n

S21?S22?…?S2n

…??…?…?…

Sn1?Sn2?…?Snn（1）

式中：S為面積；n為濟(jì)南市土地利用類型數(shù)；i，j（i，j=1，2，…，n）分別是土地類型變化前后所代表的土地類型；Sij為i與j地類間轉(zhuǎn)換的面積量［28］。

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

本文以修正后的中國生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)當(dāng)量表為依據(jù)［6，29］，結(jié)合濟(jì)南市社會經(jīng)濟(jì)發(fā)展情況，對濟(jì)南市單位面積糧食產(chǎn)量創(chuàng)造的經(jīng)濟(jì)價(jià)值進(jìn)行修正。鑒于建設(shè)用地的生態(tài)系統(tǒng)服務(wù)價(jià)值的特殊性，本文以謝高地等修正后的當(dāng)量因子表為依據(jù)，將建設(shè)用地系數(shù)設(shè)為0。

通過查詢《濟(jì)南市2021年統(tǒng)計(jì)年鑒》可以看出，1980—2020年濟(jì)南平均糧食單產(chǎn)為5130.49 kg·hm2；查詢《全國農(nóng)產(chǎn)品成本收益資料匯編2019》獲取2018年山東地區(qū)粳稻、小麥、玉米、大豆4種主要糧食作物的價(jià)格，計(jì)算得到4種糧食作物的平均價(jià)格為2.846元·kg1，并將其作為濟(jì)南市平均糧食價(jià)格，根據(jù)“一個(gè)標(biāo)準(zhǔn)生態(tài)系統(tǒng)服務(wù)經(jīng)濟(jì)價(jià)值等效系數(shù)是單位面積農(nóng)田的食物生產(chǎn)經(jīng)濟(jì)價(jià)值的1/7”［30］，核算出濟(jì)南市一個(gè)生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子為2085.91元·hm2，其次通過與修正后的當(dāng)量表相乘得出濟(jì)南生態(tài)系統(tǒng)單位面積服務(wù)價(jià)值系數(shù)表，最后以此核算濟(jì)南市ESV（表1）。

1.3.3?生態(tài)貢獻(xiàn)率

生態(tài)貢獻(xiàn)率表示不同土地利用類型在某一時(shí)段內(nèi)的服務(wù)價(jià)值變化量對區(qū)域總服務(wù)價(jià)值變化量的影響。公式（2）為：

SKT=|ΔESVKT|∑6K=1|ΔESVKT|（2）

式中：SKT代表的是地類K在時(shí)間段T內(nèi)的生態(tài)貢獻(xiàn)率；ΔESVKT代表的是地類K在時(shí)間段T內(nèi)的ESV變化量［31］。

2?結(jié)果與分析

2.1?濟(jì)南土地利用/覆被變化分析

2.1.1?土地利用類型結(jié)構(gòu)變化分析

1980—2020年間，濟(jì)南市土地利用類型中，耕地面積占比最大，維持在國土總面積的55%～63%之間，主要分布于濟(jì)南市北部、西南和東南等區(qū)域；林地、草地、建設(shè)用地面積占比在8%～20%之間，林地、草地主要分布于中部山地區(qū)域，建設(shè)用地主要分布于中部地區(qū)以及零散分布于耕地之間；水域面積占比較小，在2%～3%之間（表2，圖2）。

40年間耕地、林地、草地、水域、建設(shè)用地和未利用地6種土地類型中，建設(shè)用地面積呈現(xiàn)持續(xù)擴(kuò)張的趨勢，耕地、林地、草地呈現(xiàn)面積均呈現(xiàn)下降的趨勢，水域面積變化不大。在研究時(shí)段內(nèi)，建設(shè)用地總共增加1069.53 km2，面積翻了一倍；其中增長最快的階段在2005—2010年期間。1980—1990年，增加了12.02 km2；1990—1995年，增加了82.53 km2；1995—2000年，增加了19.01 km2；2000—2005年，增加了108.84 km2；2005—2010年，增加了722.59 km2；2010—2015年，增加了83.64 km2；2015—2020年，增加了40.90 km2。相較于1980年，2020年耕地面積減少691.14 km2，2005—2010年期間耕地面積減少最多；其中1980—1990年耕地增加了42.25 km2，1990—1995年減少了155.46 km2，1995—2000年增加了58.37 km2，2000—2005年減少了136.74 km2，2005—2010年減少了373.45 km2，2010—2015年減少了84.98 km2，2015—2020年減少了41.09 km2。林地在1980—2020年期間，面積減少了84.17 km2，2005—2010年期間林地面積減少最多；其中1980—1990年林地面積減少了34.57 km2，1990—1995年面積又增加了14.04 km2，1995—2000年減少了12.60 km2，2000—2005年增加了4.79 km2，2005—2010年減少了55.18 km2，2010—2015年減少了1.49 km2，2015—2020年增加了0.84 km2。1980—2020年期間，草地面積減少了217.44 km2，2005—2010年期間草地面積減少最多；其中1980—1990年草地面積減少了5.78 km2，1990—1995年面積又增加了13.17 km2，1995—2000年減少了13.45 km2，2000—2005年減少了3.41 km2，2005—2010年減少了207.19 km2，2010—2015年減少了1.88 km2，2015—2020年增加了1.10 km2。相比于1980年，2020年水域面積減少了15.16 km2，水域面積整體變化不大。

2.1.2?土地利用轉(zhuǎn)移矩陣分析

從土地利用轉(zhuǎn)移矩陣（表3）可以發(fā)現(xiàn)，1980—2020年期間，耕地和建設(shè)用地互相流轉(zhuǎn)頻繁，其他地類間流轉(zhuǎn)較少。耕地流出部分主要轉(zhuǎn)移為建設(shè)用地，占流出總面積的90.36%；耕地流入部分主要來源于建設(shè)用地和草地，分別占流入總面積的39.46%和26.87%。

從土地利用轉(zhuǎn)移矩陣空間分布（圖3）來看，LUCC變化最為頻繁的區(qū)域主要位于城鎮(zhèn)、居民地附近，是耕地和建設(shè)用地間的轉(zhuǎn)換。在城鎮(zhèn)化快速發(fā)展階段，建設(shè)用地面積的迅速擴(kuò)張必將擠占其他地類，而濟(jì)南市主要以平原耕種區(qū)域?yàn)橹鳎虼烁乇唤ㄔO(shè)用地占用；在進(jìn)入高質(zhì)量發(fā)展階段后，國家嚴(yán)格控制耕地紅線、違法建筑拆除復(fù)墾等系列舉措，耕地在一定程度尚有所增加。

2.2?濟(jì)南生態(tài)系統(tǒng)服務(wù)價(jià)值變化

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

1980—2020年期間，濟(jì)南市生態(tài)系統(tǒng)服務(wù)總價(jià)值呈先下降后上升再下降的趨勢（表4），總價(jià)值從1980年的237.75億元下降到1990年的234.85億元，1995年上升至235.38億元，2020年下降至214.52億元。1980—2020年，濟(jì)南市生態(tài)系統(tǒng)服務(wù)價(jià)值總體下降23.23億元，減幅約10%，主要是由于對生態(tài)系統(tǒng)服務(wù)價(jià)值無影響的建設(shè)用地面積增加，而耕地、林地、草地等面積減少，導(dǎo)致濟(jì)南市生態(tài)系統(tǒng)服務(wù)總價(jià)值量下降。

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

通過對濟(jì)南市1980—2020年單項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值變化分析發(fā)現(xiàn)，每項(xiàng)ESV均有不同程度的減少（表5），尤其在2005—2010年，下降幅度最大。在研究期內(nèi)，供給服務(wù)價(jià)值下降了2.98億元，減幅為9.82%，1980—2000年期間波動幅度較小，主要在于該時(shí)段經(jīng)濟(jì)發(fā)展緩慢，城鎮(zhèn)化率不高，耕地、林地、草地面積變化小，生產(chǎn)糧食和原材料能力較好；2005—2010年下降幅度較大，減少了1.87億元，主要原因在于該時(shí)段內(nèi)濟(jì)南市城鎮(zhèn)化進(jìn)程加快，建設(shè)用地面積快速增加，導(dǎo)致其他地類面積均有所減少，進(jìn)而促使生產(chǎn)糧食和原材料的能力有所下降；2010—2020年降幅較小，主要是由于經(jīng)濟(jì)發(fā)展模式逐漸從粗放式向高質(zhì)量高效益模式轉(zhuǎn)變，嚴(yán)格保護(hù)耕地，實(shí)現(xiàn)耕地占補(bǔ)平衡，同時(shí)大力保護(hù)生態(tài)環(huán)境，從而使得該時(shí)段內(nèi)供給服務(wù)價(jià)值整體平穩(wěn)。

在研究時(shí)段內(nèi)，濟(jì)南市的調(diào)節(jié)服務(wù)價(jià)值呈衰減趨勢，1980—2005年期間下降幅度小，減少了1.85億元，下降幅度為1.47%；2005—2010年時(shí)段內(nèi)下降幅度大，減少了9.46億元，下降幅度為7.65%；2010—2020年較為穩(wěn)定，波動小，減少了0.62億元，下降幅度為0.54%。調(diào)節(jié)服務(wù)價(jià)值變化受耕地、林地、草地的影響較大，這主要由3種地類在生態(tài)系統(tǒng)中的功能和特性所決定的，同時(shí)3者面積占比較大所造成的。但是在研究時(shí)段內(nèi)，3者的面積持續(xù)減少，導(dǎo)致了在氣候、水文等方面的調(diào)節(jié)能力有所下降，從而導(dǎo)致調(diào)節(jié)服務(wù)價(jià)值的降低。

支持服務(wù)價(jià)值主要反映保持土壤和微生物多樣性的情況，耕地、林地和草地在該方面起到了重要作用，因此支持服務(wù)價(jià)值的變化受三者變化的影響較大。1980—2020年間，支持服務(wù)價(jià)值呈現(xiàn)持續(xù)減少的趨勢，但各階段的減少幅度存在差異。1980—2005年，支持服務(wù)價(jià)值下降了1.53億元；2005—2010年期間，減少了4.99億元；2010—2020年期間，減少了0.63億元。其主要原因是耕地、林地和草地在研究時(shí)段內(nèi)都持續(xù)減少，導(dǎo)致1980—2020年期間內(nèi)支持服務(wù)價(jià)值降低。

濟(jì)南市文化服務(wù)價(jià)值受林地、水域影響較大，但在1980—2020年間，林地和水域面積持續(xù)下降，使其美學(xué)景觀和文化服務(wù)價(jià)值均持續(xù)降低。主要是由于大量的林地和水域流轉(zhuǎn)為耕地、建設(shè)用地等其他地類，導(dǎo)致兩者面積持續(xù)減少，直接導(dǎo)致文化服務(wù)價(jià)值的下降。

2.3?LUCC對生態(tài)系統(tǒng)服務(wù)價(jià)值變化的貢獻(xiàn)率

根據(jù)表6的結(jié)果表明，1980—2005年、2005—2010年、2010—2020年、1980—2020年4個(gè)不同時(shí)段內(nèi)，耕地對生態(tài)系統(tǒng)服務(wù)價(jià)值變化的貢獻(xiàn)率最大，尤其在2010—2020年間高達(dá)78.20%，主要由于研究區(qū)域耕地被人類活動大量占用，導(dǎo)致其生態(tài)系統(tǒng)服務(wù)價(jià)值變化最大；其次是林地，1980—2005年間對ESV的貢獻(xiàn)率最高，為28.38%。1980—2005年間，耕地對ESV的貢獻(xiàn)率最大，其他土地利用類型對ESV的貢獻(xiàn)率從大到小依次為林地、水域、草地、未利用地；2005—2010年，耕地對ESV的貢獻(xiàn)率最大，其他土地利用類型對ESV的貢獻(xiàn)率從大到小依次為草地、林地、水域、未利用地；2010—2020年，耕地對ESV的貢獻(xiàn)率最大，其他土地利用類型對ESV的貢獻(xiàn)率從大到小依次為水域、林地、草地、未利用地；1980—2020年，耕地對ESV的貢獻(xiàn)率最大，其他土地利用類型對ESV的貢獻(xiàn)率從大到小依次為草地、林地、水域、未利用地。

3?結(jié)論與討論

本文以濟(jì)南市為研究區(qū)域，采用謝高地當(dāng)量因子法，探討了1980—2020年間研究區(qū)內(nèi)的土地利用變化情況，并以此為基礎(chǔ)對濟(jì)南市的生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行了核算，最后通過分析不同土地利用類型的生態(tài)貢獻(xiàn)率掌握了ESV變化的主要影響因素。得出以下3點(diǎn)結(jié)論：

（1）在研究時(shí)段內(nèi)，研究區(qū)在土地利用結(jié)構(gòu)上，主要以耕地為主，其次為建設(shè)用地、林地、草地；從變化趨勢上看，1980—2020年間建設(shè)用地面積呈現(xiàn)持續(xù)增加的趨勢，在1980—2000年階段，耕地、林地、草地、水域面積呈現(xiàn)波動式變化，在2000—2020年階段處于持續(xù)下降的趨勢。

（2）1980—2020年，從生態(tài)系統(tǒng)服務(wù)總價(jià)值來看，整體呈先下降后上升再下降的趨勢；從單項(xiàng)服務(wù)價(jià)值來看，供給服務(wù)價(jià)值、調(diào)節(jié)服務(wù)價(jià)值、支持服務(wù)價(jià)值和文化服務(wù)價(jià)值均有不同程度的減少，尤其在2005—2010年，下降幅度最大。

（3）生態(tài)貢獻(xiàn)率結(jié)果表明，耕地對濟(jì)南市生態(tài)系統(tǒng)服務(wù)價(jià)值變化影響最大，主要由于研究區(qū)內(nèi)受城市發(fā)展的影響，導(dǎo)致耕地面積大幅減少，其次影響較大的是林地、草地和水域。

本文在宏觀尺度上對濟(jì)南市的生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行分析，研究結(jié)果在一定程度上可作為未來城市發(fā)展提供決策依據(jù)。但受土地利用數(shù)據(jù)空間分辨率、農(nóng)作物市場價(jià)格波動等因素的影響，可能對研究成果產(chǎn)生一定影響。下一步作者將從提高土地利用數(shù)據(jù)分辨率、修正市場價(jià)格的角度，盡可能減少二者對研究結(jié)果的影響。

參考文獻(xiàn)：

［1］?國家統(tǒng)計(jì)局，中國統(tǒng)計(jì)年鑒［M］．北京：中國統(tǒng)計(jì)出版社，1999.

［2］?國家統(tǒng)計(jì)局，中國統(tǒng)計(jì)年鑒［M］．北京：中國統(tǒng)計(jì)出版社，2022.

［3］?關(guān)梅，張文信，蔣海明，等.山東省土地利用變化的生態(tài)系統(tǒng)服務(wù)交叉敏感性評價(jià)［J］.中國農(nóng)業(yè)大學(xué)學(xué)報(bào)，2022，27（6）：192203.

［4］?謝高地，張彩霞，張雷明，等.基于單位面積價(jià)值當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價(jià)值化方法改進(jìn)［J］.自然資源學(xué)報(bào)，2015，30（8）：12431254.

［5］?馬麗娜，張飛云，翟玉鑫，等.1980—2020年新疆土地利用變化下生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空演變分析［J］.干旱區(qū)地理，2023，46（2）：253263.

［6］?Costanza R ，Arge，Groot R D，et al.The value of the world's ecosystem services and natural capital［J］.Nature，1997，387（15）：253260.

［7］?Assesment M E.Ecosystems and human wellbeing： synthesis［J］.Physics Teacher，2005，34（9）：534538.

［8］?雷軍成，王莎，汪金梅，等.土地利用變化對尋烏縣生態(tài)系統(tǒng)服務(wù)價(jià)值的影響［J］.生態(tài)學(xué)報(bào)，2019，39（9）：30893099.

［9］?DONG G，XU E，ZHANG H .Spatiotemporal Variation of Driving Forces for Settlement Expansion in Different Types of Counties［J］.Sustainability，2015，8（1）：39.

［10］?劉冰，王靖偉，吳福榮，等.黃河三角洲景觀格局變化及生態(tài)風(fēng)險(xiǎn)評價(jià)［J］.山東國土資源，2023，39（5）：3945.

［11］?劉冰，王海明，吳明陽，等.曹縣黃河故道國家濕地公園整合優(yōu)化后景觀格局變化［J］.甘肅林業(yè)科技，2023，48（1）：5458.

［12］?劉金勇，孔繁花，尹海偉，等.濟(jì)南市土地利用變化及其對生態(tài)系統(tǒng)服務(wù)價(jià)值的影響［J］.應(yīng)用生態(tài)學(xué)報(bào)，2013，24（5）：12311236.

［13］?魏佳豪，溫玉玲，龔志軍，等.近30年鄱陽湖濱岸緩沖帶土地利用變化及生態(tài)系統(tǒng)服務(wù)價(jià)值［J］.生態(tài)學(xué)報(bào)，2022，42（22）：92619273.

［14］?何治學(xué)，王磊，羅君，等.宜賓市近20年土地利用變化及其生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空效應(yīng)［J］.西華師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2023，44（4）：350357.

［15］?謝高地，張彩霞，張昌順，等.中國生態(tài)系統(tǒng)服務(wù)的價(jià)值［J］.資源科學(xué)，2015，37（9）：17401746.

［16］?盧辭，李明鴻.黃河三角洲城市化與生態(tài)系統(tǒng)服務(wù)價(jià)值的耦合關(guān)系［J］.國土與自然資源研究，2023（6）：2026.

［17］?高雯娜，吳玉琪，吳可欣，等.基于生態(tài)系統(tǒng)服務(wù)價(jià)值的土地利用結(jié)構(gòu)優(yōu)化：以清流縣為例［J］.赤峰學(xué)院學(xué)報(bào)（自然科學(xué)版），2023，39（10）：6266.

［18］?王兆峰，楊釗穎，張海燕.基于生態(tài)系統(tǒng)服務(wù)價(jià)值的洞庭湖生態(tài)經(jīng)濟(jì)區(qū)生態(tài)補(bǔ)償標(biāo)準(zhǔn)與策略優(yōu)化［J］.湖南師范大學(xué)自然科學(xué)學(xué)報(bào)，2023，46（5）：6371.

［19］?周夏青.海岸帶土地生態(tài)系統(tǒng)服務(wù)價(jià)值評價(jià)及治理保護(hù)對策研究：以威海市為例［J］.山東國土資源，2023，39（9）：7278.

［20］?趙穎，王飛.資源型省份土地利用變化及其對生態(tài)系統(tǒng)服務(wù)價(jià)值的影響［J］.中國農(nóng)學(xué)通報(bào)，2015，31（17）：127133.

［21］?楊光宗，呂凱，李峰.基于格網(wǎng)尺度的南昌市土地利用變化及生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空相關(guān)性分析［J］.中國土地科學(xué)，2022，36（8）：121130.

［22］?劉金勇，孔繁花，尹海偉，等.濟(jì)南市土地利用變化及其對生態(tài)系統(tǒng)服務(wù)價(jià)值的影響［J］.應(yīng)用生態(tài)學(xué)報(bào)，2013，24（5）：12311236.

［23］?王敏，張聰，惠洪寬.濟(jì)南市的土地利用變化及其生態(tài)服務(wù)價(jià)值研究［J］.齊魯師范學(xué)院學(xué)報(bào)，2017，32（4）：108114.

［24］?代靜，王小燕，杜文強(qiáng).濟(jì)南市土地利用變化及生態(tài)服務(wù)價(jià)值時(shí)空演變分析［J］.環(huán)境保護(hù)科學(xué)，2019，45（2）：615.

［25］?高登.大都市區(qū)土地利用變化及其生態(tài)系統(tǒng)服務(wù)價(jià)值響應(yīng)：以濟(jì)南市為例［J］.中國資源綜合利用，2021，39（12）：7274.

［26］?呂樂婷，張杰，彭秋志，等.東江流域景觀格局演變分析及變化預(yù)測［J］.生態(tài)學(xué)報(bào)，2019，39（18）：68506859.

［27］?李薇，談明洪.太行山區(qū)不同坡度NDVI變化趨勢差異分析［J］.中國生態(tài)農(nóng)業(yè)學(xué)報(bào)，2017，25（4）：509519.

［28］?翟浩然，唐新明，王光輝，等.中國超大城市建設(shè)用地遙感監(jiān)測與時(shí)空格局［J］.遙感信息，2019，34（4）：3947.

［29］?謝高地，甄霖，魯春霞，等.一個(gè)基于專家知識的生態(tài)系統(tǒng)服務(wù)價(jià)值化方法［J］.自然資源學(xué)報(bào)，2008（5）：911919.

［30］?鄧楚雄，鐘小龍，謝炳庚，等.洞庭湖區(qū)土地生態(tài)系統(tǒng)的服務(wù)價(jià)值時(shí)空變化［J］.地理研究，2019，38（4）：844855.

［31］?封建民，郭玲霞，李曉華.漢中市土地利用時(shí)空變化及其對生態(tài)系統(tǒng)服務(wù)價(jià)值的影響［J］.水土保持研究，2020，27（1）：275282.

[STHZ]Influence of Land Use Changes on the Value of Ecosystem Services in Ji'nan from 1980 to 2020

ZHANG Menghong1， JIANG Peng2， ?ZHAO Li1， ZHAO Jing1， LIU Weiguo3， JI Dong1， SHI Hai1， MOU Gang1， WANG Jing1

（1. Rizhao Bureau of Natural Resources and Planning， Shandong， Rizhao 276800， China; 2. No.1 Exploration Brigade of Shandong Coalfield Geological Bureau， Shandong Qingdao 266520， China; 3. Rizhao Central Vitality Zone Development Service Center Shandong Rizhao 276800， China）

Abstract：Faced the problem of resource shortage， unreasonable land utilization types and degradation of ecological systems， exploring the impact of land use change on the value of ecosystem services is of great significance for formulating appropriate policies and measures on territorial spatial planning， natural resource management and ecological civilization construction. Selecting Ji'nan city as the study area， the change of land use structure and transfer among land use types from 1980 to 2020 has been studied， the ecosystem service value has been calculated， and ecological contribution rate of different land use types to the change of ecosystem service value has been analyzed. It is showed that： ①The construction land area shows a trend of continuous increase within the whole period. The area of cultivated land， forest land， grassland and water area show a fluctuating change from 1980 to 2000， but are in a continuous decline between the year 2000 and 2020. ②From 1980 to 2020， the total service value of ecosystem show a general trend of down-up-down. To be specific， four main service values， such as supply， regulation， support and culture all decrease to varying degrees， and reach the largest degree from 2005 to 2010. ③Service value change of ecosystem in Ji'nan city is mostly affected by the cultivated land， which reduces tremendously in the study area due to the urban development. It will cause a significant reduction in arable land area， followed by significant impacts on forests， grasslands， and water bodies.

Key words：Land use change; ecosystem service value; ecological contribution rate; Ji'nan city

收稿日期：20231021；修訂日期：20231202；編輯：王敏作者簡介：張夢虹（1991—），女，山東日照人，工程師，主要從事測繪地理信息相關(guān)工作；Email：466370172@qq.com?*通訊作者：姜鵬（1990—），男，山東日照人，工程師，主要從事測繪地理信息工作；Email：penglove_1224@163.com