地中海盆地五种 Anthidium 羊毛蜂的物候趋势与未来气候适宜性：非平稳性、采样偏差与空间验证

《Insects》：Phenological Trends and Projected Climatic Suitability of Five Wool Carder Bees of the Genus Anthidium Across the Mediterranean Basin

　　

气候变化可能改变蜜蜂活动的时间以及蜜蜂栖息地的气候适宜性，增加植物—传粉者错配风险，但此类影响鲜少针对地中海蜜蜂进行评估。研究人员分析了5种羊毛蜂（Anthidium）的4103条分布记录（1980—2024），以评估首飞时间趋势，并利用由5个通用环流模型（general circulation models, GCMs）驱动的集合物种分布模型（ensemble species distribution models, SDMs）预测2081—2100年气候适宜区的变化。研究人员在4个物种中检测到非平稳物候动态。仅A. florentinum表现出稳健的首飞时间提前，为每十年10.0天，且经错误发现率（False Discovery Rate, FDR）校正后仍显著（q < 0.001），而A. loti的提前不显著（q = 0.082）。仅分析伊比利亚记录时，A. florentinum的提前仍持续存在。在SSP5-8.5排放情景下，15个国家中有6个满足数据充分性标准；其中，气候适宜采样地点的最大净下降预测出现在西班牙（?15.3个百分点）、摩洛哥（?13.7）和葡萄牙（?12.5），而意大利、法国和克罗地亚显示小幅净增加。在空间区块交叉验证（spatial block cross-validation）下，模型表现低于随机交叉验证，但按曲线下面积（area under the curve, AUC）仍可接受（0.75–0.81）。在满足数据充分性标准的国家中，摩洛哥、西班牙和葡萄牙超过了复合脆弱性指数（Composite Vulnerability Index, CVI）的参考阈值，该指数综合了预测栖息地下降、物候偏移和预测增温。这些结果表明存在不均衡的、国家特异性的脆弱性，并强调在评估地中海传粉者气候风险时需考虑非平稳性、采样偏差和空间验证。

研究背景与问题：昆虫物候（phenology，生活史事件发生时间）受光周期和温度驱动，温度升高可加速发育并导致季节性发生提前。蜂类（Hymenoptera: Apoidea）是野生和栽培植物的重要传粉者，气候变暖可能破坏植物—传粉者物候同步，且植物与动物的物候响应可能日益不对称。除物候外，气候变化还会重塑传粉者地理分布；但针对同一分类群，时间与分布两个维度很少被联合评估。地中海盆地升温速度高于全球许多区域，是气候变化热点和生物多样性热点，野生蜂多样性高，但关于该区域蜂类物候响应的知识仍很有限。已有发生记录在物种、年份和国家间样本量与可靠性差异很大，且部分物种存在分类合并问题，因此需要同时评估趋势稳健性、采样偏差、非平稳性和空间验证。研究人员据此选择5种较易识别且具有近几十年历史数据的羊毛蜂（Anthidium），在同一分析框架中检验首飞时间是否提前、趋势是否随时间稳定、与春季生长度日的关系，并预测本世纪末气候适宜栖息地变化。

研究人员开展的研究与主要结论：研究人员整合1980—2024年地中海盆地5种Anthidium的4103条发生记录，采用年内第10百分位日序（P₁₀ DOY）表征首飞起点，结合趋势检验、结构断点、气候—物候关系、集合物种分布模型（SDM，species distribution model）和复合脆弱性指数（CVI，Composite Vulnerability Index），评估物候趋势与2081—2100年气候适宜性。结果显示，仅A. florentinum具有经FDR校正后稳健的首飞提前（?10.0天/十年，q < 0.001）；A. loti s.l.为提示性趋势；其余3种无统计显著趋势但点估计方向一致。4/5物种呈非平稳动态，单一直线趋势推断受限。栖息地预测中，满足数据充分性标准的6国里，西班牙、摩洛哥和葡萄牙净下降最大，法国、意大利和克罗地亚小幅净增加；CVI显示摩洛哥、西班牙和葡萄牙超过近似高风险参考阈值，葡萄牙仅为边际超过。研究意义在于首次对Anthidium属开展整合物候—分布的风险评估，并强调非平稳性、采样偏差和空间验证在地中海传粉者气候风险评估中的重要性。

