耕地是一个开放的“社会—经济—生态”系统，随着社会经济发展阶段的演变，呈现出功能多样化、利用集约化和管理模式复杂化的动态特征。这一演变不仅受自然资源、政策导向和市场需求的共同驱动，还反映了耕地资源配置效率与社会经济适应能力的提升。耕地利用隐性形态是耕地利用显性形态背后的深层次耕地属性，涵盖土地质量、投入强度、产出能力、功能属性、产权模式和经营方式等核心内容^[41]，是对耕地社会—经济—生态适应能力的综合反映。上述属性及内涵揭示了耕地利用隐性形态的特征，但主要停留在定性阐述层面，缺乏与定量测度间的有效衔接。

本文针对丘陵山地区耕地利用特性，从“要素—组织—功能”的思路出发，基于核心内容动态关联性，将耕地利用隐性形态内涵归纳为投入形态、组织适应性形态和产出形态三个维度(图 1)，以便更好地揭示耕地利用隐性形态的结构特征及内在规律，并进行定量测度。其中，投入形态指耕地利用系统中与自然资源发生相互作用的农业生产“要素”的数量、质量和结构，是耕地利用活动的基础条件，决定了系统运行的潜力；组织适应性形态指在不同自然和社会经济条件下耕地对多样化生产模式和经营方式的适应能力，是资源整合和管理优化的核心，体现系统论中的“组织”，强调其对资源整合的动态调节作用；产出形态体现投入要素与耕地资源结合后的价值实现，包括耕地利用的经济、社会、生计和生态价值等的数量、质量和结构，是对产出能力和土地功能的融合重构，体现系统论中的“功能”，反映耕地利用对多目标需求的响应。这一划分避免了隐性形态内容的碎片化，增强了维度之间的系统性和逻辑关联性，揭示了耕地利用隐性形态的系统特征，为量化分析隐性形态及其演变规律提供了清晰框架。

耕地利用隐性形态投入、组织适应性和产出3个维度并非独立存在，而是通过耦合作用形成动态调节的系统，其耦合关系主要体现为：

1) 投入与产出之间的联动机制是优化资源配置的基础。耕地投入的数量与质量决定产出效益，而产出的数量、质量与结构变化通过反馈作用优化投入结构与模式。合理的投入配置能够提升资源利用效率，促进经济、社会与生态效益的协同发展。

2) 投入与组织适应性之间的匹配机制是产出优化的关键支撑。充足的资源供给和合理的投入模式为组织适应性优化提供了资源条件和调整空间；组织适应性通过资源整合和管理优化，动态调整投入模式，使资源配置更具针对性和灵活性，提升投入要素的利用效率，增强耕地利用的区域适应性。

3) 产出与组织适应性之间的互促机制是提高耕地利用效能的关键驱动力。组织适应性通过优化耕地利用模式，促使产出结构更加合理，提高耕地利用整体效益；产出数量、质量与结构可反作用于组织适应性，推动管理模式优化，增强耕地系统对市场需求和环境变化的适应能力。

综上，投入、组织适应性和产出形态之间的耦合关系构成了耕地利用隐性形态的核心动力系统，其动态联动机制为耕地利用隐性形态的优化和可持续发展提供了重要理论支撑。通过加强三者的动态联动和区域适配性调控，能有效促进耕地利用的现代化与可持续发展。

耕地利用隐性形态的耦合协调关系在不同发展阶段呈现不同特征(图 1)。耦合是指投入、组织适应性和产出三维形态之间的相互作用和联动性，是耕地利用隐性形态水平提升的核心动力；协调则是耦合关系逐步发展的结果，体现为三维之间的匹配和平衡。

1) 初级阶段：低耦合、弱联动。

耕地利用以粗放型模式为主，投入结构单一且效率较低，投入对产出的支撑能力弱，产出效益有限，组织适应性缺乏有效调控能力，资源配置缺乏优化机制，整体耕地利用效率低。此阶段三维形态联动性较弱，耦合协调状态表现为严重失调或濒临失调。

2) 中级阶段：耦合增强、联动机制逐步形成。

随着社会经济发展、技术进步和政策引导，投入结构逐渐优化，产出效益显著提升，组织适应性初步发挥调节作用，资源配置效率开始改善。此阶段，三维形态之间的耦合关系明显增强，联动机制逐渐形成，耕地利用开始向集约化、高效化方向发展，耦合协调状态表现为初级协调或中级协调。

