1 引言

农作物种植结构的时空变化信息是农业地理和农业可持续发展的重要研究内容，特别是在研究农业生态系统对陆地碳循环贡献，评价全球变化对区域农业生产影响，分析农作物空间格局动态变化机制和建立模拟模型等方面可作为基础信息，对国家农业生产管理、国家粮食安全预警和评估等理论与实践问题有重要意义。

农作物种植结构变化以各种农作物种植面积变化为基础数据，以作物播种面积占总播种面积的比例表征，其历史变化规律可为调整粮食种植结构和指导农业适应气候变化提供科学规划、合理布局的依据。

农作物种植结构的时空分布受自然资源条件和市场经济行为影响，具有空间集聚性和时序动态性等特征，因此及时、有效地获取农作物种植结构的时空信息尤为重要。

目前获取农作物种植结构信息的数据源有两类，一是基于国家统计部门的逐级统计上报和抽样估测数据集，另外则是运用遥感技术手段进行作物面积调查。遥感技术因其时效性和空间精度高，越来越多的中、高分辨率遥感数据地被应用于作物种植结构快速调查中，可及时准确地获取全国农业种植结构信息，遥感获取种植结构技术的进步，为实时监测农业种植结构的变化提供了有效手段，然而遥感监测的种植结构尚难在宏观层面进行长时序、大尺度上的变化分析。

尽管基于行政区的统计数据在时效性和空间分辨率方面较差，但统计数据在区域农业生产研究和大尺度时空变化分析方面具有显著的优势，时间序列长和数据高度综合对于掌握国家尺度上的农作物种植结构时空变化规律具有不可替代性，是分析历史时期农业种植结构变化规律的基础数据。

改革开放前，中国的粮食种植结构是长期执行计划经济体制以及错综复杂的国际国内政治经济形势综合作用的结果；改革开放后农作物的种植结构随着经济发展进行了调整和优化，但改革开放以来中国的农作物种植结构变化特征尚缺少研究。

经济的快速发展导致中国第一产业的比重不断降低，与之同期调整的还有农作物种植结构。中国农业发展虽然取得了巨大成就，但进一步的发展也面临一系列的农业结构调整问题，特别是农业总产量与结构不合理，部分农业种植区面临环境污染问题，作物缺少连茬导致病虫害加重，对中国新一轮的农作物结构调整都带来不小压力，新一轮的农作物结构调整亟待总结过去30 年的种植结构变化规律。

因此，本研究旨在基于全国农作物统计数据库，分析农作物种植结构变化规律，为国家宏观层面的种植结构调整和粮食安全保障提供科学基础。

2 数据和方法

2.1 数据源

（1） 县级农业统计数据库。

整理农业部种植司发布的农作物数据库，获得了1980-2011 年逐年的农作物总产量、面积和单产数据库。由于历史时期行政边界多次变更，研究以2001 年国家基础地理信息中心发布的县级行政边界为基础，依据历史时期行政边界变化，将全国的农作物数据库的统计单元整理合并为2341 个县，不包含台湾、香港和澳门地区。形成了一套涵盖32 年，2341 个县的长时序空间数据集。作物种类涵盖38 种，最终按照农业统计数据库的统计要求合并为11 类：水稻、小麦、玉米、大豆、薯类、油类、麻类、糖料、棉花、蔬菜、水果。

（2） 种植结构类型。

种植结构类型组合采用作物种植面积占所有作物类型种植面积的百分比组合确定，种植结构类型组合可分两种情况：

① 采用作物种植面积占所有作物类型种植面积的百分比超过30%和前三位的组合确定，一般组合不会超过3 种作。例如：某县级行政单元的作物种植面积中前三位只有水稻的种植面积占所有作物种植面积的比例超过30%，则认为该县的种植结构类型为单一水稻型；以此类推，有2 类或3 类作物超过30%，则以超过的作物进行组合，如水稻—玉米型，水稻—小麦—大豆型。

② 所有作物比例都不超过30%时，以前三位作物为组合，例如某县的水稻种植面积为25%，小麦为24%，大豆为19%，玉米为17%，则可将该县的种植结构类型组合为水稻—小麦—大豆型。作物类型的组合考虑作物类型占比和大小排序，前后两者情况的3 种作物组合可能存在比例上差距，但排序上如果一致可认为是同一种作物种植结构组合类型。

2.2 研究方法

2.2.1 种植结构类型变化趋势为考察全国尺度上的作物种植结构类型演变趋势，定义了一个种植类型组合丰富度公式，其含义表示的是某一年份的种植结构类型数量与研究时段内全部出现的类型比例，公式表示如下：

