旅游型鄉村景觀格局變化及生態系統服務價值響應
——以江蘇省無錫市太湖國家旅游度假區為例

柳迪子， 杜守帥， 王晨旭

(1.江南大學 設計學院, 江蘇 無錫 214122; 2.北京師范大學 地理科學學部, 北京 100875)

全球60%的生態系統正在或已經退化，其中人類活動是主要誘因之一[1]。近年來，鄉村旅游迅速發展，旅游建設已經成為鄉村社會—生態系統中重要活動類型，促進產業轉型、經濟增收，為鄉村振興起到重要作用。然而，旅游業對生態資源具有較強的依賴性，盲目地建設對生態環境造成壓力，引發植被破壞、環境污染、物種減少等生態問題[2-3]，勢必也會對景觀格局、生態系統服務造成影響。目前對旅游型鄉村的研究不足，景觀格局變化過程、生態系統服務響應機制尚不明晰，造成景觀規劃、生態建設缺乏依據。因此，亟需對旅游型鄉村展開研究，探究景觀格局及生態系統服務變化過程，分析二者之間的影響機制。這對于旅游型鄉村景觀規劃、資源配置具有現實意義，進而推動鄉村可持續發展。

景觀格局變化是景觀生態學研究的核心內容，是自然、人為因素共同作用的結果，也是人類活動作用于生態系統的直接表現[4]。景觀格局指數能夠高度濃縮景觀格局信息，是反映景觀結構組成、空間配置特征的簡單量化指標[5]，是常用的分析方法。目前，相關研究已取得豐碩的成果：研究區域涵蓋了不同的地貌特征類型，如山地丘陵、平原、流域等[6-7]，總結了景觀格局變化規律和驅動因素。近年來，更關注到景觀格局與生態過程之間的關系，如：范雅雙等[8]探討水質對流域景觀格局的響應，陳羚等[9]分析了森林火險與景觀空間特征的關系，這些對相關研究具有重要意義。生態系統服務價值(ecosystem services value， ESV)是將生態系統惠益于人類社會價值進行量化，是生態系統監測管理、資源合理配置的重要依據，已成為國際生態學、經濟學、社會學研究熱點[10-11]。1997年，Costanza等[12]發表關于全球ESV評估的文章，Daily[13]梳理了生態系統服務概念、發展歷程和評估方法，成為后續研究的基礎。2001年聯合國千年生態系統評估將服務類型劃分為供給、調節、支持、文化服務，被廣泛應用[14-16]。國內研究起步較晚，早期代表性研究如：歐陽志云等[17]探討了中國生態系統的間接經濟價值，謝高地等[18]在Costanza等的研究基礎上，制定“我國陸地生態系統單位面積生態系統服務價值當量表”，對國內研究起到推動作用。近些年，相關研究取得較大進展，包括不同類型、尺度區域ESV時空變化[14]、氣候和土地利用變化的影響機制[15]、與人類福祉關系等[16]，更強調人類活動與生態系統服務的相互作用，例如，陳婧祎等[19]分析ESV與社會經濟發展的關系，楊雪婷等[16]識別了影響居民福祉的主導生態系統服務因子。此外，生態系統服務權衡、供需關系也是研究的熱點問題[20]。然而，對于生態系統服務與景觀格局關系研究仍處于起步階段，前人更多從景觀面積變化展開研究，忽視景觀結構的作用[21]。并且相關研究更偏重于大中尺度，如流域[22]、市縣域等[14]，對鄉村地區的聚焦明顯不足，難以將研究成果直接應用于鄉村景觀規劃，有必要展開村域及村域群尺度的研究。

江蘇省無錫市太湖國家旅游度假區是1992年批建的國家級旅游度假區，核心區域馬山半島由6個行政村組成，是典型的旅游型鄉村。本文以馬山半島為例，選取1984—2020年5期遙感影像進行景觀分類，計算土地利用動態度、轉移矩陣、景觀格局指數，揭示景觀格局變化過程?；诋斄恳蜃臃ㄔu估ESV，并利用敏感性指數比較各類景觀對ESV的貢獻。此外，通過ESV與景觀格局指數的相關性、回歸分析，識別影響ESV的關鍵格局特征，建立定量模型揭示二者相互作用機制，進而提出ESV增效策略，為旅游型鄉村景觀規劃、生態資源配置提供科學依據。

