利用地球化学探索考古源区溯源的算术方法——对Pearce等人评论的回复

《Journal of Archaeological Science: Reports》：Arithmetic methods for exploring archaeological source provenance using geochemistry – A reply to Pearce et al.

　　

研究人员 defending 将算术方法作为与其他数据类型一并用于源区溯源调查的探索性工具；defending 在硅结砾岩(silcrete)溯源中使用微量元素/Zr比值；承认Pearce等人(2026)提出的三项方法学改进；论证了使用双变量地球化学图解的局限性，转而支持多元分析。研究人员感谢Pearce等人(2026)对其与Nash于2026年提出的评估石制品与潜在源区露头间地球化学相似性的新算术方法所作的批判性评价。该评论超出对2026年论文的评议，还涉及研究人员此前关于巨石阵(zarsen)巨石及石片的研究，其框架与论证路径存在问题，故需更广范围的回应。研究人员指出评论中存在对其工作的误读（如将其方法歪曲为可孤立替代岩石学等证据的“万能钥匙”，即稻草人论证），并澄清其动机并非为既往比值使用“辩护”，而是源于巨石阵数据集对比解释的自然延伸。研究人员接纳三项有用建议：低浓度微量元素因舍入误差对算术计算影响放大，后续使用应引入额外误差容限；认证标准物质(CRM)宜选用近硅质材料以更稳健估计精度；可按分析不确定度对Ei/Zr比值进行逆权重处理。研究人员论证硅结砾岩中硅化稀释效应不影响非Si元素相对比值，故Ei/Zr优于浓度数据；指出双变量图解在Ustica黑曜岩及巨石阵祭坛石(Altar Stone)案例中具选择性偏差，而几何均值ΔEi/Zr(%)可同时整合多因素，支持作为探索性补充。研究人员反驳Pearce等人自身亦曾单凭地球化学定源、使用PCA等变换并聚合U-Pb锆石年代学数据，与其倡导的“总岩石学(total petrography)”标准不自洽。最终研究人员重申算术方法仅提供地球化学相似性度量，须与岩石学、矿物学及考古背景联用。

研究背景方面，考古源区溯源长期依赖岩石学、矿物学与地球化学多证据交汇，其中Ixer(1997)提出的“总岩石学(total petrography)”主张矿物成分、结构、化学与（如有）年代学全对齐。然而实际田野中，受遗产保护、采样破坏性与成本限制，研究人员常仅获全岩地球化学数据，尤其对高硅质硅结砾岩(silcrete，SiO₂常>85 wt.%)而言，非Si微量元素浓度极低且受硅化稀释与重矿物纹层干扰严重。Ciborowski和Nash于2026年提出一种基于几何均值ΔE_i/Zr(%)的算术方法，用以量化石制品与露头间地球化学相似性，引发Pearce等人(2026)评论，指其敏感于舍入、忽略兼容元素、滥用稀土元素(REE)、脱离岩石学等。本回复论文发表于《Journal of Archaeological Science: Reports》，旨在纠正误读、澄清动机并吸纳合理改进。

关键技术方法方面，研究人员采用文献既发表数据集（含巨石阵zarsen露头与Stone 58 Phillips’ Core、Preseli Hills七处辉绿岩露头与Stone 62岩芯、Ustica黑曜岩制品与Sardinia/Lipari/Palmarola源区、加拿大Marathon与Biscotasing镁铁质岩墙群ICP-MS数据，以及Bevins等人(2023)pXRF巨石阵Altar Stone与碎屑样2010K.240数据）作再分析；核心算法为将各非Si微量元素浓度E_i除以Zr得E_i/Zr比值，取制品与露头均值比之几何均值百分差ΔE_i/Zr(%)；辅以双变量图解重绘、U-Pb锆石年龄直方图分解、Spearman相关计算，以检验评论方“总岩石学”实践内部一致性。

研究结果按原小节归纳如下。

