A Structural-State Effective-Friction Closure for Oceanic-Plate Subduction Interfaces

技术摘要：一种面向大洋板块俯冲界面的结构状态有效摩擦闭合模型

问题陈述
由于大洋板块俯冲界面处于埋藏、水化、热演化且几何形状复杂的环境中，其摩擦阻力无法通过直接测量获得。该界面由锁定的粗糙体（asperities）、蠕变斑块、流体、沉积物以及蚀变矿物组成的异质混合物构成。因此，表观阻力并非简单的实验室摩擦对，而是一种“有效”的场尺度现象。直接的观测挑战在于，如何开发一种有效的场尺度闭合公式，将不可达的结构状态投影到热流、有效摩擦系数、地壳形变耦合及地震矩释放等可观测变量上。本研究测试了是否可以将此前应用于其他不可达大规模摩擦问题的“结构状态摩擦架构”有效地投影到俯冲观测数据中。

研究方法
本研究采用两阶段观测审计，以测试结构状态闭合模型与现实世界数据的匹配度：

1. 热流与热模型审计：

   • 编译了一个文献编码约束数据集，汇总了不同俯冲带中基于热流和热模型推导出的有效摩擦系数 ($\mu'$)。
   • 这些约束通常表明有效摩擦值较低（范围在 0.03 至 0.07 之间），显著低于干燥实验室环境下的 Byerlee 摩擦系数 ($\mu \approx 0.6$)。
   • 将一个简化的结构状态有效摩擦闭合模型拟合到这些区间约束中，并将其与恒定有效摩擦基准进行对比。

2. USGS 地震矩审计：

   • 通过 USGS FDSN API 下载了 1900 年至今（M $\ge$ 7）的地震目录数据，涵盖 12 个代表性俯冲区域，共计 715 个事件。
   • 计算了每个区域的总地震矩 ($M_0$)，作为长期粘滑（stick-slip）矩释放的粗略投影。
   • 针对该矩数据测试了多种投影模型，包括：
     • 仅收敛速度。
     • 仅热状态。
     • 简单的“结构加载指数”。
     • 结合了速度、热状态和板块倾角（slab dip）的修正模型。
   • 使用留一法（leave-one-out, LOO）交叉验证来评估区域迁移性，并对性能进行评估。

结构状态闭合模型
所提出的闭合模型将有效摩擦 ($\mu_{eff}$) 定义为板块收敛速度 ($V$) 和降低后的热-流体弱化状态 ($\Theta$) 的函数。该模型构建了一个四层架构：

1. 物理可行集： 由俯冲力学机制规定（例如，$0 \le \mu_{eff} \le \mu_{Byerlee}$）。
2. 结构坐标： 一个综合变量 ($\Phi$)，整合了锁定、粗糙度、流体含量、蚀变矿物、蠕变及几何形状。
3. 结构响应： 代表系统对结构坐标响应的“生产损失”（production-loss）形式。
4. 投影： 将状态映射到可观测量（热流、耦合、慢滑移、地震矩）。

主要结果

• 热流审计： 结构状态闭合模型成功覆盖了所有编译的有效摩擦区间（0.03–0.07）。其产生的区间损失（0.00003544）略小于最佳恒定有效摩擦基准（0.00004225）。两种模型均正确拒绝了 0.6 的干燥实验室摩擦值，证实了俯冲观测值代表的是一种受流体和温度调制的弱有效系数。
• 地震矩审计：
  • 仅速度模型的失效： 单纯的收敛速度是总地震矩的较差预测因子（LOO RMSE = 0.8293）。例如，汤加-克马德克带具有较高的收敛速度，但其总矩却低于秘鲁-智利或日本等速度较慢的区域。
  • 简单结构指数的失效： 简单的结构加载指数表现不佳（LOO RMSE = 0.7210），这表明该公式不能被视为一种通用的仅基于速度的摩擦定律。
  • 热状态-几何形状的成功： 结合了热状态和板块几何形状的模型表现显著更好。“速度-热状态-倾角”投影模型在非目录模型中实现了最低的 LOO RMSE（0.6005），优于仅速度和仅热状态的基准模型。“结构交互”（structural_interaction）模型产生的误差则更高（LOO RMSE = 0.8876）。

意义与主张
本文保持了审慎的研究范围，明确声明其并不旨在取代热力学俯冲模型、速率-状态摩擦理论或地震循环模拟。相反，它提供了一种特定架构方法的初步验证。

其最强有力的结论如下：

1. 兼容性： 俯冲有效摩擦与结构状态闭合模型是兼容的。该模型成功覆盖了热流研究中观察到的低摩擦区间。
2. 投影的必要性： 该闭合模型的观测投影 必须 包含热-流体状态和板块几何形状。仅基于收敛速度或简单结构指数的公式不足以解释区域地震矩释放。
3. 架构优于通用系数： 研究表明，俯冲摩擦的“通用性”不在于一个固定的数值系数，而在于一种架构：硬性的俯冲物理可行状态 $\rightarrow$ 结构坐标 $\rightarrow$ 生产损失响应 $\rightarrow$ 热-流体-几何投影 $\rightarrow$ 观测值。

综上所述，本研究提供了证据，证明结构状态有效摩擦闭合模型是一个可行的框架，前提是必须通过热-流体和几何投影项进行修正。外部的普适性尚未得到证实，目前的证据支持将其视为一种初步层级的闭合模型，而非确定的通用定律。