基于MCCA的多情景土地利用演變及其ESV變化研究

李楚瑜 徐偉銘 王開家 楊慧 邵爾輝

摘 要：為了探究多情景下土地利用演變與生態(tài)系統(tǒng)服務價值（ecosystem service value，ESV）的特征及兩者間的差異性，以龍巖市為例，基于2005、2010、2015和2020年的4期土地利用數(shù)據(jù)，運用混合元胞自動機（mixed-cell cellular automata，MCCA）模型模擬2035年龍巖市土地利用格局，并采用當量因子法計算分析了歷年以來和未來情景下ESV。結果顯示：2005—2020年龍巖市林地和耕地面積總體上呈下降趨勢，分別減少了10 590.33 hm2、4 183.63 hm2，水域面積小幅度上升，建設用地面積增長迅速；龍巖市ESV在2005—2020年整體上呈現(xiàn)出下降態(tài)勢，累計減少4.15億元；在城市發(fā)展和耕地保護情景下，2035年龍巖市ESV較2020年分別損失了6.52億元和4.30億元；而生態(tài)保護情景下，ESV預測為927.43億元，較2020年上升了0.67億元，這是由于水域面積的增加使ESV得到了有效提高。因此，生態(tài)保護情景是符合未來發(fā)展趨勢的理想土地利用空間格局。

關鍵詞：土地利用演變；MCCA模型；ESV；龍巖市

中圖分類號：F301.24

文獻標志碼：A

土地利用變化是全球環(huán)境和可持續(xù)發(fā)展的研究重點^［1］。近年來，隨著我國城市化與工業(yè)化進程加快，各類城市活動對土地的需求也越來越大^［2］，導致土地利用發(fā)生了劇烈變化，形成了資源掠奪和生態(tài)環(huán)境破壞的粗放式土地利用模式，國土空間格局受到嚴重影響^［3-4］。由于土地利用的矛盾日益加劇，對水循環(huán)、土壤更新等生態(tài)過程^［5-6］產生了一定影響，其通過生態(tài)系統(tǒng)服務價值（ecosystem service value，ESV）來體現(xiàn)和表征^［7-8］。因此，在當前城市化的過程中，研究土地利用演變及其生態(tài)系統(tǒng)服務價值的變化特征，對于優(yōu)化國土空間格局，促進地區(qū)可持續(xù)發(fā)展具有重要意義^［9-10］。

ESV最早由 COSTANZA等^［11］提出，謝高地等^［12］在此基礎上進一步提出了適合中國ESV的計算方法。此后，我國學者開展了系統(tǒng)性研究工作。戴文遠等^［13］基于“三生空間”對福州新區(qū)土地利用功能轉型及生態(tài)服務價值進行研究，彭文甫等^［14］基于1996—2016年的土地利用數(shù)據(jù)探究了對生態(tài)價值服務的影響，王航等^［15］開展土地利用及景觀演化與ESV的相關性研究。總體來看，以上的研究從歷史階段的土地利用演化和ESV的變化特點出發(fā)，但ESV的測算是動態(tài)的，對不同歷史階段生態(tài)系統(tǒng)服務價值的時空演化特征進行分析，存在著一定的滯后^［16］，很難對未來的城市生態(tài)系統(tǒng)服務狀況做出精確的分析。雖然已有部分研究模擬分析了未來土地利用格局^［17］及ESV動態(tài)特征^［18-19］，但在未來不同情景下，針對政策驅動等方面引發(fā)的ESV空間差異性的研究仍然較少，當前國家致力于推進美麗中國生態(tài)文明建設，提出了一系列生態(tài)保護政策，在目前的政策導向下，未來不同發(fā)展情景下的ESV的演變特點迫切需要進行更深入的研究。此外，目前預測未來土地利用變化模型中，TOBLER^［20］最早將元胞自動機（cellular automata，CA）模型應用于地理建模之中，在此之后出現(xiàn)了CA-Markov^［21］、CLUE-S^［22］、 FLUS^［23］和PLUS^［24］等模型。以上模型屬于傳統(tǒng)的CA模型^［25］，均獲得了預期的模擬結果，但也存在以下局限性：元胞內用地類型不一定具有單一性，傳統(tǒng)CA模擬時的元胞為純元胞，即每個元胞內只考慮了一種土地利用類型，忽略了實際元胞內多種土地利用類型混合出現(xiàn)的現(xiàn)象^［26］；其次，大多數(shù)模型采用定性的轉換規(guī)則開展研究區(qū)的模擬工作，不適用于混合元胞的轉換規(guī)則，忽略了亞元胞尺度上的定量轉換規(guī)則。而LIANG等^［27］提出的混合元胞自動機（mixed-cell cellular automata，MCCA）模型考慮了土地利用模式的復雜性，設計了一種與傳統(tǒng)CA模型在元胞、元胞空間、元胞狀態(tài)和轉換規(guī)則的元胞鄰域方面有較大改進的模型，該模型彌補了其他模型元胞內單一土地結構的缺陷，實現(xiàn)了在亞元胞尺度上由定性模擬到定量模擬的跨越。因此，MCCA模型能更精細地模擬未來情景下土地利用演變趨勢。綜上，探究未來不同情景土地利用與ESV演變趨勢，為區(qū)域生態(tài)環(huán)境保護與資源優(yōu)化管理提供科學理論基礎，對完善生態(tài)環(huán)境服務的理論體系具有重要的現(xiàn)實指導作用。

