基于MODIS数据的中国耕地高中低产田空间分布格局

第70卷第5期2015年5月

地理学报

ACTAGEOGRAPHICASINICA

Vol.70,No.5May,2015

基于MODIS数据的中国耕地高中低产田

空间分布格局

冀咏赞1,2，闫慧敏1，刘纪远1，匡文慧1，胡云锋1

（1.中国科学院地理科学与资源研究所，北京100101；2.中国科学院大学，北京100049）

摘要：人口的持续增长和食物消费水平的快速提升使得中国粮食自给问题越发受到关注。后

备土地资源补给能力的不足和城市化过程对优质耕地的占用使得耕地资源“开源”和“节流”均存在一定困难，因此，提高耕地资源利用效率、提升耕地生产能力成为当前中国农业发展的根本策略，清晰地掌握全国高中低产田在空间上的分布区域成为国家中低产田改造与高标准农田建设规划的必要前提。本研究应用空间分辨率500m的MODIS遥感数据和光能利用率模型，在空间上清晰地估算现实农田生产力的基础上，结合高分辨率的耕地分布数据和耕作制度区划信息，探索新的高中低产田划分方法，掌握中国高中低产田的空间分布格局。该方法既能够体现耕地资源条件的区域分异规律，也能表达耕地现实生产能力的空间差异性，并且有效克服了利用统计数据进行高中低产田划分时以县为统计单元导致的县内高中低产田格局不明确的问题。基于该方法划分的高中低产田分别占全国耕地面积的20.66%、39.56%以及39.78%。其中，低产田约有3/4分布于丘陵山地区；高产田则53%分布于平原区。高产田面积最大的五个省均位于黄淮海区域内，其面积总和占全国高产田面积的41.75%。耕地面积位列全国前三的黑龙江省、四川省和内蒙古自治区其高产田面积占比均不足15%。关键词：粮食安全；光能利用率模型；农田生产力；高中低产田DOI:10.11821/dlxb201505008

1引言

人口和粮食消耗的不断增长给农业和自然资源资源带来了极大的压力，需要更高的粮食产量来保证未来全球的粮食安全[1]。目前，农业用地约占陆地面积的38%，是占地面积最大的土地利用方式[2-3]；而农业发展导致生物多样性下降、温室气体排放增多、环境污染等诸多问题，目前已成为导致全球环境退化的主要原因[1,4-5]。因此，从长远发展来看，全球农业面临着在保障粮食供给的同时减弱其对全球环境的负面影响这一挑战[1,4-9]。

目前许多耕地的粮食产量远低于其所处气候条件下的耕地生产潜力[1,5]，提升这些地区的耕地现实生产能力、提高耕地资源利用效率对应对挑战非常重要。研究表明，通过使用新技术、提升水肥利用效率、采用先进的农业管理手段等方法能够使耕地增产，大幅度提高粮食供给量[1,4-5,10-11]，同时保证环境影响最小化。中国作为世界上最大的粮食生

收稿日期：2014-07-25;修订日期：2014-11-28

基金项目：中国科学院科技服务网络计划(KFJ-EW-STS-019);中国科学院重点部署项目(KSZD-EW-Z-021-02);国家自

然科学基金重点项目(41430861)[Foundation:ScienceandTechnologyServiceNetworkInitiative,No.KFJ-EW-STS-019;TheKeyProgramoftheChineseAcademyofSciences,No.KSZD-EW-Z-021-02;NationalNaturalScienceFoundationofChina,No.41430861]

作者简介：冀咏赞(1989-),女,硕士研究生,主要从事土地利用变化与生态系统生产力研究。

E-mail:[email protected]

通讯作者：闫慧敏(1974-),女,副研究员,主要从事土地利用变化及其生态环境效应研究。E-mail:[email protected]

766-778页

产国和消费国，其粮食安全问题不仅关系到国家的稳定和发展，而且会对国际粮食供应

状况产生影响。中国农业发展的资源约束条件正日益严峻[10]，实现农业的持续稳定增长，根本出路在于加快农业科技创新，提高资源利用效率①，实现农业高产高效与资源永续利用协调兼顾。

为保护、支持农业发展，改善农业生产基本条件，提高农业综合生产能力和综合效益，自1988年开始进行大面积的以中低产田改造为主的农业综合开发，并于2009年起开始启动高标准农田建设工程。从1988年到2012年，开发范围涉及全国2045个县（市、区、旗）和222个国有农（牧、林）场，累计改造中低产田、建设高标准农田6.61亿亩②。

2013年初国务院批复了《国家农业综合开发高标准农田建设规划》，到2020年，全国将改造中低产田、建设高标准农田4亿亩。为了实现这一目标，需要在空间上清晰地掌握全国高中低产田的分布区域。

