快速城鎮(zhèn)化對(duì)杭州灣南岸慈溪市土地利用及生態(tài)系統(tǒng)服務(wù)價(jià)值的影響

黃日鵬，李加林?，段義斌，馮佰香，何改麗

(1.寧波大學(xué)地理與空間信息技術(shù)系，315211, 浙江寧波；2.浙江省海洋文化與經(jīng)濟(jì)研究中心，315211，浙江寧波；3.寧?？h發(fā)展與改革局，315699，浙江寧海)

人類生存發(fā)展所需要的資源與空間環(huán)境最終來源于自然生態(tài)系統(tǒng)。生態(tài)系統(tǒng)服務(wù)功能即指通過生態(tài)系統(tǒng)的結(jié)構(gòu)、過程和功能，從而得到生命支持的產(chǎn)品與服務(wù)[1-2]。R. Costanza等[3]在前人研究的基礎(chǔ)上，對(duì)全球16種生物類群的17種生態(tài)系統(tǒng)的公益價(jià)值進(jìn)行估算并繪制全球生態(tài)系統(tǒng)的平均公益價(jià)值表。謝高地等[4]基于文獻(xiàn)調(diào)研、專家知識(shí)、統(tǒng)計(jì)資料和遙感監(jiān)測(cè)等數(shù)據(jù)源，通過模型運(yùn)算和地理信息空間分析等方法，制定出中國(guó)生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子表，在土地利用變化研究方面[5-8]得到廣泛的應(yīng)用。近幾十年來，國(guó)內(nèi)外學(xué)者都對(duì)生態(tài)系統(tǒng)服務(wù)功能研究非常重視，其已然成為生態(tài)科學(xué)研究上的熱點(diǎn)與前沿。

隨著人類社會(huì)經(jīng)濟(jì)的快速發(fā)展，城鄉(xiāng)建設(shè)用地不斷擴(kuò)展，大量占用其他用地，盡管帶來較大的經(jīng)濟(jì)效益，但也導(dǎo)致水土資源流失、生態(tài)系統(tǒng)服務(wù)價(jià)值下降等問題。從長(zhǎng)遠(yuǎn)來看，這可能會(huì)影響城鎮(zhèn)水土保持平衡與區(qū)域社會(huì)經(jīng)濟(jì)的可持續(xù)發(fā)展。不少學(xué)者對(duì)這方面已經(jīng)做了相當(dāng)深入細(xì)致的研究：張修峰等[9]以肇慶仙女湖為例，評(píng)估城市湖泊退化過程中水生態(tài)系統(tǒng)服務(wù)價(jià)值；岳書平等[10]從土地利用變化的角度出發(fā)，選取東北樣帶為研究區(qū)，運(yùn)用GIS和遙感技術(shù)分析近30年來不同類型區(qū)土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響；李屹峰等[11]采用空間顯式的生態(tài)系統(tǒng)服務(wù)功能評(píng)估軟件InVEST中的“產(chǎn)水量”“土壤保持”“水質(zhì)凈化” 模型，研究流域土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)功能的影響；曾杰等[12]采用相關(guān)分析方法分析武漢城市圈不同地域的時(shí)序變化與空間分異特征，探討武漢城市圈生態(tài)系統(tǒng)服務(wù)價(jià)值與不同土地利用類型變化之間的關(guān)聯(lián)特征。但目前仍少有專門針對(duì)快速城鎮(zhèn)化所導(dǎo)致的生態(tài)系統(tǒng)服務(wù)功能變化研究[13]。筆者選取位于杭州灣南岸社會(huì)經(jīng)濟(jì)發(fā)展水平較高、城鎮(zhèn)化速度較快的慈溪市作為研究區(qū)，對(duì)快速城鎮(zhèn)化影響下的生態(tài)服務(wù)價(jià)值變化展開研究，以期為城鎮(zhèn)化建設(shè)及區(qū)域生態(tài)系統(tǒng)服務(wù)功能保育提供科學(xué)依據(jù)。

圖1 研究區(qū)示意圖Fig.1 Schematic diagram of the study area

1 研究區(qū)概況

