研究使用了核桃举肢蛾实际点位数据、生物学气候数据和海拔数据。核桃举肢蛾实际点位数据来自文献收集及实地调研，为避免区域差异引起的生态位狭窄，此部分数据选用中国范围，剔除密集点位后获得41个实际分布点位。生物学气候数据由2000—2020年中国地表气候月数据集的月温度和降雨量计算而得，包括代表年趋势（如年均温、年降水量）、季节性（如温度或降水量的年较差）、极端或限制性（如最冷或最暖月温度、湿季或干季降水量）的19个因子。海拔数据从中国科学院地理空间数据云（http://www.gscloud.cn）的DEM高程数据中提取。

常用的模型准确性评价指标很多，例如精准度（Precision）、召回率（Recall）、混淆矩阵（Confusion Matrix）、受试者工作特征曲线（ROC）及受试者工作特征曲线下面积（AUC）等。相较其他评价指标，AUC是一种整体上的衡量方法，不需要手动设置阈值，其在样本不均衡情况下仍能合理评价模型准确性，研究遂选用此方法对模型进行准确性评估。AUC的取值范围为0~1，数值越接近1说明预测效果越好。

参照赵金鹏^[10]等对物种分布可能性（Probability, P）的划分标准，将核桃举肢蛾的适生区分为四个等级，即高适生区（P≥0.66）、中适生区（0.33≤P<0.66）、低适生区（0.05≤P<0.33）和不适生区P<0.05）。

秦巴山区西起青藏高原东缘，东至华北平原西南部，区域范围包括河南、湖北、重庆、四川、陕西和甘肃六省市的75个县（区），其中河南10个、湖北7个、重庆5个、四川15个、陕西29个和甘肃9个，国土总面积约为22.5平方公里。

基于核桃举肢蛾实际点位数据、生物学气候数据和海拔数据构建核桃举肢蛾潜在地理分布初始预测模型，依据贡献率和置换重要性筛选出8个影响核桃举肢蛾潜在地理分布的主导环境因子（表1）。从表1中可知，8个主导环境因子依次为最湿月降水量、最干季平均温度、最冷月最低温度、海拔、最冷季降水量、等温性、最干季降水量和最冷季平均温度（按照贡献率大小排序），累计贡献率和置换重要性分别为93.3%和90.1%，可以用来描述核桃举肢蛾环境需求特点，故选取以上8个因子作为影响核桃举肢蛾分布的主导环境因子建立最终预测模型。

从核桃举肢蛾对主导环境因子的响应情况（图2）可知，对于8个主导环境因子，核桃举肢蛾存在概率均是先随着主导环境因子值的增大而升高，达到一定峰值后，存在概率随着主导环境因子值的增加而降低。以存在概率0.33及以上为研究对象，最湿月降水量、最干季平均温度、最冷月最低温度、海拔、最冷季降水量、等温性、最干季降水量和最冷季平均温度的范围分别为110.5~272.7 mm、−4.8~8.8 ℃、−13~4 ℃、90.5~2154.5 m、8.5~97.9 mm、24.2~33.8、15.5~87.9 mm和−4.7~7.8 ℃，最佳值为186.09 mm、1.9 ℃、−3.7 ℃、353.1 m、28.2 mm、27.5、36.1 mm和1.6 ℃。

图  2  核桃举肢蛾预测模型的ROC曲线图

Figure 2.  ROC curve of the prediction model for A. hetaohei

研究利用测试数据集对构建的预测模型进行了准确性检验，并绘制了受试者工作特征（ROC）曲线图（图1）。从图1可知，重复测试的AUC平均为0.948，标准偏差为0.008。根据模型准确性评判标准，预测模型能较好模拟核桃举肢蛾在秦巴山区的分布。

图  1  不同主导环境因子状态下核桃举肢蛾的存在概率

Figure 1.  Existence probability of A. hetaohei existence under different dominant environmental factors

