楊人豪，楊慶媛，陳伊多，曾 黎

(西南大學(xué)地理科學(xué)學(xué)院，重慶 400715)

0 引言

耕地是農(nóng)業(yè)生產(chǎn)發(fā)展最重要的資源，如今正面臨數(shù)量減少[1]、質(zhì)量下降[2]等嚴(yán)重問題。城市建設(shè)發(fā)展占用、生態(tài)退耕等原因引起耕地?cái)?shù)量迅速減少[3]，耕地質(zhì)量因耕地污染[4]、生態(tài)破壞[5]等不斷下降，甚至出現(xiàn)無法耕種的現(xiàn)象。研究耕地資源數(shù)量及其壓力變化，探討變化背后的原因等成為地理學(xué)中的研究熱點(diǎn)[6-8]。 國外學(xué)者多關(guān)注耕地環(huán)境污染、工業(yè)侵占耕地等威脅并尋求解決問題的方法[9-11]； 國內(nèi)學(xué)者更多關(guān)注耕地?cái)?shù)量變化及驅(qū)動因子[12-13]，耕地變化與城鎮(zhèn)化[14]、耕地資源與經(jīng)濟(jì)發(fā)展關(guān)系[15]、糧食安全[16]及影響因素[17]等問題。其中在耕地資源與經(jīng)濟(jì)發(fā)展關(guān)系的研究方面，學(xué)者運(yùn)用脈沖響應(yīng)函數(shù)[18]、Granger因果關(guān)系檢驗(yàn)法[19]、經(jīng)濟(jì)與環(huán)境協(xié)調(diào)度模型[20]等方法探討兩者的關(guān)系，取得了較豐富的成果，為后續(xù)研究提供了理論與方法的借鑒。文章根據(jù)統(tǒng)計(jì)數(shù)據(jù)，辨析1978～2014年河北省耕地壓力指數(shù)變化情況，進(jìn)而運(yùn)用空間計(jì)量方法探索耕地壓力與經(jīng)濟(jì)發(fā)展的地理空間關(guān)聯(lián)與耦合規(guī)律，以期在糧食產(chǎn)量、人口數(shù)量雙增長的背景下探究河北省長期以來的耕地壓力狀況，為區(qū)域制定農(nóng)業(yè)政策、轉(zhuǎn)變耕地利用方式來協(xié)調(diào)經(jīng)濟(jì)發(fā)展與耕地壓力之間的矛盾提供依據(jù)。

1 研究區(qū)概況與數(shù)據(jù)來源

1.1 研究區(qū)概況

河北省地處蒙古高原、燕山和太行山山地、華北平原三者交界處，地勢西北高、東南低，包括高原、山地、丘陵、盆地、平原等多種地形，幅員18.88萬km2，河北省屬溫帶大陸性季風(fēng)氣候，冬夏均溫相差達(dá)20℃，年均降水量484.5mm，年無霜期81～204d。2014年全省總?cè)丝跒? 384萬人，全年地區(qū)生產(chǎn)總值2.942 1萬億元，人均生產(chǎn)總值3.998 4萬元，全年糧食播種面積633.2萬hm2，糧食產(chǎn)量3 360.2萬t。

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

研究中涉及的耕地面積、人口數(shù)量、糧食產(chǎn)量、糧食播種面積、農(nóng)作物播種面積等數(shù)據(jù)主要來源于《河北省統(tǒng)計(jì)年鑒》(1978～2015)、《河北農(nóng)村統(tǒng)計(jì)年鑒》(1978～2015)、《河北省國民經(jīng)濟(jì)和社會發(fā)展統(tǒng)計(jì)公報(bào)》及各市統(tǒng)計(jì)年鑒(1978～2015)等。

2 研究方法

2.1 最小人均耕地面積與耕地壓力指數(shù)

最小人均耕地面積是指保障一定區(qū)域內(nèi)食物安全所需保護(hù)的耕地?cái)?shù)量底線，即在一定區(qū)域范圍內(nèi)，一定食物自給水平和耕地綜合生產(chǎn)能力條件下，為了滿足每個人正常生活的食物消費(fèi)所需的耕地面積[21]。

(1)

式(1)中，Smin為最小人均耕地面積(hm2)；β為糧食自給率(%);Gr為人均年糧食需求量(kg/人);p為糧食單產(chǎn)(kg/hm2);q為糧食播種面積占總播種面積的百分比(%);k為復(fù)種指數(shù)。

耕地壓力指數(shù)是指人均所需最小耕地面積與實(shí)際人均耕地面積之比[21]，可通過耕地壓力指數(shù)的大小判斷區(qū)域耕地—糧食—人口系統(tǒng)中耕地資源緊張程度。

(2)

