土地作为人地关系地域系统中的关键要素，深刻影响着社会经济的转型和发展。同时，社会经济的转型发展也带来了土地利用方式的变革，进而导致区域土地利用转型^[1]。土地利用转型指在特定的社会经济变革过程中，区域土地利用由一种形态（含显性形态和隐性形态）转变为另一种形态的过程^[2]。其中土地利用隐性形态主要涉及土地质量、产权、经营方式、投入和产出等多重属性^[3]。

土地利用强度是土地利用隐性形态的重要属性，既可以反映出人类对土地的干预程度，又可以呈现出特定时期内区域社会经济发展状态。土地利用强度的变化是土地利用形态变化的内在动力之一^[4]，因此，基于土地利用强度来研究土地利用隐性形态的变化是一种有效的尝试。当前，针对土地利用强度的内涵和外延^[5-6]、定量测度、格局演化^[7]及动力机制^[8]研究较多，而针对区域土地利用强度的变化与粮食生产的关系研究较少。基于以上分析，本文选取我国粮食主产区（黄淮海地区）与粮食生产密切相关的耕地和村庄建设用地^[9]利用强度的变化，表征区域土地利用隐性形态的转型，在土地利用隐性形态定量测度的基础上，划分土地利用隐性形态变化与粮食产量的耦合类型，在实证研究的基础上探讨基于土地利用强度变化的区域土地利用转型与粮食生产的耦合关系，用以指导区域粮食生产系统的优化与调控。

黄淮海地区占我国国土面积的4.3%，2010年黄淮海地区耕地面积、人口数量和粮食产量占全国的比重分别为20.6%、22.5%和30.8%，其中小麦和玉米产量占全国总产量的比重约为70%和40%^[10]。依照地形特征和农业生产特点该地区又可以划分4个亚区，分别为燕山太行山山麓平原区（以下简称山麓平原区）、冀鲁豫低洼平原区（以下简称低洼平原区）、黄淮平原区和山东丘陵区^[11]（图 1）。黄淮海地区作为我国耕地面积广、人口密度大和粮食产量高的“三高”地区特征日渐明显^[12]。同时，该地区又是我国城乡快速转型发展地区，耕地、村庄建设用地等与粮食生产和农民生活密切相关的土地利用形态快速转型。

图 1 黄淮海地区概况 Figure 1　Region profile of Huang-Huai-Hai Plain

图选项

本文研究过程中使用的分县统计数据来源于《中国县（市）社会经济统计年鉴》及各地（市）年鉴，统计数据主要包括化肥使用折纯量、有效灌溉面积、农作物播种面积、农业机械总动力、乡村人口数量、乡村从业人员数量、农业劳动力数量等指标，部分评价指标市辖区没有统计数据将使用地级市平均数补充。土地利用数据来源于国土资源部，该数据采用第二次全国土地利用调查分类体系，数据根据Landsat TM和ETM传感器经几何校正后合成的真彩色影像监督分类获得。经过市辖区合并处理后，本文研究的黄淮海地区共有358个县级行政单元。

本文以土地利用强度的变化表征区域土地利用隐性形态的转型。研究表明，耕地与村庄建设用地的变化是乡村地区土地利用转型的主体^[9]，因而，本文构建了耕地利用强度和村庄建设用地利用强度，用以定量刻画黄淮海地区与粮食生产密切相关的土地利用隐性形态转型。耕地利用强度（LUI_F）评价指标体系的选择以粮食生产相关的“投入”指标作为评价依据。黄淮海地区粮食生产系统中地均化肥投入（P_fer）^[10-14]、灌溉系统（P_irri）^[13]、复种指数C_ropi和地均农业机械（P_am）的使用作为评价耕地利用强度的指标。而影响粮食生产的其他要素如自然环境和资本等要素，未纳入测度指标。村庄建设用地指乡村地区与农民生产和生活相关的建设用地。城镇化进程中，农业劳动力非农就业比例不断上升，村庄建设用地利用格局不断演化，村庄建设用地在粮食生产中的功能也发生了较大变化。村庄建设用地利用强度（LUI_R）的评价指标亦考虑到与粮食生产密切相关的指标，本文以乡村人口密度（R_popd）、农业劳动力比重（P_arl）、村庄建设用地比例（P_rc）和乡村人均建设用地（P_ccl）作为评价指标。评价指标中除乡村人均建设用地面积对土地利用强度的作用为负效应外，其他指标的作用均为正效应。评价指标权重的确定参考了已有的研究成果^[15]，在咨询国内相关领域专家意见的基础上，利用熵值法^[16]确定，即考虑了数据样本客观性的同时兼顾专家判定，各评价指标的权重如表 1。

