International Journal of Ecology
Vol.06 No.02(2017), Article ID:20677,12 pages
10.12677/IJE.2017.62006

A Dynamic Evaluation of Agro-Ecosystem at Local Levels—Using Anyang City, Henan Province 2008-2014 as an Example

Feihu Gao1, Xiaoshuai Zhao2, Xiaoguang Zhang1,3, Shimei Li4, Yicheng Wang1*

1College of Resources and Environment, Qingdao Agricultural University, Qingdao Shandong

2Zouping County Environmental Protection Bureau, Zouping Shandong

3State Key Laboratory of Soil and Sustainable Agriculture, Chinese Academy of Sciences, Nanjing Jiangsu

4College of Landscape Architecture and Forestry, Qingdao Agricultural University, Qingdao Shandong

*通讯作者。

Received: Apr. 5th, 2017; accepted: May 19th, 2017; published: May 26th, 2017

ABSTRACT

The research on agro-ecosystem evaluation so far has largely been focusing on large-scale areas such as north China or even the whole China territory; little attention has been paid to local city or county levels, less to dynamic changes of agro-ecosystem values over time. This has become an inconvenience when it comes to make a policy about agro-ecosystem due to the lack of relevant data or information. In this paper, we evaluated the positive and negative effects of the agro-eco- system on a local city level, using Anyang city of Henan province as an example, and we did a dynamic analysis on the changes of agro-ecosystem value for the years of 2008-2014. We carried out the evaluations of positive and negative effects by using modified methods used in other studies, based on data collected from official Henan Yearbooks of these years. The positive effects include agricultural production, social security, soil conservation, soil carbon fixation, photosynthesis oxygen release, air cleaning, and nutrients cycling, and the negative effects include underground water consumption, greenhouse gas emission, mulching film pollution, pesticide pollution, and fertilizer loss. Results showed that the positive effects value of Anyang agro-ecosystem increased about 1 billion yuan per year, reaching 76 billion yuan in 2014. Agricultural production value is the biggest part, about 33 percent of the total positive value, while the nonproduction values, including soil carbon fixation, nutrients cycling and photosynthesis oxygen release, were the major value component, about 67% of the total positive value. Negative value showed a similar trend of increase, reaching 8.3 billion yuan in the year of 2014. The major component of negative values was underground water consumption, which was more than 72% of total negative effect value, followed by greenhouse gas emission which is about 13%. The gap between positive and negative values, namely the net value of agro-ecosystem, increased over these years, reaching about 67.6 billion in 2014. These findings provided important base data and information to decision-makings related to agricultural ecosystem management in Anyang city; it also offers a methodological framework for agro-ecosystem evaluation in other local areas in China.

Keywords:Agro-Ecosystems, Positive Effect, Negative Effect, Dynamic Evaluation, Anyang City

地方农田生态系统价值动态评价——以2008~2014年河南省安阳市为例

郜飞虎1,赵晓帅2,张晓光1,3,李士美4,王宜成1*

1青岛农业大学资源与环境学院,山东 青岛

2山东省邹平县环保局,山东 邹平

3土壤与农业可持续发展国家重点实验室,江苏 南京

4青岛农业大学园林与林学院,山东 青岛

收稿日期:2017年4月5日;录用日期:2017年5月19日;发布日期:2017年5月26日

摘 要

目前农田生态系统价值评价较多关注大尺度空间区域,地方上的动态评价很少,县市等地方政府制定农业资源管理等政策时缺少数据支撑。本文使用文献调查法获得河南省安阳市2008~2014有关数据,基于农田生态系统正负效应评价法核算了安阳市这七年农田生态系统提供的正效应和负效应价值,分析了动态变化,其中正效应包括农产品生产、社会保障、土壤保持、土壤固碳、光合作用释氧、净化大气以及维持营养物质循环,负效应包括地下水资源消耗、温室气体排放、地膜污染、农药损失和化肥流失。结果显示,安阳市农田生态系统正效应价值以约10 × 108元/年的速度增长,2014年接近760 × 108元,其中农产品生产的价值最大,约占33%,而土壤固碳等非生产性价值是主体,占总价值近67%;负效应价值基本趋势也是逐年增长,2014年约为83×108元,负效应以地下水消耗为主,占72%以上,其次是温室气体排放,约占13%。正负效应之差呈增大趋势,2014年为676 × 108元。论文为安阳市相关农业政策的制定提供了重要的数据和信息基础,为我国其他地方农田生态系统价值评价提供了参考范例。