（1）式中：Rt表示的是某一年种植结构类型的比例，其值域为(0, 1]，值越大则种植结构类型越丰富；mt是是某一年种植结构类型的种类数量；mmax是研究时段内所有种植结构组合类型的种类数量。

运用ArcGIS 10.0 软件平台将种植结构类型的空间分布作图展示，并统计归纳种植结构的空间分布特征。根据时序上的种植结构类型，分析其变化。

2.2.2 种植结构变化趋势分析为研究近30 年来中国农业种植结构变化的状况以及该结构的变化速率在空间上的差异，对全国2341 个县级行政单元内的某种作物对应的比例按年份进行一元回归分析，得到斜率系数（S），用来表示某种作物的比例状况的变化趋势。斜率为负值，表示该县的某种作物呈减少趋势；斜率为正值区表示该县的某种作物呈增加趋势，如果通过0.05 的显著性水平（p < 0.05），则认为某种作物变化比例增加或减小趋势显著。斜率计算公式如下：

（2）式中：S 表示某类作物种植比例的变化趋势线性斜率；t 为研究时段内的年份；n 为研究时段内的年份数；Xtj表示的是某年j 作物的种植比例。为分析30 年来作物类型的种植比例变化，将S的空间分布用图展示出。

2.2.3 种植结构变化的空间聚集采用局部空间自相关分析（LISA） 方法测度农业种植结构是否存在变化趋势高度集聚。局部Moran's I 指数公式：

（3）式中：LISAi表示的是种植比例变化趋势的局部自相关指数；s 为某类作物种植比例的变化趋势线性斜率的标准差；xi为i 省某类作物种植比例的变化趋势线性斜率；xj为j 某类作物种植比例的变化趋势线性斜率；x 为某类作物种植比例的变化趋势线性斜率的平均值；wij为空间权重，表示省级单元i 与j 的临近关系，这里根据邻接标准来度量；当i 省与j 省相邻时，权重wij取1，否则取0。高高（HH） 集聚的显著性LISA图可表示农业种植结构变化趋势在县域层面上的连片分布。

3 结果分析

3.1 种植结构类型的年际变化

1980-2011 年，按照农作物种植面积占30%或者低于30%的按作物种植面积的前三位组合的标准，中国的种植结构类型数为182 种。2000 年以前，国家尺度的种植结构类型数量介于21~54 种类之间，类型丰富度在0.18；2002 年后，种植结构类型组合数量介于84~106 之间，类型丰富度指数增长至0.52；种植结构类型的调整表明2000-2002 年期间，中国种植结构发生了重大调整，类型丰富度显著增加，更为多样的种植结构逐渐替代简单的种植结构（图1）。

从调整类型来看，主要是单一型种植结构被组合型种植结构替代，从以粮食作物种植结构为主的种植格局转变为粮食与经济作物组合的格局。

中国的主要种植结构类型有单一水稻型、水稻—小麦组合型、水稻—玉米组合型、单一小麦型、小麦—玉米组合型、单一玉米型、单一蔬菜型、单一水果型等等。图1b 列举了中国三大作物及其典型的组合表征的农作物种植结构类型的时序变化趋势。

三大粮食作物的占优单一型种植结构呈逐年递减趋势，单一水稻型的县数比例由最高值37.5%（1985 年） 降至19.0%（2008 年），单一小麦型由23.3%（1985 年） 降至8.8%（2011 年），单一玉米型相对稳定保持在10%~14%左右；小麦—玉米的组合型在2002 年前维持在15%~20%之间，而2002 年之后降到10%以下；水稻—小麦型也有类似趋势，2002 年前在6%左右，2002 年后降为2%以下；水稻—玉米型在2002 年前在5%左右，而之后低于1%。

图1c 展示的是经济作物的单一种植结构型的时序变化趋势。单一蔬菜型、单一薯类型、单一水果型和单一油料型在2002年前比例非常低，2002年后基本都有超过2%。

3.2 种植结构类型的区域分异

表1 表示的是种植结构比例数据变化，1980-2011 年中国出现的种植面积前10 位的主要种植结构类有16 种。

1980 年时主要的种植结构类型以粮食作物为主，全国82.7%的县级农业种植结构是水稻、小麦、玉米及其相互的组合型号，前10 位的种植结构类型占全国总县数的93.3%；

1990 年，前10 位的比例略升至95.0%，相比1990 年，单一蔬菜型和单一水果型的种植县数替代单一大豆型和单一棉花型进入种植结构前10 位；

2000 年，全国的种植结构尚未发生根本性变化，但前10 位县数已降低了近10%，即有209 个县通过调整种植作物类型，改变了种植结构比例；

