师银芳, 赵 军, 李传华

(1.西北师范大学 地理与环境科学学院， 甘肃 兰州 730070； 2.甘肃省遥感重点实验室, 甘肃 兰州 730070)

张掖市典型绿洲农田生态系统生产力时空分布特征

师银芳1,2, 赵 军1, 李传华1

(1.西北师范大学 地理与环境科学学院， 甘肃 兰州 730070； 2.甘肃省遥感重点实验室, 甘肃 兰州 730070)

[目的] 以甘肃省张掖市甘州区典型绿洲为例，研究绿洲农田生产力的时空格局及其全球变化对农田生产力的影响，为区域农业资源管理决策与规划提供决策依据。 [方法] 利用MOD17A3 NPP产品、气象数据及甘州区农业统计等数据，结合GIS空间分析技术开展研究。 [结果] 甘州区不同农田类型净初级生产力(NPP)均值最小值均出现在2001年，最大值出现在2007年。甘州区农田净初级生产力年均值变化范围为109.3～420.5 g/(m2·a)，其分布特征与农田类型的空间分布基本一致；同一农田类型，高海拔区农田生产力小于海拔较低区，即农田生产力的大小是由农田类型及海拔高度共同决定的。NPP变化趋势存在明显的空间异质性，农田生产力增加的面积占总面积的40.9%，主要是旱地和高海拔地区的水浇地。海拔1 800 m以上的旱地和水浇地年降水量的增加趋势远大于平均气温，其NPP呈明显的增加趋势。 [结论] 近10 a来甘州区农田生态系统生产力存在明显的时空差异，整个旱地和水浇地年均NPP呈增加趋势，水田年均NPP呈下降趋势。

农田生态系统; NPP; MOD17A3; 时空特征; 典型绿洲

文献参数： 师银芳, 赵军, 李传华.张掖市典型绿洲农田生态系统生产力时空分布特征[J].水土保持通报，2017,37(2)：120-125.DOI:10.13961/j.cnki.stbctb.2017.02.017； Shi Yinfang， Zhao Jun， Li Chuanhua. Spatial-temporal Distribution Characteristics of Farmland Ecosystem Productivity of Typical Oasis in Zhangye City[J]. Bulletin of Soil and Water Conservation, 2017,37(2)：120-125.DOI:10.13961/j.cnki.stbctb.2017.02.017

全球气候变化是人类迄今面临的最大的环境问题，也是21世纪人类面临的最复杂的挑战之一。而中国政府长期以来高度重视的粮食安全及农业生产与全球变化密切相关[1-3]。因此，开展农田生态系统生产力的时空格局及特征研究，不仅有利于对农田承载力进行估算，为农田的合理利用和科学管理提供依据，而且认识农田生态系统生产力变化的主导控制因素对农业可持续发展和粮食安全对策的中长期规划具有重要的意义。农田生态系统是在人类干预和控制下形成的人工生态系统，受人类干扰最为严重，与人类日常生活联系也最为密切[4]。净初级生产力(net primary productivity， NPP)是能够以统一的尺度标准体现生态系统生产力，能够避免以作物产量衡量农田生产力时由于农业机构调整、作物品种变化等干扰，是衡量农田生态系统生产力的一个很好指标[5]。而农田生态系统净初级生产力代表着农田生态系统固定大气CO2的能力，是陆地生态系统碳循环的重要组成部分[6-8]。随着遥感技术的不断发展，利用遥感技术进行农田NPP的研究已成为一种趋势。国内很多学者利用遥感数据对农田生产力进行了研究，主要包括3个方面：一方面利用统计数据来研究农田NPP，如莫宏伟等[9-11]利用农业统计数据对农牧交错区及区域尺度的农田生产力的时空格局及影响因子进行研究，表明农田存在明显的区域差异。另一方面利用遥感数据和NPP模型来模拟农田NPP，如闫慧敏等[3]利用NOAA/AVHRR遥感数据分析了中国农田生产力变化的空间格局及地形控制的作用，表明农田生产力降低与地形起伏度密切相关；刘文超等[12]采用GLO-PEM模型和VPM模型模拟了陕北地区农田生产力，并分析了耕地变化对农田生产力的影响。第三，直接利用NPP产品来研究农田NPP的变化及其影响因素，如罗玲等[13]利用MOD17A3 NPP产品对松嫩平原农田生产力时空特征与影响因素进行研究，表明松嫩平原水田与旱田总体均呈微弱的下降趋势，降水是影响研究区农作物生产的主要因素；朱峰等[14]利用MODIS NPP产品、土地利用数据及气象数据研究了东北地区农田NPP的时空特征及变化规律。综上可知，借助遥感手段研究农田生态系统生产力已成为一种趋势，而利用MODIS NPP数据产品研究农田生产力已得到广泛应用。因此，为了更全面的了解全球变化背景下农田生态系统生产力状况，在众多学者的研究基础上，以典型绿洲甘州区为例，利用MOD17A3 NPP数据，结合GIS空间分析技术，研究绿洲农田生产力的时空格局及其全球变化对农田生产力的影响，为区域农业资源管理决策与规划提供决策依据。