关键词 :农田生态系统,正效应,负效应,动态评价,安阳市

Copyright © 2017 by authors and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

1. 引言

农田生态系统是依靠自然要素(土地、光照、温度、水分等)和人为投入(种子、化肥、农药、灌溉、机械等),利用农田生物与非生物环境之间以及生物种群之间的关系进行食物、纤维和其它农产品生产的半自然生态系统 [1] 。农田生态系统是陆地生态系统的重要组成部分,占全球陆地面积10%以上,自1997年Costanza等人在《自然》杂志发表文章对全球生态系统服务价值进行评价以来,农田生态系统服务价值评价也受到越来越多的研究 [2] [3] [4] 。

现在人们已经认识到,农田生态系统除了粮食蔬菜等食物生产功能,还有土壤保持、净化大气、物质循环等正效应,同时也产生水资源消耗、温室气体排放、农药污染等负效应 [5] [6] 。忽略这些正效应或负效应都会导致其价值评价失于偏颇,长期以来过分注重食物生产功能而忽略其他正负效应,导致地下水位下降、农业面源污染、物种栖息地破坏等生态环境问题就是明证。因此,在农田生态系统价值评估中全面考虑正负效应,评价其综合价值,已成为农田生态系统价值评价新的研究热点。

目前农田生态系统价值评价多集中于大尺度空间,如华北平原、黄淮海地区甚至全中国,较少关注地方如县市尺度,更少关注动态变化,而县市这些行政单元通常是地方农业政策的具体制定者,至少是实际的管理者。在这些尺度上评价农田生态系统价值,可为地方相关政策如农业资源配置、农业生态补偿等的制定提供重要的理论基础和数据支撑。这提出了建立地方农田生态系统价值评价方法体系和动态评价的需要。

安阳市是河南省重要的农产品生产基地,有“豫北粮仓”之称,农业资源管理、农业生态补偿等有关政策的制定需要相关基础数据的支持,但目前针对该市农田生态系统服务评价的研究几近空白,这种情况国内其它县市也相似。有鉴于此,笔者借鉴文献中的正负效应评价指标体系 [7] [8] ,使用文献调查法获得安阳市2008~2014年相关数据,计算了这几年的农田生态系统价值,对其动态变化进行了分析。研究结果可为该市农业资源管理等政策的制定提供基础数据和信息支持,为国内其他地方农田生态系统价值评价提供方法借鉴。

2. 材料与方法

2.1. 安阳市概况

安阳市位于河南省最北部,东经113˚37'至114˚58',北纬35˚12'至36˚22',北邻河北邯郸,西依山西长治,南接鹤壁、新乡,东连濮阳,总面积5599 km2。地势西高东低,地形复杂多样,平原、山地、丘陵、泊洼分别占总面积的53.8%、29.7%、10.8%、5.7%。气候为典型的暖温带半湿润大陆性季风气候,四季分明,历年平均气温14.1℃。极端最高气温43.2℃,极端最低气温−21.7℃。全年平均降雨量为606.1 mm。

安阳是国家历史文化名城,七朝古都。2015年全市总人口445.78万,其中常住人口401万。安阳是河南省重要的农产品生产基地,主产小麦、玉米、棉花、油菜籽、芝麻等,是黄淮平原典型的粮食高产区,农业部规划的粮食、棉花、油料优势种植区域,国家确定的全国优质小麦生产基地市。2015年粮食种植面积38.25 × 104 ha,粮食总产量227.4 × 104 t。

2.2. 研究方法

本文基于以往研究,结合黄淮平原高产农田生态系统的实际情况,从正效应和负效应两方面构建安阳市农田生态系统价值评价方法与指标体系,其中正效应包括农产品生产、社会保障、土壤保持、土壤固碳、光合作用释氧、净化大气、维持营养物质循环共七项,负效应包括地下水资源消耗、温室气体排放、地膜污染、农药污染、化肥流失共五项;使用的价值评价方法包括市场价值法、替代成本法、机会成本法、影子价格法、防护费用法共5种 [7] [8] [9] ,具体见表1

本文所使用数据除特别说明外均来自河南省统计网2009~2015年统计年鉴。

Table 1. Evaluation methods and indicators for positive and negative effects of agro-ecosystem in Anyang City

