王芳，張路路

（1.滄州市地產(chǎn)開發(fā)服務(wù)中心，河北 滄州 061000；2.河北經(jīng)貿(mào)大學(xué)旅游學(xué)院，河北 石家莊 050061）

耕地是糧食安全和農(nóng)業(yè)可持續(xù)發(fā)展的保障，也是體現(xiàn)鄉(xiāng)村景觀風(fēng)貌和傳承農(nóng)耕文化的重要載體[1]。耕地作為農(nóng)業(yè)經(jīng)營最主要的元素，已演變成為一個自然-社會-經(jīng)濟-生態(tài)復(fù)合系統(tǒng)，兼具農(nóng)業(yè)生產(chǎn)、生態(tài)服務(wù)和景觀文化等多種屬性[2]。受經(jīng)濟發(fā)展、社會生活方式改變和居民消費結(jié)構(gòu)升級驅(qū)動等影響，耕地呈現(xiàn)出功能多元化和復(fù)合化的趨勢[3]。長期以來，人們側(cè)重于對耕地單一生產(chǎn)功能的認識和利用，忽略了其他功能，甚至產(chǎn)生了割裂或者對立現(xiàn)象[4]。在鄉(xiāng)村一二三產(chǎn)業(yè)融合背景下，亟待充分顯化并協(xié)調(diào)耕地的多種功能以確保其可持續(xù)利用。

自20 世紀70 年代以來，在全球環(huán)境問題和傳統(tǒng)農(nóng)業(yè)升級改造的驅(qū)動下，耕地多功能成為農(nóng)業(yè)和地理學(xué)等領(lǐng)域關(guān)注的重要熱點之一[5，6]。許多學(xué)者就耕地多功能的內(nèi)涵[7，8]、多功能構(gòu)成[9，10]、指標(biāo)體系的構(gòu)建[11]、多功能評價[12，13]、時空演變特征[7，14]以及功能影響因素[1]等領(lǐng)域進行了研究，有效地促進了耕地的保護與管理。整體來看，耕地功能體現(xiàn)了耕地系統(tǒng)提供滿足人類生存與發(fā)展各種需求的程度和能力[15]，具體可分為生產(chǎn)、生活和生態(tài)功能[16]，具有顯著的層次性和階段性特征，與社會經(jīng)濟發(fā)展和利益相關(guān)者聯(lián)系緊密[4]。此外，耕地各功能間還體現(xiàn)出主從性和一定的排他性[4]，多功能間此消彼長的拮抗和相互增益的協(xié)同關(guān)系日益復(fù)雜[17，18]。部分學(xué)者在空間上探討了耕地多功能權(quán)衡與協(xié)同[6]，但缺乏對地區(qū)耕地多功能之間長時間序列的耦合互動關(guān)系及其演化規(guī)律的分析，對于多功能之間相互促進和脅迫的關(guān)系研究尚不多見[2]。耕地各子系統(tǒng)之間的協(xié)調(diào)性直接關(guān)系著耕地系統(tǒng)的健康安全和可持續(xù)發(fā)展[2]，而動態(tài)耦合模型是測度系統(tǒng)之間和內(nèi)部子系統(tǒng)相互作用程度及協(xié)調(diào)配合情況的重要模型，已廣泛應(yīng)用于生態(tài)與城鎮(zhèn)化等領(lǐng)域[19]，并取得了豐碩成果，為耕地多功能協(xié)同關(guān)系研究提供了借鑒。

滄州市作為河北沿海地區(qū)率先發(fā)展增長極和京津冀區(qū)域新的增長極，目前正處于快速城鎮(zhèn)化、工業(yè)化發(fā)展階段。在國家大力打造環(huán)首都現(xiàn)代農(nóng)業(yè)科技示范帶和農(nóng)村一二三產(chǎn)業(yè)融合政策的導(dǎo)向下，耕地生產(chǎn)、生活與生態(tài)功能亟待融合和轉(zhuǎn)型。鑒于此，本研究以滄州市耕地資源為例，借助2005～2015 年的面板數(shù)據(jù)，通過構(gòu)建指標(biāo)體系與熵值法定量評價耕地多功能，運用動態(tài)耦合模型揭示滄州市耕地生產(chǎn)、生活和生態(tài)功能間的協(xié)調(diào)演變特征，為提升耕地資源整體利用效應(yīng)和在鄉(xiāng)村發(fā)展中的競爭力提供參考。

