基于聚类分析的甘肃省静态与动态农业生态服务价值时空分布和变化趋势

刘阳 张威 章高森 张昺林 姜右锦 张松林

摘要：以生態服务功能价值当量因子法估算了2008—2017年甘肃省农业生态服务静态价值构成和不同类型生态服务功能静态价值，在此基础上，考虑到农业生态服务的空间动态性，根据聚类分层分析不同服务功能的权重比赋加优先指数估算了各项生态服务功能的动态价值构成，并对其时空变化与时空分布进行了探究。结果表明，2008—2017年，甘肃省农业生态服务功能动态价值逐年增长，不同功能动态价值贡献差异较大，食物生产与其他服务功能动态价值相关性最低;在此期间，甘肃省农业生态服务静态总价值达21 873.4亿元，而动态总价值为18 345.8亿元，相比于静态价值，各项动态功能价值均降低。对其动态价值变化趋势进行冷热点预测发现，甘肃省农业生态环境变化不稳定，各服务功能构成不均衡，且弹性较差，恢复力较弱。甘肃省亟需开展农业生态服务功能建设，大力发展农业经济，积极建设生态农业，以保持农业生态平衡，推动农业生态系统可持续发展和提高农业生态系统服务功能价值。

关键词：农业生态服务价值;静态;动态;生态平衡;甘肃省

中图分类号：F062.2 文献标识码：A

文章编号：0439-8114（2020）07-0099-10

DOI：10.14088/j.cnki.issn0439-8114.2020.07.021

Abstract： The composition of agricultural ecological services static value and different types of ecological service functions static value in Gansu province from 2008 to 2017 was estimated by the ecological service function value equivalent factor method. In consideration of the spatial dynamics character of agricultural ecological services，the dynamic value composition of various ecological services was estimated by the weight ratio assignment of different service functions based on clustering and stratification and priority index. And the temporal-spatial distribution of the dynamic value of agro-ecological services in Gansu province was studied. The results showed that the dynamic value of agricultural ecological service functions in Gansu province increased by years from 2008 to 2017. Different dynamic value contributions varied greatly. Grain conservation and other service functions had less correlation with the dynamic value. Total static value of agricultural ecological services in Gansu province from 2008 to 2017 reached 2 187.34 billion yuan. The total dynamic value was 18 34.58 billion yuan. Compared with static value，the value of each dynamic function reduced. The results of variation trend of dynamic value predicted with hots-colds pots found that changes in the agro-ecological environment in Gansu province were instability. The composition of service functions was unbalanced，as well as the elasticity was poor and restorative force was weak. To sum up，Gansu province is in urgent need of the construction of agricultural ecological service functions vigorously developing agricultural economy，actively constructing ecological agriculture，in order to maintain the balance of agricultural ecology，promote the sustainable development of agricultural ecosystems，and increase the value of agricultural ecosystem services.

Key words：value of agricultural ecosystem services; static state; dynamic state; ecological balance; Gansu province

农业生态系统是陆地生态系统不可或缺的一部分。随着全球城市化进程的发展，农业生态系统已从半自然生态系统逐渐演变为以提供农产品为主的集约化生产系统[1]。这是由于对农业生态系统自然生产力与经济生产力之间平衡关系的认识不够深入[2]，人类在追求农业经济效益的过程中破坏了农业生态系统结构和功能的完整性[3]。农业资源分配和利用不均可能导致农业生态系统能量流动不畅、物质循环受阻和信息传递滞后，最终导致农业生态系统瘫痪，使原有的社会性、高产性和波动性变差[4]。因此，深入研究农业生态服务功能及其价值可以为农业生态系统可持续发展和生态服务功能完备、均衡性保持提供科学依据。

