江蘇省綠地生態(tài)系統(tǒng)服務價值時空動態(tài)變化特征

徐可依 馮陽 陳妍 吳榛

摘要：江蘇省作為我國城鎮(zhèn)化快速發(fā)展的區(qū)域之一，土地利用變化劇烈，其綠地空間格局及生態(tài)服務價值受建設活動影響顯著。探明江蘇省綠地空間格局及其生態(tài)服務價值時空變化特征，可為綠地空間景觀格局規(guī)劃與國土空間管控提供科學依據(jù)。以江蘇省為研究區(qū)，以2000—2020年中5年為1期的土地利用數(shù)據(jù)作為基礎，運用土地利用動態(tài)度和利用土地利用轉移矩陣分析研究區(qū)土地利用類型的動態(tài)變化特征；通過對相應景觀格局指數(shù)的研究來反映綠地系統(tǒng)的結構與生態(tài)過程的聯(lián)系；通過生態(tài)系統(tǒng)服務價值評價方法對江蘇省綠地生態(tài)系統(tǒng)的服務價值進行評估，由此得到江蘇省生態(tài)系統(tǒng)服務價值時空動態(tài)變化特征。結果表明，研究期間江蘇省綠地中有林地和高覆蓋度草地為主要優(yōu)勢類型，綠地總面積呈下降趨勢，降幅為20.30%；江蘇省綠地景觀格局指數(shù)變化明顯，綠地景觀分布不均勻，均衡程度被破壞，多樣性下降，總體趨于破碎化；江蘇省綠地ESV總體呈下降趨勢，研究期內共減少15.75%；綠地高ESV地區(qū)主要分布在江蘇省北部和西南部，2020年江蘇省綠地系統(tǒng)ESV最高的城市是南京市；總體下降最多的城市是鹽城市。江蘇省綠地及其生態(tài)系統(tǒng)服務價值時空變化呈現(xiàn)明顯的地域差異。研究結果可為科學規(guī)劃綠地空間，優(yōu)化綠地景觀格局，促進綠地生態(tài)系統(tǒng)服務價值水平提升提供依據(jù)。

關鍵詞：綠地空間；景觀格局；生態(tài)系統(tǒng)服務價值；時空變化特征；江蘇省

中圖分類號：TU985.1? 文獻標志碼：A

文章編號：1002-1302（2024）08-0233-09

收稿日期：2023-07-07

基金項目：江蘇高校“大學素質教育與數(shù)字化課程建設”專項課題（編號：2020JDKT024）；江蘇省社會科學基金（編號：23SHC002）；江蘇省自然科學基金（編號：BK20201359）；江蘇省高等學校自然科學研究面上項目（編號：19KJB220001）。

作者簡介：徐可依（2000—），女，江蘇南京人，碩士研究生，主要從事風景園林規(guī)劃設計研究。E-mail：375636797@qq.com。

通信作者：陳 妍，碩士，講師，主要從事風景園林規(guī)劃設計研究。E-mail：chenyan@njtech.edu.cn。

