唐秀美，潘瑜春※，劉 玉

(1.北京農(nóng)業(yè)信息技術(shù)研究中心，北京 100097； 2.國家農(nóng)業(yè)信息化工程技術(shù)研究中心，北京 100097；3.農(nóng)業(yè)部農(nóng)業(yè)信息技術(shù)重點(diǎn)實(shí)驗(yàn)室，北京 100097)

0 引言

耕地資源是人類生存和發(fā)展不可替代的重要資源。在我國，耕地資源承載著保證糧食安全和社會穩(wěn)定的功能，因此，我國政府高度重視耕地保護(hù)，將其列入基本國策，采取“世界上最嚴(yán)格的耕地保護(hù)政策”[1]。耕地資源是一種公眾資源，人類對耕地的需求呈現(xiàn)多樣性、層次性、區(qū)域性特點(diǎn)，在不同的區(qū)域耕地價值需求主體側(cè)重點(diǎn)不同，耕地價值的體現(xiàn)也隨之有所側(cè)重[2]。目前許多學(xué)者從不同角度對耕地價值進(jìn)行分類[3-5]，并對耕地的生態(tài)價值、社會保障價值、經(jīng)濟(jì)價值進(jìn)行理論探討和定量化研究[6-8]。2000年以來，隨著全球城市化進(jìn)程的加快，生態(tài)系統(tǒng)服務(wù)稀缺性變得越來越突出，耕地生態(tài)系統(tǒng)的多功能性受到了空前的重視，國外有關(guān)農(nóng)地價值，特別是耕地生態(tài)價值的研究成為熱點(diǎn)[9-11]。目前耕地生態(tài)價值評估中，采用的方法主要包括價值當(dāng)量修正法、替代市場法和條件價值評估法等方法[13]。已有研究中，耕地生態(tài)功能類型多樣，未有公認(rèn)的分類體系，且測算方法多樣，多以耕地面積為基礎(chǔ)，計(jì)算不同年份之間的生態(tài)價值差異，而對同一區(qū)域內(nèi)部空間之間的差異分析研究較少。基于此，該研究結(jié)合北京市耕地所處位置和影響因素，構(gòu)建了適合區(qū)域的耕地生態(tài)價值評估模型，以期合理測算耕地生態(tài)功能，分析區(qū)域內(nèi)部生態(tài)價值差異和變化規(guī)律，從而為區(qū)域耕地保護(hù)和多功能利用奠定基礎(chǔ)。

1 研究區(qū)概況與數(shù)據(jù)基礎(chǔ)

1.1 研究區(qū)概況

北京市位于華北平原西北隅，北緯39°38′～41°51′，東經(jīng)115°25″～117°30′，東西寬約160km，南北長約170km，東南距渤海約150km，轄16個區(qū)。北京是中國的首都，是全國的政治經(jīng)濟(jì)文化中心，其土地利用呈現(xiàn)明顯的圈層分布特點(diǎn)。北京市耕地資源稀缺，發(fā)揮耕地的多功能性對于北京市可持續(xù)發(fā)展和建設(shè)宜居城市有重要意義。

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

數(shù)據(jù)來源為中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心的土地利用數(shù)據(jù)，該數(shù)據(jù)利用陸地衛(wèi)星TM的數(shù)字圖像作為原始數(shù)據(jù)源，解譯分類的定位誤差為小于兩個像元(60m)。數(shù)據(jù)時期為1980年、1995年、2000年和2013年，數(shù)據(jù)格式為Shpfile格式。該研究根據(jù)研究目的提取了耕地作為研究對象。

2 北京市耕地生態(tài)價值分類與評估方法

耕地資源的生態(tài)價值是指耕地在調(diào)節(jié)氣候、凈化與美化環(huán)境、防止水土流失、維護(hù)物種多樣性等方面所具有的功能價值[1]。結(jié)合已有研究和北京市的實(shí)際情況，該研究中耕地生態(tài)價值包括氣體調(diào)節(jié)功能、凈化環(huán)境功能、涵養(yǎng)水源功能、固土保肥功能和營養(yǎng)物質(zhì)循環(huán)功能，并對不同的功能進(jìn)行了細(xì)分，在計(jì)算單位面積耕地的生態(tài)價值基礎(chǔ)上，分析不同區(qū)域耕地生態(tài)價值空間分布狀況及變化規(guī)律。

