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

Resumen Técnico: Un Cierre de Fricción Efectiva de Estado-Estructural para Interfaces de Subducción de Placas Oceánicas

Planteamiento del Problema
La resistencia friccional de las interfaces de subducción de placas oceánicas no puede medirse directamente a escala completa debido a que la interfaz se encuentra enterrada, hidratada, en evolución térmica y con una geometría compleja. La interfaz comprende una mezcla heterogénea de asperezas bloqueadas, parches de deslizamiento (creep), fluidos, sedimentos y minerales alterados. En consecuencia, la resistencia aparente no es un simple par de fricción de laboratorio, sino un fenómeno de campo "efectivo". El desafío observacional inmediato es desarrollar una fórmula útil que opere como un cierre de campo efectivo, proyectando estados estructurales inaccesibles sobre observables como el flujo de calor, los coeficientes de fricción efectiva, el acoplamiento geodésico y la liberación de momento sísmico. Este estudio pone a prueba si una "arquitectura de fricción de estado-estructural", aplicada previamente a otros problemas de fricción a gran escala e inaccesibles, puede proyectarse significativamente sobre las observaciones de subducción.

Metodología
El estudio emplea una auditoría observacional de dos vertientes para probar el cierre de estado-estructural contra datos del mundo real:

1. Auditoría de Flujo de Calor y Modelo Térmico:

   • Se compiló un conjunto de datos de restricciones codificadas de la literatura, resumiendo los coeficientes de fricción efectiva ($\mu'$) derivados de modelos de flujo de calor y térmicos en diversas zonas de subducción.
   • Estas restricciones indican típicamente valores de fricción efectiva bajos (que oscilan entre 0.03 y 0.07), significativamente menores que la fricción de Byerlee de laboratorio seco ($\mu \approx 0.6$).
   • Se ajustó un cierre de fricción efectiva de estado estructural reducido a estos intervalos de restricción y se comparó con una línea base de fricción efectiva constante.

2. Auditoría de Momento Sísmico del USGS:

   • Se descargaron datos de catálogos de terremotos (M $\ge$ 7) desde 1900 hasta el presente para 12 regiones de subducción representativas mediante la API FDSN del USGS, obteniendo 715 eventos.
   • Se calculó el momento sísmico total ($M_0$) para cada región como una proyección aproximada de la liberación de momento de deslizamiento-deslizamiento (stick-slip) a largo plazo.
   • Se probaron varios modelos de proyección contra estos datos de momento, incluyendo:
     • Velocidad de convergencia por sí sola.
     • Estado térmico por sí solo.
     • Un simple "índice de carga estructural".
     • Un modelo corregido que combina velocidad, estado térmico y buzamiento (dip) de la placa.
   • El rendimiento se evaluó mediante validación cruzada de dejar uno fuera (leave-one-out, LOO) para evaluar la transferibilidad regional.

El Modelo de Cierre de Estado-Estructural
El cierre propuesto define la fricción efectiva ($\mu_{eff}$) como una función de la velocidad de convergencia de la placa ($V$) y de un estado de debilitamiento térmico-fluido reducido ($\Theta$). El modelo está estructurado como una arquitectura de cuatro niveles:

1. Conjunto Físicamente Factible: Impuesto por la mecánica de subducción (ej., $0 \le \mu_{eff} \le \mu_{Byerlee}$).
2. Coordenada Estructural: Una variable compuesta ($\Phi$) que integra el bloqueo, la rugosidad, el contenido de fluidos, los minerales alterados, el deslizamiento (creep) y la geometría.
3. Respuesta Estructural: Una forma de pérdida de producción que representa la respuesta del sistema a la coordenada estructural.
4. Proyección: Mapeo del estado hacia observables medibles (flujo de calor, acoplamiento, deslizamiento lento, momento sísmico).

Resultos Clave

• Auditoría de Flujo de Calor: El cierre de estado estructural se sitúa con éxito dentro de todos los intervalos de fricción efectiva compilados (0.03–0.07). Produce una pérdida de intervalo ligeramente menor (0.00003544) en comparación con la mejor línea base de fricción efectiva constante (0.00004225). Ambos modelos rechazan correctamente el valor de fricción de laboratorio seco de 0.6, confirmando que los observables de subducción representan un coeficiente efectivo débil, modificado por fluidos y temperatura.
• Auditoría de Momento Sísmico:
  • Fallo de Velocidad Única: La velocidad de convergencia por sí sola es un mal predictor del momento sísmico total (LOO RMSE = 0.8293). Por ejemplo, la zona Tonga-Kermadec tiene una alta velocidad de convergencia pero un momento total menor que zonas más lentas como Perú-Chile o Japón.
  • Fallo del Índice Estructural Simple: Un índice de carga estructural simple funcionó deficientemente (LOO RMSE = 0.7210), demostrando que la fórmula no puede tratarse como una ley de fricción universal basada solo en la velocidad.
  • Éxito Térmico-Geométrico: Los modelos que incorporan el estado térmico y la geometría de la placa funcionaron significativamente mejor. La proyección "velocidad-térmica-buzamiento" logró el RMSE de LOO más bajo (0.6005) entre los modelos de no catálogo, superando a las líneas base de velocidad-solo y térmica-solo. El modelo de "interacción_estructural" arrojó un error aún mayor (LOO RMSE = 0.8876).

Significancia y Reivindicaciones
El artículo mantiene un alcance modesto, declarando explícitamente que no pretende reemplazar los modelos termomecánicos de subducción, las teorías de fricción de tasa y estado, ni las simulaciones del ciclo de terremotos. En su lugar, ofrece una validación a nivel piloto de un enfoque arquitectónico específico.

Las conclusiones más sólidas respaldadas son:

1. Compatibilidad: La fricción efectiva de subducción es compatible con un cierre de estado estructural. El modelo cubre con éxito los bajos intervalos de fricción observados en los estudios de flujo de calor.
2. Necesidad de Proyección: La proyección observable de este cierre debe incluir el estado térmico-fluido y la geometría de la placa. Una fórmula basada únicamente en la velocidad de convergencia o en un índice estructural simple es insuficiente para explicar la liberación de momento sísmico regional.
3. Arquitectura sobre Coeficiente Universal: La investigación sugiere que el aspecto "universal" de la fricción de subducción no es un coeficiente numérico fijo, sino una arquitectura: estado factible de subducción dura $\rightarrow$ coordenada estructural $\rightarrow$ respuesta de pérdida de producción $\rightarrow$ proyección térmica-fluido-geométrica $\rightarrow$ observables.

En resumen, el estudio proporciona evidencia de que un cierre de fricción efectiva de estado-estructural es un marco viable para las interfaces de subducción, siempre que se corrija mediante términos de proyección térmica-fluida y geométrica. La universalidad externa permanece sin probar, y la evidencia actual respalda un cierre a nivel piloto en lugar de una ley universal definitiva.