朱從謀，李武艷，杜瑩瑩，許紅衛(wèi)，王 珂

浙江省耕地多功能價(jià)值時(shí)空變化與權(quán)衡-協(xié)同關(guān)系

朱從謀1，李武艷2，杜瑩瑩3，許紅衛(wèi)1，王 珂1※

（1. 浙江大學(xué)環(huán)境與資源學(xué)院，杭州 310058；2. 浙江財(cái)經(jīng)大學(xué)土地與城鄉(xiāng)發(fā)展研究院，杭州 310018；3. 浙江工商大學(xué)公共管理學(xué)院，杭州 310018）

耕地多功能價(jià)值的時(shí)空演變與權(quán)衡-協(xié)同關(guān)系測度對(duì)于深化耕地可持續(xù)利用和管理具有重要意義。該研究以經(jīng)濟(jì)快速發(fā)展地區(qū)浙江省為例，采用價(jià)值量化方法對(duì)2000、2010和2015年浙江省耕地多功能進(jìn)行價(jià)值評(píng)估并分析其時(shí)空變化特征。運(yùn)用Spearman秩相關(guān)系數(shù)法、雙變量空間自相關(guān)模型以及可拓展隨機(jī)性環(huán)境影響評(píng)估（Stochastic Impacts by Regression on Population，Affluence and Technology，STIRPAT）模型探究耕地多功能之間的權(quán)衡-協(xié)同關(guān)系及影響因素。結(jié)果表明：1）2000—2015年浙江省耕地多功能總價(jià)值下降，地均耕地多功能價(jià)值呈現(xiàn)先下降后上升的趨勢。其中氣體調(diào)節(jié)、水源涵養(yǎng)和社會(huì)保障功能價(jià)值下降明顯，食物生產(chǎn)和美學(xué)景觀功能價(jià)值增加；耕地多功能總價(jià)值在空間上呈現(xiàn)北高南低的分布格局，2000—2015年大部分縣市耕地多功能總價(jià)值均有不同程度的下降，西南山地丘陵區(qū)耕地多功能總價(jià)值有所提升。2）浙江省耕地多功能之間主要表現(xiàn)為協(xié)同關(guān)系，各項(xiàng)功能之間的協(xié)同-權(quán)衡關(guān)系存在空間異質(zhì)性；2000—2015年，浙江省耕地多功能協(xié)同關(guān)系總體減弱，在空間上主要表現(xiàn)為高值協(xié)同區(qū)減少。3）城鎮(zhèn)居民可支配收入和地均農(nóng)業(yè)機(jī)械總動(dòng)力對(duì)耕地多功能總價(jià)值有負(fù)面影響，農(nóng)村居民可支配收入增加有助于耕地多功能總價(jià)值的提升。該研究成果可為科學(xué)劃定耕地利用與保護(hù)區(qū)，促進(jìn)耕地多功能的協(xié)同利用和提升耕地資源價(jià)值提供科學(xué)依據(jù)。

土地利用；模型；耕地；權(quán)衡；協(xié)同；多功能；浙江省

0 引 言

耕地系統(tǒng)是在自然基礎(chǔ)上經(jīng)人工改造而形成的集約化生產(chǎn)系統(tǒng)，具有典型的多功能特性[1]。耕地不僅提供著全世界66%的糧食等物質(zhì)供給，同時(shí)還具有包括氣候調(diào)節(jié)、環(huán)境凈化、維持生物多樣性與水源涵養(yǎng)等生態(tài)功能，保障農(nóng)民就業(yè)、維護(hù)國家糧食安全等社會(huì)功能以及保持農(nóng)耕文化、提供開敞空間景觀等文化功能[2-3]。隨著城鎮(zhèn)化的快速推進(jìn)，耕地資源的稀缺性日益顯現(xiàn)，人們對(duì)耕地的認(rèn)知和需求從單一的生產(chǎn)功能逐漸向系統(tǒng)化的多元功能轉(zhuǎn)變，耕地資源的非市場價(jià)值進(jìn)一步凸顯[4]。由于耕地功能類型的多樣性、空間分布不均衡性以及人類利用的選擇性，耕地各功能間存在此消彼長的權(quán)衡與相互促進(jìn)或抑制的協(xié)同關(guān)系[5-6]。在耕地利用和保護(hù)過程中考慮耕地多元功能之間的權(quán)衡-協(xié)同關(guān)系將有效規(guī)避耕地利用負(fù)向效應(yīng)，實(shí)現(xiàn)耕地利用效益最大化。2018年中國中央一號(hào)文件提出“大力開發(fā)農(nóng)業(yè)多種功能，建設(shè)一批設(shè)施完備、功能多樣的休閑觀光園區(qū)等”，這表明耕地多功能利用成為提升農(nóng)業(yè)發(fā)展質(zhì)量、促進(jìn)鄉(xiāng)村產(chǎn)業(yè)振興的重要途徑。因此，正確、全面地認(rèn)識(shí)和評(píng)估耕地多功能價(jià)值并掌握耕地多功能之間權(quán)衡-協(xié)同演變規(guī)律，對(duì)于促進(jìn)耕地可持續(xù)利用和推進(jìn)鄉(xiāng)村振興具有重要意義。

