蔡夢(mèng)卿，黃 璐，嚴(yán)力蛟

（1. 浙江大學(xué) 生命科學(xué)學(xué)院，浙江 杭州 310058；2. 浙江大學(xué) 新農(nóng)村發(fā)展研究院，浙江 杭州 310058）

聯(lián)合國(guó)千年生態(tài)系統(tǒng)評(píng)估表明：超過(guò)60%的全球生態(tài)系統(tǒng)服務(wù)呈退化或過(guò)度利用趨勢(shì)，自然因素和由人類(lèi)活動(dòng)引起的土地利用變化是導(dǎo)致生態(tài)服務(wù)功能退化的重要驅(qū)動(dòng)因素[1]。水資源供給是重要的生態(tài)系統(tǒng)服務(wù)之一，水的可用性在空間和時(shí)間上的變化可以顯著影響任何生態(tài)系統(tǒng)的功能和生產(chǎn)力[2]，對(duì)生態(tài)系統(tǒng)和區(qū)域經(jīng)濟(jì)的可持續(xù)發(fā)展具有重要作用[3?4]。有研究[5?6]表明：氣候變化和土地利用變化是引起產(chǎn)水量變化的驅(qū)動(dòng)因素。氣候變化可以通過(guò)改變地區(qū)的降水量、氣溫和蒸散量來(lái)影響產(chǎn)水量[7]。PESSACG等[8]和LANG等[3]發(fā)現(xiàn)降水量變化會(huì)導(dǎo)致產(chǎn)水量的顯著差異。土地利用變化可能通過(guò)改變?nèi)霛B速率和蒸散速率影響水量[9?10]。建設(shè)用地的增加會(huì)提高產(chǎn)水量[11?12]。但目前氣候因素和土地利用變化對(duì)產(chǎn)水量影響的貢獻(xiàn)或相關(guān)程度仍存在爭(zhēng)議。顧晉飴等[13]在太湖流域的研究發(fā)現(xiàn)：降水量較土地利用/覆被更顯著地驅(qū)動(dòng)產(chǎn)水量的變化。竇攀烽等[14]利用情景分析發(fā)現(xiàn)氣候變化對(duì)寧波地區(qū)產(chǎn)水量影響更為顯著。HU等[15]通過(guò)主成分分析和相關(guān)性分析發(fā)現(xiàn)景觀(guān)格局對(duì)洞庭湖地區(qū)產(chǎn)水量的驅(qū)動(dòng)作用大于氣候因素。如何有效綜合評(píng)估氣候因素和土地利用變化對(duì)產(chǎn)水量的影響仍未形成相對(duì)統(tǒng)一的結(jié)果。InVEST 模型產(chǎn)水量模塊基于水量平衡原理，以柵格為單元定量估算不同景觀(guān)類(lèi)型的水源供給能力，在評(píng)估與水相關(guān)的生態(tài)系統(tǒng)服務(wù)時(shí)較為可靠[16?18]。隨著城鎮(zhèn)化的進(jìn)程加塊，城市結(jié)構(gòu)、內(nèi)部特征差異和城市生態(tài)系統(tǒng)服務(wù)供給等出現(xiàn)差異[19?20]，并隨城市化梯度差異表現(xiàn)[21]。利用梯度分析研究空間布局及其時(shí)空變化過(guò)程可從空間上識(shí)別和評(píng)估土地利用變化的影響[22?23]。杭州市是浙江省政治、經(jīng)濟(jì)、文化中心，近年來(lái)城市綠地侵占與恢復(fù)并存，產(chǎn)水量的時(shí)空變化特征具有復(fù)雜性和不確定性，深刻影響著城市可持續(xù)發(fā)展和居民福祉。本研究利用遙感數(shù)據(jù)和觀(guān)測(cè)數(shù)據(jù)，在InVEST模型支持下，評(píng)估杭州市2000—2015年產(chǎn)水量，分析氣候變化和土地利用變化與年產(chǎn)水量的相關(guān)性，以期為提高城市水資源保護(hù)和利用、促進(jìn)城市生態(tài)系統(tǒng)服務(wù)發(fā)展提供參考。

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

1.1 研究區(qū)概況

