1. Introduction
The Eurasian–Greenland continental break-up during the Paleocene–Eocene transition marked the culmination of a continental rift system and the development of Volcanic Rifted Margins (Saunders et al. Reference Saunders, Fitton, Kerr, Norry, Kent, Mahoney and Coffin1997; Meyer, van Wijk & Gernigon, Reference Meyer, Van Wijk, Gernigon, Foulger and Jurdy2007). Up to the end of the Danian stage, the early geodynamic history of this continental break-up is marked by widespread extensional basin formation and localized minor magmatism. A first peak of volcanic activity around 60 Ma caused continental flood basalt eruptions and formed numerous large central volcanoes whose eroded remnants now crop out as igneous centres in West and East Greenland, NE Ireland and NW Scotland (Emeleus & Bell, Reference Emeleus and Bell2005). A second ‘pulse’ of high magma eruption rates at c. 55 Ma coincided with the final break-up of the continents and the initiation of seafloor spreading, and emplaced the volcanic seaward-dipping wedges (Saunders et al. Reference Saunders, Fitton, Kerr, Norry, Kent, Mahoney and Coffin1997). The Volcanic Rifted Margins igneous formations of the North Atlantic Igneous Province (NAIP) extend over an area of 1.3 × 106 km2 and cover nearly the whole length of the conjugate NE Atlantic margins (Saunders et al. Reference Saunders, Fitton, Kerr, Norry, Kent, Mahoney and Coffin1997; Eldholm et al. Reference Eldholm, Gladczenko, Skogseid, Planke and Nøttvedt2000).
Different Volcanic Rifted Margins structural zones characterize distinct time–space episodes of the rift-to-drift evolution history (Fig. 1) (Eldholm et al. Reference Eldholm, Gladczenko, Skogseid, Planke and Nøttvedt2000; Menzies et al. Reference Menzies, Klemperer, Ebinger, Baker, Menzies, Klemperer, Ebinger and Baker2002; Coffin & Eldholm, Reference Coffin, Eldholm, Selley, Cocks and Plimer2005): (a) continental flood basalt areas mark the late stage of a continental rift; (b) the seaward-dipping wedges are extrusive rock units from the rift-to-drift transition; (c) High Velocity Lower Crustal bodies presumably reflect massive intrusive complexes near the Continent Ocean Boundary. Unfortunately, even with state-of-the-art drilling techniques, the High Velocity Lower Crustal bodies will remain an unreachable drilling target for the time being. The specific structure of these tectono-magmatic crustal zones depends on the interaction of lithospheric and asthenospheric properties and dynamics of rifting and magma supply during the formation of a Volcanic Rifted Margins system.
Figure 1. Schematic tectono-magmatic division of a Volcanic Rifted Margin (VRM) system. VRMs characteristically comprise three zones: Continental Flood Basalt (CFB) areas, Seaward-Dipping Reflector Sequences (SDRS) and High Velocity Lower Crustal bodies (HVLC) next to the Continent Ocean Boundary (COB). UCC, MCC, LCC – Upper, Middle and Lower Continental Crust (modified from Eldholm et al. Reference Eldholm, Gladczenko, Skogseid, Planke and Nøttvedt2000, and Menzies et al. Reference Menzies, Klemperer, Ebinger, Baker, Menzies, Klemperer, Ebinger and Baker2002).
The Paleocene magmatism of the Isle of Rum in NW Scotland forms an integral part of the pre-break-up NAIP continental-flood-basalt tectono-magmatic zone. The igneous rocks of Rum comprise early felsic magmas, that are cross-cut by later plutonic ultrabasic to basic Layered Suite (Emeleus, Reference Emeleus1997). The felsic magmas are considered to originate from melting of gneissic crust (Troll, Donaldson & Emeleus, Reference Troll, Donaldson and Emeleus2004) by ascending voluminous basic mantle melts. Such associations of ‘early felsic–later mafic’ rocks have also been recovered in cores of ocean drilling campaigns at the North Atlantic seaward-dipping wedges, for example, at the Vøring Plateau (Norwegian Sea, Hole 642E, ODP Leg 104: Eldholm, Thiede & Taylor, Reference Eldholm, Thiede and Taylor1987). In Hole 642E, crustal melts (the ‘Lower Series’, 133 m drilled, total thickness not known) are covered by a 770 m thick sequence of tholeiitic basalts (‘Upper Series’). Rum and the Vøring Plateau, however, represent different zonal and temporal evolutionary stages of formation of Volcanic Rifted Margins, and a comparison of crustal melting in these two locations can contribute to a better understanding of the interactions between ascending mantle melts and crust during continental break-up. Geochemical studies of crustal melts formed during continental break-up have other merits as well. The radiogenic isotopic compositions and trace element patterns of felsic crustal melts and mixtures of basic and felsic magmas are potentially very useful to identify the crustal components involved in melting and mixing, which in turn may provide information on the depth of melting and the make-up of inaccessible crustal layers at the regional or local scale. This paper focuses on a study of the samples from the Isle of Rum. Detailed results of samples from Vøring Plateau, ODP Leg 104, Hole 642E are reported elsewhere (Meyer et al. Reference Meyer, Hertogen, Pedersen, Viereck-Götte and Abratis2009).
2. Geological setting and previous work
The varying igneous rock types observed and described from the Isle of Rum have greatly contributed to the development of modern views of the magmatic, as well as the temporal and structural evolution of the British Palaeogene Igneous Province (Bailey, Reference Bailey1945; Emeleus, Reference Emeleus1997; Emeleus & Bell, Reference Emeleus and Bell2005; Troll, Donaldson & Emeleus, Reference Troll, Donaldson and Emeleus2004; Troll et al. Reference Troll, Nicoll, Emeleus and Donaldson2008). Rum is located on the Hebridean Terrane in NW Scotland (Fig. 2). The igneous centre of Rum was emplaced at around 60.5 Ma (Hamilton et al. Reference Hamilton, Pearson, Thompson, Kelley and Emeleus1998), most likely in a relatively short period of time (about 500 000 years) (Troll et al. Reference Troll, Nicoll, Emeleus and Donaldson2008). An elliptical ring fault, 12 km across, bounds the igneous centre and marks the remnants of a caldera rim. Major volcanic activity on the Isle of Rum started with the eruption of felsic magmas consisting of extrusive rhyodacitic ignimbrites (about 100 m thick; total eruptive volume estimated at 10 km3: Troll, Emeleus & Donaldson, Reference Troll, Emeleus and Donaldson2000), and various dacitic to rhyolitic shallow-level intrusions that are all located close to the ring fault. The units of this early felsic phase of activity are cross-cut by basic and ultrabasic layered intrusions (gabbroic and peridotitic rocks) that form the famous Rum ultrabasic Layered Suite (LS), which is effectively a cumulate intrusion (Emeleus, Reference Emeleus1997). Small (< 1 m down to 10 cm width) picritic dykes (Upton et al. Reference Upton, Skovgaard, McClurg, Kirstein, Cheadle, Emeleus, Wadworth and Fallick2002) were intruded at a late stage in the evolution of the Rum basic–ultrabasic complex. Pre-dating both phases is the intrusion of a very coarse gabbro, now preserved as large blocks within some of the early felsic intrusions.
Figure 2. Geological sketch map of the Isle of Rum (modified from Emeleus & Bell, Reference Emeleus and Bell2005).
The Palaeogene volcanic and sub-volcanic rocks of Rum intrude an Archaean Lewisian gneiss basement overlain by Proterozoic (Torridon Group) and Mesozoic sediments (Emeleus, Reference Emeleus1997). The Archaean Lewisian Gneiss Complex has experienced several episodes of high-grade metamorphism producing, in this part of Scotland, geochemically distinct lower crustal granulite-facies and upper crustal amphibolite-facies gneisses during its long and often complex geological history (Dickin, Reference Dickin1981; Weaver & Tarney, Reference Weaver and Tarney1981; Thompson et al. Reference Thompson, Dickin, Gibson and Morrison1982; Kinny, Friend & Love, Reference Kinny, Friend and Love2005; see Park, Stewart & Wright, Reference Park, Stewart, Wright and Trewin2002, for a full review of the Lewisian). The Lewisian gneiss occurrences on Rum are confined to small outcrops within the Main Ring Fault zone. They are amphibolite-facies gneiss and comprise basic amphibolite gneiss through tonalitic to granodioritic varieties (Emeleus, Reference Emeleus1997). Compared to the mainland banded Lewisian gneisses, most of the gneisses on Rum have undergone intense thermal metamorphism as a consequence of the Paleocene magmatism (Tilley, Reference Tilley1944). The overlying Torridon Group and Mesozoic rocks on Rum do not generally show such thermal metamorphism, apart from a few spectacular cases of contact metamorphism close to magmatic plugs and dykes (Holness, Reference Holness1999; Holness & Isherwood, Reference Holness and Isherwood2003).
3. Samples and analytical methods
Most of the samples from the Isle of Rum were collected in the Northern Marginal Zone (NMZ), which is a remnant of the rim and infill of a caldera. The Northern Marginal Zone has been studied in depth by Troll, Donaldson & Emeleus (Reference Troll, Donaldson and Emeleus2004). A few samples from the Southern Mountains Zone (SMZ) have been included for comparison. Four samples from ‘Unit 9’ of the Layered Suite and one sample of a picrite dyke have been analysed to shed light on the relationship between the earliest felsic volcanic products on Rum and the Layered Suite. Three Lewisian gneiss samples from Rum have been selected to serve as analogues for the isotopic and trace element composition of local Lewisian crust.
The sampled igneous rocks include rhyodacites (UDPN; iRDP1; iRDP2; RDP12, Am-iRDP-05), microgranites (LLG; PG3; R-WG-3), dacites (NMZ-AM-7; SMZ-AM-5), basic margins to rhyodacite (SMZ-iRDPBM) and to the dacite (NMZ-AM-Bm). In addition, two pre-Layered Suite basic porphyry dykes (R-FDP-03-M; R-FDP-03-C) from the vicinity of the Layered Suite and two pre-Layered Suite gabbros from within the Am Màm Breccia (R-AMGB-1; R-AMGB-2) were analysed. The reader is referred to Figure 1 of Troll, Donaldson & Emeleus (Reference Troll, Donaldson and Emeleus2004) for detailed information on the local geology and volcanic stratigraphy. Peridotite P-U9, troctolite/allivalite A-U9, and pyroxene–gabbro PxG-UG are samples from Unit 9 of the Layered Suite. Picrite M9 is a sample from a late-stage, thin (10–20 cm wide) olivine–phyric mafic dyke (Upton et al. Reference Upton, Skovgaard, McClurg, Kirstein, Cheadle, Emeleus, Wadworth and Fallick2002). The gneiss dataset comprises the samples R-GN-1, SR321B and SMZ-AMG-01.The provenance and locality (UK Grid coordinates) of the samples are listed in Table 1. Petrographical descriptions and major-element and mineral compositions of this sample set and similar rocks from Rum can be found in Emeleus (Reference Emeleus1997) and Upton et al. (Reference Upton, Skovgaard, McClurg, Kirstein, Cheadle, Emeleus, Wadworth and Fallick2002).
Table 1. Results of isotope analysis and trace element analysis of samples from the Isle of Rum
LS – Rum Layered Suite; NMZ – Northern Marginal Zone; SMZ – Southern Mountains Zone
nd – not determined; ( ) – semiquantitative result. Trace element concentrations listed as μg/g (ppm).
Isotopic ratios are age corrected to 60.5 Ma.
Samples were prepared for chemical analysis at Trinity College, Dublin. Powders were ground in a Tema® tungsten carbide mill, but were unfortunately severely contaminated with W, Co and Ta. Results for these elements are hence not reported in the data table.
Trace elements were determined by a microwave acid digestion procedure and Inductively Coupled Plasma Quadrupole Mass Spectrometry (ICP-MS) at Leuven University (J. Mareels, unpub. Ph.D. thesis, Univ. Leuven, 2004). Sample powder aliquots of 100 mg were dissolved in a mixture of 6 ml 29M HF and 2 ml 14M HNO3 in high-pressure PTFE vessels in a Milestone Ethos900® microwave digestion system. The vessels were kept at 180 °C for 30 minutes, and at 210 °C for 10 minutes. After addition of 1 ml of In and Re or Tl monitor solutions, the sample solutions were brought to near-dryness to evaporate the largest part of the HF acid and to expel Si as volatile SiF4. Residues were taken up in 1 ml of 11M HClO4 and 1 ml of a saturated HBO3 solution in water, and evaporated again. All evaporations were carried out in a Milestone Vacuum Scrubber® unit attached to the microwave digestion system. Finally, the residues were taken up in 10 ml of 14M HNO3 by heating in closed vessels at 120 °C for 30 minutes. The solutions were quantitatively transferred to a PMP volumetric flask of 250 ml and diluted with ultrapure 18.2 MΩ water. HF and HNO3 acids were purified by sub-boiling in a PTFE distillation unit. The perchloric acid was commercially available ultrapure acid. Calibration was performed by analysing international reference silicate rocks BCR-2, AGV-1 and BE-N, and by analysing control solutions prepared from commercially available stock solutions. Microgabbro reference rock PM-S was analysed as a control sample. All reference rocks were digested and measured in the same way as rock samples from Rum. The solutions were measured with a HP-4500® ICP-Mass Spectrometer, equipped with a standard Babington nebulizer. Judging from the cross-analysis of the reference rocks, the analysis of rock PM-S, and the duplicate analysis of the picrite sample M9 (Table 2), accuracy is estimated to be better than 10% relative for most elements. Instrumental precision was variable depending on concentration level.
Table 2. Trace element analyses of picrite M9 and reference rock PM-S: comparison with literature data
Trace element concentrations listed as μg/g (ppm). nd – not determined; ( ) – semiquantitative result.
Isotopic ratios of Sr and Nd were determined at the Scottish Universities Environmental Research Centre (SUERC), East Kilbride, Scotland. Samples were accurately weighed into PFA screw-top beakers (Savillex®) and were dissolved using ultra-pure reagents in a HF–HNO3–HCl digestion sequence. Sr was separated in 2.5M HCl using Bio-Rad® AG50W X8 200–400 mesh cation exchange resin. A rare-earth-element (REE) concentrate was collected by elution of 3M HNO3. Nd was separated in a ‘cocktail’ of acetic acid (CH3COOH), methanol (CH3OH) and nitric acid (HNO3) using Bio-Rad® AG1×8 200–400 mesh anion exchange resin. Total procedure blanks for Sr and Nd were less than 0.5 ng. In preparation for mass spectrometry, Sr samples were loaded onto single Re filaments with phosphoric acid. The Sr samples were analysed on a VG Sector 54–30® multiple collector mass 200 spectrometer. A 88Sr intensity of 1V (1 × 10−11 A) ± 10% was maintained. The 87Sr/86Sr ratio was corrected for mass fractionation using 86Sr/88Sr = 0.1194 and an exponential law. The mass spectrometer was operated in the peak-jumping mode with data collected as 15 blocks of 10 ratios. NBS987 gave 0.710257 ± 18 (2 SD, n = 14). Nd samples were analysed on a VG Sector 54–30® instrument. 143Nd/144Nd ratios were measured with a 144Nd beam of 1V (1 × 10−11 A). Twelve blocks of 10 ratios were collected in the peak jumping mode and corrected for mass fractionation using an exponential law and 146Nd/144Nd = 0.7219. Repeat analyses of the internal laboratory standard (JM) gave 143Nd/144Nd = 0.511511 ± 9 (2 SD, n = 21). All isotope ratios were age corrected to 60.5 Ma according to the time of igneous eruption and emplacement of the Rum Igneous Centre (Hamilton et al. Reference Hamilton, Pearson, Thompson, Kelley and Emeleus1998; Troll et al. Reference Troll, Nicoll, Emeleus and Donaldson2008). The isotopic data are listed in Table 1.
4. Geochemistry of the various units
4.a. Sr and Nd isotopic composition
Figure 3 shows the variation of Sr and Nd isotopic compositions in the Rum samples; all isotopic ratios are age corrected to 60.5 Ma. This figure is the framework for the further discussion of the petrogenetic relationships among the samples. The felsic components of the Rum Marginal Zone have high 87Sr/86Sr (0.711–0.715) and low 143Nd/144Nd ratios (0.511–0.513), and all plot close to the composition of Lewisian amphibolite gneisses. The basic margins of the rhyodacites and the pre-Layered Suite gabbros are positioned on an apparent mixing line of rhyodacites and primitive mantle components of the mafic Layered Suite. The basic dykes in the Marginal Zone are enriched in radiogenic Sr (87Sr/86Sr = 0.7046–0.7067) relative to the most primitive LS end-member, but have only marginally lower 143Nd/144Nd ratios. This might point to hydrothermal alteration with Sr-rich and REE-poor fluids.
Figure 3. 143Nd/144Nd v. 87Sr/86Sr isotope diagram of Rum igneous rocks and Lewisian crust. The Lewisian granulite crust estimate is from Dickin (Reference Dickin1981). Two possible models are outlined. (1) The rhyodacites are binary mixtures between melts derived from Lewisian gneiss and from Lewisian granulite. (2) The rhyodacites are melts from a gneissic crust that is heterogeneous with respect to 87Sr/86Sr. Numbers along mixing lines are mass fractions (in %) of end-members. In both cases, the dacites, pre-Layered Suite gabbros and basic margins of rhyodacites are mixtures between rhyodacitic melts and mafic mantle melts. BSE represents the ‘Bulk Standard Earth’ estimate. The analytical errors are smaller than the size of the symbols. The data point for the Rum ‘primitive mantle melts’ is picrite B62/2 of Upton et al. (Reference Upton, Skovgaard, McClurg, Kirstein, Cheadle, Emeleus, Wadworth and Fallick2002). All isotopic data are age corrected to 60.5 Ma.
All the rhyodacites and two out of the three microgranites analysed have a Nd isotopic composition indistinguishable from that of the regional crust and high Sr-isotopic ratios typical of middle crustal material. Since the REE and Sr budget of the rhyodacites and microgranites is predominantly of crustal origin, it is conceptually straightforward to assume that these rocks were formed by anatexis of crust, as argued by Troll, Donaldson & Emeleus (Reference Troll, Donaldson and Emeleus2004). An interesting aspect of Figure 3 is the position of the rhyodacites on an apparent mixing line between melts derived from Lewisian amphibolite gneiss and melts from Lewisian granulites. If the rhyodacites are such melt mixtures, interpretation schemes for Rum felsic magmas have to be refined to allow for up to 30% admixture of granulite, a component previously not identified as a crustal contaminant in studies of the Rum ultrabasic Layered Suite (Palacz, Reference Palacz1985; Palacz & Tait, Reference Palacz and Tait1985; Tepley & Davidson, Reference Tepley and Davidson2003). However, such a conclusion might be premature. The data shown in Figure 3 can also be explained without invoking a granulite mixing component, if the Lewisian amphibolite gneiss is not seen as one single composition but as a group of rocks that is heterogeneous with respect to 87Sr/86Sr ratios.
Modelling of crustal contamination of mantle magmas and mixing of crustal melts often relies heavily, if not exclusively, on isotopic correlation diagrams of the 87Sr/86Sr, 143Nd/144Nd and 206Pb/204Pb–207Pb/204Pb–208Pb/204Pb systems. The primary aim of this exercise is the identification of the end-components involved in the mixing process. This is usually reasonably reliable, since the isotopic ratios are not significantly affected by processes such as partial melting and fractional crystallization. However, this strength can also become a shortcoming, because a petrogenetic model cannot be complete without specifying all the processes that determined the range of absolute abundances and variation trends of major and trace elements. To distinguish among the two models outlined in Figure 3, isotope data alone do not suffice. It is necessary to take into account and to model other geochemical properties of the samples, such as the trace element variation trends displayed in Figures 4 to 10.
Figure 4. (a) REE patterns of samples from the Rum Layered Suite (LS). Data for picrite dyke B62/2 are from Upton et al. (Reference Upton, Skovgaard, McClurg, Kirstein, Cheadle, Emeleus, Wadworth and Fallick2002). (b) REE patterns of ultramafic rocks modelled as adcumulates after 50% fractional crystalization of a liquid with REE abundances of picrite B62/2. Parameters of model are listed in Table 3. REE abundances are normalized to C1 chondrite values (Sun & McDonough, Reference Sun, McDonough, Saunders and Norry1989).
4.b. The Layered Suite
The samples from the ultrabasic Layered Suite (LS), which have been studied as representative for the mantle component on Rum, are a collection of plutonic rocks and picritic basalts. The Sr and Nd isotopic ratios (Fig. 3) are among the most ‘primitive’ (or ‘least contaminated’) values measured in the Rum Layered Suite (Palacz, Reference Palacz1985; Upton et al. Reference Upton, Skovgaard, McClurg, Kirstein, Cheadle, Emeleus, Wadworth and Fallick2002). The REE pattern of picrite M9 (Fig. 4a) is virtually identical to that reported by Upton et al. (Reference Upton, Skovgaard, McClurg, Kirstein, Cheadle, Emeleus, Wadworth and Fallick2002). This particular picrite dyke is strongly olivine–phyric cumulate rock (~ 60% olivine macrocrysts), which is clearly reflected in the high Ni and Cr concentrations and the overall low REE abundances. A much closer candidate for a primitive liquid composition is picrite B62/2 studied by Upton et al. (Reference Upton, Skovgaard, McClurg, Kirstein, Cheadle, Emeleus, Wadworth and Fallick2002), who regarded this sample as representative of a melt composition. Its REE pattern is included in Figure 4a for reference and can be considered as a good approximation (± 10%) of the REE composition of the parental liquid of the Layered Series. The concave-upward REE pattern with a maximum enrichment at Nd, and a slight depletion of La with respect to Ce and Nd, is similar to the primitive ‘Transitional-to-Enriched’ type tholeiitic basalt that is a common component of the North Atlantic Paleocene Igneous Province.
Although the petrogenesis of the Layered Suite is not the topic of the present paper, a first-order modelling of the cumulate rocks is presented in the next paragraph to show that late-stage basic liquids similar to B62/2 are viable analogues for the melts that formed the Layered Suite and that were the heat engine for crustal melting on Rum.
The very low REE abundances of the troctolite (allivalite) A-U9 and peridotite P-U9 attest to the truly adcumulus nature of the samples and allow for very little trapped intercumulus liquid. REE patterns calculated with a Rayleigh fractionation model of a liquid such as picrite B62/2 are shown in Figure 4b. The cumulus proportions and partition coefficients (calculated from “Geochemical Earth Reference Model – GERM' Database; http://earthref.org) are listed in Table 3. The results shown in Figure 4b are for aggregate cumulate compositions after 50% crystallization, without any trapped intercumulus liquid. The following equation was used (Neumann, Mead & Vitaliano, Reference Neumann, Mead and Vitaliano1954):
![\begin{equation}
[{\rm C}_{{\rm average}\,{\rm solid}} /{\rm C}_{{\rm initial\ liquid}} ] = [1/{\rm F}].[1 - (1 - {\rm F})^{\rm D} ]\end{equation}](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180419074659605-0201:S0016756809006244:S0016756809006244_eqnU1.gif?pub-status=live){kind=small}
where C = concentration of trace element i; F = degree of solidification (F = 0.5); Di = ΣKijXj = bulk solid/liquid distribution coefficient of trace element i; Kij = solid/liquid partition coefficient of trace element i for mineral j; Xj = crystallization proportion of mineral j (with ΣXj = 1). Bytownitic plagioclase dominates the REE patterns in troctolite A-U9 (positive Eu-anomaly; slight enrichment of LREE), while periodotite P-U9 and gabbro PxG-U9 show the imprint of cumulus clinopyroxene (higher HREE and lower LREE). No trapped intercumulus liquid is required to explain the higher REE contents of the pyroxene gabbro.
Table 3. Parameters of the fractional crystallization model
The model applies to the adcumulus rocks from the Layered Suite (Fig. 4) and the pre-LS gabbros (Fig. 6). LS – Rum Ultrabasic Layered Suite; NMZ – Northern Marginal Zone; PLAG – plagioclase; CPX – clinopyroxene; OLIV – olivine. Partition coefficients are a generic set calculated from data in the ‘Geochemical Earth Reference Models (GERM) Data Base' (http://earthref.org).
4.c. The basic members of the Rum Marginal Zones
The two basic porphyry dykes (samples R-FDP-03-M and R-FDP-03-C) from the contact aureole of the Layered Suite in the Northern Marginal Zone have a primitive Nd-isotopic composition (143Nd/144Nd = 0.5131; Fig. 3). The enhanced 87Sr/86Sr ratios (0.7046–0.7067) could be due to hydrothermal activity, as these dykes are heavily altered. However, the prominent La–Ce–Pr–Nd depletion, as compared to picrites B62/2 and M9 (Fig. 5c), is a primary magmatic feature as the REE are immobile during alteration. This is borne out by the smooth and virtually identical patterns of the two analysed samples. The basic dykes most likely represent one member from a range of primitive tholeiitic melts. Similar variations in REE patterns have been observed in successive units of, for example, the Upper Tholeiitic Series drilled during ODP Leg 104 on the Vøring Plateau (Viereck et al. Reference Viereck, Taylor, Parson, Morton, Hertogen, Morton and Parson1988, Reference Viereck, Hertogen, Parson, Morton, Love, Gibson, Eldholm, Thiede and Taylor1989).
Figure 5. Chondrite normalized REE patterns of (a) rhyodacites, (b) dacites and (c) basic members of the Rum marginal zones. The patterns of the basic margins of the rhyodacite intrusions are included in both (a) and (b) to show the difference between the rhyodacites and dacites. Pattern of picrite B62/2 is repeated in (c) for comparison with patterns of Figure 4. REE abundances are normalized to C1 chondrite values (Sun & McDonough, Reference Sun, McDonough, Saunders and Norry1989).
The ‘pre-Layered Series’ gabbro intrusions (R-AMGB-1 and R-AMGB-2) have REE patterns typical of plagioclase–cumulate plutonic rocks (overall-low REE; positive Eu anomalies; Fig. 5c). These gabbros do not appear to be related to the basic porphyry dykes, because the Sr and Nd isotopic ratios differ greatly (87Sr/86Sr = 0.710; 143Nd/144Nd = 0.512). Moreover, the lack of a depletion of La relative to Nd on chondrite-normalized plots is also an indication that these early gabbros did not crystallize from liquids with a marked La–Ce depletion such as the basic dykes. The similarity of isotopic composition in the gabbros and the basic margins to (rhyo)dacitic intrusions rather suggests that the gabbros crystallized from liquids that were similar to those that were quenched at the margins. Figure 6b shows that the REE patterns can be modelled quite closely by 50% crystallization of a 70% plagioclase–30% clinopyroxene adcumulate from the margin liquid, using the partition coefficients listed in Table 3, and the fractionation equation listed in Section 4.c.
Figure 6. (a) REE patterns of basic margins of (rhyo)dacitic intrusions and the pre-Layered Suite gabbros. (b) REE pattern of pre-Layered Suite gabbro modelled as a 70% plagioclase–30% clinopyroxene adcumulate from a liquid similar to that of the basic margins of the rhyodacites. REE abundances are normalized to C1 chondrite values (Sun & McDonough, Reference Sun, McDonough, Saunders and Norry1989).
4.d. The Lewisian amphibolite gneiss samples
Three Lewisian gneiss samples from the Isle of Rum were also analysed. Although the isotopic ratios (Fig. 3) and the shape of the REE patterns (Fig. 7a) are very much alike, there are significant differences in absolute abundances of trace elements, as expected for a petrographically heterogeneous group of rocks. Most Lewisian gneisses cropping out on Rum carry signs of thermal metamorphism (Emeleus, Reference Emeleus1997). Fluid loss during reheating could be the main reason for the consistent depletion of the trace elements Rb, Cs, Ba and U in sample SR321B, which was taken from a boulder in the conglomerate that was shed from the dome of the early Rum volcano. In the ensuing discussion, this sample is not considered as representatitve for the composition of the average ‘Lewisian gneiss upper crust’ beneath Rum. The results of samples R-GN-1 and SMZ-AMG-01 will be used to estimate the range of possible gneiss compositions (indicated by shaded ellipses in Figs 9, 10, 12, 14 and 15). Figure 7a also displays the REE pattern of an ‘average worldwide granulite crust’ (Rudnick & Gao, Reference Rudnick, Gao and Rudnick2003) that has Nd isotopic ratios similar to Lewisian amphibolite gneiss averages. However, the Sr isotopic composition of the granulite is much less radiogenenic (87Sr/86Sr = 0.7025; Fig. 3) due to evolution in an environment with a low time-integrated Rb/Sr ratio. The average worldwide granulite is generally depleted in LREE and mobile trace elements relative to Lewisian amphibolite gneiss from Rum.
Figure 7. REE patterns of (a) Lewisian amphibolite gneiss and (b) microgranite samples from the Rum Marginal Zone. The granulite pattern is an estimate for an ‘average’ granulite after Rudnick & Gao (Reference Rudnick, Gao and Rudnick2003). REE abundances are normalized to C1 chondrite values (Sun & McDonough, Reference Sun, McDonough, Saunders and Norry1989).
4.e. The dacites, rhyodacites and microgranites
The dacites and rhyodacites have high 87Sr/86Sr (0.711–0.715) and low 143Nd/144Nd (0.511–0.513) ratios that converge to the values of the Hebridean Lewisian amphibolite gneiss from Rum (Fig. 3). The dacites and rhyodacites have similar REE patterns, but the absolute REE abundances are consistently higher in the rhyodacites than in the dacites (Fig. 5). The rather small variation of abundances indicates that the rhyodacitic and dacitic magmas were homogeneous with respect to REE abundances, in spite of the small though significant variability of Sr isotopic ratios.
The absence of a prominent negative Eu anomaly in a rhyodacite sample led Troll, Donaldson & Emeleus (Reference Troll, Donaldson and Emeleus2004) to conclude that the composition of dacites and rhyodacitic magmas had not been greatly affected by fractional crystallization involving plagioclase. The new analyses show that the absence of a negative Eu anomaly is a persistent feature, and supports their conclusion. The isotopic data and the REE patterns are most readily understood if the dacites and rhyodacites are seen as melts generated by anatexis of crustal material. The most obvious candidate for the parent crustal material similar to the Lewisian amphibolite gneiss was sampled at Rum itself. Figure 8 shows the modelled REE pattern of partial melts of a plagioclase–alkali feldspar–hornblende–quartz gneiss with REE contents of the assumed average Rum amphibolite gneiss (average of samples R-GN-1 and SMZ-AMG-01). The parameters of the non-modal batch melting model (Shaw, Reference Shaw1970) and input data are summarized in Table 4. Although the match between the calculated pattern and the rhyodacite pattern is not perfect, the important point is made that anatexis of gneiss can produce melts with REE abundances similar to those of rhyodacites, provided that melting proceeded to an extent where little plagioclase or alkali-feldspar was left in the residue, that is, up to at least 50%. With little constraints on the amounts of water released during melting, it is difficult to decide whether such a high degree of melting is realistic or not. However, one can assume that the temperature of the ascending high-Mg picritic mantle melts that caused crustal anatexis substantially exceeded the solidus temperature of most amphibolite gneiss assemblages. The presence of a negative Eu anomaly in the calculated REE pattern is admittedly a reason for concern about the validity of the model. However, one should take into account that the theoretical equilibrium batch melting equation assumes that part of the Eu2+ released in the melt phase by melting of plagioclase + feldspar is repartitioned instantaneously into the residual plagioclase + feldspar fraction. If allowance is made for slow back-diffusion of Eu2+ into residual solids, one can expect a smaller negative Eu-anomaly in the partial melts than predicted by the equilibrium model.
Figure 8. Modelled REE pattern of 50% partial melts from a model Lewisian amphibolite gneiss (average of samples R-GN-1 and SMZ-AMG-01). A satisfactory match to the rhyodacite patterns can be obtained for degrees of partial melting in excess of 50%. Parameters of the model are listed in Table 4. REE abundances are normalized to C1 chondrite values (Sun & McDonough, Reference Sun, McDonough, Saunders and Norry1989).
Table 4. Parameters of the partial melting modelling of Lewisian amphibolite gneiss (Fig. 8)
PLAG – plagioclase; KFSP – alkali-feldspar; HBL – hornblende; QTZ – quartz; X0– starting modal abundance (mass %) of the minerals; p – fractional contribution of each mineral to the melt.
Partition coefficients are a generic set calculated from data in the ‘Geochemical Earth Reference Models (GERM) Data Base’ (http://earthref.org).
The microgranitic intrusions of the Northern Marginal Zone (Fig. 7b) have isotopic ratios and REE patterns that overlap with those of the dacites and rhyodacites. The greater variability of absolute REE abundances is not unexpected for rocks that did not rapidly congeal, but passed through a slower cooling and solidification process. Substantial variation of trace-element concentrations on the hand-specimen or even outcrop scale can result from local fluctuations in ratios of cumulate minerals to interstitial liquid, and from mobility of water-rich residual fluids. In the further discussion it is assumed that these microgranites are plutonic equivalents of the dacites and rhyodacites (Dunham, Reference Dunham1968).
4.f. Ranges of other incompatible trace elements
So far the discussion has focused on the REE elements. Figures 9 and 10 show the ranges and variation trends of other incompatible trace-elements as a function of LREE content, as represented by the Ce data.
Figure 9. Trace element variation trends of Rum igneous rocks. The shaded ellipse is an estimate of the concentration range of average Lewisian gneiss. Gneiss sample SR321B falls outside this range. Granulite data (star) from Rudnick & Gao (Reference Rudnick, Gao and Rudnick2003).
The elements Zr, Th, Ba and Rb (Figs 9, 10) follow the general trend of the LREE. Concentrations in the rhyodacites are higher than in the Lewisian gneiss, while dacites have similar or slightly higher contents than the gneiss. The contents of the three microgranites are scattered, but overlap with those of the rhyodacites and dacites. The rhyodacites form a rather homogeneous group, and for most elements the average contents of rhyodacites and dacites are distinctly different. The basic margins to (rhyo)dacites are intermediate between the rhyodacites and the picritic liquids of the Layered Series.
Figure 10. Variation of Cs, Rb and Ce in igneous rocks from Rum. The consistently low Cs contents of the rhyodacites are remarkable.
The trace element Cs (Fig. 10) quite unexpectedly shows a different behaviour. In silicate systems this element is present as the Cs+ species, which is the cation with the largest ionic radius. This element shows the greatest contrast in abundance between pristine mantle material, primary mantle melts (< 0.1 ppm), lower crust (0.3 ppm) and especially upper crust (4.9 ppm) (Rudnick & Gao, Reference Rudnick, Gao and Rudnick2003). The range in concentrations spans 2.5 orders of magnitude. The strong enrichment in the crust, and as a corollary, the strong depletion in depleted mantle, is primarily due to the large size of the Cs+ ion, which severely limits substitution in most rock-forming minerals, including the common accessories (apatite, zircon, etc.) There is, however, ready substitution for K+ in micas and, to a lesser extent, in K-feldspar. One would expect that Cs would be very strongly enriched in the liquid during melting of a gneiss rock, unless micas behaved as refractory residual phases during anatexis. The latter assumption is unlikely to be correct.
It is highly improbable that the low Cs concentrations in the rhyodacites might be due to analytical error. The data for the rhyodacites are very consistent, and have been obtained in two different analytical runs. Furthermore, the duplicate analyses of picrite M9 agreed well, at a concentration level ten times lower than in the rhyodacites. It is equally improbable that Cs was lost from the rhyodacites by hydrothermal alteration. Minor alteration of volcanic glass or K-feldspar would produce secondary phases (sericite, zeolites, etc.) that are as good or even better host phases for Cs+ as compared to the original magmatic phases. Removal of Cs+ to the extent observed in the rhyodacites would require extensive alteration in a regime of high water/rock ratios. However, this is ruled out by the petrographical characteristics of rhyodacites, for example, preservation of glass shards and delicate fiamme structures (Troll, Donaldson & Emeleus, Reference Troll, Donaldson and Emeleus2004).
5. Crustal melting and mixing involving Cs–Rb depleted Lewisian amphibolite gneiss
5.a. Indications of Cs–Rb depletion of amphibolite gneiss
Figure 11 shows the enrichment of lithophile trace elements in rhyodacites and dacites with respect to the assumed average Rum Lewisian amphibolite gneiss composition (average of samples R-GN-1 and SMZ-AMG-01). For most of the elements the enrichment factors of the rhyodacites are as expected for substantial crustal anatexis (~ 30–50%), taking into account that some lithophile trace elements (e.g. heavy REE) become less incompatible in felsic silicate systems (Mahood & Hildreth, Reference Mahood and Hildreth1983). The high Ba content and lack of a negative Eu anomaly in the rhyodacites indicate that little alkali-feldspar remained in the residual source. The enrichment factors of the dacites are generally smaller than for the rhyodacites. This either points to a higher degree of melting of gneiss, melting of another gneiss component, or mixing with a mafic component (see position of samples in Fig. 3).
Figure 11. Enrichment or depletion factors of trace elements in rhyodacites and dacites relative to the adopted Lewisian amphibolite gneiss average. Note the lack of a Rb enrichment and in particular the substantial depletion of Cs.
On the basis of known geochemical behaviour, one would have expected that the highly incompatible elements Cs and Rb were enriched by at least a factor of 2.5 in the rhyodacites relative to the amphibolite gneiss. This is obviously not the case. Figure 11 further shows that the depletion of Cs in the rhyodacites is not such an aberration after all, because the enrichment factor of Rb is also much smaller than expected for a highly incompatible element. The depletion of Cs lies neatly on a trend formed by the series Ba–Rb–Cs and can therefore be considered as systematic.
Conceptually, there are two ways to explain the low concentrations of the Rb and Cs in the rhyodacites. In a first model, one can view the rhyodacites as mixtures of partial melts of low-Cs Lewisian granulites with partial melts of Lewisian amphibolite gneiss similar to the analysed samples R-GN-1 and SMZ-AMG-01 from Rum. Such a mixing model is outlined in Figure 12b. It demonstrates that no mixture of melts derived from amphibolite gneiss and granulite is actually able to reproduce the Cs–Rb variation trend of the rhyodacites and microgranites. A similar conclusion holds with respect to the Cs–Ce variation trend. However, one could argue that the rather flat Cs–Rb trend is brought about by protracted fractional crystallization of a low-Cs mixed magma that was originally positioned on the mixing line close to the granulite melt end-member. Such an interpretation meets with serious difficulties, because a dominant contribution of melts from granulites to the original melt mixture is not consistent with the position of the rhyodacites close to the amphibolite gneiss end-member in the 143Nd/144Nd v. 87Sr/86Sr diagram (Fig. 3). Protracted fractional crystallization, with plagioclase and alkali feldspar as the main crystallizing phases, would also generate residual melts with negative Eu-anomalies, which is not the case. Finally, one could altogether question whether straightforward magma mixing is the appropriate mechanism, and argue that the geochemical trends were established by a more complex ‘Assimilation Fractional Crystallization’ (AFC) process. It is, however, very unlikely that partial melts of lower crustal granulites would have a high enough enthalpy to decompose and assimilate middle-to-upper crustal amphibolite gneiss rock to the extent that the contaminant is dominating the Sr-isotopic composition of the end-product. Here again, the argument about the absence of negative Eu anomalies is pertinent.
Figure 12. (a) Variation of Cs and Rb in igneous rocks from Rum; sample symbols as in Figures 9 and 10. (b) The mixing line between partial melts of amphibolite gneiss and of granulite does not reproduce the Cs–Rb variation trend of the rhyodacites and dacites. A two-fold enhancement of Cs and Rb in the crustal melts was assumed. (c) Semi-quantitative depletion model of amphibolite gneiss to explain the low Cs and Rb contents of rhyodacites and dacites, which are considered as melts of variably depleted amphibolite crust. If the depletion was an old event, it can account for the variable 87Sr/86Sr ratios of the rhyodacites (see Fig. 3).
A second model to explain the low concentrations of the Rb and Cs in the rhyodacites is based on the assumption that the source of the rhyodacites was similar to the amphibolite gneiss sampled on Rum, except for having been depleted in Rb and Cs prior to the melting event. This model is outlined by the evolution trends shown in Figure 12c. The required depletion for Cs is from ~ 2 ppm to 0.4 ppm. For Rb a smaller depletion from an average value of 60 ppm to 40 ppm is necessary. Partial melting of ~ 50% would enhance the concentrations to the levels measured in the rhyodacites (~ 0.8 ppm Cs and ~ 80 ppm Rb). Taken at face value, Figure 12 implies that the dacites were formed by melting of lesser Cs and Rb depleted Lewisian gneiss sources. This second model has the advantage of conceptual simplicity, and it readily accounts for the absence of negative Eu anomalies and for the convergence of isotopic composition rhyodacites to the amphibolite gneiss ratios. The model will be examined further in the following Sections.
5.b. Timing and causes of the Cs–Rb depletion
An important aspect of the model is the timing of the assumed depletion process. If it took place concurrently with or within a few millions of years prior to the crustal melting that gave rise to the rhyodacites, the depletion event would have had no consequences for the isotopic composition of the rhyodacites. In order to explain the spread of the 87Sr/86Sr ratios and the generally lower ratios of the rhyodacites relative to the amphibolite gneiss, one has to consider two different explanations. (1) The Cs–Rb depleted source regions of the rhyodacite melts were already heterogeneous with respect to 87Sr/86Sr, and had evolved in an environment with a lower time-integrated Rb/Sr ratio than did the analysed amphibolite gneiss sampled from Rum. (2) The rhyodacites are mixtures of melts derived from Lewisian granulites and from Cs–Rb depleted Lewisian amphibolite gneiss having 87Sr/86Sr ratios within the same range as the analysed amphibolite gneiss from Rum; the amphibolite gneiss source contributes more than ~ 70% by mass to the melt mixtures (see Fig. 3).
The situation is different when one assumes that the Cs–Rb depletion took place many tens or even hundreds of millions of years before Paleocene times. Figure 13 shows the effect of an ancient Cs–Rb depletion event on the Sr-isotopic evolution of the potential amphibolite gneiss source regions of the rhyodacites. The graph shows the change back in time of the 87Sr/86Sr ratio of two amphibolite gneiss formations, having 87Rb/86Sr = 1 (a common value for Hebridean amphibolites) and 87Sr/86Sr ratios of 0.717 and 0.715 at 60 Ma ago (as for the analysed samples from Rum; Fig. 3). In a second step, one can readily calculate the expected 87Sr/86Sr ratios of those two gneiss formations at 60 Ma ago, assuming they had lost between 25% and 50% of the Rb at a given moment in the past. The range of estimated percentage loss of Rb is based on the observed (~ 1.4) and expected (~ 2.3) enrichment of Rb in rhyodacites relative to average amphibolite gneiss (see Fig. 11). Figure 13 indicates that the depleted gneiss formations had at 60 Ma similar 87Sr/86Sr ratios to the rhyodacites, if the depletion event took place somewhere between 400 Ma and 600 Ma before the Paleocene. This time frame obviously makes sense, because it links the inferred Cs–Rb depletion event to Caledonian metamorphism. While Figure 13 is no proof of a Caledonian age of the depletion event, it none the less shows that this time frame and sequence of events is consistent with Sr-isotopic evolution and Cs–Rb data. Admittedly, this time frame is not readily reconciled with the evidence that rocks cropping out in the Hebridean Terrane have not been affected by Caledonian regional metamorphism at c. 450 Ma. However, Rum is situated at a distance of only 15 km to the west of the Moine Thrust Belt, the assumed western limit of Caledonian regional metamorphism. In the light of evidence (Kelley, Reddy & Maddock, Reference Kelley, Reddy and Maddock1994) for major movements at c. 430 Ma even in the Outer Isles Thrust zone about 100 km to the west of the Moine thrust, it is not inconceivable that the Caledonian Moine overthrusting initiated faulting and high-temperature fluid flow (Kerrich, La Tour & Willmore, Reference Kerrich, La Tour and Willmore1984) in deeper crustal zones at the eastern edge of the Hebridean Terrane.
Figure 13. Conceptual model of Sr-isotopic evolution of amphibolite gneiss that underwent variable degree of Rb (and Cs) depletion in the past. Two cases are presented. In the first, the undepleted amphibolite gneiss attains a 87Sr/86Sr ratio of 0.717 at 60 Ma, comparable to the high range of amphibolite gneiss sampled at Rum. If such a sample underwent 25% Rb loss c. 560 Ma ago, it would attain a value 87Sr/86Sr ratio of 0.7152 at 60 Ma; with 50% Rb loss, the ratio would grow to 0.7135. In the second case, the undepleted amphibolite gneiss has a 87Sr/86Sr ratio of 0.715 at 60 Ma (low range of amphibolite gneiss on Rum). The depleted samples would cover a 87Sr/86Sr range of 0.7113 to 0.7132. The 87Sr/86Sr values of depleted samples at 60 Ma corresponds to the observed range in felsic melts from Rum.
The other Palaeozoic magmatic episode that affected the Western Scottish Highlands area is the intrusion of Carboniferous–Permian dyke swarms of tholeiitic to alkaline composition (Baxter, Reference Baxter1987), but it is improbable that this event caused regional thermal metamorphism at mid-crustal depths. Although it is commonly assumed that Rb (and by inference Cs) is mobilized by medium- to high-grade regional metamorphism in a suprasubduction environment, the evidence for it is largely indirect and based on such information as resetting of Rb/Sr clocks and the generally low Rb content of granulite facies rocks. However, there is a substantial amount of data from experimental studies of the element mobility from subducted sediments and oceanic crust (e.g. Manning, Reference Manning2004; Zack & John, Reference Zack and John2007; Spandler, Mavrogenes & Hermann, Reference Spandler, Mavrogenes and Hermann2007). These experiments demonstrate that Cs and Rb are readily mobilized by water-rich fluids at subsolidus conditions. For completeness, it can be noted that a metamorphic event even at 400–600 Ma ago will not affect the 143Nd/144Nd system, because Nd and Sm contents are not noticeably changed by metamorphism and the time interval is too short compared to the 106 Ga half-life of 147Sm.
5.c. Modelling of the trace element and isotope variation trends
A model for the melting and mixing processes on Rum that takes account ancient Cs–Rb depletion of amphibolite gneiss has one important advantage over other models mentioned above, namely, it accounts for the range 87Sr/86Sr ratios of rhyodacites without invoking a contribution of melts from Lewisian granulites. It would indeed be somewhat unsatisfactory to invoke a substantial mixing component that apparently left no trace as a contaminant of samples from the Rum ultrabasic Layered Suite (Palacz, Reference Palacz1985; Palacz & Tait, Reference Palacz and Tait1985; Tepley & Davidson, Reference Tepley and Davidson2003). In addition, an ancient Cs–Rb depletion provides a causal link between the lower Cs concentrations and the lower 87Sr/86Sr ratios of rhyodacites as compared to the amphibolite gneiss samples. The model is further evaluated in Figures 14 and 15 for its success in consistently reproducing isotopic ratios of Sr and Nd, as well as the absolute concentrations of a typical incompatible lithophile trace element such as Ce. In both figures the enrichment factor of Rb and Ce during crustal melting is set equal to 2, that is, for approximately 50% partial melting (see Fig. 8). The model explores melt mixing processes only. Arguments against the feasibility of an ‘Assimilation Fractional Crystallization’ (AFC) model for the the origin of rhyodacites and dacites were given in Section 5.a. Admittedly, AFC could be an acceptable process to model geochemical variations in the basic dykes and early pre-Layered Suite gabbros. However, in view of the petrographical evidence for co-mingling of mafic and felsic magmas in Rum (Troll, Donaldson & Emeleus, Reference Troll, Donaldson and Emeleus2004), melt mixing was assumed to be the dominant process for all samples.
Figure 14. Modelled Ce v. 87Sr/86Sr evolution trend based on the ‘older Cs–Rb depletion’ model shown in Figure 12c. The granulite component is based on data of Dickin (Reference Dickin1981) and Rudnick & Gao (Reference Rudnick, Gao and Rudnick2003). All isotopic data are age-corrected to 60.5 Ma. The rhyodacites are modelled as 50% melts from a partially Rb depleted Lewisian source; the data allow up to 10% mixing with mafic mantle melts. The dacites, basic margins and the parent magma of the pre-Layered Suite gabbros are mixtures between mantle melts and Lewisian gneiss melts; see Figure 6 for the modelling of gabbro adcumulates. The numbers along the mixing lines are the mass fraction (in %) of the mafic melt in the hybrid rocks.
The position of the rhyodacites and dacites in the Ce v. 87Sr/86Sr diagram (Fig. 14) can be explained by melting of a gneiss that underwent variable Rb depletion. The rhyodacites might be mixed with ~ 10% basic mantle melts, the dacites with 30 to 45%. The basic margins of the (rhyo)dacites and the pre-Layered Suite gabbro, after correction for cumulate processes (Fig. 6), also plot on the theoretical mixing curves. They can be modelled as mixtures of ~ 70% basic mantle melts and ~ 30% (rhyo)dacitic melts.
The Ce v. 143Nd/144Nd diagram (Fig. 15) is actually not affected by the Cs–Rb depletion event(s) discussed above. The estimates of crustal or mantle melt contribution to the rhyodacites, dacites, pre-Layered Suite gabbros and basic margins are consistent with the estimates deduced from the Ce–87Sr/86Sr diagam (Fig. 14). Figure 15 clearly shows the differences between the rhyodacites and dacites. The dacites seem to be derived from a less Rb–Cs depleted crustal host rock (Fig. 12) and contain a greater contribution of basic melts. The latter model also offers an explanation for the lower overall REE abundances in the dacites (Fig. 5).
Figure 15. Modelled Ce v. 143Nd/144Nd evolution trend of the of Rum igneous rocks. This diagram does not depend on the gneiss depletion model. The granulite component is based on data of Dickin (Reference Dickin1981) and Rudnick & Gao (Reference Rudnick, Gao and Rudnick2003). Two mixing lines between crustal melts from amphibolite gneiss and mafic mantle melts are plotted to show the effect of a range of Ce contents of the gneiss crust and the mantle melts. The mixing estimates are consistent with the values derived in Figure 14. Note that mixing with a melt from a granulite source is not required to explain the position of the dacitic melts in this diagram.
6. Comparison of Rum with other units of the North Atlantic Volcanic Rifted Margins system
Crustal anatexis and the formation of mixtures between mantle and crust melts, either by magma mixing or advanced ‘Assimilation Fractional Crystallization’ (AFC), appear to be an integral part of igneous processes accompanying the development of Volcanic Rifted Margins. Figure 16 summarizes the 87Sr/86Sr–143Nd/144Nd data for three cases of the North Atlantic Volcanic Rifted Margins: Rum (this work), Vøring Plateau (Viereck et al. Reference Viereck, Taylor, Parson, Morton, Hertogen, Morton and Parson1988; Meyer et al. Reference Meyer, Hertogen, Pedersen, Viereck-Götte and Abratis2009) and SE Greenland (Fitton et al. Reference Fitton, Larsen, Saunders, Hardarson and Kempton2000). Available data for Skye define a trend that largely coincides with the field of the SE Greenland Lower Series (Carter et al. Reference Carter, Evensen, Hamilton and O'Nions1978; Dickin, Reference Dickin1981; Thirlwall & Jones, Reference Thirlwall, Jones, Hawkesworth and Norry1983) and Middle Series (Gibson, Reference Gibson1991). Lewisian average granulite gneiss is from Dickin (Reference Dickin1981). The composition of the ‘Moine’ meta-psammite is an average of data from Geldmacher et al. (Reference Geldmacher, Haase, Devey and Garbe-Schönberg1998); it is included as representative for the isotopic composition of medium-grade metamorphic upper continental crust material.
Figure 16. Comparison of the Rum Nd and Sr isotope ratios with available data from North Atlantic Volcanic Rifted Margin units from ODP Leg 152, SE Greenland margin, and from ODP Leg 104, Vøring Plateau, Norwegian Sea. See text for data sources. Symbols of Rum samples are as in Figure 15.
The isotopic characteristics and major-element chemistry of the Rum rhyodacites and the peraluminous rhyodacites from the Vøring Plateau are an indication that these rocks are essentially melts from crustal material, with perhaps a minor (< 15%) admixed component of mantle melts. Intermediate rocks (dacites, andesites, gabbros, etc.) represent variable mixtures of crust and mantle derived materials. Only a small fraction of the most mafic picrites, basalts and gabbros represent genuine mantle melts that largely escaped contamination with crustal material.
Rum, Vøring Plateau and SE Greenland represent different spatial/temporal formation stages of the North Atlantic Volcanic Rifted Margin system, and show a remarkably diverse record of mantle–crust interaction processes. The geochemistry of the igneous products is an important source of information on the nature, composition and approximate depth where the main interaction between ascending mantle melts and crust took place. In the case of SE Greenland, crustal anatexis and/or magma mixing of the Lower and Middle Series material most likely proceeded at different depths (granulitic, gneissic). Felsic melts on Rum were largely formed by anatexis of Lewisian amphibolite gneisses, though the inferred Cs–Rb depletion discussed above points to melting of gneisses that were affected by a higher grade of metamorphism than most outcropping gneisses. The indication that on nearby Skye, granulites (Carter et al. Reference Carter, Evensen, Hamilton and O'Nions1978; Thirlwall & Jones, Reference Thirlwall, Jones, Hawkesworth and Norry1983) and Lewisian gneisses (Gibson, Reference Gibson1991) have been involved is a reminder that local structure of the crust exerts a substantial control on mantle/crust interaction. Presumably, tectonic and structural factors had a great impact on the nature of mantle/crust interaction on the Vøring Plateau. The seaward-dipping wedges were emplaced on the edge of a thinned and stretched extended zone of the continental crust. Therefore, it may be in line with expectations that the crustal component involved in melting at the Vøring Plateau was akin to the psammitic/pelitic Moine schists that represent higher crustal levels in the North Atlantic realm.
The distinct geochemical trends seen in isotopic compositions are also observed for incompatible trace elements. A good illustration is the variation patterns displayed in the Th–Hf/3–Nb/16 ternary diagram (Fig. 17), which is a geochemical analogue of the Th–Hf/3–Ta diagram previously used by Thompson et al. (Reference Thompson, Morrison, Dickin, Gibson and Harmon1986) to establish the nature and extent of crustal contamination. Indeed, the Th-apex is the locus of low-to-medium-grade metamorphic crustal rocks, which in a continental break-up environment can either serve as the parent material in crustal anatexis or the contaminant in assimilation–fractional crystallization. The samples from Rum and the Vøring Plateau form an array that starts at the Th-depleted corner with mantle melts and ends at the Th-rich corner with crustal melts. The variation within the SE Greenland series is more complex, but there is an indication that the crustal component involved was a granulite or gneiss that had lost some Th (together with Rb and Cs) during high-grade metamorphism.
Figure 17. Th–Hf/3–Nb/16 diagram for the representative Rum samples and available data from the seaward-dipping reflector sequences of SE Greenland (Fitton et al. Reference Fitton, Larsen, Saunders, Hardarson and Kempton2000) and the mid-Norwegian Vøring margin (Meyer et al. Reference Meyer, Hertogen, Pedersen, Viereck-Götte and Abratis2009). Included in the Rum plot are the average granulite compositions (star symbol) proposed by Wood (Reference Wood1980) and Rudnick & Gao (Reference Rudnick, Gao and Rudnick2003). Felsic rocks from Rum and the Vøring margin extend into the Th-apex, which represents the Upper Crustal end-member. This is not the case for the SE Greenland margin samples. Fields for N-MORB and E-MORB are indicated; the dashed ellipse encloses the field of oceanic island basalts and continental alkaline and tholeiitic basic rocks (Wood, Reference Wood1980).
Figures 15 and 16 are interesting for other reasons as well. The most primitive and uncontaminated basaltic mantle components from the North Atlantic Volcanic Rifted Margins system have quite homogeneous 87Sr/86Sr–143Nd/144Nd isotopic compositions, and are in this respect close to the composition of typical North Atlantic MORB. This similarity is underscored by the Th–Hf/3–Nb/16 diagrams. The parent mantle melts in the three cases all plot in the field of primitive N-MORB with some spillover into the E-MORB field, particularly for the Vøring Plateau Upper Series basalts. There is not much of a close affinity with ‘within plate’ continental basalts.
7. Conclusions
87Sr/86Sr and 43Nd/144Nd isotopic compositions together with trace element data reinforce the conclusion of Troll, Donaldson & Emeleus (Reference Troll, Donaldson and Emeleus2004) that the early rhyodacitic volcanic rocks on the Isle of Rum are largely a product of crustal anatexis by ascending mantle melts. Mixing between rhyodacite magma and mafic mantle magmas gave rise to a series of variously mixed dacites, basic dyke margins and gabbros. The parent material of the rhyodacitic melts are Lewisian amphibolitic gneisses. However, the geochemistry of the strongly incompatible and mobile trace elements Cs and Rb indicates that the amphibolite gneiss parent material had experienced partial loss of Cs and Rb prior to the melting event at c. 60.5 Ma ago. The timing of the depletion cannot be determined exactly, but a Caledonian (c. 450–550 Ma) age for the Rb depletion can account for the heterogeneity of 87Sr/86Sr ratios of the rhyodacites and their lower ratios as compared to amphibolite gneiss sampled at Rum. There is no compelling need to invoke an additional lower-crustal granulitic component to account for the isotopic and trace element variations.
The late-stage picritic melts on Rum are close analogues for the mantle melts that caused crustal anatexis and were the parent melts of the Rum Layered Suite. These mantle melts have close affinities to slightly-depleted to slightly-enriched North Atlantic MORB.
Crustal anatexis is a common process in the rift-to-drift evolution during continental break-up and the formation of Volcanic Rifted Margins systems. Similar associations of ‘early felsic–later mafic’ volcanic rocks have been recovered from the now-drowned seaward-dipping wedges on the shelf of SE Greenland and on the Vøring Plateau (Norwegian Sea). The felsic–mafic rock associations of the Isle of Rum are an excellent onshore analogue to help to understand the igneous construction of Volcanic Rifted Margins systems. The igneous rocks from Rum, SE Greenland and the Vøring Plateau show geochemical differences that result from differences in the regional crustal composition and the depth at which crustal anatexis took place, but reflect a mafic MORB-type magma end-member that is remarkably similar in all three cases.
Acknowledgements
The Scottish Natural Heritage granted permission to VRT to conduct fieldwork and sampling on the Isle of Rum. Samples from Unit 9 were kindly provided by B. O'Driscoll. VRT acknowledges financial support by Science Foundation Ireland and the Royal Irish Academy. Trace element analyses were supported by EUROMARGINS CRP-01-LEC-EMA13F. RM and JH acknowledge funding by FWO-Vlaanderen (project G.0008.04) and by the government of Luxembourg (BFR05/133). The paper benefited from a constructive review by S.A. Gibson and comments from K. Goodenough.