1 研究区概况

研究区选择在甘肃省河西走廊中部的甘州区，位于东经100°6′—100°52′，北纬38°32′—39°24′。南依祁连山，北邻内蒙古阿拉善右旗，东连金昌，西接酒泉。东西长约65 km，南北宽约98 km，总面积约4 240 km2。年平均降水量约115 mm，蒸发量2 047 mm，日照时数3 085 h左右，平均海拔约1 474 m。甘州区属温带大陆性气候，地势平坦，土地肥沃，水源丰富，日照充足，气候温和，是典型的绿洲农业和大型灌溉农业区，主要以水浇地为主，因此，合理的水肥配合是提高产量的有效措施[15]。甘州区主要盛产小麦、玉米、豆类及各种蔬菜，是全国重点商品粮基地之一。

2 数据与研究方法

2.1 数据来源及处理方法

研究数据主要包括：(1) NPP数据 MOD17A3影像数据，空间分辨率为1 km×1 km，利用MRT对数据进行拼接、裁剪、投影等预处理； (2) 土地利用数据 张掖市2007年土地利用/土地覆盖数据集及土地调查资料，来源于国家自然科学基金委员会寒区旱区科学数据中心，根据其得到甘州区农田类型数据； (3) 气象数据 张掖市及周边气象站的日平均气温、降水量数据，来源于中国气象科学数据共享服务网； (4) DEM数据 ASTER GDEM数据，其空间分辨率为30 m(http:∥www.nasa.gov/)； (5) 农业统计数据 收集2001—2010年甘州区农业统计数据，来源于张掖市统计年鉴资料。

2.2 研究方法

(1) 甘州区农田生产力分级方法。由于各农田类型作物产量的差异，采用像元尺度的距平分析法[16-17]计算不同农田类型NPP的距平百分率，计算公式为：

(1)

(2) 农田生产力年际变化趋势的计算。以像元为基础，利用一元线性回归趋势线分析法模拟2001—2010年甘州区农田生态系统NPP每个栅格像元的变化趋势，用来反映不同时期NPP变化趋势的空间特征。计算公式[18]为：

(2)

式中：θ——甘州区农田NPP的趋势斜率；n——监测时间段的年数； NPPi——第i年NPP的平均值。θ>0

表示NPP呈增加趋势，反之则表示NPP是减少的。

3 结果分析

3.1 甘州区农田生态系统NPP验证

根据甘州区农业统计数据，估算甘州区农田NPP，并对NPP模型模拟结果进行验证。根据不同农作物收获部分的含水量、收获指数、产量及播种面积，得到甘州区农田NPP，其转换方法[17,19]为：

(3)

式中：Yi——第i种农作物的产量；MCi——第i种农作物的含水率；HIi——第i种农作物的收获指数； Ai——第i种农作物的播种面积； 0.9——作物收获指数的调整系数； 0.45——将NPP转换为植物地上生物量的碳转换系数。参考国内外相关研究成果，主要农作物收获指数和含水率[20-21]详见表1。

表1 主要农作物的收获指数和含水量

采用上述方法，应用2001—2010年甘州区农业统计数据计算得到的NPP与遥感估算的NPP进行比较，结果如图1所示。由公式(3)可知，NPP是农作物产量和面积的函数，与作物产量成正比，与农作物播种面积成反比。可以看出，由统计数据估算得到的NPP与由遥感模型估算的农田NPP之间存在较好的拟合关系，R2为0.467。考虑到只用农田NPP验证，数据较少，为保证验证的可靠性，分别统计水田及其他农田类型的NPP，R2均达到0.36以上，拟合效果较好。综上表明MODIS遥感模型估算结果能够较准确地反映甘州区农田生产力的实际情况。

图1 MODIS NPP模拟值与由统计数据计算的NPP比较

3.2 甘州区农田生态系统生产力的时间变化