綠地是指城市中自然或人為被植物覆蓋的地塊和用地，包括林地和草地［1］。本研究的綠地類型按照LUCC分類系統(tǒng)分類。綠地可以幫助改善生態(tài)、保護環(huán)境和美化景觀，是實現(xiàn)城市可持續(xù)發(fā)展的重要組成部分，具有重要的生態(tài)價值。城鎮(zhèn)化進程加快，大量林地、草地被建設用地所侵占，使得綠地系統(tǒng)被破壞，生態(tài)問題面臨重大挑戰(zhàn)［2］，人們更加意識到綠地在生態(tài)系統(tǒng)和城市環(huán)境中的重要性。深入分析區(qū)域綠地景觀格局及綠地生態(tài)服務價值時空動態(tài)演變特征，可為科學保護綠地系統(tǒng)、改善城市環(huán)境、促進環(huán)境的可持續(xù)發(fā)展提供有效支撐［3］。生態(tài)系統(tǒng)服務價值是指根據(jù)自然生態(tài)系統(tǒng)的功能、結構和過程提供的產品和服務為人類所帶來的價值。生態(tài)系統(tǒng)服務價值的評估是生態(tài)功能區(qū)劃、環(huán)境保護等的重要基礎和依據(jù)，是全球經濟的重要組成部分［4-6］，是景觀生態(tài)學、風景園林學以及生態(tài)經濟學等學科研究的重點問題之一［7］?？茖W評估生態(tài)系統(tǒng)服務價值，可為生態(tài)環(huán)境保護和國土空間規(guī)劃提供參考。隨著人們對綠地價值的日漸重視，綠地生態(tài)系統(tǒng)服務價值的研究逐漸成為各學科當前研究的焦點之一。自Costanza在1997年發(fā)文對全球綠地系統(tǒng)服務價值總量進行評估之后，國際上相繼開展了有關研究。Semenyuk等對不同區(qū)域的森林生態(tài)系統(tǒng)服務價值產生的環(huán)境影響和變化進行了評估和研究［8-9］；Admasu等分析了土地利用及覆被變化對不同流域生態(tài)系統(tǒng)服務價值的影響［10-11］。Najmuddin等對不同國家影響城市生態(tài)系統(tǒng)服務價值的相關土地變化進行了探究［12-13］。我國對于生態(tài)系統(tǒng)服務價值的相關研究起步較晚，但也取得了很多進展。鄒麗麗等采用不同的研究方法對不同地區(qū)的灘涂和濕地的生態(tài)系統(tǒng)服務價值進行了分析［14-15］。馬偉波等分析了不同區(qū)域城市群的相關生態(tài)系統(tǒng)服務價值的時空演變以及驅動［16-17］。謝高地等對我國生態(tài)系統(tǒng)服務價值進行評估并分析了我國生態(tài)系統(tǒng)服務價值的發(fā)展趨勢［18-19］。但是，目前我國對于省域等大尺度及其各行政區(qū)的綠地系統(tǒng)生態(tài)服務價值時空動態(tài)變化特征的研究還比較少。因此，本研究對江蘇省綠地生態(tài)系統(tǒng)服務價值的空間變化進行分析，旨在為生態(tài)系統(tǒng)保護和城鎮(zhèn)管理、規(guī)劃及建設提供科學依據(jù)。江蘇省內多個地級市為國家園林城市，但省域范圍內綠地空間分布不均衡。近幾十年來，城鎮(zhèn)化快速發(fā)展導致綠地空間被動受到城市建設用地擠壓，生態(tài)系統(tǒng)被破壞，需要深入分析江蘇省綠地生態(tài)服務價值時空動態(tài)演變特征，以期為國土空間規(guī)劃和生態(tài)空間保護提供科學依據(jù)。本研究以江蘇省為例，運用土地利用動態(tài)度、土地利用轉移矩陣、景觀格局指數(shù)與生態(tài)服務價值評價等方法，通過對2000—2020年江蘇省綠地景觀格局及生態(tài)系統(tǒng)服務價值時空動態(tài)演變特征的研究和分析，發(fā)現(xiàn)江蘇省綠地ESV總體呈下降趨勢，綠地景觀分布不均勻，均衡程度被破壞，多樣性下降，總體趨于破碎化；通過探討整個省域與局部市域的綠地生態(tài)系統(tǒng)服務價值空間分異及變化趨勢，為江蘇省省域及各地市間的綠地景觀格局優(yōu)化、生態(tài)空間保護、可持續(xù)發(fā)展和相關規(guī)劃提供科學依據(jù)和理論指導。

1 研究地區(qū)與研究方法

1.1 研究區(qū)概況

江蘇省位于30°45′～35°20′N、116°18′～121°57′E（圖1），省域面積約為10.72萬km2，與山東省、安徽省、浙江省和上海市接壤。江蘇省地勢平坦，以平原和低山丘陵為主，屬于亞熱帶季風氣候，年均降水量約為1 200 mm，雨量充沛，年均氣溫為14～23? ℃，四季分明。江蘇省內湖泊、綠地眾多，內含多個國家園林城市，其中水域面積占省域總面積的11%以上，全省林地覆蓋率約為26.4%，建成區(qū)的綠化覆蓋率約為45%，人均公共綠地面積為13.7 m2。近年來，由于江蘇省處于優(yōu)越的地理區(qū)位，區(qū)域城鎮(zhèn)化快速發(fā)展，造成人地矛盾問題日益突出，劇烈的人類活動致使區(qū)域綠地生態(tài)產生明顯變化。