3) 高级阶段：高度耦合、良性循环。

投入要素的精准配置、产出的多样化结构与组织适应性的动态调节能力共同作用，三维形态相互作用机制不断完善，耦合关系达到高度联动状态，形成良性循环，耕地利用效率全面提升。此阶段耕地利用隐性形态展现出高度灵活性和适应性，三维耦合协调状态表现为高级协调。

不同区域的耦合演进存在显著差异。平原地区依托优越的自然条件和完善的基础设施，耦合演进速度更快。而丘陵山地区由于地形复杂、耕地破碎等因素，资源配置效率较低，三维形态的耦合协同优化受到限制。需基于本框架的系统性、动态性和空间适配性等特征，通过优化资源配置模式和提升投入效率，克服自然条件的不利影响，推动耦合关系的提升。未来，随着农业智能化、绿色化和精准化技术的推广，耕地利用隐性形态三维耦合关系将进一步优化，推动耕地资源的高效利用和可持续发展。

重庆市位于中国西南部，东与湖北省接壤，南邻湖南省和贵州省，西邻四川省，北邻陕西省，处于东部经济发达区与西部资源丰富区的交会地带(图 2)，辖区面积8.24×10⁴ km²，下辖38个区县，实施“一区两群”发展战略，即主城都市区(包括21个区县)、渝东北三峡库区城镇群(包括11个区县)和渝东南武陵山区城镇群(包括6个区县)。重庆市地形以山地和丘陵为主，山地面积占76%，丘陵占22%，河谷平坝区仅占2%。受复杂地形影响，耕地资源呈现“总量有限、分布零散、质量梯度差异显著”的特征。重庆市第三次国土空间调查显示，全市耕地总面积为1.87×10⁶ hm²，其中，15°以上的坡耕地占耕地总量的35%，25°以上的坡耕地以及石漠化地区、林区、河道周边等不稳定耕地占耕地总量的17%。总体来看，重庆市耕地零散破碎，多为“巴掌田”和“鸡窝地”，耕地利用率低，产业规模化发展受限。2022年，全市粮食总产量为10.73×10⁶ t，同比下降1.8%。城镇化进程对耕地利用格局产生深远影响。截至2022年，重庆市常住人口约3213.34万人，常住人口城镇化率为70.96%，高于全国平均水平，但全市不同区域的城镇化水平差距大。以重庆市为研究区，深入研究丘陵山地区耕地利用隐性形态的内部耦合变化规律及其驱动力，在对山区耕地可持续发展的研究中具有典型代表意义。

研究区涵盖重庆市除渝中区外的37个区县(渝中区无耕地，不纳入研究范围)，研究单元为区(县)级尺度，以2010、2015、2020、2022年为研究时间节点，涉及数据主要来自《重庆市统计年鉴2011》^[42]、《重庆市统计年鉴2016》^[43]、《重庆市统计年鉴2021》^[44]、《重庆市统计年鉴2023》^[45]以及相关区域的国民经济与社会发展统计公报和土地利用遥感数据(表 1)。

基于耕地利用隐性形态概念内涵及相关研究^[39]，结合区域特性与数据可得性，从投入、产出与组织适应性3个维度选取指标，保证指标的综合性与包容性，合理规避共线性，构建耕地利用隐性形态复合评价指标体系(表 2)，综合表征耕地利用隐性形态。

对表 1中的部分指标说明如下：耕地的社会保障价值主要体现在其农产品供给能力上，尤其是粮食生产能力，它直接关系到区域粮食安全与社会稳定，故单位面积耕地的粮食产量可作为衡量其社会保障价值的重要产出指标；耕地的生态价值产出具有多维特征，其中植被覆盖及其生长状况是生态功能的重要组成部分^[39]，NDVI(归一化植被指数)可有效表征植被的生长状态和覆盖度，进而反映耕地生态系统的健康状况，衡量耕地利用的可持续性；选取从业者生计价值表征农户实际收入，衡量耕地对当地农户的生计保障程度；耕地整合度采用平均地块面积衡量，反映耕地的分散或集中程度，表征耕地的规模化经营潜力；耕地肥力是耕地基础地力形成的基础^[46]，耕地表层有机碳含量是耕地肥力的核心指标^[47]，可综合反映土壤保水保肥能力，适用于区域尺度研究；碳排放强度指标用来衡量耕地利用的环境影响及其可持续生产水平，反映农业生产过程中碳排放的强度及对生态系统的潜在影响。