在现有的研究中，高中低产田的划分方法可归纳为以下三种：一是地力等级法，以耕地实地调研为基础，把耕地障碍因素作为划分耕地地力等级的关键指标，也就是将不存在和基本上不存在限制因素的耕地划归高产田，再根据限制因素的限制程度划分中低产田[13]；二是平均产量法[12-13]，以粮食平均单产为中值，上下各浮动一定的比重（一般为15%左右）或数值（一般为50kg），分别作为中产田的上限和下限；三是潜力产出率法[13]，根据耕地现实产出水平与其最大生产潜力的比值大小来划分高中低产田。这些方法有其各自的优点，但同时也存在相应的问题，地力等级法需要系统细致的野外调研工作，对各种耕地限制因素进行整合评判和分析，需要大量的人力物力和时间，难以实现耕地质量的快速动态监测与评估，不适于进行大尺度耕地质量评价；平均产量法依托于农业统计数据中的粮食产量，难以避免农业结构调整和作物品种变化等因素的干扰，同时每年一次的统计数据作为现阶段全国尺度上反映实际农田生产力的重要信息源，其最小空间单元是县，难以对特定地块耕地质量做出评价；潜力产出率法则需要对耕地最大生产潜力进行估算，计算复杂，过程繁琐。因此，使用新的数据与方法实现大尺度耕地质量的及时监测与有效评估对于加强耕地质量建设和合理利用土地具有重要意义。

遥感技术多星种、多传感器、多分辨率共同发展的局面为各个时空尺度下土地利用变化、农作物长势监测等方面的研究提供了数据基础[14]，同时生产力模型的发展使其越来越多地被应用于不同尺度不同类型生态系统生产力模拟及生产力对气候和土地利用变化的响应研究[15-18]。遥感技术的进步与农田生产力遥感模型的发展为高中低产田新划分方法的建立创造了可能。生态系统净初级生产力（NetPrimaryProductivity,NPP）是绿色植被在单位面积、单位时间内所累积的有机物数量，它直接反映耕地的现实生产能力且为各种作物生产能力提供一个统一的衡量标准。随着NPP遥感估算方法日趋成熟，现已可以作为监测高中低产田的有效且可行的度量指标。VPM模型是基于MODIS数据估算NPP的重要方法之一，较其他同类模型，该模型能够模拟农业多熟种植区的生态系统生产力，与通量观测和农业统计数据的对比研究证明在中国估算农田生态系统生产力的可靠性高于同类模型[16,19]。本研究基于2010年全国耕地成分栅格数据集和空间分辨率500m、时间分辨率8天的MODIS遥感数据，应用VPM模型结合中国耕作制度区划，探索中国高中低产田遥感监测方法，掌握中国高中低产田的分布规律与特点。

①陈锡文:走中国特色农业现代化道路,http://cpc.people.com.cn/GB/64162/82819/114926/114927/6840448.html(2008/01/30)

②王光坤:农业综合开发战略调向“退减保调转”,http://www.farmer.com.cn/xwpd/btxw/201312/t20131206_920451.htm(2013/12/06)

2数据与方法

2.1基于MODIS数据的生态系统生产力遥感模型AGRO-VPM

VPM（VegetationPhotosynthesisModel）模型[20]是一个基于卫星遥感数据进行生态系统生产力估算的光能利用率模型。它首先利用Monteith方程的原理对生态系统总初级生产力（GrossPrimaryProductivity，GPP）进行模拟，进而根据自养呼吸与GPP的比例关系计算NPP。VPM模型估算GPP可以表示为：

（1）GPP=εg×FPARchl×PAR

式中：εg为光能利用率（μmolCO2/μmolPPFD）；FPARchl表示被植被光合部分所吸收的光合有效辐射（PAR）比例；PAR表示光合有效辐射（μmolphotosyntheticphotonfluxdensity,PPFD）。

其中，生态系统的εg由最大光能利用率、温度、水分以及叶片物候期共同决定，即：

εg=ε0×Tscalar×Wscalar×Pscalar（2）

式中：ε0为最大光能利用率（μmolCO2/μmolPPFD），通过应用CO2涡度通量塔测量的生态系统CO2净交换（NEE）和incidentPAR（μmol/m2/sPPFD）数据由米氏方程计算求得[16]。Tscalar、Wscalar和Pscalar分别是温度、水分以及叶片物候期对最大光能利用率的调节系数。以上参数的获取及计算方法详见KalfasJL[21]一文。

根据对小麦、向日葵、高粱、豆类等作物呼吸作用和生长速率全过程实验结果，作物生长过程中自养呼吸与光合作用生成GPP的比例是不变的，呼吸消耗与GPP的比例系数为0.42[22-25]。

AGRO-VPM模型在农田生态系统中得到了充分的验证，中国学者将其应用到东北春玉米和退化草地[26]以及禹城冬小麦—夏玉米轮作农田生态系统[19]中，均取得较好的模拟效果。2.2数据

（1）耕地数据

本研究中使用的2010年耕地数据来源于中国科学院地理科学与资源研究所“全国耕地资源遥感监测与耕地保护对策研究”项目组，该数据是利用LandsatTM/ETM进行影像解译，通过实地调查和精度验证生成1:10万矢量图的基础上统计生成。每个栅格对应的值为该栅格内耕地面积所占的万分比，其分辨率为1km。

（2）中国耕作制度区划

由于中国自然及气候条件的地域差异性极大，而依据热量、水分、地貌等环境指标划分的中国耕作制度区划使得每个区域内部具有自然条件和种植制度的相对均一性，因此利用耕作制度区划对全国进行分区，进而基于净初级生产力对各个耕作区进行高中低产田划分可以有效地避免全国不同区域因条件差异导致的耕地质量不可比的问题。中国耕作制度区划[27]根据热量条件、水分、地貌及社会经济条件、耕作制与熟制特点、作物类型，将中国耕作制度分成3个带，12个一级区，38个二级区。分区时兼顾了自然条件与社会经济条件，基本上保持了县级行政单元的完整性（图1）。