杭州灣位于中國(guó)浙江省東北部，西起浙江海鹽縣澉浦鎮(zhèn)和上虞區(qū)之間的曹娥江收閘斷面，東至揚(yáng)子角到鎮(zhèn)海角連線，與舟山、北侖港海域?yàn)猷彙：贾轂车靥幈眮啛釒暇?，屬季風(fēng)型氣候，四季分明，冬夏稍長(zhǎng)，春秋略短。平均年日照時(shí)間2 038 h，年日照比例為47%。年平均氣溫16 ℃，雨量充足，年平均降水量1 272.8 mm，平均年徑流總量5.122億m3，降水高峰月為9月，平均占全年降水量的14%。冬季盛行西北至北風(fēng)，夏季盛行東到東南風(fēng)，全年以東風(fēng)為主，100年內(nèi)未發(fā)生大的自然災(zāi)害。慈溪市位于杭州灣南岸，是長(zhǎng)江三角洲經(jīng)濟(jì)帶滬、杭、甬經(jīng)濟(jì)金三角的中心地帶，區(qū)位和交通優(yōu)勢(shì)十分明顯。在杭州灣大橋及舟山跨海大橋的建成通車后，其區(qū)域優(yōu)勢(shì)更加明顯。慈溪是2017年全國(guó)中小城市綜合實(shí)力百?gòu)?qiáng)縣市第5名。2017年12月，當(dāng)選中國(guó)工業(yè)百?gòu)?qiáng)縣。

2 研究方法

2.1 數(shù)據(jù)來源及土地分類

本研究數(shù)據(jù)主要來源于慈溪市2005、2010和2015年3期Landsat TM/OLI遙感影像，軌道號(hào)為118-39。采用ENVI5.2遙感影像處理軟件對(duì)遙感數(shù)據(jù)進(jìn)行大氣校正、幾何精校正、假彩色合成和圖像拼接等數(shù)據(jù)預(yù)處理，然后運(yùn)用慈溪市域界線進(jìn)行影像裁剪，得到研究區(qū)影像數(shù)據(jù)。再對(duì)研究區(qū)3期遙感影像進(jìn)行土地利用類型的目視解譯和人機(jī)交互解譯，得到不同時(shí)期的土地利用數(shù)據(jù)，3年解譯精度均超過0.9，符合研究判別的精度要求。根據(jù)GB/T 21010—2017《土地利用現(xiàn)狀分類》，結(jié)合研究區(qū)實(shí)際情況，將研究區(qū)分為建設(shè)用地、水田、旱地、水體、林地、灘地、未利用地7個(gè)大類。

2.2 土地利用動(dòng)態(tài)度和利用強(qiáng)度

為較為直觀的反映出土地利用類型變化的程度與速度，筆者采用土地利用類型動(dòng)態(tài)度[14]來描述土地類型的變化速度與程度。

土地利用強(qiáng)度能夠顯示出土地利用的廣度和深度，不僅反映了土地利用中土地本身的自然屬性，而且反映人類因素與自然環(huán)境因素的綜合效應(yīng)。筆者根據(jù)王秀蘭等[15]提出的土地利用程度的綜合分析方法，將土地利用程度按照土地自然綜合體在社會(huì)因素影響下的自然平衡狀態(tài)分為若干等級(jí), 不同等級(jí)賦予不同指數(shù)(表1)[16]。

表1 土地利用強(qiáng)度分級(jí)表

2.3 生態(tài)系統(tǒng)服務(wù)價(jià)值

