Kimberlit-Eruptionen, angetrieben durch Plattenfluss und Subduktionswinkel

Jul 12, 2023

Wissenschaftliche Berichte Band 13, Artikelnummer: 9216 (2023) Diesen Artikel zitieren

55 Zugänge

16 Altmetrisch

Details zu den Metriken

Kimberlite entstehen aus thermochemischen Auftrieben, die Diamanten an die Oberfläche der Erdkruste transportieren können. Der Großteil der auf der Erdoberfläche erhaltenen Kimberlite brach vor 250 bis 50 Millionen Jahren aus und wurde auf Veränderungen der Plattengeschwindigkeit oder Mantelplumes zurückgeführt. Diese Mechanismen können jedoch nicht das Vorhandensein starker Subduktionssignaturen erklären, die in einigen Kimberliten der Kreidezeit beobachtet wurden. Dies wirft die Frage auf, ob es einen Subduktionsprozess gibt, der unser Verständnis des Zeitpunkts von Kimberliteruptionen vereinheitlicht. Wir entwickeln eine neuartige Formel zur Berechnung des Subduktionswinkels basierend auf Grabenwanderung, Konvergenzrate, Plattendicke und -dichte, um den Zustrom von Plattenmaterial in den Mantel mit dem Zeitpunkt von Kimberlit-Eruptionen in Verbindung zu bringen. Wir stellen fest, dass Subduktionswinkel in Kombination mit Spitzen im Slab-Fluss Impulse von Kimberlit-Eruptionen vorhersagen. Hohe Subduktionsgeschwindigkeiten des Plattenmaterials lösen einen Mantelrückfluss aus, der fruchtbare Reservoire im Mantel stimuliert. Diese konvektiven Instabilitäten transportieren die von der Platte beeinflusste Schmelze in einem Abstand vom Graben zur Oberfläche, der dem Subduktionswinkel entspricht. Unsere Deep-Time-Slab-Dip-Formulierung bietet zahlreiche potenzielle Anwendungen, darunter die Modellierung der tiefen Kohlenstoff- und Wasserkreisläufe und ein verbessertes Verständnis subduktionsbedingter Mineralvorkommen.

Kimberlite sind mafische Vulkangesteine, die aus dem Erdmantel ausbrechen und die Wirtsgesteine ​​der meisten Diamanten sind1. Kimberlite kommen auf jedem Kraton vor und wurden seit 3 ​​Ga2 sporadisch eingelagert, doch die größte Anzahl der heute auf der Erde erhaltenen Kimberliteruptionen entstand in den letzten 250 bis 50 Millionen Jahren, hauptsächlich in Afrika und Nordamerika3. Während die Verteilung von Kimberliten mit den Rändern großer Provinzen mit geringer Scherwellengeschwindigkeit (LLSVPs)4 und Änderungen in der Winkelgeschwindigkeit der Platte3 in Verbindung gebracht wurde, erklären diese weder die Häufigkeit von Kimberliteruptionen noch angereicherte radiogene Isotopensignaturen, die auf eine Subduktion hinweisen Plattenkomponente in einigen Kimberlitpopulationen der Kreidezeit1,5. Es wurde vorgeschlagen, dass die steile Subduktion der ozeanischen Lithosphäre in den Mantel einen starken Mantelrückfluss und Magmatismusimpulse auslöst6. Doch trotz des theoretischen Zusammenhangs zwischen Vulkanausbrüchen und hohen Plattenflussraten7 haben Schwierigkeiten bei der Schätzung des Volumens und des Subduktionswinkels der ozeanischen Lithosphäre, die in alten Subduktionszonen recycelt wird, jede Korrelation mit Kimberliteruptionen vereitelt. Frühere Versuche, den Neigungswinkel von Subduktionsplatten zu charakterisieren, verwendeten eine multivariate Analyse der Eigenschaften der Subduktionszone, um nach Korrelationen zwischen Schlüsselparametern zu suchen8,9,10,11,12,13,14. Allerdings eignen sich diese Ansätze vor allem für die Reproduktion des heutigen Gesteinsabfalls und sind nur begrenzt auf Subduktionszonen anwendbar, die im Laufe der tiefen geologischen Zeit rekonstruiert wurden. Unter Verwendung eines aktuellen Modells zur Rekonstruktion tektonischer Platten15 und von Modellen der Plattenkühlung16,17,18 untersuchen wir hier die Abschätzung der Plattenneigung anhand einfacher plattenkinematischer Parameter, die die meisten Subduktionszonen auf der ganzen Welt charakterisieren, um die mögliche Rolle steil subduzierender Platten bei der Kontrolle zu untersuchen Kimberlit-Ausbrüche in Afrika und Nordamerika.

Die Tiefen der subduzierenden Platten, die aus dem Slab2-Modell19 ermittelt wurden und mit Kimberliten überlagert sind, sind in den letzten 250 Millionen Jahren ausgebrochen20. Grabensegmente haben in (i) Ozeanien die folgenden Abkürzungen: Ton, Tonga; Ker, Kermadec; NH, Neue Hebriden; Sol, Salomon; (ii) Südostasien: PNG, Papua-Neuguinea; Summe, Sumatra; Mar, Marianen; IZB, Izu-Bonin; Ryu, Ryūkyū; Mann, Manila; Ph, philippinisch; (iii) Asien: Mak, Makran; SJ, Südjapan; NJ, Nordjapan; Kur, Kurile; (iv) Europa: Hel, Helenic; Cal, Kalabrien; (v) Nordamerika: Al, Aleuten; Cas, Kaskaden; (vi) Mittelamerika: Mexiko, Mexiko; MAM, Mittelamerika; LAT, Kleine Antillen; (vii) Südamerika: EG, Ecuador; SA, Südamerika; SC, Südchile; SSW, Süd-Sandwich. Ungefähre Grabenstandorte, Grenztypen, Namen und entsprechende Abkürzungen sind in Tabelle S1 aufgeführt. Weiße Bereiche weisen auf Regionen mit nichtozeanischer Kruste hin; dicke rote Linien zeigen mittelozeanische Rücken an; Dünne rote Linien zeigen Transformationsgrenzen an. Die Karte wurde mit Cartopy21 erstellt.

Der Abwärtszug von Platten in den Erdmantel ist die größte treibende Kraft in der Plattentektonik22, wobei die Biegung der subduzierenden Platte eine wichtige Rolle bei der Modulation des Ausmaßes des Plattenzuges spielt, der auf die Oberfläche übertragen wird, um die Plattenbewegung anzutreiben23. Seit den 1980er Jahren8 haben zahlreiche Studien die Zusammenhänge zwischen dem Subduktionswinkel und einer Reihe von Parametern untersucht, darunter die Dauer der Subduktion, die Konvergenzrate, die Art der darüber liegenden Platte, die Grabenlänge und andere8,9,10,11,12,13, 14. Eine kürzlich durchgeführte Studie kehrte zu einer multivariaten Regression zurück und bekräftigte die seit langem vertretene Ansicht, dass die Plattenneigung von der Subduktionsdauer und in geringerem Maße vom Alter der absteigenden Platte und davon, ob die darüber liegende Platte kontinentaler oder ozeanischer Natur ist, beeinflusst wird13. Während diese Parameter für die Annäherung an die heutige Plattenneigung nützlich sind, ist die Abgrenzung von Subduktionszonensegmenten subjektiv und zeitabhängig, was eine Schwierigkeit darstellt, wenn versucht wird, diese Plattenneigungsbeziehungen in tiefe geologische Zeit mit sich ständig weiterentwickelnden Subduktionsgrenzen zurückzuverfolgen. Beispielsweise werden Daten über die Dauer oder den erneuten Beginn der Subduktion mit zunehmendem Alter immer spärlicher und Plattenreorganisationen machen es schwierig, die Langlebigkeit der Subduktion zu quantifizieren. Ebenso ist die Unterteilung von Plattensegmenten im Laufe der Zeit höchst subjektiv und ein Bruch in der Topologie der Subduktionszone kann zu stark unterschiedlichen Schätzungen der Plattenneigung führen. Um diese Herausforderungen zu meistern, untersuchen wir Plattenneigungskorrelationen mit einer multivariaten Analyse der Plattenrheologie und kinematischen Parameter, die aufgrund der Optimierung von Plattenmodellen für Nicht-Nettorotation über tiefe geologische Zeiträume hinweg weniger mehrdeutig sind24. Mit diesem Ansatz kann die Schätzung der Plattenneigung auf sich entwickelnde Subduktionszonen angewendet werden, die sich soweit zurück erstrecken, wie es eine bestimmte Plattenrekonstruktion zulässt.