式(2)中，K為耕地壓力指數(shù);Sa為實(shí)際人均耕地面積(hm2/人)。K<1時(shí)，耕地?zé)o明顯壓力；K=1時(shí)，耕地瀕臨壓力變化的邊界；K>1時(shí)，耕地壓力明顯，若不采取糧食增產(chǎn)措施，可能爆發(fā)區(qū)域糧食危機(jī)。

2.2 地理聯(lián)系率

地理事物在地理空間分布上存在聯(lián)系，可用地理聯(lián)系率刻畫兩種地理事物聯(lián)系的程度[22]。借鑒該方法表征耕地壓力指數(shù)與經(jīng)濟(jì)發(fā)展水平的關(guān)聯(lián)程度。

(3)

式(3)中，G為地理聯(lián)系率;Ki為i地區(qū)的耕地壓力指數(shù);K0為全域內(nèi)的耕地壓力指數(shù);Pi為i地區(qū)的經(jīng)濟(jì)發(fā)展水平;P0為全域內(nèi)的經(jīng)濟(jì)發(fā)展水平。當(dāng)Ki、K0的比值與Pi、P0的比值在空間上分布較為一致時(shí)，G值較高，表明兩者的地理聯(lián)系率較高； 當(dāng)Ki、K0的比值與Pi、P0的比值在空間上分布差異較大時(shí)，G值較低，表明兩者的地理關(guān)聯(lián)性弱。

2.3 重心模型

該文借鑒來自于物理學(xué)領(lǐng)域中的重心理論，將其引入地理學(xué)研究中，某行政區(qū)域某種屬性點(diǎn)重心可借助各次級行政區(qū)的某種屬性和地理坐標(biāo)來表征[23]。

(4)

式(4)中,當(dāng)Mi為各市的耕地壓力指數(shù)或經(jīng)濟(jì)發(fā)展水平時(shí)，則(Xt，Yt)為第t年全省耕地壓力指數(shù)或經(jīng)濟(jì)發(fā)展水平的重心。

為反映重心偏移幅度，引入重心空間偏移距離d。

(5)

式(5)中，(Xti，Yti)，(Xtj，Ytj)分別為第t年的第i和第j個重心坐標(biāo);d為兩個重心坐標(biāo)的空間距離。

2.4 地理集中度與耦合指數(shù)

地理集中度是表明某種屬性在地域上集中程度的指標(biāo)，可用其衡量某一要素在空間分布的狀況[24]。

(6)

式(6)中，Ri為i地區(qū)耕地壓力指數(shù)或經(jīng)濟(jì)發(fā)展水平的地理集中指數(shù);Ai為i地區(qū)的耕地壓力指數(shù)或經(jīng)濟(jì)發(fā)展水平；Ti為i地區(qū)的國土面積。

耕地壓力和經(jīng)濟(jì)發(fā)展的空間耦合關(guān)系可用兩者的地理集中度的比值來表征[25]，計(jì)算公式：

Ii=Ri耕地壓力/Ri經(jīng)濟(jì)發(fā)展水平

(7)

式(7)中，Ii為i地區(qū)的耦合指數(shù);Ri耕地壓力為i地區(qū)耕地壓力的地理集中度；Ri經(jīng)濟(jì)發(fā)展水平為i地區(qū)經(jīng)濟(jì)發(fā)展水平的地理集中度。

3 結(jié)果與分析

3.1 河北省耕地壓力指數(shù)變化

耕地占用和人口增長的結(jié)果必然是耕地壓力的增加，耕地壓力指數(shù)能很好地反映區(qū)域耕地資源的壓力水平。人均年糧食需求量(Gr)取值的不同直接影響耕地壓力水平的指征。眾多研究者對該變量取值為300～400kg/人[26-27]，且世界糧農(nóng)組織提出400kg/人即可滿足正常營養(yǎng)水平。當(dāng)前全國人均年糧食消費(fèi)水平400～500kg/人[28]，考慮現(xiàn)代飲食結(jié)構(gòu)的調(diào)整，把近年人們增加的蔬菜、肉類消費(fèi)部分進(jìn)行折算[29]，將1978～1990年、1991～2000年、2001～2014年Gr分別取值為380kg/人、400kg/人、420kg/人，糧食自給率為100%，計(jì)算結(jié)果如圖1。

圖1 河北省人均最小耕地面積、實(shí)際人均耕地面積及耕地壓力指數(shù)變化

3.1.1 耕地壓力指數(shù)的時(shí)序變化