近些年來，我國(guó)研究者基于我國(guó)的現(xiàn)實(shí)情況，提出適用于我國(guó)研究的價(jià)值當(dāng)量因子法以評(píng)估中國(guó)陸地生態(tài)系統(tǒng)服務(wù)價(jià)值。筆者采用謝高地等[17]的研究所得的生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子表，將生態(tài)系統(tǒng)服務(wù)劃分為食物生產(chǎn)、原料生產(chǎn)、水資源供給、氣體調(diào)節(jié)、氣候調(diào)節(jié)、凈化環(huán)境、水文調(diào)節(jié)、土壤保持、維持養(yǎng)分循環(huán)、生物多樣性和美學(xué)景觀等11種服務(wù)功能，進(jìn)而求得慈溪市不同生態(tài)系統(tǒng)單位面積的生態(tài)服務(wù)價(jià)值。根據(jù)謝高地等[17]的研究方法，將1個(gè)標(biāo)準(zhǔn)單位生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子定義為慈溪市每hm2農(nóng)田的年平均自然糧食產(chǎn)量的凈利潤(rùn)量，慈溪市的糧食產(chǎn)量主要以谷物、豆類、薯類為主。我們以谷物、豆類、薯類3類主要糧食總產(chǎn)量與其相應(yīng)的糧食單價(jià)來計(jì)算慈溪市農(nóng)田生態(tài)系統(tǒng)的糧食作物總產(chǎn)值，糧食作物總產(chǎn)值除去播種面積得單位面積糧食作物產(chǎn)值，再考慮在沒有人力投入的自然生態(tài)系統(tǒng)提供的經(jīng)濟(jì)價(jià)值是現(xiàn)有單位面積農(nóng)田提供的食物生產(chǎn)服務(wù)經(jīng)濟(jì)價(jià)值的1/7[18]，求得慈溪市1個(gè)生態(tài)服務(wù)價(jià)值量因子的經(jīng)濟(jì)價(jià)值量為1 991.21元/hm2，得到慈溪市不同生態(tài)系統(tǒng)單位面積的生態(tài)服務(wù)價(jià)值(表2)。最后采用R. Costanza等[3]提出的生態(tài)服務(wù)價(jià)值分析模型，并加以適當(dāng)修訂，計(jì)算慈溪市的生態(tài)服務(wù)價(jià)值，公式如下：

ESVak=SaVCak;

(1)

(2)

(3)

式中：ESVak分別表示該年份第a類土地k項(xiàng)服務(wù)功能價(jià)值，元；ESVa第a類土地的生態(tài)系統(tǒng)服務(wù)價(jià)值，元；ESV為該地區(qū)生態(tài)系統(tǒng)服務(wù)總價(jià)值，元；Sa表示a類型的土地面積，hm2；VCak表示a土類第k種生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)，元/hm2。

表2 慈溪市不同土地利用類型單位面積生態(tài)服務(wù)價(jià)值

注：不同土地類型盡量與謝高地等[17]的當(dāng)量因子表中二級(jí)地類相對(duì)應(yīng)，林地視作針闊混交林，灘地視作濕地，未利用地視作裸地。Notes：Woodland is regarded as mixed forest of coniferous and broad leaved land， tidal flat regarded as wetland, and unused land regarded as bare land for mapping the type of land use to the 2nd grade land in the equivalent factor table by G. Xie[17]

3 結(jié)果

3.1 杭州灣南岸慈溪市土地利用時(shí)空變化特征

結(jié)合表3和圖2可知，在3個(gè)時(shí)期慈溪市各土地利用空間格局情況具有較高的一致性。占地較多的土地利用類型是建筑用地、旱地、林地、灘地，而水田、水體與未利用所占面積相對(duì)較小。不同的土地類型空間分布呈帶狀，其中灘地分布在慈溪市外圍的北部到東部，接壤杭州灣，2015年面積為225 km2，占總面積的17%。灘地以南旱地與建設(shè)用地依次分布，2015年二者合占地為776 km2，合占總面積的58%。林地主要分布于建設(shè)用地以南，面積達(dá)163 km2。水體主要以水庫(kù)或坑塘的形式分布于山地與建設(shè)用地、旱地與灘地之間，少部分以河流網(wǎng)的形式分布于中部平坦地區(qū)；水田主要分布于山地平原接壤地、水分充足的地方；未利用地面積最少，2015年僅1.23 km2。

表3 慈溪市生態(tài)系統(tǒng)服務(wù)價(jià)值變化與城鎮(zhèn)化水平、土地利用強(qiáng)度的關(guān)系

Notes: Ecosystem service value is abbreviated as ESV below.

圖2 杭州灣南岸慈溪市2005、2010和2015年土地利用類型示意圖Fig.2 Schematic diagram of land use types of Cixi county on the south shore of Hangzhou Bay in 2005, 2010, and 2015