1 研究區概況

馬山半島位于江蘇省無錫太湖國家旅游度假區內(120°05′—120°14′N，31°36′—31°45′E)，三面太湖環繞，陸地面積約3 300 hm2。區域屬太湖流域丘陵地貌，自然資源豐富，森林覆蓋率達70%，海拔在5～266 m之間。馬山半島日照充沛、氣候溫潤，年均氣溫約17.4 ℃，日照時數1 613.6 h，降水量1 030.5 mm，屬亞熱帶季風氣候。區域自1992年批建旅游度假區至今已有近30 a的發展歷程，匯集各種旅游產業，如景區“靈山勝境”和“拈花灣”、酒店民宿、療養基地、采摘園等，地產企業也投資建設大批別墅住宅。據統計，島內現存6個行政村，19個居民點，是典型的旅游型鄉村。

2 材料與方法

2.1 數據來源及預處理

本研究以批建年份1992年為時間節點，受制于數據獲取限制，最終選取1984，1991，2001，2010，2020年5期Landsat遙感影像，均攝制于5—9月植物生長季、無云量干擾。數據源于美國地質勘探局(https:∥earthexplorer.usgs.gov/)，分辨率為30 m × 30 m。首先，利用ENVI 5.3對遙感影像預處理，包括：輻射定標、大氣校正、影像裁切等。此外，運用ENVI中最大似然法對研究區進行監督分類，分類標準參考《土地利用現狀分類》標準及謝高地等[23]對生態系統的二級分類結果，將景觀類型確定為：建設用地、闊葉林地、水田地、灌草從地、裸地、水系、其他用地。分類后，結合實地調查及高分辨率衛星圖像比對、修正，運用Ground Truth ROIs模塊進行精度驗證，得出Kappa系數大于0.92，證明分類結果可靠。另外，相關社會經濟數據源于《無錫市統計年鑒》《馬山志》等資料集。

2.2 馬山半島土地利用及景觀格局變化分析

(1) 土地利用動態度。土地利用動態度分析包括單一土地利用動態度和綜合土地利用動態度，該方法可表征區域某一時段景觀變化的劇烈程度。單一土地利用動態度描述某類景觀在某時段變化的劇烈程度，而綜合土地利用動態度反應整個區域在某時段變化的劇烈程度[24]，公式如下：

(1)

式中：K為單一土地利用動態度,即研究時段中某一景觀類型的變化速率;Ua，Ub分別為該地類在研究起止時刻的面積;T為研究時長。

(2)

式中：Lc為綜合土地利用動態度; LUi為研究起止時刻第i類景觀類型面積; ΔLUi-j為研究時段第i類景觀類型轉為非i類景觀類型面積的絕對值;T為研究時長。

(2) 景觀類型轉移矩陣。轉移矩陣能夠揭示研究期間景觀類型之間的轉移情況，該方法在土地利用變化研究中廣泛應用[25]?；贓NVI5.3中的變化檢測模塊，計算出1984—1991，1991—2001，2001—2010，2010—2020年4個時段及1984—2020年完整研究時段的轉移矩陣。根據計算結果繪制?；鶊D，圖中延伸分支的寬度對應景觀類型轉移規模大小，景觀類型轉移情況更直觀。

(3) 景觀格局指數分析。景觀格局指數是定量描述景觀格局演變及生態過程的重要方法，本文選取以往研究中廣泛應用的指數[6-7]，運用Fragstats 4.2計算。景觀格局指數包括景觀水平和類型水平，從規模、形態、連通性、破碎度等多方面表征景觀格局特征。景觀水平包括：斑塊數目(NP)、最大斑塊指數(LPI)、蔓延度(CONTAG)、面積加權平均分維數(以下簡稱“分維數”，FRAC_AM省略為FRAC)、聚集度(AI)、分離度(SPLIT)、香濃多樣性(SHDI)、香濃均勻度(SHEI)；類型水平包括：斑塊面積(CA)、最大斑塊指數(LPI)、斑塊密度(PD)、邊緣密度(ED)、分維數(FRAC)、凝聚度(COHESION)。