表1. 安阳市农田生态系统正负效应评价方法与指标体系

3. 结果与分析

3.1. 农田生态系统服务正效应

3.1.1. 农产品生产

安阳市主要农作物是小麦、玉米、大豆、油料、棉花、蔬菜和红薯。利用公式(1)计算出安阳市2008~2014年这些农产品生产价值(表2)。

3.1.2. 社会保障

安阳市最低生活保障标准按每人350元计,我国农业生产中隐性失业比例按46.8%计 [10] 。由公式(2)计算出安阳市农田生态系统社会保障价值(表3)。

3.1.3. 土壤保持

安阳市现实土壤侵蚀模数为1950 t∙km−2∙a−1,潜在侵蚀模数根据前人研究计算得到8448 t∙km−2∙a−1 [11] 。由公式(3)计算出安阳市的土壤保持量。本市的土壤层厚度按20 cm,土壤容重1.41 g∙cm−1,土地年收益按我国农业生产年均收益2.05 × 104元∙hm−2计,由公式(4)计算出安阳市农田减少土地废弃价值。

该市土壤库中N、P、K含量分别为0.99 g/kg、22.6 mg/kg、138 mg/kg [12] ,N、P主要通过施用化肥磷酸二铵(2400元∙t−1)获得,K通过施用氯化钾(2200元∙t−1)获得,磷酸二铵的N、P含量分别为14%和15%,氯化钾含K量为50%,由公式(5)计算出农田生态系统的保持土壤养分价值。两者之和即为安阳市农田生态系统土壤保持价值(表4)。

3.1.4. 土壤固碳

基于肖玉等人的研究 [13] ,本文使用36.73 t∙hm−2∙a−1为安阳市农田生态系统固碳速率;农田生态系统固定CO2成本按752.90元∙t−1计 [7] 。由公式(7)计算出安阳市农田生态系统土壤固碳价值(表5)。

3.1.5. 光合作用释氧

由公式(8)计算得到安阳市各类农作物的净初级生产力。工业制氧价格按1000元∙t−1计,由公式(9)计算出安阳市农田生态系统释氧价值(表6)。

3.1.6. 净化大气

农田生态系统净化大气的功能主要是对有毒气体的吸收和对粉尘的吸附,其中有毒气体我们关注最常见的SO2。农田生态系统SO2吸收能力本文用12.50 kg∙hm−2∙a−1,滞尘能力为2.1 t∙hm−2∙a−1 [14] 。SO2治理费和降尘清理费参考国家林业局发布的《森林生态系统服务功能评估规范(LY/T1721-2008)》,分别为1.20和0.15元∙kg−1。由公式(11)和(12)分别计算安阳市农田生态系统的SO2净化价值和降尘价值,求和得到净化大气总价值(表7)。

3.1.7. 营养物质循环

由公式(13),根据安阳市各类农作物产量、农作物植株内的N、P、K 含量以及各类养分的市场价格,计算得到安阳市农田生态系统维持营养物质循环的价值(表8)。

3.1.8. 正效应价值汇总

以上计算的各项正效应价值汇总于表9图1绘出了2008与2014年的对比(简明起见,其它年份未绘出)。安阳市农田生态系统服务价值最大的是农产品生产,历年大约在230~254 × 108元之间,占总价值约33%;其次是土壤固碳,在205 × 108元以上,占总价值近30%;然后是营养物质循环和光合作用释氧,都在124 × 108元以上,占总价值约18%。这四项是主要的价值组成,占总价值每年都在97%以上。七项

Table 2. Value of agricultural products (unit: 100 million yuan)

表2. 农产品生产功能价值(单位:108元)

Table 3. Value of social security (unit: 100 million yuan)

表3. 社会保障价值(单位:108元)

Table 4. Value of soil conservation (unit: 100 million yuan)

表4. 土壤保持功能价值(单位:108元)

Table 5. Value of soil carbon fixation (unit: 100 million yuan)

表5. 土壤固碳价值(单位:108元)

Table 6. Value of photosynthesis oxygen release (unit: 100 million yuan)

表6. 光合作用释氧价值(单位:108元)

Table 7. Value of air cleaning (unit: 100 million yuan)

表7. 净化大气价值(单位:108元)

Table 8. Value of nutrients cycling (unit: 100 million yuan)

表8. 维持营养物质循环功能价值(单位:108元)