在時(shí)間維度上，2005—2015灘地的面積急劇減少，其中2005—2010年面積減少92.94 km2，2010—2015年又減少了59.79 km2，10年年均減少15.3 km2，土地動(dòng)態(tài)度達(dá)4.04%；近些年來人們大量的灘地圍墾是其主要原因。同樣面積減少的林地在2005—2010年減少30 km2，2010—2015年則變化不大。與此相對(duì)應(yīng)，旱地和建設(shè)用地的面積都呈現(xiàn)不同程度的增加，其中建設(shè)用地最為明顯，前5年增加了72.42 km2，后5年增加了40.53 km2，10年土地動(dòng)態(tài)度達(dá)到50%。旱地前5年增加了61.27 km2,后5年反而減少了將近5 km2。水體面積平穩(wěn)增長(zhǎng)，水田呈現(xiàn)先減后增的波動(dòng)，但無論是面積增大的地類還是面積縮減的地類，2010—2015年的土地利用變化程度普遍比2005—2010年階段更為和緩。3年的土地利用強(qiáng)度分別為241、261和273。這一趨勢(shì)也表明城鎮(zhèn)化雖然持續(xù)進(jìn)行，但2010—2015年速度相對(duì)放緩。這可能與產(chǎn)業(yè)的結(jié)構(gòu)優(yōu)化升級(jí)、人們保護(hù)生態(tài)的意識(shí)增強(qiáng)有關(guān)。

3.2 快速城鎮(zhèn)化對(duì)慈溪市土地利用的影響

為量化分析快速城鎮(zhèn)化對(duì)慈溪市的生態(tài)系統(tǒng)服務(wù)價(jià)值的影響，筆者基于3期的土地利用矢量數(shù)據(jù)，統(tǒng)一以2015年的海岸線為準(zhǔn)，通過ArcGIS 10.2的空間分析功能對(duì)各年土地利用矢量數(shù)據(jù)進(jìn)行疊加分析，得到3期的土地利用轉(zhuǎn)移矩陣。

分析轉(zhuǎn)移矩陣可知，土地利用各類型間的多向轉(zhuǎn)化導(dǎo)致了土地利用類型面積在空間分布上的變化。在2005—2010年期間，慈溪市分別有38.28 km2旱地、13.18 km2水田、7.59 km2林地、4.94 km2水體、27.54 km2灘地、1.99 km2未利用地轉(zhuǎn)換為建設(shè)用地。2010年到2015年期間又有18.93 km2旱地、1.85 km2水田、1.48 km2林地、1.95 km2水體、16.85 km2灘地、0.03 km2未利用地轉(zhuǎn)換為建設(shè)用地。從中可以看出，其建設(shè)用地主要由旱地、灘地轉(zhuǎn)變而來，在2005—2010年間兩者轉(zhuǎn)換為建設(shè)用地面積合計(jì)為65.82 km2，占新增建設(shè)用地面積的70.37%，在2010—2015年間兩者轉(zhuǎn)變?yōu)榻ㄔO(shè)用地面積為35.78 km2，占87.05%。從2005—2015年，土地利用由水田轉(zhuǎn)換為建設(shè)用地的地區(qū)主要分布于杭州灣南岸東南部的掌起鎮(zhèn)、范市鎮(zhèn)、三北鎮(zhèn)與慈溪市外圍。中部地區(qū)各鎮(zhèn)每年都有少量旱地轉(zhuǎn)換為建筑用地，這是中部各鎮(zhèn)城鎮(zhèn)化擴(kuò)張的體現(xiàn)。北部的寧波杭州灣新區(qū)與慈東經(jīng)濟(jì)開發(fā)區(qū)這2個(gè)開發(fā)區(qū)周邊灘地大面積的轉(zhuǎn)變?yōu)榻ㄔO(shè)用地。杭州灣大橋的建造使得這一地區(qū)具有顯著地地理區(qū)位優(yōu)勢(shì)，加速了此地區(qū)城鎮(zhèn)化的進(jìn)程。