“多功能”概念最早來源于農(nóng)業(yè)領(lǐng)域，2000年以來，多功能領(lǐng)域的研究逐漸拓展到景觀管理[7]、生態(tài)系統(tǒng)服務(wù)[8]、土地利用以及鄉(xiāng)村發(fā)展等領(lǐng)域[9-10]。從國外研究看，學(xué)者多從村域、農(nóng)場尺度對(duì)農(nóng)業(yè)或土地多能性進(jìn)行分析，主要集中在農(nóng)業(yè)和土地多功能性概念辨析與分類、耕地多功能評(píng)價(jià)與模式識(shí)別、農(nóng)業(yè)多功能價(jià)值評(píng)估等[11-12]。近年來，隨著中國城鄉(xiāng)重構(gòu)與轉(zhuǎn)型速度加快，耕地多功能利用與管理受到高度重視。國內(nèi)學(xué)者的研究主要集中在耕地多功能內(nèi)涵[1,3]、功能評(píng)價(jià)[13-14]、功能管理等方面[15]。在耕地多功能分類方面，多數(shù)學(xué)者從土地生產(chǎn)、生態(tài)和社會(huì)3類基礎(chǔ)功能進(jìn)行細(xì)化和拓展[14,16]。隨著社會(huì)經(jīng)濟(jì)不斷發(fā)展，耕地文化功能價(jià)值逐漸受到學(xué)者和社會(huì)關(guān)注[17]。在耕地多功能評(píng)價(jià)研究中，目前常用的耕地多功能度量方法包括功能指標(biāo)體系表征法、物質(zhì)量法、能值法和價(jià)值量法[9,18-19]。其中，價(jià)值量法從土地價(jià)值角度綜合測算耕地多功能價(jià)值，可以直觀地反映耕地功能價(jià)值結(jié)構(gòu)、大小和演變趨勢，有利于全面認(rèn)識(shí)耕地資源的總體價(jià)值，對(duì)真實(shí)反映資源稀缺和資源價(jià)值的土地要素價(jià)格形成機(jī)制以及制定耕地保護(hù)的經(jīng)濟(jì)補(bǔ)償政策具有一定借鑒意義[20-21]。隨著中國生態(tài)文明建設(shè)戰(zhàn)略的實(shí)施，面臨日益稀缺的耕地資源，科學(xué)認(rèn)知并有效協(xié)同與權(quán)衡耕地利用的多功能性對(duì)于耕地保護(hù)和生態(tài)治理具有重要意義。目前有關(guān)耕地多功能權(quán)衡-協(xié)同的研究主要是借鑒生態(tài)系統(tǒng)服務(wù)權(quán)衡-協(xié)同研究的分析框架[22-23]。近年來，生態(tài)系統(tǒng)服務(wù)之間的關(guān)系研究逐漸從定性到定量分析、線性到非線性關(guān)系識(shí)別進(jìn)行拓展，其中定量方法包括相關(guān)性分析、線性回歸、均方根誤差以及耦合協(xié)調(diào)度等多種分析方法[22,24-26]，對(duì)研究耕地多功能權(quán)衡-協(xié)同關(guān)系具有重要借鑒意義。但目前有關(guān)快速發(fā)展地區(qū)耕地多功能之間交互關(guān)系的空間差異及動(dòng)態(tài)變化研究仍較少[6]，探討快速發(fā)展地區(qū)耕地多功能價(jià)值變化的社會(huì)經(jīng)濟(jì)驅(qū)動(dòng)因素的分析欠缺?？赏卣闺S機(jī)性環(huán)境影響評(píng)估（Stochastic Impacts by Regression on Population，Affluence and Technology，STIRPAT）模型是一個(gè)應(yīng)用廣泛、較為成熟的環(huán)境壓力評(píng)價(jià)模型，可定量分析各類社會(huì)經(jīng)濟(jì)因素對(duì)生態(tài)環(huán)境壓力造成的非比例影響[27]。在快速城鎮(zhèn)化進(jìn)程中，耕地多功能及其價(jià)值易受到城鄉(xiāng)社會(huì)經(jīng)濟(jì)發(fā)展的廣泛影響，采用STIRPAT模型可進(jìn)一步量化社會(huì)經(jīng)濟(jì)因素對(duì)耕地系統(tǒng)壓力的影響。

浙江省是中國東部沿海發(fā)達(dá)省份，耕地資源相對(duì)稀缺，快速工業(yè)化和城鎮(zhèn)化進(jìn)程導(dǎo)致耕地規(guī)模減少和生態(tài)質(zhì)量下降，亟待充分發(fā)揮耕地多功能，提升農(nóng)業(yè)發(fā)展質(zhì)量。該研究以浙江省67個(gè)縣市區(qū)為分析單元，運(yùn)用價(jià)值量化模型定量評(píng)估耕地多功能價(jià)值，分析浙江省2000—2015年耕地多功能價(jià)值時(shí)空變化。并運(yùn)用Spearman秩相關(guān)分析方法和雙變量空間自相關(guān)模型探討耕地多功能之間權(quán)衡-協(xié)同關(guān)系的空間格局及其動(dòng)態(tài)變化，同時(shí)采用STIRPAT模型分析影響耕地多功能價(jià)值演變的社會(huì)經(jīng)濟(jì)驅(qū)動(dòng)因素。研究結(jié)果可為了解浙江省耕地多功能價(jià)值的時(shí)空演變及制定區(qū)域耕地多功能利用和管理政策提供科學(xué)依據(jù)。

1 研究區(qū)概況與研究方法

1.1 研究區(qū)概況

浙江省（118°01'～123°10'N，27°02'～31°11'E）位于中國東南沿海地區(qū)，陸域面積10.55萬km2，地形地貌復(fù)雜多樣，地勢由西南向東北傾斜（圖1），是中國陸域面積最小的省份之一。浙江省素有“七山一水二分田”之說，耕地主要分布在浙北杭嘉湖、寧紹平原、浙中金衢盆地以及浙東沿海區(qū)，但由于人口眾多，人均耕地僅400 m2左右，不及全國人均耕地水平的一半[28]。同時(shí)，浙江省是中國經(jīng)濟(jì)最活躍的省份之一，2015年常住人口城鎮(zhèn)化率達(dá)到65.8%，人均GDP達(dá)到4.29萬元?？焖俚某擎?zhèn)化、工業(yè)化背景下，浙江省耕地面積快速減少，據(jù)統(tǒng)計(jì)數(shù)據(jù)顯示，2000—2015年浙江省耕地共減少21.21萬hm2，且主要集中在浙北平原和沿海地區(qū)。城鎮(zhèn)空間擴(kuò)張、糧食安全保障、城鄉(xiāng)居民高質(zhì)量生活需求等對(duì)耕地利用和管理提出了嚴(yán)峻挑戰(zhàn)。

圖1 研究區(qū)概況

1.2 研究方法

1.2.1 耕地關(guān)鍵功能識(shí)別

作為一種半自然半人工的復(fù)合系統(tǒng)，耕地的本質(zhì)是滿足人類生存發(fā)展的需求[5]。在耕地功能分類中，目前大部分學(xué)者將其劃分為生產(chǎn)、生態(tài)、社會(huì)和文化功能[13]。結(jié)合研究區(qū)現(xiàn)狀和數(shù)據(jù)可獲取性，本研究選取目前受學(xué)者較多關(guān)注的8種耕地功能進(jìn)行研究[21,29]，各功能的分類、估算方法及數(shù)據(jù)來源見表1，其中耕地?cái)?shù)據(jù)來自于生態(tài)環(huán)境部提供的Landsat TM影像（30 m）解譯數(shù)據(jù)。耕地的生產(chǎn)功能指通過耕地利用為社會(huì)提供農(nóng)產(chǎn)品供給，主要表現(xiàn)為耕地的食物生產(chǎn)；生態(tài)功能則從耕地生態(tài)系統(tǒng)服務(wù)角度出發(fā)，主要包括支持、調(diào)節(jié)等生態(tài)服務(wù)，同時(shí)，長期以來由于人類過度利用耕地進(jìn)而造成農(nóng)業(yè)污染，也給人類社會(huì)和其他生態(tài)系統(tǒng)帶來了負(fù)面影響[19]，因此，該研究將耕地生態(tài)功能劃分為支持服務(wù)功能（氣體調(diào)節(jié)與環(huán)境凈化）和調(diào)節(jié)服務(wù)功能（水源涵養(yǎng)與維持生物多樣性）以及農(nóng)業(yè)污染負(fù)向功能[30]；耕地社會(huì)功能主要表現(xiàn)為耕地利用保障農(nóng)民最低生活水平的能力[31]，本研究主要關(guān)注耕地生活保障功能；耕地文化功能則體現(xiàn)為耕地為人類提供美學(xué)景觀的功能[32]。在識(shí)別耕地食物生產(chǎn)、氣體調(diào)節(jié)、環(huán)境凈化、水源涵養(yǎng)、生物多樣性、農(nóng)業(yè)污染、社會(huì)保障和美學(xué)景觀8種功能基礎(chǔ)上，采用價(jià)值量化模型定量評(píng)估耕地各項(xiàng)功能價(jià)值。

1.2.2 耕地多功能價(jià)值估算

1）食物生產(chǎn)功能價(jià)值