2.3 馬山半島生態系統服務價值響應分析

(1) 生態系統服務價值評估。運用“單位面積價值當量因子法”評估馬山半島1984—2020年ESV，參照謝高地等[23]修正后的單位面積ESV當量。查閱《無錫市統計年鑒》得到2020年無錫市谷物及其他作物農業產值(1.84×109元)、糧食作物總面積(7.95×104hm2)，根據公式(3)計算出單位面積農田自然糧食產量的經濟價值(3 300.11元/hm2)，記為標準當量。值得說明的是，自然糧食產量的經濟價值需排除人力投入產生的價值，謝高地等[18]的研究確定每hm2農田自然糧食產量價值是當年實際糧食市場價值的1/7，該系數被廣泛認可[14-15]。根據公式(4)，得到單位面積各類景觀的各單項ESV(表1)。最后，根據公式(5)計算出區域各類景觀ESV及總量。

(3)

式中：Ea為單位農田自然糧食產量的經濟價值(元/hm2);i為作物種類;mi為第i種糧食作物的面積(hm2);pi為第i種糧食作物的平均價格(元/kg);qi為第i種糧食作物單產(kg);M為糧食作物總面積;n為糧食作物種類。

Eij=eijEa

(4)

式中：Eij為j景觀類型的i生態服務功能的單位面積ESV;eij為j類生態系統的i生態服務功能相對標準當量的當量因子。

ESV=∑Aj×VCj

(5)

式中：ESV為生態系統服務價值總量;Aj為第j類景觀類型的面積;VCj為第j類景觀類型單位面積ESV。

表1 馬山半島單位面積各景觀類型各單項生態系統服務價值ESV 元/hm2

(2) 敏感性指數分析。引入經濟學中敏感性指數(coefficient of sensitive， CS)是為了核實區域ESV隨時間變化對某類景觀單位面積ESV(VC)的依賴程度，并可驗證VC準確性[26]。CS越大，ESV對該類景觀VC變化越敏感，說明其ESV貢獻度越高。若CS<1，認為ESV缺乏彈性，計算結果可信；若CS>1，則認為ESV具有彈性，計算結果準確度低。本文將研究區各類景觀VC分別上下調整50%，來評價ESV對該類景觀VC變化的敏感性，公式如下：

(6)

式中：CS為敏感性指數; ESVi為初始ESV總量; ESVj為調整后ESV總量;VCik為第k種景觀類型的初始ESV系數;VCjk為第k種景觀類型調整后ESV系數。

(3) 相關性分析。景觀格局指數與ESV的相關性分析是探究二者關系的重要步驟[21]，幫助識別對ESV影響顯著的景觀格局特征。以1984—2020年馬山半島的景觀格局指數及ESV為數據源，運用SPSS 23.0進行Pearson相關性檢驗，當顯著性p<0.05，說明該組變量存在相關性。

(4) 回歸分析。為進一步探究景觀格局變化對ESV的作用機制，需在ESV與景觀格局指數之間建立定量模型。將相關性表現顯著的景觀格局指數逐一與ESV進行線性回歸，建立回歸模型，并根據回歸模型中自變量的標準回歸系數(調整后R2)，進行解釋度排序，以此評價景觀格局指數對ESV作用力大小[27]。標準回歸系數消除了不同變量之間所取單位不同的影響，其絕對值大小反映X對Y的影響程度[28]。

3 結果與分析

3.1 馬山半島景觀類型及演變特征

(1) 景觀面積變化特征。馬山半島土地利用景觀類型分為：建設用地、闊葉林地、水田地、灌草叢、裸地、水系和其他用地。由圖1和表2可知，闊葉林地是區域主要景觀類型，占35%～51%；其次為水田地和水域，分別占12.7%～17.3%和11.7%～16.8%；裸地面積最小，僅占0.3%～6.1%。從各類景觀面積變化情況來看：①總體來說，建設用地、灌草叢增加，闊葉林地、水系、裸地和其他用地縮減，水田地基本持平；②建設用地增幅最大，且持續增加，不斷向周邊景觀類型侵占，其中1991—2010年增速最快；灌草叢增幅次之，其中1991—2001年增幅最大；③闊葉林地縮減最嚴重，1984—2010年減至最低值，2010—2020年少量回升；④水系次之，1984—2001年減至最低值，2001—2020年少量回升；裸地的縮減幅度最大。