隨著農(nóng)業(yè)技術(shù)進(jìn)步和人們食物結(jié)構(gòu)的調(diào)整，河北省人均最小耕地面積Smin從1978年的0.15hm2/人震蕩下降至2014年的0.08hm2/人(圖1)，實(shí)際人均耕地面積Sa在37年中基本持續(xù)保持緩慢下降趨勢，從1978年的0.13hm2/人下降至2014年的0.09hm2/人，而耕地壓力指數(shù)K波動較大，在1992～2003年期間明顯呈“U”形，之后基本呈下降趨勢，直到2014年底，耕地壓力指數(shù)K為0.91。圖1中K值在1995年以前除1990年為0.99外都保持大于1，在1978～1995年期間我國農(nóng)業(yè)處于發(fā)展階段，糧食單產(chǎn)量還有待提高，人口數(shù)量增長速度較快，糧食需求量大，耕地壓力較大。在1995～2000年K值處于0.89～0.99之間，期間由于國家長久以來給予三農(nóng)問題的高度關(guān)注，耕地資源得到有效保護(hù)，糧食產(chǎn)量不斷提高，加之1982年以來實(shí)施計(jì)劃生育政策使人口數(shù)量增長得到一定控制，耕地壓力有所下降。在2001～2009年耕地壓力又有所回升，期間河北省城市建設(shè)發(fā)展較快，更多的農(nóng)村青年勞動力進(jìn)入城市，農(nóng)村中空心村、耕地撂荒等現(xiàn)象更加嚴(yán)重。2010年以來K值一直保持在1以下，期間國家對進(jìn)一步擴(kuò)大種糧直補(bǔ)、良種補(bǔ)貼和農(nóng)資綜合補(bǔ)貼的規(guī)模，并且出臺高標(biāo)準(zhǔn)基本農(nóng)田建設(shè)規(guī)劃，對于耕地保護(hù)、糧食安全等給予了更高的關(guān)注。

3.1.2 耕地壓力指數(shù)的空間變化

根據(jù)式(1)和(2)為保持時(shí)間間隔的一致性，利用河北省11市數(shù)據(jù)分別計(jì)算1980年、1990年、2000年、2010年的耕地壓力指數(shù)，并使用ArcGIS 9.3把各市4個年份的耕地壓力指數(shù)進(jìn)行可視化表達(dá)(圖2)，大部分地級市的耕地壓力指數(shù)均從較高水平(1.20≤K≤1.60)下降至較低水平(0.00≤K≤1.20)。1980年河北省僅石家莊市、唐山市的耕地壓力指數(shù)小于1.00，其余9市均大于1； 1990年唐山市、廊坊市、石家莊市、衡水市耕地壓力指數(shù)小于1.00，其余7市耕地壓力指數(shù)位于1.00～1.60之間； 2000年已有超過一半地級市的耕地壓力指數(shù)小于1.00，包括邯鄲市、保定市、石家莊市、邢臺市、滄州市、衡水市，河北省土地侵蝕最為嚴(yán)重的張家口市耕地壓力指數(shù)超過2.00，與其長期以來土地干旱、沙化現(xiàn)象有關(guān)，其余5市除秦皇島市外耕地壓力指數(shù)基本位于1.00～1.20之間； 2010年除了唐山市、承德市、秦皇島市、張家口市耕地壓力指數(shù)大于1.00，其余7市耕地壓力指數(shù)均低于1.00。河北省各市耕地壓力的下降與國家對農(nóng)業(yè)的持續(xù)關(guān)注密切相關(guān)。近年來建設(shè)的渤海糧倉科技示范工程、寧晉糧食綜合示范基地等大幅度提升了河北省糧食生產(chǎn)能力，進(jìn)一步緩解了耕地壓力。

圖2 河北省耕地壓力指數(shù)的空間演變

3.2 河北省經(jīng)濟(jì)發(fā)展變化

河北省經(jīng)濟(jì)發(fā)展水平較快， 2014年GDP達(dá)2.942 1萬億元，城鎮(zhèn)居民人均可支配收入2.414 1萬元； 農(nóng)村居民人均可支配收入1.018 6萬元[30]。2014年河北省GDP與1978年相比增長了近160倍，人均GDP增長了近109倍，GDP和人均GDP都實(shí)現(xiàn)了從百元向萬元的巨大飛躍(圖3)。

作為中國重要的棉糧產(chǎn)區(qū)，河北省不僅種植大量小麥、玉米、高粱等糧食作物，也種植許多棉花、油料、麻類、甜菜等經(jīng)濟(jì)作物， 2014年第一產(chǎn)業(yè)占GDP的11.7%； 改革開放以來，河北省大力發(fā)展工業(yè)，形成以新能源、汽車、電氣、煤炭等為主要方向的優(yōu)勢產(chǎn)業(yè)，已建成4個國家級高新技術(shù)產(chǎn)業(yè)開發(fā)區(qū)， 2014年實(shí)現(xiàn)工業(yè)增加值1.333 1萬億元，第二產(chǎn)業(yè)占GDP的51.1%； 隨著經(jīng)濟(jì)發(fā)展水平不斷提高，娛樂業(yè)、服務(wù)業(yè)等第三產(chǎn)業(yè)也蓬勃發(fā)展，第三產(chǎn)業(yè)占GDP的37.2%[30]。