2.1 氣體調(diào)節(jié)功能

耕地生態(tài)系統(tǒng)大氣調(diào)節(jié)功能價值的評價，首先要計(jì)算耕地地塊生態(tài)系統(tǒng)各種農(nóng)作物的年凈生物量[14]，在此基礎(chǔ)上，根據(jù)光合作用和呼吸作用方程式得知，生態(tài)系統(tǒng)每生產(chǎn)1.00g植物干物質(zhì)能固定1.62g CO2，釋放1.20g O2[15]，據(jù)此，采用固碳制氧成本法計(jì)算耕地大氣調(diào)節(jié)價值。

E調(diào)=ECO2+EO2

(1)

ECO2=PiTMC1

(2)

EO2=PiVC2

(3)

Pi=SP

(4)

其中，E調(diào)為耕地大氣調(diào)節(jié)功能價值；ECO2為耕地固碳價值；EO2為耕地釋氧價值；Pi為區(qū)域耕地第一凈生產(chǎn)力總量；T為1.00g植物干物質(zhì)固定CO2重量，取值為1.62g；M為CO2折算成碳的系數(shù)為3/11；C1為碳釋放成本價值，取值為638.65元/t；C2為釋氧的成本價，以工業(yè)制氧成本400元/t計(jì)算；V為1.00g植物干物質(zhì)釋放O2重量，取值為1.20g；S為耕地面積(hm2)；P為耕地單位面積第一凈生產(chǎn)力，取值為4.829t/(hm2·年)[16]。據(jù)此，計(jì)算1hm2耕地的大氣調(diào)節(jié)功能價值為：

E調(diào)=SPTMC1+SPVC2=1hm2×4 829kg×1.62kg×3/11×0.64元/kg+1hm2×4 829kg×1.20kg×0.4元/kg=1 365.47+2 317.92=3 683.39元/hm2

2.2 凈化環(huán)境功能

耕地凈化環(huán)境功能包括吸收有害氣體、凈化大氣、消納廢棄物等。該研究計(jì)算耕地凈化環(huán)境價值主要考慮凈化大氣和消納畜禽糞便兩個方面。

E凈=E氣+E消

(5)

式(5)中，E凈為耕地凈化環(huán)境價值；E凈為凈化大氣功能價值；E消為耕地消納廢棄物價值。

2.2.1 凈化大氣功能

耕地對大氣的凈化功能主要是指耕地可以依靠自身特殊的結(jié)構(gòu)和功能，通過吸收、過濾、阻隔和分解等生理生化過程[17]。凈化大氣功能主要計(jì)算耕地減少空氣中硫化物、氮化物、鹵素等有害物質(zhì)的含量和滯塵所實(shí)現(xiàn)的價值。該研究運(yùn)用恢復(fù)費(fèi)用法計(jì)算耕地吸收SO2、HF、NOx和滯塵能力的價值[14]。

E氣=ESO2+EHF+ENOX+Ez

(6)

ESO2=SU1C3

(7)

EHF=SU2C4

(8)

ENOX=SU3C5

(9)

EZ=SU4C6

(10)

其中，ESO2為耕地吸收SO2的價值；EHF為耕地吸收HF的價值；ENOX為耕地吸收NOx的價值；Ez為耕地滯塵的價值； 結(jié)合已有研究，確定主要的參數(shù)值[18-19]，U1為單位面積耕地資源吸收SO2量，取值為45kg/(hm2·年)； C3為治理SO2單位成本，取值為3 000元/t；U2為單位面積耕地資源吸收HF的量，取值為0.38kg/(hm2·年)；C4為治理HF單位成本，取值為600元/t；U3為單位面積耕地資源吸收NOx的量，取值為33.8kg/(hm2·年)；C5為治理NOx單位成本，取值為1.6萬元/t；U4為單位面積耕地資源的滯塵能力的量，取值為0.95kg/(hm2·年)，C6為治理粉塵單位成本，取值為170元/t。據(jù)此，計(jì)算1hm2耕地的凈化大氣功能價值為：