使用ArcGIS10.2軟件將慈溪市土地利用類型矢量圖轉(zhuǎn)換為柵格圖片，再采用Fragstats軟件計(jì)算3個(gè)時(shí)期景觀水平格局指數(shù)。結(jié)果顯示：城鎮(zhèn)化過程導(dǎo)致建設(shè)用地與其他土地利用類型的景觀格局變化明顯。在2005—2015年段，建設(shè)用地的斑塊數(shù)量從713減少至685，同時(shí)，旱地、灘地的斑塊數(shù)量則呈現(xiàn)逐年上升的趨勢(shì)，這是建設(shè)用地?cái)U(kuò)張融合的同時(shí)分割旱地與灘地的結(jié)果。占地面積相對(duì)較少的林地、水田、水體斑塊數(shù)量從2005—2015年都先是降低而后增加。在平均斑塊面積這一指標(biāo)上，10年間建設(shè)用地面積增大了56%。灘地的平均斑塊面積大幅減少，2005年單個(gè)斑塊平均面積為9.18 km2，2010年減少至6.29 km2，2015年僅有3.32 km2。該指數(shù)的大幅縮減是圍墾灘地的結(jié)果，圍墾一方面使灘地總面積大幅縮減，另一方面大的灘地斑塊遭到割裂，斑塊數(shù)量增加，使得灘地平均斑塊面積驟減。景觀邊界密度能直觀地反映景觀或景觀類型邊界被割裂的程度，同時(shí)反映景觀的破碎化程度。10年間，旱地、建設(shè)用地、灘地、水體的邊界密度在人類活動(dòng)的影響下，都有不同程度的增加，最為明顯的是建筑用地與水體，分別增加了26.37%與30.5%，表明城鎮(zhèn)化過程中建設(shè)用地的擴(kuò)張使土地破碎化程度增加。

3.3 杭州灣南岸慈溪市的生態(tài)系統(tǒng)服務(wù)價(jià)值變化

圖3 慈溪市不同土地利用類型的生態(tài)系統(tǒng)服務(wù)價(jià)值及其變化Fig.3 Ecosystem service value of different land use types and its changes in Cixi county

由圖3可知，隨著杭州灣南岸慈溪市土地利用的轉(zhuǎn)變，該地區(qū)的生態(tài)系統(tǒng)服務(wù)總價(jià)值與構(gòu)成也相應(yīng)發(fā)生變化。總體而言慈溪市生態(tài)服務(wù)價(jià)值呈現(xiàn)下降趨勢(shì)，2005年總生態(tài)服務(wù)價(jià)值達(dá)65.1億元，2010年減少為59.6億元，2015年又減少至55.1億元。慈溪市各類土地地類中，灘地所具有的生態(tài)系統(tǒng)服務(wù)價(jià)值明顯高于其他地類，2005年灘地生態(tài)系統(tǒng)服務(wù)價(jià)值為39.2億元，占地面積為28.4%的灘涂占總生態(tài)系統(tǒng)服務(wù)價(jià)值的60.2%，可見灘地對(duì)于慈溪市大氣調(diào)節(jié)、水文調(diào)節(jié)、土壤物質(zhì)循環(huán)、維持生物多樣性等方面的重要程度。近些年來的圍墾灘涂確實(shí)帶來了不少的經(jīng)濟(jì)利益，但在無形中也造成巨大的生態(tài)損失。如何權(quán)衡利弊，進(jìn)行科學(xué)可持續(xù)的開發(fā)建設(shè)將是當(dāng)下甚至于未來杭州灣南岸地區(qū)發(fā)展急需解決的問題。水體的生態(tài)服務(wù)價(jià)值量?jī)H次于灘地。在整體生態(tài)服務(wù)價(jià)值量下降的情況下，水體的價(jià)值量不降反升，從2005年的13.2億元到2010年的18.4億元，上升至2015年的20.1億元。林地與旱地緊跟其后，其原因有所不同，林地雖占地面積不算大，但單位面積價(jià)值量大，而旱地占地面積大彌補(bǔ)了單位面積的低價(jià)值量，兩者2015年的生態(tài)系統(tǒng)服務(wù)價(jià)值量分別為7.53億元和3.49億元。水田的價(jià)值量呈現(xiàn)波折變化，3個(gè)時(shí)間段分別占總生態(tài)服務(wù)價(jià)值的1.1%、0.8%和1.2%。建筑用地與未利用地提供的價(jià)值量可以忽略。

從不同年份的價(jià)值變化來看，價(jià)值量最大的灘地對(duì)整體價(jià)值量的變化影響巨大。2005—2010年減少9.6億元，2010—2015年又減少6.2億元，10年間灘地的生態(tài)服務(wù)價(jià)值共減少15.8億元。除灘地外變化量最大的是水體，水體生態(tài)服務(wù)價(jià)值在2005—2010年、2010—2015年分別增加5.2億元和1.7億元。林地生態(tài)服務(wù)價(jià)值10年間減少1.5億元，占變化總量的14.69%。旱地與水田相對(duì)較為平穩(wěn)，旱地的生態(tài)服務(wù)價(jià)值從2005年的3.0億元到2010年的3.5億元，增加0.5億元，在2010—2015年又減少394萬元；水田的生態(tài)服務(wù)價(jià)值2005—2010年減少0.2億元，2010—2015年增加0.15億元。