食物生產(chǎn)功能是耕地系統(tǒng)的基礎(chǔ)功能。浙江省種植業(yè)產(chǎn)出主要包括水稻、小麥、玉米、甘薯、大豆、油菜、棉花、茶葉和蔬菜等9種農(nóng)作物，占農(nóng)作物種植總面積的90%以上。該部分價(jià)值測算采用市場價(jià)值法[19]，計(jì)算式為

式中CP為縣市的耕地食物生產(chǎn)價(jià)值，元；P是農(nóng)作物的市場價(jià)格，元/kg；C是農(nóng)作物種植成本，元/kg；y是作物單位面積產(chǎn)出，kg/hm2；s是作物的種植面積，hm2。

表1 耕地8種功能分類及定量估算方法

2）氣體調(diào)節(jié)功能價(jià)值

氣體調(diào)節(jié)功能是農(nóng)作物將吸收的CO2通過光合作用轉(zhuǎn)換成O2，調(diào)節(jié)大氣中的CO2和O2的平衡。根據(jù)已有研究成果[19]，1 kg植物可固定1.63 kg CO2和釋放1.2 kg O2。該部分運(yùn)用成本法進(jìn)行價(jià)值測算，計(jì)算式為

式中AR是縣市的耕地氣體調(diào)節(jié)價(jià)值，元；NPP是年均作物總凈生產(chǎn)力，kg；1是工業(yè)固碳成本，為260.9元/t，2是氧氣釋放成本，為376.47元/t，采用造林成本法獲取，參考鄒月等[33]研究成果；h是作物的經(jīng)濟(jì)系數(shù)，%；r是農(nóng)作物的含水率，%，參考孔凡靖等[19]研究成果。

3）環(huán)境凈化功能價(jià)值

農(nóng)作物在成長過程中具有較為顯著的環(huán)境清潔作用，主要是對(duì)SO2、NO2、HF的吸收和滯塵作用[21]。其計(jì)算式為

式中EP是縣市的耕地環(huán)境凈化價(jià)值，元；Q是單位面積吸收的各種廢氣的質(zhì)量，其中吸收SO2為2.45 kg/hm2，吸收NO2為33.31 kg/hm2，吸收HF為0.33 kg/hm2，吸收粉塵為1 500 kg/hm2，參考Yu等[31]研究成果；C是各種廢氣的平均治理成本，其中SO2治理成本為0.6元/kg，NO2治理成本為0.6元/kg，HF治理成本為0.9元/kg，粉塵治理成本為0.17元/kg，參考Yu等[31]研究成果；是農(nóng)作物總種植面積，hm2。

4）水源涵養(yǎng)功能價(jià)值

耕地水源涵養(yǎng)功能以植物分流截留以及土壤儲(chǔ)水為主。采用綜合蓄水能力來測算耕地系統(tǒng)水源涵養(yǎng)量[6]，運(yùn)用影子工程法來評(píng)價(jià)農(nóng)業(yè)水源涵養(yǎng)價(jià)值，計(jì)算式為

式中WC是縣市的耕地水源涵養(yǎng)價(jià)值，元；是降水量，mm；s是農(nóng)作物的種植面積，hm2；k是農(nóng)作物的截留率，%，參考元媛等[21]研究成果；r是第種土壤面積，hm2；m是第種土壤的粗孔隙率（非毛管孔隙度），%，參考朱慶瑩等[6]研究成果；n是第種土壤的深度，mm；3是人工蓄水成本，元/m3。

5）生物多樣性功能價(jià)值

耕地生態(tài)系統(tǒng)中的生物多樣性是產(chǎn)生并維持其他生態(tài)系統(tǒng)服務(wù)功能的基礎(chǔ)，對(duì)維持生態(tài)系統(tǒng)結(jié)構(gòu)穩(wěn)定與服務(wù)的可持續(xù)供應(yīng)具有重要意義。本研究借鑒謝高地等[32]研究成果對(duì)耕地生物多樣性功能價(jià)值進(jìn)行測算，計(jì)算式為

式中BI是縣市的耕地維持生物多樣性功能價(jià)值，元；指的是有關(guān)耕地生物多樣性當(dāng)量因子，其中水田的生物多樣性價(jià)值當(dāng)量因子1取值為0.21，旱地生物多樣性價(jià)值當(dāng)量因子2取值為0.13；1和2分別為水田和旱地面積，hm2；E指耕地食物生產(chǎn)的單位市場價(jià)值，元/hm2；是耕地總面積，hm2。

6）農(nóng)業(yè)污染負(fù)功能價(jià)值

農(nóng)業(yè)污染指的是在耕地生產(chǎn)過程中由于過量使用化肥、農(nóng)藥而造成的對(duì)農(nóng)田生態(tài)系統(tǒng)破壞。未被土地轉(zhuǎn)化為肥力或植物吸收的化肥、農(nóng)藥部分是直接經(jīng)濟(jì)損失，同時(shí)農(nóng)藥化肥在土壤、糧食中殘留會(huì)降低糧食產(chǎn)品質(zhì)量，導(dǎo)致食物產(chǎn)出損失。據(jù)此，采用化肥、農(nóng)藥利用率、糧食產(chǎn)量損失率和市場價(jià)格法測算農(nóng)業(yè)污染負(fù)價(jià)值[19]。具體如下：

式中AP指縣市的耕地農(nóng)業(yè)污染負(fù)價(jià)值；是化肥、農(nóng)藥使用總量，t；指化肥、農(nóng)藥的利用率，%，為由于化肥農(nóng)藥而造成的食物產(chǎn)量損失率，%，參考鄧明君等[34]研究成果；指化肥、農(nóng)藥的市場價(jià)格，元/kg。

7）社會(huì)保障功能價(jià)值

耕地的社會(huì)保障功能是指耕地為農(nóng)民提供生活保障，維護(hù)社會(huì)經(jīng)濟(jì)正常發(fā)展。社會(huì)保障價(jià)值以農(nóng)戶最低生活保障水平進(jìn)行間接評(píng)估[31]。計(jì)算式如下：

式中SS是縣市的耕地社會(huì)保障價(jià)值，元；1是農(nóng)業(yè)人口數(shù)量，人；是各縣（市、區(qū)）最低生活保障標(biāo)準(zhǔn)，元/人；1是農(nóng)村人均可支配收入，元；2是城鎮(zhèn)人均可支配收入，元。

8）美學(xué)景觀功能價(jià)值

耕地美學(xué)景觀功能是為人們提供審美享受的功能。本研究基于謝高地等[32]已有研究，結(jié)合不同年份農(nóng)業(yè)發(fā)展的實(shí)際情況，測算耕地美學(xué)景觀價(jià)值。同時(shí)運(yùn)用區(qū)域經(jīng)濟(jì)發(fā)展指標(biāo)對(duì)不同年度價(jià)值進(jìn)行修正。計(jì)算式如下：

式中CL是縣市的耕地美學(xué)景觀價(jià)值，元；是修正因子，通過縣域人均GDP與全省人均GDP的比值獲??；指的是有關(guān)耕地美學(xué)景觀的價(jià)值當(dāng)量因子；E指耕地食物生產(chǎn)的單位價(jià)值，元/hm2，其計(jì)算公式參照式（7）。