20世纪90年代以来，随着生态系统各项服务功能评价方法和评估策略的不断完善[5-7]，国内外学者通过对农业生态服务静态价值的评估和预测有效保护了农业生态系统[8-10];且大多数学者侧重于生态系统服务的土地利用类型和区域性影响研究[11]，尤其是草地和森林生态服务价值[12]。然而，由于农业生态服务价值评估的数据欠缺，不具有完整性、连续性和可变性[13]，导致国内外农业生态服务研究多局限于单一年份的静态式价值变化的研究[14]，对不同区域的随时间、生物和环境因子等动态因子影响下的动态价值变化的研究还极为少见[15]，甘肃省等西北干旱、半干旱地区亦然。聚类分层分析方法在经济学、生物学、地理学和生态学等各方面的分析应用较为成熟[16]。在生物学方面，通过对动植物及其基因进行聚类分析，探究其种群结构间的关系[17];在地理学方面，主要对不同地理因子进行聚类分析，探究不同地域地理因子间相互影响的相似性等[18]。本研究对农业生态服务功能进行聚类分层分析，旨在聚类10年内各项服务功能静态价值的相似性，从而对其静态价值贡献率进行分层，找出主要服务价值和负反馈服务价值。

甘肃省是“丝绸之路”中国西北地区的交通枢纽，农业以传统的二元种植结构为主。由于过度追求农业生态系统的经济性，忽视农业生态服务功能的完备性和均衡性，影响了甘肃省农业生态系统的可持续发展。本研究在较为成熟的静态价值研究基础上，对各项服务价值进行优先指数的构建[19]，由于生态系统服务功能所产生的价值随着外部价值形态和内部价值结构的变化而变化[20]，故在构建优先指数的基础上对农业生态系统的各项服务功能进行聚类分层加权。结合农业生态服务动态价值变化，应用地统计学进行空间预测以及冷热点空间变化特征的探究，发现甘肃省农业生态服务在各项服务功能上表现出动态服务价值的差异性，可以更加深入和靶向性地针对某一生态服务功能进行调整和保护，从而更加准确地保持和保护农业生态系统的稳定性，同时提高了农业生态服务的完整性和价值创造力。因此，对甘肃省农业生态服务价值的时空变化规律进行研究可以为甘肃省农业生态服务功能评估和生态农业可持续发展提供科学依据。

1 研究区域与方法

1.1 研究区域概况

甘肃省位于中国西北地区，拥有除海洋以外的所有地貌;总面积约42.59万km2（截至2018年底），包括12个地级市、2个自治州;人口达到209.95万（截至2019年5月）处于北纬32°31′ —42°57′、东经92°13′ —108°46′ ;年平均气温为2～10 ℃，年平均降水量为250～420 mm;是世界“丝绸之路经济带”的重要组成部分，也是“丝绸之路”中国西北地区的交通枢纽。

1.2 研究方法

1.2.1 数据来源 2008—2017年，甘肃省不同农业土地利用类型面积数据来自相应年份的《中国统计年鉴》。甘肃省3种主要农作物小麦、玉米和马铃薯及其播种面积数据来自2008—2017年的《甘肃省统计年鉴》《中国统计年鉴》和《甘肃农村年鉴》。根据甘肃省各市（县）物价局监测数据、《中国物价年鉴》和居民分类消费价格指数，计算小麦和玉米平均价格。根据《中国物价年鉴》和居民分類消费价格指数计算马铃薯平均价格。

1.2.2 研究方法 参照Costanza等[5]的生态系统服务估算价值模型，谢高地等[21]将中国的生态系统服务价值分为9类，并在对中国200多位生态学者问卷调查的基础上制订了中国的生态系统服务价值当量因子表（表1）。

生态系统服务价值当量因子1是指达到全国平均产量的农田每公顷每年自然粮食产量的经济价值，主要表现为生态系统产生的生态服务相对贡献大小的潜力;其他当量因子值是指生态系统产生的该生态服务相对于农田食物生产服务的贡献大小。

1.2.3 分析方法

1）单位面积农田经济作物生产服务价值的计算。由谢高地等[21]对生态服务功能价值当量因子的定义可推知，1个生态服务价值当量因子的经济价值量等于当年全国平均粮食单产价值的1/7。由统计年鉴可知，甘肃省农田主要种植的粮食作物为马铃薯、玉米和小麦，故采用三者计算单位农田粮食作物生产服务的价值。参照谈存峰等[22]的研究结果，选用的计算公式为：

式（1）中，GESV为单位面积农田经济作物生产服务功能的价值（元/hm2）;Si为第i种粮食作物的面积（hm2）;Pi为第i种粮食作物的均价（元/t）;Ai为第i种粮食作物的单产（t/hm2）;M为n种粮食作物的总面积（hm2）。