表 1　土地利用隐性形态评价指标体系 Table 1　Evaluation index system of land use recessive morphology

表选项

为了消除评价指标体系源数据量纲对评价的影响，对源数据利用极值法进行归一化处理：功效为正的指标使用公式（1），功效为负的指标使用公式（2）：

(1)

(2)

式中：T_n^k为第n个样本的k指标归一化后的数值，X_n^k为第n个样本k指标的原始数据，X_max和X_min分别代表第k个指标的最大值和最小值。

在数据标准化的基础上利用指标权重与归一化之后的数据相乘求和的方法定量评价土地利用隐性形态转型，耕地利用强度的计算利用公式（3），村庄建设用地利用强度计算利用公式（4）：

(3)

(4)

式中：m_j和n_p分别为二者各评价指标的权重，T_ij和T_ip为i县第j和第p个指标归一化后的数值。

利用耕地利用强度和村庄建设用地利用强度，结合县域粮食产量，基于ArcGIS10.2中自然间断点（Jenks）分类方法，分别把3个序列数据由低到高分为“低”（1）、“中”（2）和“高”（3）3组，利用3伊3矩阵交叉方法，构建“低-低-低”（111）…“高-高-高”（333）共27类耕地利用强度-村庄建设用地利用强度-粮食产量耦合（FRGC）类型（图 2）。在研究过程中利用信息图谱^[17]的方式，实现在复杂耦合类型模式下的格局演化研究。

图 2 耕地利用强度-村庄建设用地利用强度-粮食产量耦合（FRGC）分类体系 Figure 2　Coupling classification system of farmland intensity-rural construction land intensity-grain yield(FRGC)

图选项

研究结果表明，耕地利用强度持续增加，村庄建设用地利用强度加速下降，粮食产量波动上升。1990—2010年 LUI_F（图 3a~图 3c）均值由0.29上升到0.4，山麓平原区耕地利用强度一直较高；黄淮平原区和低洼平原区增长较快，成为耕地高投入区；山东丘陵区耕地利用强度较低。期间，LUI_R均值由0.52下降到0.45，并呈现加速减少的趋势（图 3d~图 3f），1990—2000年村庄建设用地利用强度呈现下降趋势，山东丘陵区保持较高水平，黄淮平原西部地区LUI_R大部在增加，而低洼平原区和苏北地区及山麓平原区村庄建设用地利用强度处于下降趋势。1990—2010年黄淮海平原地区粮食产量由1.06亿t上升到1.65亿t，1990年、2000年和2010年粮食产量占全国总产量的比重分别为24.4%、24%和30.1%，在经历了2000年附近短暂的下降后，黄淮海地区粮食生产在全国的地位不断加强。黄淮平原区和低洼平原区是粮食增产的主要地区（图 3g~图 3i），而山麓平原区和山东丘陵区粮食产量在下降。

图 3 1990—2010年LUIF、LUIR和粮食产量时空格局演化 Figure 3　Spatio-temporal pattern of LUIF, LUIR and grain yield during 1990—2010

图选项