圖3 河北省GDP和人均GDP變化 圖4 河北省第一、二、三產(chǎn)業(yè)產(chǎn)值占GDP比重變化

表1 河北省耕地壓力指數(shù)重心、GDP重心、人均GDP重心移動距離

年份耕地壓力指數(shù)重心(km)GDP重心(km)人均GDP重心(km)1980～199022 954 038 21991～200044 928 369 812001～201011 9816 579 27

3.3 河北省耕地壓力與經(jīng)濟(jì)發(fā)展的重心軌跡演變

GDP和人均GDP能夠較好反映區(qū)域?qū)嶋H經(jīng)濟(jì)發(fā)展水平，利用式(3)計(jì)算河北省耕地壓力指數(shù)分別與GDP、人均GDP的地理聯(lián)系率G1和G2，進(jìn)而利用模型尋找河北省1980年、1990年、2000年和2010年的耕地壓力指數(shù)重心、GDP重心、人均GDP重心，測算40年來每一次重心移動的距離。

3.3.1 耕地壓力指數(shù)與GDP的重心軌跡擬合

1980年、1990年、2000年和2010年的河北省耕地壓力指數(shù)與GDP的地理聯(lián)系率G1均超過了90，分別為95.28、94.90、94.76和95.36，表明耕地壓力與經(jīng)濟(jì)發(fā)展實(shí)現(xiàn)的GDP產(chǎn)值的地理分布較為一致，在空間上存在較高的關(guān)聯(lián)。利用式(4)計(jì)算1980年、1990年、2000年和2010年耕地壓力指數(shù)重心和GDP重心變化(圖5)。40 年期間河北省的耕地壓力指數(shù)重心和GDP重心移動軌跡大體一致，均向東北方向偏移。在1980～1990年，耕地壓力指數(shù)重心和GDP重心移動位置均位于保定市，從重心移動距離和方向來看，耕地壓力指數(shù)重心向東北方向移動，且移動距離較大，超過20km，而GDP重心向西北方向發(fā)生4.03km的移動； 1991～2000年期間耕地壓力指數(shù)重心和GDP重心均發(fā)生東移跨越市行政界線，前者依然保持較長距離的移動至廊坊市，后者緊跟后向東移動至滄州市，此時(shí)耕地壓力指數(shù)重心和GDP重心不再處于同一市域； 2001～2010年耕地壓力指數(shù)重心和GDP重心繼續(xù)分別位于廊坊市和滄州市，并不斷向彼此移動靠近。40 年來耕地壓力和GDP的在移動方向和軌跡的一致性，一定程度上更加表明GDP與耕地壓力密切相關(guān)。

圖5 河北省耕地壓力指數(shù)重心與GDP重心移動軌跡 圖6 河北省耕地壓力指數(shù)重心與人均GDP重心移動軌跡

3.3.2 耕地壓力指數(shù)與人均GDP的重心軌跡擬合

1980年、1990年、2000年和2010年的河北省耕地壓力指數(shù)與人均GDP的地理聯(lián)系率G2比G1更高，均超過95，分別為97.58、98.41、97.80 和96.42，表明耕地壓力與人均GDP在地理空間上的分布也存在較高的關(guān)聯(lián)性。利用公式(4)計(jì)算1980年、1990年、2000年和2010年耕地壓力指數(shù)重心和人均GDP重心變化(圖6)，兩者均大致向東北方向移動。40年期間耕地壓力指數(shù)重心和人均GDP重心都在保定市和廊坊市之間移動，但耕地壓力指數(shù)重心總體移動距離較大，總計(jì)移動距離超過79km，呈現(xiàn)“增大—減少”的震蕩變化，并且在方向上出現(xiàn)“先東北，后東南”的移動軌跡； 而人均GDP重心的移動距離基本保持不變，呈現(xiàn)“東北—東南—東北”的東移態(tài)勢。在1980～2000年期間兩者分別沿自己的移動軌跡向東北偏移，在2001～2010年向彼此方位靠近，也反映經(jīng)濟(jì)水平的提高，人們收入的增加與區(qū)域內(nèi)耕地壓力關(guān)系密切。

3.4 耕地壓力與經(jīng)濟(jì)發(fā)展的空間關(guān)聯(lián)