作為閩西部快速城鎮(zhèn)化區(qū)域，龍巖市用地結構變化迅速，同時隨著革命老區(qū)現(xiàn)代化發(fā)展的推進，導致用地結構失衡和功能紊亂，使得人地之間的矛盾愈發(fā)突出。本文基于生態(tài)型現(xiàn)代化幸福龍巖的愿景和新時代支持革命老區(qū)發(fā)展的目標，結合考慮混合元胞特點的MCCA模型，通過對龍巖市多情景下土地利用與ESV的動態(tài)變化特點進行探討，旨在為龍巖市城鎮(zhèn)用地布局科學化、城鎮(zhèn)可持續(xù)發(fā)展提供參考。

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

1.1 研究區(qū)概況

龍巖市位于福建省西部，地處閩粵贛三省交界，跨新羅區(qū)、永定區(qū)、長汀縣、上杭縣、連城縣、武平縣、漳平市四縣兩區(qū)一市，界于115°51′E ～117°45′E， 24°23′N～26°02′N 之間，主要為亞熱帶季風氣候，歷年平均降水量為1 442.1 mm，擁有2個國家級自然保護區(qū)，3個國家森林公園，生物資源豐富，為典型的山地城市。龍巖市是先進制造業(yè)基地和全國重要的紅色城市，生態(tài)文明建設基礎較好，有著“三山為屏、四水串城”的總體格局，作為福建省打造21世紀海上絲綢之路核心區(qū)的內陸樞紐之一，可為福建省自貿區(qū)擴區(qū)提供支撐，未來發(fā)展情景廣闊。研究區(qū)位置示意圖如圖1所示。

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

本文所采用的各期土地利用數(shù)據(jù)（2005年、2010年、2015年和2020年）來源于中科院資源環(huán)境科學與數(shù)據(jù)中心，并按照中國土地利用/覆蓋數(shù)據(jù)的三級分類系統(tǒng)將土地利用類型重分類為耕地、草地、林地、水域、建設用地和未利用地六大類，空間分辨率為30 m。DEM數(shù)據(jù)、基礎地理數(shù)據(jù)分別來源于地理空間數(shù)據(jù)云和OSM開放數(shù)據(jù)平臺。生態(tài)保護紅線數(shù)據(jù)與永久基本農田數(shù)據(jù)來源于龍巖市自然資源局。由于各類數(shù)據(jù)來源不一，將坐標系統(tǒng)一為WGS_1984_UTM_Zone_50N。各生態(tài)系統(tǒng)單位面積服務價值所用的糧食產量、播種面積和價格數(shù)據(jù)來源于《龍巖市統(tǒng)計年鑒》和《中國農村統(tǒng)計年鑒》。

2 研究方法

2.1 MCCA模型