1.2.3 耕地多功能權(quán)衡-協(xié)同測度方法

Spearman秩相關(guān)系數(shù)可定量描述2個(gè)變量之間的關(guān)聯(lián)程度與方向[24]。本研究采用Spearman秩相關(guān)分析探討8種耕地功能兩兩之間的相關(guān)關(guān)系，來分析不同功能間的權(quán)衡與協(xié)同關(guān)系。當(dāng)相關(guān)系數(shù)大于0，且通過顯著性檢驗(yàn)時(shí)，則認(rèn)為2種功能之間存在協(xié)同關(guān)系，反之，則表現(xiàn)為權(quán)衡關(guān)系。

為探討多個(gè)變量間之間的空間相關(guān)性，Anselin等[35]在單變量空間自相關(guān)基礎(chǔ)上提出雙變量空間自相關(guān)指數(shù)，該方法可對(duì)不同要素的空間分布的關(guān)聯(lián)性特征進(jìn)行定量分析，已成為研究不同地理要素間空間關(guān)系的一種有效方法[24]。據(jù)此，本研究采取相關(guān)公式為

1.2.4 驅(qū)動(dòng)因素分析

本研究采用可拓展隨機(jī)性環(huán)境影響評(píng)估模型（STIRPAT）分析耕地多功能價(jià)值變化的社會(huì)經(jīng)濟(jì)驅(qū)動(dòng)因素。STIRPAT模型是一個(gè)多變量非線性模型，可提供一個(gè)人類社會(huì)經(jīng)濟(jì)因素對(duì)環(huán)境影響的因果分析框架[36]。借鑒STIRPAT模型的標(biāo)準(zhǔn)形式，本研究參考以往研究成果[19,29]，構(gòu)建分析社會(huì)經(jīng)濟(jì)因素對(duì)耕地多功能價(jià)值影響的擴(kuò)展模型。具體表達(dá)式為

式中為各縣市地均耕地多功能價(jià)值，考慮到浙江省耕地空間分布差異明顯，利用單位面積數(shù)值可消除耕地多功能價(jià)值的面積依賴性；1為各縣市農(nóng)業(yè)人口數(shù)量，1為各縣市人均GDP，1和2分別為各縣市城鎮(zhèn)人均可支配收入和農(nóng)村人均可支配收入，為各縣市地均農(nóng)業(yè)機(jī)械化水平，為各縣市第二產(chǎn)業(yè)比重，為各縣市城鎮(zhèn)化率，1～7分別為各自變量的彈性系數(shù)。通過對(duì)式（12）的變量采取對(duì)數(shù)形成新的模型，運(yùn)用SPSS.22進(jìn)行逐步回歸分析。模型中社會(huì)經(jīng)濟(jì)指標(biāo)數(shù)據(jù)來源于《浙江省統(tǒng)計(jì)年鑒（2001—2016）》和各地級(jí)、縣級(jí)統(tǒng)計(jì)年鑒。

2 結(jié)果與分析

2.1 浙江省耕地多功能價(jià)值時(shí)空變化

由表2可知，浙江省耕地各類功能價(jià)值結(jié)構(gòu)存在顯著差異，其中耕地的水源涵養(yǎng)、氣體調(diào)節(jié)以及食物生產(chǎn)功能價(jià)值占比較大。2000—2015年，浙江省耕地多功能價(jià)值總體呈下降趨勢，由2000年的848.31億元下降至2015年的701.74億元，在2010年以后下降趨勢放緩。然而，浙江省地均耕地多功能價(jià)值經(jīng)歷了先快速減少后逐步提升的過程，地均耕地多功能價(jià)值由2000年的31 206.44元/hm2增長到2015年的32 060.48元/hm2。在耕地各類功能價(jià)值演變中，氣體調(diào)節(jié)功能價(jià)值減少幅度最大，達(dá)到84.22億元，且主要集中在2000—2010年。食物生產(chǎn)功能價(jià)值增加8.5%，美學(xué)景觀功能價(jià)值增長最快，2000—2015年增加了了237.7%。2000-2015年間浙江省耕地農(nóng)業(yè)污染下降了6.5%，表明該地區(qū)耕地受化肥農(nóng)藥污染情況有所緩解。其他功能值呈現(xiàn)下降趨勢：環(huán)境凈化、水源涵養(yǎng)、生物多樣性、社會(huì)保障分別下降31.47%（3.4億元）、15.08%（40.7億元）、47.94%（6.63億元）、31.49%（48.3億元）。

在時(shí)空分布上，浙江省耕地多功能價(jià)值呈現(xiàn)明顯的空間差異（圖2）。耕地多功能價(jià)值高的縣市主要分布在浙北平原和浙東南沿海平原地區(qū)，這些地區(qū)水網(wǎng)密布、雨熱同期，是傳統(tǒng)的農(nóng)業(yè)主產(chǎn)區(qū)。耕地多功能價(jià)值低的縣市則集中在浙西南山地區(qū)，這些地區(qū)耕地資源少、破碎度高，交通不便，農(nóng)業(yè)發(fā)展水平相對(duì)較低。2000—2015年，浙江省大部分縣市耕地多功能價(jià)值呈下降趨勢，其中浙北平原區(qū)杭州市和浙中盆地區(qū)金華市耕地多功能價(jià)值減少最多（圖2d），浙西南山地區(qū)（麗水、衢州等地）提升明顯。

2.2 浙江省耕地多功能權(quán)衡-協(xié)同關(guān)系

浙江省2000和2015年耕地4類主要功能兩兩之間的相關(guān)關(guān)系如表3所示。整體上，浙江省耕地多功能之間存在顯著的協(xié)同關(guān)系，但各類耕地功能之間的協(xié)同關(guān)系差異較大，其中生產(chǎn)與生態(tài)功能、社會(huì)與文化功能協(xié)同關(guān)系較強(qiáng)，生產(chǎn)與文化功能協(xié)同關(guān)系相對(duì)較弱。隨著城鎮(zhèn)化的不斷推進(jìn)，2000—2015年耕地主要功能的協(xié)同關(guān)系總體減弱。其中社會(huì)與文化功能協(xié)同關(guān)系下降最大，而生產(chǎn)與文化功能、生態(tài)與文化功能協(xié)同關(guān)系有所提升。

表2 2000—2015年浙江省耕地多功能價(jià)值測算

圖2 2000—2015年浙江省各市縣耕地多功能價(jià)值時(shí)空分布

表3 浙江省2000和2015年耕地多功能相關(guān)關(guān)系