（3）MODIS数据

VPM模型的输入数据为增强型植被指数EVI（EnhancedVegetationIndex）以及地表水分指数LSWI（LandSurfaceWaterIndex），利用空间分辨率500m、时间分辨率8天的MODIS地表反射率产品（MOD09A1）（http://www.edc.usgs.gov/）计算可得。

EVI考虑了土壤背景对植被指数的影响，对气溶胶等残留做了进一步校正[28]。

注：引自刘巽浩《中国耕作制度》。0级带分界线将全国分为一熟区、二熟区和三熟区，1-12为12个一级区代码，在

图中用相同颜色标示，将一级区细分得到38个二级区，在图中用数字标出。二级区代码对应名称在表1中列出。

图1中国耕作制度区划图

Fig.1ThefarmingsystemzoningmapofChina

ρnir-ρred

（3）

nirredblue式中：G=2.5，C1=6，C2=7.5，L=1；ρnir、ρred和ρblue分别表示近红外波段、红波段和蓝波段的反射率。

LSWI通过利用短波红外（SWIR）波段对水分含量的敏感性来进行计算。

ρ-ρ

LSWI=nirswir（4）

nirswir

式中：ρnir和ρswir分别表示近红外和短波红外波段的反射率。

温度数据由来自国家气象局的中国气象科学数据共享服务网（http://www.cma.gov.cn/2011qxfw/2011qsjgx/）提供的2000-2008年的气象站点气温数据经ANUspline软件进行样条插值得到，是与遥感数据时空分辨率（8天，500m）匹配的栅格数据；PAR数据为利用MODIS地表反射率产品反演得到的与遥感数据时空分辨率（8天，500m）匹配的栅格数据。

2.3高中低产田划分方法

2.3.1高中低产田定义本研究中，中低产田是指目前的产出水平远未达到所处的自然和社会经济条件下应有的生产能力，具有较大增产潜力的耕地；高产田是指不存在或较少存在制约农业生产的限制因素，生产能力较高的耕地。

2.3.2VPM模型模拟利用VPM模型模拟2000-2008年全国耕地净初级生产力，为了消除气候波动以及农业管理措施变化等因素的影响，计算9年NPP的均值作为高中低产田划分的数据基础。

耕地净初级生产力大体呈现出从西北向东南逐渐增加的趋势，这与中国水热条件的

EVI=G×

变化相一致。南方三熟区耕地分布虽不集中，但产能最高；黄淮海地区耕地分布广泛且连续，大部分耕地产能较高；四川盆地和东北部分地区产能次之。值得注意的是新疆西北部和西部有少量耕地产能很高。宁夏和甘肃中部、陕西北部为全国耕地产能较低区域。2.3.3高中低产田产能区间界定在确定高中低产田标准过程中，为了降低遥感信息中混合像元对划分结果的影响，首先提取耕地面积比例大于等于50%的栅格所对应的净初级生产力用于划分标准的制定，获取各个耕作区内NPP值的频率分布直方图。

耕作区存在四种典型的频率分布直方图（图2），每种均包含一个主要的正态分布，该正态分布前后存在一定的噪声。由于噪声的存在会影响每个耕作区内高中低产田划分标准的确定，因此将正态分布前后两个拐点对应的NPP值作为分界值，落在正态分布分界值以内的栅格则为有效栅格。定义分界值为NPPa和NPPb（NPPa

对耕地进行高中低产田划分的依据实质上是按照平均分配原则将耕地分为三类，低产田上限标准为NPPa+1/3×NPPdif，高产田下限标准为NPPb–1/3×NPPdif。但是，考虑到随着中低产田改造的实施，低产田面积比例减少，中高产田面积比例增加，因此将低、中、高产田所占比例进行适当调整，调整后的占比分别为30%、35%和35%。即高中低产田标准的划分界限为：

低产田上限标准=NPPa+NPPdif×30%高产田下限标准=NPPb–NPPdif×35%

以此为标准对各个耕作区内所有耕地进行高中低产田划分。各耕作区内有效栅格NPP范围及高中低产田划分标准见表1

。

图2四种典型的频率分布直方图示例

Fig.2Histogramsoftypicalfrequencydistribution

表1

Tab.1

耕作区有效栅格NPP范围及高中低产田划分标准(gC/m2/a)

mediumandlowyieldcropland(gC/m2/a)

NPPrangeofeffectivegridsineachfarmingareaandcriteriafortheclassificationofhigh,

低产田

上限[***********][***********][***********][***********][***********][***********]675835

高产田下限[***********][***********][***********][***********][***********][***********]209381168

提升潜力[***********][***********][***********][***********][***********][***********]33

序号[***********][***********][***********][1**********]38

二级区代码及名称

11-藏东南川西谷地喜凉作物一熟区12-海北甘南高原喜凉作物一熟轮歇区

21-后山坝上晋北高原山地半干旱喜凉作物一熟区22-陇中青东宁中南黄土丘陵半干旱喜凉作物一熟区31-辽吉西蒙东南冀北半干旱喜温一熟区32-黄土高原东部易旱喜温一熟区33-晋东半湿润易旱一熟填闲区

34-渭北陇东半湿润易旱冬麦一熟填闲区41-大小兴安岭山麓岗地凉温作物一熟区42-三江平原长白山地温凉作物一熟区43-松嫩平原喜温作物一熟区