本文采用MCCA模型模擬未來土地利用空間分布。其模型結構主要包括2個部分，土地利用成分發(fā)展?jié)摿τ嬎愫蛠喸叨认峦恋乩米兓M的CA 模型^［28］。技術路線如圖2所示。

2.1.1 模型介紹

1）土地利用成分發(fā)展?jié)摿?/p>

首先，基于起始年份和目標年份，通過隨機森林回歸（random forest regression，RFR）算法訓練出元胞內各土地利用成分面積占比的變化量與驅動因素之間的關系；其次，根據(jù)訓練關系預測得到各土地利用類型的開發(fā)潛力。其中，RFR有效克服了多重空間變量之間共線性的相關問題。此外，混合元胞內土地利用成分的面積占比可能增加或減少。由于混合元胞的土地利用成分是相互約束的，一個成分的增加意味著其他成分的減少，因此我們只需要關注土地利用成分增加的部分，避免重復計算。其表達式為：

R_k=RFR_train （Ys_i，k，F(xiàn)s_i）（1）

P_i，k=RFR_predict（R_k，F(xiàn)_i）（2）

式中：i為混合元胞；k為土地利用類型；s為樣本閾值；R_k為驅動因子與混合元胞內地類覆蓋比例變化的關系；RFR_train（）為RFR的訓練過程；Ys_i，k為混合元胞的樣本數(shù)據(jù)集；Fs_i為各驅動因素的樣本數(shù)據(jù)集；P_i，k為各土地利用成分的開發(fā)潛力；RFR_predict（）為RFR預測過程；F_i為驅動因素數(shù)據(jù)集。

2）亞元胞尺度下的土地利用變化模擬的CA模型

亞元胞，指基于混合元胞內的土地利用類型占比進行模擬，即每個元胞可能存在多種地類。MCCA模型的土地利用變化模擬中，各混合元胞中土地利用成分的面積比例是否增加，以及增加的量，在獲得土地利用成分發(fā)展?jié)摿?、用地需求以及自適應慣性系數(shù)后，由一組數(shù)量轉換規(guī)則估算轉為增加的土地利用成分和其他用地成分改變的數(shù)量，其競爭機理與FLUS模型類似。該輪盤競爭機制中，每個扇區(qū)代表一種地類，各種地類通過輪盤與其他地類競爭，每個元胞內土地利用類型的個數(shù)為競爭的輪數(shù)，每一輪與一個地類競爭，若土地利用類型k獲勝，則該地類數(shù)量會增加，其他地類相應減少。土地利用類型k的總變化概率可以表示為

Ot_i，k=P_i，k×Ωt_i，k×Dt_k（3）

式中：t為迭代次數(shù)；Ot_i，k為土地利用成分的總變化概率；Ωt_i，k為混合元胞的鄰域效應；Dt_k為對土地利用類型k的未來需求的反饋，是自適應慣性系數(shù)。當模擬的土地利用成分數(shù)值等于目標數(shù)值時，MCCA將輸出k層模擬結果，每層為一個地類的空間分布。

經(jīng)多次實驗，MCCA模型的回歸樹數(shù)目設為100，采樣率設為5%，鄰域為3，步長為1。

2.1.2 情景設定

圍繞“建設閩西南生態(tài)型現(xiàn)代化城市，打造有溫度的幸福龍巖”的建設愿景，以龍巖市2005—2020 年土地利用變化特征及地區(qū)未來發(fā)展規(guī)劃為依據(jù)，根據(jù)城市的實際情況，設置城市發(fā)展、耕地保護、生態(tài)保護3個情景，模擬未來土地利用需求，如表1所示。

2.1.3 轉換成本矩陣

在模擬未來多情景下的土地利用需求時，根據(jù)設計的情景，設置不同地類轉換為其他地類的難度，見表2。其中：a為耕地；b為林地；c為草地；d為水域；e為建設用地；f為未利用地；1表示地類之間可以轉換，0表示地類之間禁止轉換。

2.2 ESV評估

1）基于糧食價格的修正

據(jù)文獻［29］的研究結果可知，1個當量因子的價值為全國平均糧食單產市場價格的1/7。考慮到糧食價格的波動性，參照2020年全國主要糧食平均價格和2005—2020年龍巖市平均糧食產量，修正方法^［30］如下：