2008—2017年甘肃省单位面积农田主要粮食作物生产服务功能价值分别为947.09、940.07、 1 015.00、1 128.62、1 313.86、1 342.28、1 377.00、 1 428.40、1 467.34、1 502.92元/hm2。除2009年略有下降外，甘肃省单位面积农田主要粮食作物生产服务功能价值随着时间延续逐年增加，两者之间呈极显著线性正相关（R2为0.966 8，P<0.01）。

2）单位面积农田经济作物生产服务功能价值的修订。中国学者认为中国生态服务价值的计算应该在Costanza等[5]学者提出的生态服务经济价值的计算原理和方法基礎上进行修正。这是因为中国生态学和经济学界认为Costanza等的研究方法能够很好地体现欧美发达国家的生态服务功能价值，对于中国这样的发展中国家，该方法使得生态服务功能价值被高估，而农业耕地的生态服务价值被低估。所以，谢高地等[21]学者在Costanz等方法的基础上根据对中国几百位生态学者的调查报告对不同生态服务功能单价进行了一定的修正。本研究根据研究范围对修正模型进行矫正后的公式为：

式（2）中，GESVi为单位面积农田经济作物生产服务功能的修订价值（元/hm2）;Tb为不同农业土地类型生态系统总生物量;B为甘肃省农业生态单位面积生物量。

3）农业生态服务功能静态价值（AESV0）计算。

2008—2017年甘肃省农业生态服务功能的单价详见表2。

4）农业生态服务功能优先指数的构建。同一生态系统服务功能在不同类型的生态系统中所占的权重不同，即同样的生态服务功能在不同类型生态系统中所发挥的价值不同;同时，随着研究区域的时空变化，其生态系统的生态服务功能价值随之变化。因此，王云等[19]提出用生态服务价值的多年变化率与其均值的比值作为优先指数进行价值计算。为了更具体地看出农业生态服务功能价值动态变化与其不同价值构成的趋向性，本研究在Govender等[16]研究的基础上，采用分层聚类的方法，通过对每项农业生态服务功能的静态价值年平均变化率来表示对该项服务需要保护的迫切程度，再利用分层聚类的方法对服务价值损失越多的生态系统类型增加更高的聚类权重比，即优先级别越高。公式为：

式（4）中，[αi]为第i种生态服务功能价值的优先指数;[vGESVi]为第i种生态服务价值的10年平均变化率;[vaGESV]为不同生态服务价值的10年平均变化率的平均值;[Ci]为第i种生态服务价值的10年聚类分析的占比权重。

5）农业生态服务功能动态价值计算。在充分考虑不同类型生态服务功能的动态价值对甘肃省农业生态系统保护与恢复导向性差异的基础上，参照王云等[19]、谈存峰等[22]的研究结果，确定其生态服务功能价值（AESV）计算模型，为：

式（5）中，A为农业用地面积（hm2）;Ei为第i种生态服务功能价值的当量因子;GESVi为修订后单位面积农田经济作物生产服务功能价值;αi为第i种生态服务功能价值的优先指数。

采用SPSS 24.0统计软件处理、分析试验数据，采用ArcGIS 10.5和Origin Pro 9.0绘图。

2 结果与分析

2.1 2008—2017年甘肃省农业生态服务功能静态价值的时空分布

10年来，甘肃省农业生态服务静态总价值达到21 873.4亿元。2008—2017年，甘肃省农业生态服务静态总价值逐年增加，分别为1 636.6亿、1 640.8亿、1 776.4亿、1 961.3亿、2 255.5亿、2 310.7亿、2 352.7亿、2 439.0亿、2 645.9亿、2 854.3亿（图1）。从图1可以看出，2009—2012年甘肃省农业生态服务静态总价值增长较快，2012年同比增长率为近10年内的最大值（15%）。

2008—2017年，甘肃省大部分地区农业生态服务静态总价值呈平稳增长的趋势，少部分地区呈缓慢增长甚至下降趋势，且不同城市间呈碎片化和不规律的空间分布格局（图2）。对于同一区域来说，10年来静态总价值一直在增加，由于局限于行政区划农业生态系统服务价值构成不同，所以不同区域间的静态价值差异性一直存在且比较明显。2011—2012年，甘肃省东部和北部（酒泉）地区的农业生态服务静态总价值明显增加，其余年份的静态总价值空间分布较不明显。

