Appalachian遗留煤矿土壤肥力限制与探索性多变量分类：酸度连续梯度与独立痕量金属维度

《Soil Systems》：Soil Fertility Constraints and Exploratory Multivariate Classification of Appalachian Legacy Coal Mine Soils

　　

Appalachian煤田（Appalachian coalfield）的遗留煤矿地保留着化学异质性土壤，但一致的多站点比较仍缺乏。本研究探讨土壤性质是否在站点间存在差异、哪些性质构成样品水平变异，以及样品是否形成不同的土壤化学类别。研究人员在来自西弗吉尼亚州和宾夕法尼亚州7个遗留煤矿站点的135个表层（0–20 cm）样品中测定了18种土壤测试性质，并使用Kruskal–Wallis检验（Holm校正）、17种性质的主成分分析（principal component analysis, PCA）和Ward聚类进行分析，同时以bootstrap重采样和替代设定作为稳健性检验。所有18种性质在站点间均存在差异，尽管硼差异与实验室批次重合。五个成分解释了69.9%的方差；酸度/盐基阳离子、痕量金属和磷有效性成分可重复（Tucker一致性≥0.87），但两个次要成分不可重复。聚类将酸性、需石灰的多数（n = 106）与近中性少数（n = 29）分开，但分离较弱（轮廓系数=0.27），仅酸性组稳定（Jaccard=0.90 versus 0.67），且该划分匹配跨站点边界的pH阈值6.0。这些土壤更适合描述为具有独立痕量金属维度的连续酸度梯度，而非离散肥力类别，为设计但不替代复垦和生物能源物种的田间试验提供基线。

一、研究背景与问题
阿巴拉契亚地区地表采煤遗留大量扰动土地和剧烈改变土壤。自1970年代末以来，阿巴拉契亚煤田约2.5 million acres（约10,000 km²，约占>150,000 km²煤田的6.5%）受地表采煤扰动；1976—2015年中阿巴拉契亚累计地表采矿足迹约5900 km²（1.5 million acres），占该区7.1%。该遗留包括1977年《地表采矿控制与复垦法》（Surface Mining Control and Reclamation Act, SMCRA）后许可复垦土地和废弃矿地（abandoned mine lands, AMLs）。仅西弗吉尼亚州联邦清单列有5478个AML问题区、9790个记录问题，其中44%记录为复垦完成、8%正在复垦、48%待复垦。许多遗留景观仍保留退化土壤，持续限制植被恢复、土地生产力和下游生态系统服务。这些土壤更多受采矿暴露地质材料而非自然成土作用影响。地表采矿翻转覆盖层剖面，使未风化岩石到达地表；酸砂岩和黄铁矿页岩可形成强酸土壤（pH低于4），碳酸盐岩层可形成碱性土壤（pH近8），中间材料可随时间漂移。酸碱核算决定石灰需求和缓冲pH，历史上是西弗吉尼亚州和宾夕法尼亚州预采矿规划与复垦许可的主要依据。矿土还常受氮磷有效性限制；矿土硝化速率可低至未扰动土壤的7%，酸性条件下磷固定增加。采矿和复垦还改变容重、有机碳和养分库。由于母质组合和复垦实践在站点间及站点内变化，遗留矿土肥力限制预计具有站点特异性，不能从单一地点外推。现有比较多集中于单一复垦区内土地利用类型或复垦年龄；跨站点比较易与复垦技术、气候、地貌、表土材料和管理混淆。描述性阈值比较揭示位置特异土壤限制，主成分和聚类分析用于识别驱动肥力变异的变量，但多变量分类对数据预处理、极端值和变量选择敏感，结果只有在重采样和替代设定下可重复才有意义。研究人员认为，尚无比较–多变量框架将站点水平比较与多变量结构稳定性评估耦合应用于阿巴拉契亚多个独立遗留煤矿地。本研究是评估中阿巴拉契亚遗留矿地用于生物质开发的更广项目中的土壤表征部分；土壤化学是任何植被恢复策略的必要初筛，但不足以确定作物适宜性或生产力。研究人员分析西弗吉尼亚州和宾夕法尼亚州7个遗留煤矿站点土壤样品，以回答：站点间土壤肥力指标是否不同；哪些土壤性质解释样品间最大变异；样品能否分组成不同且可重复的土壤化学类别。研究人员明确评估多变量结果对转换、变量选择、极端值和重采样的稳健性，评估站点内空间依赖，并将结果解释为可指导但不替代田间评价的土壤化学基线。