44-辽河平原丘陵温暖作物一熟填闲区51-河套、河西灌溉一熟填闲区52-北疆灌溉一熟填闲区53-南疆、东疆绿洲二熟一熟区

61-燕山太行山山前平原水浇地套复二熟旱地一熟区62-黑龙港缺水低平原水浇二熟旱地一熟区63-鲁西北豫北低平原水浇地粮二熟棉一熟区64-山东丘陵水浇地二熟旱坡地花生棉花一熟区65-黄淮平原南阳盆地旱地水浇地两熟区66-汾渭谷地水浇地二熟旱地一熟二熟区67-豫西丘陵山地旱坡地一熟水浇地二熟区71-秦巴山区旱地二熟一熟兼水田二熟区72-川鄂湘黔低高原山地水田旱地二熟兼一熟区73-贵州高原水田旱地二熟一熟区74-云南高原水田旱地二熟一熟区

75-滇黔边境高原山地河谷旱地一熟二熟水田二熟区81-江淮平原麦稻两熟兼旱三熟区

82-鄂豫皖丘陵平原水田旱地两熟兼早三熟区91-盆西平原水田麦稻两熟填闲区92-盆东丘陵低山水田旱地两熟三熟区101-沿江平原丘陵水田早三熟二熟区102-两湖平原丘陵水田中三熟二熟区111-浙闽丘陵山地水田旱地三熟二熟区112-南岭丘陵山地水田旱地三熟二熟区113-滇南山地旱地水田二熟兼三熟区121-华南低平原晚三熟区

122-华南沿海西双版纳台南二熟三熟与热作区

NPPa[***********][***********][***********][***********][***********][***********]450550

NPPb[***********][***********][***********][***********][***********][***********][***********]500

NPPdif[***********][***********][***********][***********][***********][***********]750950

3结果与分析

高中低产田的划分标准体现了由耕地自然生产潜力所决定的耕地产能的区域差异，低产田上限和高产田下限总体上的规律为三熟区最高、二熟区居中、一熟区最低（图3）。水热充沛地区的低产田实际产能会高于干旱寒冷区的高产田产能，在38个二级耕作区中，低产田上限最高的耕作区为滇南山地旱地水田二熟兼三熟区（113），其上限标准达840gC/m2/a，高于27个二级耕作区的高产田下限标准。

图3高中低产田划分标准变化趋势(横轴对应表1中38个二级区序号)

correspondstothesequencenumberofTab.1)