式中：V為龍巖市單位面積農田糧食作物的市場價值，元/hm2；P為全國主要糧食作物平均單價，元/kg；Q為研究區(qū)糧食作物單位面積產量，kg/hm2。

2）ESV系數(shù)的修正

采用2015年修正的符合中國等發(fā)展中國家實際的ESV當量因子表^［12］，并依據(jù)龍巖市的生態(tài)系統(tǒng)服務類型，制定了單位面積ESV當量，建設用地設為0。單位面積ESV系數(shù)為

V_i=C_iV i=1，2，…，n（5）

式中：V_i為第i類土地的ESV系數(shù)，元/hm2；C_i為第i類土地單位面積生態(tài)系統(tǒng)價值當量。

對單位面積ESV系數(shù)進行修正，結果如表3所示。

3）ESV計算

利用ESV分析模型計算龍巖市各時期的ESV（V_ES），計算公式^［31］為

V_ES=∑A_iV_i（6）

式中，A_i為第i類土地的面積，hm2。

3 結果與分析

3.1 龍巖市土地利用及其ESV變化

3.1.1 土地利用演變分析

通過4期土地利用數(shù)據(jù)，分析研究區(qū)2005—2020年間土地利用變化數(shù)據(jù)，如表4和圖3所示。龍巖市以林地為主，水域主要位于其西南方向，耕地圍繞建設用地分布。從表4和圖3可以看出：龍巖市的林地占據(jù)約73%的土地面積，從2005年的73.62%至2020年的73.06%，林地總體上呈下降趨勢，共減少了10 590.33 hm2；耕地面積明顯減少，2005—2015年間面積占比呈現(xiàn)下降趨勢，到2020年耕地比例有所緩和并開始上升，面積累計減少了4 183.63 hm2；水域在這15年間增長較緩，面積增加了310.45 hm2，其比例在2010年上升后，直到2020年無明顯變化，這也證明了龍巖市一直致力于保護好四大流域的重大舉措取得了成效；建設用地2005—2020年的面積持續(xù)增長，共計增加了 12 080.34 hm2，其中2005—2010年間增長面積最大，表明龍巖市的建設用地在這一時期的擴張速度是最快的；草地總體呈增長趨勢，2005—2010年間增長面積相對較多，僅次于建設用地。未利用地占比小，總體上面積變化不明顯，減少約1.83%。綜上，林地與耕地面積的下降主要是源于建設用地的擠壓和侵占。

3.1.2 ESV變化分析

根據(jù)式（6）可得到龍巖市2005—2020年ESV的變化，如表5所示。從表5可知：龍巖市2005、2010、2015、2020年ESV分別為930.91、927.77、927.95、926.76億元，15年間ESV共計減少了4.15億元，總體呈下降趨勢；在ESV總量中，林地對其貢獻最大，約占價值總量的83%，這與龍巖市高森林覆蓋率有關；15年間林地的ESV明顯降低，共計減少了5.84億元，而草地ESV的增加量最大，共計增加了1.15億元?？傮w來看，2005—2010年龍巖市的ESV下降較快，2010—2015年小幅度上升，但2015—2020年又呈下降趨勢。

3.2 未來不同情景下土地利用及其ESV變化

3.2.1 不同情景下土地利用格局變化

1）精度檢驗

LIANG等^［27］和蔣小芳等^［28］分別采用空間分辨率在250～1 250 m和250～500 m的混合像元驗證MCCA模型的模擬精度，證明當空間分辨率為250 m時，土地利用模擬精度最高。本文以龍巖市為研究區(qū)，較上述研究范圍更大，鑒于此，采用250 m的空間分辨率模擬未來結果圖。本文基于2005年土地利用柵格圖模擬2020年的土地利用圖，采用OA、Kappa系數(shù)和mcFoM評估MCCA模型的預測精度，其中mcFoM與FoM類似，模擬結果愈接近1愈準確^［32］。采用2020年各地類的實際結果和MCCA模型計算的模擬結果，分別與2005年土地利用數(shù)據(jù)進行對比，得到15年間土地利用成分變化情況，變化程度為-1～1，如圖4所示。變化程度值為正表示地類的面積在擴張，為負表示減少，為0則表示地類沒有變化。由圖4可以發(fā)現(xiàn)，各地類的變化區(qū)域基本吻合，模擬效果較好。OA、Kappa系數(shù)和mcFoM的計算結果分別為0.95、0.89和0.222，說明模擬結果可靠。