通过静态总价值（图2）分析发现，甘肃省西北部、中部和南部农业生态服务静态总价值增长迅速，中部温带大陆性气候的城市经济区至西南部（甘南藏族自治州）大陆性气候中间的农牧过渡地区[23]呈先增长后减小的趋势，同比静态总价值变幅为

-4.1%～10.5%，说明城市经济区与农牧过渡区农业生态服务静态总价值的差距较大，经济-农牧过渡区农业生态服务功能的变化不稳定。虽然中部经济带区的农业生态服务静态总价值逐年增长，但总体水平低于西南部农牧区与经济-农牧过渡区。

由于甘肃省特殊的东南-西北走向的狭长地形，2008—2017年，甘肃省农业生态服务静态总价值的时空分布总体表现为西北高、东南低，这可能是因为不同区域的农业生态系统的构成不同导致服务功能的价值贡献不同。甘肃省西部毗邻祁连山山脉，丰富的自然资源条件[24-26]为其农业服务功能发挥作用提供了基础;东部紧邻阿拉善荒漠，其静态总价值虽整体低于西部，但是呈稳定增长的趋势，说明不同区域农业生态服务功能静态总价值的占比和贡献率可变。

2.2 不同农业生态服务功能静态与动态价值分析

2.2.1 不同农业生态服务功能静态价值分析 生态系统服务功能是否正常发挥事关农业生态系统的平衡、安全与可持续发展。2008—2017年，甘肃省不同农业生态服务功能的静态价值均具有增长趋势，但各自对甘肃省农业生态服务功能静态10年总价值的均值贡献占比不同，依次为土壤形成与保护（18.85%）>水源涵养（16.25%）>废物处理（14.47%）>生物多样性与保护（14.19%）>气体调节（13.08%）>气候调节（11.97%）>原材料生产（7.92%）>食物生产（4.80%）>娱乐文化（2.89%），这与谈存峰等[22]对2008年甘肃省兰州市农业生态服务静态价值的研究结果类似。

从图3可以看出，甘肃省农业生态系统土壤形成与保护功能静态价值的同比增长率与农业生态服务功能静态总价值的同比增长率相似，2012年达到峰值，2012—2014年急劇下降。水源涵养功能与废物处理的静态价值同比增长率变化相类似，整体呈平缓增长趋势。2008—2017年，气体调节、气候调节和原材料生产的同比增长率变化相似。气体调节和气候调节对生命支持系统分别具有辅助和支撑作用，对生物种群的生存与发展有较大影响[27]。

近10年来生物多样性与保护对甘肃省农业生态服务年均静态总价值的贡献为310.41亿元，与谢高地等[21]对青藏高原每年陆地生态系统中生物多样性与保护对生态系统服务价值的贡献（1 494.5亿元）相比少79.23%，说明甘肃省农业生态系统中的生物多样性少、食物链简单、生物对农业生态系统服务功能的静态价值贡献较低，因此，应保持农业生态系统中的生物平衡[28]，以促进甘肃省农业生物多样性的发展。娱乐文化功能的静态价值占比最低，同比增长率波动较大，整体呈下降趋势（图3）。

从图4可以看出，不同农业生态服务功能近10年静态价值的变化趋势仅与原材料生产呈显著相关（P<0.05），说明在静态价值转换的模型计算中只有原材料生产的变化对农业生态功能具有显著反映，而其他服务功能对农业生态功能的贡献不明显。这也可能是导致近10年静态价值变化平稳的最主要原因。但是对近10年不同农业服务功能静态价值进行热图分析（图4a）发现，不同服务功能静态价值时间变化层次明显，并且同一服务功能静态价值年际变化也十分明显。对于不同的生态服务功能需要动态分析的方法，要充分考虑到农业生态区各项服务功能时空分布上的变化特征[20]。然而，对静态价值进行聚类分析（图4b）发现，不同服务功能之间变化呈现出其静态价值随时间变化不明显。这更进一步说明农业生态服务价值不仅需要考虑静态变化，更要考虑到其他生态环境因子的影响。