根据2001—2010年甘州区NPP数据及土地利用数据，得到不同类型农田的NPP，结果如图2所示。甘州区农田由水田、水浇地和旱地组成，其中水浇地面积占农田总面积的97%；2001—2010年水浇地的年均NPP变化范围为280.9～334.2 g/(m2·a)，10 a平均NPP为311.5 g/(m2·a)；水田和旱地面积远小于水浇地，近10 a水田的年均NPP为338.6 g/(m2·a)，略高于水浇地，旱田的NPP变化范围为100.7～201.0 g/(m2·a)，均小于水浇地和水田的多年NPP最小值。从变化趋势看，2001—2010年甘州区水浇地和水田NPP均值的各年变化步调一致，就时间角度而言，水浇地NPP呈上升趋势，水田NPP有轻微的下降趋势。近10 a来水田NPP均值均大于水浇地和旱地NPP，不同农田类型NPP均值最小值均出现在2001年，最大值出现在2007年。

图2 甘州区2001-2010年不同农田类型NPP均值变化

为了进一步分析各年农田生产力的变化，分别统计2001—2010年农田NPP的最大值、最小值及均值，结果如表2所示。从时间角度看，甘州区农田NPP最大值、最小值及均值变化存在明显的差异；2001—2010年农田NPP最小值出现在2005年，其值为58.1 g/(m2·a)；农田NPP最大值出现在2004年，其值为504.3 g/(m2·a)。2001—2010年农田NPP极差变化范围为275.4～423.1 g/(m2·a)，极差最大值出现在2003年，极差最小出现在2010年。2001年农田NPP均值最小，其值为281.34 g/(m2·a)，农田NPP均值2007年达到最大，为335.30 g/(m2·a)，这与各类型农田NPP均值变化一致。

表2 典型绿洲甘州区农田NPP的时间变化 g/(m2·a)

综上分析可知，甘州区各类型农田，面积最大的为水浇地，其次为旱地；农田NPP总量贡献最大的为水浇地，水田次之。2001—2010年，甘州区不同农田类型NPP均值的最小值均出现在2001年，最大值出现在2007年；近10 a来甘州区旱地和水浇地年均NPP呈增加趋势，水田年均NPP呈下降趋势。农田NPP极差变化范围较大，极差最大值和最小值分别出现在2003和2010年，其值分别为423.1，275.4 g/(m2·a)。

3.3 甘州区农田生态系统生产力的空间变化

根据甘州区2001—2010年农田NPP数据，得到农田NPP均值空间分布图(图3)。从图3中可知，整体而言，近10 a甘州区农田NPP年均值变化范围为109.3～420.5 g/(m2·a)，其分布特征与农田类型的空间分布基本一致。就各农田类型而言，水浇地主要分布在海拔1 400～1 700 m，其NPP大于220 g/(m2·a)；少量分布在海拔1 800 m以上，其NPP值也相对较小；旱地集中分布在甘州区南部地区，海拔1 800 m以上，NPP值小于190 g/(m2·a)；水田主要分布在海拔1 400～1 450 m，其NPP值较大，大于300 g/(m2·a)。

为了进一步分析农田NPP的空间分布差异，根据2001—2010年甘州区农田NPP数据，采用距平分析法，得到不同农田类型距平百分率，并根据距平百分率将农田分为5类：低产农田(<-20%)、较低产农田(-20%～-10%)、中产农田(-10%～10%)、较高产农田(10%～20%)、高产农田(>20%)，得到其空间分布图(图3)。由图3可知，中产和较高产农田所占面积最大，分别占农田总面积的25%和30%；较低产农田所占的面积最少，仅占0.6%；高产农田占农田总面积的17%，这与甘州区农田以水浇地为主是相一致的。

图3 典型绿洲甘州区农田生产力空间分布及分级

综上可知，2001—2010年甘州区农田NPP年均值变化范围为109.3～420.5 g/(m2·a)，其分布特征与农田类型的空间分布基本一致。水田和大部分水浇地主要分布在海拔较低地区，NPP值较大；旱地和少量水浇地分布在高海拔地区(>1 800 m)，NPP值较小，这表明农田生产力的大小是由农田类型及海拔高度共同决定的。

3.4 农田生产力与气候因子变化趋势分析