Table 9. Positive effects of agro-ecosystem in Anyang 2008-2014 (unit: 100 million yuan)

表9. 安阳市农田生态系统2008~2014年正效应价值(单位:108元)

括号中数字是该项在当年总价值中的百分比。Numbers in brackets are value percentage in the year.

Figure 1. Comparison of positive effects of agro-ecosystem in Anyang between 2008 and 2014

图1. 安阳市农田生态系统2008与2014年各项正效应价值比较

价值除社会保障外,其它六项都逐年递增,这主要归因于逐年增长的农作物种植面积和总产量。社会保障价值受到保障人数等复杂社会因素的影响有所起伏,约3 × 108元,占总价值不到0.5%。另外值得注意的是,土壤固碳、营养物质循环等这些非生产性价值约占总价值的67%,且近三年来有增加趋势,显示安阳市农田生态系统在生产农产品的同时,更对维护当地自然生态系统平衡起到巨大作用。

3.2. 农田生态系统负效应

3.2.1. 地下水资源消耗

安阳市地处黄淮平原,年降水量较少,境内无大河或水库,农田灌溉主要依靠地下水。单位面积灌溉成本以黄淮海地区平均灌溉成本2000元∙hm−2计 [15] ,由河南省统计局得到安阳市历年农田灌溉面积,由公式(14)得到安阳市地下水资源消耗费用(表10)。

3.2.2. 温室气体排放

北方旱地排放的温室气体主要是CO2和N2O。基于前人研究,安阳市农田生态系统CO2排放量取0.57 t∙hm−2∙a−1,N2O排放量折算为CO2是1.37 t∙hm−2∙a−1,CO2排放成本取752.90元∙t−1 [16] 。由公式(15)计算出安阳市2008~2014年温室气体排放成本(表11)。

3.2.3. 地膜污染

地膜污染成本主要表现为影响土壤物理结构和生物活动从而导致农作物减产。根据已有研究,地膜残留造成的农作物减产损失率为0.5% [7] 。由公式(16)计算出安阳市地膜污染成本(表12)。

3.2.4. 农药损失

根据农业部最新数据,我国农药利用率为36.6% [17] ,按此数值计算安阳市农药流失,农药价格按农药市场信息河南省近期农药平均价格37.55元∙kg−1,由公式(18)计算安阳市农药流失成本。据研究,滥用农药导致农作物减产率约1% [7] ,由公式(19)计算粮食减产造成的经济损失。两者之和为农药损失成本(表13)。

3.2.5. 化肥流失

按农业部最新数据,我国化肥利用率为35.2% [17] ,化肥平均价格按2300 元∙t−1计,由河南省统计网得到安阳市2008~2014年化肥使用量,由公式(20)计算出安阳市化肥流失成本(表14)。

3.2.6. 负效应成本与构成

以上计算的各项负效应成本总结于表15图2绘出了2008与2014年成本对比。可以看出,地下水消耗是负效应最主要的部分,每年都在60 × 108元左右,占总成本72%以上。安阳市地处华北平原,降水量较小,也无大面积的地上水体,灌溉只能以地下水为主,所以安阳市要确保农业可持续发展,一个主要的工作重心是解决灌溉用水问题,一方面保护地下水资源,另一方面要改进灌溉技术,减少水资源消耗。温室气体排放居第二位,每年11 × 108元左右,占总成本约13%,且呈增长趋势,主要是由于大量施用化肥特别是氮肥所致。另外地膜和农药的大量施用也导致这两项成本呈逐年增长趋势。

3.3. 安阳市农田生态系统正负效应动态分析

图3绘出了安阳市农田生态系统2008~2014年正负效应及其差值的变化。正效应明显大于负效应,每年大约都是负效应的9倍左右。正负效应之差(即综合效应)逐年增长,由2008年的624 × 108元增加到2014年的676 × 108元,这主要归因于农田生态系统生物量的增加以及由此带来的农产品生产、土壤固碳、营养物质循环等价值的增长,而负效应部分则相对增加较小。但前面分析已发现,负效应中70%以上是

Table 10. Cost of underground water consumption (unit: 100 million yuan)

表10. 地下水资源消耗费用(单位:108元)

Table 11. Cost of greenhouse gases release (unit: 100 million yuan)

表11. 温室气体排放成本(单位:108元)

Table 12. Cost of mulching plastic film pollution (unit: 100 million yuan)

表12. 地膜污染成本(单位:108元)