通過測度耕地壓力指數(shù)與GDP、人均GDP的重心和地理聯(lián)系率是從整體上反映耕地壓力與經(jīng)濟(jì)發(fā)展的耦合程度，而地理集中度可以從區(qū)域內(nèi)各地級市的角度反映兩者的耦合關(guān)系。運(yùn)用式(6)(7)計(jì)算河北省11市耕地壓力、GDP、人均GDP的地理集中指數(shù)和耦合指數(shù)(表2)。

表2 河北省耕地壓力指數(shù)與GDP、人均GDP的地理集中度與耦合指數(shù)

地區(qū)1980年1990年2000年2010年RGR1R2I1I2RGR1R2I1I2RGR1R2I1I2RGR1R2I1I2唐山市0 252 252 040 110 121 161 901 620 610 721 262 542 160 500 581 413 072 090 460 68邯鄲市1 682 191 870 770 901 561 851 400 841 121 181 691 210 700 981 251 821 230 691 02保定市0 750 950 610 801 250 781 070 620 731 270 591 000 560 591 050 670 860 850 780 80承德市0 420 210 361 981 160 490 210 352 301 400 450 150 262 891 720 520 210 342 511 54廊坊市2 411 341 701 801 422 442 073 441 180 712 581 893 341 370 772 711 922 271 411 20石家莊市0 741 961 540 380 480 852 001 380 420 620 772 271 470 340 520 981 990 820 491 20秦皇島市2 291 202 931 910 782 211 312 951 690 753 211 252 812 571 143 431 102 633 111 31邢臺市1 631 041 091 581 491 441 171 081 221 320 981 111 000 880 971 050 900 901 161 16張家口市0 590 290 392 071 530 510 310 411 621 231 010 230 314 383 290 680 240 372 791 84滄州市1 631 011 041 601 561 451 070 971 361 501 101 251 180 870 930 871 521 970 570 44衡水市2 590 811 253 192 071 701 071 561 591 091 011 191 820 850 561 080 821 621 310 66 注：RG、R1、R2分別為耕地壓力指數(shù)、GDP、人均GDP的地理集中度，I1、I2分別為耕地壓力指數(shù)與GDP、人均GDP的耦合指數(shù)

3.4.1 耕地壓力指數(shù)、GDP、人均GDP的地理集中度分異

1980～2010年區(qū)間，耕地壓力指數(shù)地理集中度的高值區(qū)均出現(xiàn)在衡水市、廊坊市、秦皇島市，表現(xiàn)為耕地面積較少，但人口數(shù)量較多，糧食供給壓力較大； 保定市、張家口市、承德市為耕地壓力指數(shù)地理集中度的低值區(qū)，最高值與最低值相差4～5倍。GDP地理集中度的最高值均出現(xiàn)在石家莊市、邯鄲市、唐山市、廊坊市，且都是河北省重要的工業(yè)城市，經(jīng)濟(jì)水平發(fā)展程度較高； 承德市、張家口市、衡水市為GDP地理集中度的低值區(qū)，最高值與最低值相差近10倍。人均GDP地理集中度的高值區(qū)出現(xiàn)在秦皇島市、唐山市、邯鄲市、廊坊市、滄州市等城市； 而低值區(qū)為承德市、張家口市、保定市等區(qū)域，最高值與最低值的差值從1980年的7倍增長到2000年的12倍，再在2010年回歸到7倍。

3.4.2 耕地壓力指數(shù)與GDP、人均GDP的耦合指數(shù)分異

耕地壓力指數(shù)與GDP的耦合指數(shù)的高值區(qū)出現(xiàn)在秦皇島市、承德市、張家口市、衡水市，說明這些區(qū)域耕地壓力的增長與GDP的增長有著更密切的聯(lián)動關(guān)系； 低值區(qū)為唐山市、石家莊市，表明兩個城市GDP的增長帶來的耕地壓力的增長幅度比其他城市小， 1980～2010年期間最高值與最低值相差從30倍降至10倍左右。

耕地壓力指數(shù)與人均GDP的耦合指數(shù)的高值區(qū)在1980年和1990年出現(xiàn)在張家口市、滄州市、衡水市，而在2000年和2010年高值區(qū)變?yōu)榍鼗蕧u市、承德市、張家口市，說明近年來多市出現(xiàn)經(jīng)濟(jì)發(fā)展影響耕地壓力水平的現(xiàn)象； 低值區(qū)在1980～2010年期間基本都出現(xiàn)在秦皇島市、石家莊市、廊坊市、唐山市等區(qū)域，最高值與最低值相差從10倍降至3倍。從河北省各市耕地壓力與GDP、人均GDP的耦合指數(shù)來看， 40年間雖然高值、低值區(qū)域有較小變化，但總體趨勢是內(nèi)部區(qū)域的耦合指數(shù)高值與低值的差距逐漸縮小，逐漸趨近于1，耕地壓力與經(jīng)濟(jì)發(fā)展水平的協(xié)調(diào)一致程度逐漸提高。