為進(jìn)一步探討浙江省耕地不同功能在空間上的權(quán)衡-協(xié)同關(guān)系，基于GeoDa軟件對(duì)4類主要功能進(jìn)行雙變量空間自相關(guān)分析，結(jié)果如圖3、圖4所示。浙江省耕地4類功能之間的全局Moran’s I指數(shù)均為正值，且都通過了1%的顯著性水平檢驗(yàn)，說明耕地4類功能之間均存在空間協(xié)同關(guān)系，與表3結(jié)論一致。在空間分布上，4類耕地功能兩兩之間的權(quán)衡-協(xié)同關(guān)系空間分布格局類似，均呈現(xiàn)出顯著的空間異質(zhì)性。具體而言，各類功能的高值協(xié)同區(qū)主要集中在浙北地區(qū)杭州市區(qū)、紹興部分縣市以及慈溪市等，低值協(xié)同區(qū)主要集中在浙西南地區(qū)麗水市及衢州部分縣市；低-高值權(quán)衡區(qū)和高-低值權(quán)衡區(qū)分別集中在高值協(xié)同區(qū)和低值協(xié)同區(qū)周邊。隨著城鎮(zhèn)化的推進(jìn)，2000—2015年浙江省耕地4類功能之間的空間權(quán)衡-協(xié)同關(guān)系逐漸減弱。其中，生產(chǎn)與生態(tài)功能、生產(chǎn)與文化功能、生態(tài)與文化功能的高值協(xié)同區(qū)和2類權(quán)衡區(qū)域減少。生產(chǎn)與社會(huì)、生態(tài)與社會(huì)、社會(huì)與文化功能的空間協(xié)同關(guān)系減弱，表現(xiàn)為高值協(xié)同區(qū)減少、低-高值權(quán)衡區(qū)增加。

2.3 浙江省耕地多功能價(jià)值變化影響因素

將浙江省2000、2010和2015年各縣市所選取的社會(huì)經(jīng)濟(jì)指標(biāo)以及地均耕地多功能價(jià)值，根據(jù)式（12）取其對(duì)數(shù)，運(yùn)用逐步回歸方法進(jìn)行線性回歸分析。在比較多個(gè)線性回歸模型結(jié)果后，選取最優(yōu)線性回歸方程，表達(dá)式為

ln=9.322-1.81ln1＋1.651ln2-0.135ln

模型R為0.648，且通過了顯著性檢驗(yàn)，說明該模型的擬合度較好，3個(gè)自變量1、2和均在0.01的水平下對(duì)因變量產(chǎn)生顯著影響。這表明城鎮(zhèn)人均可支配收入（1）、農(nóng)村人均可支配收入（2）以及地均農(nóng)業(yè)機(jī)械總動(dòng)力（）是引起浙江省各市縣地均耕地多功能價(jià)值（）空間差異的重要因素。而其他因素如農(nóng)業(yè)人口（1）、人均GDP（1）、第二產(chǎn)業(yè)比重（）以及城鎮(zhèn)化率（）對(duì)因變量影響并不顯著。具體而言，城鎮(zhèn)人均可支配收入和地均農(nóng)業(yè)機(jī)械總動(dòng)力的彈性系數(shù)為負(fù)（-1.810，-0.135），表明城鎮(zhèn)居民可支配收入或地均農(nóng)業(yè)機(jī)械動(dòng)力每增加1%，地均耕地多功能價(jià)值會(huì)相應(yīng)降低1.81%或0.135%。農(nóng)村居民可支配收入的彈性系數(shù)為正（1.651），說明農(nóng)村居民收入的增加可有效提升耕地多功能價(jià)值。

圖3 2000年浙江省耕地4類功能雙變量空間局部自相關(guān)集聚分布

圖4 2015年浙江省耕地4類功能雙變量空間局部自相關(guān)集聚分布

3 討 論

快速城市化不僅降低了耕地多功能價(jià)值，也改變了耕地各功能之間的權(quán)衡-協(xié)同關(guān)系。耕地多功能的時(shí)空變化及其權(quán)衡-協(xié)同關(guān)系研究旨在實(shí)現(xiàn)不同區(qū)域、不同層次的耕地多功能的協(xié)同發(fā)展，在充分認(rèn)識(shí)耕地多功能時(shí)空變化、交互關(guān)系及驅(qū)動(dòng)機(jī)制基礎(chǔ)上，尋求區(qū)域經(jīng)濟(jì)發(fā)展與耕地多功能保護(hù)之間的平衡[15]。價(jià)值量化方法是評(píng)估耕地多功能的重要手段，基于耕地多功能視角全面評(píng)估耕地資源價(jià)值可進(jìn)一步豐富耕地保護(hù)補(bǔ)償研究。該研究結(jié)果顯示，浙江省2000—2015年耕地多功能總價(jià)值呈現(xiàn)下降趨勢，這與湖南、重慶等中西部耕地資源功能價(jià)值呈上升趨勢的結(jié)論有所不同[10,19]，這主要是由于中西部地區(qū)耕地利用強(qiáng)度不斷增強(qiáng)，農(nóng)作物播種面積不降反升。而對(duì)于快速發(fā)展的東部沿海地區(qū)而言，耕地占用、農(nóng)業(yè)非糧化等現(xiàn)象普遍存在[37]，導(dǎo)致耕地多功能價(jià)值下降。然而，浙江省地均耕地多功能價(jià)值呈現(xiàn)出先降低后提升的趨勢，表明城鎮(zhèn)化進(jìn)程中浙江省耕地多功能利用效益在逐步顯現(xiàn)和提升。農(nóng)業(yè)活動(dòng)帶來的耕地利用負(fù)效應(yīng)不容忽視，但由于耕地負(fù)功能定量方法復(fù)雜，本研究采用的農(nóng)業(yè)污染價(jià)值量化方法可能存在重復(fù)估算的不足，但其結(jié)果并不影響研究結(jié)論。未來可進(jìn)一步在耕地正負(fù)功能內(nèi)涵界定、定量方法優(yōu)化等方面開展深入研究。

另外，本研究顯示耕地多功能之間主要表現(xiàn)為協(xié)同關(guān)系，這種協(xié)同關(guān)系在快速城鎮(zhèn)化過程中逐漸減弱，這與朱慶瑩等[6]研究結(jié)論一致。同時(shí)，耕地不同功能的權(quán)衡-協(xié)同關(guān)系空間分異特征明顯。因此，可根據(jù)不同地區(qū)耕地功能權(quán)衡-協(xié)同關(guān)系及其變化制定耕地利用和保護(hù)差異化策略。具體而言，可將耕地功能高-高值協(xié)同區(qū)劃為“耕地優(yōu)勢功能區(qū)”，主要分布在浙北平原縣市；耕地功能低-低值協(xié)同區(qū)劃為“耕地全面優(yōu)化區(qū)”，主要集中在浙西南山地縣市；耕地功能低-高值和高-低值權(quán)衡區(qū)容易受到周邊地區(qū)耕地利用的影響，可劃為“耕地重點(diǎn)優(yōu)化區(qū)”。通過加強(qiáng)縣域間合作，采用工程修復(fù)技術(shù)改良耕地本底質(zhì)量，運(yùn)用經(jīng)濟(jì)、政策等手段優(yōu)化耕地利用模式，促進(jìn)耕地多功能合理利用和保護(hù)。

4 結(jié) 論

本研究以經(jīng)濟(jì)快速發(fā)展地區(qū)浙江省為例，從生產(chǎn)、生態(tài)、社會(huì)和文化方面選取8種典型功能，定量評(píng)估了2000—2015年浙江省耕地多功能價(jià)值，并運(yùn)用相關(guān)系數(shù)法、雙變量空間自相關(guān)模型以及STIRPAT模型揭示了耕地多功能之間的權(quán)衡-協(xié)同關(guān)系的時(shí)空格局演變及影響因素，研究結(jié)果可為經(jīng)濟(jì)快速發(fā)展地區(qū)劃定耕地利用與保護(hù)區(qū)，促進(jìn)耕地多功能協(xié)同利用和推進(jìn)農(nóng)業(yè)可持續(xù)發(fā)展提供科學(xué)參考。主要結(jié)論如下：