2002 年，是中国农业种植结构的重大转折时期，以粮食作物为主的种植结构发生了根本性变化，前10 位中已种植粮食作物为主要种植类型的县数降至50.7%，单一蔬菜型县数上升至第4 位，占9.1%；

2011 年，中国前10位的县数总量为1699 个，只占72.6%，相比1990 年减少了524 个县（22.4%），粮食作物与经济作物的混合结构出现，例如水稻—蔬菜组合型。农作物种植结构的调整，不仅是大政策的体现。

2002 年以后，中国的种植结构趋于丰富，既与国家农业政策的调整有关，也与农户的微观选择有关，例如在大城市周边地区形成的粮食与蔬菜作物组合型或是粮食与水果作物组合型，就是城市化带来的经济作物需求。

图2 展示的是5 个时期中国县级行政单元内的农作物种植结构类型的空间分布及其变化。

从种植结构类型变化上看，1980-2011 年中国63.8%的县市（1494 个） 调整了种植结构，主要改变方向为降低主要粮食作物产量的种植比例，增加经济作物和蔬菜瓜果类作物的种植比例；空间分布来看，改变了种植结构类型的县主要分布于中国沿海和西部地区，未改变种植结构的县主要分布于中国中部地区（图2f）。

只有近1/3 的县级未改变农业种植结构，其中单一水稻型385 个，单一小麦型105 个，单一玉米型160 个，小麦—玉米组合型130 个，尽管如此这些县中的主导作物类型比例还是有不同程度的降低趋势。

1980-2011 年间，虽然以单一水稻型为典型特征的种植结构县数比例在持续下降，水稻依然是第一大种植作物；空间分布上，1980 年中国单一水稻型作物主要分布在中国南方的11 个省市，然而到2011 年分布省份无变化，但各省内都有不少县调整了单一水稻型的种植结构，以长三角、珠三角和闽三角的城市化地区最为典型。

单一玉米型则主要分布在东北至西南的东三省、内蒙古东部至广西等地，从1980-2011 年在同一线上有所外扩。

1980 年单一小麦型则主要分布在山东新疆北部、甘肃、宁夏、内蒙中部、河南、安徽北部和山东南部等地，2011 年则相对萎缩至河南、安徽北部和山东南部局部区域。

1980 年小麦—玉米组合型主要分布在山东中北部、河北南部、河南北部、山西、陕北和新疆中南部、而到2011年则萎缩至山东中北部、河北南部、河南北部、山西南部等地。

3.3 作物种植比例变化的分布

图3 和表2 表示的是1980-2011 年中国农作物种植比例的年际变化趋势。全国县级尺度上，三大粮食作物所占的农作物种植比例呈减少趋势，47%的水稻、61%的小麦和29.6%的玉米种植县的种植比例呈显著减少趋势（p < 0.05）。

结合种植结构类型变化看（图2），粮食作物由水稻为主的格局调整为水稻、小麦和玉米的种植区域共存的格局；其他作物的种植比例呈现空间上的显著增加。

水稻种植比例显著减少的地区主要集中在福建、广东和浙江等沿海省份，改变了其原始的种植结构类型；传统的水稻种植区的江西、湖南、安徽南部、四川盆地和重庆等地也有减少，但减少程度相对于沿海较轻，该地区主要的农业种植结构调整有别于沿海省份，基本上不改变种植结构类型，但降低种植比例。

小麦种植比例显著减少的地区主要集中在黄土高原地区包括甘肃南部、宁夏、陕北高原至山西一带，这些地区的是种植结构类型也发生了显著调整；此外华北平原的传统小麦种植区的种植比例也显著减少，但未改变种植结构类型。

玉米种植比例显著减少的地区在空间上覆盖中国的第二阶梯，在高原区形成北东—西南走向，从东北三省至西南部云贵高原的条带地区，可以称之为种植比例“玉米减少带”（图3b）；其中种植结构类型变化的地区主要分布在云南北部、四川南部、陕西南部、湖北西部、京津冀地区和东北平原区；而种植比例减少包含约650 个县级行政单位，覆盖从黑龙江北部至云南西南部的15个省市。

与之相应，经济作物如大豆、蔬菜、薯类、水果和油料作物则呈增加趋势，50.3%的大豆、87.7%的蔬菜、59.5%的薯类、70.7%的水果和69.5%的油料种植县的种植比例增加趋势显著（p < 0.05）（图2）。