3.4.3 耕地壓力與經(jīng)濟(jì)發(fā)展的耦合類型分異

圖7 河北省耕地壓力與經(jīng)濟(jì)發(fā)展耦合類型演變

4 結(jié)論與討論

該文利用耕地壓力指數(shù)模型測算河北省1978～2014年耕地壓力指數(shù)，利用地理聯(lián)系率、重心模型等方法擬合了河北省耕地壓力指數(shù)重心與GDP、人均GDP重心演變軌跡，進(jìn)而測算耕地壓力、GDP、人均GDP的地理集中度和耦合指數(shù)，分析其時(shí)空耦合類型。

(1)1978～2014年間河北省耕地壓力指數(shù)呈現(xiàn)下降趨勢。伴隨著農(nóng)業(yè)技術(shù)條件的進(jìn)步，糧食產(chǎn)量不斷提高，河北省的人均最小耕地面積逐漸減少，耕地壓力指數(shù)也隨之降低，從1978年的1.16下降至2014年的0.91，并在2010年之后持續(xù)保持低于1的態(tài)勢，出現(xiàn)因?yàn)楦孛娣e導(dǎo)致糧食供應(yīng)不足現(xiàn)象的可能性較低； 由于耕地產(chǎn)能、人口數(shù)量、城市發(fā)展等影響，河北省11市耕地壓力區(qū)域差異明顯，其中耕地壓力指數(shù)小于1的數(shù)量也從1980年的僅2個增加到2010年的7個，耕地壓力在空間分布上出現(xiàn)較大面積的降低。1978～2014年間河北省經(jīng)濟(jì)發(fā)展速度較快。改革開放以來，河北省大力發(fā)展工業(yè)，拓展高新技術(shù)等優(yōu)勢產(chǎn)業(yè)，GDP和人均GDP增長迅速，經(jīng)濟(jì)總量位居全國前列。期間河北省的第一產(chǎn)業(yè)占GDP的比重從28.52%降至11.7%，第二產(chǎn)業(yè)占GDP的比重持續(xù)保持在50%左右，第三產(chǎn)業(yè)蓬勃發(fā)展，在2014年已占GDP的37.30%，其產(chǎn)業(yè)結(jié)構(gòu)發(fā)生了較合理的改變。

(2)河北省耕地壓力與經(jīng)濟(jì)發(fā)展存在密切的聯(lián)系。1978～2014年期間河北省耕地壓力指數(shù)與GDP、人均GDP的地理聯(lián)系均較高，到2010分別達(dá)到95.36和96.42。隨著區(qū)域內(nèi)GDP和人均GDP的迅速增長，最小耕地面積和耕地壓力指數(shù)不斷減少，正是耕地保護(hù)與經(jīng)濟(jì)發(fā)展之間矛盾統(tǒng)一的牽制關(guān)系。經(jīng)濟(jì)發(fā)展推動耕作技術(shù)的革新、加速農(nóng)業(yè)發(fā)展的進(jìn)程，可以更好地利用耕地得到更高的糧食產(chǎn)量，并且經(jīng)濟(jì)發(fā)展之后可利用更多的社會經(jīng)濟(jì)資源開展土地整理整治項(xiàng)目，恢復(fù)因挖損、塌陷等造成的耕地?fù)p失； 但經(jīng)濟(jì)發(fā)展可能帶來侵占部分優(yōu)質(zhì)良田、污染耕地等現(xiàn)象。1978～2014年期間耕地壓力指數(shù)重心與GDP重心、人均GDP重心的移動軌跡保持較高的一致性，均向東北方向發(fā)生移動。就移動幅度來看，耕地壓力指數(shù)重心的偏移距離和角度大于GDP重心、人均GDP重心，說明耕地壓力指數(shù)受到影響的可能性更大，移動更具敏感性，而區(qū)域內(nèi)經(jīng)濟(jì)處于穩(wěn)定增長發(fā)展時(shí)期，更加具有穩(wěn)定性； 就重心偏移的演變路徑而言，GDP重心和人均GDP重心后于耕地壓力指數(shù)重心，且在2010年出現(xiàn)三者彼此集中的現(xiàn)象，說明在經(jīng)濟(jì)發(fā)展和城市建設(shè)過程中，確實(shí)占用一部分耕地，經(jīng)濟(jì)水平提高的同時(shí)，耕地壓力也隨之升高。