二、主要关键技术方法（不超过250字）
研究人员采集西弗吉尼亚州和宾夕法尼亚州7个遗留煤矿地0–20 cm表层土壤，最终获得135个样品，站点样本数为Marion 20、Webster 23、Glen Rogers Rd 19、Turkey Creek East 13、Monongalia 20、Raleigh 20、Fayette 20；样品来自六个县，属于Appalachian Plateaus生理区。测定18种土壤测试性质。统计采用Kruskal–Wallis检验（Holm校正）、Dunn–Bonferroni两两比较、ICC(1)（intraclass correlation coefficient，组内相关系数）、Moran’s I和Mantel检验；对17种性质进行主成分分析（PCA），先winsorize和z标准化，并行分析确定成分数，varimax旋转，Tucker一致性系数比较替代设定，bootstrap评估变异。聚类采用Ward最小方差连接，k-means调整Rand指数（ARI）比较，Jaccard相似度、预测强度、留一站分析和替代设定评估稳定性。

三、研究结果
3.1 描述统计：研究人员总结135个样品18种性质。pH和缓冲pH相对变异较低；Ni、Ca、Zn变异系数高于175%，磷饱和比（phosphorus saturation ratio, PSR）和P高于125%。Ca绝对范围最宽，从强酸性Fayette样品的22 mg kg^?1（pH 4.6）到近中性Webster样品的26,160 mg kg^?1（pH 7.4），极端值经实验室数据核对，反映单个样品而非站点整体。Turkey Creek East pH最不酸且变异最大，若干样品近中性至碱性；Fayette有机质、PSR、P、K突出而pH最低；B模式跟踪实验室批次而非地理。多数性质站点内离散大，站点间中心趋势差异与大量重叠共存。

3.2 分布假设：站点内Shapiro–Wilk检验显示，除Al外，每种性质至少在1个站点偏离正态或呈退化分布；126个性质×站点组合中84个非正态或退化。Brown–Forsythe检验显示18种性质中12种站点间方差不齐。按预设规则无性质符合ANOVA，故全部用Kruskal–Wallis检验。合并样品统计也非正态、主要右偏，支持第4阶段使用Spearman相关。

3.3 站点间差异：所有18种性质在未校正和Holm校正后均站点间差异。非批次影响性质中，Na、有机质、Mn效应最大，其次P、Cu、PSR、Zn、Al；pH和缓冲pH也大效应，电导率（electrical conductivity, EC）、Ca、Fe、Ni中等。B效应很大但与实验室批次重合，不解释为站点效应。Dunn–Bonferroni识别94个显著站点对，其中83个不涉及B；Fayette和Marion分别参与36和34个。秩ICC从0.08（Ni）到0.62（Na），表明同站点样品中至强相似，外推到更广站点总体时有效独立观测数远小于样品数。站点采样范围1.8 ha（Raleigh）至7.7 ha（Webster, Bolair），最近邻间距中位数18–41 m。125个站点内Moran’s I检验中30个显示正空间自相关（p<0.05），约偶然期望5倍，11个经错误发现率控制后仍显著，集中于Marion和Monongalia的酸度相关性质；Mantel检验在Marion（r=0.34，p=0.002）和较弱在Glen Rogers Rd（r=0.17，p=0.042）显示样品越近整体土壤化学越相似，其他5站点无此关系。