大豆的种植比例增加区主要分布在中国东北的黑龙江和吉林省南部的少量县，而在中国其他地区都有不同程度的增加；蔬菜种植比例显著增加的地区主要集中在沿海地区，包括京津冀、长三角和珠三角等大城市周边，与城市化需求密切相关，还有一部分分布于山东半岛，山东与江苏交界，广西东北部等县市；薯类种植比例显著增加的地区主要分布于内蒙中部、宁夏南部、甘肃西南以及西南五省市的交界山区；水果种植比例显著增加的地区主要分布于山东烟台、北京周边、陕北黄土高原、新疆和东南沿海的广东、福建等地，油料种植比例显著增加的区域集中分布于中国的西藏南部、新疆北部，青海东部以及长江中下游的湖北至安徽沿岸县市。

此外，棉花、糖料和麻类是区域性种植作物，全国棉花有42.2%的棉花，糖料有55.6%的糖料以及24.2%的麻类种植县存在种植比例变化趋势，其中30.7%的棉花和24.4%的糖料呈增加趋势，13.6%的麻类呈减少趋势，三类作物的分布呈现区域性，棉花的显著增加区主要分布于新疆，华北平原和长江中下游平原；麻类主要分布于长江中下游沿岸；而糖料则零星分布于云南、广西、广东和福建等地。

3.4 城市地区的种植结构变化

为了解城市化对区域性种植结构变化的影响，将中国的城市地区按照市辖区的范围提取了660 个城市，共717 个县区的农业种植结构变化趋势，分析城市地区种植结构变化趋势。

总体趋势表明，城市及其周边地区的种植结构变化表现出两个特征，

一是农作物种植面积在过去30 年间迅速萎缩，耕地大量转为城市建设用地，粮食作物用地和林地改为经济作物用地，种植结构显著改变，90%的县区范围内农作物种植面积减少，粮食作物面积迅速减少，而水果和蔬菜类种植面积迅速增加；

二是种植结构由早期的粮食作物转向蔬菜，水果为主经济作物；68.6%的县区发生了种植结构变化，以经济作物中的单一蔬菜型和单一水果型为主的单一型种植结构随着城市扩展而迅速增加，1980 年只有3 个县区以蔬菜型种结构植为主，2011 年则有54 个县区级形状单位以蔬菜为主；1980 年只有13 个以水果为主要种植类型的县区，2011 年为29 个。

224 个县市未发生种植结构类型变化，都是粮食作物为主，112 个是单一水稻型，3 个水稻—玉米型；6 个水稻—小麦型；20个小麦型，48个小麦—玉米型；31个玉米型。

3.5 种植结构变化的空间集聚

运用局部自相关方法，对11 种作物进行空间集聚分析，获得了1980-2011 年中国种植结构的空间集聚特征（图4），1300 个县级行政单位具有作物的空间集聚效应；其中，水稻的高度增长聚集区占全国县级行政单位的2.86%，主要分布在中国东北地区的吉林和黑龙江的松嫩平原和黑龙江的三江平原地区；

小麦的高度集聚区占5.64%，主要分布在东北的辽宁和吉林中部地区，此外在广西地区有显著的连片集聚区。

玉米种植结构增长的连片区占6.11%，分布在中国北方的甘肃黑河流域和宁夏，内蒙高原中部附近；大豆的4.53%连片集聚在内蒙古东北，黑龙江和吉林的东南等地；1.62%的麻类集聚区主要分布在湖南，江西和湖北的交界区；棉花的高级聚集区占7.77%，主要分布在新疆全境，山东北部、天津和河北东南以及河南东部部分县市；8.24%的蔬菜类聚集区主要分布在京津冀地区、上海、浙江、福建的沿海县市和珠江三角洲的外围城市如肇庆、清远、从化和惠州等地；12%的薯类种植比例增加的集聚区主要分布在中国的西部省份，包括从内蒙中部，宁夏南部，甘肃南部、四川山区、重庆周边及其与之接壤的贵州、云南等地；10%的水果种植比例集聚区主要分布在北京周边、陕北黄土高原、新疆中南部、江西赣南、福建南部、广州的粤东粤西两轴和广西南部；1.41%的糖料主要分布在贵州南部、广西南部和广东的雷州半岛附近；9.35%的油料作物主要分布在青海东部，西藏南部，湖北东部、湖南北部、河南南部、安徽南部和江西北部的集中区域。

4 结论与讨论

4.1 结论

（1） 国家尺度上的种植结构类型变化，中国的种植结构调整从2002 年起发生重大变化，类型丰富度显著增加，多元种植结构逐渐替代单一型种植结构。三大粮食作物的占优单一种植结构类型呈逐年递减趋势，其中既有三大粮食作物的组合型结构混合因素，又有三大粮食作物种植面积转换为其他作物种植面积的变化因素。