1）2000—2015年，浙江省耕地多功能價(jià)值呈下降趨勢，地均耕地多功能價(jià)值表現(xiàn)出先下降后上升的趨勢。其中氣體調(diào)節(jié)、水源涵養(yǎng)和社會(huì)保障價(jià)值下降最明顯，食物生產(chǎn)和美學(xué)景觀價(jià)值有所上升。在空間分布上，耕地多功能價(jià)值較高地區(qū)主要分布在浙江北部平原和東南沿海地區(qū)，低值區(qū)主要集中在西南山地丘陵區(qū)。2000—2015年，浙江省大部分縣市耕地多功能價(jià)值呈現(xiàn)不同程度的下降趨勢，浙西南地區(qū)耕地多功能價(jià)值有所提升。

2）浙江省耕地多功能之間存在顯著協(xié)同關(guān)系，各項(xiàng)功能間的權(quán)衡-協(xié)同關(guān)系存在空間異質(zhì)性。耕地各功能高值協(xié)同區(qū)分布在杭州、紹興等縣市，低值協(xié)同區(qū)集中在浙西南麗水、衢州等地，低-高和高-低值權(quán)衡區(qū)分布在協(xié)同區(qū)周邊。2000—2015年，浙江省耕地多功能協(xié)同關(guān)系逐漸減弱，在空間上主要表現(xiàn)為高值協(xié)同區(qū)減少。

3）地均耕地多功能價(jià)值受到城鄉(xiāng)居民可支配收入和農(nóng)業(yè)機(jī)械化水平的顯著影響。其中城鎮(zhèn)居民可支配收入和地均農(nóng)業(yè)機(jī)械總動(dòng)力與地均耕地多功能價(jià)值呈顯著負(fù)相關(guān)，農(nóng)村居民可支配收入與地均耕地多功能價(jià)值呈顯著正相關(guān)關(guān)系。

[1] 姜廣輝，張鳳榮，孔祥斌，等. 耕地多功能的層次性及其多功能保護(hù)[J]. 中國土地科學(xué)，2011，25(8)：42-47. Jiang Guanghui, Zhang Fengrong, Kong Xiangbin, et al. The difference levels and the protection of multi-functions of cultivated land[J]. China Land Science, 2011, 25(8): 42-47. (in Chinese with English abstract)

[2] Marsden T, Sonnino R. Rural development and the regional state: Denying multifunctional agriculture in the UK[J]. Journal of Rural Studies, 2018, 24(4): 422-431.

[3] 宋小青，歐陽竹. 耕地多功能內(nèi)涵及其對(duì)耕地保護(hù)的啟示[J]. 地理科學(xué)進(jìn)展，2012，31(7)：859-868. Song Xiaoqing, Ouyang Zhu. Connotation of multifunctional cultivated land and its implications for cultivated land protection[J]. Progress in Geography, 2012, 31(7): 859-868. (in Chinese with English abstract)

[4] OECD. Multifunctionality in agriculture: Evaluating the degree of jointness, policy implications[R]. Paris: Organization for Economic Cooperation and Development (OECD), 2008.

[5] 楊雪，談明洪. 近年來北京市耕地多功能演變及其關(guān)聯(lián)性[J]. 自然資源學(xué)報(bào)，2014，29(5)：733-743. Yang Xue, Tan Minghong. Changes and relationships of arable land functions in Beijing in recent years[J]. Journal of Natural Resources, 2014, 29(5): 733-743. (in Chinese with English abstract)

[6] 朱慶瑩，胡偉艷，趙志尚. 耕地多功能權(quán)衡與協(xié)同時(shí)空格局的動(dòng)態(tài)分析：以湖北省為例[J]. 經(jīng)濟(jì)地理，2018，38(7)：143-153. Zhu Qingying, Hu Weiyan, Zhao Zhishang. Dynamic analysis on spatial-temporal pattern of trade-offs and synergies of multifunctional cultivated land: Evidence from Hubei Province[J]. Economic Geography, 2018, 38(7): 143-153. (in Chinese with English abstract)

[7] 彭建，呂慧玲，劉焱序，等. 國內(nèi)外多功能景觀研究進(jìn)展與展望[J]. 地球科學(xué)進(jìn)展，2015，30(4)：465-476. Peng Jian, Lv Huiling, Liu Yanxu, et al. International research progress and perspectives on multifunctional landscape[J]. Advances in Earth Science, 2015, 30(4): 465-476. (in Chinese with English abstract)

[8] De groot R S, Wilson M A, Boumans R M. A typology for the classification, description and valuation of ecosystem functions, goods and services[J]. Ecological Economics, 2002, 41(3): 393-408.

[9] 陳睿山，蔡運(yùn)龍，嚴(yán)祥，等. 土地系統(tǒng)功能及其可持續(xù)性評(píng)價(jià)[J]. 中國土地科學(xué)，2011，25(1)：8-15. Chen Ruishan, Cai Yunlong, Yan Xiang, et al. The functions of land system and its sustainability assessment[J]. China Land Science, 2011, 25(1): 8-15. (in Chinese with English abstract)

[10] 譚雪蘭，安悅，蘇洋，等. 長株潭地區(qū)農(nóng)業(yè)功能的時(shí)空變化特征及發(fā)展策略研究[J]. 地理科學(xué)，2018，38(5)：708-716. Tan Xuelan, An Yue, Su Yang et al. Spatio-temporal dynamics of agricultural function in Changsha-Zhuzhou- Xiangtan urban agglomerations[J]. Scientia Geographica Sinica, 2018, 38(5): 708-716. (in Chinese with English abstract)

[11] Aubry C, Ramamonjisoa J, Dabat M, et al. Urban agriculture and land use in cities?: An approach with the multi-functionality and sustainability concepts in the case of Antananarivo(Madagascar)[J]. Land Use Policy. 2012, 29: 429-439.

[12] Van der Ploeg J D, Laurent C, Blondeau F, et al. Farm diversity, classification schemes and multifunctionality[J]. Journal of Environmental Management. 2009, 90: 124-131.

[13] 范業(yè)婷，金曉斌，項(xiàng)曉敏，等. 蘇南地區(qū)耕地多功能評(píng)價(jià)與空間特征分析[J]. 資源科學(xué)，2018，40(5)：980-992. Fan Yeting, Jin Xiaobin, Xiang Xiaomin, et al. Evaluation and spatial characteristics of arable land multifunction in southern Jiangsu[J]. Resources Science, 2018, 40(5): 980-992. (in Chinese with English abstract)

[14] 辛蕓娜，孔祥斌，鄖文聚. 北京大都市邊緣區(qū)耕地多功能評(píng)價(jià)指標(biāo)體系構(gòu)建：以大興區(qū)為例[J]. 中國土地科學(xué)，2017，31(8)：77-87. Xin Yunna, Kong Xiangbin, Yun Wenju. Design and application of multi-functional evaluation index system for cultivated land in metropolitan fringe of Beijing: A case study in Daxing District[J]. China Land Science, 2017, 31(8): 77-87. (in Chinese with English abstract)