运用熵值法确定指标权重并采用极值处理法消除指标量纲的影响。正向指标标准化公式为：

负向指标标准化公式为：

其中：x_j^′为第j项指标的标准化值；x_j为第j项指标的原始值；x_max为第j项指标的最大值；x_min为第j项指标的最小值。

引入“耕地利用隐性形态指数”，以实现耕地利用隐性形态的量化比较与分析。指数取值区间为[0, 1]，数值越大表示该指标的作用越突出。区域耕地利用综合隐性形态指数值为3个维度耕地利用隐性形态指数值均值，各维度耕地利用隐性形态指数值的计算公式如下：

式中：N_i，d为第i个研究单元耕地利用隐性形态维度d的指数值；n为特定隐性形态准则层下的指标个数；L_ij为第i个研究单元的第j个指标标准化后的值；W_j为第j个指标的权重。

耦合是两个及两个以上的系统之间相互作用、相互依赖的动态过程^[23]。本文借助耦合协调度模型，探究耕地利用隐性形态投入—组织适应性—产出3个维度的耦合作用关系。具体公式如下：

其中：C为耦合度，取值范围为[0, 1]，其值越大说明各维度之间相互影响作用越强烈，P_i、R_i、E_i分别为投入、组织适应性和产出隐性形态的指数值，i∈{1，2，…，37}，n表示维度个数；T为耕地利用隐性形态的综合协调指数，α、β、γ分别为各维度的权重系数，本文取α=β=γ= 1 3；D为耦合协调度，取值范围为[0, 1]，表征各功能之间的互动关系，数值越大表明耦合协调水平越好。为了更好地体现耦合协调水平，参照已有研究^[48]，采用均值分段法，将耕地利用隐性形态投入—组织适应性—产出的耦合协调度D值划分为5种类型(表 3)。

核密度估计方法的优点在于不依赖数据分布的假设，也无须建立参数模型，能够通过连续的密度曲线直观展示变量的阶段性和动态变化特征。本文采用核密度估计方法对耕地利用隐性形态投入—组织适应性—产出耦合协调度进行分析，以刻画其时序特征与演进规律。具体公式如下：

式中：n为样本个数；H为带宽；x_i为样本观测值；x为样本均值；K(·)为核函数，具体采用高斯函数进行估计。核密度曲线的波峰数量反映了极化特征，波峰的高度体现了差异特征，曲线的拖尾厚度表示高(低)值区域的数量，而拖尾长度则揭示了高低值区域的演化特征。

在解析具有显著空间差异的地理归因和空间分异特征方面，地理探测器的稳健性和精确性已被学术界广泛验证^[48]并运用于自然和经济社会等领域中各种现象的驱动力和影响因素研究^[49]，由因子探测、交互探测、风险探测和生态探测4个基础模块组成。本文利用因子探测模块研究耕地利用隐性形态投入—组织适应性—产出耦合协调度演化过程中，多维因子X对耦合协调度Y的解释程度，为制定科学合理的政策提供依据，具体公式如下：

式中：q∈[0, 1]，为该因子对耦合协调度的解释力，值越大表明该因子对耕地利用隐性形态投入—组织适应性—产出耦合协调度的影响程度越强，反之越弱；L为X或Y的分层；h为耕地利用隐性形态投入—组织适应性—产出耦合协调度的等级；N和N_h分别为整个区域和第h级的样本个数；σ²和σ_h²分别为整个区域和第h级的耦合协调度。

基于前期数据计算结果，利用ArcGIS软件分析重庆市2010、2015、2020和2022年耕地利用投入、组织适应性和产出形态的时空分异结果(图 3)。

2010-2022年间，重庆市耕地利用投入、组织适应性和产出形态发生了显著变化，呈现出不同的空间分化态势。2010年，重庆市大部分研究单元的耕地利用投入和产出形态值集中在0.2~0.4，无明显分化。组织适应性形态值集中在0.4~0.6，呈西高东低格局，具体为主城新区水平高于其他区域。2015年，产出形态分化显著，呈西高东低空间分化，中心城区与主城新区的产出优化是产生分异的主要原因。2020年，重庆市各区域的投入形态持续增加，出现东高西低的显著空间分化，产出形态值则持续上升，仍为西高东低格局。至2022年，全市已建成高标准农田约8.78×10⁵ hm²，显著提高了耕地的粮食产能，单位面积粮食产量增加了10%至20%，并且具有节水、节电、节肥、节药的效果，对传统投入(化肥、农药等)依赖下降，促进了投入结构的优化与调整，大幅提升了组织适应性形态水平，呈现出较强的空间优化态势。