2.2.2 不同农业生态服务功能的动态价值评估 采用聚类分层分析法将近10年来不同农业生态服务功能静态价值进行矩阵二维化，可以明显看出，虽然同一年份不同服务功能之间变化不明显，但从时间变化可以看出2008—2011、2012—2014、2015—2017这3个不同时间段静态价值热图表达明显不同，聚类系数同比变化较大，为28.6%～73.7%（图4b）。对不同农业生态服务功能而言，娱乐文化和原材料生产、食物生产以及其他6项服务功能在这3个年份变化阶段聚类分层明显不同，对农业生态服务功能价值反映各占权重的25%。气候调节、气体调节和生物多样性与保护归为一类，这与生态系统中生物多样性与气候调节和气体调节有直接关系[27]。由图4b的服务功能聚类系数，依据式4计算得出不同农业生态服务功能优先指数（表3）[19]。

由优先指数计算结果（表3）可知，土壤形成与保护和食物生产优先指数较高，娱乐文化与水源涵养优先指数较低，其余变化范围较小，说明对于农业生态服务功能价值贡献而言，土壤形成与保护和食物生产功能价值对其贡献较高。

由不同服务功能的优先指数，利用式（5）估算出甘肃省2008—2017年农业生态服务动态总价值为18 345.8亿元，低于静态总价值。通过估算发现，不同农业生态服务功能的动态价值与静态价值相比，不同年份、不同服务功能价值变化规律不同。分析发现，食物生产相较于其他服务功能对农业生态服务功能价值有直接影响。甘肃省农业生态服务动态价值对水源涵养功能有较强的依赖性，可能与近年来甘肃省对本省水资源现状的认识比较到位，采取了积极的抗旱节水措施，对水资源进行了有效地利用与保护有关[27]。

不同生态服务功能动态价值相关性分析表明，食物生产与其他服务功能的动态价值相关性较低（图5a），结合其动态价值聚类分析（图5b），说明其对甘肃省农业生态服务功能的影响较大，且具有单一性。这可能与甘肃省近年来主要粮食作物产量不稳定、经济价值浮动较大[29]有关，所以要积极均衡各类型粮食作物的种植，粮食作物要多样化，作物种植技术要科学化，产量要保持稳定化。气候调节与气体调节的相关性系数较高，为0.996 8（图5a），表明这两种服务功能对区域性气候变化以及空气质量有一定的影响，导致的直接结果是影响到农业作物的生长和产量，间接地对大气质量有一定的影响，有可能对农业生态服务功能产生负作用[2]，也可能对生态环境产生一定破坏。通过对不同农业生态服务功能静态和动态价值进行分析发现，虽然不同农业生态服务功能价值逐年上升，但是不同服务价值对总价值产生的贡献价值不同，甚至有些服务价值在一定年份产生负影响，说明甘肃省各类型生态服务价值时空分布不稳定。通过聚类分析不同服务动态价值对总价值的贡献值，可以更好地调控贡献值较低的服务功能，发掘其潜在的服务价值，更好地利用在农业生态环境，维持农业生态环境的稳定。

2.3 甘肃省农业生态服务功能的动态价值时空分布与评估

通过图5b甘肃省农业生态服务功能静态价值进行聚类发现，2008、2011、2014、2017年4个年份静态价值热图表达明显不同，说明这4个年份不同服务价值功能变化较大。所以，选择2008、2011、2014、2017年的生态服务动态价值进行克里金插值分布预测，发现甘肃省西部和东部农业生态服务功能变化较大。从图6可以看出，甘肃省西部酒泉地区农业生态服务动态总价值变化与分布不稳定，近年来酒泉地区逐渐呈聚集增长和发散型降低;甘肃省东部和东南部地区也呈现出相类似的斑块状变化，说明甘肃省近年来农业生态服务功能时空分布不稳定，发展不均衡。这可能与甘肃省地形狭长、气候多变、农业用地分布不均匀、年际变化较大有一定的关系[30]。同时，在促进城市化发展的基础上也要提高农业用地的利用率，减少建筑用地的占有率[31]，不仅可以增加食物生产，也可以加大不同农业生态服务功能对农业生态系统的补偿作用，对农业生态环境的有效利用有很好的预测和指示作用。

