贫困地理学的研究对象是“贫困地域系统”，也是人地关系在贫困地域的具体表现形式^[23]. 根据地域类型分类，贫困地域系统可划分为“乡村贫困地域系统”“城镇贫困地域系统”“城乡交错贫困地域系统”3种类型. 国内学者率先将贫困地域系统定义为：在特定贫困地域内，由地理区位、生态环境、资源禀赋、经济水平、社会福利等要素相互作用构成的，具有一定功能和结构的复杂系统. 将贫困地域系统分解为3个子系统，即根据人地关系地域系统中的“人”“地”“业”要素因子，划分为以“人”为核心要素的社会子系统，以“地”为核心要素的环境子系统，以“业”为核心要素的经济子系统^[9]. 在系统科学理论、人地关系地域系统理论、空间贫困和多维贫困理论、乡村地域系统理论等基础理论的指导下，本文将乡村贫困地域系统概念初步界定为：在特定时空范围的乡村地域，由“人”“地”“业”要素构成的系统内部存在不同程度的维度剥夺或耦合关系失调，导致区域内部各类资本要素供给不足、结构低效或动力失调，在自然、社会和经济方面处于绝对或相对贫困状态的乡村地域系统.

根据上述对乡村贫困地域系统的内涵挖掘和基本概念界定，乡村贫困地域系统多维结构可被表征为“人”“地”“业”3类核心要素所衍生的关联子系统及多维结构. “人”“地”“业”分别衍生对应人口子系统、土地子系统和产业子系统. 3个子系统互相关联、互相影响、互相制约，共同影响区域乡村贫困地域系统演化方向、进程和状态.

1) 以“人”为核心要素的人口子系统，表征了贫困地区因自然资源环境不足以支撑人口高质量发展，或因社会资源、公共政策等供给不足，导致个体社会福利缺失，反映了贫困与社会资源及福利供给之间的关系，并直接反映区域社会贫困演化状态.

2) 以“地”为核心要素的土地子系统，表征了贫困地区在地理区位、资源禀赋、生态环境等方面的劣势，反映了贫困与自然地理环境之间的关系，并直接反映区域自然贫困程度，同时也是区域社会贫困和经济贫困的根本原因.

3) 以“业”为核心要素的产业子系统，表征了乡村贫困地域在产业基础、产业结构、产业发展、要素流通、市场开发潜力等方面的不足，主要反映贫困与区域经济发展之间的关系，并直接反映区域经济贫困程度.

“人—地—业”各子系统按照“基底—支撑—潜力”的递进式逻辑，进一步解构为：人力资源本底、人力资源保障、人力资源开发、土地资源本底、土地生态承载、土地功能开发、产业空间区位、产业收益保障、产业发展潜力共9个维度(图 1).

构建科学合理的评价指标体系是开展乡村贫困地域系统演化研究的重要基础. 本文以目标性与代表性、科学性与简明性、时效性与地域性、可量化与可获取为评价指标遴选原则，参照近年来有关人地关系地域系统、贫困地理学、发展地理学、乡村地理学、人口地理学等相关研究中的评价指标设置，结合研究区实地调查可获取统计数据或已有的相关数据产品，按照人口、土地、产业3个子系统的多维结构和“基底—支撑—潜力”递进式逻辑，每个维度选取不少于20个候选指标，共计选取66项候选指标构成了评价指标候选集合，采用极差法对评价指标进行标准化处理.

1) 人口子系统评价指标候选集. 人力资源本底维度主要遴选了能够体现村域人口数量规模和机构特征的基底型指标；人力资源保障维度主要遴选了能够体现村域人口发展状态和基本保障水平的支撑型指标；人力资源开发维度主要遴选了能够体现村域人力资源开发前景和转型动力的潜力型指标.

2) 土地子系统评价指标候选集. 土地资源本底维度主要遴选了能够体现村域土地资源特征和土地利用状态的基底型指标；土地生态承载维度主要遴选了能够体现村域土地对生态环境承载作用及能够产出生态资产价值的支撑型指标；土地功能开发维度主要遴选了能够体现村域土地多功能价值和转型开发利用前景的潜力型指标.

3) 产业子系统评价指标候选集. 产业空间区位维度主要遴选了能够体现村域经济贫困地理特征的基底型指标；产业收益保障维度主要遴选了能够体现村域产业发展状态和产业发展优势的支撑型指标；产业发展潜力维度主要遴选了能够体现村域产业发展韧性和转型潜力的指标.

