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

Technisches Resümee: Ein struktureller Zustands-Effektiv-Reibungs-Abschluss für ozeanische Platten-Subduktionsgrenzflächen

Problemstellung
Der Reibungswiderstand von ozeanischen Platten-Subduktionsgrenzflächen kann aufgrund der Tatsache, dass die Grenzfläche vergraben, hydriert, thermisch evolvierend und geometrisch komplex ist, nicht direkt im vollen Maßstab gemessen werden. Die Grenzfläche besteht aus einem heterogenen Gemisch aus arretierten Asperitäten, kriechenden Patches, Fluiden, Sedimenten und veränderten Mineralen. Folglich ist der scheinbare Widerstand kein einfaches Labormaterial-Reibungspaar, sondern ein „effektives“ feldskaliges Phänomen. Die unmittelbare Beobachtungsherausforderung besteht darin, eine nützliche Formel zu entwickeln, die als effektiver feldskaliger Abschluss fungiert und den unzugänglichen strukturellen Zustand auf Observablen wie Wärmefluss, effektive Reibungskoeffizienten, geodätische Kopplung und seismisches Momentenfreisetzungsverhalten projiziert. Diese Studie untersucht, ob eine „strukturelle Zustands-Effektiv-Reibungs-Architektur“, die zuvor auf andere unzugängliche großskalige Reibungsprobleme angewendet wurde, sinnvoll auf Subduktionsbeobachtungen projiziert werden kann.

Methodik
Die Studie verwendet ein zweigliedriges Beobachtungs-Audit, um den strukturellen Zustands-Abschluss gegen reale Daten zu testen:

1. Wärmefluss- und thermisches Modell-Audit:

   • Es wurde ein Datensatz aus literaturkodierten Randbedingungen zusammengestellt, der die effektiven Reibungskoeffizienten ($\mu'$) zusammenfasst, die aus Wärmefluss- und thermischen Modellen über verschiedene Subduktionszonen hinweg abgeleitet wurden.
   • Diese Randbedingungen deuten typischerweise auf niedrige effektive Reibungswerte hin (im Bereich von 0,03 bis 0,07), die signifikant niedriger sind als die trockene Labor-Byerlee-Reibung ($\mu \approx 0,6$).
   • Ein reduzierter struktureller Zustands-Effektiv-Reibungs-Abschluss wurde an diese Intervall-Randbedingungen angepasst und mit einer konstanten effektiven Reibungs-Baseline verglichen.

2. USGS Seismisches Moment-Audit:

   • Erdbebenkatalogdaten (M $\ge$ 7) von 1900 bis heute wurden für 12 repräsentative Subduktionsregionen über die USGS FDSN API heruntergeladen, was 715 Ereignisse ergab.
   • Das gesamte seismische Moment ($M_0$) wurde für jede Region als grobe Projektion der langfristigen Stick-Slip-Momentenfreisetzung berechnet.
   • Verschiedene Projektionsmodelle wurden gegen diese Momentendaten getestet, einschließlich:
     • Konvergenzgeschwindigkeit allein.
     • Thermischer Zustand allein.
     • Ein einfacher „struktureller Belastungsindex“.
     • Ein korrigiertes Modell, das Geschwindigkeit, thermischen Zustand und Plattenabtauchwinkel (Slab Dip) kombiniert.
   • Die Leistung wurde mittels Leave-One-Out (LOO) Kreuzvalidierung bewertet, um die regionale Übertragbarkeit zu beurteilen.

Der Strukturelle Zustands-Abschluss-Modell
Der vorgeschlagene Abschluss definiert die effektive Reibung ($\mu_{eff}$) als Funktion der Plattenkonvergenzgeschwindigkeit ($V$) und eines reduzierten thermisch-fluiden Schwächungszustands ($\Theta$). Das Modell ist als vierstufige Architektur strukturiert:

1. Physischer zulässiger Bereich: Aufgelegt durch die Subduktionsmechanik (z. B. $0 \le \mu_{eff} \le \mu_{Byerlee}$).
2. Strukturelle Koordinate: Eine zusammengesetzte Variable ($\Phi$), die Verhakung (Locking), Rauheit, Fluidgehalt, veränderte Minerale, Kriechen und Geometrie integriert.
3. Strukturelle Antwort: Eine Produktionsverlust-Form, die die Reaktion des Systems auf die strukturelle Koordinate darstellt.
4. Projektion: Abbildung des Zustands auf messbare Observablen (Wärmefluss, Kopplung, langsames Gleiten, seismisches Moment).

Hauptergebnisse

• Wärmefluss-Audit: Der strukturelle Zustands-Abschluss fällt erfolgreich innerhalb aller zusammengestellten effektiven Reibungsintervalle (0,03–0,07). Er liefert einen etwas geringeren Intervallverlust (0,00003544) als die beste konstante effektive Reibungs-Baseline (0,00004225). Beide Modelle lehnen den trockenen Laborreibungswert von 0,6 korrekt ab, was bestätigt, dass Subduktions-Observablen einen schwachen, durch Fluide und Temperatur modifizierten effektiven Koeffizienten repräsentieren.
• Seismisches Moment-Audit:
  • Versagen der Geschwindigkeit allein: Die Konvergenzgeschwindigkeit allein ist ein schlechter Prädiktor für das gesamte seismische Moment (LOO RMSE = 0,8293). Beispielsweise hat die Tonga-Kermadec-Zone eine hohe Konvergenzgeschwindigkeit, aber ein geringeres Gesamtmoment als langsamere Zonen wie Peru-Chile oder Japan.
  • Versagen des einfachen strukturellen Index: Ein einfacher struktureller Belastungsindex schnitt schlecht ab (LOO RMSE = 0,7210), was zeigt, dass die Formel nicht als universelles reines Geschwindigkeits-Reibungsgesetz behandelt werden kann.
  • Erfolg von Thermisch-Geometrisch: Modelle, die den thermischen Zustand und die Slab-Geometrie einbeziehen, schnitten signifikant besser ab. Die „Geschwindigkeit-Temperatur-Dip“-Projektion erreichte den niedrigsten LOO RMSE (0,6005) unter den Nicht-Katalog-Modellen und übertraf die Speed-only- und Thermal-only-Baselines. Das „structural_interaction“-Modell lieferte einen noch höheren Fehler (LOO RMSE = 0,8876).

Bedeutung und Ansprüche
Die Arbeit hält einen bescheidenen Umfang ein und stellt explizit klar, dass sie nicht beansprucht, thermomechanische Subduktionsmodelle, Rate-and-State-Reibungstheorien oder Erdbebenzyklus-Simulationen zu ersetzen. Stattdessen bietet sie eine Pilot-Validierung eines spezifischen architektonischen Ansatzes.

Die am stärksten gestützten Schlussfolgerungen sind:

1. Kompatibilität: Die effektive Reibung der Subduktion ist mit einem strukturellen Zustands-Abschluss kompatibel. Das Modell deckt erfolgreich die in Wärmeflussstudien beobachteten niedrigen Reibungsintervalle ab.
2. Notwendigkeit der Projektion: Die beobachtbare Projektion dieses Abschlusses muss den thermisch-fluiden Zustand und die Slab-Geometrie enthalten. Eine Formel, die nur auf Konvergenzgeschwindigkeit oder einem einfachen strukturellen Index basiert, reicht nicht aus, um die regionale seismische Momentfreisetzung zu erklären.
3. Architektur über universellem Koeffizienten: Die Forschung legt nahe, dass der „universelle“ Aspekt der Subduktionsreibung nicht ein fester numerischer Koeffizient ist, sondern eine Architektur: harter physisch möglicher Subduktionszustand $\rightarrow$ strukturelle Koordinate $\rightarrow$ Produktionsverlust-Antwort $\rightarrow$ thermisch-fluide-geometrische Projektion $\rightarrow$ Observablen.

Zusammenfassend lässt die Studie darauf schließen, dass ein struktureller Zustands-Effektiv-Reibungs-Abschluss ein tragfähiger Rahmen für Subduktionsgrenzflächen ist, sofern er durch thermisch-fluide und geometrische Projektionsterme korrigiert wird. Eine externe Universalität bleibt unbewiesen, und die vorliegenden Evidenzen stützen eher einen Pilot-Abschluss als ein definitives universelles Gesetz.