1990—2010年耕地利用强度-村庄建设用地利用强度-粮食产量耦合类型特征呈现稳中有变的时空格局。1990—2010年LUI_F较低（111~133）的区域主要集中于山东丘陵区、低洼平原东部和黄淮平原南部地区（图 4a~图 4c），LUI_R和粮食产量较低的区域主要集中在环渤海湾地区；山东丘陵区和山麓平原部分地区LUI_R和粮食产量较高。耦合类型中LUI_F居中（211~233）的区域主要分布在低洼平原北部和山东丘陵东部地区，并在黄淮平原地区离散分布（图 4d~图 4f）；这些县域中LUI_R和粮食产量较低的区域主要集中在低洼平原北部地区，而LUI_R和粮食产量较高的地区主要集中在黄淮平原区。黄淮海地区LUI_F较高（311~333）的地区主要集中在山麓平原大部、低洼平原南部及苏北和豫北地区，并逐渐在冀鲁豫交界地区形成明显的集聚区（图 4g~图 4i），而山麓平原北部地区逐渐退出该耦合类型；FRGC类型中LUI_F、LUI_R和粮食产量均较高的地区逐渐向冀鲁豫交界地区集中，而苏北地区和山麓平原地区此类县域逐渐减少，石家庄市附近地区耦合类型逐渐由“333”类型为主演化为“322”类型为主。

图 4 1990—2010年耕地利用强度-村庄建设用地利用强度-粮食产量耦合（FRGC）时空演化过程 Figure 4　Spatio-temporal evolution process of LUIF, LUIR and GY(FRGC)during 1990—2010

图选项

土地利用转型背景下，不同地类土地利用强度相互影响，共同推动区域粮食生产格局的转型。转型初期，耕地和村庄建设用地利用强度均增加，推动粮食生产能力提升。城镇化背景下，农业剩余劳动力析出，农业劳动力数量相应减少，就业结构和收入结构不断改善，农户农业生产的资本性投入能力增强，耕地利用强度提高^[18]。农民改善居住条件的需求增加，农村宅基地开始向外扩张，村庄建设用地的比例增加，这一阶段，农业生产收入比例仍占主要地位，农民的兼业情况主要呈现为季节性和短距离外出务工或本地非农就业为主。乡村人口总体增长状态下，农业劳动力比重下降不明显，村庄建设用地利用强度增加。这一阶段耕地利用强度与村庄建设用地利用强度普遍上升，粮食生产能力得到提升，粮食安全格局稳定（图 5a）。

图 5 土地利用转型与粮食生产耦合关系 Figure 5　Coupling relationship modes between land use transitions and grain production

图选项

转型中期，耕地利用强度继续增加和村庄建设用地利用强度开始下降，区域粮食安全保障能力下降。农业劳动力转移过程中，由于人口增加和经济因素导致的村庄建设用地面积呈现快速的外延扩张，人均建设用地面积迅速提高（图 5b）。同时，由于大量人口外流，村庄宅基地出现废弃和空置现象，农村“空心化”^[19]问题突出。此外，部分农业劳动力实现资本原始积累，彻底脱离农业生产，农村宅基地需求下降。农业劳动力的兼业化比例提高，农业劳动力数量快速下降，农业劳动力比重下降，乡村人口密度也随之下降。城镇化快速推进过程中，城镇建设用地和工矿用地需求增加，占用耕地和征用村庄用地的比例提高，村庄建设用地在建设用地中的比例不断下降。以上多种因素共同推动了村庄建设用地利用强度的下降。随着农户收入结构的改善，农户继续追加耕地投入的能力迅速提升，但是纯农户的比例不断减少^[20]，并且农户的年龄结构趋于“老弱化”，大规模增加资本性投入需求下降，耕地投入主要以省工性投入为主，地均化肥使用量不断增加。农户耕地流转的趋势开始出现，且比例不断提升。而耕作条件较好的平原农区随着农业机械化的大力推广，粮食播种面积不降反升（如研究时段内的黄淮平原区），灌溉比例也相应提升，而丘陵及耕地质量不高地区的耕地利用形态会出现弃耕或撂荒^[21]的现象，耕地利用强度区域格局差异明显。这一阶段耕地利用强度仍保持较高水平，粮食产量较稳定。但同时也应该看到粮食生产要素投入产出的效率不断下降，粮食生产成本不断增加，粮食生产所带来的资源环境问题不断凸显^[22]。耕地资源被建设占用的比例迅速提高（图 6），村庄建设用地利用强度不断下降，粮食生产的保障能力降低。