本文综合运用“变异系数分析(CV)+方差膨胀系数(VIF)”复合方法，去除部分共线性指向指标和贡献度相对不足的指标，最后遴选保留了31项评价指标，通过CRITIC(Criteria Importance Though Intercrieria Correlation)赋权方法测算得到各项评价指标权重，进而构建了乡村贫困地域系统演化发展评价指标体系(表 1).

为科学表征乡村贫困地域系统“人—地—业”演化状态，本文基于发展地理学视角，参考运用生计多边形法测度乡村贫困地域系统多维发展指数^[24-25]，测度过程考量了不同维度指标之间的相互影响和相互制约作用. 按照人力资源本底、人力资源保障、人力资源开发、土地资源本底、土地生态承载、土地功能开发、产业空间区位、产业收益保障、产业发展潜力共9个方面及相应评价指标，构建“加权生计多边形”乡村贫困地域多维发展指数测度模型，并结合本文评价指标维度分解情况，运用生计9边形模型测度乡村贫困地域系统演化发展指数(图 2).

O为原点，其中原点O与各顶点间的直线为评价维度，原点O与各顶点方向直线实测线段长度(O-HB，O-HS，O-HD……O-IE)为该子维度得分，O-HB，O-HS，O-HD……O-IE分别以A，B，C……I表示，任意两个子维度之间夹角为α(α=360°/9)，则9边形“A-I”的面积即可被定义为某区域乡村贫困地域系统演化发展指数. 在9个维度中，因排序组合方式不同，导致其面积计算结果也不同. 因此，本文采用全排序法取9个维度得分两两相乘后的加总值，作为某区域乡村贫困地域演化发展指数. 全排序多边形测度模型符合乡村贫困地域的系统化和复杂性特征. 计算公式为：

式(1)、式(2)中，MRPI_p为某评价单元乡村贫困地域系统演化发展指数；A为人力资源本底维度演化发展指数得分；w为某评价指标权重；y为评价指标标准化值；B-I维度演化发展指数得分测度方法同上.

冷热点空间分析是在给定一组加权变量要素的基础上，借助ArcGIS PRO 3.0中Getis-Ord Gi分析工具，识别研究区内具有空间统计显著性的热点和冷点. z得分(临界值)和p值(显著性水平)均用于统计显著性的度量，逐要素判断是否拒绝零假设. 如果某要素z得分高且p值小，则表征高值空间聚类；如果某要素z得分低并为负数且p值小，则表征低值的空间聚类. z得分趋向越高或越低，则空间聚类密集度越大. 若z得分接近于0，则表示不存在明显的空间聚类. 为了更加直观地反映研究区某一时间节点变量要素在空间上的冷、热点分布情况，本文选取2013年、2016年、2019年和2021年4个时间节点，以乡村贫困地域系统演化综合指数测度结果为变量要素，分期识别具有统计显著性多维减贫指数的高值(热点)和低值(冷点)空间聚类区，可以更加直观地研判乡村贫困地域系统不同时点的演化特征.

时空立方体是依托ArcGIS PRO 3.0平台，聚焦数据的时间序列分析、空间与时间集成模式分析，是对时空数据进行可视化探测研究. 本文以研究区各村域2013-2021年乡村贫困地域系统演化综合指数为分析变量，通过制作变量数据聚合点创建乡村贫困地域系统演化时空立方体，使时间截点要素在空间上进行聚合，刻画分析时空演化状态和模式. ①分期变量与数据：将2013-2021年各村域贫困地域系统演化综合指数得分按照半年期进行平均切分，共得到17期数据；②生成时空立方体：借助ARCGISpro 3.0平台，以村域X，Y坐标位置创建时空聚合点，通过时空聚合点分别创建乡村贫困地域演化综合指数时空立方体；③设置参数：将创建的各类立方体输入新兴时空热点分析工具，以各时空立方体Count字段(事件个数之和)作为分析变量，鉴于村域地理分布的不规则性，以“K—最邻近”作为空间关系概念化选项，空间领域数设置为8，领域时间步长设置为6，定义全局窗口为整个立方体. ④结果输出：运行生成相应类别新兴时空热点栅格图层，主要内容包括新增热点、连续热点、加强热点、持续热点、逐渐减少的热点、新增冷点、连续冷点、加强冷点、持续冷点、逐渐减少的冷点等16种时空演化模式.

