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

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

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

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

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

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

0 引 言

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

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

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

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

1.1 研究區(qū)概況

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

圖1 研究區(qū)概況

1.2 研究方法

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

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

1.2.2 耕地多功能價值估算

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

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

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

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

2）氣體調(diào)節(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。該部分運用成本法進行價值測算，計算式為

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

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

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

式中EP是縣市的耕地環(huán)境凈化價值，元；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)功能價值

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

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

5）生物多樣性功能價值

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

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

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

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

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

7）社會保障功能價值

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

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

8）美學景觀功能價值

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

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

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

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

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

1.2.4 驅(qū)動因素分析

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

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

2 結(jié)果與分析

2.1 浙江省耕地多功能價值時空變化

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

在時空分布上，浙江省耕地多功能價值呈現(xiàn)明顯的空間差異（圖2）。耕地多功能價值高的縣市主要分布在浙北平原和浙東南沿海平原地區(qū)，這些地區(qū)水網(wǎng)密布、雨熱同期，是傳統(tǒng)的農(nóng)業(yè)主產(chǎn)區(qū)。耕地多功能價值低的縣市則集中在浙西南山地區(qū)，這些地區(qū)耕地資源少、破碎度高，交通不便，農(nóng)業(yè)發(fā)展水平相對較低。2000—2015年，浙江省大部分縣市耕地多功能價值呈下降趨勢，其中浙北平原區(qū)杭州市和浙中盆地區(qū)金華市耕地多功能價值減少最多（圖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)功能、社會與文化功能協(xié)同關(guān)系較強，生產(chǎn)與文化功能協(xié)同關(guān)系相對較弱。隨著城鎮(zhèn)化的不斷推進，2000—2015年耕地主要功能的協(xié)同關(guān)系總體減弱。其中社會與文化功能協(xié)同關(guān)系下降最大，而生產(chǎn)與文化功能、生態(tài)與文化功能協(xié)同關(guān)系有所提升。

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

圖2 2000—2015年浙江省各市縣耕地多功能價值時空分布

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

為進一步探討浙江省耕地不同功能在空間上的權(quán)衡-協(xié)同關(guān)系，基于GeoDa軟件對4類主要功能進行雙變量空間自相關(guān)分析，結(jié)果如圖3、圖4所示。浙江省耕地4類功能之間的全局Moran’s I指數(shù)均為正值，且都通過了1%的顯著性水平檢驗，說明耕地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)化的推進，2000—2015年浙江省耕地4類功能之間的空間權(quán)衡-協(xié)同關(guān)系逐漸減弱。其中，生產(chǎn)與生態(tài)功能、生產(chǎn)與文化功能、生態(tài)與文化功能的高值協(xié)同區(qū)和2類權(quán)衡區(qū)域減少。生產(chǎn)與社會、生態(tài)與社會、社會與文化功能的空間協(xié)同關(guān)系減弱，表現(xiàn)為高值協(xié)同區(qū)減少、低-高值權(quán)衡區(qū)增加。

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

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

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

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

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

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

3 討 論

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

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

4 結(jié) 論

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

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

2）浙江省耕地多功能之間存在顯著協(xié)同關(guān)系，各項功能間的權(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）地均耕地多功能價值受到城鄉(xiāng)居民可支配收入和農(nóng)業(yè)機械化水平的顯著影響。其中城鎮(zhèn)居民可支配收入和地均農(nóng)業(yè)機械總動力與地均耕地多功能價值呈顯著負相關(guān)，農(nóng)村居民可支配收入與地均耕地多功能價值呈顯著正相關(guān)關(guān)系。

[1] 姜廣輝，張鳳榮，孔祥斌，等. 耕地多功能的層次性及其多功能保護[J]. 中國土地科學，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)涵及其對耕地保護的啟示[J]. 地理科學進展，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]. 自然資源學報，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é)同時空格局的動態(tài)分析：以湖北省為例[J]. 經(jīng)濟地理，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]. 地球科學進展，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] 陳睿山，蔡運龍，嚴祥，等. 土地系統(tǒng)功能及其可持續(xù)性評價[J]. 中國土地科學，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è)功能的時空變化特征及發(fā)展策略研究[J]. 地理科學，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è)婷，金曉斌，項曉敏，等. 蘇南地區(qū)耕地多功能評價與空間特征分析[J]. 資源科學，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ū)耕地多功能評價指標體系構(gòu)建：以大興區(qū)為例[J]. 中國土地科學，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] 宋小青，歐陽竹. 中國耕地多功能管理的實踐路徑探討[J]. 自然資源學報，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ū)耕地多功能評價及時空演變分析[J]. 西南大學學報，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] 靳文娟，邊振興，魏忠義，等. 沈陽城市周邊耕地多功能評價[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ù)評估：以河北省欒城縣為例[J]. 資源科學，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ù)價值時空變化特征及其驅(qū)動因素分析[J]. 中國生態(tài)農(nóng)業(yè)學報，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ū)耕地非市場價值評估[J]. 地理學報，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ù)功能正負效應綜合評價[J]. 生態(tài)學雜志，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] 郝夢雅，任志遠，孫藝杰，等. 關(guān)中盆地生態(tài)系統(tǒng)服務(wù)的權(quán)衡與協(xié)同關(guān)系動態(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ù)時空變化及權(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è)工程學報，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]. 中國土地科學，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] 陳占明，吳施美，馬文博，等. 中國地級以上城市二氧化碳排放的影響因素分析：基于擴展的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è)工程學報，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] 李景剛，歐名豪，張效軍，等. 耕地資源價值重建及其貨幣化評價：以青島市為例[J]. 自然資源學報，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ù)價值變化及其影響因子的回歸分析[J]. 生態(tài)學報，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] 謝高地，張彩霞，張雷明. 基于單位面積價值當量因子的生態(tài)系統(tǒng)服務(wù)價值化方法改進[J]. 自然資源學報，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] 鄒月，周忠學. 西安市景觀格局演變對生態(tài)系統(tǒng)服務(wù)價值的影響[J]. 應用生態(tài)學報，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] 鄧明君，鄧俊杰，劉佳宇. 中國糧食作物化肥施用的碳排放時空演變與減排潛力[J]. 資源科學，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ù)價值驅(qū)動因素與空間異質(zhì)性分析[J]. 農(nóng)業(yè)機械學報，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

朱從謀，李武艷，杜瑩瑩，等. 浙江省耕地多功能價值時空變化與權(quán)衡-協(xié)同關(guān)系[J]. 農(nóng)業(yè)工程學報，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

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

朱從謀，博士生，主要研究方向為土地利用與規(guī)劃管理Email：congmouzhu1993@163.com

王珂，博士，教授，博士生導師，主要從事資源環(huán)境遙感監(jiān)測、評價、規(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