利用甘州区2001—2010年农田NPP数据，根据式(2)计算多年NPP变化趋势，得到农田生产力变化趋势空间分布图(图4)。从图4中可知，从空间格局整体分布而言，农田NPP呈增加和减少2种趋势，其变化率的均值为0.31 g/(m2·a)，旱地NPP增加趋势较为显著。局部空间变化而言，农田生产力增加的面积占总面积的40.9%，主要是低产和较低产农田，即旱地和高海拔地区的水浇地；NPP呈轻微减少趋势的农田主要是较高产的水浇地，海拔位于1 500～1 800 m。

图4 典型绿洲甘州区农田生产力变化趋势

为了进一步分析其产生的原因，根据2001—2010年甘州区及周边气象台站的平均气温和降水量数据，利用IDW+DEM插值，得到各年平均气温和年降水量栅格数据，然后计算平均气温和年降水量的空间变化趋势，并分析NPP变化趋势与年平均气温和年降水量变化的关系。近10 a年平均气温和年降水量均呈增加趋势，存在明显的空间异质性，其平均增加速率分别为0.08 ℃/a，4.85 mm/a；海拔1 800 m以上的旱地和水浇地年降水量的增加趋势远大于平均气温的，即其NPP增加趋势显著；而中低海拔地区的水浇地年降水年平均温度和年降水量均呈显著增加趋势，但其NPP有降低趋势。本文只考虑这2个气候因子，要了解农田生产力变化趋势及其与气候的关系，应综合考虑各气候因子之间的相互作用，以及考虑其他因素如灌溉方式和施肥措施等的影响，这将有待于进一步深入探讨和研究。

4 讨论与结论

(1) 2001—2010年旱地和水浇地年均NPP呈增加趋势，水田年均NPP呈下降趋势；甘州区各类型农田，面积最大的为水浇地，其次为旱地；农田NPP总量贡献最大的为水浇地，水田次之。

(2) 2001—2010年，不同农田类型NPP均值的最小值均出现在2001年，最大值出现在2007年；农田NPP极差变化范围较大，极差最大值和最小值分别出现在2003和2010年，其值分别为423.1，275.4 g/(m2·a)。

(3) 2001—2010年甘州区农田NPP年均值变化范围为109.3～420.5 g/(m2·a)，其分布特征与农田类型的空间分布基本一致，农田生产力的大小是由农田类型及海拔高度共同决定的。就农田NPP变化趋势的空间格局而言，农田NPP呈增加和减少2种趋势，其变化率的均值为0.31 g/(m2·a)，旱地NPP增加趋势较为显著。

(4) 本研究仅以典型绿洲甘州区为研究区，利用MOD17A3 NPP数据对其农田生产力的时空变化及其与气候因素的影响进行研究，并未考虑灌溉条件、施肥措施、化肥及农膜使用量等人为能量投入对农田生产力的影响；多年农田生产力变化的主导因素是气候因素还是人为能量投入在后续研究中仍需深入探讨分析，这对绿洲农田种植结构的合理化及可持续性发展具有重要意义。而借助遥感手段，准确估算像元尺度农田生产力的空间分布，对农田生态系统的碳汇、碳循环等的研究具有重要的意义，对农业整体素质和综合能力的提升具有重要意义。

[1] 肖国举,张强,王静.全球气候变化对农业生态系统的影响研究进展[J].应用生态学报,2007,18(8):1877-1885.

[2] 曹明奎.气候变化挑战农业系统和食物安全[J].科学时报,2004-05-13(1).

[3] 闫慧敏,刘纪远,曹明奎.中国农田生产力变化的空间格局及地形控制作用[J].地理学报,2007,62(2):171-180.

[4] 王留芳.农业生态学[M].陕西 西安:陕西科学技术出版社,1994.

[5] 王宗明,张柏,宋开山,等.松嫩平原作物现实生产力空间分布影响因素分析[J].干旱区资源与环境,2007,21(7):85-91.

[6] 赵其国,钱海燕.低碳经济与农业思考[J].生态环境学报,2009,18(5):1609-1614.

[7] Prince S D, Haskrtt J, Steinger M, et al. Net primary production of U. S. Midwest croplands from agricultural harvest yield data[J]. Ecological Applications, 2001,11(4):1194-1205.

[8] 师银芳,赵军.张掖市农田生态系统碳源/汇时空分布特征[J].水土保持通报,2015,35(4):283-287.

[9] 莫宏伟,刘礼英,任志远,等.农牧交错区耕地净第一性生产力动态变化研究:以陕西榆林市榆阳区为例[J].干旱地区农业研究,2007,25(1):15-20.