核桃举肢蛾在秦巴山区分布较广，总适生区面积达19.63×10⁴ km²。根据当前情景下核桃举肢蛾在秦巴山区的潜在适生面积（表2）和潜在分布图（图3），核桃举肢蛾高适生区主要集中在秦巴山区的中部，河南、四川和重庆等6省市均有分布，其中甘肃省面积最小（0.13×10⁴ km²），主要分布在甘肃与陕西交界的地区，包括两当县、徽县、成县、武都县、康县和文县；陕西省面积最大（4.8×10⁴ km²），除周至县和太白县外，陕西南部的17个县（区）均有分布。不适生区面积约为1.6×10⁴ km²，主要集中在秦巴山区的西部和东南部，即四川的北川、平武两县，甘肃的宕昌、文、礼三县及河南的淅川和镇平两县。

图  3  当前情景下核桃举肢蛾在秦巴山区的潜在分布图

Figure 3.  Potential suitable distribution of A. hetaohei in Qinba Mountain area under current situation

目前对分类预测模型评价指标有很多，例如精确率、召回率、F值和AUC值（ROC曲线）。相较于其他方法，AUC因不受诊断阈值的影响，且能提供所有阈值范围上的性能评价结果，被公认为生态位分类预测模的最佳评价指标。基于选取的主导环境变量重复测试的AUC平均为0.948，标准偏差为0.008，模型预测结果达到“极好”标准, 说明建立的预测模型能够很好地模拟核桃举肢蛾在秦巴山区分布。

物种与环境之间的关系是研究物种生态学和空间分布的重要方面之一，对核桃举肢蛾对主导环境因子的响应情况进行了分析，绘制了相应的响应曲线。研究结果表明，最湿月降水量、最干季平均温度、最冷月最低温度、海拔、最冷季降水量、等温性、最干季降水量和最冷季平均温度8个环境变量是主导核桃举肢蛾在秦巴山区分布的环境因子，且对于8个主导环境因子，核桃举肢蛾存在概率均是先随着主导环境因子值的增大而升高，达到一定峰值后，存在概率随着主导环境因子值的增加而降低。查阅相关文献，各地区的核桃举肢蛾普遍以老熟幼虫在树冠下1~3厘米深的表土层中、石块下及树缝中、枯果内结茧越冬。越冬幼虫普遍于 5 月上旬开始出土化蛹，5 月中旬或6月上旬为出土化蛹盛期，而成虫主要发生在 5 月下旬至7 月下旬，8 月上旬或8月底第2代老熟幼虫进入脱果盛期，9 月至翌年 3 月以老熟幼虫结茧越冬^[11-14]。由此可见，最冷月最低温度、最冷季降水量和最冷季平均温度关系着核桃举肢蛾是否能够安全过冬，而最湿月降水量、最干季平均温度和最干季降水量则关系着核桃举肢蛾是否能够发育繁殖。研究进一步分析表明核桃举肢蛾中适生区以上区域的最湿月降水量、最干季平均温度、最冷月最低温度、最冷季降水量、最干季降水量和最冷季平均温度范围分别为110.5~272.7 mm、−4.8~8.8 ℃、−13~4℃、8.5~97.9 mm、15.5~87.9 mm和−4.7~7.8 ℃，最佳值为186.09 mm、1.9 ℃、−3.7 ℃、28.2 mm、36.1 mm和1.6 ℃。

核桃举肢蛾在秦巴山区分布较广，总适生区面积达19.63×10⁴ km²。高适生区主要集中在秦巴山区的中部，其中甘肃省面积最小，约为0.13×10⁴ km²，包括两当县、徽县、成县、武都县、康县和文县；陕西省面积最大，约为4.8×10⁴ km²，除周至县和太白县外，陕西南部的17个县（区）均有分布。由此可看出，核桃举肢蛾在秦巴山区危险系数极高, 有关部门应引起高度重视, 严防核桃举肢蛾造成重大损失。昆虫的潜在地理分布不仅受气候和海拔因素的影响，也受地区的经济、人口等因素影响，但研究未考虑此类情况；同时，在统计核桃举肢蛾适生面积过程中，未扣除水体、城镇等核桃举肢蛾非生长区。以上两个问题，不可避免地会影响预测结果的准确性，未来研究将重点关注上述两个问题，进一步提高预测模型结果的准确性。