《Ecological Indicators》：Weighted Jaccard similarity: A threshold-free multivariate indicator of correlation-network reorganization across coral, freshwater, and climate systems
编辑推荐：
生态临界点的检测始终是环境监测的核心挑战。基于临界减速（Critical Slowing Down, CSD）理论的经典早期预警信号仅作用于单变量时间序列，在结构层面对多变量相关性重组缺乏敏感性。加权Jaccard（Weighted Jaccard, WJ）相
生态临界点的检测始终是环境监测的核心挑战。基于临界减速（Critical Slowing Down, CSD）理论的经典早期预警信号仅作用于单变量时间序列，在结构层面对多变量相关性重组缺乏敏感性。加权Jaccard（Weighted Jaccard, WJ）相似性指数已被提议作为此类重组的无阈值指标。本研究针对前期投稿中审稿人提出的关于独立同分布（Independent Identically Distributed, IID）置换推断对存在空间及时间依赖性的生态数据适用性的关切进行了修订。研究人员在底物匹配的结构保持零假设下重新分析了三组经典生态临界点数据集：全球珊瑚礁白化（全球珊瑚白化数据库，41,361条记录，36个环境变量，1980–2020年）、美国湖泊富营养化（美国环保署国家湖泊评估2007、2012、2022年，23个水化学变量，三次独立调查跨越15年）以及全球气候遥相关网络（19个月尺度的美国国家海洋和大气管理局指数，1950–2024年）。在所有三种底物中，WJ值均具稳健性和可重复性；其统计显著性取决于结构保持零假设所应用的聚类尺度，而该聚类尺度本身即作为参数化测量轴发挥作用。在珊瑚礁系统中，结构重组信号在所有窗口与方法组合下均经受住了最精细空间零假设的检验（站点级聚类K≈10,000个站点）。在湖泊系统中，信号在标准的9区域环保署生态区聚类下减弱，但在更精细的分辨率下趋于稳定；营养梯度WJ特征在所有三次独立调查中均得到重复。在气候系统中，信号在6个月块置换下消失，且在多变量相位随机化替代数据检验下与零假设无法区分。三种底物指纹维度（连续-离散Type 1间隙、符号处理Type 2间隙、底物投影Type 6间隙以及直接成对水平符号翻转率）在不同底物间具有正式可区分性（非对称Tversky Type 3间隙的Kruskal-Wallis p<10?3）。在有意义的相关幅度（|r|≥0.2）下，符号保持具有跨底物的普遍性（符号翻转率<3%）。研究结果将WJ定位为具有明确推断边界条件的生态监测底物指纹指标，而非普遍显著的临界点探测器。WJ的适用条件以及依赖结构被恰当建模的底物特异性聚类尺度，均从应用该指标的同一数据集中经验性地导出。
生态系统的临界点检测是当代环境科学的核心议题。当渐进的环境变化触发系统结构与功能的快速重组时，生态系统可能发生突变的非线性转变。此类转变涵盖从浅水湖泊浊度变化到珊瑚礁相态转换，再到大规模气候系统重组的多重尺度。基于临界减速理论的经典早期预警信号通过监测单变量时间序列的方差与自相关上升来预测系统趋近分岔点，然而这类指标在结构上对多变量重组不敏感——即变量间耦合关系的重新配置可能构成系统压力的互补性标志。例如，珊瑚礁对热胁迫的响应涉及温度异常、加热周数、热胁迫频率、风驱动环流及生物覆盖度等多因素的协同变化；湖泊富营养化伴随着养分输入、浊度、叶绿素、电导率及地球化学缓冲之间耦合关系的重组；气候系统突变则表现为热带太平洋振荡、大气环流模态与大西洋多年代际变率之间遥相关耦合的改变。因此，仅依赖单变量监测的框架可能错失集中于多变量耦合结构的信号。
相关网络比较方法通过构建不同条件下的相关矩阵并量化其相似性来弥补这一空白。传统方法以固定阈值（通常为绝对相关性的前5%）将矩阵二值化后使用Jaccard指数比较，但该方法存在两大局限：阈值由研究者主观设定且结果对选择敏感；二值化丢弃了幅度信息，相关系数0.95与0.10若均高于阈值则被等同处理。加权Jaccard相似性指数通过对连续相关幅度向量运算避免了这些局限，其取值范围为0（完全相异的幅度剖面）至1（完全相同的剖面），具有无阈值构造特性且不受变量数量缩放的影响。
前期投稿针对珊瑚礁白化、湖泊富营养化和气候遥相关三种经典生态临界点数据集应用WJ，在假设观测独立同分布的置换检验下报告了显著的相关网络重组。审稿人提出了两项实质性关切：其一，IID置换不适用于空间聚类数据（珊瑚礁站点、湖泊）和时间序列相关数据（月尺度气候指数）；其二，基于|r|的未签名WJ将r=+0.8与r=?0.8等同处理，因而无法识别相关符号反转——即配对水平上最大可能的重组形式。
本研究针对上述关切扩展了分析框架，以底物匹配的结构保持零假设替代IID置换作为主要检验方法：对珊瑚礁采用站点聚类置换，对湖泊采用多尺度生态区聚类置换，对气候系统采用块置换结合多变量相位随机化替代数据。同时，系统变化聚类尺度以产生WJ显著性的粒度轨迹，该轨迹本身即作为底物依赖性尺度的测量。此外，研究加入了WJ家族的符号感知扩展（非对称Tversky、符号分裂、底物投影）以及直接成对水平符号翻转率测量，形成可进行跨底物正式区分性检验的底物指纹维度。研究的实证贡献框架相应调整：不再声称统一的临界点检测，而是报告WJ信号在不同底物中经受或未能经受结构保持零假设推断的特异性边界条件、正式区分三种系统的底物指纹，以及各底物中依赖性运作的经验尺度。
研究所用数据包括：全球珊瑚白化数据库的41,361条调查记录（1980–2020年，36个分析变量），按时间划分为基线期（1980–2009）、升级期（2010–2015）、大规模白化期（2016–2017）和后期（2018–2020）；美国环保署国家湖泊评估的三次独立调查数据（2007、2012、2022年，23个水化学参数），按叶绿素-a浓度分为贫营养、中营养、富营养和超富营养四种状态；美国国家海洋和大气管理局的19个月度气候指数（1950–2024年，171对相关性），按已确立的气候突变边界分为基线期（1950–1976）、后1977期（1977–1997）、后1998期（1998–2009）和近期加速期（2010–2024）。
研究构建了相关向量并计算了三种WJ公式：基于|r|的未签名WJ、基于符号分裂向量的符号感知WJ，以及基于前5%相关配对的二值化Jaccard。研究还计算了四种配对家族间隙测量：Type 1（连续-离散间隙，未签名WJ减二值化Jaccard）用于定位沿值分布维度的发散；Type 2（符号处理间隙，未签名WJ减符号感知WJ）量化符号驱动的不相似性贡献；Type 3（非对称-对称间隙，非对称Tversky覆盖度之差）指示结构重组的方向性（相对基线的架构扩张或收缩）；Type 6（底物投影间隙，Pearson未签名WJ减Spearman未签名WJ）指示线性结构与秩单调结构在底物配对关系中的相对贡献。此外，研究补充了直接成对水平符号翻转率测量，并按最小绝对相关幅度分层检验。
主要技术方法包括：底物匹配的结构保持零假设生成（珊瑚礁站点聚类置换、湖泊多尺度生态区聚类置换、气候系统多尺度块置换及多变量相位随机化替代数据），各生成1,000次置换；聚类粒度轨迹分析，将聚类尺度作为参数化测量轴报告WJ z分数的轨迹；跨底物正式检验（Kruskal-Wallis检验底物来源标签的可区分性，Bootstrap 95%置信区间）；以及成对差异的Fisher z变换与Benjamini-Hochberg错误发现率校正。分析采用Python进行，全部代码公开可得。
**珊瑚礁白化结果**：未签名WJ显示，与1980–2009基线相比，相关网络相似性在后续三个窗口中渐进下降。在IID置换下所有比较均达p<0.0001显著性；在站点聚类置换下（K=9,388至11,712个站点），所有窗口与方法组合仍达p<0.0001显著性，表明珊瑚礁结构重组信号经受住了最精细空间零假设的检验。聚类粒度轨迹呈非单调形态：在空间粒度较粗时（K≤481）观测WJ超过聚类置换零假设（z分数+1.4至+3.4），仅在站点级粒度时观测WJ低于零假设（z=?6至?14），这一非单调形状构成了珊瑚礁底物指纹。Type 6底物投影间隙在所有三个比较中均为正值，表明Pearson检测到比Spearman更多的相似性，与珊瑚环境关系中的非线性单调结构一致。Type 3非对称Tversky间隙在三个窗口中符号变化（升级期正值对应架构扩张，大规模白化期负值对应架构收缩，后期混合）。
**湖泊富营养化结果**：未签名WJ的营养梯度特征在2007、2012、2022年三次独立调查中均得到重复，跨调查标准差（约0.03）比效应量小一个数量级。在标准环保署9聚合生态区尺度（K=9）聚类置换下，大多数比较未达显著性；但聚类粒度轨迹显示清晰单调模式：随粒度变精细z分数加深（Spearman ρ=?0.522，p<10
?6）。在HUC2流域聚类（K≈18–21）时中位数z分数为?2.0（边缘显著），在HUC8聚类（K≈300–1150）时为?2.8。Type 3非对称Tversky间隙在所有18个湖泊比较中均为负值（符号检验18/18，p<10
?5），与富营养化态相比问题的-效基线的系统性耦合丧失一致，潭泊是三种底物中Type 3间隙方向唯一具有正式一致性的底物。Type 6底物投影间隙为负值，与湖泊水化学关系中的异常值影响线性结构一致。同一营养状态跨调查的建筑稳定性WJ值（0.67–0.71）表明审核的营养相关性网络在15年间具有中度稳定性。
**气候遥网络结果**：在IID置换下所有三个时代比较均达p<0.0001显著性。但在底物匹配块置换下，信号在块大小6个月时消失：Pearson z分数轨迹从块大小1个月的?4.95降至6个月的?0.97，之后围绕零值波动（Spearman ρ=+0.780，p<10
?8）。在多变量相位随机化替代数据检验下，三个时代比较均未达显著性（z分数+0.10至+0.71，p=0.53–0.75）。Type 6底物投影间隙接近零，与月尺度气候指数间的近似线性关系一致。Type 3间隙在三个时代中符号变化（后1977期负值，后1998期和近期正值）。
**直接成对水平符号翻转率结果**：在有意义相关幅度（|r|≥0.2）下，符号保持具有跨底物普遍性。符号翻转率为：珊瑚礁0.14%、湖泊1.23%、气候1.04%，均拒绝50%随机率（p<0.001），但跨底物间无显著差异（Kruskal-Wallis p=0.70）。这一结果 empirically 回应了审稿人关于WJ对符号反转不敏感的关切：在这三种生态底物中，有意义相关幅度下的符号翻转普遍罕见。
**跨底物正式检验结果**：Type 6间隙在三种底物间可区分（Kruskal-Wallis H=6.80，p=0.033），各成对Bootstrap 95%置信区间互不重叠（珊瑚礁+0.045、湖泊?0.090、气候?0.004），分别对应非线性单调、异常值影响线性和近似线性结构。Type 3间隙区分性更强（H=17.37，p=1.7×10
?4），湖泊与珊瑚礁和气候均可区分。符号翻转率不可区分（p=0.70）。推断间隙（IID z减结构保持z）在气候与空间底物间可区分（气候?5.1 vs 珊瑚礁?19.4/湖泊?23.2），与空间底物相比气候的依赖性结构更易被零假设模型化。
讨论部分首先重新框架了实证贡献，将前期声称的"统一临界点检测"调整为报告"WJ信号在不同底物特异性推断条件下的存废"。这一调整在方法上更为严谨：WJ作为可用指标保留于三种底物中，而依赖性恰当建模的底物特异性聚类尺度现从同一数据集中经验导出，而非假设为研究者设定的粒度。
聚类粒度作为测量轴是本研究的方法论核心。历史应用中聚类尺度需由研究者论证，本研究将其视为参数化测量轴。湖泊的单调深化轨迹（ρ=?0.522）揭示了标准K=9生态区尺度仅为连续体上信号可检测但较弱的一点；气候的单调松弛轨迹（ρ=+0.780）揭示了月指数自相关尺度约为6个月——这些均从数据本身经验提取而非假设。珊瑚礁的非单调轨迹（粗粒度正z分数，细粒度负z分数）则构成其特有的底物指纹，与湖泊和气候的单调模式形成对比。
底物指纹的解读表明，Type 6间隙的三种签名可解释：珊瑚礁的非线性单调结构、湖泊的异常值影响线性结构、气候的近似线性结构。Type 3间隙中湖泊的统一收缩方向（18/18负值）最为一致，与浅水湖泊动力学的双稳态框架相符：富营养化态伴随着相比浮游-富营养基线的系统性变量耦合丧失，而非简单的耦合重分布。
直接符号翻转率的普遍低发生率（<3%）表明，在这些底物中生态相关网络重组主要通过网络内主要经度变化而非符号反转，使未签名WJ的符号盲点在经验上成为非关键问题。
各底物的生态解读如下：珊瑚礁分析与前认识一致，WJ渐进下降反映2016–2017全球白化事件后环境-生物耦合结构的累积重组，信号对站点聚类零假设的稳健性表明其可检测于站点间环境变异之上。湖泊分析中三次独立调查的营养梯度 signature 重复提供了强经验支持，表明其为体系特性而非调查特定伪影，聚类粒度轨迹揭示监控项目应报告轨迹而非依赖单一标准生态区检验。气候分析呈现最清晰的推断边界案例：月尺度可检测，但吸收于≥6个月的自相关结构中，这并非否定既有文献中后1976气候突变的存在，而是明确了WJ指标适用的时间聚合尺度。
研究将WJ置于早期预警信号的更广泛背景下进行讨论。WJ捕获的信号与经典CSD指标不同：CSD检测单个变量趋近分岔时的恢复速率减慢，WJ检测多变量耦合结构的重组，两者互补而非竞争。WJ属于非CSD类别的结构性相似度量，不假设转变类型（ fold、Turing、噪声诱导或速率诱导转变）。
研究承认了若干局限性：气候分析未通过多变量相位随机化替代检验，应视为底物特异性推断边界而非方法失败；未签名WJ的符号盲点已通过Layer 2I直接测量弥补；聚类粒度框架依赖多尺度结构的可用性；结构保持零假设内部的处理仍待空间块自举等更严谨方法完善；相关性结果不隐含因果性。WJ对监测项目的实用优势在于无需阈值选择、不受变量数量影响、可从标准多变量数据集计算，且配对家族分解提供多个底物指纹维度。
研究结论如下：加权Jaccard相关网络相似性在三种经典生态临界点系统中作为多变量生态重组指标具有可用性，但其显著性底物特异性地依赖于结构保持零假设的条件。珊瑚礁信号经受最精细空间零假设，湖泊信号在标准生态区尺度减弱但在更精细尺度稳定，气候信号在约6个月自相关尺度下消失。聚类尺度作为参数化测量轴的轨迹提供底物特异性依赖尺度诊断。Type 3和Type 6签名在跨底物间正式可区分。有意义相关幅度下的符号翻转率普遍低（<3%）。WJ应定位为具有明确推断边界条件的生态监测底物指纹指标，而非普遍显著的临界点探测器。其适用条件从应用该指标的同一数据集中经验导出，解决了审稿人关于推断严谨性的实质性关切。
译文（研究结论部分）：
加权Jaccard相关网络相似性作为多变量生态重组指标，在底物匹配的结构保持零假设下，于三种经典临界点系统——珊瑚礁白化、湖泊富营养化和气候遥相关突变——中仍具可用性Semantic。该信号在珊瑚礁中经受最精细空间零假设的检验，在湖泊中于标准9区域环保署生态区零假设下减弱（并在更精细尺度趋于稳定），在气候中因时间自相关而消失（经验自相关尺度约为6个月）。结构保持的聚类尺度是参数化测量轴，其轨迹提供底物特异性依赖尺度诊断。四种配对家族间隙测量（Type 1连续-离散、Type 2符号处理、Type 3非对称Tversky、Type 6底物投影）及一种直接成对测量（Layer 2I符号翻转率）产生底物指纹；其中Type 3和Type 6签名在三种底物间具有正式可区分性（Kruskal-Wallis p<0.05及p<10
?3）。在有意义相关幅度下，符号翻转率于三种底物中均普遍低（|r|≥0.2时<3%）。该框架将WJ定位为具有明确推断边界条件的生态监测底物指纹指标，而非普遍显著的临界点探测器。WJ的适用条件从应用该指标的同一数据集中经验导出，通过扩展方法论而非修补单一推断步骤，回应了审稿人关于推断严谨性的实质性关切。