[10] 苏本营,张璐,陈圣宾,等.区域农田生态系统生产力的时空格局及其影响因子研究:以山东省为例[J].生态环境学报,2010,19(9):2036-2041.

[11] 张智全,黄高宝,李广.陇东耕地净第一性生产力与生态服务价值分析[J].中国沙漠,2011,31(6):1516-1520.

[12] 刘文超,颜长珍,秦元伟,等.近20年陕北地区耕地变化及其对农田生产力的影响[J].自然资源学报,2013,28(9):1373-1382.

[13] 罗玲,王宗明,宋开山,等.2000—2006年松嫩平原农田生产力时空特征与影响因素研究[J].农业系统科学与综合研究,2010,26(4):468-474.

[14] 朱峰,刘志明,王宗明,等.东北地区农田净初级生产力时空特征及其影响因素分析[J].资源科学,2010,32(11):2079-2084.

[15] 马强,宇万太,沈善敏,等.旱地农田水肥效应研究进展[J].应用生态学报,2007,18(3):665-673.

[16] 张树文,张养贞,李颖,等.东北地区土地利用/覆盖时空特征分析[M].北京:科学出版社,2006.

[17] 国志兴,王宗明,刘殿伟,等.基于MOD17A3数据集的三江平原低产农田影响因素分析[J].农业工程学报,2009,25(2):152-155.

[18] Stow D, Daeschner S, Hope A, et al. Variability of the seasonally integrated normalized difference vegetation Index across the north slope of Alaska in the1990s[J]. International Journal of Remote Sensing, 2003,24(5):1111-1117.

[19] Lobell D B, Hicke J A, Asner G P, et al. Satellite estimates of productivity and light use efficiency in the United States agriculture,1982-1998[J]. Global Change Biology, 2002,8(8):722-735.

[20] 张福春,朱志辉.中国作物的利用指数[J].中国农业科学,1990,23(2):83-87.

[21] 苏本营,张璐,陈圣宾,等.区域农田生态系统生产力的时空格局及其影响因子研究:以山东省为例[J].生态环境学报,2010,19(9):2036-2041.

Spatial-temporal Distribution Characteristics of Farmland Ecosystem Productivity of Typical Oasis in Zhangye City

SHI Yinfang1,2， ZHAO Jun1， LI Chuanhua1

(1.CollegeofGeographyandEnvironmentScience,NorthwestNormalUniversity,Lanzhou,Gansu730070,China; 2.KeyLaboratoryofRemoteSensing,Lanzhou,Gansu730070,China)

[Objective] Taking the typical oasis of Ganzhou District of Zhangye City, Gansu Province as an example, the temporal and spatial pattern of farmland productivity in that oasis and the effect of global change on the farmland productivity were researched, to provide decision-making basis for regional agricultural resource management and planning. [Methods] Using the data such as MOD17A3 NPP product, meteorological data and agricultural statistics data, the research was conducted by GIS spatial analysis technology. [Results] The minimum and maximum values of different types of net primary productivity(NPP) in Ganzhou District appeared in 2001 and 2007 respectively. The range of annual average farmland NPP in Ganzhou District was 109.3～420.5 g/(m2·a)， and its spatial distribution was approximately consistent with the farmland. For same type of farmland, its productivity of high elevation area was less than that of low elevation. All of which indicated that the agricultural productivity was determined by both of farmland type and altitude. The NPP trend had an obvious spatial heterogeneity. Area of increased agricultural productivity covered 40.9% of the total area, among which, the mainly increased types were dry land and irrigated cropland in high altitude region. The dry land and irrigated cropland above 1 800 m had a higher increasing trend of annual precipitation than that of average temperature, which made the NPP a significant increasing trend. [Conclusion] The ecosystem productivity of Ganzhou District from 2001 to 2010 had obvious spatial and temporal differences, the annual NPP of dry land and irrigated cropland exhibited an increasing trend, while the paddy NPP showed a decreasing trend.

farmland ecosystem productivity; NPP; MOD17A3; temporal and spatial variation; typical oasis

2016-06-12

2016-09-06

甘肃省遥感重点实验室基金项目“基于综合顺序分类理论的青藏高原区潜在植被GIS模拟的不确定性研究”； 国家自然科学基金项目“基于GIS的石羊河流域生态系统服务与权衡研究”(41661084)。

师银芳(1985—),女(汉族),甘肃省靖远县人,博士,讲师,研究方向为环境遥感与GIS应用。E-mail:amandasyf@126.com。

A

1000-288X(2017)02-0120-06

S181， S157