E凈=ESO2+EHF+ENOX+Ez=1hm2×45kg/(hm2·年)×3元/kg+1hm2×0.38kg/(hm2·年)×0.6元/kg+1hm2×33.8kg/(hm2·年)×16元/kg+1hm2×0.95kg/(hm2·年)×0.17元/kg=135+0.228+540.80+0.16=676.19元/hm2

2.2.2 消納廢棄物功能

耕地對畜禽糞便的消納降解可以減少畜禽糞便污染環(huán)境，該研究中主要考慮耕地對氮、磷的消納能力[20]，確定耕地消納廢棄物的計(jì)算方法：

E消=S×(∑Vj×Cj)

(11)

式(11)中，Vj為耕地年消納氮、磷的量，參考已有研究[21-22]，北京市耕地年消納氮、磷的限量確定為200和40kg/hm2，北京市有機(jī)肥的施用比例確定為30%，即北京市耕地消納氮、磷的量取值為60和12kg/hm2；Cj為j種肥料的價格，根據(jù)中國化肥網(wǎng)的統(tǒng)計(jì)數(shù)據(jù)，尿素的平均價格為1 450元/t，過磷酸鈣的價格為1 500元/t，氯化鉀的價格為1 900元/t，根據(jù)純元素的折算率，氮為46%，磷為17%，鉀為55%，確定全氮、全磷、全鉀的價格分別為3 152元/t、8 823元/t、3 455元/t[16]，據(jù)此，計(jì)算1hm2耕地的消納廢棄物價值為：

E消=1hm2×(60kg/hm2×3.15元/kg+12kg/hm2×8.82元/kg)=294.84元/hm2

1hm2耕地的凈化環(huán)境的價值為：

E凈=676.19元/hm2+294.84元/hm2=971.03元/hm2

2.3 涵養(yǎng)水源功能

耕地涵養(yǎng)水源功能主要指耕地對降水的截留、吸收和貯存，將地表水轉(zhuǎn)換為地表徑流或地下水的作用[16]。該研究根據(jù)水庫工程的蓄水成本(替代工程法)確定耕地的水源涵養(yǎng)功能。

2.3.1 涵養(yǎng)水源量

以耕地區(qū)域水量平衡法計(jì)算耕地涵養(yǎng)水源量。

W水=(G-Z)×S=G×(1-n)×S

(12)

式(12)中，W水為涵養(yǎng)水源量(t)；G為平均降水量(mm/年)，取值為北京市年平均降水量585mm；Z為平均蒸散量(mm/年)；n為耕地資源上作物年蒸騰量占總降雨的比例，參考已有研究，取值為0.84[23]。據(jù)此，計(jì)算1hm2耕地的涵養(yǎng)水源量為：

W水=626mm/年×(1-0.84)×1hm2=936t

2.3.2 涵養(yǎng)水源價值

耕地涵養(yǎng)水源的價值根據(jù)水庫工程的蓄水成本(替代工程法)確定。

E涵=W水×C庫

(13)

式(13)中，E涵為耕地涵養(yǎng)水源的價值(元)；W水為涵養(yǎng)水源量(t/年)；C庫為水庫造價(元/t)，根據(jù)中國水利年鑒的水庫庫容造價數(shù)據(jù)，結(jié)合歷年價格指數(shù)[24]，確定2015年的北京市水庫庫容造價為6.48元/m3，采用公式(13)計(jì)算涵養(yǎng)水源價值。據(jù)此，計(jì)算1hm2耕地的涵養(yǎng)水源價值為：

E涵=936t×6.48元/m3=6 065.28元/hm2

2.4 固土保肥功能

耕地中農(nóng)作物覆蓋在土壤上，一方面可以減少土壤侵蝕量，另一方面可以保持因土壤侵蝕而使泥沙淤積在水庫、江河、湖泊的間接經(jīng)濟(jì)價值。因此，耕地的固土保肥功能有保持土壤肥力、減少河流湖泊泥沙淤積和固定肥力(化肥和有機(jī)肥)3個方面[24]。研究采用替代市場法和影子工程法計(jì)算耕地減少這3方面損失的價值。

E固=E侵+E淤+E肥+E有

(14)