浙江省杭州市位于中國(guó)東南沿海北部、長(zhǎng)江三角洲南沿 (29°11′～30°33′ N，118°21′～120°30′ E)，屬于亞熱帶季風(fēng)氣候，2019年全市平均氣溫約18.1 ℃，年降水量為1 647.7 mm。夏季高溫多雨，冬季低溫干燥。杭州市地形復(fù)雜多樣，東部為浙北平原地區(qū)，西部為浙西丘陵地帶，京杭大運(yùn)河和錢(qián)塘江穿城而過(guò)。截至2019年，全市下轄10區(qū)2縣1市，分別為上城區(qū)、下城區(qū)、拱墅區(qū)、西湖區(qū)、江干區(qū)、濱江區(qū)、蕭山區(qū)、余杭區(qū)、臨安區(qū)、富陽(yáng)區(qū)、桐廬縣、淳安縣和建德市。

1.2 數(shù)據(jù)來(lái)源及處理

本研究所采用的數(shù)據(jù)包括Landsat遙感影像數(shù)據(jù)、氣象觀(guān)測(cè)數(shù)據(jù)、土壤數(shù)據(jù)和地形高程數(shù)據(jù)。采用中國(guó)科學(xué)院地理科學(xué)與資源研究所提供的長(zhǎng)三角城市群2000、2005、2010、2015年共4期土地利用數(shù)據(jù)為基礎(chǔ)數(shù)據(jù)，將地類(lèi)劃分合并為耕地、林地、草地、水體、建設(shè)用地和未利用土地等6類(lèi)。氣象數(shù)據(jù)由國(guó)家氣象科學(xué)數(shù)據(jù)中心提供的浙江省基準(zhǔn)地面氣象觀(guān)測(cè)站及自動(dòng)站共計(jì)18個(gè)站點(diǎn)的2000—2015年觀(guān)測(cè)數(shù)據(jù)集，去除異常值后在ArcGIS平臺(tái)下進(jìn)行投影，使用空間分析工具中的克里金(Kringing)插值分析得到2000、2005、2010和2015年的降水量、潛在蒸散量分布圖，其中潛在蒸散量按Penman-Monteith (PM)公式計(jì)算得到，實(shí)際蒸散量由InVEST模型產(chǎn)水模塊運(yùn)行計(jì)算得到。利用中科院寒區(qū)旱區(qū)數(shù)據(jù)中心提供的中國(guó)土壤類(lèi)型圖提取杭州市土壤根系深度、物理組成、有機(jī)質(zhì)等信息，在ArcGIS平臺(tái)柵格數(shù)據(jù)空間疊加運(yùn)算下得到土壤有效含水量(AWC)[24]。采用中國(guó)科學(xué)院地理所提供的數(shù)字高程模型(DEM)計(jì)算得到流域矢量數(shù)據(jù)。所有數(shù)據(jù)進(jìn)行柵格化處理，轉(zhuǎn)為同一投影的柵格圖像，空間坐標(biāo)系為Krasovsky_1940_Albers，像素單元大小為 30 m×30 m。

1.3 研究方法

1.3.1 產(chǎn)水量及其計(jì)算 利用InVEST模型的產(chǎn)水模塊，考慮實(shí)際蒸散量，建立年降水量與Budyko熱液耦合平衡的關(guān)系，估算杭州市不同土地利用類(lèi)型的產(chǎn)水功能。其核心算法是運(yùn)用水量平衡法結(jié)合氣候、地形和土地利用類(lèi)型計(jì)算得出每個(gè)柵格的產(chǎn)水量。產(chǎn)水量為區(qū)域內(nèi)各柵格的降水量減去實(shí)際蒸散量。產(chǎn)水量計(jì)算公式如下：

式(1)中：Yx為柵格x的年產(chǎn)水量(mm)；TAEx為柵格x的年實(shí)際蒸散量(mm)；Px為柵格x的降水量(mm)。TAEx根據(jù)Budyko曲線(xiàn)近似得到，公式[25]為：

式(2)～(5)中：Rx為柵格x的Budyko干燥度指數(shù)，是潛在蒸散量(TPEx)與降水量(Px)的比值。CKx為植被蒸散系數(shù)，在不同的植被類(lèi)型下該值是不同的，表示不同發(fā)育期植物潛在蒸散量與參考蒸散量(TE0x)的比值，具體取值參考模型默認(rèn)范例。TAEx通過(guò)TE0x直接計(jì)算，并由TE0x和Px共同決定其取值。wx是土壤可利用水量與降水量的比值。根據(jù)DONOHUE等[26]的定義：