1 研究區(qū)概況

滄州市地處環(huán)渤海灣經(jīng)濟區(qū)，位于華北平原東部黑龍港流域。屬暖溫帶大陸性季風(fēng)氣候，雨量集中但分布不均，干濕季明顯，四季分明。滄州市是河北省糧、棉、油集中產(chǎn)區(qū)，耕地面積78.69 萬hm2，優(yōu)質(zhì)耕地少，后備耕地資源匱乏[20]，土壤鹽堿化問題突出。隨著耕地投入水平增加，耕地集約利用水平呈波動上升趨勢[21]，但耕地生態(tài)安全水平亦不容樂觀[22]。隨著農(nóng)業(yè)產(chǎn)業(yè)化經(jīng)營的不斷擴大，經(jīng)營主體由家庭承包逐漸向農(nóng)民合作社、家庭農(nóng)場、農(nóng)業(yè)龍頭企業(yè)等具有法人資格的新型經(jīng)營主體轉(zhuǎn)變，加速了耕地資源的多樣化利用與城鄉(xiāng)轉(zhuǎn)型發(fā)展。近年來，隨著國家沿海發(fā)展戰(zhàn)略和京津綠色農(nóng)副產(chǎn)品保障基地定位的提出，滄州市農(nóng)業(yè)發(fā)展面臨著向都市農(nóng)業(yè)轉(zhuǎn)型升級的機遇，亟需厘清耕地多功能特征和相互關(guān)系可為耕地多功能轉(zhuǎn)型升級提供支撐。

2 數(shù)據(jù)來源與研究方法

2.1 數(shù)據(jù)來源

本研究涉及耕地數(shù)據(jù)和社會經(jīng)濟數(shù)據(jù)，耕地數(shù)據(jù)主要來源于1996～2016 年《河北省土地調(diào)查統(tǒng)計年鑒》，社會經(jīng)濟數(shù)據(jù)來源于1996～2016 年《河北農(nóng)村統(tǒng)計年鑒》及滄州市國民經(jīng)濟和社會發(fā)展統(tǒng)計公報。

2.2 研究方法

2.2.1 耕地多功能表征指標(biāo)體系構(gòu)建 將耕地社會-經(jīng)濟-生態(tài)復(fù)合系統(tǒng)分解為有限的、可被量化的功能指標(biāo)，是多功能評價的重要環(huán)節(jié)。借鑒“三生”功能理論和相關(guān)研究[16，23～25]，將耕地多功能劃分為生產(chǎn)、生活和生態(tài)功能，參考已有研究成果和當(dāng)?shù)赝恋亍⑥r(nóng)業(yè)、規(guī)劃等方面專家意見，建立了耕地多功能表征指標(biāo)體系（表1）。

表1 滄州市耕地多功能評價指標(biāo)體系Table 1 Evaluation index system of cultivated land multi-function in Cangzhou

生產(chǎn)功能是耕地生產(chǎn)活動所帶來的農(nóng)產(chǎn)品供給及促進區(qū)域農(nóng)業(yè)經(jīng)濟、國民經(jīng)濟增長的能力[9]。選擇糧食單產(chǎn)、耕地墾殖率、耕地單位年收益和耕地國民經(jīng)濟貢獻度作為表征指標(biāo)。其中，糧食單產(chǎn)直接反映耕地系統(tǒng)生產(chǎn)力大小；耕地墾殖率反映耕地在區(qū)域土地開發(fā)利用中的比重，體現(xiàn)了耕地系統(tǒng)結(jié)構(gòu)和生產(chǎn)功能的穩(wěn)定性；耕地單位年收益反映出耕地經(jīng)濟產(chǎn)出水平和生產(chǎn)力的高低；耕地國民經(jīng)濟貢獻度反映了通過農(nóng)業(yè)生產(chǎn)和農(nóng)產(chǎn)品產(chǎn)出供給促進農(nóng)業(yè)經(jīng)濟與國民經(jīng)濟增長的能力。

