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

Sintesi Tecnica: Una Chiusura per l'Attrito Effettivo dello Stato Strutturale per le Interfacce di Subduzione delle Placche Oceaniche

Problematica
La resistenza d'attrito delle interfacce di subduzione delle placche oceaniche non può essere misurata direttamente su scala reale poiché l'interfaccia è sepolta, idratata, termicamente evolutiva e geometricamente complessa. L'interfaccia comprende un mix eterogeneo di asperità bloccate, patch di creep, fluidi, sedimenti e minerali alterati. Di conseguenza, la resistenza apparente non è una semplice coppia di attrito di laboratorio, ma un fenomeno di "attrito effettivo" su scala di campo. La sfida osservativa immediata consiste nello sviluppare una formula utile che operi come una chiusura efficace su scala di campo, proiettando stati strutturali inaccessibili su osservabili quali il flusso di calore, i coefficienti di attrito effettivo, l'accoppiamento geodetico e il rilascio del momento sismico. Questo studio testa se una "architettura dell'attrito a stato strutturale", precedentemente applicata ad altri problemi di attrito su larga scala inaccessibili, possa essere significativamente proiettata sulle osservazioni di subduzione.

Metodologia
Lo studio impiega un audit osservativo a due rami per testare la chiusura dello stato strutturale rispetto ai dati reali:

1. Audit del Flusso di Calore e del Modello Termico:

   • È stato compilato un dataset di vincoli codificati dalla letteratura, che riassume i coefficienti di attrito effettivo ($\mu'$) derivati da modelli di flusso di calore e termici in varie zone di subduzione.
   • Questi vincoli indicano tipicamente bassi valori di attrito effettivo (compresi tra 0,03 e 0,07), significativamente inferiori al valore di attrito di Byerlee a secco ($\mu \approx 0,6$).
   • Una chiusura dell'attrito effettivo a stato strutturale ridotta è stata adattata a questi vincoli intervallari e confrontata con una linea di base a coefficiente di attrito costante.

2. Audit del Momento Sismico USGS:

   • I dati del catalogo sismico (M $\ge$ 7) dal 1900 ad oggi sono stati scaricati per 12 regioni di subduzione rappresentative tramite l'API FDSN dell'USGS, ottenendo 715 eventi.
   • Il momento sismico totale ($M_0$) è stato calcolato per ogni regione come una proiezione grossolana del rilascio di momento stick-slip a lungo termine.
   • Sono stati testati vari modelli di proiezione rispetto a questi dati di momento, tra cui:
     • La sola velocità di convergenza.
     • Lo stato termico da solo.
     • Un semplice "indice di carico strutturale".
     • Un modello corretto che combina velocità, stato termico e inclinazione della lastra (slab dip).
   • Le prestazioni sono state valutate utilizzando la validazione incrociata leave-one-out (LOO) per valutare la trasferibilità regionale.

Il Modello di Chiusura dello Stato Strutturale
La chiusura proposta definisce l'attrito effettivo ($\mu_{eff}$) come funzione della velocità di convergenza delle placche ($V$) e di uno stato di indebolimento termico-fluido ridotto ($\Theta$). Il modello è strutturato su un'architettura a quattro livelli:

1. Insieme Fisico Fattibile: Imposto dalla meccanica di subduzione (es. $0 \le \mu_{eff} \le \mu_{Byerlee}$).
2. Coordinata Strutturale: Una variabile composita ($\Phi$) che integra bloccaggio, rugosità, contenuto di fluidi, minerali alterati, creep e geometria.
3. Risposta Strutturale: Una forma di perdita di produzione che rappresenta la risposta del sistema alla coordinata strutturale.
4. Proiezione: Mappatura dello stato su osservabili misurabili (flusso di calore, accoppiamento, scivolamento lento, momento sismico).

Risultati Chiave

• Audit del Flusso di Calore: La chiusura dello stato strutturale ricade con successo all'interno di tutti gli intervalli di attrito effettivo compilati (0,03–0,07). Produce una perdita di intervallo leggermente inferiore (0,00003544) rispetto alla migliore linea di base a coefficiente di attrito costante (0,00004225). Entrambi i modelli rifiutano correttamente il valore di attrito di laboratorio a secco di 0,6, confermando che le osservabili di subduzione rappresentano un coefficiente effettivo debole, modificato da fluidi e temperatura.
• Audit del Momento Sismico:
  • Fallimento della Velocità-Solo: La sola velocità di convergenza è un cattivo predittore del momento sismico totale (LOO RMSE = 0,8293). Ad esempio, la zona Tonga-Kermadec ha un'alta velocità di convergenza ma un momento totale inferiore rispetto a zone più lente come il Perù-Cile o il Giappone.
  • Fallimento del Semplice Indice Strutturale: Un semplice indice di carico strutturale ha performato scarsamente (LOO RMSE = 0,7210), dimostrando che la formula non può essere trattata come una legge di attrito universale basata solo sulla velocità.
  • Successo Termico-Geometrico: I modelli che incorporano lo stato termico e la geometria della lastra hanno performato significativamente meglio. La proiezione "velocità-termica-inclinazione" ha ottenuto il LOO RMSE più basso (0,6005) tra i modelli non basati su catalogo, superando le linee di base basate solo su velocità o solo su temperatura. Il modello "structural_interaction" ha prodotto un errore ancora maggiore (LOO RMSE = 0,8876).

Significato e Rivendicazioni
Il documento mantiene un ambito modesto, dichiarando esplicitamente di non pretendere di sostituire i modelli termomeccanici di subduzione, le teorie di attrito rate-and-state o le simulazioni del ciclo sismico. Invece, offre una validazione a livello pilota di un approccio architettonico specifico.

Le conclusioni più solidamente supportate sono:

1. Compatibilità: L'attrito effettivo di subduzione è compatibile con una chiusura a stato strutturale. Il modello copre con successo i bassi intervalli di attrito osservati negli studi sul flusso di calore.
2. Necessità di Proiezione: La proiezione osservabile di questa chiusura deve includere lo stato termico-fluido e la geometria della lastra. Una formula basata esclusivamente sulla velocità di convergenza o su un semplice indice strutturale è insufficiente per spiegare il rilascio del momento sismico regionale.
3. Architettura sopra il Coefficiente Universale: La ricerca suggerisce che l'aspetto "universale" dell'attrito di subduzione non è un coefficiente numerico fisso, ma un'architettura: stato di subduzione fattibile duro $\rightarrow$ coordinata strutturale $\rightarrow$ risposta di perdita di produzione $\rightarrow$ proiezione termica-fluida-geometrica $\rightarrow$ osservabili.

In sintesi, lo studio fornisce prove del fatto che una chiusura dell'attrito a stato strutturale è un quadro viabile per le interfacce di subduzione, a condizione che sia corretta dai termini di proiezione termica-fluida e geometrica. La universalità esterna rimane non dimostrata, e l'evidenza attuale supporta una chiusura a livello pilota piuttosto che una legge universale definitiva.