镇沅县位于云南省西南部，隶属普洱市，地处哀牢山与无量山之间，东经100°21′-101°31′、北纬23°24′-24°22′. 全县总面积4 223 km²，山区面积占97.7%. 全县辖8镇1乡109个村，总人口21.3万人. 2001年被认定为国家级扶贫开发重点县，2013年以来镇沅县累计减少贫困人口10 582户，37 548人，3个贫困乡镇、76个贫困村全部脱贫出列，2019年镇沅县退出国家级贫困县序列. 镇沅县地理环境、自然资源、人口民族、贫困特征等要素均具有典型的区域代表性，是开展西南山区乡村贫困地域系统时空演化格局研究的理想实证区.

实证研究区土地利用矢量数据主要来源于《云南省普洱市镇沅县2014年度土地利用变更数据》《云南省普洱市镇沅县2017年度土地利用变更数据》(镇沅县国土资源局)；归一化植被指数(NDVI)、植物净初级生产力(NPP)数据来源于美国航空航天局(NASA)数据发布中心的MODIS 13A2数据，原格式为HDF格式，将数据经过拼接、转投影为Albers双纬线圆锥投影，转格式为tif格式，最后经裁剪得到研究区域NDVI，NPP数据；DEM数据从USGS网站(http://earthexplorer.usgs.gov)下载，数据的空间分辨率为30 m. 海拔、坡度和坡向数据通过ArcGIS PRO 3.0中的spatial_analysis模块生成. 村域调查数据主要涉及人口结构、农业生产、基础设施、产业结构、经济收入等相关数据，数据时序为2013-2021年，其中2013-2017年数据为2018年镇沅县实地调查收集获得，2018-2021年数据是在原有数据结构基础上进行的补充调查结果.

根据全排列生计9边形模型计算得到镇沅县乡村贫困地域系统演化发展指数测度结果显示，2013-2021年镇沅县乡村贫困地域系统演化状态整体呈现出较快的发展变化态势，其中2013-2019年呈现出快速上升态势，2020-2021年呈现出增长滞缓并趋于稳定态势. 通过对2013-2021年109个评价单元乡村贫困地域系统演化综合指数的中位数、平均数对比分析可见，2013-2015年测度结果中的中位数高于平均数，说明各村域减贫效果差异性较为显著，减贫政策措施的普惠作用尚未充分发挥；从2016年起，测度结果中的平均数开始高于中位数，且平均数与中位数之差呈现出逐年增大态势，反映出研究区各村域减贫水平差异性逐步缩小，减贫政策措施的普惠作用逐步显现(图 3).

采用自然断点法将2013-2021年镇沅县各评价单元乡村贫困地域系统演化发展指数均值划分为5个等级，分别为低水平区(0~1.32]、较低水平区(1.32~1.87]、中等水平区(1.87~2.45]、较高水平区(2.45~3.02]和高水平区(3.02~4.22]. 从各等级类型区占比来看，低水平区占比为9.17%、较低水平区占比为25.69%、中等水平区占比为26.61%、较高水平区占比为26.61%、高水平区占比为11.93%. 从分析结果看，镇沅县乡村贫困地域系统演化水平呈现出明显的空间异质性，且总体趋于中等水平. 指数高值区主要分布在县域中部的恩乐镇恩乐村、团结村、民江村、玻烈村和中南部的按板镇部分村域；指数低值区主要分布在县域东北部的九甲镇、东部的和平镇和西部的振太镇等地(图 4).

乡村贫困地域系统演化综合指数可以有效表征区域综合减贫水平. 演化综合指数分布在置信区间“-3到3”中的村域，置信度为99%的统计显著性；置信区间“-2到2”中的村域，置信度为95%的统计显著性；置信区间“-1到1”中的村域，置信度为90%的统计显著性；置信区间0的村域，不具有统计显著性. 热点区域主要表征某一时间点相关村域减贫综合水平与其他村域相比相对较高，冷点区域表征某一时间点相关村域综合减贫水平与其他村域相比相对滞缓，但不能表征该区域减贫水平出现持续性增长或衰退. 结合相关村域的实际情况，通过对镇沅县2013年、2016年、2019年和2021年乡村贫困地域系统演化静态时空冷热点分布图(图 5)进行分析可知：① 4个时间节点县域冷热点空间分布相对稳定且呈现出显著的空间关联性和集聚特征，热点空间在不同时点的显著程度较为稳定，冷点空间显著程度按照时间由远及近逐步减弱；②热点空间呈现出中心集聚特征，主要集中在人口密度相对较大、耕地资源良好、距县域中心较近或产业相对发达的区域，如恩乐镇恩乐村人口密度大，民江村、团结村距离县政府距离近，按板镇红星村、文牍村人均耕地资源较为丰富；③冷点空间呈现出分散多点集聚特征，主要集中在人口密度相对较低、距离县域中心相对较远、耕地资源较为贫瘠的范围，且多处于哀牢山、无量山等坡度起伏度较大的区域，主要分布在和平镇、者东镇、九甲镇、振太镇的部分村域.

本文通过乡村贫困地域系统演化时空立方体生成了镇沅县2013-2021年乡村贫困地域系统时空演化模式. 通过分析可知，镇沅县各村域乡村贫困地域系统时空动态演化主要呈现出连续热点、加强热点、连续冷点和持续冷点4种模式.

1) 连续热点模式. 涉及评价单元数量69个，占比为63.30%. 连续热点刻画了相关变量在各个时点呈现出连续线性的较快演化状态，结合乡村贫困地域系统演化综合指数测度结果进行分析，显示该时空演化模式表征了镇沅县乡村贫困地域系统整体上经历了较快且连续稳定的减贫正向演化过程.