分析3个维度的构成指标可知，重庆市耕地利用隐性形态整体水平显著提升，但各维度分异与结构性矛盾并存，具体表现为：

1) 投入形态：传统要素流失快，现代化要素补充不足。

2010-2022年，投入形态呈现“先优化后衰退”的阶段性变化，其核心矛盾为传统投入要素(劳动力、化肥)的快速流失与现代化要素(机械动力、能源)的补充滞后。机械动力投入在2010-2022年间增长率为2.8%，提升缓慢。受社会经济发展水平和地形限制，渝东北、渝东南地区机械化推进速度和水平显著低于主城都市区；农村用电量年均增长率缓慢上升，但空间上仍呈主城新区提升速度快于渝东北、渝东南地区，农业现代化仍存在区域差异。劳动力、养分与农药投入呈下降趋势。主城都市区劳动力投入年均下降约3.5%，渝东南地区劳动力投入下降较慢，对劳动力和农药化肥等依赖性仍较强。劳动力的大量非农就业转移及对施用农药化肥等的控制是投入形态后期表现欠佳的主要原因。

2) 组织适应性形态：土地整治驱动下的区域差异收敛。

组织适应性形态的提升主要得益于重庆市土地整治、高标准农田建设下耕地整合度及耕地质量的显著提高。随着城镇化进程的加速，乡村劳动力转移缓解了耕地资源压力，加上城镇化带动的交通改善、资金和信息流入、农产品需求扩张和农业生产的集约化经营，促进了现代化与规模化生产。土地资源整合和集约利用带来的管理效率提升，进一步优化了耕地利用的组织适应性。2022年，重庆市高标准农田建设面积达8.78×10⁵ hm²，推动耕地整合度提升(如巴南区的耕地平均地块面积从7.60 hm²/块增至8.72 hm²/块)，土地整治政策虽缓解了主城都市区的“细碎化-低效化”矛盾，但对渝东北、渝东南地区的辐射效应有限。

3) 产出形态：经济-生态功能协同性增强。

产出形态的优化滞后于投入结构调整，2010-2022年经济价值产出年均增长约5.2%，生态价值(NDVI)提升约12.6%。但后期受市场波动和资源约束影响，2020年后部分区域(如巫溪县)出现“高投入-低产出”的情况，反映山区耕地系统韧性不足。

总体来看，2010-2022年，重庆市耕地利用三维隐性形态变化呈现出显著的时空不均衡性。与投入形态和组织适应性形态相比，产出形态的提升相对滞后。在空间格局上，投入和产出形态的区域分化加剧，中心城区和主城新区的优势明显，而渝东北、渝东南地区受人口流失、耕地本底条件和基础设施等因素影响，发展相对滞后，且内部各维度发展水平参差不齐。

投入形态是三维隐性形态优化的基础，其变化趋势直接影响产出形态和组织适应性形态的变化。重庆市耕地利用投入形态呈“传统要素流失加速、现代化要素增长缓慢”的特征。主城都市区依托政策干预和技术进步，投入结构不断优化，现代化要素占比提升。而渝东北、渝东南地区受社会经济和自然条件制约，传统投入要素流失加剧，现代化要素补充不足，导致投入形态结构性失衡，制约了三维隐性形态的整体优化。因此，需采取差异化的优化策略，以改善耕地投入形态，提高耕地利用的整体效能。

根据重庆市耕地利用隐性形态投入—组织适应性—产出耦合协调度测算结果，借助核密度估计方法，绘制2010-2022年重庆市耕地利用隐性形态耦合协调度演变核密度曲线图，刻画2010-2022年耦合协调度的时序演变特征(图 4)。

2010-2022年，重庆市耕地利用隐性形态投入—组织适应性—产出三维耦合协调度呈现持续提升的趋势，核密度曲线在横轴的投影不断向右延展，反映整体耦合协调水平逐年提高。