2）2035年不同情景預測結果

MCCA模型對2035年地類結構進行空間可視化，土地利用類型的覆蓋比例為0～1，結果如圖5所示。在城市發(fā)展情景下，相較于2020年，建設用地擴張十分明顯，例如上杭縣古田鎮(zhèn)，在原有基礎上繼續(xù)擴張，成團發(fā)展，表明在不受各項保護政策的影響下，城鎮(zhèn)化發(fā)展迅速。具體來講，建設用地的面積增加了9 189.17 hm2，林地和草地分別減少了5 108.67 hm2和2 109.22 hm2，說明林地與草地受到建設用地的侵占較為突出。

耕地保護情景下，耕地較2020年增加了331.15 hm2，建設用地的面積增加了6 868.23 hm2。耕地在基本農田保護紅線的政策下保障了糧食安全，得到了有效的利用和保護，在長汀北部、連城北部和上杭南部分布相對集中。隨著區(qū)域經(jīng)濟的不斷發(fā)展，城市建設用地持續(xù)擴張，城鎮(zhèn)化進程必然會導致其他地類的萎縮，僅有耕地得到了有效的控制。

生態(tài)保護情景下，除了草地和未利用地，其他用地面積都有所增長，是符合未來發(fā)展趨勢的理想土地利用空間格局。在此情景下，建設用地、林地、水域、耕地面積分別增長了1 456.49、602.47、360.14、131.66 hm2，而草地和未利用地分別下降了1 587.65 hm2、963.11 hm2，說明草地和未利用地逐漸轉化為建設用地、林地、水域以及耕地。另外，林地和耕地被建設用地侵占的程度也有所降低，在一定程度上，有利于維護地區(qū)的生態(tài)安全，有利于地區(qū)的可持續(xù)發(fā)展。

3.2.2 不同情景下生態(tài)服務價值變化

2035年，在城市發(fā)展、耕地保護和生態(tài)保護情景下，ESV分別為920.24億元、922.46億元和927.43億元，如圖6所示。相比于2020年，龍巖市ESV只有在生態(tài)保護情景下提高了0.67億元，在城市發(fā)展情景和耕地保護情景下均呈現(xiàn)減少趨勢，分別損失了6.52億元和4.30億元。在城市發(fā)展情景下，生態(tài)系統(tǒng)的服務價值損失最大。由于經(jīng)濟的快速發(fā)展和城鎮(zhèn)化的逐步推進，需要占用耕地和林地來擴大建設用地的面積，而建設用地的增多導致生態(tài)系統(tǒng)的服務價值降低，這對龍巖市今后的發(fā)展具有重要的參考價值，必須適時調整發(fā)展戰(zhàn)略，優(yōu)化國土空間布局。在耕地保護情景下，ESV比城市發(fā)展情景多2.22億元，說明永久基本農田的劃定有效地控制了耕地的數(shù)量，一定程度上遏制了建設用地侵占耕地。在生態(tài)保護情景下，與城市發(fā)展和耕地保護情景相比，ESV分別增加7.19億元和4.97億元，說明林地作為龍巖市最大的用地類型，對生態(tài)保護起著非常重要的作用，一定程度上也延緩了建設用地的擴張。在政策的干預下，生態(tài)保護紅線對生態(tài)的發(fā)展和構建幸福龍巖有至關重要的作用。

3.2.3 不同情景下土地利用變化對生態(tài)服務價值影響差異