1.2 數(shù)據(jù)來源與預處理

本研究基礎數(shù)據(jù)選取江蘇省2000年、2005年、2010年、2015年、2020年的土地利用分類數(shù)據(jù)，來源于中國科學院資源環(huán)境科學與數(shù)據(jù)中心（www.resdc.cn/）。本研究參照現(xiàn)有的LUCC（land-use and land-cover change）分類體系，運用ArcGIS軟件將江蘇省的土地利用類型劃分成耕地、林地、草地、建設用地、水域、未利用地6個一級大類。同時，根據(jù)研究所需將林地和草地這2類綠地進行二次分類，其中林地劃分為有林地、疏林地、灌木林、其他林地4個二級子類，草地劃分為高覆蓋度草地、中覆蓋度草地、低覆蓋度草地3個二級子類（圖2）。

1.3 研究方法

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

土地利用動態(tài)度反映某種土地類型在一段研究期內的面積變化情況［20］。本研究利用土地利用動態(tài)度來反映各綠地類型面積的動態(tài)變化情況［21-22］。計算公式為

K=Sb-SaSa×1T×100%。

式中：K表示某類土地類型的動態(tài)度；Sa表示某類土地類型研究初期的面積；Sb表示某類土地類型研究末期的面積；T表示總研究時長。

1.3.2 土地利用轉移矩陣

土地利用轉移矩陣是通過對某研究區(qū)不同時期土地利用狀況的變化分析求得的一個二維矩陣。利用土地利用轉移矩陣可直觀揭示各類型土地利用的移動方向以及變化數(shù)量，有利于分析研究區(qū)土地利用類型的動態(tài)變化特征［23］。本研究利用土地轉移矩陣分析2000—2020年江蘇省綠地變化狀況及特征，對江蘇省綠地景觀類型與其他土地利用類型間相互轉換的情況和趨勢進行分析。

1.3.3 景觀格局指數(shù)

景觀格局是指大小、形狀不同的景觀要素在景觀空間上的分布，其影響并決定著各種生態(tài)過程。景觀格局指數(shù)是分析生態(tài)過程變化的常用指標，這種指數(shù)可以較直觀與客觀地數(shù)量化景觀格局，便于對景觀格局進行描述，可以更好地解釋景觀功能與景觀結構的相關聯(lián)系［24］。本研究通過對相應景觀格局指數(shù)的分析來反映江蘇省綠地系統(tǒng)的結構與生態(tài)過程的聯(lián)系，選取斑塊數(shù)量（NP）、形態(tài)指數(shù)（LSI）、聚集度（AI）、斑塊密度（PD）、香農多樣性指數(shù)（SHDI）、香農均勻度指數(shù)（SHEI）和蔓延度（CONTAG）等7個綠地景觀格局指標，并運用Fragstats 4.2軟件進行相關景觀格局指數(shù)計算［25-27］。

1.3.4 綠地生態(tài)系統(tǒng)服務價值評價

生態(tài)系統(tǒng)服務價值評價是指對生態(tài)系統(tǒng)中所含的不同生態(tài)資源的價值進行評估。學者把生態(tài)系統(tǒng)服務價值的評估方法分為替代市場法和假想市場法，前者是指利用生態(tài)服務功能的相關產品的市場價格進行評估，后者是指利用旅行費、防護費等間接價值進行假想評估。本研究采用替代市場法對江蘇省綠地生態(tài)系統(tǒng)服務價值進行評估。

本研究根據(jù)江蘇省土地利用以及研究目的，按照LUCC分類把各綠地分為有林地、灌木林、疏林地、其他林地、高覆蓋度草地、中覆蓋度草地以及低覆蓋度草地。將各類型綠地的單位面積ESV與謝高地等的基礎當量表中綠地系統(tǒng)二級分類對應［28］，分別對應針闊混交、灌木、闊葉、針葉、灌草叢、草甸以及草原。測算出的江蘇省綠地生態(tài)系統(tǒng)服務價值各綠地單位面積價值見表1。綠地生態(tài)系統(tǒng)服務價值估算公式如下［29］。

ESV=∑pi×Li。