(2) 土地利用動態度。據表3可知，馬山半島綜合土地利用動態度呈上升趨勢，2001年以后明顯高于之前，2010—2020年動態最劇烈，1991—2001年最穩定。1984—1991年，裸地動態最劇烈，其次為水田地、灌草叢、水系。1991—2001年，大部分景觀類型動態度增加，灌草叢、其他用地最劇烈，建設用地、水田地次之。2001—2010年，多數景觀類型動態度下降，建設用地、其他用地最活躍。2010—2020年，裸地、水田地、闊葉林地動態劇烈，而建設用地趨于穩定。

表2 1984-2020年馬山半島不同景觀類型面積 hm2

表3 1984-2020年馬山半島土地利用動態度 %

(3)景觀轉移過程。根據?；鶊D可知(圖2a—2d)：1984—1991年，水系明顯向水田地、其他用地轉移，部分景觀類型向裸地轉移，而裸地主要向建設用地轉移。1991—2001年，闊葉林地大規模向灌草叢、水田地、裸地轉移，裸地向建設用地、水田地轉移；2001—2010年，闊葉林地持續向灌草叢、水田地轉移，建設用地主要由水田地、裸地轉移；2010—2020年，闊葉林地大規模由水田地、灌草叢轉移，裸地向建設用地、水田地轉移。從整體轉移情況來看(圖2e)：各類景觀之間轉換頻繁，闊葉林地主要向灌草叢、水田地、建設用地轉移；建設用地主要向闊葉林地、水田地、灌草叢、裸地侵占。

3.2 馬山半島景觀格局變化特征

(1) 類型水平景觀格局指數。根據景觀格局指數增降情況(圖3)可知，研究期間，建設用地形態破碎化、復雜化，并趨于聚集；闊葉林地優勢性減弱，形態破碎、復雜化；水田地和灌草叢破碎化，并且邊緣復雜化；裸地形態結構簡化，趨于聚集；其他用地形態破碎、復雜化，趨于聚集；水系幾乎保持穩定。

(2) 景觀水平景觀格局指數。由圖4可知，研究期間，區域景觀破碎化、結構復雜化，連通性和優勢度均減弱。具體來看，1984—1991年，斑塊趨于聚集、連通性加強，景觀豐度下降、斑塊形態趨于復雜；1991—2001年，各景觀指數增降幅度最大，斑塊趨于離散，形態變得復雜、破碎，連通性和優勢性均減弱；2001—2010年各項指數持續之前的態勢，但變化速率明顯降低；2010—2020年，各項指數與1991—2010年的趨勢相反，景觀的連通性、優勢度增強，結構簡化，斑塊趨于聚集。

3.3 馬山半島生態系統服務價值變化及敏感性分析

(1) 生態系統服務價值變化。根據1984—2020年馬山半島不同時期各類景觀的各單項ESV及其變化情況(表4—5)來看，1984—2020年，ESV總量呈下降趨勢，累計減少7.08×105元。具體來看：1984—1991年，ESV縮減5.01×105元，其中水系ESV縮減最嚴重；1991—2001年，ESV持續縮減5.03×105元，闊葉林地ESV縮減最嚴重，灌草叢ESV明顯增加；2001—2010年，ESV顯著回升，主要因為水系ESV明顯增加；2010—2020年，ESV少量回升，闊葉林地ESV增加顯著，但水系ESV縮減抵消了部分闊葉林地增量。

從各類景觀ESV比重來看(表4)，水系與闊葉林地ESV遠高于其他景觀類型，具體排序為：水系>闊葉林地>灌草叢>水田地>裸地。從各類景觀ESV的增降情況來看，1984—2020年，水系ESV縮減最嚴重，縮減6.63×105元；闊葉林地次之，縮減1.26×105元；灌草叢增加8.12×104元；水田地、裸地基本持平。從各單項ESV來看(表5)：各單項ESV均減少。其中調節服務損失最嚴重，縮 減6.17×105元，占總損失的87.2%。由此可見，旅游業發展對水系、闊葉林地ESV影響較大，且調節服務衰退最嚴重。