2) 加强热点模式. 涉及评价单元数量5个，占比为4.59%，且均分布于同一乡镇. 加强热点刻画了相关变量在时间序列上呈现出显著的突变性演化状态，结合乡村贫困地域系统演化综合指数测度结果进行分析，显示该时空演化模式表征了相关村域经历了显著快速的减贫正向演化过程.

3) 连续冷点模式. 涉及评价单元数量1个，占比为0.92%，主要分布在哀牢山脉附近. 连续冷点刻画了相关变量在各个时点呈现出连续缓慢的演化状态. 结合测度结果分析，显示该村域减贫综合指数变化幅度较小，且前期呈现出小幅度波动性变化，表征该区域乡村贫困地域系统经历了连续缓慢的减贫演化过程.

4) 持续冷点模式. 涉及评价单元数量10个，占比为9.17%，空间上呈现出分散集聚特征，主要分布在哀牢山、无量山附近. 连续冷点刻画了相关变量在各个时点呈现出强度相近的持续缓慢演化状态. 结合测度结果进行分析，显示该时空演化模式表征了相关村域乡村贫困地域系统经历了状态稳定的持续性缓慢减贫演化过程.

另外，未检测到模式涉及评价单元数量24个，占比为22.02%，表征相关评价单元乡村贫困地域系统时空演化状态与其他村域相比不具有显著性.

“空间贫困陷阱”是贫困地理学的重要开创性理论，Daimon^[13]认为“空间贫困陷阱”是由于地理位置偏远或迁移成本过高所导致的一种持续性贫困状态，这种持续性贫困状态可以理解为是在地理空间上与周边区域相比的一种空间相对贫困状态，其内涵更多聚焦于贫困人口分布的空间集聚性，而非“人—地—业”多维空间贫困瞄准. 本文引入乡村贫困地域“人—地—业”系统演化水平与模式等要素，进一步拓展“空间贫困陷阱”内涵特征，将人口、土地、产业要素在地理空间上呈现出的持续性多维相对贫困状态作为识别空间贫困陷阱的标准.

乡村贫困地域系统演化水平分级刻画了各评价单元乡村贫困地域系统演化的绝对状态，其评价得分可以从一定程度上反映区域多维贫困状态和减贫效果，但无法直观地体现区域减贫的时空动态效果. 基于时空立方体的时空冷热点模式分析刻画了各评价单元乡村贫困地域系统演化的相对变化状态，冷热点模式识别可以有效表征区域减贫动态特征，但无法直观地体现这种动态特征与减贫目标呈现正相关还是负相关. 基于上述分析，本文构建了涵盖乡村贫困地域系统演化水平和时空演化模式的“空间贫困陷阱”识别刻画模式，按照不同程度等级组合将空间贫困陷阱划分为3个等级(图 7).