式(14)中，E固為耕地土壤固土保肥價值；E侵為耕地減少土壤侵蝕價值；E淤為耕地減少河流湖泊泥沙淤積價值；E肥為耕地固定肥力(元)；E有為減少有機(jī)質(zhì)損失價格。

2.4.1 土壤侵蝕總量

耕地土壤侵蝕總量的計(jì)算公式為：

W土=S×M

(15)

式(15)中，W土為耕地減少土壤侵蝕總量(t)；M為耕地侵蝕模數(shù)[t/(hm2·年)]，根據(jù)有效林地的侵蝕模數(shù)[25]，結(jié)合北京市水土保持公報，北京市耕地侵蝕模數(shù)平均值為25t/(hm2·年)[26]。

W土=1hm2×25 t/(hm2·年)=25t

2.4.2 減少侵蝕價值

將由計(jì)算出的耕地減少土壤侵蝕總量除以土地耕作層的平均厚度，即得耕地減少土地資源損失的面積WA。

(16)

式(16)中，WA為減少土地資源損失面積(hm2)；W土為耕地減少土壤侵蝕總量(t)；T為土壤的平均密度，取值為1.3g/cm3；L為土地耕作層的平均厚度(m)，根據(jù)已有研究，取值為0.5 m[16]。

E侵=WA×C耕

(17)

以耕地生產(chǎn)用地的年平均收益，作為耕地減少廢棄土地的機(jī)會成本，計(jì)算耕地減少土壤侵蝕總量的價值。E耕為耕地減少土壤侵蝕總量的價值，C耕為耕地的年平均收益，根據(jù)北京市統(tǒng)計(jì)年鑒及中國農(nóng)村統(tǒng)計(jì)年鑒數(shù)據(jù)，確定北京市2015年耕地平均收益為3.2萬元/hm2。據(jù)此，計(jì)算1hm2耕地的減少侵蝕價值為：

E侵=25×103kg/1.3×103kg/m3/0.5m×3.2元/m2=123.08元/hm2

2.4.3 減少淤積價值

按照我國主要流域的泥沙運(yùn)動規(guī)律，一般土壤侵蝕流失的泥沙有24%淤積于水庫、江河、湖泊，這部分泥沙直接造成了水庫、江河、湖泊蓄水量的下降，在一定程度上增加了干旱、洪澇災(zāi)害發(fā)生的機(jī)會，另有33%滯留， 37%入海[27]。該文僅考慮淤積于水庫、江河湖泊的24%，即每年減少泥沙淤積的經(jīng)濟(jì)價值。按照采用影子工程法計(jì)算減少河流湖泊泥沙淤積的價值：

E淤=W土×24%×C庫

(18)

據(jù)此，計(jì)算1hm2耕地的減少淤積的價值為：

E淤=25t×24%×6.48元/t=38.88元/hm2

2.4.4 固定肥力價值

固定肥力價值主要根據(jù)減少土壤流失中的土壤肥力價值計(jì)算。其中，有機(jī)質(zhì)損失的價值采用耕地可增加的薪柴的費(fèi)用確定。

E肥=∑Rj×Cj×W土

(19)

E有=W有×T有×C有

(20)

式(19)(20)中，Rj為單位侵蝕物中第j種養(yǎng)分元素的含量(g/kg)； 結(jié)合已有研究，確定北京市土壤中有機(jī)質(zhì)、氮、磷、鉀含量為：有機(jī)質(zhì)15.3g/kg，氮0.9g/kg，琳0.8g/kg，鉀21g/kg[28-29]；T有為薪材轉(zhuǎn)化成有機(jī)質(zhì)的比例為2： 1；C有為薪材的平均價格取值為211元/t。據(jù)此，計(jì)算1hm2耕地的固土保肥價值為：

E肥=25 000kg×(0.9 g/kg×10-3×3.15元/kg+0.8 g/kg×10-3×8.82元/kg+21 g/kg×10-3×3.46元/kg)=2 063.78元/hm2

E有=25 000kg×15.3g/kg×10-3×2×0.211元/kg=161.42元/hm2

E固=123.08+38.88+2 063.78+161.42=2 387.16元/hm2