生活功能反映了耕地通過食物生產(chǎn)保障特定區(qū)域內(nèi)人口糧食需求以及承載農(nóng)村剩余勞動力的能力[14]。利用人均糧食供給量、糧食作物面積占比來衡量糧食保障安全；選取反映鄉(xiāng)村農(nóng)業(yè)從業(yè)人員比重的地均勞動力承載指數(shù)、反映農(nóng)業(yè)機械對勞動替代能力的農(nóng)業(yè)機械化水平2 項指標(biāo)來衡量就業(yè)保障情況。

生態(tài)功能指維護農(nóng)村土壤健康和生態(tài)系統(tǒng)平衡的能力[26]。耕地生態(tài)功能不僅與自身多樣性有關(guān)，還受到人類施用化肥等外部行為的制約。選擇農(nóng)田生態(tài)系統(tǒng)多樣性、單位耕地化肥負荷和單位耕地農(nóng)藥負荷指標(biāo)進行衡量。農(nóng)田生態(tài)系統(tǒng)多樣性計算公式[9]為：

式中，D 為農(nóng)田系統(tǒng)多樣性指數(shù)；bi為各類農(nóng)作物播種面積占總播種面積比例；n 為研究區(qū)的農(nóng)作物種類數(shù)，結(jié)合滄州市耕地種植實際[27]，選取小麥、玉米、棉花、油料作物、豆類、谷子6 類典型農(nóng)作物進行計算。

2.2.2 耕地多功能評價模型 采用綜合加權(quán)求和模型對耕地功能進行評價。按照0.20 步長，從小到大將耕地功能分為非常低、較低、中等、較高、非常高5 個級別。

式中，S 為耕地功能值；xij為第i 年j 指標(biāo)經(jīng)標(biāo)準(zhǔn)化處理后的值；w 為指標(biāo)權(quán)重，考慮到耕地功能變化的復(fù)雜性和系統(tǒng)性，采用熵值法確定權(quán)重[9]。

為揭示耕地功能類型的多樣性與數(shù)量的均勻度，采用辛普森倒數(shù)指數(shù)法來計算。具體計算公式為：

式中，Ti為研究區(qū)第i 年的辛普森倒數(shù)指數(shù)，其值越大，表明耕地功能類型越多樣，數(shù)量越均勻；Sij為研究區(qū)第i 年第j 項耕地功能值；n 為研究區(qū)耕地功能的種類數(shù)。

2.2.3 耕地多功能協(xié)調(diào)模型 耕地資源系統(tǒng)處于連續(xù)動態(tài)變化之中，與自然、社會和經(jīng)濟要素間經(jīng)過不停地物質(zhì)、能量、信息的交流與交換以實現(xiàn)系統(tǒng)演變[13]，因此，子系統(tǒng)之間存在著相互促進或脅迫的耦合互動關(guān)系。系統(tǒng)論認為有序度是反映系統(tǒng)協(xié)調(diào)性水平的重要參數(shù)[15]，借鑒動態(tài)耦合度函數(shù)模型[14]，通過測算系統(tǒng)功能水平來分析耕地系統(tǒng)的協(xié)調(diào)性。

假設(shè)耕地生產(chǎn)、生活和生態(tài)功能分別是S1、S2、S3，都是關(guān)于時間t 的連續(xù)函數(shù)Si=f（t），i=1，2，3?？赏ㄟ^對Si=f（t）求導(dǎo)，得到子系統(tǒng)的演化速度VSi=dSi/dt。耕地系統(tǒng)發(fā)展演化遵循“S”型發(fā)展變化趨勢[16]，子系統(tǒng)之間相互作用，多功能之間的脅迫關(guān)系也存在一定的周期性變化規(guī)律。以生產(chǎn)、生態(tài)功能為例，其演化速度分別為VS1、VS3，兩者構(gòu)成了子系統(tǒng)V=f（VS1，VS3）。由于V 的變化是由VS1、VS3引起的，通過V 的變化情況研究耕地功能之間的協(xié)調(diào)性作用關(guān)系。將2 個子系統(tǒng)的協(xié)調(diào)度表示成VS1和VS3的夾角，其公式為：