Wir extrahieren aktuelle Plattenneigungsdaten aus dem Slab2-Modell19, das die Tiefe und Geometrie abtauchender Platten auf der ganzen Welt anhand von Erdbebentiefen und Tomographiemodellen schätzt (Abb. 1) und kombinieren diese mit nahezu aktuellen kinematischen Platteneigenschaften aus einem aktuellen Plattenrekonstruktionsmodell25 mit pyGPlates26. Der Neigungswinkel \(\theta\) der nach unten gerichteten Platte wird in mehreren Tiefenintervallen innerhalb des Grabens gemessen (Abb. 2). Der durchschnittliche Neigungswinkel \(\theta _{\textrm{av}}\) ist einfach das arithmetische Mittel aller vom Slab2-Modell aufgelösten Tiefenintervalle orthogonal zum Graben,

Dabei ist d das Tiefenintervall und n die Gesamtzahl der Tiefenintervalle, die vom Slab2-Modell aufgelöst werden können. Der gesamte Arbeitsablauf zur Berechnung des Neigungswinkels aus dem Slab2-Modell sowie die Plattenrheologie und die kinematischen Parameter der subduzierenden Platte sind auf GitHub (https://github.com/brmather/Slab-Dip) offen verfügbar. Die statistisch signifikantesten Kombinationen von Parametern, die empfindlich auf den aktuellen Neigungswinkel subduzierender Platten reagieren, werden in den folgenden Abschnitten beschrieben.

Schematische Darstellung von Platteneinbrüchen und Subduktionsparametern. Die Konvergenzrate \(v_c\) ist die Summe der Geschwindigkeit der subduzierenden Platte \(v_s\) und der Geschwindigkeit der übergeordneten Platte \(v_o\). Wir gehen davon aus, dass die Geschwindigkeiten der darüber liegenden Platte und des Grabens gleich und entgegengesetzt sind (\(v_o = -v_t\)). \(v_s\) und \(v_o\) sind positiv zum Graben und jeder Vektor ist orthogonal zum Graben. \(v_{\textrm{hsp}}\) ist die halbe Ausbreitungsrate an mittelozeanischen Rücken, die proportional zur Rate des flüchtigen Zuflusses (\(q_v\)) in die Platte ist, und \(\theta\ ) ist der Neigungswinkel der subduzierenden Platten, der in verschiedenen Intervallen berechnet wird.

Frühere Studien zur multivariaten Analyse von Subduktionskoeffizienten fanden keinen statistisch signifikanten Zusammenhang zwischen der Neigung der Platte und dem Alter der Subduktion der ozeanischen Lithosphäre8,9,13. Wir finden jedoch eine signifikante Beziehung, bei der die Plattenneigung proportional zur Konvergenzgeschwindigkeit \(v_c\) (\(P=0,39\)) und der Dicke der nach unten gerichteten Platte \(h_{\textrm{ Platte}}\) (\(P=0,59\)) (Abb. 3). Die Plattendicke wurde aus der thermischen Entwicklung der ozeanischen Lithosphäre vorhergesagt.

wobei \(\sqrt{\kappa t}\) aus Modellen der Plattenkühlung stammt16,27. Diese Modelle beschreiben die Verdickung der thermischen Grenzschicht als Funktion des Alters des Meeresbodens, die bei einem konstanten thermischen Diffusionskoeffizienten von \(\kappa = 1\) mm\(^2\)/ eine maximale Dicke von etwa 80 Myr erreicht. S. Wir verwenden ein Plattenmodell der ozeanischen Lithosphärenkühlung, das nachweislich eine optimale Anpassung an Tiefen- und Wärmeflussdaten liefert, um die Dicke der ozeanischen Lithosphäre zu berechnen, die in Gräben auf der Erde recycelt wird16. Das Produkt aus der Plattendicke \(h_{\textrm{Platte}}\) und der Konvergenzgeschwindigkeit \(v_c\), integriert entlang eines Grabensegments, ergibt die volumetrische Rate des lithosphärischen Recyclings oder den „Plattenfluss“. Wir beproben Altersgitter des Meeresbodens an Grabengrenzen, um die Dicke der subduzierenden ozeanischen Lithosphäre zu berechnen.

Korrelationen zwischen Plattenneigung und Konvergenzgeschwindigkeit, Plattendicke und Ausbreitungsgeschwindigkeit. 2D-Histogramme veranschaulichen die Wahrscheinlichkeitsdichte zwischen den einzelnen Parametern mit Plattenneigung für alle Subduktionszonensegmente, die in Schritten von 0,5 Grad mosaikiert sind (a–c); Subduktionszonensegmente werden für jede Subduktionszone gruppiert, wobei Whisker eine Standardabweichung vom Mittelwert (d–f) angeben. \(P_r\) ist der Wert des Pearson-Korrelationskoeffizienten und P ist der p-Wert.

Der zweite wichtige zu berücksichtigende Parameter ist die Geschwindigkeit des Grabenvorschubs oder -rückzugs. Es wird allgemein angenommen, dass das Zurückrollen der Platte zu einer Subduktion mit geringem Winkel und einer Stagnation der Platte an der oberen und unteren Mantelgrenze führt12,28,29, was teilweise erklärt, warum geringere Subduktionswinkel häufiger in Subduktionszonen zwischen zwei ozeanischen Platten beobachtet werden, wo es zu einem Rückzug der Gräben kommt häufiger als in Subduktionszonen zwischen Ozean und Kontinent10. Die Geschwindigkeit der Grabenwanderung, \(v_t\), wird relativ zum Bezugssystem des Mantels berechnet und kann mit der Konvergenzgeschwindigkeit, \(v_c\), verglichen werden, um verschiedene Subduktionsdynamiken zu charakterisieren:

Wenn \(v_c = -v_t\), wird die gesamte Konvergenzrate auf Rollback aufgeteilt;

Wenn \(v_t = 0\) ist der Graben stationär und die Geschwindigkeit der subduzierenden Platte ist gleich der Konvergenzrate (\(v_s = v_c\));

Wenn \(v_t > 0\), schreitet der Graben in Richtung der Subduktion voran.

Ein dritter Parameter, den wir berücksichtigt haben, ist die flüchtige Anreicherung der subduzierenden Platte. Eine erhöhte Häufigkeit flüchtiger Komponenten wie Wasser und Kohlenstoff verbessert die Kopplung zwischen der subduzierenden Platte und der darüber liegenden Platte30. Die Quelle der Anreicherung flüchtiger Stoffe liegt an mittelozeanischen Rücken, wo flüchtige Schmelze durch Kanäle in der neu entstehenden ozeanischen Lithosphäre zirkuliert31. Die Sequestrierungsrate flüchtiger Stoffe in der ozeanischen Lithosphäre, \(q_v\), ist proportional zur Ausbreitungsrate des Meeresbodens, \(v_{\textrm{hsp}}\)31. Wir beproben Gitter der Ausbreitungsrate des Meeresbodens, die mit den Arbeitsabläufen in32 erstellt wurden, an Grabengrenzen auf die gleiche Weise, wie wir Altersgitter interpolieren, um die Plattendicke \(h_{\textrm{plate}}\) zu berechnen. Wir stellen fest, dass \(v_{\textrm{hsp}}\) eine starke negative Korrelation mit der Plattenneigung aufweist (\(P=0,34\), Abb. 3c,f). Flüchtige Stoffe werden der Platte auch durch andere Prozesse hinzugefügt, beispielsweise durch hydrothermale Veränderung der oberen Kruste33, Sepertinisierung während der ultralangsamen Ausbreitung34 und während der Biegung und Rissbildung der Platte vor dem Eintritt in die Subduktionszone35. Die oben skizzierte Ausbreitungsgeschwindigkeitsbeziehung stellt jedoch die Beziehung dar Die häufigste Quelle flüchtiger Stoffe, die weltweit in den meisten Gräben subduziert wird.

Ein vierter Parameter, den wir berücksichtigten, war die mittlere Dichte der subduzierenden Platte. Der lithosphärische Auftrieb kann anhand des Plattenalters, der thermischen Struktur, der Krustendicke und der Erschöpfung abgeschätzt werden36. Die Kruste ist aufgrund ihrer relativ geringen Dichte (\(\sim\) 2900 kg/m\(^3\))37 die Hauptquelle für positiven Auftrieb in der ozeanischen Lithosphäre. Unter der Annahme moderner Manteltemperaturbedingungen weist die 10–20 Millionen Jahre alte Lithosphäre einen negativen Auftrieb auf36 und wird mit zunehmendem Alter immer negativer. Ein 60 Millionen Jahre langer Abschnitt der ozeanischen Lithosphäre ist etwa 79,4 km dick16, davon sind 7 km Kruste, was einer mittleren Dichte von \(\rho _{\textrm{av}} = 3278\) kg/m\(^3 entspricht \). Während der Auftrieb ozeanischer Platten möglicherweise nicht ausreicht, um allein eine Subduktion auszulösen38, verändern seitliche Änderungen in der Dichtestruktur der absinkenden Platte die Auftriebskraft in etablierten Subduktionszonen, was sich auf die Plattenneigung auswirken kann. Solche positiven Auftriebsanomalien sind mit ozeanischen Hochebenen verbunden, die häufig Subduktionszonen verstopfen39 oder zur Subduktion flacher Platten führen40.

Unter Berücksichtigung all dieser rheologischen und kinematischen Beziehungen für heutige Subduktionszonen haben wir eine Nearest-Neighbor-Regression angewendet, um den Neigungswinkel \(\theta _{\textrm{av}}\) subduzierender Platten vorherzusagen. Diese Regression implementiert einen kd-Baum, um k Nachbarn mit dem kürzesten euklidischen Abstand vom Trainingsdatensatz \(X_{\textrm{train}}\) zum Testdatensatz \(X_{\textrm{test}}\) effizient zu suchen. und verwendet den gewichteten Mittelwert der entsprechenden Plattenneigungen, um die Testplattenneigung vorherzusagen.

Dabei ist \(d_k\) der euklidische Abstand zwischen den Trainings- und Testdaten für k nächste Nachbarn, \(d_k = \Vert X_{\textrm{train}} - X_{\textrm{test}} \Vert _k\). Unter Verwendung einer Teilmenge der aktuellen Konfiguration von Subduktionsgrenzen und Platteneinbrüchen, die aus dem Slab2-Modell erhalten wurden, als Testdatensatz vergleichen wir die Leistung der Platteneinbruchvorhersage mit dem Trainingsdatensatz (Abb. 4). Der Trainingswert (\(R^2\)-Wert) misst die Übereinstimmung zwischen dem Trainings- und dem Testdatensatz, der für \(k=1\) maximal 1 beträgt. Dies liegt daran, dass zum aktuellen Zeitpunkt eine exakte Übereinstimmung zwischen dem Testdatensatz und dem Trainingsdatensatz besteht. Wenn k größer wird, werden mehr Nachbarn in den Durchschnitt einbezogen, was den \(R^2\)-Wert und die Kreuzvalidierungsbewertung verringert. Durch die Kombination dieser beiden Metriken wird die Leistung des Schätzers so bewertet, dass das Problem nicht über- oder unterangepasst ist. Wir entscheiden uns für \(k=5\), wo ein optimaler Kompromiss zwischen der Kreuzvalidierung und den Trainingsergebnissen besteht.

Kreuzvalidierungs-Score und Trainings-Score (\(R^2\)-Wert) für den Nearest-Neighbours-Regressionsalgorithmus zur Vorhersage des Neigungswinkels subduzierender Platten. Der Algorithmus wird unter Verwendung einer zunehmenden Anzahl von nächsten Nachbarn (k) getestet, die in der Berechnung aus Gleichung (1) verwendet werden. 3. Schattierte Bereiche geben eine Standardabweichung vom Mittelwert an (durchgezogene Linien). Der Kreuzvalidierungs-Score wird aus Gründen der visuellen Klarheit als negative Zahl angezeigt, um den optimalen Kompromiss mit dem Trainings-Score vergleichen zu können.

Wir haben ein flexibles objektorientiertes Python-Paket entwickelt, um die Plattenneigung auf der Grundlage heutiger Plattenrheologie und kinematischer Parameter abzuschätzen (https://github.com/brmather/Slab-Dip). Der Standard-Regressionsalgorithmus und der Trainingsdatensatz wurden oben beschrieben. Der Benutzer kann jedoch auch seinen eigenen Trainingsdatensatz und einen beliebigen Regressionsalgorithmus definieren, der im Scikit-Learn-Python-Paket gebündelt ist, um einen maßgeschneiderten Schätzer für die Neigung der Bodenplatte zu erstellen. Im GitHub-Repository stehen dem Benutzer mehrere Beispiele als Jupyter-Notebooks sowie Installationsanweisungen zur Verfügung. Die von uns gewählte Regression der nächsten Nachbarn ist allgemein und stellt einen robusten Zusammenhang zwischen den heutigen Eigenschaften subduzierender Platten und denen im Laufe der Zeit her. Obwohl es möglicherweise nicht auf bestimmte geodynamische Kontexte anwendbar ist, die sich von den heutigen unterscheiden, fasst es einige der Haupttreiber der Subduktion zusammen, die aus der plattentektonischen Theorie vorhergesagt werden, wie z. B. die Beziehung zwischen dem Alter des Meeresbodens, der Plattendicke, der Grabenwanderung usw. und Konvergenzgeschwindigkeit. Dabei wird berücksichtigt, dass die Plattensaugkraft, die mit älteren dichten Platten verbunden ist, die stark in den Erdmantel absinken, die Geschwindigkeit der tektonischen Platten an der Erdoberfläche antreibt und eine Zunahme des Ausmaßes des Zurückrollens zu einer Abflachung der nach unten gerichteten Platten führt. gehende Platte.

Es ist wichtig zu beachten, dass die von uns formulierte Beziehung das Einfallen der Platte bis zur maximalen Tiefe vorhersagt, die vom Slab2-Modell aufgelöst wird. Die Flugbahn der absteigenden Platte in den Erdmantel kann von den vorhergesagten Plattentiefen abweichen, da sie auf Viskositätskontraste im Erdmantel trifft, was manchmal zu Plattenstagnation41,42 und möglicherweise zur Plattenverankerung43 führt. Diese Dynamik wird von unserer Regression nicht erfasst, was ihre Anwendung auf bestimmte Subduktionszonen über einen nominellen Zeitbereich ausschließen kann. Dennoch gilt unsere Slab-Dip-Formulierung für die meisten Subduktionszonen weltweit und bietet den Vorteil, dass sie über tiefe geologische Zeiträume hinweg angewendet werden kann, indem tektonische Plattenrekonstruktionen verwendet werden, um die Flugbahn subduzierender Slabs in den Mantel vorherzusagen. Dies kann verwendet werden, um den Abstand zwischen Gräben und Vulkanbögen, das Auftreten flacher Plattensubduktion und das Recycling ozeanischer Lithosphäre in tiefere Teile des Mantels abzuschätzen.

Kimberlite sind Vulkangesteine, die schnell aus dem Erdmantel aufsteigen und weltweit in Kratonen eingelagert sind2. Kimberlit-Eruptionen wurden mit Mantelauftrieben aus dem LLSVP4 und Dehnungstektonik im Zusammenhang mit lithosphärischer Entladung44 oder Änderungen der Plattengeschwindigkeit3 in Verbindung gebracht. Diese Mechanismen erklären jedoch nicht die starken subduzierten Slab-Signaturen, die in Kimberliten aus der Kreidezeit in Afrika, Brasilien und Nordamerika aufgrund erhöhter Strontium-Isotopen-Verhältnisse beobachtet wurden1,5. Während afrikanische Kimberlite einen statistisch signifikanten Zusammenhang mit der Entfernung zu LLSVPs aufweisen, besteht bei Kimberliten in Nordamerika kein solcher Zusammenhang (Abb. 5). Im Gegensatz dazu wurden Kimberliteruptionen in Nordamerika mit einer flachen Subduktion der Farallon-Platte während der Laramid-Orogenese in Verbindung gebracht45. Hier könnte Magma durch wasserdurchströmtes Dekompressionsschmelzen der Mantelübergangszone46 erzeugt und durch subduktionsinduzierten Mantelrückfluss47 nach oben transportiert worden sein. Ein hoher Slab-Fluss wurde zuvor mit der Häufigkeit von Vulkanausbrüchen in Verbindung gebracht, bei denen Mantelauftriebe durch die Subduktion großer Mengen ozeanischer Lithosphäre in den Mantel vorangetrieben werden48. Subduktionsbedingte Mantelaufschwünge wurden mit der Bildung der 260 Ma großen magmatischen Provinz Emeishan im Südwesten Chinas49 und dem kenozoischen Vulkanismus im Nordosten Chinas50 in Verbindung gebracht. Um die Rolle der Subduktion bei der Erzeugung unterschiedlicher Kimberlitpopulationen in Afrika und Nordamerika in Einklang zu bringen, haben wir die vertikalen und horizontalen Komponenten des Plattenflusses mithilfe unserer Formulierung des Plattenabfalls im vorherigen Abschnitt getrennt und die globalen Subduktionsgrenzen mithilfe von pyGPlates26 rekonstruiert. Wir verwendeten ein von 15 modifiziertes 170-Myr-Plattenmodell mit deutlich verbesserten Subduktionszonengrenzen entlang des Westrandes Nordamerikas und einer verbesserten Auflösung der karibischen Platte51 (siehe Methoden). Die Standorte der Kimberlite (\(\le\) vor 170 Ma wurden mithilfe einer Zusammenstellung von Kimberliten mit Eruptionsalter52, die wir in ein dreifach verfeinertes ikosoedrisches Netz umwandelten,53 bis zu ihrer Eruptionszeit rekonstruiert, um Duplikate und geografische Stichprobenverzerrungen zu vermeiden. Anschließend haben wir die Kimberlite in (i) nordamerikanische und (ii) afrikanische Populationen unterteilt. Zusammengenommen machen diese Kimberlitpopulationen 91 % des globalen Kimberlitdatensatzes aus, der in den letzten 250–50 Millionen Jahren ausgebrochen ist und als „Kimberlitblüte“3 bezeichnet wird. Beim Vergleich des Zeitpunkts von Kimberlit-Eruptionen mit den Geschwindigkeiten des Plattenflusses berücksichtigen wir nur Gräben, in denen die Subduktion in Richtung der Kimberlit-Populationen erfolgt. Wir finden eine starke Korrelation zwischen dem hohen Plattenfluss entlang des westlichen Randes Nord- und Mittelamerikas, der mit der Subduktion der Farallon-Platte verbunden ist, und sowohl afrikanischen als auch nordamerikanischen Kimberlitpopulationen (Abb. 6a). Der niedrige Subduktionswinkel (30–35\(^\circ\)), der aus unserer Plattenneigungsanalyse an rekonstruierten Subduktionszonen vorhergesagt wurde, deutet darauf hin, dass sich Platten über mehr als 1.000 km aus dem Graben erstrecken würden, bevor sie die 660 km lange Mantelübergangszone kreuzen, und eindringen könnten tiefer in den unteren Erdmantel. Eine schnelle Subduktion von Plattenmaterial kann zu einem Mantelrückfluss führen35, aus dem Mantelauftriebe Kimberliteruptionen auslösen können. In den folgenden Abschnitten untersuchen wir, wie der subduktionsinduzierte Mantelrückfluss aufgrund hoher Plattenflussraten mit Kimberliteruptionen in Afrika und Nordamerika zusammenhängen könnte.

Räumlich-zeitlicher Zusammenhang zwischen Kimberliteruptionen und LLSVP-Grenzen. (a) Scherwellengeschwindigkeitsanomalie aus dem SMEAN2-Tomographiemodell54, überlagert mit der räumlichen Verteilung von Kimberliten (Diamanten), rekonstruiert zu ihrem Eruptionsalter. Wir definieren LLSVP-Grenzen als die 1 % langsame Kontur des 2800 km tiefen Schnitts aus der Geschwindigkeitsschwankung, die durch die rote gestrichelte Linie dargestellt wird. Als Referenz werden aktuelle Küstenlinien hinzugefügt. Die Karte wurde mit Cartopy21 erstellt. (b) Kumulative Dichtefunktion des Abstands zwischen LLSVP-Grenzen und Kimberlitpopulationen in Afrika und Nordamerika mit zufälligen kontinentalen Standorten zum heutigen Tag und rekonstruiert auf 170 Ma, als der Pangäa-Superkontinent entstand. Afrikanische Kimberlite weisen eine statistisch signifikante Beziehung zu den LLSVP-Grenzen auf, während dies bei nordamerikanischen Kimberliten nicht der Fall ist.

Es besteht eine starke Korrelation zwischen dem Slab-Fluss und der Häufigkeit von Kimberliteruptionen während des Höhepunkts afrikanischer Kimberliteruptionen zwischen 120 und 130 Ma (Abb. 6b). Die Subduktion hielt am Westrand des Superkontinents Pangäa während seiner Entstehung an, bis vor etwa 120 Millionen Jahren die Spaltung zwischen Afrika und Südamerika begann (Abb. 7). Zwischen 160 und 120 Ma korrelieren zwei Spitzen bei Kimberlit-Eruptionen mit Impulsen hoher Plattenströmung von der subduzierenden Farallon-Platte unter Amerika, die 30–35 \(^\circ\) in den Erdmantel eintauchte. Der größere Höhepunkt der Kimberliteruptionen bei 120 Ma entspricht einem Plattenfluss von 60 km\(^3\)/Jahr. Ein zweiter Höhepunkt afrikanischer Kimberliteruptionen trat zwischen 80 und 90 Ma auf, was mit einem zweiten Impuls im Plattenfluss (bis zu 80 km\(^3\)/Jahr) und einem Maximum der Plattengeschwindigkeit (6 cm/Jahr) korreliert Die Ausbreitungsgeschwindigkeit des Meeresbodens nimmt zwischen Afrika und Südamerika zu (Abb. 6b). Während gezeigt wurde, dass mit dem LLSVP verbundene Mantelplumes einen erheblichen Teil der kratonischen Lithosphäre in Afrika erodiert haben55 und möglicherweise mit einigen Kimberlit-Eruptionen in Zusammenhang stehen4, erklärt dies weder die Subduktionssignaturen in Kimberliten noch den Zeitpunkt ihres Ausbruchs.

Wir schlagen vor, dass ein Reservoir recycelter Platten den Mantel aufgrund der allgegenwärtigen Subduktion während des Aufbaus von Pangäa besetzt, was zu einem Dehydrierungsschmelzen des darüber liegenden Mantels führt, da in den kalten Platten mitgeführtes Wasser freigesetzt wird56. Dann, wenn Pangäa aufzubrechen beginnt, treibt die schnelle Subduktion von Plattenmaterial in einem geringen Winkel den Mantelrückfluss aus diesem fruchtbaren Mantelreservoir voran, was zu Kimberlitausbrüchen führt. Da subduzierende Platten die tiefe Mantelstruktur beeinflussen und möglicherweise einen verstärkten Plume-Fluss an den Rändern des afrikanischen LLSVP7 auslösen, könnte dies die Lieferung von durch Platten beeinflusster Schmelze vom obersten unteren Mantel an die Oberfläche beschleunigen, was einer erhöhten Häufigkeit von Kimberlit-Eruptionen bei 120 entspricht Ma. Dies ähnelt möglicherweise einem Prozess, der zur Bildung der 260 Millionen Jahre alten magmatischen Provinz Emeishan beiträgt, in der der recycelte Paläo-Tethys-Ozean unter Südwestchina großflächige Mantelauftriebe aus Tiefen von 410–660 km auslöst49. Der zweite Impuls bei Kimberlit-Eruptionen vor 80–90 Ma kann aufgrund der großen Entfernung von der nächsten Subduktionszone nach der Öffnung des Südatlantiks nicht ohne weiteres mit dem Plattenfluss in Verbindung gebracht werden (Abb. 7). Stattdessen wäre ein Anstieg der afrikanischen Plattengeschwindigkeit erforderlich Setzen Sie mehr kratonische Lithosphäre den Mantelauftrieben aus, die mit dem LLSVP verbunden sind, und erhöhen Sie dadurch die Häufigkeit von Kimberlit-Eruptionen3,4 (Abb. 6b).

Zusammenhang zwischen Kimberlit-Eruptionsdichte, Plattenfluss und Plattengeschwindigkeit über 170 Millionen Jahre. (a) Impulse in kombinierten Kimberlit-Eruptionen korrelieren mit Perioden hoher Plattenströmung der nach Osten einfallenden Farallon-Platte am westlichen Rand Nordamerikas. Der aus unserer Analyse vorhergesagte niedrige Neigungswinkel der Subduktion (30–35\(^\circ\)) wird durch eine starke Seitwärtskomponente des Plattenflusses angezeigt. (b) der erste Höhepunkt afrikanischer Kimberlit-Eruptionen korreliert mit einem hohen Slab-Fluss zwischen 120 und 130 Ma, als sich Pangäa ausbreitete; Der zweite Höhepunkt bei 80–90 Ma lässt sich aufgrund der größeren Entfernung zwischen Amerika und Afrika ab der Öffnung des Südatlantiks eher durch einen Anstieg der afrikanischen Plattengeschwindigkeit3 als durch den Plattenfluss erklären. (c) Pulse in nordamerikanischen Kimberlit-Eruptionen korrelieren eng mit Pulsen im Plattenfluss, der entlang nach Osten abfallender Subduktionszonen innerhalb eines Radius von 2500 km um den nächsten Kimberlit-Ausbruch in Nordamerika integriert ist, ohne dass eine Korrelation mit der Plattengeschwindigkeit beobachtet wurde.

Entwicklung des Plattenneigungswinkels in globalen Subduktionszonen von 160 Ma bis heute, überlagert mit der räumlichen Verteilung von Kimberliteruptionen. Ein großer Höhepunkt des Plattenflusses und der Kimberlit-Eruptionsdichte tritt bei 120 Ma während des Auseinanderbrechens des Gondwana-Superkontinents auf, und ein zweiter (kleinerer) Höhepunkt bei 80 Ma ist mit der Entfernung der flachen Farallon-Platte von der Basis der darüber liegenden kontinentalen Lithosphäre im Norden verbunden Amerika. Weiße Regionen weisen auf nichtozeanische Kruste hin, graue Regionen auf heutige Küstenlinien. Karten wurden mit Cartopy21 erstellt. Eine vollständige Zeitreihe von 170 Ma bis 0 Ma ist als Animation in Film S1 verfügbar.

A second population of kimberlite eruptions occurred between 110 and 40 Ma while North America migrated westward during the opening of the North Atlantic Ocean. It has been proposed that the dehydration of hydrous minerals stored within the flat-subducting Farallon plate promoted magmatism and kimberlite generation approximately 1500 km from the nearest trench45, however, geodynamic models suggest that flat subduction inhibits arc magmatism as the release and convection of fluids from the slab are obstructed by the asthenospheric wedge57. From our reconstructions of slab dip, the average dip angle along the western margin of North America varies between 30 and 36\(^\circ\) and the slab flux predicts the peaks and troughs in kimberlite eruption frequency between 110 and 40 Ma (Fig. 6c). Slab dip is spatially and temporally variable along North American subduction boundaries during the Laramide period, which has been attributed to the flat subduction of the Shatsky Rise conjugate on the northernmost section of the Farallon plate40. Its subduction predicts the distribution of magmatic and amagmatic zones in North America. From 95 to 60 Ma, the subduction of relatively young seafloor (5–50 Ma) combined with subduction of the buoyant conjugate Shatsky Rise leads to flat slab subduction beneath central USA58 (Fig. 7). The distribution of kimberlite eruptions during this period are focused in Canada and the south of North America on either side of the conjugate (Fig. 8). Abrupt changes in subduction angles could be accommodated by slab tears adjacent to the Arizona–New-Mexico magmatic belt57. It is likely that melts associated with the dehydration of recycled slab material in the mantle transition zone were delivered to the surface through subduction-induced return flow47. Removal of the flat Farallon slab from the base of overriding continental lithosphere at 50 Ma2.3.CO;2 (1995)." href="/articles/s41598-023-36250-w#ref-CR59" id="ref-link-section-d64025640e2584">59 würde den Mantelrückfluss weiter stimulieren und zu ausgedehnteren Kimberliteruptionen führen, die in der ehemals amagmatischen Zone im Zentrum der USA auftreten (Abb. 6c).

Entwicklung der Flachplattensubduktion entlang Nordamerikas. (a) Die flache Slab-Subduktion der Farallon-Platte wird durch den Shatsky Rise verursacht, der bei 95 Ma in den Graben eintritt. (b) Ein Großteil des Shatsky-Konjugats wird bis 70 Ma verbraucht. Zu diesem Zeitpunkt überdeckt der Farallon einen Großteil der zentralen USA und erzeugt eine tektonische Hebung, die mit der Laramid-Orogenese verbunden ist. (c) Bei 50 Ma wird die Farallon-Platte von der Basis der darüber liegenden kontinentalen Lithosphäre entfernt, was zu ausgedehnten Kimberlit-Eruptionen führt. Schwarze Polygone zeigen rekonstruierte große magmatische Provinzen an, graue Polygone zeigen abgeleitete LIP-Konjugate an, Pfeile zeigen die absolute Plattengeschwindigkeit an. Karten wurden mit Cartopy21 erstellt.

Die Dichotomie der Kimberlit-Eruption zwischen afrikanischen und nordamerikanischen Populationen ist wichtig, da der Rückflussmechanismus im Mantel zwar konsistent ist, die Stimulation der Quellregionen im Mantel jedoch unterschiedlich ist. In Afrika hat der mit subduzierenden Platten verbundene Mantelrückfluss wahrscheinlich Aufschwünge an den Rändern des LLSVP angeregt und ein fruchtbares Mantelreservoir belebt, das aus sinkenden Plattenresten stammt, die beim Zusammenbau von Pangäa übrig geblieben sind, was die in diesem Kimberlit aus der Kreidezeit beobachteten Subduktionssignaturen erklärt Bevölkerung. Der zweite Impuls afrikanischer Kimberlit-Eruptionen vor 80–90 Ma ist wahrscheinlich mit einem Anstieg der Plattengeschwindigkeit verbunden, da Afrika über mit dem LLSVP3 verbundene Aufschwünge wandert. Mittlerweile werden Kimberliteruptionen in Nordamerika durch den Rückfluss des oberen Erdmantels in Regionen neben der flachen Subduktion des Shatsky Rise und innerhalb der Region, die von der Laramid-Orogenese nach der Entfernung der flachen Platte unter Nordamerika betroffen ist, vorangetrieben. Wichtig ist, dass beide Kimberlitpopulationen mit der schnellen Subduktion der Farallon-Platte in einem geringen Winkel unter Amerika zusammenhängen, was darauf hindeutet, dass diese Platte eine wichtige Rolle bei der Förderung des durch Auftrieb verursachten Vulkanismus spielt. Ob die Farallon-Platte mit einer hohen Anreicherung flüchtiger Stoffe wie H\(_2\)O durchdrungen ist, was die Solidustemperatur senkt und teilweises Schmelzen und Magmabildung fördert45, oder ob es bereits vorhandene fruchtbare Mantelreservoirs stimuliert, ist unklar. Dennoch unterstreicht diese Studie die Bedeutung der Subduktion für die Entstehung von Kimberliten und stellt frühere Vorstellungen in Frage, dass Kimberlite hauptsächlich durch Mantelplumes erzeugt werden.

Der Neigungswinkel der subduzierenden ozeanischen Lithosphäre ist ein Schlüsselparameter, der die Mantel- und Kontinentaldynamik in Subduktionszonen charakterisiert. Wir schlagen einen einfachen Rahmen zur Vorhersage des Plattenabfalls anhand der Dicke der nach unten gerichteten Platte, der Konvergenzrate, der Geschwindigkeit der Grabenwanderung, der Dichte und der flüchtigen Anreicherung der Platte vor. Die Anwendung dieses Rahmenwerks auf Plattenrekonstruktionen liefert neue Einblicke in die Dynamik vergangener Subduktionszonen, die räumliche Verteilung des Bogenvulkanismus über tiefe geologische Zeiträume und das Schicksal subduzierter Platten. Mithilfe dieses Vorhersagerahmens rekonstruieren wir den Plattenneigungswinkel von Subduktionszonensegmenten über die letzten 170 Millionen Jahre, um Impulse bei Kimberliteruptionen zu erklären. Hohe Subduktionsraten stimulieren die Rückströmung des Mantels, was das teilweise Schmelzen und die Entstehung von Magma fördert. Kimberlit-Eruptionen in Afrika und Nordamerika stehen im Zusammenhang mit der Subduktion der Farallon-Platte unter Amerika. In Afrika zeigen Spitzen bei Kimberlit-Eruptionen eine starke Korrelation mit einem hohen Plattenfluss während der Anfangsstadien des Aufbrechens des Superkontinents und hohen Plattengeschwindigkeiten (bis zu 6 cm/Jahr), wenn Afrika über das LLSVP wandert. In Nordamerika führt die Subduktion des Shatsky Rise-Konjugats von 95–50 Ma zu einer flachen Subduktion unterhalb der zentralen USA, verbunden mit der Laramid-Orogenese, die den Magmatismus auf die Ränder der flachen Platte beschränkt, bis sie sich danach von der Basis der Lithosphäre entfernt 50 Ma. Unsere Ergebnisse unterstreichen die wichtige Rolle des Subduktionswinkels bei der Modulation des horizontalen Plattenflusses und damit der Verteilung und des Zeitpunkts des Vulkanismus. Dies hilft, die Dichotomie der Kimberlit-Populationen in Afrika und Nordamerika zu erklären und hat wichtige Auswirkungen auf die Bildung alter Kimberlite in Australien, Indien und Südamerika, die durch Rekonstruktionen der Plattenneigung durch tiefe geologische Zeiträume erklärt werden können.

Das in dieser Arbeit verwendete Plattenmodell wurde gegenüber einem kürzlich veröffentlichten Modell15 wie folgt modifiziert. Subduktionszonen, denen keine Bewegungen zugewiesen wurden, waren zunächst im Laufe der Zeit stationär; Diesen wurden neue Bewegungen relativ zum globalen Plattenkreis zugewiesen, um einen moderaten Rückzug des Grabens zu modellieren und gleichzeitig die tomografischen Einschränkungen einzuhalten. Die Orcas-Platte wurde vor 170–130 Ma in zwei separate Platten aufgeteilt, entsprechend ihrer Konfiguration nach 130 Ma, um Divergenzen an den Plattengrenzen auszugleichen. Die karibische Platte wurde außerdem vor 140–120 Ma durch ein neues Back-Arc-Ausbreitungszentrum gespalten, um die Bildung der großen magmatischen Provinz der Karibik an einem sich ausbreitenden Rücken zu ermöglichen51. Schließlich wurde das absolute Plattenbewegungsmodell mithilfe eines iterativen Optimierungsworkflows24 eingeschränkt. Ein ZIP-Archiv mit Plattenrekonstruktionsdateien zur Verwendung in GPlates ist auf Zenodo verfügbar (doi.org/10.5281/zenodo.5769002).

Die während der aktuellen Studie generierten und/oder analysierten Datensätze sind im Zenodo-Repository verfügbar, https://doi.org/10.5281/zenodo.5831990. Eine Reihe von Jupyter-Notebooks mit einem Python-Workflow zur Berechnung der Plattenneigung mithilfe eines Plattenrekonstruktionsmodells ist auch bei Zenodo über https://doi.org/10.5281/zenodo.5831990 und GitHub (https://github.com/brmather/) verfügbar. Slab-Dip). Die Plattenrekonstruktionssoftware GPlates und pyGPlates sind kostenlos verfügbar unter www.gplates.org/download.

Woodhead, J. et al. Kimberlite offenbaren die 2,5 Milliarden Jahre währende Entwicklung eines tiefen, isolierten Mantelreservoirs. Natur 573, 578–581. https://doi.org/10.1038/s41586-019-1574-8 (2019).

Artikel ADS CAS PubMed Google Scholar

Shirey, SB & Richardson, SH Beginn des Wilson-Zyklus bei 3 Ga, dargestellt durch Diamanten aus dem subkontinentalen Mantel. Wissenschaft 333, 434–436. https://doi.org/10.1126/science.1206275 (2011).

Artikel ADS CAS PubMed Google Scholar

Tappe, S., Smart, K., Torsvik, T., Massuyeau, M. & de Wit, M. Geodynamik von Kimberliten auf einer abkühlenden Erde: Hinweise auf die plattentektonische Entwicklung und tiefe volatile Zyklen. Planet Erde. Wissenschaft. Lette. 484, 1–14. https://doi.org/10.1016/j.epsl.2017.12.013 (2018).

Artikel ADS CAS Google Scholar

Torsvik, TH, Burke, K., Steinberger, B., Webb, SJ & Ashwal, LD Diamanten, die von Plumes aus der Kern-Mantel-Grenze entnommen wurden. Natur 466, 352–355. https://doi.org/10.1038/nature09216 (2010).

Artikel ADS CAS PubMed Google Scholar

Tappe, S., Graham Pearson, D., Kjarsgaard, BA, Nowell, G. & Dowall, D. Eingabe der Mantelübergangszone in den Kimberlit-Magmatismus in der Nähe einer Subduktionszone: Ursprung der anomalen Nd-Hf-Isotopensystematik am Lac de Gras, Kanada . Planet Erde. Wissenschaft. Lette. 371–372, 235–251. https://doi.org/10.1016/j.epsl.2013.03.039 (2013).

Artikel ADS CAS Google Scholar

Faccenna, C. et al. Durch Subduktion ausgelöste magmatische Impulse: eine neue Klasse von Plumes? Planet Erde. Wissenschaft. Lette. 299, 54–68. https://doi.org/10.1016/j.epsl.2010.08.012 (2010).

Artikel ADS CAS Google Scholar

Cao, X., Flament, N., Bodur, Ö. F. & Müller, RD Die Entwicklung der basalen Mantelstruktur als Reaktion auf die Aggregation und Ausbreitung von Superkontinenten. Wissenschaft. Rep. 11, 1–16. https://doi.org/10.1038/s41598-021-02359-z (2021).

Artikel CAS Google Scholar

Jarrard, RD Beziehungen zwischen Subduktionsparametern. Rev. Geophys. 24, 217–284. https://doi.org/10.1029/RG024i002p00217 (1986).

Artikel ADS Google Scholar

Lallemand, S., Heuret, A. & Boutelier, D. Über die Beziehungen zwischen Plattenneigung, Back-Arc-Spannung, absoluter Bewegung der oberen Platte und Krustennatur in Subduktionszonen. Geochem. Geophys. Geosyst.https://doi.org/10.1029/2005GC000917 (2005).

Artikel Google Scholar

Syracuse, EM & Abers, GA Globale Zusammenstellung von Variationen in der Plattentiefe unter Bogenvulkanen und deren Auswirkungen. Geochem. Geophys. Geosyst.https://doi.org/10.1029/2005GC001045 (2006).

Artikel Google Scholar

Schellart, WP, Freeman, J., Stegman, DR, Moresi, L. & May, D. Entwicklung und Vielfalt von Subduktionszonen, gesteuert durch die Plattenbreite. Nature (London) 446, 308–311 (2007).

Artikel ADS CAS PubMed Google Scholar

MacDougall, JG, Kincaid, C., Szwaja, S. & Fischer, KM Der Einfluss von Plattenneigungsvariationen, Lücken und Rollback auf die Mantelkeilströmung: Erkenntnisse aus Fluidexperimenten. Geophys. J. Int. 197, 705–730. https://doi.org/10.1093/gji/ggu053 (2014).

Artikel ADS Google Scholar

Hu, J. & Gurnis, M. Subduktionsdauer und Plattenabfall. Geochem. Geophys. Geosyst. 21, 1–34. https://doi.org/10.1029/2019GC008862 (2020).

Artikel Google Scholar

Schellart, WP Kontrolle des Alters und der Größe der Subduktionszone bei der Subduktion flacher Platten. Vorderseite. Erdwissenschaft. 8, 1–18. https://doi.org/10.3389/feart.2020.00026 (2020).

Artikel Google Scholar

Clennett, EJ et al. Eine quantitative tomotektonische Plattenrekonstruktion im westlichen Nordamerika und im östlichen Pazifikbecken. Geochem. Geophys. Geosyst. 21, 1–25. https://doi.org/10.1029/2020GC009117 (2020).

Artikel Google Scholar

Grose, CJ Eigenschaften der ozeanischen Lithosphäre: Überarbeitete Vorhersagen des Plattenkühlungsmodells. Planet Erde. Wissenschaft. Lette. 333–334, 250–264. https://doi.org/10.1016/j.epsl.2012.03.037 (2012).

Artikel ADS CAS Google Scholar

Richards, FD, Hoggard, MJ, Cowton, LR & White, NJ Neubewertung der thermischen Struktur der ozeanischen Lithosphäre mit überarbeiteten globalen Inventaren der Kellertiefen und Wärmeflussmessungen. J. Geophys. Res. Feste Erde 123, 9136–9161. https://doi.org/10.1029/2018JB015998 (2018).

Artikel ADS Google Scholar

Richards, F., Hoggard, M., Crosby, A., Ghelichkhan, S. & White, N. Struktur und Dynamik des ozeanischen Lithosphären-Asthenosphären-Systems. Physik. Planet Erde. Inter.https://doi.org/10.1016/j.pepi.2020.106559 (2020).

Artikel Google Scholar

Hayes, GP et al. Slab2, ein umfassendes Subduktionszonen-Geometriemodell. Wissenschaft 362, 58–61. https://doi.org/10.1126/science.aat4723 (2018).

Artikel ADS CAS PubMed Google Scholar

Giuliani, A. & Pearson, DG Kimberlite: Von der tiefen Erde zu Diamantenminen. Elemente 15, 377–380. https://doi.org/10.2138/GSELEMENTS.15.6.377 (2019).

Artikel CAS Google Scholar

Met Office. Cartopy: eine kartografische Python-Bibliothek mit einer Matplotlib-Schnittstelle. Exeter, Devon, https://doi.org/10.5281/zenodo.7430317 (2023).

Conrad, CP & Lithgow-Bertelloni, C. Die zeitliche Entwicklung der Plattenantriebskräfte: Bedeutung der Plattensaugkraft gegenüber der Plattenzugkraft während des Känozoikums. J. Geophys. Res. Feste Erde 109, 1–14. https://doi.org/10.1029/2004JB002991 (2004).

Artikel Google Scholar

Capitanio, FA, Morra, G. & Goes, S. Dynamik der Plattenbiegung am Graben und Platte-Platten-Kopplung. Geochem. Geophys. Geosyst.https://doi.org/10.1029/2008GC002348 (2009).

Artikel Google Scholar

Tetley, MG, Williams, SE, Gurnis, M., Flament, N. & Müller, RD Beschränkung absoluter Plattenbewegungen seit der Trias. J. Geophys. Res. Feste Erde 124, 7231–7258. https://doi.org/10.1029/2019JB017442 (2019).

Artikel Google Scholar

Müller, RD et al. Ein globales Plattenmodell, das die Deformation der Lithosphäre entlang großer Rifts und Orogene seit der Trias berücksichtigt. Tektonik 38, 1884–1907. https://doi.org/10.1029/2018TC005462 (2019).

Artikel ADS Google Scholar

Müller, RD et al. GPlates: Aufbau einer virtuellen Erde durch tiefe Zeit. Geochem. Geophys. Geosyst. 19, 2243–2261. https://doi.org/10.1029/2018GC007584 (2018).

Artikel ADS Google Scholar

Parsons, B. & Sclater, JG Eine Analyse der Variation der Bathymetrie des Meeresbodens und des Wärmeflusses mit dem Alter. J. Geophys. Res. 82, 803–827. https://doi.org/10.1029/JB082i005p00803 (1977).

Artikel ADS Google Scholar

Capitanio, FA, Stegman, DR, Moresi, LN & Sharples, W. Kontrolle der oberen Platte auf tiefe Subduktion, Grabenwanderungen und Deformationen an konvergenten Rändern. Tektonophysik 483, 80–92. https://doi.org/10.1016/j.tecto.2009.08.020 (2010).

Artikel ADS Google Scholar

Schepers, G. et al. Der Vorstoß der südamerikanischen Platte und der erzwungene Rückzug der Andengräben waren die Treiber für vorübergehende flache Subduktionsepisoden. Nat. Komm. 8, 1–9. https://doi.org/10.1038/ncomms15249 (2017).

Artikel CAS Google Scholar

Gonzalez, CM, Gorczyk, W. & Gerya, TV Dekarbonisierung subduzierender Platten: Erkenntnisse aus der petrologisch-thermomechanischen Modellierung. Gondwana Res. 36, 314–332. https://doi.org/10.1016/j.gr.2015.07.011 (2016).

Artikel ADS CAS Google Scholar

Keller, T., Katz, RF & Hirschmann, MM Flüchtige Stoffe unter mittelozeanischen Rücken: Tiefes Schmelzen, kanalisierter Transport, Fokussierung und Metasomatismus. Planet Erde. Wissenschaft. Lette. 464, 55–68. https://doi.org/10.1016/j.epsl.2017.02.006 (2017).

Artikel ADS CAS Google Scholar

Seton, M. et al. Ein globaler Datensatz des heutigen Krustenalters der Ozeane und der Ausbreitungsparameter des Meeresbodens. Geochem. Geophys. Geosyst. 21, e2020GC009214. https://doi.org/10.1029/2020GC009214 (2020).

Artikel ADS CAS Google Scholar

Dutkiewicz, A., Müller, RD, Cannon, J., Vaughan, S. & Zahirovic, S. Sequestrierung und Subduktion von Tiefseekarbonat im globalen Ozean seit der frühen Kreidezeit. Geologie 47, 91–94. https://doi.org/10.1130/G45424.1 (2019).

Artikel ADS CAS Google Scholar

Merdith, AS, Atkins, SE & Tetley, MG Tektonische Kontrollen der Kohlenstoff- und Serpentinitspeicherung in der subduzierten oberen ozeanischen Lithosphäre in den letzten 320 Ma. Vorderseite. Erdwissenschaft. 7, 1–23. https://doi.org/10.3389/feart.2019.00332 (2019).

Artikel ADS Google Scholar

Faccenda, M., Gerya, TV & Burlini, L. Tiefe Plattenhydratation, hervorgerufen durch biegebedingte Schwankungen des tektonischen Drucks. Nat. Geosci. 2, 790–793. https://doi.org/10.1038/ngeo656 (2009).

Artikel ADS CAS Google Scholar

Korenaga, J. Entstehung und Entwicklung der Plattentektonik auf der Erde: Theorien und Beobachtungen. Ann. Rev. Planet Erde. Wissenschaft. 41, 117–151. https://doi.org/10.1146/annurev-earth-050212-124208 (2013).

Artikel ADS CAS Google Scholar

Afonso, JC, Ranalli, G. & Fernàndez, M. Dichtestruktur und Auftrieb der ozeanischen Lithosphäre erneut untersucht. Geophys. Res. Lett.https://doi.org/10.1029/2007GL029515 (2007).

Artikel Google Scholar

Gurnis, M., Hall, C. & Lavier, L. Sich entwickelndes Kräftegleichgewicht während der beginnenden Subduktion. Geochem. Geophys. Geosyst.https://doi.org/10.1029/2003GC000681 (2004).

Artikel Google Scholar

Knesel, KM, Cohen, BE, Vasconcelos, PM & Thiede, DS Schnelle Änderung der Drift der australischen Platte dokumentiert Kollision mit dem Ontong-Java-Plateau. Natur 454, 754–757. https://doi.org/10.1038/nature07138 (2008).

Artikel ADS CAS PubMed Google Scholar

Liu, L. et al. Die Rolle der ozeanischen Plateau-Subduktion bei der Laramid-Orogenese. Nat. Geosci. 3, 353–357. https://doi.org/10.1038/ngeo829 (2010).

Artikel ADS CAS Google Scholar

Mao, W. & Zhong, S. Stagnation der Platte aufgrund einer Schicht mit verringerter Viskosität unter der Mantelübergangszone. Nat. Geosci. 11, 876–881. https://doi.org/10.1038/s41561-018-0225-2 (2018).

Artikel ADS CAS Google Scholar

King, SD, Frost, DJ & Rubie, DC Warum kalte Brammen in der Übergangszone stagnieren. Geologie 43, 231–234. https://doi.org/10.1130/G36320.1 (2015).

Artikel ADS Google Scholar

Chen, YW, Wu, J. & Suppe, J. Ausbreitung der Nazca-Subduktion entlang der Anden nach Süden. Natur 565, 441–447. https://doi.org/10.1038/s41586-018-0860-1 (2019).

Artikel ADS CAS PubMed Google Scholar

Jelsma, H., Barnett, W., Richards, S. & Lister, G. Tektonische Lage von Kimberliten. Lithos 112, 155–165. https://doi.org/10.1016/j.lithos.2009.06.030 (2009).

Artikel ADS CAS Google Scholar

Currie, CA & Beaumont, C. Stehen diamanthaltige Kimberlite aus der Kreidezeit mit der Subduktion bei geringem Winkel unter dem Westen Nordamerikas in Zusammenhang? Planet Erde. Wissenschaft. Lette. 303, 59–70. https://doi.org/10.1016/j.epsl.2010.12.036 (2011).

Artikel ADS CAS Google Scholar

Kjarsgaard, BA, Heaman, LM, Sarkar, C. & Pearson, DG Der mittelkreidezeitliche Kimberlitkorridor Nordamerikas: Nasses, kantengetriebenes Dekompressionsschmelzen einer tiefen Mantelquelle vom OIB-Typ. Geochem. Geophys. Geosyst. 18, 2727–2747. https://doi.org/10.1002/2016GC006761 (2017).

Artikel ADS CAS Google Scholar

Chen, Y. et al. Vereinbarkeit seismischer Strukturen und Kimberlit-Magmatismus aus der späten Kreidezeit im Norden von Alberta, Kanada. Geologie 48, 872–876. https://doi.org/10.1130/G47163.1 (2020).

Artikel ADS CAS Google Scholar

Mather, BR et al. Intraplattenvulkanismus, ausgelöst durch Ausbrüche im Plattenfluss. Wissenschaft. Adv. 6, 1–8. https://doi.org/10.1126/sciadv.abd0953 (2020).

Artikel Google Scholar

He, C. Abdrücke der subduzierten ozeanischen Paläo-Tethys-Lithosphäre auf Diskontinuitäten im oberen Erdmantel und die Bildung der großen magmatischen Provinz Emeishan. Geophys. J. Int. 231, 1298–1308. https://doi.org/10.1093/gji/ggac251 (2022).

Artikel ADS Google Scholar

Yang, J. & Faccenda, M. Intraplattenvulkanismus, der aus der aufsteigenden wasserhaltigen Mantelübergangszone stammt. Natur 579, 88–91. https://doi.org/10.1038/s41586-020-2045-y (2020).

Artikel ADS CAS PubMed Google Scholar

García-Reyes, A. & Dyment, J. Struktur, Alter und Ursprung der Karibischen Platte entschlüsselt. Planet Erde. Wissenschaft. Lette. 571, 117100. https://doi.org/10.1016/j.epsl.2021.117100 (2021).

Artikel CAS Google Scholar

Giuliani, A. et al. Kimberlit-Entstehung aus einer gemeinsamen karbonatreichen Primärschmelze, die durch Assimilation des lithosphärischen Mantels verändert wurde. Wissenschaft. Adv. 6, 1–10. https://doi.org/10.1126/sciadv.aaz0424 (2020).

Artikel CAS Google Scholar

Moresi, L. & Mather, B. Stripy. Ein Python-Modul zur (eingeschränkten) Triangulation in kartesischen Koordinaten und auf einer Kugel. Journal of Open Source Software 4, 1410. https://doi.org/10.21105/joss.01410 (2019).

Artikel ADS Google Scholar

Jackson, MG, Konter, JG & Becker, TW Ursprüngliches Helium, das von den heißesten Mantelwolken mitgerissen wird. Natur 542, 340–343. https://doi.org/10.1038/nature21023 (2017).

Artikel ADS CAS PubMed Google Scholar

Celli, NL, Lebedev, S., Schaeffer, AJ & Gaina, C. Afrikanische kratonische Lithosphäre, geschnitzt von Mantelwolken. Nat. Commun.https://doi.org/10.1038/s41467-019-13871-2 (2020).

Artikel PubMed PubMed Central Google Scholar

Regier, ME et al. Der in Superdeep-Diamanten aufgezeichnete Kohlenstoffkreislauf von der Lithosphäre bis zum unteren Erdmantel. Natur 585, 234–238. https://doi.org/10.1038/s41586-020-2676-z (2020).

Artikel CAS PubMed Google Scholar

Axen, GJ, van Wijk, JW & Currie, CA Basale Lithosphäre des kontinentalen Mantels, verschoben durch flache Plattensubduktion. Nat. Geosci. 11, 961–964. https://doi.org/10.1038/s41561-018-0263-9 (2018).

Artikel ADS CAS Google Scholar

Currie, CA & Copeland, P. Numerische Modelle der Farallon-Plattensubduktion: Erstellen und Entfernen einer flachen Platte. Geosphäre 18, 476–502. https://doi.org/10.1130/GES02393.1 (2022).

Artikel Google Scholar

Humphreys, E. D. Post-Laramide removal of the Farallon slab, western United States. Geology 23, 987. 2.3.CO;2">https://doi.org/10.1130/0091-7613(1995)023<0987:PLROTF>2.3.CO;2 (1995).

2.3.CO;2" data-track-action="article reference" href="https://doi.org/10.1130%2F0091-7613%281995%29023%3C0987%3APLROTF%3E2.3.CO%3B2" aria-label="Article reference 59" data-doi="10.1130/0091-7613(1995)0232.3.CO;2">Artikel ADS Google Scholar

Referenzen herunterladen

Diese Studie wurde vom AuScope Simulation, Analysis & Modeling-Knoten unterstützt, der von der australischen Regierung im Rahmen der National Collaborative Research Infrastructure Strategy, NCRIS, finanziert wurde. Wir danken dem Australian Research Council für die Finanzierung durch den Zuschuss DP200100966 (MS).

EarthByte Group, School of Geosciences, The University of Sydney, Sydney, 2006, Australien

Ben R. Mather, R. Dietmar Müller, Christopher P. Alfonso, Maria Seton und Nicky M. Wright

Sie können diesen Autor auch in PubMed Google Scholar suchen

Sie können diesen Autor auch in PubMed Google Scholar suchen

Sie können diesen Autor auch in PubMed Google Scholar suchen

Sie können diesen Autor auch in PubMed Google Scholar suchen

Sie können diesen Autor auch in PubMed Google Scholar suchen

BM konzipierte die Idee, entwickelte die Python-Tools, generierte die Abbildungen und schrieb das Manuskript. RDM konzipierte die Idee und Anwendungen der Slab-Dip-Formulierung zur Erklärung von Kimberlit-Eruptionen. CPA verbesserte die tektonische Rekonstruktion in Nordamerika und stellte Altersraster und Raster zur Ausbreitungsrate des Meeresbodens bereit. MS konzipierte die Experimente und verbesserte das Manuskript. NMW lieferte die Rekonstruktion großer magmatischer Provinzen.

Korrespondenz mit Ben R. Mather.

Die Autoren geben an, dass keine Interessenkonflikte bestehen

Springer Nature bleibt neutral hinsichtlich der Zuständigkeitsansprüche in veröffentlichten Karten und institutionellen Zugehörigkeiten.

Ergänzende Informationen 1.

Open Access Dieser Artikel ist unter einer Creative Commons Attribution 4.0 International License lizenziert, die die Nutzung, Weitergabe, Anpassung, Verbreitung und Reproduktion in jedem Medium oder Format erlaubt, sofern Sie den/die Originalautor(en) und die Quelle angemessen angeben. Geben Sie einen Link zur Creative Commons-Lizenz an und geben Sie an, ob Änderungen vorgenommen wurden. Die Bilder oder anderes Material Dritter in diesem Artikel sind in der Creative Commons-Lizenz des Artikels enthalten, sofern in der Quellenangabe für das Material nichts anderes angegeben ist. Wenn Material nicht in der Creative-Commons-Lizenz des Artikels enthalten ist und Ihre beabsichtigte Nutzung nicht gesetzlich zulässig ist oder über die zulässige Nutzung hinausgeht, müssen Sie die Genehmigung direkt vom Urheberrechtsinhaber einholen. Um eine Kopie dieser Lizenz anzuzeigen, besuchen Sie http://creativecommons.org/licenses/by/4.0/.

Nachdrucke und Genehmigungen

Mather, BR, Müller, RD, Alfonso, CP et al. Kimberlit-Eruptionen, angetrieben durch Plattenfluss und Subduktionswinkel. Sci Rep 13, 9216 (2023). https://doi.org/10.1038/s41598-023-36250-w

Zitat herunterladen

Eingegangen: 25. Oktober 2022

Angenommen: 30. Mai 2023

Veröffentlicht: 06. Juni 2023

DOI: https://doi.org/10.1038/s41598-023-36250-w

Jeder, mit dem Sie den folgenden Link teilen, kann diesen Inhalt lesen:

Leider ist für diesen Artikel derzeit kein Link zum Teilen verfügbar.

Bereitgestellt von der Content-Sharing-Initiative Springer Nature SharedIt

Durch das Absenden eines Kommentars erklären Sie sich damit einverstanden, unsere Nutzungsbedingungen und Community-Richtlinien einzuhalten. Wenn Sie etwas als missbräuchlich empfinden oder etwas nicht unseren Bedingungen oder Richtlinien entspricht, kennzeichnen Sie es bitte als unangemessen.