从曲线重心位置来看，曲线重心持续向右迁移，大致由0.65迁至0.80附近，耕地利用隐性形态三维耦合协调度逐年提高。从波峰变化来看，2015年和2020年的主峰显著高于2010年，耦合协调度的分布逐渐趋于集中；然而，2022年的主峰高度有所下降，表明区域差异有所增大。从波峰数量来看，2010年曲线呈单峰分布，表明区域间协调度水平较为均衡，主峰位于0.6~0.7区间，集中于中等协调水平；2015年后，逐渐形成“一主一副”双峰分布，高协调度区域增加的同时，部分区域仍处于较低水平，两极分化趋势显现。从曲线拖尾特征看，2010-2020年曲线左侧拖尾明显，中等协调度区域占比较大；2022年右侧拖尾延长，高耦合协调度区域占比进一步提升。

整体来看，重庆市耕地利用隐性形态的三维耦合协调度在不断提高。但因区域社会经济发展与耕地利用水平的非均衡性，区域间耦合协调水平的非均衡性加剧。部分区域(如主城新区)，耕地利用效率和可持续发展能力增强，协调度持续提升，而渝东北、渝东南部分区县三维协调度提升缓慢，仍存在低协调度区域。

基于2010、2015、2020和2022年耕地利用隐性形态投入、组织适应性和产出形态的耦合协调度计算结果，利用ArcGIS软件分析重庆市耕地利用隐性形态三维耦合协调度的空间演化特征(图 5)。

2010-2022年，重庆市耕地利用隐性形态的耦合协调度整体呈逐步提升态势，但区域间协调水平呈现组团式极化分布特征，主城都市区和渝东北地区协调水平较高，部分区域仍存在濒临失调或严重失调现象。具体空间演变特征如下：

1) 2010年，中级协调区域占主体，局部区域呈低水平协调。

2010年，重庆市大多数区域耕地利用隐性形态的耦合协调度处于初级协调或中级协调水平，中级协调区域占比达56.75%，为研究区的主体。部分区域协调水平较低，如江北区和渝北区，因投入和产出水平偏低(＜0.3)，尽管组织适应性较高(＞0.5)，但系统内部失衡导致整体协调水平处于濒临失调状态。巫溪县则表现为“高投入-低产出”的结构性失衡，区域投入要素结构不合理，利用效率较低，难以实现高效产出。

2) 2015年，中级协调区域扩大，高级协调区域初步形成。

2015年，重庆市耦合协调水平显著提升，中级协调区域占比增至78.38%，主城都市区协调水平提升明显。其中，荣昌区和永川区凭借三维均衡发展，升级为高级协调区域。渝北区因组织适应性形态优化投入要素结构，带动产出效率提升，耦合协调水平从濒临失调改善至初级协调。然而，江北区和巫溪县依然处于濒临失调状态，但表现形式和成因各异：江北区受城市化挤压影响，耕地利用投入和产出水平持续下降，组织适应性的提升难以弥补资源投入不足的影响；巫溪县则因自然条件和市场限制，高投入难以转化为高效产出，组织适应性作用有限，其失调体现了山区农业区“高投入-低产出”的结构性矛盾。

3) 2020年，中心城区协调度下降，区域差异显著。

2020年，重庆市的中级协调区域比例维持在78.38%，但内部空间结构变化显著，呈现“中心塌陷，边缘优化”格局。中心城区部分区域协调度由中级降至初级，农业生产功能进一步边缘化，仅江北区从濒临失调状态升级为初级协调的优化状态。巫溪县通过土地整治优化耕地本底条件，组织适应性显著提升，协调水平从濒临失调升级至初级协调。整体上，区域间协调度差异加剧，快速城市化对中心城区的耕地利用带来了显著冲击：耕地资源减少、投入维度衰退，同时农业产出功能弱化，导致中心城区农业生产功能逐步边缘化，仅部分区域因政策干预实现小幅提升。

4) 2022年，高级协调区域增加，局部协调性下降明显。

2022年，高级协调区域占比上升至45.86%，集中于主城都市区和渝东北地区，得益于长期高水平的组织适应性对投入与产出的协同优化作用。中级协调区域占比45.95%，仍为主体。然而，部分区域协调度退化明显。江北区受限于有限的耕地资源和发展空间，投入与产出进一步减少，耦合协调度退化至严重失调状态，耕地利用系统内部失衡加剧。

总体来看，2010-2022年，重庆市耕地利用隐性形态的耦合协调度优化过程呈现两大特征：一是中级协调区域长期占主导，高级协调区域在后期显著增加，特别是在主城都市区和渝东北地区，长期高水平的组织适应性优化了投入与产出的协同关系；二是协调度的区域差异日益显著，主城都市区和渝东北地区协调水平较高，而部分区域受自然、经济和政策等多重因素影响，协调度退化或失衡。