Table 13. Cost of pesticide loss (unit: 100 million yuan)

表13. 农药损失成本(单位:108元)

Table 14. Cost of fertilizer loss (unit: 100 million yuan)

表14. 化肥流失成本(单位:108元)

Table 15. Negative effects of agro-ecosystem in Anyang 2008-2014 (unit: 100 million yuan)

表15. 安阳市农田生态系统2008~2014年负效应成本(单位:108元)

括号中数字是该项在当年总成本中的百分比。Numbers in brackets are cost percentage in the year.

Figure 2. Comparison of negative effects of agro-ecosystem in Anyang, 2008-2014

图2. 安阳市农田生态系统2008与2014年各项负效应成本比较

Figure 3. Dynamics of positive and negative effects and their differences (i.e. net effect)

图3. 正负效应及其差值(综合效应)的历年变化

地下水消耗,而黄淮海地区地下水正越来越成为稀缺资源,地下水位正以每年1~3 m的速度下降 [15] 。由此可以预计,农业灌溉使用地下水的成本有可能迅速升高,导致正负效应的差距迅速缩小。这再一次突显了该地区地下水资源保护和节水灌溉的重要性。

4. 总结和讨论

4.1. 总结

本文核算了安阳市农田生态系统正负效应的价值,主要发现总结如下:

1) 正效应价值以大约10 × 108元/年的速度增长,2014年接近760 × 108元,其中农产品生产价值约占33%,其次分别是土壤固碳(28%)、营养物质循环(18%)、光合作用释氧(18%)、土壤保持(1.7%)、社会保障(0.4%)和净化大气(0.3%),这些非生产性价值占总价值约67%。当地农田生态系统在提供农产品价值的同时,为自然生态系统平衡提供更大的价值。

2) 负效应价值基本趋势也是逐年增加,2014年是83 × 108元多,其中主要是地下水消耗,占总成本72%以上,其次是温室气体排放(13.5%)、化肥流失(8.5%)、农药污染(4.0%)和地膜污染(1.5%)。安阳市农业要实现可持续发展,必须确保足够的灌溉用水,重点可从两方面开展工作,一是保护地下水资源,二是采用节约型灌溉,减少水资源消耗。

3) 正效应与负效应之差即农田生态系统净效益逐年增长,2014年是676 × 108元,主要归因于逐年增长的生物量以及由此带来的农产品生产、土壤固碳、营养物质循环等价值的增长。

4.2. 讨论

从正、负效应两方面考察农田生态系统价值是近年来生态农业领域的研究热点之一,本文选取华北平原一个典型的粮食高产区--河南省安阳市进行了案例研究,为安阳市以及国内其他县市农业资源管理、农业生态补偿等相关政策的制定提供了基础数据和信息参考。农田生态系统价值评价也可从其它角度展开,如从生命周期评价角度评价农田生态系统的环境影响 [18] ,用能值分析法评估农田生态系统服务价值 [19] ,基于碳足迹的农田生态系统评价 [20] 等,多种方法的存在为方法比较和择优提供了可能,这值得在将来研究中进一步探讨。

本文正负效应各项都比较容易判别和计算,但并不是农田生态系统价值的全部。选取了小麦、玉米等7种主要农作物而忽略其它农产品,数据获取和处理过程不可避免会出现误差,动态评价忽略物价变化,这些都会给评价结果造成一定影响,但短期内不会太显著,如果对长期(比如几十年)的生态系统服务价值进行比较则应当考虑物价变化的影响,必要时包括其它农产品也会提高评估精度。

农田生态系统还有其它如文化功能、审美价值等正效应,也有物种单一、栖息地破碎、地下水污染等负效应。其中有些能够定量评价,如用人力资本法评价农业生产区地下水硝酸盐污染造成的损失 [21] ,但通常难以定量化或者数据获取比较困难,本文没有包括在计算之中,随着方法的改进和数据的完善,有待以后在农田生态系统价值评价中包括这些方面。

基金项目

土壤与农业国家重点实验室开放基金项目(Y20160007);山东省高校科技发展计划项目(J14LF04)。

文章引用

郜飞虎,赵晓帅,张晓光,李士美,王宜成. 地方农田生态系统价值动态评价——以2008~2014年河南省安阳市为例
A Dynamic Evaluation of Agro-Ecosystem at Local Levels—Using Anyang City, Henan Province 2008-2014 as an Example[J]. 世界生态学, 2017, 06(02): 49-60. http://dx.doi.org/10.12677/IJE.2017.62006