式(6)中：CAWx為柵格x的土壤有效含水量(mm)，Z為季節(jié)常數(shù)，基于杭州市自然地理特征及相關(guān)參考文獻(xiàn)將其取值為6.5[27]。本研究年產(chǎn)水總量為杭州市所有柵格單元產(chǎn)水量的總和。

1.3.2 梯度分析 以城市原點(diǎn) (紫薇園，30°15′00″ N，120°10′12″ E)為城市中心、以 5 km 為間距，建立覆蓋杭州市主城區(qū)的10個(gè)同心圓梯度帶[28]。為消除因面積不同產(chǎn)生的結(jié)果差異，采用ArcGIS空間分析工具Zonal對(duì)每個(gè)梯度帶要素進(jìn)行空間分布均值統(tǒng)計(jì)分析。

1.3.3 影響因素相關(guān)性分析 選取年降水量 (Pre)、年實(shí)際蒸散量 (AET)和年均氣溫(Tem)作為影響年產(chǎn)水量的氣象因子指標(biāo)，根據(jù)主成分分析(PCA)從密度、面積、邊緣、形狀、多樣性和聚散性6個(gè)方面選取斑塊密度(PD)、最大斑塊指數(shù)(LPI)、邊緣密度(ED)、景觀(guān)形狀指數(shù)(LSI)、香農(nóng)多樣性指數(shù)(SHDI)和蔓延度(CONTAG)等6個(gè)景觀(guān)格局指數(shù)作為影響年產(chǎn)水量的景觀(guān)格局指標(biāo)。對(duì)研究區(qū)進(jìn)行1 km×1 km網(wǎng)格化處理并去除邊緣值，分別構(gòu)建氣象因子和景觀(guān)格局指數(shù)網(wǎng)格數(shù)據(jù)集，其中景觀(guān)格局指數(shù)由土地利用類(lèi)型數(shù)據(jù)通過(guò)Fragstats 4.2軟件計(jì)算得到。利用SPSS 25中Pearson分析計(jì)算得到年產(chǎn)水量與氣象因子、景觀(guān)格局指數(shù)的相關(guān)性。

2 結(jié)果與分析

2.1 年產(chǎn)水量的年間變化

由表1可知：2000—2015年間杭州市年產(chǎn)水總量與平均值均呈先下降后上升趨勢(shì)。就總量來(lái)看，2005年最小，為 9.81×109m3，比 2000年下降了 23.24%，2015年達(dá)到最大值 (24.48×109m3)；15 a平均增長(zhǎng)率為6.10%。就柵格單元平均值來(lái)看，最小值為 582.42 mm (2005年)，最大值為1 457.24 mm(2015年)，4 a相對(duì)誤差為5.90%，2000—2015年平均增長(zhǎng)率為6.17%。

從不同土地利用類(lèi)型來(lái)看，2000—2015年間各類(lèi)型土地柵格單元年產(chǎn)水量均呈先下降后上升趨勢(shì)(表1)。建設(shè)用地高于其他土地利用類(lèi)型，最大值為1 796.25 mm (2015年)，水體產(chǎn)水量低于其他土地利用類(lèi)型，最小值為214.30 mm (2005年)。2000—2015年間，各類(lèi)型土地年產(chǎn)水量增長(zhǎng)率從大到小依次為水體(181.72%)、未利用土地(95.98%)、林地(90.10%)、耕地(86.97%)、建設(shè)用地(79.42%)、草地(77.78%)。

表1 2000—2015年杭州市的產(chǎn)水量Table 1 Water yield of Hangzhou during 2000?2015

2.2 年產(chǎn)水量的空間變化

從空間區(qū)域上來(lái)看，2000—2015杭州市年產(chǎn)水量呈現(xiàn)東北高、西南低的分布特征 (圖1)。高值區(qū)主要分布在東北部的主城區(qū)，且呈聚集擴(kuò)張趨勢(shì)，柵格單元多年平均產(chǎn)水量為1 042.68～1 254.85 mm；其他區(qū)縣市(富陽(yáng)區(qū)、臨安區(qū)、桐廬縣、淳安縣、建德市)相對(duì)較小，為933.57～1 009.06 mm (圖2)。

圖1 杭州市 2000—2015年產(chǎn)水量空間分布示意圖Figure 1 Spatial distribution of water yield in Hangzhou in 2000－2015

圖2 杭州市各縣區(qū)年產(chǎn)水量變化趨勢(shì)Figure 2 Changes of water yield of counties and districts in Hangzhou