系統(tǒng)的協(xié)調(diào)性水平判定主要依據(jù)子系統(tǒng)的交互作用的強弱程度。按照夾角的值來進行判斷，在一個演化周期內(nèi)，整個系統(tǒng)將經(jīng)歷低級協(xié)調(diào)共生（-90°～0°）、協(xié)調(diào)發(fā)展（0°～90°）、極限發(fā)展（90°～180°）和螺旋式上升（-180°～-90°）4 個階段[19]。但系統(tǒng)并不是完全按上述4 個階段演化，特別是子系統(tǒng)之間出現(xiàn)拮抗，脅迫作用越來越明顯，系統(tǒng)的惡化程度超過其承載力閾值時，若不采取保護措施緩和矛盾，系統(tǒng)發(fā)展將會出現(xiàn)停滯甚至倒退。4 個發(fā)展階段總是相對存在的，某一條件的改變會導(dǎo)致新一輪協(xié)調(diào)關(guān)系的發(fā)展演化。假設(shè)VSi表示耕地3 個功能變化速率，根據(jù)相關(guān)研究[18]，可以將復(fù)合系統(tǒng)協(xié)調(diào)性狀態(tài)分為8 種類型（表2）。

3 結(jié)果與分析

3.1 滄州市耕地資源多功能評價結(jié)果分析

通過模型計算獲得1995～2015 年滄州市耕地功能值（圖1）。耕地綜合功能呈現(xiàn)上升趨勢，但生產(chǎn)、生活和生態(tài)功能的演變趨勢有所不同。

3.1.1 耕地生產(chǎn)功能 耕地生產(chǎn)功能保持著平穩(wěn)快速的增長趨勢。生產(chǎn)功能由較低狀態(tài)向較高的狀態(tài)轉(zhuǎn)變，1995 年功能值僅0.18；2015 年功能值達0.73，增幅305.56%。1995～2004 年耕地生產(chǎn)功能值的增幅低于2005～2015 年，這與國家農(nóng)業(yè)宏觀政策緊密相關(guān)。1995～2004 年受“三農(nóng)”問題影響和制約，農(nóng)產(chǎn)品生產(chǎn)功能的增速較緩[7]。該階段研究區(qū)內(nèi)糧食單產(chǎn)和耕地墾殖率均呈現(xiàn)波動性變化，糧食單產(chǎn)由3.03 t/hm2增至3.42 t/hm2，平均值為3.82 t/hm2；耕地墾殖率由55.99%調(diào)整至52.91%，平均值為55.12%；耕地單位年收益由0.75 萬元/hm2持續(xù)提升至1.34 萬元/hm2。2004 年以后，耕地生產(chǎn)功能呈現(xiàn)恢復(fù)性增長趨勢，糧食單產(chǎn)持續(xù)提升，由2005年的3.52 t/hm2提高至2015 年的6.05 t/hm2，耕地單位年收益則由1.46 萬元/hm2提升至4.8 萬元/hm2。1995～2015 年耕地的國民經(jīng)濟貢獻度由34.75%下降至10.48%，表明近年來耕地的經(jīng)濟創(chuàng)收功能明顯弱化。

表2 滄州市耕地多功能協(xié)調(diào)性狀態(tài)類型評判標(biāo)準(zhǔn)Table 2 Evaluation criteria of cultivated land multi-functional coordination state types in Cangzhou