Fig.3Changetrendsofcriteriafortheclassificationofhigh,mediumandlowyieldcropland(Thehorizontalaxisdigital

高产田下限与低产田上限的差值意味着各耕作区内低产田改造的提升潜力，为了掌握各个区域的产能提升潜力，根据高产田下限与低产田上限的差值将各耕作区划分为250四类产能提升类型区。其中黄土高原东部易旱喜温一熟区（32）、辽河平原丘陵温暖作物一熟填闲区（44）、盆西平原水田麦稻两熟填闲区（91）和松嫩平原喜温作物一熟区（43）四个耕作区提升潜力最低；北疆灌溉一熟填闲区（52）、浙闽丘陵山地水田旱地三熟二熟区（111）、华南低平原晚三熟区（121）、云南高原水田旱地二熟一熟区（74）、滇南山地旱地水田二熟兼三熟区（113）和华南沿海西双版纳台南二熟三熟与热作区（122）六个耕作区提升潜力最高（表1）。

中国高、中、低产田面积分别占全国耕地面积的20.66%、39.56%以及39.78%。高产田主要分布在华北平原、长江中下游平原、四川盆地、吉林省中部以及辽宁省东北部，这些区域地势平坦，水资源充沛，为农作物生长提供了良好的条件；高产田在陕西中部、宁夏北部、甘肃中部及南部、新疆西北部地区也有少量分布，这些地区气候条件相对较差，但河流的流经为水资源的供给、农业的发展提供了保障，使其相较同一耕作区内其他耕地生产力更高。中产田集中分布于东北平原、四川盆地、河南省南部、安徽省中北部，这些地区多位于高产田分布区附近，水土条件也比较适宜发展农业，中产田在湖南、江西、广东、广西等地区也有零散分布。低产田则多分布于黄土高原、黄河入海口处、云贵高原等自然和气候条件较恶劣地区（图4）。

除热量、水分等已体现在中国耕作制度区划中的因素之外，地形条件也是导致耕地生产力区域差异的因素之一。地形起伏度是较好的区域尺度地形指标，能够反映地表的起伏变化，用地面一定距离范围内最大的高程差表示[29]。本研究基于1km的数字高程模型，以5km×5km为适宜尺度分析窗口[30]，提取中国区域地形起伏度，将小于30m的地区定义为平原，大于30m的地区定义为丘陵山地[29]。平原区耕地约占全国耕地面积的40%，丘陵山地区耕地约占全国耕地面积的60%。平原区中中产田占比最高，为48.14%，其次为高产田，占27.95%，低产田仅占平原区耕地的23.91%

；丘陵山地区中低

产田面积最大，占比为

49.15%，其次为中产田，占34.61%，高产田仅占丘陵山地区的16.24%（图5a）。

在平原区，中、高产田的比例皆高于全国平均水平，低产田比例则远远低于全国平均水平，而在丘陵山地区，低产田占比接近50%，高于全国平均值39.78%，中、高产田比例则低于全国平均状况。就三类农田的地形组成结构来看，约3/4的低产田分布于丘陵山地

图4全国高中低产田分布图

Fig.4Thedistributionmapofhigh,mediumandlowyieldcroplandin

China

图5不同地形区以及不同等级耕地的组成结构

Fig.5Thecompositionstructureofdifferentterrainareasandcroplands

区，中产田在平原区和丘陵山地区的分布大致相当，而高产田中约有53%分布于平原区，47%位于丘陵山地区（图5b）。

高产田面积排名全国前五的省份分别为河南省、山东省、江苏省、河北省和安徽省，均位于黄淮海区域内，其次为湖北省、吉林省和四川省。高产田在这8个省份内集中连片分布，占全国高产田面积的57.3%。中低产田在黑龙江省、四川省、内蒙古自治区、河南省和河北省分布较多，占全国中低产田面积的34.88%。河南省、河北省和山东省耕地面积大，高产田和中低产田均有广泛分布（表2）。

全国共有12个省（自治区）内高产田比例高于全国平均水平，江苏省更是以接近50%的高产田比例位居全国首位，同时省内低产田比例不足15%，其次分别为河南省、山东省、湖北省、安徽省、河北省、吉林省、新疆维吾尔族自治区、台湾省、辽宁省、海南省以及广东省，省/自治区内高产田分别占34.69%、32.86%、31.73%、28.80%、27.82%、26.40%、25.18%、23.47%、22.96%、21.17%和20.78%。低产田比例高于全国平均水平的有17个省（直辖市、自治区）

，其中西藏自治区低产田比例全国最高

表2

各省高中低产田面积(万km2)

行政区江西省甘肃省广西浙江省重庆市湖南省福建省贵州省宁夏海南省青海省台湾省上海市西藏天津市北京市

高产田面积

0.890.870.570.560.550.500.440.420.270.180.160.150.070.060.060.05

中产田面积

1.331.772.041.161.322.000.501.250.510.440.220.310.250.070.350.23

Tab.2Theareasofhigh,mediumandlowyieldcroplandineveryprovince(tenthousandkm2)

行政区河南省山东省江苏省河北省安徽省湖北省吉林省四川省新疆陕西省辽宁省黑龙江省云南省内蒙古山西省广东省

高产田面积

3.723.353.252.682.282.121.961.611.481.481.461.411.071.030.990.90

中产田面积低产田面积

4.474.952.514.264.332.523.155.112.091.892.948.781.614.062.161.50

2.541.901.002.711.302.042.325.312.303.791.965.614.176.162.951.94

低产田面积

2.183.902.511.041.973.511.243.301.040.250.430.180.120.330.27

0.20

（71.8%），而高、中产田各只占不足15%，其次分别为贵州省、云南省、甘肃省、湖南省、福建省、宁夏回族自治区、内蒙古自治区、陕西省、青海省、重庆市、江西省、广西壮族自治区、山西省、广东省、四川省和北京市，省/市/自治区内低产田分别占66.48%、60.85%、59.64%、58.37%、56.98%、56.95%、54.76%、52.99%、52.57%、51.29%、49.47%、49.00%、48.41%、44.70%、44.13%和41.20%（图6）。

图6各省高中低产田所占耕地比例

Fig.6Theproportionsofhigh,mediumandlowyieldcroplandineachprovince

4结论与讨论

4.1结论

本研究基于农田生产力遥感模型划分的高中低产田分别占全国耕地面积的20.66%、39.56%以及39.78%。约3/4的低产田分布于丘陵山地区，丘陵山地区的高产田主要分布在吉林中东部、甘肃中部及南部、新疆西北部、四川及重庆东北部、湖北东部等地区；

而平原区仍有近一半（48.14%）的耕地为中产田，高产田仅占27.95%，平原区的低产田

主要分布在黄河入海口处、湖南和江西中北部以及黑龙江和吉林的东部等地区。

位于黄淮海粮食主产区的江苏省、河南省、山东省、安徽省和河北省高产田面积总和占全国高产田面积的41.75%，位于长江中下游平原的湖北省高产田面积占全国高产田面积的5.78%，位于东北平原的吉林省高产田面积占全国高产田面积的5.36%。

中低产田面积占全省耕地总面积比例高于85%的省份有贵州省、湖南省、黑龙江省、天津市、内蒙古自治区、北京市、广西壮族自治区、甘肃省、四川省、西藏自治区和重庆市；低产田面积比例最高的三个省份为西藏自治区、贵州省和云南省，低产田面积占全省耕地总面积的60%以上。耕地面积位列全国前三的黑龙江省、四川省和内蒙古自治区（表2），其高产田面积占比均不足15%，中产田面积分别为8.78万km2、5.11万km2和4.06万km2，占全省耕地总面积的55.57%、42.46%和36.08%；低产田面积分别为5.61万km2、5.31万km2和6.16万km2，占全省耕地总面积的35.48%、44.13%和54.76%。黑龙江省和内蒙古自治区的中低产田中分别有61.77%和26.33%分布于平原区，四川省的中低产田中有88.01%分布于丘陵山地区。4.2讨论

在已有的高中低产田划分研究中，《我国中低产田分布及粮食增产潜力研究》（1988）采用地力等级法将耕地分为七个等级，之后按照等级进行高中低产田划分，得到高产田占全国耕地比例为32.16%，中产田占32.90%，低产田占34.94%[13]；《中国耕地资源及其开发利用》（1992）采用平均产量法划分的高中低产田分别占比21.54%、37.24%和41.22%[13,31]；潜力产出率法划分的高中低产田分别占比34.92%、41.95%和23.14%[13]。本研究中的高中低产田所占比例与《中国耕地资源及其开发利用》中的结果较为相近，但与另外两种划分结果相比本研究高产田所占比例低于另外两种方法。比较发现，《中国耕地资源及其开发利用》和本文的划分方法相对更为保守，在全国东北区（黑吉辽）、华北区（京津冀鲁豫）、黄土高原区（晋陕甘）、西北干旱区（新蒙宁）、青藏高原区（青藏）、长江中下游区（苏皖鄂湘赣浙沪）、华南区（粤桂闽琼台）和西南区（川滇渝黔）等八大区中，三个区的耕地以中产田为主，五个区的耕地以低产田为主，而潜力产出率法的划分结果中有两个区的耕地以高产田为主，四个区的耕地以中产田为主，仅有两个区的耕地以低产田为主。划分结果出现差异的一个原因是每一种划分方法均无法避免人为主观因素的干扰，但这是中低产田划分方法共有的不足之处，也是导致本方法划分结果不确定性的原因之一。已有方法多以县为最小调查单元，高中低产田在栅格级别上的空间分布格局不明确，但是就高中低产田在全国八大区的分布情况来看，各个方法的划分结果均能够反映区域尺度的耕地质量状况：高产田所占比例高于全国平均水平的区域为水热条件较好的华北区和长江中下游区，低产田所占比例远高于全国平均水平的区域为不太适宜发展农业的黄土高原区和西北干旱区。

农业的不断发展带来了诸多环境问题，城市化进程的加快导致农业生产资源受限，因此提高资源利用效率是实现农业持续稳定增长、保障国家粮食安全的必由之路。在空间上清晰地识别全国高中低产田的分布区域是进行农业综合开发中低产田改造的前提。遥感技术的不断进步和生态系统生产力遥感模型的日益成熟为基于净初级生产力划分高中低产田创造了条件。基于遥感的生态系统生产力模型能够清晰表达耕地NPP的空间格局，很好地解决了利用统计数据进行高中低产田划分时以县为统计单元导致的县内高中低产田格局不明确的问题。同时，遥感信息综合反映了自然条件和管理因素对农作物生长的影响。耕作制度区划的引入使高中低产田的划分基于耕地所在耕作区的自然和社会

条件，因此耕地质量并非由耕地净初级生产力值在全国范围内的高低决定，而是由其在相应耕作区内的高低决定。这样可以使得划分入某一类的耕地虽然生产力高低不同，但均能反映其在当前的气候和自然条件下是否需要进行改造及改造后增产潜力的高低。目前，黄土高原东部易旱喜温一熟区（32）、辽河平原丘陵温暖作物一熟填闲区（44）、盆西平原水田麦稻两熟填闲区（91）和松嫩平原喜温作物一熟区（43）四个耕作区的低产田提升潜力最低，小于100gC/m2/a；北疆灌溉一熟填闲区（52）、浙闽丘陵山地水田旱地三熟二熟区（111）、华南低平原晚三熟区（121）、云南高原水田旱地二熟一熟区

（74）、滇南山地旱地水田二熟兼三熟区（113）和华南沿海西双版纳台南二熟三熟与热作区（122）六个耕作区的低产田提升潜力最高，大于250gC/m2/a，其余28个耕作区的低产田提升潜力均在100~250gC/m2/a之间。

本研究应用空间分辨率500m的MODIS数据估算耕地生产力，为了消除混合像元对高中低产田划分标准的影响，在确定各个耕作区的划分标准时使用了耕地占栅格面积大于等于50%的栅格对应的NPP数据，使得农田为栅格内主要的土地利用方式，但对于耕地所占面积比例低的栅格，该方法并不能降低混合像元影响产生的结果不确定性。因此，更为精细的区域中低产田改造规划需要更高分辨率的生产力数据的支持。此外，耕地产能的稳定程度体现了耕地在外界条件（主要是气候因素）突变的情况下保持其原有生产力的能力，不断提高高产稳产耕地的数量是实现粮食供给持续增长、保障国家粮食安全的有效途径。随着高标准农田建设步伐的推进，如何有效检测高标准农田的分布将是支撑未来现代农业开发规划亟待解决的技术之一。

参考文献(References)

[1]FoleyJA,RamankuttyN,BraumanKA,etal.Solutionsforacultivatedplanet.Nature,2011,478(7369):337-342.

[2]FoodandAgricultureOrganizationoftheUnitedNations(FAOSTAT).http://faostat.fao.org/site/567/default.aspx#ancor

(accessed,March2011).

[3]RamankuttyN,EvanAT,MonfredaC,etal.Farmingtheplanet:1.Geographicdistributionofglobalagriculturallands

intheyear2000.GlobalBiogeochemicalCycles,2008,22(1).

[4]FoleyJA,DeFriesR,AsnerGP,etal.Globalconsequencesoflanduse.Science,2005,309(5734):570-574.

[5]TilmanD,BalzerC,HillJ,etal.Globalfooddemandandthesustainableintensificationofagriculture.Proceedingsof

theNationalAcademyofSciences,2011,108(50):20260-20264.

[6]SeufertV,RamankuttyN,FoleyJA.Comparingtheyieldsoforganicandconventionalagriculture.Nature,2012,485

(7397):229-232.

[7]GodfrayHCJ,BeddingtonJR,CruteIR,etal.Foodsecurity:Thechallengeoffeeding9billionpeople.Science,

2010,327(5967):812-818.

[8]RobertsonGP,SwintonSM.Reconcilingagriculturalproductivityandenvironmentalintegrity:Agrandchallengefor

agriculture.FrontiersinEcologyandtheEnvironment,2005,3(1):38-46.

[9]CassmanKG,DobermannA,WaltersDT,etal.Meetingcerealdemandwhileprotectingnaturalresourcesand

improvingenvironmentalquality.AnnualReviewofEnvironmentandResources,2003,28(1):315-358.

[10]ShiQuanhong,WangHong,ChenFu,etal.Thespatial-temporaldistributioncharacteristicsandyieldpotentialof

medium-lowyieldedfarmlandinChina.ChineseAgriculturalScienceBulletin,2010,26(19):369-373.[石全红,王宏,陈阜等.中国中低产田时空分布特征及增产潜力分析.中国农学通报,2010,26(19):369-373.]

[11]ZhangLin,ZhangFengrong,JiangGuanghui,etal.Potentialimprovementofmedium-lowyieldedfarmlandand

guaranteeoffoodsafetyinChina.ResearchofAgriculturalModernization,2005,26(1):22-25.[张琳,张凤荣,姜广辉等.我国中低产田改造的粮食增产潜力与食物安全保障.农业现代化研究,2005,26(1):22-25.]

[12]WuBingfang,DuXin,MengJihua,etal.Aremotesensingmethodtoclassifyhigh,mediumandlowyieldcroplandand

itsapplicationinHuang-Huai-HaiBasin//ModernAgriculturalDevelopmentandNationalFoodSecurity&Huang-Huai-HaiModernAgriculturalDevelopmentStrategyForumProceedings,2009:104-108.[吴炳方,杜鑫,蒙继华,等.高中低产田遥感监测方法及在黄淮海流域的应用//现代农业发展与国家粮食安全暨黄淮海现代农业发展战略高峰论坛

论文集,2009:104-108.][13]LinPengsheng.Studyonthedistributionandpossibleproductionincrementofmediumandlow-yieldfarmlandin