2.5 營養(yǎng)物質(zhì)循環(huán)功能

耕地生態(tài)系統(tǒng)的營養(yǎng)物質(zhì)循環(huán)主要在生物庫、凋落物庫和土壤庫之間進(jìn)行。對于農(nóng)田生態(tài)系統(tǒng)來說，因其凋落物量較小，生物與土壤之間的養(yǎng)分交換過程是最主要的過程[30]，該研究只考慮土壤庫和生物庫。該文對參與評價的生物庫和土壤庫中的營養(yǎng)元素僅考慮含量相對較多的氮、磷、鉀[31]。

E營=E生+E土

(21)

(22)

(23)

其中，E營為耕地營養(yǎng)物質(zhì)循環(huán)功能價值；E生為生物庫參與營養(yǎng)元素循環(huán)的價值；E生為生物庫中營養(yǎng)物質(zhì)循環(huán)的總價值；Cnj為第j類農(nóng)產(chǎn)品生物質(zhì)中含氮的百分比；Cpj為第j類農(nóng)產(chǎn)品生物質(zhì)中含磷的百分比；Ckj為第j類農(nóng)產(chǎn)品生物質(zhì)中含鉀的百分比；Pn、Pp、Pk分別對應(yīng)于氮、磷、鉀的市場價格；E土為土壤庫中營養(yǎng)物質(zhì)循環(huán)的價值；Mj為第j類農(nóng)產(chǎn)品土壤庫總量；Snj為第j類農(nóng)產(chǎn)品土壤庫中含氮的百分比；Spj為第j類農(nóng)產(chǎn)品土壤庫中含磷的百分比；Skj為第j類農(nóng)產(chǎn)品土壤庫中含鉀的百分比； 根據(jù)已有研究[15-16]，生物庫中氮、磷、鉀的營養(yǎng)元素含量比例取平均值，分別為3.09%， 0.74%和3.28%，土壤庫中氮、磷、鉀的百分含量換算后為0.09%、0.08%和2.10%；fN、fP、fK分別為氮、磷、鉀在土壤中的周轉(zhuǎn)率，其值分別為0.08、0.01和0.01[32]。據(jù)此計(jì)算1hm2耕地營養(yǎng)物質(zhì)循環(huán)價值為：

E營=1hm2×4 829kg×(3.09%×3.15元/kg+0.74%×8.82元/kg+3.28%×3.46元/kg)+1×104m2×0.5m×1.3×103kg/m3×(0.09%×0.08×3.15元/kg+0.08%×0.01×8.82元/kg+2.10%×0.01×3.46元/kg)=1 333.24+6 655.74=7 988.98元/hm2

圖1 北京市單位面積耕地生態(tài)價值組成情況

3 結(jié)果分析

3.1 北京市單位耕地生態(tài)價值

通過對北京市耕地生態(tài)價值的計(jì)算結(jié)果，分析耕地生態(tài)價值組成情況，北京市耕地生態(tài)價值總值為2.109 6萬元/hm2，其中，價值最高的為營養(yǎng)物質(zhì)循環(huán)功能，為7 988.98元/hm2，其次是涵養(yǎng)水源功能，為6 065.28元/hm2，氣體調(diào)節(jié)功能是3 683.39元/hm2，固土保肥功能為2 387.16元/hm2，最低的為凈化大氣功能，價值為971.03元/hm2。

表1 北京市耕地生態(tài)價值變化

年份耕地面積(萬hm2)氣體調(diào)節(jié)功能(億元)凈化環(huán)境功能(億元)涵養(yǎng)水源功能(億元)固土保肥功能(億元)營養(yǎng)物質(zhì)循環(huán)功能(億元)總價值(億元)1980580858214056435231387464012254199544687816464342710106735709427200048689717934732953116238901027220133645261343354221187029127690

3.2 生態(tài)總價值變化