重庆市耕地利用隐性形态投入—组织适应性—产出耦合协调度水平受多因素综合影响。根据地理探测器分析结果，各影响因子对耦合协调度的解释力(q值)存在显著差异(图 6)。

依据q值在不同研究时间点的表现，可将影响因素分为4类(表 4)：主要因素(在3个及以上研究节点q≥0.5，长期占据主导地位)、次要因素(在两个研究节点q≥0.5，且无强烈波动，具有阶段性影响)、一般因素(仅在一个研究节点q≥0.5，影响力有限)和其他因素(在所有研究节点q＜0.5，影响力较弱)。结合丘陵山区耕地利用特性，丘陵山区受地形限制，机械化推进速度较缓慢，因此，机械动力投入(X₁)不作为主要因素，而归为次要因素。生态价值(X₈)2010、2015、2020和2022年的q值分别为0.990 1、0.003 8、0.039 6、0.582 4，存在异常波动，因此不作为次要因素，归为其他因素。

主要因素长期稳定影响耦合协调水平，包括能源投入(X₅)、从业者生计价值(X₉)、耕地地势平坦度(X₁₃)和耕地肥力(X₁₅)4项。

1) 能源投入(X₅)。反映农业现代化的投入水平及技术支持作用。2010-2020年q值均大于0.5。然而，2022年q值下降至0.32，表明其影响力逐步被其他现代化因素取代，体现了农业投入结构的优化升级。

2) 从业者生计价值(X₉)。该因素是农民耕地经济效益的核心指标，是驱动农户优化耕地利用的重要动力。其q值在2010-2020年保持在0.6以上，2022年降至0.48，仍具较大影响。高生计价值产出的农户可通过组织适应性更强的方式(如合作化经营、机械化生产)优化资源配置，提高土地利用效率。

3) 耕地地势平坦度(X₁₃)。2010-2022年，耕地地势平坦型q值均＞0.68，明显高于平原地区^[27]，表明地形是丘陵山地区耕地系统演化中的重要约束因素。重庆地处丘陵山地区，地形复杂且起伏大，限制了机械化的推广空间。因此，耕地坡度对耕地投入、组织适应性和产出的协调发展具有重要影响。

4) 耕地肥力(X₁₅)。耕地肥力是决定耕地的生产能力的核心自然属性。良好的耕地肥力能够提供更高的有机碳含量，有利于提高耕地产出效益并稳定农作物生产能力，其对耦合协调度的影响稳定而深远。

次要因素在特定阶段对耕地利用影响较大，包括机械动力投入(X₁)、耕地相对丰富性(X₁₀)、耕地区位邻近性(X₁₁)和耕地整合度(X₁₄)4项。

1) 机械动力投入(X₁)。农业机械动力替代传统劳动力，显著提高农业生产效率，促进了耕地要素的集约化配置。2015年后，随着农业技术进步和集约化生产的推进，机械动力对劳动力流失的替代效应显现，大幅提升了农业产出效率，优化了耕地要素的配置结构。主城新区的永川区等以农业机械化推进为基础，建设机械化生产基地、改善宜机化生产条件，大幅提升了耕地利用系统的组织适应性和资源利用效率。

2) 耕地相对丰富性(X₁₀)。人均耕地面积是衡量劳动力与耕地资源匹配度的重要指标，较高的耕地相对丰富性能为农业从业者提供更好的资源支撑，推动现代农业组织模式的推广与优化。

3) 耕地区位邻近性(X₁₁)。耕地与市场的距离直接影响其经济价值和产出方式。主城新区的耕地由于接近市场，利用效率较高，远离城市的渝东北、渝东南地区受交通不便影响较大。随着交通条件的改善，这一影响逐渐减弱。

4) 耕地整合度(X₁₄)。耕地整合度反映了耕地空间布局的细碎化程度和资源整合效率，影响规模化经营潜力。整合度高的耕地更适合规模化经营和机械化作业，有助于提高耕地资源利用效率。主城新区的耕地整合度高，而渝东北和渝东南地区耕地较为细碎，限制了规模化、机械化的推广。与耕地地势平坦度和耕地肥力相比，耕地整合度对隐性形态的影响更多体现为对耕地资源整合或组织效率的间接影响。