China[D].Beijing:ChineseAcademyofAgriculturalSciences,2008.[林鹏生.我国中低产田分布及增产潜力研究

[D].北京:中国农业科学院,2008.]

[14]LiFolin,LiBenxun,CaoWeixing.Statusandprospectofcrop-yieldassessbyremotesensing.JournalofYunnan

AgriculturalUniversity,2005,20(5):680-684.[李佛琳,李本逊,曹卫星.作物遥感估产的现状及其展望.云南农业大学学报,2005,20(5):680-684.]

[15]LiuWenchao,YanChangzhen,QinYuanwei,etal.Croplanddynamicsanditsinfluenceontheproductivityinrecent20

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

[16]YanHuimin,LiuJiyuan,HuangHeqing,etal.Impactsofcroplandtransformationonagriculturalproductionunder

urbanizationandGrainforGreenProjectinChina.ActaGeographicaSinica,2012,67(5):579-588.[闫慧敏,刘纪远,黄河清,等.城市化和退耕还林草对中国耕地生产力的影响.地理学报,2012,67(5):579-588.]

[17]HuangMei,JiJinjun,PengLili.Theresponseofvegetationnetprimaryproductivitytoclimatechangeduring1981-

2000intheTibetanPlateau.ClimaticandEnvironmentalResearch,2008,13(5):608-616.[黄玫,季劲钧,彭莉莉.青藏高原1981-2000年植被净初级生产力对气候变化的响应.气候与环境研究,2008,13(5):608-616.]