式中：ESV表示研究區(qū)域綠地生態(tài)系統(tǒng)服務價值總量，元/年；pi表示第i類綠地的生態(tài)系統(tǒng)服務單位面積價值，元/（hm2·年）；Li表示研究區(qū)域內第i類綠地面積，hm2。本研究將ESV劃分為極低（ESV<1.5萬元/年）、低（1.5萬元/年≤ESV＜10.0萬元/年）、中（10.0萬元/年≤ESV＜60.0萬元/年）、高（60.0萬元/年≤ESV＜400.0萬元/年）、極高（ESV≥400.0萬元/年）等5個等級。

2 ?結果與分析

2.1 江蘇省土地利用變化

2000—2020年江蘇省土地利用類型中耕地占比最大，主要分布在蘇北、蘇中平原和濱海地區(qū)，但耕地面積逐年呈減小趨勢，總降幅達11.77%。建設用地在2000—2010年呈快速增加趨勢，由 1 465 484.40 hm2 增加到2 003 813.91 hm2，增幅為 36.73%；2010—2020年增速逐漸放緩，由 2 003 813.91 hm2 增加到2 149 048.26 hm2，增幅為7.25%。林地面積在2000—2005年小幅增加，由337 542.30 hm2增加到337 934.97 hm2；但在2005—2010年快速減小，由337 934.97 hm2減小到312 680.88 hm2，降幅為7.47%；2010—2020年緩慢減小，由312 680.88 hm2減小到306 821.25 hm2，減小了1.87%。草地在2005—2010年快速減少，由130 641.75 hm2減少到80 607.69 hm2，減幅高達38.30%；但在2010—2015年出現(xiàn)增加，由 80 607.69 hm2 增加到 92 871.18 hm2，增加15.21%。林地和草地面積在2005—2010年的快速減少是由于此期間城鎮(zhèn)屬于快速發(fā)展階段，這一階段綠地被大量建設用地和耕地取代，但在2010—2015年政府為防止建設用地盲目擴張破壞生態(tài)環(huán)境，深化了綠地保護，但在多重影響因素下，耕地、林地、草地仍在不同程度地減少。水域面積除2010—2015年小幅減小外，一直在穩(wěn)定增加，為保護生態(tài)環(huán)境，江蘇省在2017年出臺了《江蘇省生態(tài)河湖行動計劃（2017—2020年）》，加強了水資源保護和水生態(tài)修復，從而使水域面積持續(xù)增加（圖3）。

由圖4可知，江蘇省綠地系統(tǒng)中有林地面積占綠地比例最大，其次為高覆蓋度草地，中覆蓋度草地和低覆蓋度草地占綠地比例較小。各類型綠地整體呈減少趨勢。其中，各類型林地在2005—2010年減少較快，分別減少1.38%、15.75%、10.92%、32.92%。各類型草地也在2005—2010年快速減少，其中高覆蓋度草地變化最明顯，減少了 49 914.90 hm2，減幅為38.79%；但在2010—2015年高覆蓋度草地面積又明顯增加，增加 12 262.41 hm2，增幅達15.57%；到2015—2020年又恢復下降趨勢。由圖4還可知，2000—2020年有林地為江蘇省綠地中的優(yōu)勢類型，其次為高覆蓋度草地，2000—2020年林地占綠地的比例從70.79%增加到77.43%，有林地占比從46.09%增加到53.49%，草地占比從29.21%減少到22.57%，高覆蓋度草地占比從28.80%減少到22.11%。江蘇省綠地總面積從476 840.21 hm2減小到396 255.51 hm2，減小了16.90%。

由圖5可知，2000—2005年除其他林地和高覆蓋度草地動態(tài)度明顯外，其余綠地類型動態(tài)度都較小。2005—2010年所有景觀動態(tài)度都呈明顯浮動，這期間高覆蓋度草地變化最大，動態(tài)度為-7.80%，其次是其他林地，動態(tài)度為-6.90%，灌木林動態(tài)度達-3.14%，疏林地動態(tài)度為-2.20%。總體來看，各類型林地的動態(tài)度比各類型草地變幅較小，無論在哪個時期，高覆蓋度草地動態(tài)度一直較大。這也表明草地容易受到人類活動的影響，在城鎮(zhèn)化過程中受到的傷害最明顯。