图7对时间变化显著的农业生态服务动态总价值进行了差值贡献率计算，即[β=ΔAESVΔTAESV]，其中，ΔAESV表示年度间农业生态服务动态总价值增减量，ΔTAESV表示总农业生态服务动态价值增减量，且β>0表示该区域农业生态服务功能正增长，β<0表示该区域农业生态服务功能负增长。从图7可以看出，甘肃省近年来农业生态服务总价值不断增长，但是各区域增长速率不同，2008—2011年农业生态服务动态价值热点比2011—2014年多，说明2008—2011年农业生态服务动态总价值要高;而2014—2017年动态价值整体热点增多，大多数β>0，说明甘肃省近年来农业生态服务动态总价值整体提高。但是，纵观2008—2017年冷热点变化，甘肃省农业生态服务动态价值发展依旧不均衡，整体生态服务动态总价值较低。

由甘肃省农业生态服务动态总价值冷热点变化分析结合各服务功能静态价值年际变化和动态价值变化分析可知，导致甘肃省近年来生态服务动态价值变化差距大、不规律的原因可能有：①酒泉等西北部地区大部分地区均为戈壁、沙漠，不同生态服务功能动态价值贡献率低，人口稀疏，农业生态环境的价值创造力偏低，利用不合理，导致了酒泉等西北部地区斑块化的不均衡发展模式[32，33]。②兰州、天水等中部到东部地区农业生态服务动态价值逐年增长，发展均衡，与这些地区为保护农业生态环境、近年来大量减少建筑用地对农业用地的占用、积极提高土地的利用率有关;东南部甘南地区整体农业生态服务动态价值水平较高，但发展不均衡，可能是该地区农业生态环境破坏严重、生态弹性较差导致各生态服务功能发展差异性较大。③食物生产对甘肃省农业生态服务动态价值的影响较大，要增加粮食用地面积，提高各区域粮食作物种植的均衡性，丰富农业用地的利用率，间接性地对农业生态系统中气候调节、气体调节、水源涵养和生物多样性与保护的价值起到增强作用。④娱乐文化的动态和静态价值均较低，像戈壁、沙漠这样的区域，农作物产量差，耕地面积少，应该根据其独特的地理风貌，发挥娱乐文化价值，因地制宜大力发展旅游业。现阶段规模化观光农业的价值大增，如青海门源的百里油菜花田。因此，应大力发展甘肃省的休闲观光旅游农业[34，35]，增加农业生态旅游的含金量[36]。

3 小结与讨论

农业生态系统为人类的生存与发展提供了物质和能量，是人类赖以生存不可或缺的生态系统之一。农业生态系统服務功能如若退化，那么人类的物质来源和能量来源都将受到影响，不仅会影响人类的健康和安全，甚至还可能会影响到人类的生存。由于生态系统服务价值具有空间性而并非静态式[20]，农业生态系统服务价值亦然，即农业生态系统各服务功能价值的变化会受到农业生态系统外部生态环境因子和内部自身变化的影响。本研究利用当量因子法对甘肃省2008—2017年农业生态服务静态价值进行估算，同时考虑到农业生态系统变化的动态性，对其各项服务功能的静态价值通过聚类分层的方法赋加静态价值权重和优先指数，从而对静态价值进行重构，估算其动态价值。甘肃省农业生态服务静态总价值为21 873.4亿元，2008—2017年呈逐年增长趋势，依次为1 636.6亿、1 640.8亿、1 776.4亿、1 961.3亿、2 255.5亿、2 310.7亿、2 352.7亿、 2 439.0亿、2 645.9亿、2 854.3亿元。2012年，甘肃省农业生态服务功能静态总价值的同比增长率最高，为15%。甘肃省2008—2017年农业生态总价值为18 345.8亿元，低于静态总价值。