关键技术方法（不超过250字）：研究人员使用全球生物多样性信息网络（GBIF）1980—2024年5种Anthidium的4103条记录，并与EBHD数据库交叉核对；研究范围覆盖地中海盆地23个候选国家，15国含确认记录，北马其顿1条记录仅用于物候和分布建模。以P₁₀ DOY表征首飞起点，采用Mann–Kendall趋势检验、Theil–Sen斜率、Chow结构断点检验和线性/广义可加模型（GAM，generalised additive model）分析；以ERA5计算春季生长度日（GDD，growing degree days），以WorldClim变量经方差膨胀因子（VIF，variance inflation factor）筛选建模；构建随机森林（RF）、梯度提升（GB）和逻辑回归（LR）的集合SDM，采用随机与1°×1°空间区块交叉验证，并用5个GCM在SSP2-4.5、SSP5-8.5下预测；计算CVI和物候债务钟（PDC，Phenological Debt Clock）。

研究结果：
3.1 五种Anthidium物候偏移：基于P₁₀时间序列的Theil–Sen趋势和Mann–Kendall检验显示，仅A. florentinum在FDR校正后显著提前（?10.0天/十年，q < 0.001）；A. loti s.l.名义显著（?12.2天/十年，p = 0.033）但FDR后不显著（q = 0.082）；A. manicatum、A. oblongatum、A. punctatum无显著趋势但均方向一致。4/5物种斜率前后半段偏差超过50%，判为非平稳。更严格年样本阈值下仅A. florentinum稳定显著；Weibull第10百分位亦支持其提前，而年均观测日期无显著趋势；A. florentinum在伊比利亚记录及嵌套西部子集中仍显示提前。
3.2 线性与GAM气候—物候响应比较：GAM在小样本（n = 11–22）下严重过拟合，交叉验证R²为负（?0.29至?59.0）；线性模型样本内R²较低（0.00–0.25）。仅A. florentinum的局地GDD—DOY偏相关名义显著（偏r = ?0.47，未校正p = 0.028），但FDR校正后不显著（q = 0.140），且对GDD设定敏感；因此推断仅基于线性模型。
3.3 群落水平SDM验证：随机交叉验证下AUC为0.810–0.836、TSS为0.480–0.556；空间区块交叉验证下AUC降至0.748–0.809、TSS降至0.374–0.458，所有模型超过AUC阈值，除A. oblongatum外均超过TSS阈值。随机交叉验证高估空间外推能力。MESS（多变量环境相似面）分析显示，SSP2-4.5和SSP5-8.5下分别有7.7%和26.0%背景格点至少在一个GCM中处于校准气候范围之外，主要位于北非、东地中海和北部亚得里亚—巴尔干区域。
3.4 A. florentinum首飞趋势短期外推至2050：将Theil–Sen趋势自2012年参考点线性外推，P₁₀由2012年DOY 164.6（6月14日）移至2030年DOY 146.6（5月27日）和2050年DOY 126.6（5月7日），提前38天；该外推仅为趋势延续的理论示意，不是生态预测。
3.5 SSP5-8.5下国家与物种净栖息地变化：在15个有确认记录国家中，6国满足数据充分性标准。物种平均净下降最大为西班牙（?15.3个百分点）、摩洛哥（?13.7）和葡萄牙（?12.5）；意大利（+1.7）、法国（+8.0）和克罗地亚（+8.2）显示小幅净增加。西班牙、葡萄牙和摩洛哥在SSP2-4.5下亦为净下降，方向一致。其余9国因记录少或基线适宜面积小被标记为仅方向性指示。
3.6 国家水平CVI：等权CVI显示摩洛哥（0.74）、西班牙（0.72）和葡萄牙（0.67）超过近似高风险参考阈值（≈0.67），葡萄牙仅边际超过；土耳其（0.70）亦超过但未满足数据充分性标准。摩洛哥、西班牙和葡萄牙高分主要由栖息地下降（D1）及其与物候偏移的交互（D2）驱动；意大利和克罗地亚较高分主要由预测增温（D3）驱动。以ERA5观测增温趋势替代D3时，同一四国超过相应阈值，国家得分高度相关（Spearman ρ = 0.96）。
3.7 VIF逐步筛选与相关矩阵：对19个WorldClim BIO变量进行逐步VIF剔除，保留8个变量（bio1、bio2、bio4、bio8、bio9、bio14、bio15、bio19），最终最大VIF为7.49（bio1）。保留变量中bio14与bio15相关性最高（r = ?0.78），但二者VIF均低于10而保留；残余共线性可能增加逻辑回归组分的系数方差。
3.8 五种Anthidium物候债务率：经FDR校正后无物种达到CONFIRMED层级。A. florentinum为PARTIAL层级：物候趋势FDR显著，而GDD敏感性仅名义显著；观测提前?10.0天/十年超过局地GDD敏感性预期提前?4.0天/十年，债务率为?6.0天/十年。其余4种为EXPLORATORY层级（债务率?3.8至?15.1天/十年），因缺乏稳健基础信号仅作示例。PDC为概念性、假设生成指标，不用于确证推断。
3.9 CVI对组分权重的敏感性：等权方案下摩洛哥、西班牙、土耳其和葡萄牙超过≈0.67参考阈值；栖息地加权方案下同一四国仍超阈值，前五名国家与等权一致；温度加权方案下土耳其第一、摩洛哥第二，仅两国超阈值，西班牙降至第三且低于阈值，葡萄牙降至第十一。摩洛哥、西班牙、土耳其和意大利在三种方案下均居前五；葡萄牙位置依赖温度组分权重。无经验证据支持其他权重时，等权方案被保留为主要且方法学上可辩护的选择。