3.1.2 耕地生活功能 耕地生活功能呈先波動性下降后上升的趨勢。1995～2001年耕地生活功能呈現(xiàn)波動下降趨勢。其中，糧食保障有所提高，人均糧食供給量由372.81 kg/人提升至448.56 kg/人，糧食作物占農(nóng)作物面積比例保持在80%左右；耕地規(guī)模減少和農(nóng)村勞動力的轉(zhuǎn)移，共同導(dǎo)致地均勞動力承載指數(shù)由8.05 人/hm2提高至8.53人/hm2；農(nóng)業(yè)機械化水平呈小幅增長，由4.50%提升至5.50%。2001～2015 年耕地生活保障功能由0.35 提高至0.7，提高幅度為100%。得益于耕地單產(chǎn)水平的提高和城鎮(zhèn)化進程加速，人均糧食供給量顯著提高。與此同時，種植結(jié)構(gòu)小幅調(diào)整，糧食作物占農(nóng)作物的面積由80%下降至78%，農(nóng)業(yè)機械化水平顯著提高，就業(yè)保障功能也逐漸提升，整體上使生活功能由較低狀態(tài)躍遷至較高狀態(tài)。

3.1.3 耕地生態(tài)功能 耕地生態(tài)功能值呈先持續(xù)下降后波動性小幅度提升的趨勢，主要分為2 個階段。第1 階段為1995～2007年，表現(xiàn)為耕地生態(tài)功能持續(xù)下降。糧食生產(chǎn)功能的提升主要得益于農(nóng)業(yè)技術(shù)和農(nóng)資投入的增長，由此也帶來了農(nóng)藥、化肥使用過量，土壤污染、板結(jié)等問題，生態(tài)功能級別由非常高下降至較低狀態(tài)。單位耕地農(nóng)藥負荷由13.77 kg/hm2增加至18.34 kg/hm2，單位耕地化肥負荷由281.05 kg/hm2增加至411.31 kg/hm2。單一化的種植結(jié)構(gòu)使得農(nóng)田生態(tài)系統(tǒng)多樣性指數(shù)由1.65 下降至1.26。第2 階段為2008～2015 年，表現(xiàn)為耕地生態(tài)功能波動性變化，并有小幅提升，處于中等級別。隨著國家“生態(tài)文明”戰(zhàn)略的實施，耕地生態(tài)問題逐漸受到重視，耕地保護逐步由“數(shù)量、質(zhì)量”向“數(shù)量、質(zhì)量、生態(tài)”一體化管護方向發(fā)展。2015 年化肥和農(nóng)藥負荷量分別較2008 年降低了12.28%和31.31%；通過土地整治，農(nóng)田生態(tài)系統(tǒng)多樣性指數(shù)降低趨勢減緩，由2008 年的1.25 下降至2015 年的1.19。

圖1 1995～2015 年滄州市耕地多功能演化Fig.1 Multi-functional evolution ofcultivated land from 1995 to 2015 in Cangzhou

3.1.4 耕地綜合功能 耕地綜合功能整體呈現(xiàn)增長趨勢，由1995 年的0.43 提高至2015 年的0.65。其中，1995～2004 年呈現(xiàn)波動性變化；2004 年以后該指標(biāo)值持續(xù)上升。結(jié)合耕地生產(chǎn)、生活和生態(tài)功能的變化來看，滄州市耕地多功能利用轉(zhuǎn)型在2004～2007 年；隨后，耕地的生產(chǎn)功能、生活功能持續(xù)穩(wěn)定上升，生態(tài)功能減弱明顯趨緩，開始沿城鄉(xiāng)居民多需求發(fā)展方向演變。

辛普森倒數(shù)指數(shù)與耕地綜合功能變化相反。1995～2007 年耕地多功能的辛普森倒數(shù)指數(shù)波動頻繁，但整體呈增加狀態(tài)；2007～2015 年辛普森倒數(shù)指數(shù)一直呈下降狀態(tài)。

綜合分析表明，滄州市耕地多功能利用轉(zhuǎn)型發(fā)生在2007 年。這與宋小青等[9]關(guān)于全國耕地多功能利用轉(zhuǎn)型研究的結(jié)論相似，后者認為全國耕地功能自2006 年開始總體增強，各類功能均衡發(fā)展態(tài)勢明顯。