不同农业生态服务静态价值构成中，食物生产和原材料生产有显著年际变化，这可能与农业生态系统以农田生态系统为主，而农田生态系统是食物生产和原材料生产的源头[37，3]。研究发现甘肃省服务功能静态总价值变化以3年为一个年度发生显著变化。不同农业生态服务静态价值通过动态转换，可以更加明显地区分不同服务功能对农业生态系统服务贡献的重要性，也反映了不同农业生态服务功能在空间上的差异性。通过不同服务功能动态价值分析发现：水源涵养、气候调节、气体调节和生物多样性与保护的动态价值对甘肃省农业生态服务贡献较高。但是，由于甘肃省自身自然条件的因素制约，水源涵养和生物多样性与保护都不是自身的优势服务功能，这说明甘肃省这2个服务功能的直接价值利用较好，并且潜在服务价值巨大。而食物生产的动态服务价值虽然变化浮动较大，整体呈上升趋势，但是其与甘肃省农业生态服务动态总价值的相关性较低，说明甘肃省对农田生态系统的利用不完善，农田生态系统与森林、草地、水体等农业生态系统在空间配置上不合理;农业生态补偿措施不完善[38];应该大力发展生态农业旅游[36]，使得农田生态系统多功能化，激活耕地的农业生态系统价值，提高粮食生产所产生的农业生态服务价值。故对于甘肃省农业生态系统而言，促进甘肃省食物生产和娱乐文化的发展对维护甘肃省农业生态系统平衡十分重要。

对于某一区域农业生态服务功能的研究，要利用不同服务功能需要静态价值分析与动态价值分析相结合的方法，才能全面反映研究区域农业生态服务功能的变化与发展趋势，给农业生态环境的保护与利用提供一定的指示和预测作用。从农业生态服务价值构成、时间变化和不同服务功能类型动态、静态价值分析等方面对甘肃省农业生态服务功能进行研究发现，社会生态平衡和农业生态平衡至关重要。一方面，在大力发展农业经济的过程中要重视农业生态平衡的作用以及农业生态服务功能的重要性;另一方面，要推进生态农业建设。转结构、提效能，大力推行生态文明建设，不断保持农业生态系统平衡，持续提高农业生态系统服务功能价值。

参考文献：

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

[2] SWINTON S M，LUPI F，ROBERTSON G P，et al. Ecosystem services and agriculture：Cultivating agricultural ecosystems for diverse benefits[J]. Ecological economics，2007，64（2）：245-252.

[3] 杨正勇，杨怀宇，郭宗香. 农业生态系统服务价值评估研究进展[J].中国生态农业学报，2009，17（5）：1045-1450.

[4] 朱玉林，李明杰. 湖南省农业生态系统能值演变与趋势[J]. 应用生态学报，2012，23（2）：499-505.

[5] COSTANZA R，DéARGE R，GROOT R，et al. The total value of the world's ecosystem services and natural capital[J]. Nature，1997，387：253-260.

[6] TURNER R K，VAN DEN BERGH J C J M，S？DERQVIST T，et al. Ecological-economic analysis of wetlands： Scientific integration for management and policy[J]. Ecological economics，2000，35（1）：7-23.

[7] JING Y，ZHANG F，HE Y，et al. Assessment of spatial and temporal variation of ecological environment quality in Ebinur Lake Wetland National Nature Reserve，Xinjiang，China[J]. Ecological indicators，2020，110：1057-1064.

[8] MAES J，EGOH B，WILLEMEN L，et al. Mapping ecosystem services for policy support and decision making in the European Union [J]. Ecosystem services，2012，1（1）：31-39.

[9] DECLERCK F A J，JONES S K，ATTWOOD S，et al. Agricultural ecosystems and their services：The vanguard of sustainability？[J].Current opinion in environmental sustainability，2016，23：92-99.

[10] BAUMBER A. Enhancing ecosystem services through targeted bioenergy support policies[J]. Ecosystem services，2017，26（8）：98-110.

[11] 張玲玲，巩 杰，张 影. 基于文献计量分析的生态系统服务研究现状及热点[J]. 生态学报，2016，36（18）：5967-5977.

[12] SUTTON P C，COSTANZA R. Global estimates of market and non-market values derived from nighttime satellite imagery，land cover，and ecosystem service valuation[J]. Ecological economics，2002，41（3）：509-527.

[13] ARROW K，SOLOW R，PORTNEY P，et al. Report of the NOAA panel on contingent valuation[J]. Federal register，1993，58（10）：4601-4614.

[14] 虞依娜，彭少麟. 生态系统服务价值评估的研究进展[J]. 生态环境学报，2010，19（9）：2246-2252.

[15] 谢高地，肖 玉，甄 霖，等. 我国粮食生产的生态服务价值研究[J]. 中国生态农业学报，2005，13（3）：10-13.

[16] GOVENDER P，SIVAKUMAR V. Application of k-means and hierarchical clustering techniques for analysis of air pollution：A review （1980—2019）[J].Atmospheric pollution research，2020，11（1）：40-56.