結(jié)合北京市1980年、1995年、2000年和2013年的耕地分布狀況，計(jì)算北京市耕地生態(tài)價值變化情況，結(jié)果如圖2、表1。1980～2013年，耕地生態(tài)價值分布呈現(xiàn)明顯的從中心城區(qū)到遠(yuǎn)郊區(qū)逐漸降低的趨勢， 1980～2013年，耕地生態(tài)總價值共減少了45.64億元，其中，營養(yǎng)物質(zhì)循環(huán)功能減少最高，達(dá)到17.28億元，其次是涵養(yǎng)水源功能，減少13.12億元，氣體調(diào)節(jié)功能減少7.97億元，固土保肥功能減少5.17億元，凈化環(huán)境功能減少2.1億元。

3.3 分區(qū)縣耕地生態(tài)價值變化

表2 北京市各區(qū)耕地生態(tài)價值變化 億元

由表2分析各區(qū)耕地生態(tài)價值的變化情況， 1980～2013年，各區(qū)縣耕地生態(tài)價值都呈現(xiàn)了明顯的下降趨勢，特別是1980～1995年，減少幅度較高。從遠(yuǎn)郊區(qū)到城區(qū)，耕地生態(tài)價值的減少幅度逐漸增加，其中，豐臺、石景山、房山、順義的耕地生態(tài)價值逐年減少； 而朝陽、海淀、門頭溝、通州、昌平、大興、懷柔、平谷、密云和延慶等區(qū)的耕地生態(tài)價值都是1980～1995年降低，1995～2000年增加， 2000～2013年又增加的趨勢。分區(qū)統(tǒng)計(jì)發(fā)現(xiàn)，從1980年到2013年，大興區(qū)的耕地生態(tài)價值減少數(shù)量最多，減少了6.01億元，其次是朝陽區(qū)，減少了5.9億元，昌平區(qū)和順義區(qū)也減少了5.63億元和5.48億元，而石景山減少的數(shù)量較少，僅為0.32億元，門頭溝減少了0.73億元，從減少比例計(jì)算，減少比例最高的為石景山區(qū)，達(dá)到98.94%，其次是朝陽區(qū)，減少比例為98.17%，延慶的減少比例最低，為15.18%，平谷區(qū)為17.91%，密云區(qū)為19.37%。從各區(qū)縣耕地不同生態(tài)功能變化來說，各區(qū)縣營養(yǎng)物質(zhì)循環(huán)功能減少最高，其次是涵養(yǎng)水源功能和氣體調(diào)節(jié)功能，固土保肥功能和凈化環(huán)境功能減少數(shù)量最低。

4 討論與結(jié)論

4.1 討論

對耕地價值的認(rèn)識經(jīng)歷了由無價到單一的經(jīng)濟(jì)價值核算，再到綜合價值核算的發(fā)展過程。近年來，由于耕地生態(tài)價值評估越來越受到學(xué)術(shù)界研究的重視，耕地生態(tài)價值本身的影響因素復(fù)雜多樣且多難以定量，學(xué)者采用各種方法試圖對其進(jìn)行合理核算，評估方法多樣?？傮w而言，耕地生態(tài)價值的核算應(yīng)該結(jié)合其所在區(qū)位、經(jīng)濟(jì)發(fā)展?fàn)顩r等合理確定類型及估算方法。該研究中結(jié)合北京市的實(shí)際情況，確定了北京市耕地5類生態(tài)價值類型和各類型中不同功能的核算方法，研究對于完善耕地生態(tài)價值核算方法，分析不同時間和區(qū)域耕地生態(tài)價值分布情況有重要意義。將來應(yīng)進(jìn)一步對耕地多功能價值估算方法進(jìn)行深入研究，特別是耕地多功能價值的時空變化規(guī)律及驅(qū)動力是研究重點(diǎn)。

4.2 結(jié)論

該文在對北京市耕地生態(tài)價值進(jìn)行定義的基礎(chǔ)上，根據(jù)北京市的實(shí)際情況，確定了北京市耕地生態(tài)價值類型，在對不同類型的生態(tài)價值模型構(gòu)建的基礎(chǔ)上，以北京市4期土地利用現(xiàn)狀圖為基礎(chǔ)，分析了不同年代北京市耕地生態(tài)價值的變化情況。在對分類型耕地生態(tài)價值評估中，完善了目前已有研究中耕地生態(tài)價值的計(jì)算方法，在氣體調(diào)節(jié)功能中，計(jì)算了固碳和釋氧價值； 在凈化環(huán)境功能中，計(jì)算了凈化大氣功能和消納廢棄物的價值； 在固土保肥功能中，計(jì)算了減少土壤侵蝕、減少淤積和固定肥力的價值。相對全面地計(jì)算了不同功能的價值類型。同時，研究彌補(bǔ)了已有研究多重視數(shù)量分析而忽略空間分布的問題，對于分區(qū)域耕地保護(hù)有重要意義。