由表2可知，轉換趨勢有明顯的階段性，2000—2020年江蘇省綠地轉出面積最大的為水域， 共轉出

536.33 km2，其次為耕地，共轉出 529.48 km2；江蘇省綠地轉入面積最大來源于水域，共轉入322.26 km2，其次為耕地，共轉入 227.26 km2。2005—2010年林地轉換最明顯，林地轉為耕地和建設用地面積較大，分別為181.19、111.01 km2，其次是未利用地；草地主要轉出為水域，這期間綠地面積共減少 753.43 km2。到2015年，林地轉為耕地和建設用地的面積大幅減小，草地主要由水域轉入，這期間草地面積有小幅增加，共增加119.09 km2。以上數(shù)據(jù)均說明江蘇省綠地最大轉出類型為水域和耕地，這是由于政府在這一時期重視對水域和耕地的保護，嚴格控制水域和耕地被侵占的現(xiàn)象，這一時期也陸續(xù)出臺了《江蘇省耕地質量管理條例》《中共江蘇省委江蘇省人民政府關于進一步加強耕地保護和改進占補平衡的實施意見》等許多耕地保護政策；綠地轉入面積最大的也是水域和耕地，這主要是由于退耕還林、還草的政策要求。2005—2010年林地轉為建設用地面積較大，是因為這一時期江蘇省內各地的建設用地無序擴張導致綠地大面積減少。

2.2 綠地景觀空間格局變化

景觀指數(shù)可以反映景觀格局及其變化。由圖6可知，相關指數(shù)變化趨勢明顯，綠地斑塊密度在2005—2010年呈明顯下降趨勢，表明這段時期景觀破碎化加劇，連通性降低，這是由于該時期建設用地無序擴張，使得綠地景觀被割裂；香農均勻度指標和香農多樣性指標變化趨勢相似，除2010—2015年小幅增長外總體呈下降趨勢，2005—2010年有明顯下降，說明這一時期江蘇省綠地系統(tǒng)景觀均衡程度被破壞，多樣性下降，綠地景觀分布不均勻；綠地斑塊數(shù)量在2000—2005年和2010—2015年有所增加，但在2005—2010年大幅減少，從3 424個減少到3 199個，減少了6.57%，這都與綠地面積快速下降有關；蔓延性只在2010—2015年小幅降低，其余時段都穩(wěn)定增長，說明研究期內江蘇省綠地系統(tǒng)的優(yōu)勢景觀形成了較好的連接性；聚集度整體呈下降趨勢，說明研究期內江蘇省綠地景觀趨于分散。形態(tài)指數(shù)只在2005—2010年小幅降低，是因為人類活動對綠地景觀的影響。2010—2015年各指標都呈小幅增長趨勢，這是由于這段時期有林地面積的增加，使得景觀破碎化程度有所降低，連通性與多樣性有所增加，但總體趨勢還是趨于破碎化。

2.3 綠地生態(tài)系統(tǒng)服務價值時間變化

根據(jù)江蘇省綠地生態(tài)系統(tǒng)服務當量因子表，計算得出2000—2020年江蘇省綠地系統(tǒng)ESV及其變化（圖7）。江蘇省綠地生態(tài)系統(tǒng)服務總價值呈減少趨勢，共減少15.75%，其原因主要與江蘇省綠地面積的快速減小有關。其中，高覆蓋度草地綠地系統(tǒng)ESV在2005—2010年下降最明顯，共下降38.79%，這與該期間快速城鎮(zhèn)化導致的高覆蓋度草地大面積被占用有關。有林地、灌木林與疏林地綠地系統(tǒng)ESV呈緩慢下降趨勢；中覆蓋度草地與低覆蓋度草地綠地系統(tǒng)ESV變化較??；其他林地在2000—2005年有增長趨勢，在2005—2010年快速下降，隨后穩(wěn)定下降。

2.4 江蘇省綠地生態(tài)系統(tǒng)服務價值空間分異演變