圖7是2020—2035年龍巖市各土地利用類型ESV變化。從圖7可以看出：引起3種情景內ESV變化的主要因素是水域面積的變化。其中，在城市發(fā)展情景下，水域的ESV減少了8.5%，變動幅度較大，且耕地保護情景下ESV的損失次之，減少4.5%。林地和草地面積的減少成為城市發(fā)展和耕地保護情景中ESV降低的重要因素。而在生態(tài)保護情景下，水域面積的增加使得ESV得到了有效提升，提高了1.08億元，一定程度上削弱了草地ESV降低的趨勢。此外，由于林地、耕地和水域轉換為建設用地的概率降低和生態(tài)保護紅線的限制，使得林地、耕地面積增加帶來的價值增長與草地價值減少表現(xiàn)出互補趨勢。相較于其他2種情景，草地在生態(tài)保護情景下ESV的減少量最小，且林地和耕地ESV呈正向增長（由于林地和耕地的基數(shù)較大，ESV上升幅度不明顯），表明在此情景下，生態(tài)環(huán)境得到改善，生態(tài)結構日趨合理化且生態(tài)系統(tǒng)合理化水平更高。

4 結論與討論

本文基于2005—2020年龍巖市土地利用與ESV時空演變特征，利用MCCA模型模擬了不同情景下2035年土地利用格局及其ESV變化差異。主要結論如下：

1）隨著城鎮(zhèn)化的高速發(fā)展，2005—2020年龍巖市林地和耕地面積總體上呈下降趨勢，分別減少了10 590.33 hm2和4 183.63 hm2，水域面積小幅度上升，建設用地面積則增長迅速，這是由于建設用地的無序擴張，對林地與耕地造成了侵占和擠壓。

2）龍巖市的生態(tài)服務價值在2005—2020年整體上呈現(xiàn)出下降的態(tài)勢，累計減少4.15億元。林地的ESV減少量最大，減少了5.84億元，而草地和水域的ESV則一直處于上升趨勢，分別增加了1.15億元、0.93億元。由此可知：城鎮(zhèn)化的迅速發(fā)展引起的土地利用改變是影響龍巖市ESV變化的重要因素；建設用地的迅速擴張，會使林地和耕地的生態(tài)結構發(fā)生變化，從而削弱了ESV。

3）在城市發(fā)展情景和耕地保護情景下，龍巖市2035年的ESV較2020年分別損失了6.52億元和4.30億元，林地和草地面積的減少成為城市發(fā)展和耕地保護情景中ESV降低的重要因素。而生態(tài)保護情景下，ESV預測為927.43億元，較2020年上升了0.67億元，這是由于水域面積的增加使得ESV得到了有效提高。

本文在分析龍巖市歷史土地利用變化規(guī)律基礎上，借助MCCA模型模擬城市發(fā)展、耕地保護與生態(tài)發(fā)展3種情景下龍巖市土地利用演變情況，進而利用當量因子法對3種情景下用地結構的生態(tài)系統(tǒng)服務能力進行價值化表達。對龍巖市而言，堅持發(fā)展生態(tài)保護理念是提高ESV的正確方向，該情景不僅推進了土地利用的協(xié)同布局，也對實現(xiàn)區(qū)域的持續(xù)健康發(fā)展起著一定作用。從研究方法上看，基于2020年的糧食價格與生物量進行估算，客觀上可能存在對未來ESV的高估或低估，未來可以嘗試采用更加精確的ESV評估方法。此外，本文僅采用250 m混合像元空間分辨率的數(shù)據(jù)進行模擬，研究尺度較單一，沒有考慮更高或更低的分辨率數(shù)據(jù)進行對比實驗，未來可進一步探究多尺度視角下土地利用演變對ESV時空演變的影響。

參考文獻：

[1]劉紀遠， 張增祥， 徐新良， 等. 21世紀初中國土地利用變化的空間格局與驅動力分析［J］. 地理學報， 2009， 64（12）： 1411-1420.

［2］ 樊杰. 主體功能區(qū)戰(zhàn)略與優(yōu)化國土空間開發(fā)格局［J］. 中國科學院院刊， 2013， 28（2）： 193-206.

［3］ 白中科， 周偉， 王金滿， 等. 試論國土空間整體保護、系統(tǒng)修復與綜合治理［J］. 中國土地科學， 2019， 33（2）： 1-11.

［4］ 黃金川， 林浩曦， 漆瀟瀟. 面向國土空間優(yōu)化的三生空間研究進展［J］. 地理科學進展， 2017， 36（3）： 378-391.