圖2 滄州市耕地綜合功能和辛普森倒數(shù)指數(shù)的變化Fig.2 Changes of comprehensive function and Simpson reciprocal index ofcultivated land in Cangzhou

3.2 滄州市耕地多功能耦合協(xié)調(diào)度

利用Matlab 7.0 軟件，對滄州市1995～2015 年耕地生產(chǎn)、生活和生態(tài)功能的變化進行非線性曲線擬合，在保證最佳擬合優(yōu)度的前提下，得到各功能的曲線方程（表3）。

對Si=f（t）函數(shù)求導(dǎo)，得到耕地生產(chǎn)功能、生活功能和生態(tài)功能的演化速度，以及各個功能之間的協(xié)調(diào)度；最后根據(jù)協(xié)調(diào)性類型判定標(biāo)準(zhǔn)劃分類型（表4）。數(shù)據(jù)顯示，盡管滄州耕地功能整體上在增加，但多功能之間的多樣化協(xié)同效應(yīng)還未得到發(fā)揮。1995～2005 年，耕地功能協(xié)同性表現(xiàn)為生態(tài)功能滯后型，即耕地生產(chǎn)功能有所增長，生活功能有所進步，而生態(tài)功能有所下降；2006～2015 年，耕地功能協(xié)同性則表現(xiàn)為生態(tài)-生產(chǎn)滯后型，即生產(chǎn)功能盡管在增長，生態(tài)功能也有小幅提升，但與生活功能的增幅不匹配。

表3 耕地多功能演化擬合曲線Table 3 Fitting curve of cultivated land multi-function evolution

表4 1995～2015 年滄州市耕地多功能演化速度Table 4 Multi-functional evolution rate of cultivated land in Cangzhou from 1995 to 2015

進一步分析單項功能之間的協(xié)同性，結(jié)果顯示（圖3），1995～2015 年滄州市耕地生產(chǎn)-生活功能經(jīng)歷了協(xié)調(diào)發(fā)展和低級協(xié)同共生的階段。耕地的生產(chǎn)功能和生活功能作為耕地的傳統(tǒng)與基本功能，受社會經(jīng)濟發(fā)展水平、政策等多因素的影響，1995～2005 年農(nóng)戶仍以家庭經(jīng)營方式種植傳統(tǒng)糧食作物為主要生計來源，糧食生產(chǎn)穩(wěn)定，農(nóng)村種植業(yè)從業(yè)人員比例和人均糧食保證率等表征的耕地生活功能較強，耕地生產(chǎn)與生活功能的協(xié)同度處于較好狀態(tài)。2006 年以來，隨著耕地投入和集約化水平的增加，耕地產(chǎn)出的邊際效應(yīng)逐漸顯現(xiàn)，生產(chǎn)功能的提升幅度開始下降，與此同時，對耕地的社會保障功能需求日益提升，使生產(chǎn)與生活功能的協(xié)同性呈現(xiàn)迅速下降而后又逐漸有所提升的趨勢。

圖3 1995～2015 年滄州市耕地生產(chǎn)-生活-生態(tài)功能協(xié)調(diào)度的變化Fig.3 Coordination degree change of cultivated land production-life-ecology function in Cangzhou from 1995 to 2015

生產(chǎn)-生態(tài)功能協(xié)調(diào)性結(jié)果顯示，2000 年以前，耕地經(jīng)濟與生態(tài)子系統(tǒng)處于低級協(xié)調(diào)共生階段。為了追求更大的產(chǎn)出效益，耕地集約化水平、利用強度和程度越來越高，以高強度的農(nóng)藥、化肥投入等作為高產(chǎn)條件，使得耕地生態(tài)功能逐漸減弱，生產(chǎn)功能對生態(tài)功能產(chǎn)生了較強的脅迫作用，反過來生態(tài)功能對生產(chǎn)功能也產(chǎn)生了明顯的約束作用。2000 年以后，生產(chǎn)與生態(tài)功能間的拮抗趨勢開始緩解。政府開始提高耕地保護意識，不斷加大對耕地科技和管理的投入力度，不可持續(xù)的農(nóng)業(yè)生產(chǎn)方式逐漸被取代。2006 年以后，耕地生產(chǎn)-生態(tài)功能協(xié)調(diào)性水平提高，二者逐步進入?yún)f(xié)調(diào)發(fā)展?fàn)顟B(tài)。

