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

Resumo Técnico: Um Fechamento de Atrito Efetivo de Estado-Estrutural para Interfaces de Subducção de Placas Oceânicas

Definição do Problema
A resistência friccional das interfaces de subducção de placas oceânicas não pode ser medida diretamente em escala total devido ao fato de a interface estar enterrada, hidratada, termicamente evoluindo e geometricamente complexa. A interface compreende uma mistura heterogênea de asperezas travadas, manchas de escorregamento (creep), fluidos, sedimentos e minerais alterados. Consequentemente, a resistência aparente não é um simples par de fricção de laboratório, mas um fenômeno de campo "efetivo". O desafio observacional imediato é desenvolver uma fórmula útil que opere como um fechamento de estado-estrutural efetivo em escala de campo, projetando estados estruturais inacessíveis sobre observáveis como fluxo de calor, coeficientes de atrito efetivos, acoplamento geodésico e liberação de momento sísmico. Este estudo testa se uma "arquitetura de atrito de estado-estrutural", previamente aplicada a outros problemas de fricção de grande escala inacessíveis, pode ser significativamente projetada sobre observações de subducção.

Metodologia
O estudo emprega uma auditoria observacional de duas vertentes para testar o fechamento de estado-estrutural contra dados do mundo real:

1. Auditoria de Fluxo de Calor e Modelo Térmico:

   • Um conjunto de dados de restrições codificadas da literatura foi compilado, resumindo coeficientes de atrito efetivos ($\mu'$) derivados de modelos de fluxo de calor e térmicos em várias zonas de subducção.
   • Essas restrições tipicamente indicam valores de atrito efetivo baixos (variando de 0,03 a 0,07), significativamente menores que o atrito de Byerlee seco de laboratório ($\mu \approx 0,6$).
   • Um fechamento de atrito efetivo de estado-estrutural reduzido foi ajustado a esses intervalos de restrição e comparado contra uma linha de base de atrito efetivo constante.

2. Auditoria de Momento Sísmico do USGS:

   • Dados de catálogos de terremotos (M $\ge$ 7) de 1900 até o presente foram baixados para 12 regiões de subducção representativas via API FDSN do USGS, resultando em 715 eventos.
   • O momento sísmico total ($M_0$) foi calculado para cada região como uma projeção grosseira da liberação de momento de stick-slip a longo prazo.
   • Vários modelos de projeção foram testados contra estes dados de momento, incluindo:
     • Velocidade de convergência isoladamente.
     • Estado térmico isoladamente.
     • Um simples "índice de carregamento estrutural".
     • Um modelo corrigido combinando velocidade, estado térmico e mergulho (dip) da placa.
   • O desempenho foi avaliado usando validação cruzada de deixar um fora (leave-one-out, LOO) para avaliar a transferibilidade regional.

O Modelo de Fechamento de Estado-Estrutural
O fechamento proposto define o atrito efetivo ($\mu_{eff}$) como uma função da velocidade de convergência das placas ($V$) e de um estado de enfraquecimento térmico-fluido ($\Theta$). O modelo é estruturado em uma arquitetura de quatro níveis:

1. Conjunto Físico Viável: Imposto pela mecânica de subducção (ex: $0 \le \mu_{eff} \le \mu_{Byerlee}$).
2. Coordenada Estrutural: Uma variável composta ($\Phi$) integrando travamento, rugosidade, conteúdo de fluido, minerais alterados, escorregamento (creep) e geometria.
3. Resposta Estrutural: Uma forma de perda de produção representando a resposta do sistema à coordenada estrutural.
4. Projeção: Mapeamento do estado para observáveis mensuráveis (fluxo de calor, acoplamento, escorregamento lento, momento sísmico).

Resultos Principais

• Auditoria de Fluxo de Calor: O fechamento de estado-estrutural situa-se com sucesso dentro de todos os intervalos de atrito efetivo compilados (0,03–0,07). Ele produz uma perda de intervalo ligeiramente menor (0,00003544) em comparação com a melhor linha de base de atrito efetivo constante (0,00004225). Ambos os modelos rejeitam corretamente o valor de atrito seco de laboratório de 0,6, confirmando que os observáveis de subducção representam um coeficiente efetivo fraco, modificado por fluidos e temperatura.
• Auditoria de Momento Sísmico:
  • Falha da Velocidade Isolada: A velocidade de convergência isoladamente é um pobre preditor do momento sísmico total (LOO RMSE = 0,8293). Por exemplo, a zona Tonga-Kermadec possui alta velocidade de convergência, mas menor momento total do que zonas mais lentas como Peru-Chile ou Japão.
  • Falha do Índice Estrutural Simples: Um simples índice de carregamento estrutural desempenhou mal (LOO RMSE = 0,7210), demonstrando que a fórmula não pode ser tratada como uma lei de atrito universal baseada apenas na velocidade.
  • Sucesso Térmico-Geométrico: Modelos incorporando o estado térmico e a geometria da placa performaram significativamente melhor. A projeção "velocidade-térmica-mergulho" alcançou o menor LOO RMSE (0,6005) entre os modelos de não-catálogo, superando as linhas de base de velocidade-apenas e térmica-apenas. O modelo "interação_estrutural" apresentou um erro ainda maior (LOO RMSE = 0,8876).

Significância e Alegações
O artigo mantém um escopo modesto, declarando explicitamente que não pretende substituir modelos termomecânicos de subducção, teorias de atrito de taxa-e-estado ou simulações de ciclo de terremotos. Em vez disso, oferece uma validação de nível piloto de uma abordagem arquitetônica específica.

As conclusões mais fortemente sustentadas são:

1. Compatibilidade: O atrito efetivo de subducção é compatível com um fechamento de estado-estrutural. O modelo cobre com sucesso os baixos intervalos de atrito observados em estudos de fluxo de calor.
2. Necessidade de Projeção: A projeção observável deste fechamento deve incluir o estado térmico-fluido e a geometria da placa. Uma fórmula baseada apenas na velocidade de convergência ou em um índice estrutural simples é insuficiente para explicar a liberação regional de momento sísmico.
3. Arquitetura sobre Coeficiente Universal: A pesquisa sugere que o aspecto "universal" do atrito de subducção não é um coeficiente numérico fixo, mas uma arquitetura: estado de subducção fisicamente viável $\rightarrow$ coordenada estrutural $\rightarrow$ resposta de perda de produção $\rightarrow$ projeção térmica-fluida-geométrica $\rightarrow$ observáveis.

Em resumo, o estudo fornece evidências de que um fechamento de atrito efetivo de estado-estrutural é um arcabouço viável para interfaces de subducção, desde que seja corrigido por termos de projeção térmica-fluida e geométrica. A universalidade externa permanece não provada, e a evidência atual apoia um fechamento de nível piloto em vez de uma lei universal definitiva.