［5］ 王永琪， 馬姜明. 基于縣域尺度珠江—西江經(jīng)濟帶廣西段土地利用變化對生態(tài)系統(tǒng)服務價值的影響研究［J］. 生態(tài)學報， 2020， 40（21）： 7826-7839.

［6］ 雷璇， 楊波， 蔣衛(wèi)國， 等. 東洞庭濕地植被格局變化及其影響因素［J］. 地理研究， 2012， 31（3）： 461-470.

［7］ 劉瑩， 耿文亮， 邵靜文， 等. “三生空間”視角下土地利用變化與生態(tài)系統(tǒng)服務價值響應： 以黃河下游地區(qū)為例［J］. 地域研究與開發(fā)， 2021， 40（4）： 129-135.

［8］ 楊清可， 段學軍， 王磊， 等. 基于“三生空間”的土地利用轉型與生態(tài)環(huán)境效應： 以長江三角洲核心區(qū)為例［J］. 地理科學， 2018， 38（1）： 97-106.

［9］ 龍花樓. 論土地利用轉型與土地資源管理［J］. 地理研究， 2015， 34（9）： 1607-1618.

［10］劉永強， 廖柳文， 龍花樓， 等. 土地利用轉型的生態(tài)系統(tǒng)服務價值效應分析： 以湖南省為例［J］. 地理研究， 2015， 34（4）： 691-700.

［11］COSTANZA R， D′ARGE R， DE GROOT R， et al. The value of the worlds ecosystem services and natural capital［J］. World Environment， 1997， 387： 253-260．

［12］謝高地， 張彩霞， 張昌順， 等. 中國生態(tài)系統(tǒng)服務的價值［J］. 資源科學， 2015， 37（9）： 1740-1746.

［13］戴文遠， 江方奇， 黃萬里， 等. 基于“三生空間”的土地利用功能轉型及生態(tài)服務價值研究： 以福州新區(qū)為例［J］. 自然資源學報， 2018， 33（12）： 2098-2109.

［14］彭文甫， 周介銘， 楊存建， 等. 基于土地利用變化的四川省生態(tài)系統(tǒng)服務價值研究［J］. 長江流域資源與環(huán)境， 2014， 23（7）： 1053-1062.

［15］王航， 秦奮， 朱筠， 等. 土地利用及景觀格局演變對生態(tài)系統(tǒng)服務價值的影響［J］. 生態(tài)學報， 2017， 37（4）： 1286-1296.

［16］李龍， 吳大放， 王芳， 等. 中國快速城市化區(qū)域生態(tài)系統(tǒng)服務價值預測及權衡研究： 以佛山市為例［J］. 生態(tài)學報， 2020， 40（24）： 9023-9036.

［17］孫瑋健， 張榮群， 艾東， 等. 基于元胞自動機模型的土地利用情景模擬與驅動力分析［J］. 農業(yè)機械學報， 2017， 48（增刊1）： 254-261.

［18］何玲， 賈啟建， 李超， 等. 基于生態(tài)系統(tǒng)服務價值和生態(tài)安全格局的土地利用格局模擬［J］. 農業(yè)工程學報， 2016， 32（3）： 275-284.

［19］歐陽曉， 賀清云， 朱翔. 多情景下模擬城市群土地利用變化對生態(tài)系統(tǒng)服務價值的影響： 以長株潭城市群為例［J］. 經(jīng)濟地理， 2020， 40（1）： 93-102.

［20］TOBLER W R. Cellular Geography［M］. Netherlands： Springer， 1979： 379-386.

［21］崔旺來， 蔡莉， 奚恒輝， 等. 基于土地利用/覆蓋變化的浙江大灣區(qū)生態(tài)安全評價及多情景模擬分析［J］. 生態(tài)學報， 2022， 42（6）： 2136-2148.

［22］趙旭， 湯峰， 張蓬濤， 等. 基于CLUE-S模型的縣域生產-生活-生態(tài)空間沖突動態(tài)模擬及特征分析［J］. 生態(tài)學報， 2019， 39（16）： 5897-5908.