從梯度變化來(lái)看，年產(chǎn)水量呈現(xiàn)隨著與城市中心距離的增加先增大后減小趨勢(shì)(表2)，最高值均出現(xiàn)在距離城市中心10 km處。對(duì)比距離城市中心5～15 km和40～50 km的年產(chǎn)水量發(fā)現(xiàn)：2000、2005、2010和2015年的差異率分別為4.50%、16.50%、9.25%和10.31%，這表明年產(chǎn)水量的城鄉(xiāng)變化在2005年前后出現(xiàn)拐點(diǎn)。

表2 2000—2015年杭州市沿梯度變化的年產(chǎn)水量Table 2 Water yield along gradient in Hangzhou during 2000?2015

2.3 土地利用和氣象因子對(duì)年產(chǎn)水量的影響

2.3.1 土地利用類(lèi)型變化 杭州市各土地利用類(lèi)型轉(zhuǎn)移矩陣表明 (表 3)：2000—2015年間土地利用變化最劇烈的為耕地與建設(shè)用地之間的轉(zhuǎn)化，耕地面積減少了2.75%，建設(shè)用地面積增加了3.67%，15 a間有3.24%耕地轉(zhuǎn)化為建設(shè)用地。結(jié)合表1可知：杭州市2000—2015年間建設(shè)用地面積比例增加顯著，導(dǎo)致杭州市2015年產(chǎn)水量顯著大于2000年產(chǎn)水量。

表3 杭州市 2000—2015年各土地利用類(lèi)型轉(zhuǎn)移矩陣Table 3 Statistics of conversion rates of landscape types in Hangzhou from 2000 to 2015

2.3.2 景觀(guān)格局變化 從景觀(guān)格局來(lái)看，2000—2015年間杭州市存在景觀(guān)明顯破碎化和異質(zhì)性(表4)。15 a間，PD先下降后上升，LPI下降1.23、ED上升1.07，表明景觀(guān)破碎化程度明顯；LSI上升3.69，表明斑塊形狀較不規(guī)則、復(fù)雜度不斷增加；CONTAG下降1.77，表明景觀(guān)連通性處于下降趨勢(shì)，反映景觀(guān)破碎化增加；SHDI由0.95上升至1.00，表明景觀(guān)異質(zhì)性增加。從梯度變化來(lái)看，隨著與城市中心距離的增加，LPI、CONTAG總體呈現(xiàn)先減小后逐步增加的趨勢(shì)，最大斑塊面積和景觀(guān)連通性先減小后增大；ED、LSI、PD、SHDI總體呈現(xiàn)先增加后逐步減小的趨勢(shì)，表明隨著與城市中心距離的增加，景觀(guān)破碎化和異質(zhì)性先增大后減小(圖3)。其中，CONTAG和SHDI的拐點(diǎn)出現(xiàn)在10 km處，表明在距離城市中心10 km處，景觀(guān)連通性最小、異質(zhì)性最大；PD、ED、LSI的拐點(diǎn)出現(xiàn)在20 km處，表明在距離城市中心20 km處，景觀(guān)破碎化最為嚴(yán)重；LPI的拐點(diǎn)則在15和30 km處，表明在距離城市中心15和30 km處，最大斑塊面積減小導(dǎo)致一定程度的景觀(guān)破碎化。因此，在距離城市中心10～20 km處存在最為顯著的景觀(guān)破碎化和異質(zhì)性。結(jié)合年產(chǎn)水量發(fā)現(xiàn)：?jiǎn)挝幻娣e產(chǎn)水量拐點(diǎn)出現(xiàn)在距離城市中心10 km處，與景觀(guān)格局的梯度變化趨勢(shì)較為一致。

圖3 杭州市各景觀(guān)格局指數(shù)沿梯度變化趨勢(shì)Figure 3 Changes of landscape pattern indices along gradient in Hangzhou

表4 2000—2015年杭州市各景觀(guān)格局指數(shù)Table 4 Statistics of landscape pattern indices in Hangzhou from 2000 to 2015

