Lithogéochimie des unités géologiques de la région du lac Serpentine
Les tableaux ci-dessous résument les caractéristiques lithogéochimiques des unités géologiques de la région du lac Serpentine. Ces unités sont décrites dans le Bulletin géologiQUE couvrant ce territoire et dans le Lexique stratigraphique du Québec. Les 138 analyses utilisées ici proviennent d’échantillons collectés lors de la campagne de cartographie de 2022. Elles ont été sélectionnées en fonction des critères suivants : une somme des oxydes majeurs comprise entre 97,5 % et 101,5 %, une perte au feu (LOI) inférieure à 3 % (à l’exception des roches ultramafiques, des roches calcosilicatées et de plusieurs roches volcaniques). Une sélection préalable des analyses avec l’indice d’altération d’Ishikawa a été faite pour séparer les analyses de roches altérées des roches non altérées (Large et al. 2001). Ces analyses ont été réalisées par le laboratoire d’Actlabs d’Ancaster, en Ontario.
Les analyses ont été soumises à un processus d’assurance et de contrôle de la qualité interne et en laboratoire. Ainsi, pour s’assurer de la justesse et de la précision des valeurs fournies par le laboratoire, la Direction de l’acquisition des connaissances géoscientifiques (DAGQ) insère régulièrement des blancs, des standards et des duplicatas. Les matériaux de référence représentent ~10 % des analyses.
La majorité des échantillons de la base de données ont été analysés pour les oxydes majeurs, les éléments en traces et les métaux. Les analyses ont été effectuées par différentes techniques en fonction des éléments, telles que la spectrométrie de masse à plasma à couplage inductif (ICP-MS), la spectrométrie d’émission optique à plasma à couplage inductif (ICP-AES) et l’activation neutronique (INAA). Pour plus de renseignements sur les techniques d’analyse et de dissolution utilisées, se référer à l’information disponible pour chaque échantillon dans SIGÉOM à la carte.
Lorsque les unités géologiques comprennent plus de 10 analyses, les profils des éléments de terres rares et multiéléments sont regroupés pour constituer des enveloppes comprenant les 25e et 75e percentiles de la population. Cette procédure a été retenue pour simplifier la visualisation d’un grand nombre de profils. Les graphiques affichent les enveloppes ainsi que la ligne médiane.
Les éléments de terres rares sont normalisés d’après les valeurs de Palme et O’Neill (2004). Les teneurs anomales, distinctives ou jugées importantes sont inscrites en caractère gras dans les tableaux.
Roches supracrustales
Unité stratigraphique ou lithologique |
Classification |
Affinité |
Environnement tectonique |
Nbre Mg |
Terres rares |
Diagramme multiélément |
Remarques |
|
---|---|---|---|---|---|---|---|---|
Roches volcaniques et lithologies associées |
||||||||
Groupe de Parent (pPpa9a)
|
Basalte subalcalin
|
Tholéiitique (Diagrammes)
|
Source : N-MORB (Diagrammes)
|
62,02 à 64,53 |
0,39 < (La/Yb)N < 0,74 0,52 < (La/Sm)N < 0,68 0,88 < (Gd/Yb)N < 1,06 0,98 < Eu/Eu* < 1,03
|
Profils montrant une pente positive; Anomalies négatives : Th et Nb (fortes); Zr, Ti et Y (faibles); Anomalies positives : Sm et Eu (faibles)
|
Enrichissement en terres rares lourdes par rapport aux terres rares légères. Pente positive. Anomalie négative en Th et Nb pouvant résulter d’une contamination crustale | |
Groupe de Parent (pPpa9b) |
Basalte |
Tholéiitique à transitionnel |
Source : E-MORB |
48,94 à 66,72 |
1,91 < (La/Yb)N < 3,32 1,11 < (La/Sm)N < 1,54 1,45 < (Gd/Yb)N < 1,63 0,96 < Eu/Eu* < 1,10 |
Profil montrant une légère pente négative; Anomalies négatives : Sm et Y (faibles) |
Pente légèrement négative. Enrichissement en terres rares légères par rapport aux terres rares lourdes |
|
Groupe de Parent (pPpa9c) |
Basalte subalcalin, andésite basaltique et andésite |
Transitionnel à calco-alcalin |
Environnement de formation variable; principalement entre les basaltes d’arc insulaire (IAB) et les E-MORB |
41,30 à 56,12 |
2,02 < (La/Yb)N < 4,72 1,11 < (La/Sm)N < 2,11 1,30 < (Gd/Yb)N < 2,02 0,94 < Eu/Eu* < 1,13 |
Pente montrant une pente négative; Anomalies négatives : Nb, Sm et Ti (faibles); Anomalies positives : La et Zr (faibles) |
Pente négative montrant un enrichissement en terres rares légères par rapport aux terres rares lourdes | |
Groupe de Parent (pPpa10a) |
Basalte subalcalin |
Tholéiitique |
N-MORB à E-MORB |
55,51 à 64,06 |
1,50 < (La/Yb)N < 1,62 1,04 < (La/Sm)N < 1,07 1,23 < (Gd/Yb)N < 1,52 1,03 < Eu/Eu* < 1,09 |
Profil montrant une faible pente négative et une légère contamination crustale montrée par une anomalie négative en Nb et Ta; Anomalies négatives : Th et Nb (fortes), Zr, Hf, Gd et Y (faibles); Anomalies positives : Eu et Ti (faibles) |
Profil de terres rares présentant un profil assez plat avec faible pente négative vers les terres rares lourdes | |
Groupe de Parent (pPpa10b) |
Basalte et basalte subalcalin |
Tholéiitique à transitionnel |
N-MORB à E-MORB |
54,86 à 63,34 |
1,85 < (La/Yb)N < 2,13 1,13 < (La/Sm)N < 1,27 1,25 < (Gd/Yb)N < 1,48 0,98 < Eu/Eu* < 1,11 |
Profil montrant une pente négative et une légère contamination crustale montrée par une anomalie négative en Nb et Ta; Anomalies négatives : Nb, Zr et Hf (fortes); Ti et Y (faibles); Anomalies positives : La, Eu, Ti et Er (faibles) |
Profil de terres rares présentant une faible pente négative vers les terres rares lourdes | |
Groupe de Parent (pPpa10c) |
Basalte, andésite basaltique et andésite |
Transitionnel à calco-alcalin |
Environnement de formation variable; principalement entre les basaltes d’arc insulaire (IAB) et les E-MORB |
34,62 à 53,94 |
2,62 < (La/Yb)N < 5,23 1,55 < (La/Sm)N < 2,70 1,21 < (Gd/Yb)N < 1,75 0.92 < Eu/Eu* < 1.06 |
Profil montrant une pente négative et une possible contamination crustale montrée par une anomalie négative en Nb et Ta; Anomalies négatives : Nb, Ta et Ti (fortes); Sm et Y (faibles); Anomalies positives : Th et La (fortes); Zr (modérée); Er (faible) |
Profil de terres rares présentant une pente négative modérée avec un enrichissement en terres rares légères | |
Groupe de Parent (pPpa12) |
Basalte subalcalin |
Indéterminée |
Environnement de formation variable; principalement entre les basaltes d’arc insulaire (IAB) et les E-MORB |
40,85 à 64,52 |
1,22 < (La/Yb)N < 1,25 1,11 < (La/Sm)N < 2,03 1,06 < (Gd/Yb)N < 1,63 0,92 < Eu/Eu* < 1,25 |
Profil montrant une pente négative et une possible contamination crustale montrée par une anomalie négative en Nb et Ta; Anomalies négatives : Nb et Ta (fortes); Zr, Hf et Y (faibles); Anomalies positives : Eu et Er (faibles) |
Patron de terres rares présentant une pente négative modérée montrée par un enrichissement en terres rares légères par rapport aux terres rares lourdes | |
Suite de Qummitaliup (pPqmp) |
Gabbro alcalin, gabbro et gabbronorite |
Principalement tholéiitique à transitionnel |
N-MORB à E-MORB |
41,90 à à 76,24 |
1,14 < (La/Yb)N < 5,72 0,78 < (La/Sm)N < 2,16 0,99 < (Gd/Yb)N < 1,82 0,89 < Eu/Eu* < 1,11 |
Profil montrant une pente modérément négative; Anomalies négatives : Nb et Th (fortes), Ta, Zr, Y et Yb (faibles); Anomalies positives : Nd, Eu et Er (faibles) |
Profil montrant une pente modérément négative avec un enrichissement plus important en terres rares légères par rapport aux terres rares lourdes | |
Roches sédimentaires |
||||||||
Groupe de Spartan (pPsp1)
|
Wacke, shale et shale ferrifère (Diagramme A)
|
|
14,24 à 51,74 |
4,60 < (La/Yb)N < 54,27 3,17 < (La/Sm)N < 5,61 1,02 < (Gd/Yb)N < 4,78 0,54 < Eu/Eu* < 0,88 |
Profil montrant une forte pente négative et une possible contamination crustale montrée par une forte anomalie négative en Nb et Ta; Anomalies négatives : Nb, Ta et Ti (fortes); Sm et Eu (modérées); Y (faible) Anomalies positives : Zr, Hf et Er (faibles) |
Patron des terres rares présentant une pente forte avec un fort enrichissement en terres rares légères et une faible anomalie négative en Eu. | ||
Complexe de Qaaneq (pPqaa1) |
Shale et shale ferrifére (Diagramme A) |
39,58 à 50,64 |
4,03 < (La/Yb)N < 13,41 3,28 < (La/Sm)N < 4,41 0,60 < (Gd/Yb)N < 1,88 0,69 < Eu/Eu* < 0.93 |
Profil montrant une forte pente négative et une possible contamination crustale montrée par une forte anomalie négative en Nb et Ta; Anomalies négatives : Nb, Ta et Ti (fortes); Eu et Sm (modérée); Y et Yb (faibles); Anomalies positives : Zr, Hf et Er (faibles) |
Patron de terres rares présentant une forte pente négative avec un fort enrichissement en terres rares légères. | |||
Complexe de Qaaneq (pPqaa2) |
Shale ferrifère à grès ferrifère (Diagramme A) |
67,25 à 77,40 |
1,37 < (La/Yb)N < 6,31 0,95 < (La/Sm)N < 2,14 1,23 < (Gd/Yb)N < 1,94 0,66 < Eu/Eu* < 1,74 |
Profil montrant une pente négative et contamination crustale suggérée par une forte anomalie négative en Nb et Ta; Anomalies négatives : Nb et Ta (fortes); Th et Ti (modérée); Zr, Hf et Y (faibles); Anomalies positives : Eu (modérée); Gd et Er (faible) |
Patron de terres rares présentant une pente négative modérée |
Roches intrusives
Unité stratigraphique ou lithologique |
Classification |
Affinité |
Environnement tectonique |
Nbre Mg |
Terres rares |
Diagramme multiélément |
Remarques |
---|---|---|---|---|---|---|---|
Roches intrusives felsiques à intermédiaires |
|||||||
Suite de Sanimuapik (pPsak1)
|
Diorite litée et granodiorite
|
Granitoïde de type I-S; métalumineux à faiblement peralumineux (Diagramme A et |
Granite d’arc volcanique
|
44,62 à 58,50 |
7,27 < (La/Yb)N < 7,72 2,75 < (La/Sm)N < 3,37 1,42 < (Gd/Yb)N < 1,87 0,79 < Eu/Eu* < 1,06
|
Profil montrant une forte pente négative et une possible contamination crustale montrée par une forte anomalie négative en Nb et Ta; Anomalies négatives : Nb, Sm et Ti (fortes à modérées); Ta et Y (faibles);
Anomalies positives : Zr et Hf (modérées); Nd, Er, Yb et Lu (faibles)
|
Profil montrant une forte pente négative attribuable à un enrichissement en terres rares légères et une faible anomalie négative en Eu |
Suite de Sanimuapik (pPsak2)
|
Diorite litée
|
Granitoïde de type I-S; métalumineux (Diagramme A et
|
Granite d’arc volcanique
|
49,86 |
(La/Yb)N = 9,42 (La/Sm)N = 2,75 (Gd/Yb)N = 2,00 Eu/Eu* = 0,97
|
Profil montrant une forte pente négative et une possible contamination crustale montrée par une forte anomalie négative en Nb et Ta;
Anomalies négatives : Th, Nb, Ta et Ti (fortes); Sm et Y (faibles); Anomalies positives : Zr et Er (faibles)
|
Forte pente négative montrée par un enrichissement en terres rares lourdes. Il n’y a qu’une seule analyse pour cette unité.
|
Suite de Niqituraaqiaruk (pPnqk1) |
Diorite, monzodiorite, gabbro alcalin et gabbro |
Granitoïde de type I-S; métalumineux (Diagramme A et |
Granite d’arc volcanique |
57,28 à 75,51 |
6,41 < (La/Yb)N < 28,44 2,13 < (La/Sm)N < 4,49 1,84 < (Gd/Yb)N < 2,85 0,84 < Eu/Eu* < 0,96 |
Profil montrant une forte pente négative et une possible contamination crustale montrée par une forte anomalie négative en Nb et Ta; Anomalies négatives : Nb, Ta, Zr et Ti (fortes); Th et Hf (modérées); Y et Yb (faibles); Anomalies positives Eu et Er (faibles) |
Forte pente négative montrée par un enrichissement plus élevé en terres rares légères qu’en terres rares lourdes |
Unité stratigraphique ou lithologique |
Classification |
Affinité |
Environnement tectonique |
Nbre Mg |
Terres rares |
Diagramme multiélément |
Remarques |
---|---|---|---|---|---|---|---|
Roches intrusives mafiques-ultramafiques et lithologies associées |
|||||||
Essaim de Franklin (nPfra)
|
Gabbro
|
Transitionnel à calco-alcalin
|
E-MORB montrant de la contamination crustale |
66.80 |
(La/Yb)N = 4,31 (La/Sm)N = 2,07 (Gd/Yb)N = 1,61 Eu/Eu* = 1,04
|
Profil montrant une pente négative modérée et une possible contamination crustale montrée par une faible anomalie négative en Nb et Ta; Anomalies négatives : Th, Nb, Ta et Sm (modérées à faibles)
|
Il n’y a qu’une seule analyse, pour cette unité. Pente négative modérée
|
Suite de Vanasse (pPvas1) |
Gabbro à olivine, gabbro, gabbronorite et diorite. |
Tholéiitique à transitionnel. |
Environnement de formation variable; E-MORB à N-MORB ainsi que basaltes de rifts continentaux (CRB) et basaltes d’île océanique (OIB) (Diagrammes A et B) | 39.53 à 84,15 |
1,46 < (La/Yb)N < 7,76 0,76 < (La/Sm)N < 2,61 1,30 < (Gd/Yb)N < 2,56 0,69 < Eu/Eu* < 1,38 |
Profil montrant une pente négative et une possible contamination crustale montrée par une forte anomalie négative en Nb et Ta; Anomalies négatives : Th, Nb, Ta, Zr, Hf et Ti (fortes à modérées); Sm, Y et Yb (faibles); Anomalies positives : Eu (modérée); Er (faible) |
Forte pente négative montrée par un enrichissement plus élevé en terres rares légères qu’en terres rares lourdes |
Suite de Vanasse (pPvas2) |
Gabbro alcalin, monzogabbro et monzodiorite |
Très variable |
Environnement de formation variable; E-MORB à N-MORB, ainsi que basaltes de rifts continentaux (CRB) et basaltes d’île océanique (OIB) (Diagrammes A et B) |
42,23 à 69,54 |
2,08 < (La/Yb)N < 8,94 1,01 < (La/Sm)N < 2,74 1,58 < (Gd/Yb)N < 2,42 0,91 < Eu/Eu* < 1,40 |
Profil montrant une pente négative et une possible contamination crustale montrée par une forte anomalie négative en Nb et Ta; Anomalies négatives : Th, Nb, Ta et Ti (fortes à modérées); Sm, Y et Yb (faibles); Anomalies positives : Eu (modérée); Er (faible) |
Forte pente négative montrée par un enrichissement plus élevé en terres rares légères qu’en terres rares lourdes |
Suite de Sanimualuk (pPsnm1) |
Gabbro et roche ultramafique |
83,71 à 91,67 |
0,32 < (La/Yb)N < 0,59 0,31 < (La/Sm)N < 0,63 1,01 < (Gd/Yb)N < 1,15 0,98 < Eu/Eu* < 1,08 |
Pente légèrement positive devant rapidement plat vers les terres rares lourdes; Anomalies négatives : Zr (modérée); Ti, Y, Yb et Lu (faibles); Anomalies positives : Sm, Dy et Er (faibles) |
Pente faiblement positive avec un enrichissement plus important en terres rares lourdes qu’en terres rares légères | ||
Suite de Sanimualuk (pPsnm2) |
Gabbro à gabbronorite |
Tholéiitique |
N-MORB |
59,29 à 71,40 |
0,30 < (La/Yb)N < 1,07 0,40 < (La/Sm)N < 0,97 0,87 < (Gd/Yb)N < 1,23 0,91 < Eu/Eu* < 2,33 |
Pente légèrement positive devant rapidement plat vers les terres rares lourdes; Anomalies négatives : Nb, Zr, Ti et Y (faibles); Anomalies positives : Ta et Sm |
Pente faiblement positive avec un enrichissement plus important en terres rares lourdes qu’en terres rares légères |
Suite d’Illuinaqtuut (pPiqt1) |
Gabbro à gabbronorite |
Tholéiitique |
En dehors des MORBs, montre de la contamination crustale |
56,79 à 79,44 |
2,64 < (La/Yb)N < 2,78 1,11 < (La/Sm)N < 1,28 1,47 < (Gd/Yb)N < 1,89 1,41 < Eu/Eu* < 1,56 |
Profil montrant une faible pente négative ainsi qu’une anomalie négative en Nb, Ta qui pourrait témoigner d’une contamination crustale. Anomalies négatives : Th, Nb, Ta, Zr et Hf (fortes); Ti (modérée); Y (faible); Anomalie positive : Eu (modérée) |
Profil montrant une pente faiblement négative avec une anomalie positive en Eu |
Suite d’Illuinaqtuut (pPiqt1b) |
Clinopyroxénite |
Tholéiitique |
MORB |
69,71 à 88,91 |
1,05 < (La/Yb)N < 1,51 0,38 < (La/Sm)N < 0,54 1,64 < (Gd/Yb)N < 2,34 0,74 < Eu/Eu* < 1,05 |
Profil montrant une pente faiblement à modérément négative ainsi qu’une anomalie négative en Nb, Ta qui pourrait témoigner d’une contamination crustale; Anomalies négatives : Th, Ta, Zr et Hf (fortes); La, Ce et Ti (modérées); Pr et Y (faibles); Anomalie positive : Eu (forte) |
Pente fortement positive montrant un enrichissement en terres rares légères et pente modérément négative montrant un appauvrissement en terres rares lourdes |
Suite d’Illuinaqtuut (pPiqt2) |
Gabbro à olivine, clinopyroxénite et gabbro foïdique. |
Transitionnel |
Environnement de formation variable; E-MORB à N-MORB, ainsi que basaltes de rifts continentaux (CRB) et basaltes d’île océanique (OIB) (Diagramme) |
62,30 à 67,08 |
3,02 < (La/Yb)N < 10,26 0,97 < (La/Sm)N < 1,28 2,03 < (Gd/Yb)N < 4,37 0,90 < Eu/Eu* < 1,10 |
Profil montrant une forte pente négative ainsi qu’une anomalie négative en Nb, Ta qui pourrait témoigner d’une contamination crustale; Anomalies négatives : Ta, Zr, Hf et Ti (fortes); Ti, Th et Nb (modérées); Y (faible); Anomalie positive : Eu (forte) |
Forte pente négative montrée par un enrichissement plus élevé en terres rares légères qu’en terres rares lourdes |
Suite d’Illuinaqtuut (pPiqt3) |
Clinopyroxénite, gabbro à olivine |
Tholéiitique à transitionnel |
Environnement de formation variable; E-MORB ainsi que basaltes de rifts continentaux (CRB) et basaltes d’île océanique (OIB) |
47,35 à 62,31 |
2,06 < (La/Yb)N < 3,55 1,04 < (La/Sm)N < 1,51 1,13 < (Gd/Yb)N < 1,95 0,97 < Eu/Eu* < 1,24 |
Profil montrant une faible pente à faiblement négative; Anomalies négatives : Th, Nb, Zr et Hf (modérées à faibles); Anomalies positives : Ta et Ti (faibles) |
Profil montrant une pente faiblement négative |
RÉFÉRENCES
Autres publications
AGRAWAL, S., 2008. Tectonic discrimination of basic and ultrabasic volcanic rocks through log-transformed ratios of immobile trace elements. International Geology Review; volume 50, pages 1057-1079. doi.org/10.2747/0020-6814.50.12.1057
BARTON, M.D., YOUNG, S., 2002. Non-pegmatitic deposits of beryllium: mineralogy, geology, phase equilibria and origin In: Beryllium: Mineralogy, Petrology and Geochemistry (Grew, E.S., editor). Reviews in Mineralogy and Geochemistry; volume 50, pages 591-691.
CABANIS, B., LECOLLE, M., 1989. The La/10-Y/15-Nb/8 diagram: A tool for distinguishing volcanic series and discovering crustal mixing and/or contamination: Comptes Rendus de l’Académie des Sciences, Série 2, Sciences de la Terre; volume 309, 20 pages.
DE LA ROCHE, H., LETERRIER, J., GRANDCLAUDE, P., MARCHAL, M., 1980. A classification of volcanic and plutonic rocks using R1-R2 diagrams and major element analyses – its relationships with current nomenclature. Chemical Geology; volume 29, pages 183-210. doi.org/10.1016/0009-2541(80)90020-0
FROST, B.R., FROST, C.D., 2008. A geochemical classification for feldspathic igneous rocks. Journal of Petrology; volume 49, pages 1955-1969. doi.org/10.1093/petrology/egn054
FROST, B.R., BARNES, C.G., COLLINS, W.J., ARCULUS, R.J., ELLIS, D.J., FROST, C.D., 2001. A geochemical classification for granitic rocks. Journal of Petrology; volume 12, pages 2033- 2048. doi.org/10.1093/petrology/42.11.2033
HERRON, M.H., 1988. Geochemical classification of terrigenous sands and shales from core or log data. Journal of Sedimentary Research; volume 58, pages 820-829. doi.org/10.1306/212F8E77-2B24-11D7-8648000102C1865D
MANIAR, P.D., PICCOLI, P.M., 1989. Tectonic discrimination of granitoids. Geological Society of America Bulletin; volume 101, pages 635-643. doi.org/10.1130/0016-7606(1989)101<0635:TDOG>2.3.CO;2
MCDONOUGH, W.F., SUN, S.S., 1995. The composition of the Earth. Chemical Geology; volume 120, pages 223-253. doi.org/10.1016/0009-2541(94)00140-4
MIDDLEMOST, A.K., 1994. Naming materials in the magma/igneous rock system. Earth-Science Reviews; volume 37, pages 215-224. doi.org/10.1016/0012-8252(94)90029-9
PALME, H., O’NEILL, H.S.C., 2004. Cosmochemical estimates of mantle composition. In: Treatise on Geochemistry. (Holland, H.D., Turekian, K.K. editors), Elsevier, Amsterdam, The Netherlands; volume 2, pages 1-38. doi.org/10.1016/B978-0-08-095975-7.00201-1
PEARCE, J.A., 2008. Geochemical fingerprinting of oceanic basalts with applications to ophiolite classification and the search for Archean oceanic crust, Lithos; volume 100, pages 14-48. doi.org/10.1016/j.lithos.2007.06.016
PEARCE, J.A., 2014. Immobile Element Fingerprinting of Ophiolites, Elements; volume 10, pages 101-108.
PEARCE, J.A., GALE, G.H., 1977. Identification of ore-deposition environment from trace element geochemistry of associated igneous host rocks. Geological Society, London; Special Publications, volume 7, pages 14-24. doi.org/10.1144/GSL.SP.1977.007.01.03.
PEARCE, J.A., HARRIS, B.W., TINDLE, A.G., 1984. Trace element discrimination diagram for tectonic interpretation of granitic rocks. Journal of Petrology; volume 25, pages 956-983. doi.org/10.1093/petrology/25.4.956
ROSS, P.-S., BÉDARD, J.H., 2009. Magmatic affinity of modern and ancient subalkaline volcanic rocks determined from trace-element discriminant diagrams. Canadian Journal of Earth Science; volume 46, pages 823-839. doi.org/10.1139/E09-054
WHALEN, J.B., CURRIE, K.L., CHAPPELL, B.W., 1987. A-Type granites: Geochemical characteristics, discrimination and petrogenesis. Contributions to Mineralogy and Petrology; volume 95, pages 407-419. doi.org/10.1007/BF00402202
WINCHESTER, J.A., FLOYD, P.A., 1977. Geochemical discrimination of different magma series and their differentiation products using immobile elements. Chemical Geology; volume 20, pages 325-343. doi.org/10.1016/0009-2541(77)90057-2
WOOD, D.A., 1980. The application of a Th–Hf–Ta diagram to problems of tectonomagmatic classification and to establishing the nature of crustal contamination of basaltic lavas of the British Tertiary volcanic province. Earth and Planetary Science Letters; volume 50, pages 11-30. doi.org/10.1016/0012-821X(80)90116-8