参考文献 (References)

  1. 1. 谢高地, 肖玉. 农田生态系统服务及其价值的研究进展[J]. 中国生态农业学报, 2013, 21(6): 645-651.

  2. 2. Costanza, R., d’Arge, R., Groot, R. D., et al. (1997) The Value of the World’s Ecosystem Services and Natural Capital. Nature, 387, 253-260. https://doi.org/10.1038/387253a0

  3. 3. Kinzig, A.P., Perrings, C., Chapin III, F.S., et al. (2011) Paying for Ecosystem Services—Promise and Peril. Science, 334, 603-604. https://doi.org/10.1126/science.1210297

  4. 4. Daily, G.C. and Matson, P.A. (2008) Ecosystem Services: From Theory to Implementation. Proceedings of the National Academy of Sciences of the United States of America, 105, 9455-9456. https://doi.org/10.1073/pnas.0804960105

  5. 5. Zhang, W., Ricketts, T.H., Kremen, C., et al. (2007) Ecosystem Services and Dis-Services to Agriculture. Ecological Economics, 64, 253-260.

  6. 6. Foley, J.A., Defries, R., Asner, G.P., et al. (2005) Global Consequences of Land Use. Science, 309, 570-574. https://doi.org/10.1126/science.1111772

  7. 7. 元媛, 刘金铜, 靳占忠. 栾城县农田生态系统服务功能正负效应综合评价[J]. 生态学杂志, 2011, 30(12): 2809- 2814.

  8. 8. 杨志新. 北京郊区农田生态系统正负效应价值的综合评价研究[D]. 北京: 中国农业大学, 2006: 15-20.

  9. 9. 刘鸣达, 黄晓姗, 张玉龙, 等. 农田生态系统服务功能研究进展[J]. 生态环境, 2008, 17(2): 834-838.

  10. 10. 孙新章, 周海林, 谢高地. 中国农田生态系统的服务功能及其经济价值[J]. 中国人口∙资源与环境, 2007, 17(4): 55-60.

  11. 11. 陈宇, 李子轩, 齐建怀, 等. 漳卫河流域山区土壤侵蚀强度变化特征分析[J]. 海河水利, 2013(6): 20-22.

  12. 12. 张晓勇. 安阳市土壤养分现状分析[J]. 河南农业, 2014(3): 24-25+27.

  13. 13. 肖玉, 谢高地, 安凯, 等. 华北平原小麦-玉米农田生态系统服务评价[J]. 中国生态农业学报, 2011, 19(2): 429- 435.

  14. 14. 田耀武, 王建东, 张世平. 郑州西南绕城高速公路不同时期生态系统服务价值研究[J]. 中南公路工程, 2006, 31(3): 81-85.

  15. 15. 隋鹏, 陈素英, 陈源泉, 等. 黄淮海地下水亏水区农业生产用水生态代价评估——以河北省栾城县为例[J]. 中国生态农业学报, 2007, 15(5): 178-181.

  16. 16. 白杨, 欧阳志云, 郑华, 等. 海河流域农田生态系统环境损益分析[J]. 应用生态学报, 2010, 21(11): 2938-2945.

  17. 17. 李慧. 我国化肥利用率提高2.2个百分点[EB/OL]. 光明日报. http://epaper.gmw.cn/gmrb/html/2015-12/03/nw.D110000gmrb_20151203_3-08.htm, 2015-12-03.

  18. 18. 梁龙, 陈源泉, 高旺盛, 等. 华北平原冬小麦-夏玉米种植系统生命周期环境影响评价[J]. 农业环境科学学报, 2009, 28(8): 1773-1776.

  19. 19. 马凤娇, 刘金铜. 基于能值分析的农田生态系统服务评估——以河北省栾城县为例[J]. 资源科学, 2014, 36(9): 1949-1957.

  20. 20. 段华平, 张悦, 赵建波, 等. 中国农田生态系统的碳足迹分析[J]. 水土保持学报, 2011, 25(5): 203-208.

  21. 21. 刘光栋, 吴文良, 刘仲兰, 等. 华北农业高产粮区地下水面源污染特征及环境影响研究——以山东省桓台县为例[J]. 中国生态农业学报, 2005, 13(2): 125-129.

期刊菜单