(3)河北省耕地壓力與經(jīng)濟(jì)發(fā)展耦合關(guān)系明顯，耦合類型隨著時(shí)間的推移發(fā)生演變。1980～2010年期間河北省各地級市的耕地壓力地理集中度與GDP地理集中度、人均GDP地理集中度在數(shù)值上的關(guān)聯(lián)度逐漸增高，各地級市的分異逐步均衡； 越來越多地級市的耕地壓力與經(jīng)濟(jì)發(fā)展的耦合指數(shù)逐漸趨近于1，在2010年底有6個城市的耦合指數(shù)在[0.90， 1.20]之間。但其中張家口市、承德市在絕大部分年度的耦合指數(shù)大于1.50，耕地壓力與經(jīng)濟(jì)發(fā)展的耦合程度非常高，說明當(dāng)?shù)貐^(qū)域在城市化和工業(yè)化發(fā)展過程中，應(yīng)進(jìn)一步加強(qiáng)耕地質(zhì)量建設(shè)，積極發(fā)展生態(tài)農(nóng)業(yè)，增強(qiáng)耕地抗災(zāi)能力，提升耕地產(chǎn)出水平，防范因耕地壓力過大而引發(fā)的糧食供給不足、耕地資源過度利用等問題。

該研究通過對耕地壓力與經(jīng)濟(jì)發(fā)展的重心移動軌跡，刻畫了河北省耕地壓力和經(jīng)濟(jì)發(fā)展的空間格局變化，為耕地壓力和經(jīng)濟(jì)發(fā)展水平關(guān)系的研究提供了新的視角，也為當(dāng)?shù)睾侠砝酶刭Y源，轉(zhuǎn)變經(jīng)濟(jì)發(fā)展方式提供了理論依據(jù)。但由于統(tǒng)計(jì)資料和收集數(shù)據(jù)的限制，借助蔡運(yùn)龍等人的觀點(diǎn)僅從糧食安全的角度去測定耕地壓力指數(shù)，并未考慮環(huán)境污染、生態(tài)破壞等因素與耕地壓力水平的關(guān)系； 在研究耕地壓力與經(jīng)濟(jì)發(fā)展的耦合關(guān)系時(shí)，僅利用GDP和人均GDP的值判斷區(qū)域經(jīng)濟(jì)發(fā)展的水平，未建立系統(tǒng)的指標(biāo)體系去測度經(jīng)濟(jì)發(fā)展的速度與質(zhì)量。今后需進(jìn)行更多的對當(dāng)?shù)馗刭Y源存在的污染、損毀等因素綜合評判區(qū)域耕地壓力水平； 對于經(jīng)濟(jì)發(fā)展水平的衡量應(yīng)根據(jù)相應(yīng)問題補(bǔ)充指標(biāo)。

[1] 許麗麗， 李寶林，袁燁城，等.2000～2010年中國耕地變化與耕地占補(bǔ)平衡政策效果分析.資源科學(xué)， 2015, 37(8): 1543～1551

[2] 徐明崗， 盧昌艾，張文菊，等.我國耕地質(zhì)量狀況與提升對策.中國農(nóng)業(yè)資源與區(qū)劃， 2016, 37(7): 8～14

[3] 封志明， 劉寶勤，楊艷昭.中國耕地資源數(shù)量變化的趨勢分析與數(shù)據(jù)重建： 1949～2003.自然資源學(xué)報(bào)， 2005, 20(1): 35～43

[4] 宋偉， 陳百明，劉琳.中國耕地土壤重金屬污染概況.水土保持研究， 2013, 20(2): 293～298

[5] 譚永忠， 吳次芳，王慶日，等.“耕地總量動態(tài)平衡”政策驅(qū)動下中國的耕地變化及其生態(tài)環(huán)境效應(yīng).自然資源學(xué)報(bào)， 2005, 20(5): 727～734

[6] 李玉平， 蔡運(yùn)龍.區(qū)域耕地—人口—糧食系統(tǒng)動態(tài)分析與耕地壓力預(yù)測——以河北省邢臺市為例.北京大學(xué)學(xué)報(bào)(自然科學(xué)版)， 2007, 43(2): 230～234

[7] 耿艷輝， 閔慶文，成升魁，等.涇河流域耕地—人口—糧食系統(tǒng)與耕地壓力指數(shù)時(shí)空分布.農(nóng)業(yè)工程學(xué)報(bào)， 2008, 24(10): 68～73

[8] 甘彩紅， 李陽兵，邵景安，等.三峽庫區(qū)腹地縣域耕地壓力研究——以奉節(jié)縣27個鄉(xiāng)鎮(zhèn)為例.資源科學(xué)， 2014, 36(7): 1365～1373

[9] Wesstr?m I，Joel A，Messing I.Controlled drainage and subirrigation-A water management option to reduce non-point source pollution from agricultural land.Agriculture Ecosystems & Environment， 2014， 198： 74～82