讨论总结与结论翻译：讨论指出，仅A. florentinum具有稳健且经FDR校正的首飞提前，其他物种趋势不显著或仅为提示性，且多数物种非平稳，说明单一直线趋势推断不适用于本数据集中的多数物种。P₁₀估计在样本稀疏年份精度下降，严格阈值分析显示唯一确证结果不依赖最稀疏年份；对其他物种进行线性外推会产生生物学上不合理的日期，因此外推仅限A. florentinum至2050年。GAM在小样本下过拟合，支持使用线性模型进行推断。集合SDM在随机交叉验证下表现可接受，但空间区块交叉验证下降，提示随机交叉验证高估空间迁移能力。栖息地净下降集中于西南地中海，法国、意大利和克罗地亚显示小幅净增益，呈向研究区北缘移动的格局。CVI高分国家主要由栖息地下降驱动，意大利和巴尔干部分国家主要由预测增温驱动；权重方案会影响葡萄牙等国的排名。研究还强调地理采样不平衡（西班牙、法国和葡萄牙合计贡献84%目标物种记录，北非记录很少）、北非和东地中海部分格点处于气候外推范围，以及A. florentinum和A. loti s.l.分类合并可能带来的偏差。物候提前本身不应直接解释为负面种群效应，研究也未确立变暖为原因，而将其描述为气候相关趋势。结论翻译：在5种地中海Anthidium物种中，A. florentinum的首飞时间表现出稳健的、经FDR校正的提前，在仅分析伊比利亚记录时仍存在；A. loti s.l.为提示性，其余物种未产生统计可靠的时间趋势。国家水平上，栖息地估计在6国满足数据充分性标准；净下降集中于西班牙、摩洛哥和葡萄牙，这些国家也超过本研究的近似CVI参考阈值，而法国、意大利和克罗地亚显示小幅净增加。证据因此指向不均衡的、物种和国家特异性的气候脆弱性，而非统一的地中海响应。尽管存在局限，本研究是Anthidium属首个整合物候—分布风险评价，将长期数据集、群落水平集合建模和多维脆弱性分析纳入同一方法框架，覆盖地中海盆地尺度。未来研究结合多种物候估计量、更大的多GCM集合、花资源数据、寄主植物物候互作、景观水平栖息地连通性和局地适应能力变量，将有助于确定所观察模式的稳健性并阐明其机制和生态后果。鉴于地中海盆地是对全球气候变化最敏感的区域之一且为全球生物多样性热点，这里记录的Anthidium蜂物候趋势和气候适宜性预测变化为该区域传粉者保护政策提供参考框架。