一般因素包括劳动力投入(X₂)、养分投入(X₃)、农户去依赖性(X₁₂)和耕地碳排放强度(X₁₆)4项，这些因素对耦合协调度的影响主要集中在2010年，随后逐渐减弱。

1) 劳动力投入(X₂)。该因素直接影响农业生产模式的选择，但随着机械化的推广，其影响逐渐减弱。2010年耕地利用对劳动力的依赖较大，而后期其影响逐渐被技术进步所取代。

2) 化肥投入(X₃)。化肥施用可以短期提升产量，但长期依赖会导致土壤质量下降。随着有机肥和生态农业的推广，其影响逐渐减弱。

3) 农户耕地去依赖性(X₁₂)。城镇化率影响农村劳动力结构与耕地利用方式。虽然其对耕地利用的长期影响是显著的，但在短期内，其作用弱于能源投入、耕地地势平坦度等因素。

4) 耕地碳排放强度(X₁₆)。该因素体现了耕地利用的生态效益。随着农业现代化和绿色发展模式的推广，其影响作用逐渐减弱。

其他因素包括农药投入(X₄)、社会保障价值(X₆)、经济价值(X₇)和生态价值(X₈)4项，对耦合协调水平的影响较弱，多为间接影响或阶段性作用。

综上，能源投入、从业者生计价值、耕地地势平坦度和耕地肥力是影响重庆市耕地利用隐性形态耦合协调度的主要因素，其作用长期稳定。机械动力投入和耕地整合度等次要因素在特定时间节点产生较强作用，但其影响具有阶段性。劳动力投入、养分投入、农户耕地去依赖性和耕地碳排放强度等一般因素随着农业生产方式变革逐渐减弱。其他因素则主要表现为间接影响，反映耕地利用系统的复杂性。整体来看，重庆市耕地利用隐性形态耦合协调水平是经济、社会和自然因素相互交织的结果，这些因素共同推动了农业现代化发展与耕地可持续利用。

本文聚焦于快速工业化与城镇化进程中区域耕地利用隐性形态的时空演化特征，基于耕地利用隐性形态的投入形态—组织适应性形态—产出形态，重构理论内涵、创新解析框架。通过综合指数模型，量化耕地利用三维隐性形态，并探析三维隐性形态时空演化规律，识别影响隐性形态耦合协调演化的驱动因素，结论如下：

1) 2010-2022年，重庆市耕地利用投入—组织适应性—产出三维隐性形态整体水平明显提高，但各维度空间依赖性明显，呈现相似组团分布特征，空间分异显著。整体来看，研究区投入形态先优化后衰退，后期的衰退主要受丘陵山地区农业机械化推进缓慢影响，未能有效弥补农业劳动力非农转移导致的要素流失。同时，化肥施用控制导致传统投入要素数量减少，而新型要素补充不同步，加剧了要素替代失衡，形成投入结构性空缺。从具体时间节点看，2010年，三维形态整体水平较低，空间分化不明显，组织适应性呈现“西高东低”格局；2015年，产出形态显著优化，中心城区与主城新区的产出水平显著高于渝东北和渝东南地区；2020年，投入形态的空间分化加剧，呈“东高西低”空间分布格局；2022年，耕地利用系统进一步优化，投入形态结构显著调整，产出形态虽有局部波动但整体提升，组织适应性大幅改善，表现出较强的区域适应能力。整体来看，中心城区和主城新区耕地资源禀赋较优，三维隐性形态表现较为突出，而渝东北和渝东南地区由于资源禀赋质量劣势性、经济发展水平弱势以及基础设施等限制，发展进程相对滞后。

2) 2010-2022年，重庆市耕地利用隐性形态的三维耦合协调度持续提升，整体水平从初级协调向高级协调迈进，但区域差异和动态变化仍较为显著。主城都市区资源流失和农业功能的弱化，使整体协调度呈现“中心塌陷，边缘改良”格局。耦合协调度的时空演化特征与三维隐性形态演化存在一定相似性：2015年，中心城区与主城新区率先达到较高协调水平；2020年，受快速城镇化影响，中心城区协调度出现下降，但渝东北地区因土地整治工程协调度有所改善；2022年，高级协调区域显著增加，主城都市区和渝东北地区表现突出。然而，部分区域(如江北区和巫溪县)仍存在耦合协调度下降或失调的现象，耕地利用系统内部结构和功能的动态平衡较为复杂，协调提升过程伴随波动和区域分化。