误读工作部分，研究人员指出Pearce等人将其方法歪曲为孤立替代岩石学的“魔法子弹”，忽视原文明确写有“非万灵药、须结合岩石学等指示”的语句，构成稻草人论证；并说明“总岩石学”在数百露头采样中不现实，算术法定位为初步探索工具。

误读动机部分，研究人员否认因与Hancock等人(2024)争论而“发明方法自辩”，重申方法源于Stonehenge数据集对比解释的自然需求，且从未与Pearce等人沟通动机。

方法改进三点，研究人员接纳：低浓度舍入放大问题，建议后续对E_i/Zr引入浓度区间误差容限；认证标准物质(CRM)应近硅质（如OREAS 21h/293等石英砂/砾岩），原Nash等人(2020)用ALS Minerals常规CRM不理想但不涉算术理论本身；可按分析不确定度逆权重各E_i/Zr，原论文为聚焦实用性未植入该系数但已警示只用高质元素。

硅结砾岩小节，研究人员论证硅化稀释与重矿物纹层改变浓度但不改非Si相对比值，故E_i/Zr优于浓度；指出Pearce等人“样品增大则变异减小”不成立，因硅化度在手标本至巨砾间仍剧变。

REE选择小节，研究人员称使用者可自定元素集；锆石非唯一载REE相（磷灰石、钛铁矿、云母、粘土、独居居石等均可），Totterdown Wood与West Woods源区正是靠全REE/Zr区分；同一E_i/Zr标尺比较无偏向。

Stone 62与Preseli小节，研究人员示明算术法算得Stone 62岩芯与Carn Goedog化学最似，但从未断言其为源；指出Bevins等人(2014)曾单凭地球化学将非斑状辉绿岩制品归Carn Goedog，评论方自悖；并指评论方自身用PCA等变换，却批算术“失真”属双重标准。

兼容元素小节，研究人员引原文Worked Example 5已演示先用不兼容元素/E_i/Zr筛岩套、再用兼容元素(Ni,Cr,Sc)/Zr定层位，Preseli例未重复此步但文中明述，非“忽略兼容元素”。

Ustica、双变量图解、总岩石学与Altar Stone小节，研究人员重绘Ustica的Th-Yb、Nb-Rb等双变量图，显示Palmarola源制品落Lipari区、Lipari源制品落Sardinia区，证明双变量选择性呈现会误导；再以Bevins等人(2023)干样pXRF数据绘Altar Stone与2010K.240的P-Zn-Mn-K-Ca-Fe-Ba-Sr图，示两样在Mn/Zn/P上系统性分离且2010K.240不落Altar Stone强相关阵列(Spearman ρ=0.85)，评其“总岩石学”未解释这些元素；进一步拆Clarke等人(2024)聚合U-Pb锆石年龄，示三个碎屑样2010K.240/MS3/FN593彼此年代分布异，五个Orcadian Basin砂岩样亦互异，单体比对无相似，质疑苏格兰Old Red Sandstone源结论。

对Pearce等人结论逐条回应中，研究人员重申：灵敏度问题不否定方法发展；原文已言无硬置信限；权重改进可接受；未只用水兼容元素；REE行为被简化；双变量图有价值但非客观；Stone 62仅示“最似”非定源；任何统计变换均“去原始性”，评论方自用PCA；动机诽谤“隐藏变异”无据，所有例题数据公开。

讨论与结论翻译部分，研究人员总结：感谢Pearce等人压力测试并指出三改进点；但其论证建基于选择性引文与误读，且自身“总岩石学”实践（单用化学定源、PCA变换、聚合异层位锆石）前后不一；双变量图若选择性展示亦具缺陷。研究人员建议读者细读Ciborowski和Nash(2026)原文并带其明示 caveat 使用方法。结论重申算术法仅给地球化学相似性度量，须与岩石学、矿物学、考古情境联用，非孤立万灵药；在硅结砾岩等高纯材料中E_i/Zr可抗硅化稀释；方法可扩展权重与误差容限；未来须以多方法交叉验证推进巨石阵等 monument 理解。