圖1 1984-2020年馬山半島土地利用景觀分類

表4 1984-2020年馬山半島各景觀類型生態服務價值(ESV)及總量變化 104元

(2) 敏感性指數分析。敏感性指數分析可比較各類景觀對ESV貢獻度差異。據表6可知，各類景觀CS均小于1，說明ESV對各類景觀VC缺乏彈性，計算結果可信。各類景觀CS差異顯著，排序為：水系>闊葉林地>灌草叢>水田地>裸地。水系CS最高，說明其ESV貢獻最大，例如，2020年水系CS為0.561，表明水系VC每變化1%，引起ESV變化0.561%，原因是水系VC遠高于其他景觀，并且面積較大；闊葉林地具有絕對的面積優勢，VC僅次于水系，因此它對ESV貢獻也較高；水田地、裸地的貢獻較低。從時間序列來看：水系由于面積縮減，ESV貢獻度在減弱；闊葉林面積雖然也在縮減，但其面積優勢顯著，其VC遠不及水系，因此在水系的貢獻度下降后，闊葉林貢獻度便凸顯；灌草叢的貢獻也在上升，水田地較穩定。

表5 1984-2020年馬山半島各單項生態服務價值(ESV)變化 104元

3.4 馬山半島ESV與景觀格局指數的相關性分析

(1) 景觀格局指數與ESV總量的相關性。ESV與景觀格局指數的相關性分析(表7)可知，1984—2020年馬山半島ESV與水系面積呈顯著正相關，與景觀分維數、建設用地分維數、凝聚度呈負相關。說明水系面積越大，區域ESV越高；景觀形態越復雜、建設用地越聚集，ESV越低。

(2) 景觀格局指數與各單項ESV的相關性。供給、調節、支持、文化服務功能價值分別與景觀格局指數進行相關性分析(表7)。

①供給服務價值與蔓延度、水系面積呈正相關，與多樣性、分維數及建設用地分維數、凝聚度呈負相關。說明景觀連通性越強、水系規模越大，供給服務越強；景觀結構越復雜、建設用地越聚集、復雜，供給服務越弱。②調節服務價值與水系面積呈顯著正相關，與分維數、建設用地凝聚度呈負相關。說明水系規模越大，調節功能越強；景觀整體形態越復雜、建設用地越聚集，調節功能越弱。③支持服務價值與蔓延度呈顯著正相關，與多樣性、分離度呈負相關。說明景觀連通性越強，支持服務越強；景觀結構越復雜、分離，支持服務越弱。另外，在類型層面，闊葉林地面積、灌草叢分維數與支持服務價值的相關性顯著，表現為：闊葉林面積越大，支持服務越強；灌草叢形態越復雜，支持服務越弱。另外，支持服務還與建設用地分維數、灌草叢面積呈負相關，建設用地形態越復雜、灌草叢面積越大，支持服務越弱。④文化服務價值與多樣性呈顯著負相關，說明景觀結構越復雜，文化服務越弱。另外，文化服務價值還與分離度、建設用地分維數、灌草叢分維數呈負相關，與蔓延度、聚集度及闊葉林面積呈正相關。這說明景觀連通性、聚集程度越高，文化服務越強；建設用地、灌草叢斑塊形態越復雜，文化服務價值越低。

3.5 馬山半島ESV與景觀格局指數的回歸分析

將與ESV相關性顯著的景觀格局指數再分別與各類ESV進行回歸分析，進一步揭示景觀格局變化對ESV的作用機制，并建立回歸模型(表8)。根據標準化回歸系數排序，列出對ESV及各單項ESV解釋度最高的前3項指數，以此識別對不同類型生態系統服務效能影響最大的景觀格局特征。