图 6 1990—2010年黄淮海地区耕地面积与村庄建设用地面积变化 Figure 6　Area changes of farmland and rural construction land inHuang-Huai-Hai Plain during 1990—2010

图选项

转型后期，耕地和村庄建设用地利用强度均下降，粮食生产高效化和集约化保障区域粮食安全。城镇化发展中后期，农业劳动力比重低，村庄建设用地整治需求旺盛，农村聚落社区化，人口居住集中化开始大力推进，村庄建设用地的比例迅速下降，村庄建设用地将实现优化重组^[23]。乡村人口密度不断降低，农业劳动力数量和比重不断下降，村庄建设用地利用强度下降。耕地利用方面由于乡村人口逐渐减少，土地流转的比例不断增加，粮食生产的规模化和集约化不断加强，可持续的农业生产系统^[24]不断优化，耕地利用强度逐步下降，并趋于稳定。农村剩余劳动力转移基本完成，新型农业经营主体成为粮食生产的主力，粮食生产格局进入新的状态（图 5c）。

城镇化进程中，我国土地利用形态与粮食生产均发生了深刻的变化，系统研究二者耦合关系、动力机制及其演化规律，对于揭示我国粮食生产系统要素转型、结构变革和功能演化，保障我国粮食安全的同时，促进乡村良性发展，优化乡村土地利用结构，实现社会良性转型具有至关重要的作用^[25]。因此，粮食生产转型过程中，应合理优化土地利用转型与粮食生产的耦合关系，加强耕地的保护力度，限制村庄建设用地的无序扩张，合理搭配乡村农业生产体系中的人口及其他生产要素，确保农户生产与生活水平稳定提升。

（1）城镇化进程中，黄淮海地区土地利用形态和粮食生产格局迅速转型。1990—2010年土地利用转型格局呈现出耕地利用强度不断增加，村庄建设用地利用强度加速减少，且时空格局区域差异明显。区域粮食生产能力经历短暂的下降后持续增强。

（2）1990—2010年耕地利用强度-村庄建设用地利用强度-粮食产量耦合类型时空格局特征明显。“低-低-低”（111）类型主要分布在环渤海湾地区，“高-高-高”（333）类型主要集聚于冀鲁豫交界地区。而“中-中-中”类型主要分布在山东半岛东部和低洼平原北部地区且不断集中。

（3）土地利用转型推动粮食生产格局不断演化。转型初期，耕地和村庄建设用地利用强度均增加，推动粮食生产能力提升。转型中期，耕地利用强度继续增加和村庄建设用地利用强度开始下降，区域粮食安全保障能力下降。转型后期，耕地和村庄建设用地利用强度均下降，粮食生产高效化和集约化保障区域粮食安全。依据土地利用转型与粮食生产不同阶段的耦合关系，制定合理的粮食生产政策能够有效保障区域粮食生产安全。

（4）本文基于土地利用强度变化初步探讨了土地利用转型与粮食生产的耦合关系，是土地利用转型对粮食生产系统影响的初步研究。而对土地利用隐性形态中涉及质量、产权和经营方式等土地利用形态尚缺乏深入的研究。本文初步构建了土地利用转型与粮食生产耦合关系模式，但对二者耦合过程的动力机制研究仍薄弱，这些内容是进一步研究的重点。