3.4 相关与主成分分析：13对性质绝对Spearman相关系数≥0.5。PSR与P关联最强部分因定义所致；pH与缓冲pH强相关，pH与Ca、Mg正相关，与石灰需求（lime requirement, LR）负相关。Zn与Na、Cu关联。有机质仅弱关联酸度和盐基阳离子性质，与Mn、Al负相关，与Cu、Zn正相关，说明其变化较大程度独立于酸度驱动化学。总体KMO=0.645，Bartlett检验显著，除Al的MSA为0.35外均≥0.51。五个成分特征值>1，解释69.9%方差；平行分析支持前4个，第5个仅勉强超过Kaiser标准。旋转后：PC1为酸度/盐基阳离子成分（pH、Ca、缓冲pH、Mg、EC正，LR负）；PC2为痕量金属成分（Zn、Ni、Cu、Na）；PC3为磷有效性成分（P、PSR、K）；PC4为Al–Fe成分（Al、Fe正，缓冲pH负）；PC5为Mn–有机质成分（Mn正，有机质负）。前3个成分高度可重复：PC1一致性≥0.96、PC3≥0.97、PC2≥0.88；bootstrap第5百分位分别为0.88、0.87、0.92。PC4和PC5较不稳定，解释暂定。

3.5 聚类分析与稳定性：Ward聚类中最终合并融合距离增幅最大，平均轮廓宽度在k=2最高为0.27。两簇为酸性Cluster 1（n=106）和近中性Cluster 2（n=29）。k-means划分96.3%一致（ARI=0.84）。稳定性显示Cluster 1稳定（平均Jaccard=0.90，从未溶解），Cluster 2仅中等稳定（Jaccard=0.67），24%重采样中溶解。k=2预测强度0.72，低于0.8阈值；k=3和k=4分别为0.49和0.35。留一站分析中90–100%保留样品归入与原分析相同簇，省略任一站仅中度改变其余样品划分（ARI=0.66–1.00），近中性组不依赖单一站点。包含B、log转换、去除最极端样品和k-means时与主划分一致性高；未转换、去除Al或EC时中等；加入7个bootstrap生成记录中度一致并扩大近中性簇。秩转换、保留成分得分、去除LR和PSR未重现酸性/近中性划分，而分离少量极端样品或产生无关划分。PC1–PC2平面上，两簇位于连续酸度梯度两端，痕量金属梯度穿过两簇。Holm校正后簇在11种性质上不同：Cluster 2有更高pH、缓冲pH、EC、Ca、Mg、Na、Zn、Cu、Mn，更低LR和P；最大效应为LR（r=0.72）、pH（0.66）、缓冲pH（0.64）、Mg（0.51）、Ca（0.50）。pH≥6.0规则与聚类94.8%一致，LR=0规则94.1%一致。簇成员与站点相关但非一一对应：Fayette无Cluster 2，Turkey Creek East 8/13、Raleigh 7/20贡献较多，其他站点均含两类型。

