Journal of Petrology Advance Access originally published online on December 3, 2004
Journal of Petrology 2005 46(3):579-601; doi:10.1093/petrology/egh089
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Oxygen and Hydrogen Isotope Stratigraphy of the Rustenburg Layered Suite, Bushveld Complex: Constraints on Crustal Contamination
1 DEPARTMENT OF GEOLOGICAL SCIENCES, UNIVERSITY OF CAPE TOWN, RONDEBOSCH 7700, SOUTH AFRICA
2 DÉPARTEMENT DE GÉOLOGIE, UMR 6524, UNIVERSITÉ JEAN MONNET, 23 RUE PAUL MICHELON, F-42023 CÉDEX 2, FRANCE
3 LABORATOIRE MAGMAS ET VOLCANS, UMR 6524, 5 RUE KESSLER, 63000 CLERMONT-FERRAND, FRANCE
4 SCHOOL OF GEOSCIENCES, UNIVERSITY OF THE WITWATERSRAND, PRIVATE BAG 3, PO WITS, 2050, SOUTH AFRICA
RECEIVED JANUARY 1, 2004; ACCEPTED OCTOBER 5, 2004
 |
ABSTRACT |
|---|
TOP
ABSTRACT
INTRODUCTIONNew
18O values for plagioclase, pyroxene and olivine, and limited whole-rock D values are presented for samples from the Rustenburg Layered Suite of the Bushveld Complex, South Africa. In combination with existing data, these provide a much more complete composite O-isotope stratigraphy for the intrusion. Throughout the layered suite, mineral 18O values indicate that the magmas from which they crystallized had 18O values that were about 7·1
, that is, 1·4 higher than expected for mantle-derived magmas, suggesting extensive crustal contamination. More limited H-isotope data suggest that the OH present within whole rocks, regardless of the degree of alteration, is of magmatic origin and not an alteration phenomenon. There appears to be no systematic change in 18O value with stratigraphic height and this requires the contamination to have taken place in a ‘staging chamber’ before emplacement of the magma(s) into the present chamber. Large amounts (30–40%) of contamination by the lower to middle crust are needed to explain these 18O values, which is in general agreement with previous estimates based on Sr- and Nd-isotope data. Alternatively, smaller amounts of contamination (
20%) by sedimentary rocks, or their partial melts, represented by the country rock can explain the data, but it is not apparent how such material could have been present at the depth of the ‘staging chamber’ in the lower to middle crust. KEY WORDS: Bushveld Complex; Rustenburg Layered Suite; oxygen isotopes; hydrogen isotopes; crustal contamination
|
INTRODUCTION |
|---|
The mafic–ultramafic component of the Bushveld Complex of South Africa is the largest such igneous intrusion on Earth (e.g. Eales & Cawthorn, 1996
) and contains some of the most important magmatic ore deposits yet discovered. The intrusion covers an area of roughly 65 000 km2 (e.g. Tankard et al., 1982) and lies almost entirely within the bounds of the sedimentary rocks of the Transvaal Supergroup (Fig. 1). The layered mafic–ultramafic rocks of the Bushveld complex have been designated the Rustenburg Layered Suite [South African Committee on Stratigraphy (SACS), 1980], henceforth abbreviated to RLS, and were emplaced at 2058·9 ± 0·8 Ma (U/Pb date on titanite, Buick et al., 2001). The RLS is most commonly subdivided using a zonal stratigraphy into a norite Marginal Zone, an ultramafic Lower Zone (LZ), an ultramafic to mafic Critical Zone (CZ), a gabbronoritic Main Zone (MZ), and a ferrogabbroic Upper Zone (UZ). The boundaries between the zones are not always defined in the same way by different researchers, but the exact position of boundaries is not of great significance to the present paper. A schematic stratigraphic column is shown in Fig. 1. Despite the large amount of existing data, no consensus has yet been reached on the petrogenesis of the Bushveld magma(s). Most workers (e.g. Davies et al., 1980; Cawthorn et al., 1981; Sharpe, 1981; Kruger, 1994; Eales & Cawthorn, 1996) are agreed that the Rustenburg Layered Suite appears to have crystallized from at least three distinct magma types, and that these magmas were affected by a significant amount of crustal contamination.
|
A large number of stratigraphic profiles have been compiled through different sections of the layered rocks; these include studies on the changes in mineral proportion and cryptic variations recorded in mineral compositions (e.g. von Gruenewaldt, 1973
; Cameron, 1978; Teigler & Eales, 1996). One of the most significant stratigraphic datasets, which has greatly influenced ideas concerning the origin and evolution of the magma(s), is that of initial Sr-isotope ratio (Hamilton, 1977; Kruger & Marsh, 1982; Harmer & Sharpe, 1985; Sharpe, 1985; Kruger, 1994). The base of the LZ records the lowest initial Sr-isotope ratios (0·7048) and thereafter the LZ shows considerable variability, reaching a high of 0·707 (Fig. 1). The CZ is also highly variable, with a sudden increase from 0·7065 to 0·7075 at the level of the Merensky Reef (Kruger & Marsh, 1982) just below the contact between the CZ and the MZ. The lower part of the MZ has variable initialSr-isotope ratios (0·7075–0·709) before becoming fairly constant (0·7085) through its upper part. At the Pyroxenite Marker (Fig. 1), just below the contact between the MZ and UZ, there is a sudden shift to lower initial Sr-isotope ratios (0·7073), which remain at a similar value throughout the UZ. The Pyroxenite Marker (von Gruenewaldt, 1973) represents an important event, because it marks a reversal in composition of both plagioclase and pyroxene, and a reversal in pyroxene mineralogy in that primary orthopyroxene reappears at this layer, whereas there is inverted pigeonite below. Together with the Sr-isotope break, these observations suggest a significant role for magma recharge. Kruger (1994) suggested that the variable Sr-isotope ratios from the LZ to the lower MZ represent an open-system ‘integration stage’ with numerous influxes of magma, whereas the upper MZ and UZ represent a closed-system ‘differentiation stage’ where the evolution of the magmas was dominated by fractional crystallization with infrequent addition of new magma or in situ contamination.
In contrast to the Rb–Sr system, there are comparatively few Sm–Nd isotope data for the Bushveld Complex. Maier et al. (2000) demonstrated that the 
Nd stratigraphy of a 4700 m section of the LZ to MZ in the western limb follows an inverse relationship to that of initial Sr-isotope ratio, and that there is a negative (although not especially strong) correlation between Nd and initial Sr-isotope ratios. Those workers concluded that the parental magmas that fed into the lower part of the intrusion had assimilated a relatively small amount of a partial melt of the crust, whereas the magmas parental to the upper part of the complex had assimilated a higher proportion of an incompatible element poor residue of that previous partial melting event.
Like Nd-isotopes, oxygen-isotope data for Bushveld mafic rocks are comparatively scarce. Previous work has shown that the 18O values of Bushveld magmas, estimated from mineral 18O values (Schiffries & Rye, 1989; Reid et al., 1993; Harris & Chaumba, 2001), are typically about 1·5 higher than the value of 5·7 expected for a mantle-derived basaltic magma (Ito et al., 1987; Eiler, 2001). Unlike Sr-isotopes, the oxygen-isotope data appear to show no systematic change with stratigraphic height. The constancy of these data was taken by those workers to suggest that the parental magmas had already assimilated a significant amount of crust before emplacement and that progressive contamination of the magma in situ did not occur to a significant degree.
At present, it is not easy to reconcile models explaining the radiogenic isotopes in terms of the influx of different magmas that had experienced variable degrees of contamination and/or different contaminants, with the lack of change in the O-isotope composition. Nevertheless, it ought to be possible to use O-isotope data in combination with existing radiogenic isotope data to produce a well-constrained model for the types of contamination process and the various contaminants involved. Oxygen isotopes have one important advantage over radiogenic isotopes in that the concentration of oxygen in the various end-members would not be expected to vary significantly. Modelling of contamination processes using oxygen isotopes produces inherently better constrained solutions than is the case for radiogenic isotopes (e.g. Sr and Nd) because the Sr and Nd concentrations are generally not known in all end-members. In the case of the Rustenburg Layered Suite, modelling of radiogenic isotopes is problematic because the rocks are cumulates. Although the initial isotope ratios in the cumulates and the liquids from which they crystallized ought to be the same, it is not a simple matter to determine the elemental concentration of Sr and Nd in the liquids based on the cumulate compositions because of the trapped liquid effect (e.g. Cawthorn, 1996). This is compounded by the problem of estimating element concentrations in the proposed contaminant, which for some zones may be a partial melt of a crustal rock.
|
PURPOSE OF STUDY |
|---|
Our first aim is to produce a more detailed O-isotope stratigraphy of the Rustenburg Layered Suite. The existing data comprise samples (n = 24) from the eastern limb (Schiffries & Rye, 1989
) with additional data for the northern limb from Harris & Chaumba (2001). Existing data for the LZ and Marginal Zone are particularly sparse (only one sample of LZ), whereas the Merensky Reef and its immediate footwall and hanging wall have been comparatively thoroughly analysed by Schiffries & Rye (1989) and Reid et al. (1993). It is clearly important to improve the density and distribution of sampling, and in this paper we combine the existing data with mineral analyses of 25 additional samples from the northern limb (Fig. 1), and 14 samples of LZ and Marginal rocks from the eastern limb. We regard the acquisition of data from the Marginal Zone and the LZ of particular importance, as they potentially provide information on the composition of the earliest magma(s) that were intruded. The resulting dataset, although composite, gives a much more complete view of the O-isotope stratigraphy.
Our second aim is to produce a model for crustal contamination that can explain both the Sr- and Nd-isotope data and the O-isotope data for the Bushveld layered rocks. In particular, it is important to explain why the radiogenic isotopes apparently vary systematically with stratigraphic height whereas O-isotopes apparently do not. Although, as discussed above, O-isotope studies permit well-constrained crustal contamination models, they are more susceptible than Sr- and particularly Nd-isotopes to change during alteration processes. The approach used in this, as in previous papers (Schiffries & Rye, 1989; Reid et al, 1993; Harris & Chaumba, 2001), has been to analyse separated minerals as opposed to whole-rock powders. The use of mineral data has several advantages over whole rocks; only fresh mineral grains are selected for analysis, and the difference in 18O value of coexisting plagioclase and pyroxene (
plagioclase–pyroxene) indicates whether or not the minerals are in oxygen-isotope equilibrium at magmatic temperatures.
It is important for this study that the effects of secondary alteration are well understood and can be eliminated as a possible cause of O-isotope variation. With this in mind, a subset of samples have been analysed for their hydrogen-isotope composition.
|
GEOLOGY AND SAMPLE SELECTION |
|---|
The Rustenburg Layered Suite (SACS, 1980
) can be divided into eastern and western limbs of approximately the same size, and a smaller northern limb (Fig. 1). The far western limb consists mainly of Marginal Zone rocks and the Bethal limb is situated under the cover of Karoo Supergroup rocks. The oxygen-isotope data of Schiffries & Rye (1989) are for samples collected predominantly in the eastern limb. Our new analyses are for LZ and Marginal Zone rocks from the Olifants River Trough (Cameron, 1978) and the Clapham Trough (Valigy, 1998) of the eastern limb. The location and a summary of the stratigraphic position of the samples are given in Fig. 1. These data are combined with new and existing data for samples from a 3 km core through the northern limb at Bellevue. The stratigraphy of the Bellevue core has been described by Knoper & von Gruenewaldt (1992), Ashwal et al. (2004) and Barnes et al. (2004), and is not repeated here. An important difference between this section through the northern limb and similar sections through the eastern and western limbs is the absence of the Pyroxenite Marker, and its associated reversal in mineral compositions and changes in pyroxene mineralogy. The core, drilled from the top of the UZ, does not extend as far as the Platreef, which is the mineralized zone found at the base of the MZ in the northern limb (Buchanan et al., 1981; Lee, 1996; Harris & Chaumba, 2001).
To relate the samples in a composite stratigraphy, the Bellevue core data and the data of Schiffries & Rye (1989) were combined using the appearance of cumulus magnetite at the UZ–MZ boundary as a common reference. The Clapham Trough samples were related to the base of the LZ, with Marginal Zone samples (i.e. stratigraphically lower than the LZ) being assigned negative height (the Marginal Zone being of the order of 220 m thick here). Two samples are from the Burgersfort area and sample 382 was taken 1 m above the contact with metasedimentary rock. The ultramafic rocks above this contact were considered by Sharpe & Hulbert (1985) to be a sill formed by ejection of an olivine-rich mush expelled from the LZ. The Olifants River Trough samples are from the middle harzburgite unit of the LZ (Cameron, 1978) and have been assigned a height of 1100 m. It should be noted that the LZ section at Clapham Trough is compressed relative to that of Olifants River Trough, suggesting that each developed as a separate ‘basin’ (e.g. Uken & Watkeys, 1997) separated by the Schwerin fold (Fig. 1).
|
PETROGRAPHY |
|---|
Modal proportions for the samples from Bellevue, Clapham Trough and Olifants River Trough sections are presented in Tables 1 and 2. The petrography and mineral chemistry of samples analysed from the Bellevue core (Fig. 1) were described by Knoper & von Gruenewaldt (1992)
and Ashwal et al. (2004). The samples include gabbros, norites, gabbronorites and anorthosites. Magnetite and Fe-rich olivine are present in some of the UZ samples and up to 10% modal quartz is present in some of the uppermost samples of the UZ. In the case of some anorthosites, it was not possible to separate sufficient pyroxene for analysis. The proportion of modal plagioclase in the analysed samples (Table 1) varies from 46 to 98% (mean 71%). Just over 400 m below the MZ–UZ boundary in the Bellevue core (at 1970·8 m), there is a feldspathic clinopyroxenite layer about 1 m thick (here termed the Pyroxenite Horizon). This horizon is not equivalent to the Pyroxenite Marker of the western and eastern limbs (Ashwal et al., 2004), as there are fundamental differences. It does not mark a reversal in mineral compositions, and it coincides with a change from primary orthopyroxene (below) to primary (now inverted) pigeonite (above), i.e. in the reverse sense to the Pyroxenite Marker (Ashwal et al., 2004). Although there is no a priori evidence to link this horizon to an input of new magma, by analogy with the western and eastern limbs, it must be close to the level where the input of new (UZ) magma occurred. Ashwal et al. (2004) suggested that up to 500 m of the uppermost MZ (including the Pyroxenite Marker), is missing from the northern limb, possibly as a result of thermal and/or mechanical erosion of the uppermost MZ by the emplacement of UZ magmas. The two lowest samples from the Bellevue core (2849·4 and 2901; Table 1) contain significant quantities of olivine (10% and 12%, respectively) and as such are atypical for the MZ, which elsewhere contains no olivine. The olivine-bearing rocks comprise four troctolitic layers found near the base of the core. The two samples analysed are from the uppermost and lowermost of these layers, which are of the order of 50 m in thickness. Ashwal et al. (2004) suggested that the primitive nature of mineral compositions in these troctolites (An70–80, En80–83, Fo75–78) is more akin to the CZ, suggesting that the troctolitic horizons might represent a sliver of CZ rocks dismembered by the intrusion of MZ magmas.
|
|
The samples from the Clapham Trough and Olifants River Trough, being from the LZ, are considerably more mafic and comprise norites, and feldspathic harzburgites and pyroxenites. The Olifants River Trough samples come from a short section that is rich in olivine (30–90%). In four samples (1352, 1532, 1535 and 1582; Table 2) it was possible to separate fresh olivine for analysis, although these olivines showed minor serpentinization along cracks and at grain boundaries. The Clapham Trough samples contain no olivine, but have much higher amounts of orthopyroxene (65–92%) than the Olifants River Trough samples. The Marginal Zone norites consist of orthopyroxene and plagioclase with variable proportions of clinopyroxene, magnetite, quartz and biotite, the latter two minerals (which are not found in the LZ, CZ or MZ rocks) suggesting some degree of local crustal contamination.
|
ANALYTICAL METHODS |
|---|
Mineral separates were prepared by hand picking clean sieved material under a binocular microscope, in some cases after initial magnetic separation. Oxygen-isotope ratios of the silicate minerals were determined at the University of Cape Town (UCT) and Université Jean Monnet (UJM) after drying powdered material in an oven at 50°C, and degassing under vacuum on conventional silicate lines at 200°C for 2 h. The silicate minerals were reacted with ClF3 (UCT) or BrF5 (UJM) and the O2 was converted to CO2 using a hot platinized carbon rod. Stable isotope ratios were measured using either a Finnigan MAT252 (UCT) or a Micromass Isoprime (UJM) mass spectrometer and are reported in the familiar
notation where = (Rsample/Rstandard – 1) x 1000 and R = 18O/16O. Duplicate splits of an internal standard (Murchison Line quartz, MQ) were run with each batch of samples in both laboratories. The 18O of MQ has been accurately determined to be 10·1 after calibration against the NBS-28 quartz standard, assuming a value for NBS-28 of 9·64 (Coplen et al., 1983). The average value obtained for MQ was used to normalize the raw data to the SMOW scale. The average difference between duplicates of MQ analysed during the course of this work was 0·11 (UCT, n = 8) and 0·20 (UJM, n = 15). These are equivalent to 1
values of 0·06 and 0·17, respectively, and represent the typical precision of the analyses. Further details of the methods employed for extraction of oxygen from silicates at UCT have been given by Vennemann & Smith (1990) and Harris & Erlank (1992); the UJM extraction procedure has been described by Gerbe & Thouret (2004). The yields of CO2 produced from each mineral were measured to confirm complete reaction. The average yield of the conventional extraction method was 98%. A smaller number of mineral separates (including all olivine samples) were analysed using laser fluorination methods at UJM using the same equipment as described by Harris et al. (2000), but with BrF5 as reagent. Replicate analyses of the Monastery garnet standard (Harris et al., 2000) suggest that the precision is comparable with that of the conventional fluorination data. Unlike the conventional analyses, which comprise many individual grains, the laser data were obtained on 1–3 individual mineral grains. The average yields for laser fluorination during the course of this work were plagioclase 97%, quartz 94%, pyroxene 97% and olivine 96%.
Hydrogen isotopes were determined at UCT using the method of Vennemann & O'Neil (1993). Whole-rock samples were degassed on the vacuum line at 200°C prior to pyrolysis. An internal water standard (CTMP; D = –9) was used to calibrate the data to the SMOW scale and a second water standard (DML; D = –300) was used to correct for scale compression (e.g. Coplen, 1993). Typical reproducibility of internal biotite standards during the period of analysis was ±2 (1). Water contents were determined either from the voltage measured on the mass 2 collector or (in the case of large samples) from the pressure measured during sample inlet using identical inlet volume to standards of known number of micromoles. Repeated measurements of water standards of known mass suggest that the typical relative error for the water content is 3%. However, it should be noted that many of the whole-rock samples analysed contain very little H2O+ and in these samples the errors might be somewhat higher. Duplicate analyses of sample 2046 gave D values of –52 and –54 and H2O+ values of 0·15 and 0·15 wt %.
|
RESULTS |
|---|
The
18O values of plagioclase and pyroxene samples from the Bellevue core are given in Table 1 and presented graphically in Figs 2–4. Both plagioclase and pyroxene show a fairly restricted range in 18O values, from 6·1 to 8·4 (mean 7·32; n = 41) and from 5·8 to 7·6 (mean 6·45; n = 32), respectively. The only exceptions are three samples (1211·91, 1560, 1745) that have plagioclase of much higher 18O value (9·1, 10·4 and 13·1, respectively), which have not been included in the average. The per mil difference () between plagioclase and pyroxene ranges from +0·6 to +1·3 (mean value +0·98; n = 31, not including the samples with abnormally high plagioclase 18O values), with two exceptions at –0·4 and +0·1 (samples 2115 and 847·42, respectively). A single olivine from one of the olivine-bearing zones at the base of the Bellevue core gave a 18O value of 6·4.
|
|
|
Data from the LZ and Marginal Zone of the eastern limb are given in Table 2. Pyroxene and olivine from the ultramafic rocks of the Olifants River Trough gave
18O values of 5·7–6·8 and 5·6–6·5, respectively. Replicate analyses were made of olivine from two samples and the difference (0·6) is somewhat larger than predicted by the normal analytical precision. This may be due to oxygen-isotope heterogeneity among olivine grains, but it is also possible that small amounts of alteration are present along cracks, which affects the 18O value of each grain to a different degree. The Marginal Zone and LZ samples from the Clapham Trough have very consistent plagioclase (mean 7·47) and pyroxene (mean 7·04) 18O values and a single Marginal Zone norite from the Schwerin fold contains plagioclase with a 18O of 7·8.
A comparison of conventional and laser fluorination data for selected samples from the Bellevue core is shown in Table 3. Although there is broad agreement between data obtained by the different methods for some samples (e.g. 1843·34 and 2849·40), the laser data sometimes differ considerably from the conventional data, particularly so for plagioclase. For the samples in Table 3, the mean 18O values for plagioclase and pyroxene by conventional analysis are 6·96 and 6·21, whereas by laser fluorination the values are 6·54 and 5·98, respectively. Thus it appears that the laser data are generally slightly lower than the conventional values. It is important to note that the conventional data represent an average of many grains, whereas the laser data often represent only one grain, having 5–15% of the mass of the sample analysed by the conventional method. Apart from analytical error, possible explanations for this apparent difference are, first, that individual mineral grains contain variable quantities of impurities and/or minor alteration phases and, second, that the 18O values of minerals are inherently heterogeneous. Variability in plagioclase 18O values within the same sample could be due to post-magmatic interaction with fluids, which is not petrographically visible. Pyroxene is more resistant to alteration, but unlike plagioclase there is the possibility of the presence of small magnetite inclusions, which are not always visible under the binocular microscope. A single magnetite was analysed, which has a much lower 18O value (–0·1, Table 3). The observed per mil difference between pyroxene and magnetite in this sample corresponds to a temperature of 515°C [using the equations of Chiba et al. (1989)]. The presence of small magnetite inclusions within pyroxene grains in the UZ might explain why the laser pyroxene analyses tend to be more variable than the conventional analyses. Variable 18O values within fresh phenocryst populations (and sometimes within individual crystals) have been recognized in volcanic rocks (e.g. Baker et al., 2000) and related to crystal accumulation during contamination. In the RLS, Prevec et al. (2004) showed that individual rocks from the Merensky and Bastard Reefs of the RLS contain minerals with variable initial Nd- and Sr-isotope ratios. It is, therefore, possible that individual minerals in the RLS have inherently heterogeneous 18O values as a result of magmatic processes. Because of the greater variability of the laser data, we have chosen to use the conventional 18O values for pyroxene and plagioclase on all plots. A more detailed study of the intra-sample variability in 18O values within RLS rocks, combined with radiogenic isotopes, is required to resolve this issue.
|
Figure 2 shows a plot of plagioclase vs pyroxene
18O values for the Bellevue core samples, Clapham Trough samples, and samples from all three zones of the layered suite analysed by Schiffries & Rye (1989). Also plotted are Merensky Reef data from Brakspruit in the western limb (E. A. Mathez & P. Agrinier, unpublished data, 2004) and the range of values for Merensky Reef footwall rocks from Impala Platinum, western limb (Reid et al., 1993). Most samples plot between the 550°C and 1150°C isotherms with a significant minority (including most of the Schiffries & Rye samples) having plagioclase–pyroxene values between 0 and 0·58. It should be noted that the three samples that have very high plagioclase 18O values are not plotted, and there are four samples that have values of plagioclase–pyroxene close to zero. The three samples for which both olivine and pyroxene have been analysed are also plotted in Fig. 2. Two of the samples show olivine–pyroxene values that are close to the predicted difference of –0·45 at 1150°C. There is no correlation (Fig. 3) between plagioclase 18O value (or pyroxene; not shown) and the modal percent plagioclase, nor is there a correlation between plagioclase–pyroxene and modal percent plagioclase.
Whole-rock hydrogen-isotope compositions and water contents for selected samples are given in Tables 1 and 2, and presented graphically in Fig. 5. The range of D values from –53 to –99 is similar to values previously obtained for the Bushveld (Mathez et al., 1994; Harris & Chaumba, 2001; Willmore et al., 2002). The Bellevue samples are notable for their relatively low whole-rock water contents (0·18–0·65 wt %) whereas some of the LZ samples have much higher water contents as a result of partial serpentinization. It should be noted that there is no correlation between water content and D value or between D value and plagioclase–pyroxene. The highly serpentinized sample 1582 with 6·86 wt % water (equivalent to about 50% serpentine) has a D value (–77) that is comparable with the D value in comparatively unaltered samples. Those samples with <0·5% H2O+ have similar average D values (–82) to those samples with >0·5% H2O+ (–76). Mathez et al. (1994) determined the D values in samples within a 40 m section intersecting the Merensky Reef at Atok (now Lebowa Platinum Mines) in the eastern limb. They found very low bulk-rock water contents (0·04–0·26 wt %, mean 0·13, n = 36) and concluded that the water resided either as structural water within pyroxene, as suggested by Bell & Rossman (1992) for mantle orthopyroxenes, or as submicroscopic phlogopite along orthopyroxene cleavage planes. The data presented in Tables 1 and 2 suggest that there is little or no difference in H-isotope composition between water in high-temperature minerals such as amphibole and biotite, and water in hydrous minerals such as serpentine formed at lower temperatures, at least in the LZ.
|
|
MAGMATIC VS HYDROTHERMAL SIGNATURES |
|---|
Because some layered intrusions (e.g. Skaergaard and Skye; Taylor & Forester, 1979
; Taylor, 1987) show shifts in 18O associated with subsolidus hydrothermal alteration (e.g. Gregory & Criss, 1986), it is important to establish whether or not the samples analysed here reflect the composition of the crystallizing cumulates, or result from subsequent fluid–rock interaction. In general, the effects of post-crystallization hydrothermal activity on plutonic rocks can be investigated by evaluating the degree of oxygen-isotope equilibrium between coexisting minerals and/or whole rocks, using so-called – plots (Gregory & Criss, 1986; Gregory et al., 1989). The most useful – diagrams plot the 18O value of a mineral that exchanges oxygen relatively rapidly vs the 18O value of a coexisting mineral that exchanges oxygen more slowly. Minerals in equilibrium in a suite of rocks are characterized by arrays that lie on an equilibrium line of constant per mil difference between the two minerals (), which are lines of constant temperature. Rock assemblages that have interacted with an external fluid will form arrays that are not parallel to these equilibrium lines because of the greater susceptibility of one of the minerals to equilibrate with the fluid.
Schiffries & Rye (1989, 1990) showed that plagioclase–pyroxene values for samples from the eastern limb of the Bushveld Complex showed no evidence for interaction with hydrothermal fluids. In Fig. 2, our new data are plotted together with those of Schiffries & Rye (1989), and it can be seen that there is a relatively high degree of internal oxygen-isotope equilibrium. The new data are generally more scattered than the existing data but the overall spreads of data are similar.
Values of plagioclase–pyroxene between 0·58 and 1·74 [from 1150 to 550°C using the plagioclase–diopside fractionation curve of Chiba et al. (1989)] can be explained by continued oxygen-isotope exchange during slow cooling (e.g. Giletti, 1986). The plagioclase samples with high 18O values have presumably been affected by post-crystallization interaction with fluids at low temperatures, which would have raised their 18O values. The two samples at the base of the Bellevue core (2849·4 and 2901), which possibly represent large xenoliths or screens of LZ or CZ material, are notable for having plagioclase 18O values slightly lower than that of pyroxene, which plot away from the main dataset, along with two LZ samples and sample 2115 in Fig. 3. These samples with negative plagioclase–pyroxene values presumably indicate O-isotope disequilibrium as a result of alteration. Sample 2849·4 has olivine–pyroxene = –0·6, which is close to the 1150°C value of 0·45 given by Chiba et al. (1989). These data suggest that, in at least some rocks, plagioclase 18O values have been lowered during alteration, in this case at high temperatures. The most important feature of the data, however, is that the plagioclase–pyroxene values of between 0·3 and 1·5 observed for the Bushveld rocks are typical of fresh gabbros worldwide (e.g. Taylor, 1968; Gregory & Criss, 1986), and imply the preservation of magmatic oxygen-isotope compositions in the vast majority of rocks of the Rustenburg Layered Suite.
The range of D values for bulk rocks and minerals of –53 to –99 obtained by previous workers has been interpreted as magmatic in origin (Mathez et al., 1994; Harris & Chaumba, 2001; Willmore et al., 2002). Harris & Chaumba (2001) suggested, on the basis of palaeomagnetically determined latitude, that meteoric water interacting with the Bushveld rocks at 2050 Ma would have had a D value of about –20. Water of this isotope composition would have produced serpentine with a D value of about –40 [assuming a temperature of 400°C, and the serpentine–water of –20 given by Suzuoki & Epstein (1976)]. The hydrogen-isotope data, therefore, do not suggest a significant role for 2050 Ma (or recent) meteoric water, even in the serpentinized LZ samples with high water content. This suggests either that the estimate for the D value of ambient meteoric water is wrong, or that the fluids responsible for serpentinization were of magmatic rather than meteoric origin.
|
OXYGEN-ISOTOPE COMPOSITION OF THE PARENT MAGMA |
|---|
It is important to note that the norites and pyroxenites of the Bushveld intrusion do not represent quenched liquid compositions, with the possible exception of the Marginal Zone rocks. It is, therefore, necessary to estimate the magma
18O value from mineral data. It was shown above (Fig. 2) that most of the analysed samples have plagioclase and pyroxene 18O values that are consistent with oxygen-isotope equilibrium at magmatic temperatures in that they plot between the closure (550°C) and crystallization (1150°C) isotherms. Let us consider the case of a bimineralic gabbro with 71% plagioclase and 29% pyroxene, as is the case for the average Bellevue sample. At the moment of crystallization from a mantle-derived basaltic magma with, for example, a 18O value of 5·70, this rock would have plagioclase with a 18O value of 5·90, pyroxene with a 18O value of 5·40, and a whole-rock 18O value of 5·76 (i.e. the cumulate has a slightly higher bulk 18O value than the magma), assuming plagioclase–pyroxene = 0·5, appropriate for plagioclase (An70) and pyroxene at 1150°C (Chiba et al., 1989) and that plagioclase–melt = +0·2 (Kyser et al., 1981). If this cumulate cooled slowly to about 750°C, a plausible closure temperature for pyroxene in a medium-grained igneous rock, plagioclase–pyroxene = 1 (Chiba et al., 1989) and the 18O values of the coexisting plagioclase and pyroxene will therefore be 6·05 and 5·05, respectively (by mass balance, assuming equal concentrations of oxygen in both minerals, with the bulk-rock 18O value remaining at 5·76). The change in pyroxene 18O value is larger than that of plagioclase because its modal abundance is less. Hence the pyroxene 18O value is now 0·65 less than that of the original magma.
The above approach cannot be used to relate mineral and magma 18O values exactly because the magnitude of the change in pyroxene and plagioclase 18O value during slow cooling is also dependent on parameters such as grain size and cooling rate (e.g. Gregory & Criss, 1986). Furthermore, the closure temperature to oxygen diffusion is not known and the rocks commonly depart significantly from bimineralic assemblages. Nevertheless, the original magma 18O value is unlikely to differ greatly from that of the bulk rock, even for rocks with extreme modal mineralogy. In Fig. 6c the bulk-rock 18O value for the Bellevue samples has been calculated from the mineral 18O values and the modal proportions, assuming that the rocks contain only plagioclase and pyroxene. In the absence of a more rigorous approach, which is not justified, these bulk-rock 18O values are assumed to approximate those of the original magmas. The lack of correlation between modal percent plagioclase and mineral 18O values (Fig. 3) supports this assumption. For rocks where only one mineral has been analysed, and for the Schiffries & Rye (1989) data, for which no modes are available, it is assumed that plagioclase–rock and pyroxene–rock are +0·35 and –0·65, respectively.
|
|
O- AND H-ISOTOPE STRATIGRAPHY |
|---|
Upper and Main Zones (Bellevue)
Although there is a certain amount of scatter in the data, particularly for plagioclase, the following general features of the oxygen-isotope stratigraphy of Bellevue are evident. Pyroxene
18O values in the Bellevue core (Fig. 4) show an overall decrease with increasing stratigraphic height, whereas plagioclase 18O values show more scatter but do not appear to show a systematic change. This feature is also seen in the data of Schiffries & Rye (1989) through a much larger stratigraphic thickness of the eastern limb of the Bushveld. In addition, between the level of the Pyroxenite Horizon and the UZ–MZ boundary (Fig. 4), there appears to be a zone with generally lower values of plagioclase–pyroxene.
As a general rule, the value of plagioclase–pyroxene would be expected to increase with stratigraphic height for two reasons.
(1) The plagioclase becomes more sodic with height (Ashwal et al., 2004; Table 3) and plagioclase–pyroxene is known to increase with decreasing anorthite content in the plagioclase (Chiba et al., 1989).
(2) It is generally understood that the crystallization temperature in the layered suite decreased with stratigraphic height (e.g. Wager & Brown 1968) and this would have resulted in an increase in plagioclase–pyroxene of the primary minerals with stratigraphic height. However, post-crystallization reaction could obscure such trends.
The data presented in Fig. 4 for the Bellevue core suggest a value of plagioclase–pyroxene of about 0·9 for the lowest analysed part of the MZ (above the olivine-bearing rocks). There is no apparent systematic change in plagioclase 18O value and because plagioclase constitutes typically 80–90% of these rocks, no change in bulk-rock 18O is implied. The predicted plagioclase–pyroxene at crystallization temperatures is 0·53 [calculated from the data of Chiba et al. (1989) assuming An70 and 1150°C]. Just below the pyroxenite horizon the plagioclase composition is An64. Plagioclase of this composition and an observed plagioclase–pyroxene of 0·9 suggest final O-isotope equilibrium at 850°C. Just above the Pyroxenite Horizon, plagioclase–pyroxene is 0·4, which implies closure to O diffusion at much higher, magmatic temperatures. The value of plagioclase–pyroxene appears to remain constant for the remainder of the MZ (Fig. 4c). In the UZ, the 18O value of plagioclase varies significantly, probably the result of interaction with fluids, but there is no indication of a systematic change. Although the data for the upper part of the UZ are scattered, values for plagioclase–pyroxene are fairly constant at about 1·3. For a plagioclase of An47 (typical for the UZ, Table 1), this corresponds to a temperature of 730°C using the Chiba et al. (1989) fractionation factors.
The rocks between the Pyroxenite Horizon and the MZ–UZ contact are much more pyroxene rich (generally about 50% pyroxene 50% plagioclase) but the change in modal proportions in these essentially bimineralic rocks should not affect plagioclase–pyroxene. Using the observed changes in modal proportions and assuming plagioclase and pyroxene 18O values of 6·2 and 7·1 (below the Pyroxenite Horizon) and 6·6 and 6·9 (above the Pyroxenite Horizon), it could be argued that there was a change in magma 18O value from 7·0 (below) to 6·8 (above) at the level of the Pyroxenite Horizon, although it should be noted that the difference is well within the analytical errors and uncertainty in estimating magma composition from mineral 18O values.
The hydrogen-isotope stratigraphy of the Bellevue core is shown in Fig. 5. There is no obvious systematic change in D with stratigraphic height, but it may be significant that the least negative D values are found in the rocks between the Pyroxenite Horizon and the MZ–UZ contact. Water contents in the MZ rocks are uniformly low (0·15–0·30 wt %), although this may, in part, be due to the fact that only the freshest looking samples were analysed. The UZ rocks have variable, but generally higher water content (up to 0·65 wt %), which is often much higher than expected given the small amount of biotite and/or amphibole usually present (Table 1). This is consistent with the occurrence of significant quantities of hydroxyl-bearing minerals as very small inclusions within, for example, pyroxene and/or the presence of water within the pyroxene structure as discussed by Mathez et al. (1994).
Lower and Marginal Zones (Clapham and Olifants River Troughs)
The Lower Zone samples (Table 2) have more varied 18O values than the rocks from elsewhere in the Bushveld Complex. Kruger (1994) documented fairly large variations in initial Sr-isotope ratio in the LZ (0·7047–0·7072), which suggest that greater variation in magma 18O than the MZ and UZ might also have existed. Both olivine and pyroxene in the Olifants River Trough samples, and pyroxene and plagioclase in the Clapham Trough samples, have variable 18O values. The Olifants River Trough samples have pyroxene 18O values of between 5·7 and 6·8