[17] ZHANG B，TANG S，YANG R，et al. Streptomyces dangxiongensis sp. nov.，isolated from soil of Qinghai-Tibet Plateau[J]. International journal of systematic and evolutionary microbiology，2019，69（9）：2729-2734.

[18] CHRISTIE O H J，ESBENSEN K，MEYER T，et al. Aspects of pattern recognition in organic geochemistry[J].Organic geochemistry，1984，6：885-891.

[19] 王 云，潘竟虎. 基于生态系统服务价值重构的干旱内陆河流域生态安全格局优化——以张掖市甘州区为例[J]. 生态学报，2019，39（10）：3455-3467.

[20] 李双成，刘金龙，张才玉，等. 生态系统服务研究动态及地理学研究范式[J]. 地理学报，2011，66（12）：1618-1630.

[21] 谢高地，鲁春霞，冷允法，等. 青藏高原生态资产的价值评估[J]. 自然资源学报，2003，18（2）：189-196.

[22] 谈存峰，王生林. 兰州农业生态系统服务功能价值实证分析[J]. 西南农业大学学报（社会科学版），2012，10（4）：62-66.

[23] 毛爱华. 乡村旅游规划在生态环境保护下的策略研究[J]. 环境与发展，2019，31（10）：200，202.

[24] 赵 军，杨 凯. 生态系统服务价值评估研究进展[J]. 生态学报，2007，27（1）：346-356.

[25] 贺俊生，龚大洁，黄棨通，等. 祁连山区天祝县陆生野生脊椎动物多样性及地理区系研究[J]. 干旱区资源与环境，2020，34（1）：148-155.

[26] QIAOZHUO W，GUOFENG Z，HUIWEN G，et al. Influence of vegetation coverage and climate environment on soil organic carbon in the Qilian Mountains[J]. Scientific reports，2019，9（1）：285-298.

[27] 花 婷，王訓明，郎丽丽，等. 甘肃省气温与降水变化趋势及其对主要流域径流量的影响[J]. 中国沙漠，2015，35（3）：744-752.

[28] 李晓鸿. 甘肃生物多样性保护探讨[J]. 中国林业，2012（14）：31.

[29] 方彦杰，张绪成，侯慧芝，等. 近20年来甘肃省黄土高原区粮食产量的影响因素及未来粮食需求分析[J]. 中国农业科技导报，2015，17（4）：165-175.

[30] CRAIG P S，GIRAUDOUX P，SHI D，et al. An epidemiological and ecological study of human alveolar echinococcosis transmission in south Gansu，China[J]. Acta tropica，2000，77（2）： 167-177.

[31] HAO Y，ZHENG S，ZHAO M，et al. Reexamining the relationships among urbanization，industrial structure，and environmental pollution in China-New evidence using the dynamic threshold panel model[J]. Energy reports，2020，6：28-39.

[32] MA L，CHEN M，FANG F，et al. Research on the spatiotemporal variation of rural-urban transformation and its driving mechanisms in underdeveloped regions：Gansu province in western China as an example[J]. Sustainable cities and society，2019，50：1016-1028.

[33] HU Y，GAO M，BATUNACUN. Evaluations of water yield and soil erosion in the Shaanxi-Gansu Loess Plateau under different land use and climate change scenarios[J]. Environmental development，2020（34）：100488.

[34] 王立安，郭晓东，李亚兵. 西部县域文化生态旅游开发构思——以甘肃省武山县为例[J]. 干旱区资源与环境，2010，24（6）：196-200.

[35] 曹 昀，陈学林. 甘肃省生态旅游资源类型及特征[J]. 安徽农业科学，2009，37（10）：4795-4797.

[36] 曹庆楼，吴忠权. 乡村振兴视域下中国农业生态旅游发展创新[J]. 环境与可持续发展，2018，43（3）：15-17.

[37] 李胜男，张建军，陈 静，等. 农田生态系统服务及其价值量化研究进展[J]. 河北农业科学，2016，20（3）：87-94.

[38] 方 斌，王雪禅，魏巧巧. 以土地利用为视角的农田生态补偿理论框架构建[J]. 东北农业大学学报，2013，44（2）：98-104.