3.4 快速城鎮(zhèn)化影響下的生態(tài)系統(tǒng)服務(wù)價(jià)值

筆者以建設(shè)用地在所有土類中所占比例作為該區(qū)域的城鎮(zhèn)化水平指數(shù)，結(jié)合慈溪市各年的生態(tài)服務(wù)價(jià)值、土地利用強(qiáng)度變化分析，得到圖4。結(jié)合圖4和表3來看，慈溪市的生態(tài)系統(tǒng)服務(wù)價(jià)值量和城鎮(zhèn)化進(jìn)程、土地利用強(qiáng)度呈現(xiàn)負(fù)相關(guān)關(guān)系。在本次研究時(shí)段內(nèi)，慈溪市的城鎮(zhèn)化水平指數(shù)從2005年的17.0%升至2010年22.4%，2015年達(dá)到25.5%，呈逐年上升的趨勢(shì)，年均變率0.85%。生態(tài)系統(tǒng)服務(wù)價(jià)值的變化趨勢(shì)與此相反，2005—2010年、2010—2015年2個(gè)時(shí)期，慈溪市的生態(tài)系統(tǒng)服務(wù)價(jià)值分別減少5.5億元和4.5億元，減幅有所縮小，10年減少將近15.4%，年均達(dá)1.5%。土地利用強(qiáng)度10年內(nèi)增加了13.3%，其中2005—2010年階段變化最為顯著，年均增加1.6%，2010—2015年年均增幅為0.9%。

圖4 各時(shí)間慈溪市生態(tài)系統(tǒng)服務(wù)價(jià)值與城鎮(zhèn)化 水平的關(guān)系Fig.4 Relationship between ecosystem service value and urbanization level in Cixi county at various periods

為量化分析快速城鎮(zhèn)化對(duì)整個(gè)區(qū)域的生態(tài)系統(tǒng)服務(wù)價(jià)值的影響，筆者將因其他類型土地轉(zhuǎn)變?yōu)榻ㄔO(shè)用地而損失的生態(tài)系統(tǒng)服務(wù)價(jià)值視作為快速城鎮(zhèn)化直接導(dǎo)致的生態(tài)系統(tǒng)服務(wù)價(jià)值流失量。由土地利用轉(zhuǎn)移矩陣與生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量表(表2)，得出表4。由表4可知：2005—2010年因快速城鎮(zhèn)化損失的總生態(tài)系統(tǒng)服務(wù)價(jià)值達(dá)4.8億元，占該時(shí)段生態(tài)系統(tǒng)服務(wù)價(jià)值流失的88%。2010—2015年損失量少于上個(gè)時(shí)段，但也達(dá)到2.5億元，占55%。從2005—2015年10年內(nèi)損失價(jià)值量達(dá)7.3億元，占73%。以上數(shù)據(jù)可以說明，快速城鎮(zhèn)化下建設(shè)用地?cái)U(kuò)張是該地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值減少的主要原因。

表4 快速城鎮(zhèn)化直接影響下的生態(tài)系統(tǒng)服務(wù)價(jià)值變化