（2） 近30 年来中国出现的种植面积前10 位的主要种植结构类型有16 种。1980 年全国82.7%的县级农业种植结构是水稻、小麦、玉米及其相互的组合型；但从2002 年起，中国农业种植结构进入调整期，以粮食作物为主的种植结构发生了根本性变化，前10 位中以种植粮食作物为主要种植类型的县数降至50.7%。

（3） 全国县级尺度上，三大粮食作物所占的农作物种植比例呈减少趋势，47%的水稻、61%的小麦和29.6%的玉米种植县的种植比例减少趋势显著（p < 0.05）。粮食作物以水稻为主的种植格局调整为水稻、小麦和玉米的种植区域共存的格局，其中玉米种植面积比例在空间上变化最为显著，在中国形成东北—西南向的“玉米减少带”，同时该区域也是中国种植结构多样性最丰富的地带。作为种植结构调整最为热点的城市地区，城市化对种植结构变化影响显著，水果和蔬菜类种植比例在城市化地区快速增加。

（4） 种植结构变化趋势的空间集聚区域。近30 年来中国种植结构在1300 个县级行政单位形成了的空间集聚效应；其中，水稻的高度增长聚集区占全国县级行政单位的2.86%；小麦占5.64%；玉米种植结构增长的连片区占6.11%；4.53%的大豆；1.62%的麻类；棉花的高度聚集区占7.77%；8.24%的蔬菜类聚集区；12%的薯类；10%的水果种植比例；1.41%的糖料；9.35%的油料作物。

4.2 讨论

（1） 与种植结构密切相关的宏观决策问题是中国粮食安全的保障程度。《国家粮食安全中长期规划纲要（2008-2020 年）》设定的2020 年粮食自给率为95%，其中谷物为100%，这意味着农作物种植结构变化与粮食供给能力密切相关，中国种植结构变化对粮食安全的影响将涉及到以下几个方面：

① 受粮食价格的影响，农户趋向于种植非粮食作物，或者撂荒土地；而主粮作物作为肉类食品的低端原料，受中国人口增长和经济增长带来的肉类食品需求影响，将原有的主粮转作饲料投放，因此未来受粮食消费习惯改变、国内外粮食价格变化、区域水土资源条件等影响，依托调整作物种植结构实现粮食增产的潜力有限，进而影响粮食安全；

② 为减少温室气体排放，有效保护生态环境，实现可持续发展，越来越多的粮食作物种植区将会因生物质能源需求上升而转种生物质能料作物，这将直接占用种植农作物的耕地资源，进而增加粮食供给保障的负荷，因此需要统筹粮食安全保障与生物质燃料种植对种植结构的调整；

③ 随着城市化发展，越来越多城市周边的农业种植以生态观光农业需求为主，同时因为物流技术的发展，增加对远程经济作物和水果蔬菜类作物的需求，而改变城市近郊和传统农业主产区的种植结构，而被置换的农业用地也可能被城市建设用地侵占，进而减少粮食供应能力。

（2） 种植结构变化的驱动因素分析。一般研究认为，种植结构变化的驱动是政策、技术进步、社会需求、经济效益、自然条件等多因素综合的作用过程[19]。也有研究选择单一作物水稻的种植面积变化与驱动因素进行定量分析表明，除了上述原因外，城市化和气候变化也是驱动种植结构调整的重要因素。

然而，本研究中种植结构变化并不是线性趋势，而是在某一阶段迅速发生突变，因而种植结构演变的驱动因素尚难以定量化，就定性的描述来说，中国种植结构调整的几大重要因素包括：

① 城市化水平提升改变了城市周边地区的种植类型，将原有粮食作物改种为水果、蔬菜类经济作物，另外直接减少耕地面积，减少种植面积；

② 不同作物的价格差异驱动种植者将经济价值低的作物改种为价值高的作物；

③ 技术进步包含两个方面，一是农业技术进步，提升了作物产量和种植面积，另外则是物流成本降低，加大了不同种植区间的作物需求交换；

④ 近30年中国农业投入的加大，特别是单位面积的农业投入提升大大增加了总的作物产量，相应地就减少了作物种植面积的需求；

⑤ 肉类消耗，增加了饲养作物牲畜所需的主粮作物的供应；

⑥ 气候变化的影响，已有研究表明近30 年中国的气候变化对作物产量和种植面

积都有温和的正效应；而在过去的50 年中，由于气候变暖造成了全国种植制度改变玉米、冬小麦和双季稻种植北界北移，也为结构变化提供了地理空间。

为此，种植结构变化的影响及驱动因素的时空分异特征将是下一步研究的重点。

参考文献略，具体可参阅原文。