1995～2015 年滄州市耕地生活-生態(tài)功能的協(xié)調(diào)性一直處于低級協(xié)調(diào)共生階段。耕地生活功能緩慢提升，但對耕地生態(tài)功能的重視度不高，耕地生態(tài)功能的降低，在直接影響農(nóng)業(yè)經(jīng)濟增長的同時，也降低了農(nóng)民的環(huán)境福祉。隨著區(qū)域整體生態(tài)環(huán)境的下降，耕地產(chǎn)生的生態(tài)價值受到重視，滄州市耕地生活-生態(tài)功能的協(xié)調(diào)性雖然有所提升，但仍處于低級協(xié)調(diào)共生階段。

4 結(jié)論與討論

本研究以滄州市耕地資源為例，以1995～2015 年面板數(shù)據(jù)為主要數(shù)據(jù)源，通過構(gòu)建耕地多功能評價指標(biāo)體系與熵值法定量評價了耕地功能，運用動態(tài)耦合模型揭示了滄州市在社會經(jīng)濟發(fā)展過程中耕地多功能的結(jié)構(gòu)變化和重組規(guī)律，體現(xiàn)了不同社會經(jīng)濟發(fā)展階段對耕地的需求，為耕地資源的合理利用指明了方向，且表明了耕地“數(shù)量—質(zhì)量—生態(tài)”三位一體管護政策的重要價值。主要結(jié)論如下：

（1）研究期間，滄州市耕地綜合功能整體上呈現(xiàn)上升趨勢，而生產(chǎn)、生活和生態(tài)功能的演變出現(xiàn)分異。耕地生產(chǎn)功能保持著平穩(wěn)快速的增長趨勢，耕地生活功能呈先波動性下降后上升的趨勢，耕地生態(tài)功能總體呈先持續(xù)下降后波動性小幅提升趨勢。滄州市耕地多功能利用轉(zhuǎn)型發(fā)生在2007 年，開始沿城鄉(xiāng)居民多需求發(fā)展方向演變。

（2）滄州市耕地多功能之間的多樣化協(xié)同效應(yīng)還未得到發(fā)揮，多功能協(xié)調(diào)前期以生態(tài)功能滯后型為主，后逐步過渡至生態(tài)-生產(chǎn)滯后型。從單項功能來看，耕地生產(chǎn)-生活功能經(jīng)歷了協(xié)調(diào)發(fā)展低級協(xié)同共生的階段，生產(chǎn)-生態(tài)功能由低級協(xié)調(diào)共生階段向協(xié)調(diào)發(fā)展?fàn)顟B(tài)演變，生活-生態(tài)功能一直處于低級協(xié)調(diào)共生階段。功能之間相互耦合協(xié)調(diào)性結(jié)果也支持了耕地保護由單純的數(shù)量保護向“數(shù)量—質(zhì)量—生態(tài)”一體化管護發(fā)展的政策。研究期間，盡管耕地的生產(chǎn)功能在逐漸提升，但由于長期忽視其生態(tài)功能，導(dǎo)致了多功能間的不協(xié)調(diào)，也充分體現(xiàn)了實施一體化管護政策的必要性和急迫性。未來在注重耕地數(shù)量保護的同時，應(yīng)該加強耕地的生態(tài)管護，提升耕地系統(tǒng)的協(xié)調(diào)性，確保耕地可持續(xù)利用。

受研究數(shù)據(jù)的制約，本文未對耕地的文化景觀功能進行探討。同時，未來還可借助模型進行耕地多功能的預(yù)測和協(xié)調(diào)性預(yù)警，為區(qū)域耕地的永續(xù)利用提供參考。