研究發(fā)現(xiàn)， 1980～2013年，耕地生態(tài)價值從遠(yuǎn)郊區(qū)到中心城區(qū)逐漸降低，且各區(qū)耕地生態(tài)價值都呈現(xiàn)了明顯的下降趨勢，這主要是因?yàn)楸本┦薪?0年來城市擴(kuò)展占用了大量的耕地，使得耕地面積大量減少，造成北京市1980～2013年總體耕地生態(tài)價值降低，在2000年前，由于北京市進(jìn)行了土地開發(fā)等活動，增加了部分耕地，所以1997～2000年耕地生態(tài)價值有所提高。而從分區(qū)看，耕地生態(tài)價值減少較高的區(qū)域位于朝陽區(qū)和大興、順義、昌平等近郊平原區(qū)，而平谷、密云、延慶等生態(tài)涵養(yǎng)區(qū)的生態(tài)價值較少，數(shù)量和比例相對低，說明北京市近30年來主要的建設(shè)用地擴(kuò)展位于北京市東南部的平原區(qū)。總體而言，北京市的耕地生態(tài)價值減少較大，且在空間上不均衡，北京市應(yīng)該進(jìn)一步加強(qiáng)耕地保護(hù)，提高耕地利用效率，并積極促進(jìn)耕地多功能應(yīng)用。

[1] 李佳，南靈.耕地資源價值內(nèi)涵及測算方法研究——以陜西省為例.干旱區(qū)資源與環(huán)境， 2010， 24(9)： 73～19

[2] 李翠珍，孔祥斌，孫憲海.北京市耕地資源價值體系及價值估算方法.地理學(xué)報， 2008， 63(3)： 321～329

[3] 蔡運(yùn)龍，霍雅勤.中國耕地價值重建方法與案例研究.地理學(xué)報， 2006， 61(10)： 1084～1092

[4] 蔡銀鶯，李曉云，張安錄.湖北省農(nóng)地資源價值研究.自然資源學(xué)報， 2007， 22(1)： 121～130

[5] 武燕麗.農(nóng)用土地資源價值測度方法研究.太原：山西農(nóng)業(yè)大學(xué)， 2005

[6] 段瑞娟，郝晉珉，張潔瑕.北京區(qū)位土地利用與生態(tài)服務(wù)價值變化研究.農(nóng)業(yè)工程學(xué)報， 2006， 22(9)： 21～28

[7] 陳麗，曲福田，師學(xué)義.耕地資源社會價值測算方法探討——以山西省柳林縣為例.資源科學(xué)， 2006， 28(6)： 86～90

[8] 李景剛，歐名豪，張效軍，等.耕地資源價值重建及其貨幣化評價——以青島市為例.自然資源學(xué)報， 2009， 24(11)： 1870～1880

[9] Christopher D.Ives，Dave Kendal.Values and attitudes of the urban public towards peri-urban agricultural land.land Use Policy， 2013， 34(12)： 80～90

[10]Sutton N J，Choc S，Armswortha P R.A reliance on agricultural land values in conservation planning alters the spatial distribution of priorities and overestimates the acquisition costs of protected areas.Biological Conservation， 2016， 194:2～10

[11]Johanna Choumert，Pascale Phélinas.Determinants of agricultural land values in Argentina.Ecological Economics， 2015， 110:134～140

[12]Ustaoglu E，Perpina Castillo C，Jacobs-Crisioni C.Economic evaluation of agricultural land to assess land use changes.Land Use Policy， 2016， 56:125～146

[13]唐秀美，陳百明，劉玉，等.耕地生態(tài)價值評估研究進(jìn)展分析．農(nóng)業(yè)機(jī)械學(xué)報， 2016， 47(9)： 256～265