通過克里金插值法，得到江蘇省各時段綠地生態(tài)系統(tǒng)服務價值空間分異圖（圖8）。其中，高綠地生態(tài)系統(tǒng)服務價值地區(qū)和極高綠地生態(tài)系統(tǒng)服務價值地區(qū)主要分布在江蘇省北部和西南地帶，是由于北部地區(qū)以平原為主，花果山風景區(qū)、海上云臺山國家森林公園等多個風景名勝區(qū)相連，大大小小的水庫周圍形成了自然綠化帶；西南地區(qū)以平原和低山丘陵為主，歷史文化悠久，自然條件優(yōu)越，擁有老山國家森林公園等生態(tài)展示區(qū)。中等綠地生態(tài)系統(tǒng)服務價值地區(qū)主要分布在東部沿海地區(qū)。低綠地生態(tài)系統(tǒng)服務價值地區(qū)主要分布在江蘇省的四周。極低綠地生態(tài)系統(tǒng)服務價值地區(qū)主要分布在江蘇省中部，這是因為該地區(qū)城鎮(zhèn)規(guī)?？焖贁U大，城鄉(xiāng)聚落高密度集中，該地區(qū)農業(yè)發(fā)達，是各種農產品生產加工出口的基地，因此中部以大片耕地和建設用地為主。2000—2005年江蘇省綠地生態(tài)系統(tǒng)服務價值沒有明顯變化，該時期較穩(wěn)定。2005—2010年是江蘇省綠地生態(tài)系統(tǒng)服務價值變化最顯著的一段時期，該時期中等綠地生態(tài)系統(tǒng)服務價值地區(qū)大面積減少，極高以及高等綠地生態(tài)系統(tǒng)服務價值地區(qū)部分轉變?yōu)橹械取⒌鸵约皹O低綠地生態(tài)系統(tǒng)服務價值地區(qū)。變化最大的是東部沿海地區(qū)，主要是因為人類活動導致的建設用地無序擴張使得大片綠地被占，從而破壞了生態(tài)系統(tǒng)。2010—2015年變化最明顯的還是東部沿海地區(qū)，該區(qū)域內極高以及高綠地生態(tài)系統(tǒng)服務價值地區(qū)從東北部轉向東南部，南部部分地區(qū)又從中等綠地生態(tài)系統(tǒng)服務價值轉變?yōu)榈鸵约皹O低綠地生態(tài)系統(tǒng)服務價值。2010—2020年變化較少，這是由于伴隨著城鎮(zhèn)化與工業(yè)化的快速推進，經歷過快速擴張后逐漸趨于穩(wěn)定，因此生態(tài)系統(tǒng)也逐漸趨于穩(wěn)定。2000—2020年江蘇省綠地生態(tài)系統(tǒng)服務價值的空間變化量最大的是東部沿海地區(qū)和北部平原地區(qū)。

由圖9可知，2020年江蘇省綠地ESV最高的城市是南京市，約占總綠地系統(tǒng)ESV的19.21%。其次是無錫市，約占總綠地ESV的12.51%?；窗彩泻托熘菔芯G地ESV在江蘇省總綠地系統(tǒng)ESV中占比一直較高，分別占12.01%和10.32%左右。連云港市、鎮(zhèn)江市、常州市綠地ESV略低于徐州市和淮安市，宿遷市、泰州市和揚州市綠地系統(tǒng)ESV較低。

由圖10可知，蘇州市和泰州市綠地ESV有所上升，其余城市都呈下降趨勢。其中，鹽城市綠地ESV減少最多，高達81.57億元；其次是徐州市，減少34.52億元，變化55.38%、27.96%；宿遷市和無錫市減少量較少，分別減少0.61億、1.34億元，變化了0.04%和0.01%。蘇州市綠地ESV的增加與2010年以后蘇州市高覆蓋度草地面積的增加有關，泰州市綠地ESV的增加是因為泰州市大部分土地為耕地，受退耕還林的影響使得林地面積增加。鹽城市綠地ESV的減少與鹽城市綠地的大面積減少有關，主要與2006年以后鹽城市濕地、居民區(qū)以及耕地的無序擴張有關。圖11顯示，2000—2020年變化率最高的城市是泰州；南通市綠地生態(tài)系統(tǒng)服務價值變化率在2005—2010年、2010—2015年都最高。江蘇省各地市均在2005—2010年期間變化率都呈負增長趨勢，但在2010—2015年部分城市呈增長趨勢，這可能與政府發(fā)布相關政策控制城市無序擴張有關；在2015—2020年增長趨勢又被抑制，說明政府未能在城鎮(zhèn)建設與綠地生態(tài)保護之間達到平衡，因此綠地保護是需要被重點關注的問題。