[10]Heathwaite A L，Quinn P F，Hewett C J M.Modeling and managing critical source areas of diffuse pollution from agricultural land using flow connectivity simulation.Journal of Hydrology， 2005， 304(1-4)： 446～461

[11]Olarewaju O E，Adetunji M T，Adeofun C O，et al.Nitrate and phosphorus loss from agricultural land：implications for nonpoint pollution.Nutrient Cycling in Agroecosystems， 2009, 85(1): 79～85

[12]王國強(qiáng)， 宋艷華.基于耕地質(zhì)量數(shù)量的河南省糧食生產(chǎn)能力研究.中國農(nóng)業(yè)資源與區(qū)劃， 2012, 33(1): 49～55

[13]劉玉， 郭麗英，劉彥隨.1980～2008年環(huán)渤海地區(qū)縣域糧食的時(shí)空動態(tài)變化及分析.農(nóng)業(yè)工程學(xué)報(bào)， 2012, 28(12): 230～236

[14]楊秉珣， 董廷旭，劉泉.四川省城鎮(zhèn)化進(jìn)程及耕地保護(hù)研究.中國農(nóng)業(yè)資源與區(qū)劃， 2016, 37(9): 149～152

[15]吳群， 郭貫成，萬麗平.經(jīng)濟(jì)增長與耕地資源數(shù)量變化：國際比較及其啟示.資源科學(xué)， 2006, 28(4): 45～51

[16]傅澤強(qiáng)， 蔡運(yùn)龍，楊友孝,等.中國糧食安全與耕地資源變化的相關(guān)分析.自然資源學(xué)報(bào)， 2001, 16(4): 313～319

[17]劉彥隨， 王介勇，郭麗英.中國糧食生產(chǎn)與耕地變化的時(shí)空動態(tài).中國農(nóng)業(yè)科學(xué)， 2009, 42(12): 4269～4274

[18]彭凌， 丁恩俊，謝德體.中國耕地?cái)?shù)量變化與耕地保護(hù)政策關(guān)系的實(shí)證分析.西南大學(xué)學(xué)報(bào)(自然科學(xué)版)， 2011, 33(11): 103～110

[19]龍開勝， 陳利根.耕地資源數(shù)量與經(jīng)濟(jì)發(fā)展和城鄉(xiāng)收入關(guān)系的計(jì)量分析.資源科學(xué)， 2007, 29(4): 139～145

[20]張夢琳， 陳利根，龍開勝.江蘇省耕地資源與經(jīng)濟(jì)發(fā)展計(jì)量與協(xié)調(diào)性分析.中國人口·資源與環(huán)境， 2009, 19(3): 82～86

[21]蔡運(yùn)龍， 傅澤強(qiáng)，戴爾阜.區(qū)域最小人均耕地面積與耕地資源調(diào)控.地理學(xué)報(bào)， 2002, 57(2): 127～134

[22]賈科利， 馬欣，張俊華.寧夏耕地壓力與社會經(jīng)濟(jì)發(fā)展耦合時(shí)空特征分析.干旱區(qū)地理， 2014, 37(4): 812～819

[23]陳志剛， 王青，黃賢金，等.長三角城市群重心移動及其驅(qū)動因素研究.地理科學(xué)， 2007, 27(4): 457～462

[24]Ellision G，Glaeser E L.Geographic concentration in U.S.Manufacturing industries：A dartboard approach.Journal of Urban Economics， 2000, 9(5): 48～67

[25]張定， 曹衛(wèi)東，朱勝清，等.安徽省物流與經(jīng)濟(jì)時(shí)空耦合研究.地域研究與開發(fā)， 2014, 3(33): 27～32

[26]王千， 金曉斌，阿依吐爾遜·沙木西，等.河北省糧食產(chǎn)量空間格局差異變化研究.自然資源學(xué)報(bào)， 2010，25(9)： 1525～1535

[27]范秋梅， 蔡運(yùn)龍.基于糧食安全的區(qū)域耕地壓力研究——以遼寧省為例.地域研究與開發(fā)， 2010, 29(5): 110～113

[28]熊昌盛， 韋仕川，欒喬林，等.基于Moran′s I分析方法的耕地質(zhì)量空間差異研究——以廣東省廣寧縣為例.資源科學(xué)， 2014, 36(10): 2066～2074

[29]唐華俊， 李哲敏.基于中國居民平衡膳食模式的人均年糧食需求量研究.中國農(nóng)業(yè)科學(xué)， 2012, 45(11): 2315～2327

[30]河北省統(tǒng)計(jì)局. 湖北省統(tǒng)計(jì)年鑒2015.北京：中國統(tǒng)計(jì)出版社， 2015