[14]崔朝偉，許學(xué)工.北京市農(nóng)用地生態(tài)系統(tǒng)服務(wù)價值評估.生態(tài)經(jīng)濟(jì)(學(xué)術(shù)版)， 2007，(2)： 338～340， 358

[15]曹志宏，郝晉珉，梁流濤.黃淮海地區(qū)耕地資源價值核算.干旱區(qū)資源與環(huán)境， 2009， 23(9)： 7～12

[16]湯進(jìn)華，陳志，朱俊成，等.武漢城市圈耕地資源生態(tài)服務(wù)價值核算.中國農(nóng)學(xué)通報， 2015， 31(4)： 237～244

[17]白楊，歐陽志云，鄭華，等.海河流域農(nóng)田生態(tài)系統(tǒng)環(huán)境損益分析.應(yīng)用生態(tài)學(xué)報， 2010， 21(11)： 2938～2945

[18]謝高地，肖玉，甄霖，等.我國糧食生產(chǎn)的生態(tài)服務(wù)價值研究.中國生態(tài)農(nóng)業(yè)學(xué)報， 2005， 13(3)： 10～13

[19]白偉.糧食生產(chǎn)中的自然資源價值研究.北京：中國農(nóng)業(yè)大學(xué)， 2005

[20]楊志新，鄭大瑋，文化.北京郊區(qū)農(nóng)田生態(tài)系統(tǒng)服務(wù)功能價值的評估研究.自然資源學(xué)報， 2005， 20(4)： 564～571

[21]潘瑜春，孫超，劉玉，等.基于土地消納糞便能力的畜禽養(yǎng)殖承載力.農(nóng)業(yè)工程學(xué)報， 2015， 31(4)： 232～239

[22]武蘭芳，歐陽竹.基于農(nóng)田氮磷收支的區(qū)域養(yǎng)殖畜禽容量分析——以山東禹城為例.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報， 2009， 28(11)： 2277～2285

[24]中國森林生態(tài)系統(tǒng)服務(wù)評估項(xiàng)目組.中國森林生態(tài)系統(tǒng)服務(wù)評估．北京：中國林業(yè)出版社， 2010

[25]田石磊，廖超英，王小翠.藍(lán)田縣森林生態(tài)系統(tǒng)服務(wù)價值的評價.西北農(nóng)林科技大學(xué)學(xué)報(自然科學(xué)版)， 2009， 37(5)： 133-138

[26]歐陽志云，王效科.中國陸地生態(tài)系統(tǒng)服務(wù)功能及其生態(tài)經(jīng)濟(jì)價值的初步研究.生態(tài)學(xué)報， 1999， 19(5)： 607-613

[27]吳剛，馮宗煒，王效科，等.黃淮海平原農(nóng)林生態(tài)系統(tǒng)N、P、K營養(yǎng)元素循環(huán)——以泡桐—小麥、玉米間作系統(tǒng)為例.應(yīng)用生態(tài)學(xué)報， 1993， 4(2)： 141～145

[28]葉文虎，魏斌，仝川，等，城市生態(tài)補(bǔ)償能力衡量和應(yīng)用.中國環(huán)境科學(xué)， 1998， 18(4)： 298～301

[29]王克如，李少昆，曹連莆，等.新疆高產(chǎn)棉田氮、磷、鉀吸收動態(tài)及模式初步研究.中國農(nóng)業(yè)科學(xué)， 2003， 36(7)： 775～780

[30]歐陽志云，趙同謙，趙景柱，等.海南島生態(tài)系統(tǒng)生態(tài)調(diào)節(jié)功能及其生態(tài)經(jīng)濟(jì)價值研究.應(yīng)用生態(tài)學(xué)報， 2004， 15(8)： 1395～1402

[31]趙同謙，歐陽志云，鄭華，等.中國森林生態(tài)系統(tǒng)服務(wù)功能及其價值評價.自然資源學(xué)報， 2004， 19(4)： 480～491

[32]魯如坤，劉鴻翔，聞大中，等.我國典型地區(qū)農(nóng)業(yè)生態(tài)系統(tǒng)養(yǎng)分循環(huán)和平衡研究Ⅳ.農(nóng)田養(yǎng)分平衡的評價方法和原則.土壤通報， 1996,(5)： 197～199