3 結論

本研究從省域角度出發(fā)，基于江蘇省近20年（2000—2020年）各類綠地數(shù)據(jù)，運用生態(tài)服務價值評價方法和空間分析法，對2000—2020年江蘇省綠地系統(tǒng)ESV的時空變化和動態(tài)變化進行分析，得出以下結論：（1）有林地和高覆蓋度草地為江蘇省綠地系統(tǒng)中的優(yōu)勢類型。城鎮(zhèn)化使江蘇省綠地面積快速減小，其中草地相比林地更容易受到人類活動的影響，綠地主要轉出為水域和耕地，這是由于近年來政府重視對耕地的保護以及快速城鎮(zhèn)化，陸續(xù)出臺《江蘇省耕地質量管理條例》《中共江蘇省委江蘇省人民政府關于進一步加強耕地保護和改進占補平衡的實施意見》等保護政策。同時，綠地主要由水域和耕地轉入，這是退耕還林等政策要求所致。（2）江蘇省綠地系統(tǒng)ESV呈穩(wěn)定下降趨勢，共減少15.75。高覆蓋度草地綠地系統(tǒng)ESV變化最大，其中江蘇省高綠地系統(tǒng)ESV主要分布江蘇省北部和西南部，中等綠地系統(tǒng)ESV主要分布在東部沿海地區(qū)，低綠地系統(tǒng)ESV主要分布在因城鎮(zhèn)規(guī)?？焖贁U大，城鄉(xiāng)聚落高密度集中，農業(yè)發(fā)達，以耕地和建設用地為主的中部地區(qū)。2005—2010年江蘇省綠地系統(tǒng)ESV變化最大，這段時期中等綠地系統(tǒng)ESV大幅減少，高綠地系統(tǒng)ESV地區(qū)向中等和低綠地系統(tǒng)ESV地區(qū)轉變，這和建設用地擴張導致的大面積綠地被侵占有關。（3）2020年，江蘇省綠地系統(tǒng)ESV最高的城市是南京市，其次為無錫市，較低的城市為泰州市和揚州市。鹽城市因濕地、居民區(qū)以及耕地的無序擴張導致綠地系統(tǒng)ESV降低最多，宿遷市和無錫市最少。江蘇省各地市在2005—2010年都呈負增長，2010—2015年因政府抑制無序擴張使得部分城市綠地系統(tǒng)ESV有所增長，但在2015—2020年這種增長趨勢又被抑制。綜上，本研究通過土地利用動態(tài)度、土地利用轉移矩陣、綠地生態(tài)系統(tǒng)服務價值評價得到相關數(shù)據(jù)，通過空間分析法揭示江蘇省2000—2020年綠地系統(tǒng)ESV空間分異，可為科學規(guī)劃綠地空間，優(yōu)化綠地景觀格局，促進綠地生態(tài)系統(tǒng)服務價值水平提升提供一定的科學依據(jù)。

參考文獻：

［1］王保忠，王彩霞，何 平，等. 城市綠地研究綜述［J］. 城市規(guī)劃匯刊，2004（2）：62-68，96.

［2］張 超，吳 群，彭建超，等. 城市綠地生態(tài)系統(tǒng)服務價值估算及功能評價——以南京市為例［J］. 生態(tài)科學，2019，38（4）：142-149.

［3］王祥榮. 生態(tài)園林與城市環(huán)境保護［J］. 中國園林，1998，14（2）：14-16.

［4］賁薛晨，余 成. 蘇州近40年土地利用與景觀格局演變特征分析［J］. 蘇州科技大學學報（自然科學版），2022，39（3）：65-73.

［5］陳 婷，楊小艷，陳龍高，等. 近30年來連云港地區(qū)土地利用/覆被的時空變化［J］. 江蘇師范大學學報（自然科學版），2015，33（4）：76-81.

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

［7］陳宏偉，許 晶，劉 娜，等. 沈陽市城區(qū)綠地生態(tài)系統(tǒng)服務價值的時空變化［J］. 應用生態(tài)學報，2018，29（10）：3391-3397.