[15] 宋小青，歐陽竹. 中國耕地多功能管理的實(shí)踐路徑探討[J]. 自然資源學(xué)報(bào)，2012，27(4)：540-551. Song Xiaoqing, Ouyang Zhu. Route of multifunctional cultivated land management in China[J]. Journal of Natural Resources, 2012, 27(4): 540-551. (in Chinese with English abstract)

[16] 何青澤，謝德體，王三，等. 重慶市北碚區(qū)耕地多功能評(píng)價(jià)及時(shí)空演變分析[J]. 西南大學(xué)學(xué)報(bào)，2020，42(1)：79-87. He Qingze, Xie Deti, Wang San, et al. Spatio-temporal evolution of cultivated land multi-function based on a fuzzy optimization model in Beibei District, Chongqing[J]. Journal of Southwest University: Natural Science Edition, 2020, 42(1): 79-87. (in Chinese with English abstract)

[17] 靳文娟，邊振興，魏忠義，等. 沈陽城市周邊耕地多功能評(píng)價(jià)[J]. 中國農(nóng)業(yè)資源與區(qū)劃，2019，40(5)：166-173. Jin Wenjuan, Bian Zhenxing, Wei Zhongyi, et al. Multi-functional evaluation of cultivated land in peri-urban areas of Shengyang City[J]. Chinese Journal of Agricultural Resources and Regional Planning, 2019, 40(5): 166-173. (in Chinese with English abstract)

[18] 馬鳳嬌，劉金銅. 基于能值分析的農(nóng)田生態(tài)系統(tǒng)服務(wù)評(píng)估：以河北省欒城縣為例[J]. 資源科學(xué)，2014，36(9)：1949-1957. Ma Fengjiao, Liu Jintong. Agricultural ecosystem services assessment based on energy analysis in Luancheng County[J]. Resources Science, 2014, 36(9): 1949-1957. (in Chinese with English abstract)

[19] 孔凡靖，陳玉成，陳慶華，等. 重慶市農(nóng)田生態(tài)服務(wù)價(jià)值時(shí)空變化特征及其驅(qū)動(dòng)因素分析[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào)，2019，27(11)：1-12. Kong Fanjing, Chen Yucheng, Chen Qinghua, et al. Temporal and spatial variation and driving factors of farmland ecological service value in Chongqing[J]. Chinese Journal of Eco-Agriculture, 2019, 27(11): 1-12. (in Chinese with English abstract)

[20] 李廣東，邱道持，王平. 三峽生態(tài)脆弱區(qū)耕地非市場價(jià)值評(píng)估[J]. 地理學(xué)報(bào)，2011，66(4)：562-575. Li Guangdong, Qiu Daochi, Wang Ping. Assessing non-market value of cultivated land in ecologically fragile areas of Three Gorges Reservoir[J]. Acta Geographica Sinica, 2011, 66(4): 562-575. (in Chinese with English abstract)

[21] 元媛，劉金銅，靳占忠. 欒城縣農(nóng)田生態(tài)系統(tǒng)服務(wù)功能正負(fù)效應(yīng)綜合評(píng)價(jià)[J]. 生態(tài)學(xué)雜志，2011，30(12)：2809-2814. Yuan Yuan, Liu Jintong, Jin Zhanzhong. An integrated assessment of positive and negative effects of high-yielding cropland ecosystem services in Luancheng County, Hebei Province of North China[J]. Chinese Journal of Ecology, 2011, 30(12): 2809-2814. (in Chinese with English abstract)

[22] Swallow B M, Shang J K, Nyabenge M, et al. Tradeoffs, synergies and traps among ecosystem services in the Lake Victoria basin of East Africa[J]. Environmental Science & Policy, 2009, 12(4): 504-519.

[23] 郝夢雅，任志遠(yuǎn)，孫藝杰，等. 關(guān)中盆地生態(tài)系統(tǒng)服務(wù)的權(quán)衡與協(xié)同關(guān)系動(dòng)態(tài)分析[J]. 地理研究，2017，36(3)：592-602. Hao Mengya, Ren Zhiyuan, Sun Yijie, et al. The dynamic analysis of trade-off and synergy of ecosystem services in the Guanzhong Basin[J]. Geographical Research, 2017, 36(3): 592-602. (in Chinese with English abstract)

[24] 鄭德鳳，郝帥，呂樂婷，等. 三江源國家公園生態(tài)系統(tǒng)服務(wù)時(shí)空變化及權(quán)衡-協(xié)同關(guān)系[J]. 地理研究，2020，39(1)：64-78. Zheng Defeng, Hao Shuai, Lv Leting, et al. Spatial-temporal change and trade-off/synergy relationships among multiple ecosystem services in Three-River-Source National Park[J]. Geographical Research, 2020, 39(1): 64-78. (in Chinese with English abstract)