根据上述空间贫困陷阱识别的逻辑对应关系，将各类乡村贫困地域演化水平等级图层与时空演化模式图层进行叠加统计分析，制作了镇沅县空间贫困陷阱分布图(图 8)，通过分析可知，镇沅县空间贫困陷阱主要分布在县域东西两端的哀牢山和无量山脉附近，且呈现出组团集聚特征. Ⅰ级空间贫困陷阱主要涉及和平镇丫口村、光山村，振太镇秀山村、帮庆村等地，上述区域农业生产条件较差，人力资源本底规模小、人口发展能力弱，土地资源本底条件差、生态承载能力弱，且产业空间区位偏远，产业收益保障低下，“人—地—业”多维结构呈现出失调状态，导致该区域减贫难度大，脱贫稳定性差，返贫风险强且制约因素复杂多样；Ⅱ级空间贫困陷阱主要涉及和平镇那洛村、麻洋村，振太镇黄力村、台头村等地；Ⅲ级空间贫困陷阱主要涉及振太镇文兴村、长安村等地. 上述区域人均土地资源和人均产业效益相对较低，人口发展潜力不足，土地功能尚未得到充分挖掘，产业发展潜力不足，导致该区域脱贫稳定性较差，存在一定的返贫风险.

本文在已有研究成果和乡村贫困地域系统理论探析的基础上，以乡村贫困地域系统“人—地—业”多维结构为基础，同时以“基底—支撑—潜力”递进式逻辑为支撑，构建了乡村贫困地域系统演化评价指标体系；设计了基于生计9边形的乡村贫困地域系统演化发展指数测度模型，得到了镇沅县2013-2021年109个评价单元的乡村贫困地域系统演化发展指数，并采用自然断点法绘制了乡村贫困地域系统时空演化水平分级图，进行了时空演化静态冷热点分析和基于时空立方体—新兴时空热点分析的动态时空演化模式识别，在此基础上还开展了基于“水平分级+演化模式”的空间贫困陷阱识别瞄准. 主要研究结论如下：

1) 镇沅县2013-2021年乡村贫困地域系统演化总体趋于中等水平，呈现出明显的时序阶段性和空间异质性特征. 时序阶段性特征主要表现为：2013-2015年为差异波动式演化阶段，2016-2019年为普惠式快速演化阶段，2020-2021年为增长滞缓趋稳阶段. 空间异质性主要表现为：时空演化热点空间呈现出中心集聚特征，主要集中在人口密度相对较大、耕地资源良好、距县域中心较近或产业相对发达区域，时空演化冷点空间呈现出分散集聚特征，主要集中在人口密度相对较低、距离县域中心相对较远、耕地资源较为贫瘠区域，且多分布于哀牢山和无量山等坡度起伏度较大的山地区域.

2) 镇沅县乡村贫困地域系统时空演化模式以连续热点模式为主，加强热点模式呈零星分布，表征了2013-2021年镇沅县乡村贫困地域系统整体上经历了较快且连续稳定的减贫演化过程，是减贫政策普惠作用在时空演化模式上的具体表现. 连续冷点模式和持续冷点模式表征了县域东西两端呈集聚分布态势，印证了区域乡村贫困地域系统演化过程深受自然地理条件制约和胁迫. 通过探索和设计乡村贫困地域时空演化模式分析方法，将“时空立方体—新兴时空热点分析”技术运用于乡村贫困地域系统时空演化模式识别，实现了多维时空乡村地域演化状态数据的叠加，可以更加直观地探究乡村贫困地域时空演化特征.

3) 通过构建涵盖乡村贫困地域系统演化水平和时空演化模式的“空间贫困陷阱”识别刻画方式，进一步瞄准了镇沅县空间贫困陷阱分布态势，相关研究方法为区域返贫风险的时空瞄准提供了决策参考.

脱贫地区在消除绝对贫困、巩固拓展脱贫成果与乡村振兴有效衔接的减贫过程中，其时空演化既具有自然周期性特征，又兼具跨周期或逆周期现象. 乡村发展是乡村地域系统循环累积与动态演化的结果，乡村贫困地域系统作为乡村地域系统的一种特殊状态，其转型发展具有时序演进和空间异质双重特征，梳理及评估其时空演化状态特征，对科学推动脱贫地区转型发展具有重要意义. 未来的研究在分析刻画乡村贫困地域系统演化发展的基础上，应进一步探究乡村贫困地域系统转型机理及科学途径，深入挖掘不同类型乡村转型发展的内生动力及协同机制，深化乡村人地系统耦合过程与情景模拟研究，充分发挥地理学科优势，重点关注全球减贫和乡村演化发展时空格局、相对贫困治理机制与成效提升策略、城乡融合与乡村振兴互馈机制等科学命题，为乡村振兴实践提供决策参考.