①ESV總量、供給服務及調節服務受水系面積、建設用地聚集度、景觀形態復雜性的影響最大，究其原因：水系生態價值極高，尤其以水文調節和水資源供給功能，幾乎主導著馬山半島的供給、調節以及生態系統綜合服務能力。隨著旅游開發強度增加，建設用地規模、密度加大，人類對景觀的擾動范圍也在擴大，侵占生態用地導致ESV降低；另外，馬山半島以丘陵地貌為主，人類對景觀的擾動主要發生在地勢平緩的山坳處，建設用地、灌草叢不斷向山坳深處蔓延，闊葉林地持續收縮(圖5)，造成景觀形態復雜化，同時高生態價值的闊葉林地也嚴重縮減。②支持服務受景觀多樣性、灌草叢形狀復雜性、闊葉林面積的影響最大，究其原因，景觀異質性很大程度受景觀基質(闊葉林地)影響，闊葉林為馬山半島提供強大的土壤保持、維持養分循環、生物多樣性等支持功能，生態價值極高，其規模決定著馬山半島整體支持服務效能；而隨著灌草叢向山林蔓延，闊葉林地被侵占，造成景觀整體的連通性下降、ESV縮減。③文化服務受景觀多樣性、建設用地和灌草叢形狀復雜性的影響最大，究其原因：景觀異質化伴隨著景觀結構復雜化、破碎化，尤其闊葉林地表現最明顯，非常不利于馬山半島整體美感；建設用地和灌草叢向山林蔓延，造成形狀復雜化，意味著闊葉林收縮嚴重(圖5)，其美學價值也嚴重降低。

圖3 1984-2020年馬山半島景觀水平景觀指數變化情況

表6 1984-2020年馬山半島各景觀類型生態服務系統價值ESV敏感性分析

表7 馬山半島生態系統服務價值與景觀格局指數的相關性

圖4 1984-2020年馬山半島景觀水平景觀指數變化情況

表8 馬山半島生態系統服務價值與景觀格局指數的回歸分析

圖5 人類干擾丘陵景觀導致景觀形態復雜化過程示意圖

4 討論與結論

4.1 討 論

(1) 馬山半島不同時期土地利用景觀格局、ESV變化的影響因素。通過梳理馬山半島地方志及相關歷史資料，探討1984—2020年4個時段景觀格局、ESV變化的影響因素：①1984—1991年，馬山半島還未發展旅游業，景觀動態相對穩定。此時最明顯的變化是，部分水系(主要為內陸坑塘)向水田地轉變，造成水系ESV下降嚴重。查閱資料了解到，20世紀80年代末，馬山半島淡水養殖效益下滑，挫傷農民“內塘養魚”的積極性，故將坑塘填蓋，投入到收益更穩定種植業使用。由此可見，旅游度假區建設之前，農業生產效益是造成景觀格局及ESV變化的主要因素；②1991—2001年是旅游度假區建設初期，據1991，2001年景觀類型圖可知，旅游建設的影響僅體現在以靈山景區為中心的小范圍區域，表現為建設用地向周邊景觀蔓延。該時期最明顯的變化是闊葉林大面積向灌草叢轉變，造成闊葉林地ESV嚴重下降。據《馬山志》記載，20世紀80年代末，馬山半島林果業發展迅速，原村集體茶廠、林果場均分至個人承包經營，山林分到農戶，激發了農民種植茶葉、楊梅等作物的積極性，這也引發灌草叢地激增、闊葉林地縮減，致使ESV下降。由此可見，此時景觀格局及ESV受旅游業、種植業發展的共同影響。③2001—2010年，靈山景區初具規模，馬山半島吸引了大批投資建設，島內新增多處度假酒店、康養基地等休閑娛樂產業，建設用地在島內散布，區域景觀破碎化。同時，靈山景區附近規劃新農村居民點，人類對闊葉林地擾動范圍進一步擴大，致使其ESV跌至最低。2007年，政府對太湖流域實施“退耕還湖”政策，沿太湖邊緣的圩田轉變為生態價值更高的水域，這對馬山半島ESV起到積極作用。④2010—2020年，隨著旅游度假區建設不斷完善，馬山半島步入穩定建設期。建設用地擴張速度明顯減緩。2010年政府頒布的林地保護政策，引導全社會嚴格保護林地、優化林地資源配置。這推動了馬山半島退耕還林，林地ESV顯著回升，也扭轉了馬山半島景觀持續破碎化的趨勢，連通性顯著增強。由此可知，隨著旅游度假區建設程度不同，景觀格局、ESV的主要影響因素也不同，農業活動、政策頒布與旅游建設共同作用于旅游型鄉村景觀格局與ESV。