四、讨论总结与结论翻译
讨论部分围绕三个研究问题总结。对于RQ1，所有18种性质在多重检验校正后站点间差异，但B差异与实验室批次重合，LR差异可能部分来自程序性因素；中至高ICC表明同站点样品非独立重复，且仅7站点，结果说明这些站点不同，但不能说明原因或推广到所有遗留煤矿地。酸度虽效应量小于Na、Mn或有机质，但决定石灰需求并定义主要多变量梯度PC1，是最具管理相关性的站点差异：pH跨度3.7–8.0，四分之三样品获得非零石灰推荐，Turkey Creek East最不酸。原生未采矿阿巴拉契亚森林土壤通常pH约4.5–6.5，一些中阿巴拉契亚森林表土平均pH 3.7–3.9；本范围同时覆盖原生森林酸性端和新鲜碳酸盐母质碱性端。持久站点间pH差异与废石残余产酸或中和潜力有关，但本研究未测酸碱核算。对于RQ2，PCA识别三个可重复梯度：PC1酸度/盐基阳离子最稳健；PC2痕量金属与酸度大体独立，可能反映覆盖层和废石痕量金属变化，但未测岩性、废石放置和金属形态，因此是未来检验假设；PC3磷有效性，所有PSR低于10%，描述供磷差异而非环境磷流失风险；Fayette有机质、P、PSR最高，但现有记录不能确定原因。PC4和PC5较不可重复，解释谨慎。有机质弱关联酸度和养分性质，负关联Mn，可能反映有机络合、植被凋落物输入或地质煤碳计入有机质，研究未区分，故高有机质不应直接读作生物土壤肥力。RQ2答案是样品水平变异多维，酸度/盐基状态、痕量金属富集和磷有效性 largely independent，酸度单独不能完整代理整体土壤状况。对于RQ3，Ward聚类产生酸性需石灰组和近中性组，但分离弱、仅酸性组重采样稳定、近中性组约四分之一bootstrap中溶解，且划分对极端值和派生酸度变量处理敏感，几乎完全由pH 6.0阈值重现。因此样品可分成酸性多数和较小、较不一致的近中性组，但更适合理解为连续酸度梯度两端而非离散土壤类别。该分类为探索性，经内部重采样和站点排除检查，但未用独立站点验证；简单pH阈值可重现，故聚类超出常规土壤测试的操作价值有限，但主要贡献是指出痕量金属富集是酸度类别未捕获的独立维度，需独立评估。簇成员与站点相关但非一一对应，站点内异质性常见，混合站点更适合更密采样，站点身份本身不足以作为全站点统一土壤处方基础。对于复垦规划和未来生物能源作物试验，本研究未包含作物建立、生物量、植物组织或作物耐受性测定，不能确定任何生物能源作物的适宜性或生产力；意义限于任何植被恢复策略必须适应的土壤化学条件及田间试验设计。三点启示：酸度是最广泛限制，田间试验应包含石灰处理或跨观测pH范围测试物种；痕量金属状态独立于酸度，应直接测量而非从pH或站点身份推断，数据集Zn 79 mg kg^?1和Cu 32 mg kg^?1低于美国EPA植物生态土壤筛选水平160和70 mg kg^?1，但5个样品（3个来自Turkey Creek East）超过Mn植物筛选水平220 mg kg^?1，筛选水平为总浓度保守基准而非Mehlich-3提取浓度，只标识试验设计中应纳入Mn耐受的地点而非已证明植物毒性；酸性与近中性在多个站点共存，试验小区应基于样品水平土壤测试定位。讨论中列举的物种如switchgrass（Panicum virgatum）、giant miscanthus（Miscanthus × giganteus）、native willows（Salix discolor和S. eriocephala）等是候选测试对象而非推荐。局限性包括：仅7站点、可及性采样非概率设计、站点内空间自相关减少有效样本量、3个样品无记录位置、Webster Ponds 3个样品太少、Turkey Creek East仅13样品；公共记录只确立监管状态而无复垦日期、方法、材料、土地利用、石灰和改良历史，站点差异联合反映这些因素；实验室批次与站点部分混淆，B有批次效应，Na和LR可能受影响；KMO=0.645，Al低MSA，PC4和PC5不稳健，仅前三个成分可靠；两簇弱分离且未独立验证；有机质可能含地质煤碳，未测阳离子交换量、物理性质、总金属、金属形态、生物指标，未来需测容重、碎石、持水能力、根区水分等，并进行多指标时间序列评估。结论翻译：本研究使用一致工作流程比较阿巴拉契亚7个遗留煤矿地土壤化学，并明确评估多变量结果稳健性。对于RQ1，所有18种土壤测试性质在多重检验校正后站点间差异，尽管硼差异反映实验室批次，且站点内样品中至强相关，所以这些差异描述这些站点而非一般遗留煤矿地。对于RQ2，研究人员将样品水平变异组织为三个可重复、 largely independent 的梯度：酸度/盐基阳离子状态、痕量金属富集和磷有效性。对于RQ3，样品可分为酸性需石灰多数和较小近中性组，但分离弱；仅酸性组稳定，划分几乎完全由pH阈值6.0重现；因此分类是对连续酸度梯度的探索性描述，而非经过验证的土壤类别集。对于复垦规划，这些结果表明酸度和痕量金属状态应在样品而非站点水平分别评估，然后再植被恢复。本研究未评价作物表现；确定特定生物能源或复垦物种是否适合这些土壤，将需要重复田间试验，在本文记录的酸度和痕量金属条件下测量建立、生物量产量和植物组织化学。