［23］楊浩， 盧新海. 基于“三生空間”演化模擬的村莊類型識別研究： 以湖南省常寧市為例［J］. 中國土地科學， 2020， 34（6）： 18-27.

［24］李安林， 周艷， 唐麗毅， 等. 怒江州土地利用模擬及生態(tài)系統(tǒng)服務價值評估： 基于PLUS模型的多情景分析［J/OL］. 中國農業(yè)資源與區(qū)劃， 2022： 1-12. （2022-04-17）［2022-09-03］. http：//kns.cnki.net/kcms/detail/11.3513.S.20220415.1200.006.html.

［25］戴爾阜， 馬良. 土地變化模型方法綜述［J］. 地理科學進展， 2018， 7（1）： 152-162.

［26］SHI B X， YANG J Y. Scale， distribution， and pattern of mixed land use in central districts： a case study of Nanjing， China［J］. Habitat International， 2015， 46： 166-177.

［27］LIANG X， GUAN Q F， CLARKE K C， et al. Mixed-cell cellular automata： a new approach for simulating the spatio-temporal dynamics of mixed land use structures［J］. Landscape Urban Plan， 2021， 205： 103960.1-103960.17.

［28］蔣小芳， 段翰晨， 廖杰， 等. 基于多模型的黑河中游甘臨高地區(qū)土地利用情景模擬［J］. 農業(yè)機械學報， 2022， 53（9）： 178-188.

［29］謝高地， 魯春霞， 冷允法， 等. 青藏高原生態(tài)資產的價值評估［J］. 自然資源學報， 2003， 18（2）： 189-196.

［30］白羽萍， 王佳斌， 劉戀， 等. 黑河流域生態(tài)系統(tǒng)服務價值預測及權衡分析： 以張掖市為例［J］. 生態(tài)學報， 2022， 42（13）： 5187-5199.

［31］王兆峰， 許靜. 土地利用演變對國家公園生態(tài)系統(tǒng)服務價值的影響： 以三江源國家公園為例［J］. 生態(tài)學報， 2022， 42（17）： 1-11.

［32］張曉瑤， 張瀟， 李冬花， 等. 城市土地利用變化對生態(tài)系統(tǒng)服務價值影響的多情景模擬： 以深圳市為例［J］. 生態(tài)學報， 2022， 42（6）： 2086-2097.

（責任編輯：周曉南）

Multi-Scenario Land Use Eevolution and ESV Change Based on MCCA： A Case Study of Longyan City

LI Chuyu^{1，2，3}， XU Weiming^{*1，2，3}， WANG Kaijia^{1，2，3}， YANG Hui^1，2，3， SHAO Erhui^1，2，3

（1.Institute of Digital China， Fuzhou University （Fujian）， Fujian 350108， China; 2.Key Laboratory of Spatial Data Mining and Information Sharing， Ministry of Education， Fuzhou University， Fujian 350002， China; 3.National and Local Joint Engineering Research Center of Geospatial Information Technology， Fuzhou University， Fujian 350002， China）

Abstract： In order to explore the characteristics of land use evolution and ecosystem service value （ESV） under multiple scenarios and the differences between the two， this paper takes Longyan City as an example and is based on the four phases of land use data in 2005， 2010， 2015 and 2020. The mixed-cell cellular automata （MCCA） model was used to simulate the land use pattern of Longyan in 2035， and the ESV was calculated and analyzed by using the equivalent factor method. The results showed that the area of forest land and cultivated land in Longyan City decreased by 10 590.33 hm2 and 4 183.63 hm2 respectively from 2005 to 2020， the water area increased slightly， and the area of construction land increased rapidly. The ESV of Longyan City showed a downward trend from 2005 to 2020， with a total decrease of 415 million yuan; under the scenario of urban development and cultivated land protection， the ESV of Longyan City in 2035 will lose 652 million yuan and 430 million yuan respectively compared with 2020. Under the ecological protection scenario， the ESV is predicted to be 92.743 billion yuan， an increase of 67 million yuan compared with 2020， which is a result of the increase in water area. Therefore， the ecological protection scenario is an ideal spatial pattern of land use in line with the future development trend.

Key words： land use evolution; MCCA model; ESV（ecosystem service value）; Longyan City