[25] 潘競虎，李真. 干旱內(nèi)陸河流域生態(tài)系統(tǒng)服務(wù)空間權(quán)衡與協(xié)同作用分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33(7)：280-289. Pan Jinghu, Li Zhen. Analysis on trade-offs and synergies of ecosystem services in arid inland river basin[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(7): 280-289. (in Chinese with English abstract)

[26] 方瑩，王靜，孔雪松，等. 耕地利用多功能權(quán)衡關(guān)系測度與分區(qū)優(yōu)化[J]. 中國土地科學(xué)，2018，32(11): 57-64.Fang Ying, Wang Jing, Kong Xuesong. Trade-off relation measurement and zoning optimization of multi-functionality of cultivated land use: A case study of Henan Province[J]. China Land Science, 2018, 32(11): 57-64. (in Chinese with English abstract)

[27] 陳占明，吳施美，馬文博，等. 中國地級(jí)以上城市二氧化碳排放的影響因素分析：基于擴(kuò)展的STIRPAT 模型[J]. 中國人口·資源與環(huán)境，2018，28(10)：45-54. Chen Zhanming, Wu Shimei, Ma Wenbo, et al. Driving forces of carbon dioxide emission for China’s cities: Empirical analysis based on extended STIRPAT Model[J]. China Population, Resources and Environment, 2018, 28(10): 45-54. (in Chinese with English abstract)

[28] 李武艷，朱從謀，王華，等. 浙江省耕地質(zhì)量多尺度空間自相關(guān)分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2016，32(23)：239-245. Li Wuyan, Zhu Congmou, Wang Hua, et al. Multi-scale spatial autocorrelation analysis of cultivated land quality in Zhejiang province[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(23): 239-245. (in Chinese with English abstract)

[29] 李景剛，歐名豪，張效軍，等. 耕地資源價(jià)值重建及其貨幣化評(píng)價(jià)：以青島市為例[J]. 自然資源學(xué)報(bào)，2009，24(11)：1870-1880. Li Jinggang, Ou Minghao, Zhang Xiaojun, et al. Reconstruction of cultivated land resources value system and its evaluation: A case study of Qingdao[J]. Journal of Natural Resources, 2009, 24(11): 1870-1880. (in Chinese with English abstract)

[30] 岳東霞，杜軍，鞏杰，等. 民勤綠洲農(nóng)田生態(tài)系統(tǒng)服務(wù)價(jià)值變化及其影響因子的回歸分析[J]. 生態(tài)學(xué)報(bào)，2011，31(9)：2567-2575. Yue Dongxia, Du Jun, Gong Jie, et al. Dynamic analysis of farmland ecosystem service value and multiple regression analysis of the influence factors in Minqin Oasis[J]. Acta Ecologica Sinica, 2011, 31(9): 2567-2575. (in Chinese with English abstract)

[31] Yu M, Yang Y J, Chen F, et al. Response of agricultural multifunctionality to farmland loss under rapidly urbanizing processes in Yangtze River Delta, China[J]. Science of the Total Environment, 2019, 666: 1–11.

[32] 謝高地，張彩霞，張雷明. 基于單位面積價(jià)值當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價(jià)值化方法改進(jìn)[J]. 自然資源學(xué)報(bào)，2015，30(8)：1243-1254. Xie Gaodi, Zhang Caixia, Zhang Leiming. Improvement of the evaluation method for ecosystem service value based on per unit area[J]. Journal of Natural Resources, 2015, 30(8): 1243-1254. (in Chinese with English abstract)

[33] 鄒月，周忠學(xué). 西安市景觀格局演變對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2017，28(8)：2629-2639. Zou Yue, Zhou Zhongxue. Impact of landscape pattern change on ecosystem service value of Xi’an City, China[J]. Chinese Journal of Applied Ecology, 2017, 28(8): 2629-2639. (in Chinese with English abstract)

[34] 鄧明君，鄧俊杰，劉佳宇. 中國糧食作物化肥施用的碳排放時(shí)空演變與減排潛力[J]. 資源科學(xué)，2016，38(3)：534-544. Deng Mingjun, Deng Junjie, Liu Jiayu. On the space-time evolution of carbon emissions and reduction potential in Chinese grain crop fertilizer application[J]. Resources Science, 2016, 38(3): 534-544. (in Chinese with English abstract)

[35] Anselin L, Syabri I, Kho Y. GeoDa: An introduction to spatial data analysis[J]. Handbook of Applied Spatial Analysis, 2010, 38(1): 73-89.

[36] 唐秀美，郝星耀，劉玉，等. 生態(tài)系統(tǒng)服務(wù)價(jià)值驅(qū)動(dòng)因素與空間異質(zhì)性分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2016，47(5)：336-342. Tang Xiumei, Hao Xingyao, Liu Yu, et al. Driving factors and spatial heterogeneity analysis of ecosystem services value[J].Transactions of the Chinese Society for Agricultural Machinery, 2016，47(5): 336-342.

[37] Su Y, Li C, Wang K, et al. Quantifying the spatiotemporal dynamics and multi-aspect performance of non-grain production during 2000-2015 at a fine scale[J]. Ecological Indicators, 2019, 101: 410-419.

Spatial-temporal change, trade-off and synergy relationships of cropland multifunctional value in Zhejiang Province, China

Zhu Congmou1，Li Wuyan2，Du Yingying3，Xu Hongwei1，Wang Ke1※

(1.,,310058,; 2.,,310018,; 3.,,310018,)

Multifunctional Agriculture (MFA) management has been increasingly important for sustainable cropland utilization in recent years. The farmland loss and changing agricultural activities have posed serious impacts on the structure and functions of cropland in the rapid urbanization. Therefore, it is necessary to explore the spatial-temporal evolution characteristics, the trade-off and synergy relationships of multifunctional cropland. In this study, taking Zhejiang Province, a representative region with rapid social-economic development in China as an example, a series of widely used quantitative methods were employed to systematically evaluate cropland multifunctional value on the county scale. According to four dimensions of cropland production, ecological, social, and cultural functions, eight main functions of cropland were identified, including food production, atmosphere regulation, environmental purification, water conservation, biodiversity, agricultural pollution, social security, and cultural leisure. Spearman's rank correlation coefficient method, bivariate local spatial autocorrelation model and Stochastic Impacts by Regression on Population Affluence and Technology (STIRPAT) model were used to investigate the trade-off and synergy relationships among cropland multifunctions, and the social-economic driving factors. The data was collected from Zhejiang Statistical Yearbook (2001-2016); digital land use/cover maps for 2000, 2010 and 2015 were obtained from the Chinese Ministry of Environmental Protection; the precipitation data came from the National Meteorological Information Center; and the soil type data of Zhejiang Province was from the second soil survey database of China. The results showed that: 1)the total value of cropland multifunctionality displayed a declined trend during 2000 and 2015, while the value of cropland multifunctionality per hectare decreased first and then gradually increased. The value of atmospheric regulation, water conservation and social security decreased significantly. Meanwhile, the value of food production remained constant, whereas, the value of cultural leisure increased rapidly. The total value of cropland multifunctionality showed different spatial characteristics, indicating high values in the north, but low in the south of Zhejiang Province. In the counties, Hangzhou City showed the highest value, while the lowest value was in Shengsi County. The cropland multifunctional value decreased in many counties to some degree, but the value increased dramatically in the southwest mountainous areas. 2) Synergy relationships were dominant among cropland multifunction, showing spatial heterogeneity. The high-high synergy areas were distributed in Hangzhou and Shaoxing cities, while the low-low synergy areas were clustered in Lishui and Quzhou cities in southwestern Zhejiang Province, where the low-high and high-low trade-off areas were distributed around the synergy areas. During 2000 and 2015, synergy relationships among cropland multifunctions totally weakened, as the decrease in the number of high-high synergy areas. 3) Social-economic system was proved to be related closely with cropland multifunctional value. Specifically, the disposable income of urban residents and the total power of agricultural machinery have negative impacts on the multifunctional value of cropland per hectare. The increase in the disposable income of rural residents can be conducive to promote the multifunctional value of cropland per hectare. It infers that the food production of cropland can be balanced with other functions, including ecological, social, and cultural functions, in order to develop multifunctional agriculture, and further to maximize the human welfare that provided by cropland production system. The findings can provide an important decision-making support to delimit cropland use and protection zone, and further to promote the overall coordinated utilization of multifunctional cropland.

land use; models; cropland; trade-off; synergy; multi-function; Zhejiang Province

朱從謀，李武艷，杜瑩瑩，等. 浙江省耕地多功能價(jià)值時(shí)空變化與權(quán)衡-協(xié)同關(guān)系[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2020，36(14)：263-272.doi：10.11975/j.issn.1002-6819.2020.14.032 http://www.tcsae.org

Zhu Congmou，Li Wuyan，Du Yingying, et al. Spatial-temporal change, trade-off and synergy relationships of cropland multifunctional value in Zhejiang Province, China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(14): 263-272. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.14.032 http://www.tcsae.org

2020-01-16

2020-07-01

國家自然科學(xué)基金面上項(xiàng)目（41971236）；浙江省基礎(chǔ)公益研究計(jì)劃（LGN18D010002）

朱從謀，博士生，主要研究方向?yàn)橥恋乩门c規(guī)劃管理Email：congmouzhu1993@163.com

王珂，博士，教授，博士生導(dǎo)師，主要從事資源環(huán)境遙感監(jiān)測、評(píng)價(jià)、規(guī)劃及其信息系統(tǒng)研制。Email：kwang@zju.edu.cn

10.11975/j.issn.1002-6819.2020.14.032

S-01

A

1002-6819(2020)-14-0263-10