建設(shè)用地對(duì)各類用地的占用，對(duì)慈溪市生態(tài)服務(wù)價(jià)值的減少有不同程度的影響。分時(shí)段而言，2005—2010年與2010—2015年2個(gè)時(shí)段建設(shè)用地所占用的各土地利用類型的面積之間比例相近，但前一時(shí)段面積數(shù)值普遍大于后一時(shí)段，即前一時(shí)段的建設(shè)用地?cái)U(kuò)張強(qiáng)度高于后一時(shí)段。由表4可以看出，在2個(gè)時(shí)段中，對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值影響最大的是建設(shè)過程中對(duì)灘地的圍墾，2005—2010年和2010—2015年減少的價(jià)值量分別為2.85億元和1.74億元，分別占據(jù)因快速城鎮(zhèn)化影響下總流失價(jià)值量的58.9%和70.8%。這是因?yàn)檎紦?jù)灘地面積大，且灘地的生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)很高。同樣，被侵占的水體面積雖然所占的比例不高，2個(gè)時(shí)段皆為5%；但由于其極高的生態(tài)系統(tǒng)服務(wù)價(jià)值，2個(gè)時(shí)段損失的生態(tài)系統(tǒng)服務(wù)價(jià)值也達(dá)到1.2億元和0.5億元，分別占25%和20%。旱地、林地生態(tài)系統(tǒng)服務(wù)價(jià)值所占比例都不高，原因有所不同。旱地被侵占的面積是最大的，但由于生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)很低，所以流失的價(jià)值量也不高，2005—2010年減少了0.31億元，2010—2015年減少了0.15億元。林地的情況相反，生態(tài)服務(wù)價(jià)值系數(shù)高而面積變化不大，2個(gè)時(shí)段林地占據(jù)的面積分別為8.1%和3.6%，損失的價(jià)值量分別為0.34億元和0.07億元。流失價(jià)值量比重最低的地類是水田與未利用地，10年內(nèi)水田轉(zhuǎn)化為建設(shè)用地而損失的價(jià)值量為0.12億元，而未利用地?cái)?shù)值可以忽略。

4 結(jié)論與討論

1)杭州灣南岸慈溪市近些年來的土地利用變化主要體現(xiàn)在灘地面積大量減少，10年減少將近40.4%，其次為林地，10年間也減少大約32 km2，灘地遭到圍墾后多數(shù)轉(zhuǎn)變?yōu)楹档?，少量轉(zhuǎn)變?yōu)樗?。旱地的面積持續(xù)增長(zhǎng)，10年共增加56 km2；建設(shè)用地面積增量最為顯著，從2005—2015年年均增加11.3 km2，其中大部分由旱地、灘地轉(zhuǎn)變而來，建設(shè)用地的擴(kuò)張使其他用地的破碎程度增加。

2)2005—2015年10年間慈溪市生態(tài)系統(tǒng)服務(wù)價(jià)值持續(xù)下降，達(dá)10.02億元，其主要的原因是灘地面積的銳減，水體的生態(tài)服務(wù)價(jià)值不降反升；慈溪市各種土地利用類型中，灘地的生態(tài)服務(wù)價(jià)值遠(yuǎn)高于其他地類，2005年灘地生態(tài)服務(wù)價(jià)值占整個(gè)慈溪市地區(qū)的60.2%，其次為水體、林地、旱地。

3)慈溪市的城鎮(zhèn)化水平與生態(tài)系統(tǒng)服務(wù)價(jià)值呈負(fù)相關(guān)關(guān)系。城鎮(zhèn)化直接影響下的建設(shè)用地?cái)U(kuò)張是生態(tài)系統(tǒng)服務(wù)價(jià)值流失的主要原因，10年間累積損失達(dá)7.3億元，占總減少量的73%，且2005—2010年各類土地生態(tài)服務(wù)價(jià)值的損失量都要高于2010—2015年。價(jià)值量損失結(jié)構(gòu)中，10年內(nèi)因?yàn)┑剞D(zhuǎn)變?yōu)榻ㄔO(shè)用地而損失的價(jià)值量達(dá)到4.6億元，占62.9%，其次為被占用的水體，價(jià)值損失量占23.6%，其余土地類型影響相對(duì)較小。

4)本研究結(jié)果的誤差主要可能來自以下幾個(gè)方面：首先是遙感影像解譯的準(zhǔn)確率，其次在計(jì)算慈溪市單位面積的生態(tài)服務(wù)價(jià)值時(shí)所采用的慈溪市統(tǒng)計(jì)年鑒、各種作物的平均市場(chǎng)價(jià)等數(shù)據(jù)精確程度有待提高。同時(shí)，本文采用的是謝高地等[17]2015年最新修訂的生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量表，在準(zhǔn)確性上有所提高；但總的來說，當(dāng)前我們?nèi)圆荒軌蛘J(rèn)清生態(tài)系統(tǒng)所有的服務(wù)的價(jià)值，對(duì)生態(tài)系統(tǒng)服務(wù)功能價(jià)值評(píng)估的體系還未完善，存在著許多問題仍有待解決，因此結(jié)果僅為最低保守值。