2.3.3 氣象因子變化 從時(shí)間上看 (表 5)：2000—2015年間，杭州市年降水量最小的是2005年，最大的是2015年。與2000年相比，2005年降水量減少183.54 mm (13.90%)，2015年 增 加696.21 mm(52.75%)。年均氣溫最低的是2000年，最高的是2015年，相差0.31 ℃。年均實(shí)際蒸散量最小值和最大值分別出現(xiàn)在2005和2010年；與2000年相比，2005年實(shí)際蒸散量減少9.21 mm (1.64%)，2010年增加2.31 mm (0.41%)。結(jié)合表1可知：降水量變化與年產(chǎn)水量趨勢(shì)一致。從空間上看，各縣區(qū)氣象因子存在差異(表6)。主城區(qū)與富陽(yáng)區(qū)年降水量為1 540～1 545 mm，其他 4區(qū)縣市均高于 1545 mm。主城區(qū)年均氣溫為 17.35～17.39 ℃，其他 5區(qū)縣市為17.34～17.38 ℃。主城區(qū)年實(shí)際蒸散量明顯低于其他5區(qū)縣市，其中最低值出現(xiàn)在下城區(qū)(289.10 mm)，最高值出現(xiàn)在淳安縣(617.13 mm)。從城鄉(xiāng)梯度來(lái)看，在距城市中心50 km內(nèi)，年降水量變化不大；隨著距離增加，年均氣溫逐漸下降，年實(shí)際蒸散量先減小后增大，最小值出現(xiàn)在距離10 km處(圖4)。結(jié)合2000—2015年產(chǎn)水量的空間分布特征和梯度變化趨勢(shì)可知：年產(chǎn)水量空間變化與年實(shí)際蒸散量表現(xiàn)相似。

圖4 杭州市氣象因子沿梯度變化趨勢(shì)Figure 4 Changes of climate factors along gradient in Hangzhou

表5 2000—2015年杭州市的氣象因子Table 5 Statistics of climate factors in Hangzhou from 2000 to 2015

表6 2000—2015年杭州市各縣區(qū)的氣象因子Table 6 Statistics of climate factors of counties and districts in Hangzhou from 2000 to 2015

2.3.4 年產(chǎn)水量與氣象因子、景觀(guān)格局指數(shù)的相關(guān)性分析 由表 7可知：2000—2015年，杭州市年產(chǎn)水量與降水量相關(guān)性最強(qiáng)，相關(guān)系數(shù)為0.959，與氣溫呈極顯著正相關(guān)(P＜0.01)，與實(shí)際蒸散量呈極顯著負(fù)相關(guān)(P＜0.01)。在景觀(guān)格局指數(shù)方面，年產(chǎn)水量與PD相關(guān)性最強(qiáng)，相關(guān)系數(shù)為0.153，與ED、LSI、SHDI呈極顯著正相關(guān)(P＜0.01)，與LPI呈顯著負(fù)相關(guān)(P＜0.05)，與CONTAG呈極顯著負(fù)相關(guān)(P＜0.01)。

表7 杭州市年產(chǎn)水量與氣象因子、景觀(guān)格局指數(shù)的Spearman相關(guān)性分析Table 7 Spearman correlation coefficients among water yield, climate factors and landscape pattern indices

3 討論

InVEST模型水量平衡模型原理認(rèn)為：降水量和實(shí)際蒸散量是估算產(chǎn)水量的2個(gè)重要參數(shù)，其中降水量作為氣象因子的主要變量，由自然條件決定，受人為因素影響較小。也有研究[15]認(rèn)為：景觀(guān)格局比氣象因子對(duì)年產(chǎn)水量的影響更大。本研究發(fā)現(xiàn)：2000—2015年杭州市年產(chǎn)水量與降水量的相關(guān)性大于其與指數(shù)，與顧晉飴等[13]、竇攀烽等[14]結(jié)果一致。氣溫與年產(chǎn)水量具有一定的相關(guān)性，而實(shí)際蒸散量由于計(jì)算方法受到氣象因子、土地利用和景觀(guān)格局的綜合影響，相關(guān)系數(shù)較小。年產(chǎn)水量與各景觀(guān)格局指數(shù)相關(guān)系數(shù)均較小，但顯著性較強(qiáng)，說(shuō)明雖然景觀(guān)格局不是影響杭州市年產(chǎn)水量變化最重要的影響因素，但仍值得關(guān)注。提示不同地區(qū)或生態(tài)系統(tǒng)間影響產(chǎn)水量的決定性因素應(yīng)該是不同的。