［8］Semenyuk O V，Bodrov K S，Stoma G V，et al. Estimating the ecosystem services value of the bitzevsky forest natural park［J］. Moscow University Soil Science Bulletin，2019，74（3）：111-117.

［9］Estoque R C，Myint S W，Wang C Y，et al. Assessing environmental impacts and change in Myanmars mangrove ecosystem service value due to deforestation （2000—2014）［J］. Global Change Biology，2018，24（11）：5391-5410.

［10］Admasu S，Yeshitela K，Argaw M. Impact of land use land cover changes on ecosystem service values in the Dire and Legedadi watersheds，central Highlands of Ethiopia：implication for landscape management decision making［J］. Heliyon，2023，9（4）：e15352.

［11］Aneseyee A B，Soromessa T，Elias E. The effect of land use/land cover changes on ecosystem services valuation of Winike watershed，Omo Gibe Basin，Ethiopia［J］. Human and Ecological Risk Assessment，2020，26（10）：2608-2627.

［12］Najmuddin O，Li Z H，Khan R，et al. Valuation of land-use/land-cover-based ecosystem services in Afghanistan—an assessment of the past and future［J］. Land，2022，11（11）：1906.

［13］Biedemariam M，Birhane E，Demissie B，et al. Ecosystem service values as related to land use and land cover changes in Ethiopia：a review［J］. Land，2022，11（12）：2212.

［14］鄒麗麗，李傳武，張云峰，等. 25年來江蘇鹽城沿海灘涂生態(tài)系統(tǒng)服務價值研究［J］. 海洋環(huán)境科學，2023，42（2）：280-285.

［15］曹 波，蔡春軼，張 勇，等. 漢石橋濕地自然保護區(qū)土地利用格局及其生態(tài)系統(tǒng)服務價值動態(tài)［J］. 濕地科學，2022，20（6）：778-784.

［16］馬偉波，楊 帆，王 楠，等. 長三角城市群地區(qū)生態(tài)系統(tǒng)服務價值時空演變及驅動因素研究［J］. 生態(tài)與農村環(huán)境學報，2022，38（11）：1365-1376.

［17］潘 悅，王 錦，李婧熹，等. 滇中城市群綠色空間生態(tài)系統(tǒng)服務價值時空演變及驅動分析［J］. 水土保持研究，2023，30（2）：352-360.

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

［19］景曉棟，田貴良，班晴晴，等. 基于文獻計量的21世紀以來我國生態(tài)系統(tǒng)服務研究現(xiàn)狀及發(fā)展趨勢［J］. 生態(tài)學報，2023，43（17）：7341-7351.

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

［21］危小建，孫顯星，謝亞娟. 基于MS PA的全國主要城市綠地形態(tài)時空變化［J］. 江蘇農業(yè)科學，2020，48（24）：260-267.

［22］何 舸，劉茂松，徐 馳，等. 江蘇泰州地區(qū)耕地占地率動態(tài)變化［J］. 生態(tài)學雜志，2010，29（10）：2016-2021.

［23］劉 瑞，朱道林. 基于轉移矩陣的土地利用變化信息挖掘方法探討［J］. 資源科學，2010，32（8）：1544-1550.

［24］張秋菊，傅伯杰，陳利頂. 關于景觀格局演變研究的幾個問題［J］. 地理科學，2003，23（3）：264-270.

［25］布仁倉，胡遠滿，常 禹，等. 景觀指數(shù)之間的相關分析［J］. 生態(tài)學報，2005，25（10）：2764-2775.

［26］高 賓，李小玉，李志剛，等. 基于景觀格局的錦州灣沿海經濟開發(fā)區(qū)生態(tài)風險分析［J］. 生態(tài)學報，2011，31（12）：3441-3450.

［27］夏成琪，毋語菲. 鹽城海岸帶土地利用與景觀空間格局動態(tài)變化分析［J］. 西南林業(yè)大學學報（自然科學），2021，41（1）：140-149.

［28］謝高地，魯春霞，成升魁. 全球生態(tài)系統(tǒng)服務價值評估研究進展［J］. 資源科學，2001，23（6）：5-9.

［29］劉瑞清，李加林，何改麗，等. 浙江省綠地系統(tǒng)時空變化特征及生態(tài)損益評價［J］. 生態(tài)學雜志，2018，37（10）：3119-3129.