[18]WangYuan,HuangMei,WangXiangrong.Impactsoflanduseandclimatechangeonagriculturalproductivityin

Shanghai.ActaScientiaeCircumstantiae,2010,30(3):641-648.[王原,黄玫,王祥荣.气候和土地利用变化对上海市农田生态系统净初级生产力的影响.环境科学学报,2010,30(3):641-648.]

[19]YanHM,FuYL,XiaoXMetal.Modelinggrossprimaryproductivityforwinterwheat-maizedoublecroppingsystem

usingMODIStimeseriesandCO2eddyfluxtowerdata.Agriculture,Ecosystems&Environment,2009,129(4):391-400.

[20]XiaoXM,HollingerD,AberJ,etal.Satellite-basedmodelingofgrossprimaryproductioninanevergreenneedleleaf

forest.RemoteSensingofEnvironment,2004,89(4):519-534.

[21]KalfasJL,XiaoX,VanegasDX,etal.Modelinggrossprimaryproductionofirrigatedandrain-fedmaizeusing

MODISimageryandCO2fluxtowerdata.AgriculturalandForestMeteorology,2011,151(12):1514-1528.

[22]GiffordRM.WholeplantrespirationandphotosynthesisofwheatunderincreasedCO2concentrationandtemperature:

long-termvs.short-termdistinctionsformodelling.GlobalChangeBiology,1995,1(6):385-396.