研究發(fā)現(xiàn)：2000—2015年間，杭州市年產(chǎn)水量總體上表現(xiàn)為先降低后增加，建設(shè)用地的年產(chǎn)水量高于其他各土地利用類(lèi)型。建設(shè)用地是研究區(qū)變化規(guī)模和速度最快的土地利用類(lèi)型，由不透水性的下墊面組成，植被覆蓋少，不利于有效攔截大量降水，同時(shí)地表蒸散發(fā)能力較弱，人類(lèi)活動(dòng)更為劇烈；建設(shè)用地的增加一定程度上引起產(chǎn)水量的增加[11]。相比之下，水體單位面積產(chǎn)水量最小。研究區(qū)水域廣布(千島湖、西湖、西溪濕地、錢(qián)塘江)，水體下滲程度小于其他土地利用類(lèi)型，但由于水面蒸散速度與強(qiáng)度高[29]，因此水體的年產(chǎn)水量也最小。這也可能是研究區(qū)生態(tài)系統(tǒng)質(zhì)量整體提升所帶來(lái)的效應(yīng)[30]。林地通過(guò)林冠層、枯枝落葉層和土壤層攔截降雨，并通過(guò)蒸騰作用散失水分。快速城市化導(dǎo)致研究區(qū)部分高質(zhì)量城市綠地?fù)p失，景觀(guān)破碎化和異質(zhì)性顯著增大，因此林地提供的產(chǎn)水量不穩(wěn)定[15]；而為建設(shè)綠色城市、生態(tài)城市新增的城市林地大部分處于幼齡階段，這可能是林地產(chǎn)水量在2000—2015年間增長(zhǎng)率大于耕地和草地的原因，也說(shuō)明了城市林地恢復(fù)及提供優(yōu)質(zhì)服務(wù)需要較長(zhǎng)時(shí)間，具有較大的提升空間。

由主城區(qū)梯度變化可知：景觀(guān)格局在距離城市中心10～20 km處破碎化程度尤其嚴(yán)重，表明2000—2015年間此區(qū)域是人類(lèi)活動(dòng)最為密切、景觀(guān)破壞最為嚴(yán)重的區(qū)域。距離城市中心10 km處的單位面積產(chǎn)水量最大，表明這是年產(chǎn)水量在梯度上變化的拐點(diǎn)；氣象因子在該區(qū)域內(nèi)變化不大，表明景觀(guān)格局的變化將重新分配城市產(chǎn)水服務(wù)，對(duì)城市生態(tài)平衡和可持續(xù)發(fā)展產(chǎn)生積極或消極影響[31]。該結(jié)果提示：快速發(fā)展的城市化進(jìn)程中，合理布局規(guī)劃城市用地[32]，可避免因大規(guī)模集中建設(shè)導(dǎo)致的生態(tài)系統(tǒng)服務(wù)功能的降低。

目前InVEST模型產(chǎn)水模塊已被廣泛應(yīng)用，與邱問(wèn)心等[33]研究相比，本研究利用Penman-Monteith公式計(jì)算年蒸散量具有一定的可行性?！逗贾菔兴Y源公報(bào)》數(shù)據(jù)顯示：杭州市2000、2005、2010和2015年實(shí)際產(chǎn)水深度均值分別為0.63、0.54、1.13和1.43 m。本研究結(jié)果對(duì)比驗(yàn)證發(fā)現(xiàn)：4個(gè)年份相對(duì)誤差為5.90%，說(shuō)明模擬結(jié)果相對(duì)可信。但其中部分輸入?yún)?shù)的選擇仍會(huì)影響結(jié)果的準(zhǔn)確性，選擇不同的經(jīng)驗(yàn)公式、空間插值方式和土地利用類(lèi)型分類(lèi)均會(huì)得到不同的結(jié)果。本研究采用研究年份的當(dāng)年氣象數(shù)據(jù)進(jìn)行模型估算，未考慮到氣象因子作用于環(huán)境過(guò)程的時(shí)間滯后性[34]；對(duì)景觀(guān)格局的分析僅考慮了景觀(guān)水平，對(duì)斑塊和類(lèi)型水平的研究有待進(jìn)一步研究。未來(lái)可通過(guò)評(píng)估其他指標(biāo)與景觀(guān)格局的空間相關(guān)性來(lái)進(jìn)一步研究生態(tài)系統(tǒng)服務(wù)與城鄉(xiāng)景觀(guān)格局的關(guān)系，為實(shí)現(xiàn)區(qū)域可持續(xù)規(guī)劃提供科學(xué)依據(jù)。