3) 重庆市耕地利用隐性形态的三维耦合协调由经济、社会和自然因素共同驱动。经济因素中，机械动力的提升增强了产出效率和组织适应性；农村用电量代表的能源投入反映了农业现代化水平，对三维形态的优化具有重要作用。社会因素中，耕地相对丰富性和城镇化率对耕地资源配置、劳动力转移和产出效率的协同作用显著，但其影响因各区域社会经济发展水平不同而存在显著差异。自然因素中，耕地地势平坦度和耕地肥力决定了耕地利用的基础条件。此外，农业技术进步、土地整治、高标准农田建设等措施在一定程度上促进了三维耦合协调度的提升，但区域发展仍存在不均衡性，需通过资源优化配置和政策干预进一步改善。

本文基于投入—组织适应性—产出三维耦合框架，揭示了快速工业化与城镇化进程中丘陵山地区耕地利用的三维隐性形态演化特征，为阶段化政策干预提供依据，并为差异化耕地管理政策设计提供量化依据。

1) 本文构建投入—组织适应性—产出三维框架，将投入和产出视作不同子系统，并叠加组织适应性，全面揭示了丘陵山地区耕地系统的复杂互动关系。相比传统投入—产出二维分析，该框架扩展了研究范畴，使其更好与区域城乡综合发展情况衔接，深化了土地利用转型理论，为耕地利用优化提供了新视角。丘陵山地区耕地管理政策应聚焦于优化投入—组织适应性—产出三维隐性形态的耦合协调，以促进农业可持续发展。其中，农业现代化与自然条件的协同是提升耦合协调度的核心驱动力。当前丘陵山地区耕地利用隐性形态面临的关键问题是要素替代失衡，即劳动力转移加剧但技术替代进程缓慢。因此，重庆市应采取差异化“分阶段干预”政策，构建“集约投入—高组织适应性—多样化产出”的优化路径。对于形态发展不均衡的区域，应重点提升劣势维度的核心指标，以引导区域三维形态的协调演化。

2) 重庆市的耕地利用隐性形态时空演化格局趋势与平原地区整体相似^[39]，均呈现出隐性形态水平与耦合协调水平逐渐优化的态势。然而，相较于平原地区，丘陵山地区的时空演变相对缓慢^[50]，且空间极化现象更为突出。受地形复杂、生态脆弱和交通不便等因素影响，技术替代进程受到制约^[28]，导致隐性形态的提升进程相对缓慢，区域差异显著。在劳动力流失或资金投入不足时，部分区域可能出现停滞或下降。丘陵山地区的地块破碎化和资源分散使得平原地区的主要规模化经营模式难以直接适用。丘陵山地区应依托土地整治、耕地宜机化改造和适地化技术应用，提高生产效率，弥补传统规模化不足。随着农业机械化水平提升，农业现代化的替代效应凸显。因此，重庆市应遵循“技术补偿地形”的山区农业路径，采取“适地化技术+精准化政策”双轮驱动发展策略。社会经济发达、基础设施完善的主城新区，具备产业带动能力强、资金投入充足、科技创新能力高等优势，可探索“耕地多功能复合利用”模式，在保障粮食安全的同时提升耕地的生态与经济价值，以缓解城市化对耕地功能的挤压；同时，通过政策引导企业参与耕地流转，推动适度规模化经营。渝东北和渝东南地区应优先推进土地整治与农田水利建设，提升耕地整合度，同时推广适地化农业机械，发展高效绿色农业，以缓解地形限制带来的高投入—低产出的结构性矛盾。

3) 本文探讨了重庆市耕地利用三维隐性形态的耦合协调水平演化及其驱动因素，揭示了耕地利用与地理环境、社会经济及城镇化进程之间的复杂适应关系，丰富了丘陵山地区耕地利用隐性形态内部系统的互动机制研究，但研究仍存在局限性：一方面，研究结论主要适用于丘陵山地区，对其他复杂地形区的适用性需进一步验证；另一方面，受数据可得性限制，未能进行更小尺度、更长时间序列的隐性形态测度。未来可结合新型数据源与方法，增强隐性形态表征的全面性与精准性，扩展研究的时空范围，覆盖不同地貌类型的典型区域，深化多尺度研究，以期更全面、深入地揭示耕地利用隐性形态的演化规律。