[23]LuoY,SimsDA,ThomasRB,etal.SensitivityofleafphotosynthesistoCO2concentrationisaninvariantfunctionfor

C3plants:Atestwithexperimentaldataandglobalapplications.GlobalBiogeochemicalCycles,1996,10(2):209-222.

[24]ChengW,SimsDA,LuoY,etal.Photosynthesis,respiration,andnetprimaryproductionofsunflowerstandsin

ambientandelevatedatmosphericCO2concentrations.GlobalChangeBiology,2000,6(8):931-941.

[25]AlbrizioR,StedutoP.Photosynthesis,respirationandconservativecarbonuseefficiencyoffourfieldgrowncrops.

AgriculturalandForestMeteorology,2003,116(1):19-36.

[26]WangZ,XiaoX,YanX.Modelinggrossprimaryproductionofmaizecroplandanddegradedgrasslandinnortheastern

China.Agric.For.Meteorol.,2010,150(9):1160-1167.

[27]LiuXunhao.ChineseCroppingSystem.Beijing:AgriculturePress,1993:12-15.[刘巽浩.中国耕作制度.北京:农业出

版社,1993:12-15.]

[28]LiWenmei,QinZhihao,LiWenjuan,etal.ComparisonandanalysisofMODISNDVIandMODISEVI.Remote

SensingInformation,2010,(6):73-78.[李文梅,覃志豪,李文娟,等.MODISNDVI与MODISEVI的比较分析.遥感信息,2010,(6):73-78.]

[29]LiuXinhua,YangQinke,TangGuoan.ExtractionandapplicationofreliefofChinabasedonDEMandGISmethod.

BulletinofSoilandWaterConservation,2001,21(1):57-60.[刘新华,杨勤科,汤国安.中国地形起伏度的提取及在水土流失定量评价中的应用.水土保持通报,2001,21(1):57-60.]

[30]ZhangWei,LiAinong.StudyontheoptimalscaleforcalculatingthereliefamplitudeinChinabasedonDEM.

GeographyandGeo-InformationScience,2012,28(4):8-12.[张伟,李爱农.基于DEM的中国地形起伏度适宜计算尺度研究.地理与地理信息科学,2012,28(4):8-12.]

[31]AgriculturalNaturalResourcesandAgriculturalRegionalPlanningInstituteofChineseAcademyofAgricultural

Sciences,NationalSoilandFertilizerStation.ChineseArableLandResourcesanditsExploitation.Beijing:SurveyingandMappingPress,1992:15.[中国农业科学院农业自然资源和农业区划研究所,农业部全国土壤肥料总站.中国耕地资源及其开发利用.北京:测绘出版社,1992:15.]

AMODISdataderivedspatialdistributionofhigh-,medium-

andlow-yieldcroplandinChina

JIYongzan1,2,YANHuimin1,LIUJiyuan1,KUANGWenhui1,HUYunfeng1

(1.InstituteofGeographicSciencesandNaturalResourcesResearch,CAS,Beijing100101,China;

2.UniversityofChineseAcademyofSciences,Beijing100049,China)

Abstract:WiththepopulationincreaseandfoodconsumptionupgradeinChina,theissueoffoodself-sufficiencyattractsmuchattentionfrombothChinesegovernmentandinternationalsociety.Underthecircumstancesofinadequatecroplandresourcessupplyandplentyofcroplandoccupiedbyurbanconstruction,improvingtheutilizationefficiencyofarablelandresourcesandincreasingthecroplandproductivityhavebecomethefundamentalstrategiesofagriculturaldevelopmentinChina.Since1988,croplandimprovementprojects(medium-yieldandlow-yieldfieldimprovementandhigh-standardcroplandconstruction)havebeenlaunchedandimplementedonalargescale,thereforeaspatiallyexplicitmapforthedistributionofhigh-,medium-andlow-yieldcroplandwasessentialforcroplandimprovementplanning.Inthisstudy,anewmethodforrecognizinghigh-,medium-andlow-yieldfieldisdevelopedbasedoncroplandproductivity,whichiscalculatedbyusingalightuseefficiencymodelandMODISdatawitha500-mresolution.Thismethodcannotonlyreflecttheregionalheterogeneityofcroplandcondition,butalsoexpressthespatialdifferencesonagridscale.Atthesametime,iteffectivelyovercomestheshortageofstatisticaldatabasedmethodinacountyunit.Theresultsshowthattheproportionofhigh-yield,medium-yieldandlow-yieldcroplandinChinais20.66%,39.56%and39.78%,respectively.About3/4oflow-yieldcroplandislocatedinthehillyandmountainousregions,while53%ofthehigh-yieldcroplandislocatedinplainarea.Thefiveprovinceswiththelargestareaofhigh-yieldcroplandareHenan,Shandong,Jiangsu,HebeiandAnhui,whicharealllocatedintheHuang-Huai-Hairegion.Thesumofthehigh-yieldcroplandareainthesefiveprovincesaccountsfor41.75%ofthenationaltotalhigh-yieldcroplandarea.InHeilongjiangprovince,SichuanprovinceandInnerMongoliaautonomousregion,wherethecroplandarearanksthetopthreeofChina,theproportionofhigh-yieldcroplandareaineachprovinceonlyaccountsfornotmorethan15%.

Keywords:foodsecurity;lightuseefficiencymodel;croplandproductivity;high-yield,medium-yieldandlow-yieldcropland