(2) 影響馬山半島ESV的關鍵景觀格局特征及增效策略。根據相關性分析、回歸分析，可總結影響馬山半島ESV的主要景觀格局特征，并提出相應的景觀格局優化策略(圖6)，實現ESV增效。其中，標準回歸系數排序靠前的指數所對應的景觀格局特征可作為主要措施，其他作為輔助措施：①生態系統服務綜合效能方面，馬山半島要控制建設用地肆意蔓延，修復山坳深處已破壞植被；同時增加水體景觀，發揮水生態系統服務能力。②供給服務方面，可通過建設水塘、疏通水系來滿足水資源補給、食物生產灌溉的需求；限制建設用地擴張，保障水田、闊葉林等景觀的供給服務。③調節服務方面，通過修復闊葉林、增加水域面積，加強空氣凈化、水文調節能力；另外，在高密度建設區域增加綠地和水景，調節微氣候。④支持服務方面，首先要阻止苗圃、種植園等形式的灌草地向山林蔓延；構建生態安全格局、生態廊道，保障闊葉林斑塊的優勢性與連通性，為生物提供優質棲息地，維持物質能量流通。⑤文化服務方面，通過培育修復山林景觀，恢復闊葉林地美學價值；整合補綴破碎斑塊，提升馬山半島的整體美感。

(3) 本文以無錫太湖國家旅游度假區馬山半島為例，分析1984—2020年景觀格局變化及ESV響應，揭示二者之間相互作用關系。研究方法來看： ①利用遙感影像進行景觀分類，結合土地利用動態度、轉移矩陣、景觀格局指數等方法，可以從規模、動態、過程、結構多方面呈現鄉村景觀格局的演變過程； ②當量因子法適用于村域群尺度ESV評估，敏感度指數分析可比較景觀類型對ESV貢獻的差異； ③基于SPSS軟件對ESV與景觀格局指數相關性、回歸分析，成功識別影響生態系統服務的主要景觀格局特征，定量表達景觀格局與ESV之間的作用機制，為景觀資源配置提供依據。

(4) 本文仍存在不足： ①受到遙感影像分辨率限制，分析結果僅能反應大體變化趨勢。獲取更高精度數據有助于捕捉細節變化，使研究結果更準確； ②本文僅研究馬山半島單一案例，缺少同類區域比較，具有一定的局限性。此外，我國地形地貌復雜多樣，本案例對丘陵地貌鄉村更具參考價值，普適性不足。后續可參考本文思路，對更多不同地貌區的旅游型鄉村進行比較研究，形成更充分的結論； ③由于鄉村社會—生態系統的復雜性，本文僅對馬山半島ESV增效提出策略性措施，具體優化方案還需深入剖析鄉村生態系統服務真實需求，解決各功能之間權衡配置，確定區域生態系統服務功能定位，這也是進一步展開鄉村生態系統服務研究的關鍵。

圖6 影響馬山半島生態系統服務的關鍵景觀格局特征及相應增效措施

4.2 結 論

(1) 隨著無錫市太湖國家旅游度假區建設發展，馬山半島土地利用景觀格局、ESV變化顯著。各類景觀轉換頻繁，闊葉林地縮減嚴重，建設用地持續擴張。景觀格局復雜化、破碎化，連通性減弱。同時，區域ESV明顯降低，各單項ESV均下降，調節服務功能衰退最嚴重。

(2) 旅游建設、農業需求、政策頒布是影響旅游型鄉村景觀格局、ESV變化的主要因素。并且在旅游業發展的不同階段，各因素的影響程度也不同。

(3) 旅游型鄉村景觀格局與ESV之間關系顯著，不同的景觀類型、格局特征對區域ESV的作用程度不同。通過相關性、回歸分析，可識別影響區域ESV的主要景觀格局特征，從而提出相應的ESV增效措施和策略，對科學配置、合理規劃景觀資源具有重要意義。