| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Journal of Petrology | Volume 43 | Number 7 | Pages 1287-1309 | 2002
© Oxford University Press 2002
Submarine Basalts of the Northern Kerguelen Plateau: Interaction Between the Kerguelen Plume and the Southeast Indian Ridge Revealed at ODP Site 1140
1DÉPARTEMENT DES SCIENCES DE LA TERRE ET DE LENVIRONNEMENT, UNIVERSITÉ LIBRE DE BRUXELLES, CP 160/02, AVENUE F. D. ROOSEVELT, 50, B-1050 BRUSSELS, BELGIUM
2DEPARTMENT OF EARTH, ATMOSPHERIC AND PLANETARY SCIENCES, MASSACHUSETTS INSTITUTE OF TECHNOLOGY, CAMBRIDGE, MA 02139, USA
Received June 20, 2001; Revised typescript accepted February 11, 2002
| ABSTRACT |
|---|
|
|
|---|
During Ocean Drilling Program Leg 183, basaltic cores were retrieved from the Northern Kerguelen Plateau (NKP) at Site 1140 on the extreme north of the plateau,
270 km north of the Kerguelen Archipelago. Amongst the six basement units recovered, five were pillow basalts with fresh glassy rims together with flow lobes. This is the first evidence for submarine eruption anywhere on the Kerguelen Plateau. Each flow unit has distinct geochemical and isotopic characteristics that span the range from Southeast Indian Ridge (SEIR) mid-ocean ridge basalt to tholeiitictransitional basalts derived from the Kerguelen mantle plume. Relationships between element abundance ratios involving incompatible elements and isotopic ratios reflect mixing of near-primary melts from an SEIR source and the Kerguelen plume with an increasing role for an incompatible element-rich end-member in the order: Unit 1 < 6
5 < 2 < 3; that is, there is no systematic downcore geochemical variation. Modeling results show that the depleted SEIR component is dominant, ranging from 9990% in Unit 1 to 7663% in Unit 3. Comparison with Kerguelen Archipelago flood basalts with similar geochemical characteristics indicates that in the northern part of the archipelago, where 2928 Ma tholeiitictransitional basalts are present, incorporation of asthenospheric mantle by the plume on its way to the surface may have been significant. At 34 Ma, when Site 1140 basalts were formed, the SEIR was only
50 km away. Apparently, the proximity of the spreading center and plume allowed mixing of geochemically diverse magmas. In contrast to older Cretaceous parts of the Kerguelen Plateau, the Cenozoic Northern Kerguelen Plateau appears oceanic in origin and there is no evidence for a component derived from continental crust in Site 1140 basalts. KEY WORDS: isotopic geochemistry; Kerguelen plume; Northern Kerguelen Plateau; plumeridge interaction; trace elements
| INTRODUCTION |
|---|
|
|
|---|
The Northern Kerguelen Plateau (NKP) was not sampled until 1998 when the Kerimis survey cruise dredged Skiff Bank, also known as Leclaire Rise (Fig. 1). These dredges recovered basalts and picritic basalts. In addition, the dredges included granites and syenites, very unusual rocks for an oceanic environment; these were interpreted as ice-rafted material (Weis et al., 2002
350 km west of the Kerguelen Archipelago, and Site 1140 at the extreme north end of the NKP,
270 km north of the Kerguelen Archipelago (Fig. 1). The objectives of drilling at Sites 1139 and 1140 were to establish the age, geochemical characteristics and eruption environment of the submarine igneous basement forming the NKP and to compare these results with the igneous basement forming other parts of the Kerguelen Plateau.
|
Because formation of the Kerguelen Plateau has been attributed to volcanism from the Kerguelen mantle plume (e.g. Frey et al., 2000a
), an important goal was to constrain magma production rates associated with the Kerguelen plume. The ages and lithologies recovered from the sites in the NKP are very different: at Site 1139, the igneous basement is formed by an
68 Ma (ArAr ages) alkaline lava series ranging from trachybasalt to trachyte and rhyolite (Kieffer et al., 2002
), whereas at Site 1140 the
34 Ma igneous basement is made of pillow lava flow units of tholeiitic basalt (Shipboard Scientific Party, 2000
; Duncan, 2002
). Drilling at Site 1140 recovered the first basaltic glass associated with volcanic products of the Kerguelen plume. This is also the first example of submarine volcanism on the Kerguelen Plateau, as all previous ODP sites have indicated subaerial eruption environments. The apparent scarcity of submarine volcanism on the Kerguelen Plateau contrasts strongly with the OntongJava Plateau in the western Pacific Ocean, which was entirely submarine (Mahoney et al., 1993
).
The Kerguelen Archipelago is interpreted as an
30 Ma to recent (Nicolaysen et al., 2000
) subaerial expression of the Kerguelen plume on the NKP. Weis et al. (2002)
showed an apparent age and magma production trend from Site 1140 across the Kerguelen Archipelago to Heard Island on the Central Kerguelen Plateau, including seamounts between the archipelago and Heard Island that may represent the hotspot track of the Kerguelen plume. The oldest flood basalts on the northwestern part of the Kerguelen Archipelago contain a component derived from depleted mantle (Yang et al., 1998
; Doucet et al., 2002
). Storey et al. (1988)
and Gautier et al. (1990)
proposed that interaction between the Kerguelen plume and a depleted mantle reservoir [i.e. the source of Southeast Indian Ridge (SEIR) mid-ocean ridge basalt (MORB)] decreased as the distance between the SEIR and the archipelago increased with time. Detailed investigation of several stratigraphically controlled flood basalt sections demonstrates that this model was oversimplified (e.g. Yang et al., 1998
; Doucet et al., 2002
). The compositions of these basalts show that there was a general decrease in the amount of partial melting with time, as indicated by increasing alkalinity of the lavas with decreasing eruption age, and the absence of a depleted component (relatively low 87Sr/86Sr and high 143Nd/144Nd), in <25 Ma flood basalts (e.g. Frey et al., 2000b
, 2002a
; Weis et al., 1998b
, 2002
; Mattielli et al., 2002
). Thus, the
34 Ma Site 1140 samples from the extreme north of the Kerguelen Plateau enable us to constrain the relative importance of interaction between the Kerguelen plume and the depleted upper mantle.
Finally, an important aspect of Kerguelen Plateau volcanism is the role of continental lithosphere. A role for continental lithosphere is evident in some of the Cretaceous basalts forming the Kerguelen Plateau (Mahoney et al., 1995
; Frey et al., 2000a
; Weis et al., 2001a; Ingle et al., 2002
), but not in the Cenozoic lavas forming the Kerguelen Archipelago (Weis et al., 1998b
; Yang et al., 1998
; Frey et al., 2000b
; Mattielli et al., 2002
). On the Kerguelen Archipelago, the only evidence for a component of continental lithosphere has been found in a few ultrabasic xenoliths (Hassler & Shimizu, 1998
; Mattielli et al., 1999
). Hence, an important goal is to assess the potential role of continental lithosphere in the petrogenesis of Site 1140 lavas.
This paper presents geochemical and isotopic data for the basalts that were recovered at Site 1140 on the NKP. Eruption ages for Site 1140 lavas (Duncan, 2002
) and for the flood basalts of the Kerguelen Archipelago (Nicolaysen et al., 2000
; Doucet et al., 2002
) show that the NKP is distinctly younger (Cenozoic) than the Cretaceous Central and Southern Kerguelen Plateau and Elan Bank. Site 1140 lavas show that mixing between SEIR depleted MORB and plume-derived magmas was important at
34 Ma at a distance of
50 km from the spreading SEIR axis, but there is no evidence for a continental component.
| GEOLOGICAL SETTING, BASEMENT LITHOLOGY AND PETROLOGY |
|---|
|
|
|---|
Site 1140 is at a water depth of 2394 m on the northern flank of the NKP at a distance of
270 km from the Kerguelen Archipelago (Fig. 1). The boundary between the NKP and the AustraliaAntarctic Basin lies
5 km north of Site 1140 and offsets basement by
400 m. The site is at a location with simple basement structure and a thin sedimentary section in a region of sufficient water depth for drilling. The top of the acoustic basement is horizontal and covered by a sediment sequence
350 m thick (Fig. 2; Shipboard Scientific Party, 2000
|
Drilling at Site 1140 penetrated 87·4 m of basement with
50% recovery. The basement is divided into six units based on interbedded sedimentary layers and changes in phenocryst abundance and type; these are five pillow basaltic units (Units 1, 2, 3, 5 and 6) and an interbedded sedimentary unit (Unit 4) of dolomite (Fig. 3). In Units 2 and 3, the pillows are small (40100 cm), but in Units 1 and 6, <1 m pillows are accompanied by large
5 m lobes (Figs 3 and 4). Unit 5 is more massive and is composed of three large lobes (>3 m thick). Unit 4 is a pale brown to orange dolomite and only 0·8 m were recovered, including upper and lower contacts. Because of the presence of thin sedimentary layers (chalk between Units 2 and 3 and dolomite between Units 5 and 6), these basaltic units are interpreted as representing separate eruptions (Shipboard Scientific Party, 2000
). Although the contact between Units 1 and 2 was not recovered, a probable sedimentary interlayer was recorded in the logging data (Fig. 3) and the two units have different magnetic polarity. Biostratigraphic and paleomagnetic constraints were used to infer eruption ages between 33·1 and 35·4 Ma (Shipboard Scientific Party, 2000
). This age has been confirmed by ArAr dating at
34 Ma (Duncan, 2002
).
|
|
Submarine eruption has produced distinctive textures. The outermost 12 cm of pillow rims are composed of dark brown glass, some of which is unaltered (Wallace, 2002
). Pillow margins are sparsely porphyritic with a very fine-grained to glassy groundmass that contains extremely fine-grained plagioclase microlites or small laths with ragged or swallow-tailed terminations. Olivine forms sub-equant, euhedral phenocrysts and, in the glassy rinds, some of the olivine phenocrysts contain numerous inclusions of unaltered glass and rare inclusions of chromite (Fig. 5). The massive interiors of pillows and lava lobes have the same mineral phases, but are more crystalline and coarser grained than the pillow margins, reflecting slower cooling rates in the flow interiors (especially in Units 5 and 6, where the maximum clinopyroxene and plagioclase grain size reaches 0·7 mm). Clinopyroxene and plagioclase are major groundmass phases in crystalline areas; minor phases are titanomagnetite (
5%) and intersertal glass. In some samples, glass has partially crystallized to cryptocrystalline clinopyroxene and opaque minerals. Sulfide of possible magmatic origin occurs as minute (<0·01 mm) inclusions in primary minerals and fresh glass. The Site 1140 basalts are unique amongst those recovered on Leg 183 in that they contain unaltered olivine and glass (Figs 4 and 5). Anorthite-rich plagioclase (An6070) is the dominant phenocryst type in all units and is accompanied by olivine (Figs 3 and 5), except in Unit 6, and by clinopyroxene in Units 2 and 6. Units 1 and 2 are petrographically distinct: Unit 1 contains <1% phenocrysts (except in pillow margins), whereas Unit 2 contains
20% phenocrysts. Many phenocrysts of plagioclase and clinopyroxene are strongly zoned.
|
Petrographically, Site 1140 basalts appear relatively fresh, except in the areas adjacent to pillow margins. The interstices between pillows are filled by Foraminifera-bearing dolomitic chalk and the glassy margins contain thin carbonate veins. In the freshest parts of Units 1 and 2, alteration is limited to replacement of interstitial glass by clay and the filling of vesicles by clay and carbonate (Shipboard Scientific Party, 2000
). Units 2, 3, 5 and 6 contain secondary pyrite in veins and as small disseminated grains that replace residual glass.
| ANALYTICAL TECHNIQUES |
|---|
|
|
|---|
Abundances of the major elements and the trace elements Sr, V, Ni, Zn, Ga, Y and Zr were determined by X-ray fluorescence (XRF) at the University of Massachusetts, Amherst (Rhodes, 1996
|
Whole-rock samples for Sr, Nd and Pb isotope analyses were selected to encompass the entire range of chemical composition using samples with relatively low loss on ignition (LOI) values (Tables 1 and 2). The chemical procedure used was similar to that described by Weis et al. (1987)
with improvements as discussed by Weis & Frey (1991)
. All of the samples were acid-leached to remove secondary or alteration phases. Total blank values were <1 ng for all three isotopic systems considered. Such values are negligible in view of the elemental concentrations in the samples. Sr and Nd isotopic compositions were measured on single Ta filaments and triple ReTa filaments, respectively, in the dynamic mode on a VG Sector 54 multicollector mass spectrometer with an internal precision better than 1 x 10-5 unless specified in Table 2. Sr isotopic ratios were normalized to 86Sr/88Sr = 0·1194 and Nd isotopic ratios, measured as metal, to 146Nd/144Nd = 0·7219. The average 87Sr/86Sr value of the NBS 987 Sr standard was 0·710272 ± 17 (2
m on 48 samples) and analyses of the Rennes Nd standard yielded 143Nd/144Nd = 0·511962 ± 9 (2
m on 56 samples). Pb isotopic compositions were measured on single Re filaments using the H3PO4silica gel technique. All Pb isotopic ratios were corrected for mass fractionation (0·12 ± 0·04% per a.m.u.) on the basis of 72 analyses of the NBS 981 Pb standard for a temperature range between 1090°C and 1150°C. Between-run precision was better than
0·1% for 206Pb/204Pb and 207Pb/204Pb, and better than
0·15% for 208Pb/204Pb. The reproducibility of our Pb analyses is reflected by the complete replicate analysis of Sample 1140A-32R4, 06 cm, which had surprisingly high 206Pb/204Pb for its position in Unit 5 (see discussion in Fig. 12, below).
|
|
|
|
| RESULTS |
|---|
|
|
|---|
Chemical composition
The five basaltic units at Site 1140 are tholeiitic basalts (Fig. 6). Two samples from Unit 5 show alkali enrichment, as well as moderate enrichment in Rb and U, relative to other samples from Unit 5 (Table 1 and Fig. 6). The LOI values are the highest in Unit 5, all >3 wt % (up to 5·3 wt % in one sample), whereas in the rest of the core, they vary between 1 and 2 wt %. Hence, we interpret the higher alkali content of the two Unit 5 samples as reflecting the presence of clays.
|
The MgO contents vary from 8·33 to 5·48 wt % and are fairly homogeneous for each unit (Table 1); Units 2 and 3 are characterized by lower MgO contents (<6·1 wt %; Fig. 7). Each flow unit consists of tholeiitic basalts with distinctive major and trace element contents (Figs 7 and 8). Unit 1 is the least enriched in incompatible elements and has the highest Ni and Cr contents (K2O 0·090·12 wt %, P2O5 0·11 wt %, Zr/Nb = 2927, La/Sm = 1·021·15, Cr
280 ppm and Ni
90 ppm), whereas Units 2 and 3 are the most enriched in incompatible elements and have the lowest Ni and Cr contents (K2O 0·320·58 wt %, P2O5 0·250·39 wt %, Zr/Nb = 12·810·6, La/Sm = 2·362·70, Cr <120 ppm and Ni <70 ppm). Units 5 and 6 have intermediate concentrations and ratios. Therefore, there are no systematic downcore variations in abundance ratios of incompatible elements (Fig. 9). Differences in the relative abundances of incompatible elements in the five basement units are apparent in a mantle-normalized plot (Fig. 8). The abundances of the highly incompatible elements, such as Ba, Th, U, Nb, Ta and light REE (LREE), increase in the order: Unit 1 < Unit 5
Unit 6 < Unit 2 < Unit 3 (Fig. 8). Units 2 and 3 are characterized by negative Sr anomalies, especially Unit 3. At a given MgO content, lavas from Units 2 and 3 have lower Sr/Nd and CaO/Al2O3; and lower abundances of Sc and Cr than lavas from Units 1, 5 and 6 (Fig. 10). Such results are consistent with an important role for plagioclase and clinopyroxene fractionation during the petrogenesis of the Unit 2 and 3 basalts (Fig. 3 and Shipboard Scientific Party, 2000
).
|
|
|
Incompatible element abundance ratios in Site 1140 lavas, such as Ce/Y and Nb/Zr, range from values that overlap with the oldest
30 Ma transitional flood basalts of the Kerguelen Archipelago (i.e. Units 2 and 3) to values similar to SEIR MORB (Fig. 11). The younger, mildly alkalic basalts in the Kerguelen Archipelago have ratios even higher than those of Units 2 and 3. In binary trace element ratio diagrams, Site 1140 pillow basalts plot on a mixing trend, intermediate between depleted mantle compositions for Unit 1 basalts and tholeiitictransitional basalts from the Kerguelen Archipelago for Units 2 and 3 (Yang et al., 1998
; Doucet et al., 2002
; Frey et al., 2002a
).
|
Isotope geochemistry
In agreement with the distinct compositional characteristics for each flow unit, Site 1140 units have distinct Sr, Nd and Pb isotopic ratios (Table 2 and Fig. 12ac). Unit 1 samples have very homogeneous measured 87Sr/86Sr (
0·70341) and 143Nd/144Nd (
0·51301). Unit 2 has higher 87Sr/86Sr (
0·70426) and lower 143Nd/144Nd (0·51282), and Unit 3 has even higher 87Sr/86Sr (
0·70444) and lower 143Nd/144Nd (0·51278). Units 5 and 6 have very comparable intermediate 87Sr/86Sr and 143Nd/144Nd (
0·70375 and
0·51294), although Unit 5 seems to have slightly higher 87Sr/86Sr (averages are 0·70378 vs 0·70372). The corrections for 34 Myr of in situ 87Rb and 147Sm decay are small [(26) x 10-5 for 87Sr/86Sr and (34) x 10-5 for 143Nd/144Nd] and the units show even more homogeneous ratios when corrected for decay. In a SrNd isotope diagram (Fig. 12a), Site 1140 basalts define a linear trend between the field of SEIR basalts not influenced by hotspots and the field for Kerguelen Archipelago flood basalts. Unit 1 basalts overlap with the SEIR field, whereas Unit 2 and 3 basalts overlap with the depleted end of the field for 2925 Ma Kerguelen Archipelago flood basalts. Units 5 and 6 have intermediate compositions and plot close to the field of St. Paul Island (Frey & Weis, 1995
; Weis et al., 2001b), i.e. in the field of SEIR basalts from the AmsterdamSt. Paul Plateau (Johnson et al., 2000
; K. Nicolaysen, unpublished data, 1999).
In PbPb isotope diagrams (Fig. 12b), Site 1140 lavas overlap the fields for Kerguelen Archipelago flood basalts and basalts forming the Ninetyeast Ridge, a hotspot track related to the Kerguelen plume (e.g. Weis & Frey, 1991
). The linear trends for Site 1140 lavas can be interpreted as reflecting mixing between an SEIR-related component and the Kerguelen plume as represented by the Mt. Crozier series, the most radiogenic Pb flood basalts of the Kerguelen Archipelago (Mattielli et al., 2002
; Weis et al., 2002
). The two analyzed 1140 samples that plot off this trend have clearly been enriched in U by alteration at or just before their emplacement at 34 Ma (see Fig. 8); after age-correction for in situ 238U decay, they plot along the linear trend. Similarly, in a 87Sr/86Sr vs 206Pb/204Pb diagram (Fig. 12c), sample 1140A-32R4 06 cm from Unit 5 has a very high measured 206Pb/204Pb relative to its 87Sr/86Sr and plots off the 1140 mixing trend. However, after correction for in situ decay, this sample falls within the linear mixing array (238U/204Pb > 60) (Table 2). Another sample from Unit 6 (Sample 114037R-3, 104108 cm Piece 2) is also slightly off the trend, but this shipboard sample has no UPb concentration data. Each of these three samples has high K2O concentrations for their MgO contents (Fig. 7) and it is likely that they have very high 238U/204Pb as a result of alteration.
| DISCUSSION |
|---|
|
|
|---|
Mixing of depleted and plume components
The linear correlation defined by Site 1140 lavas in isotopic ratio plots and the range of variation in <90 m of core is remarkable (Fig. 12). For example, correlation coefficients for the calculated initial ratios are 0·992 for 87Sr/86Sr143Nd/144Nd, 0·92 for 208Pb/204Pb206Pb/204Pb, 0·64 for 207Pb/204Pb206Pb/204Pb and 0·96 for 87Sr/86Sr206Pb/204Pb. If unit averages are used, the respective correlation coefficients are 0·99, 0·96, 0·86 and 0·98. The correlations defined by Site 1140 basalts in PbPb isotopic plots (Fig. 12b) extend from the field for SEIR N-MORB (Unit 1) into the field defined by flood basalts of the Kerguelen Archipelago (Unit 3). Such trends are consistent with two-component mixing between a depleted MORB component and a plume-derived component and are also observed in HfNd isotopic diagrams (Mattielli et al., 2000
The Nb/Y vs Zr/Y diagram developed by Fitton et al. (1997)
to distinguish lavas derived from the Iceland mantle plume from North Atlantic MORB is also useful in evaluating mixing components in the Kerguelen plume system. In this plot (Fig. 13), the position of the different basalt units of Site 1140 is also consistent with mixing components represented by N-MORB from the SEIR and the Kerguelen mantle plume: Unit 1 plots within the field of SEIR basalts, Units 5 and 6 plot on the line separating the plume field from the MORB field, as defined for Iceland and also valid for Kerguelen, and Units 2 and 3 plot in the plume field and overlap with the tholeiitictransitional basalts of the Kerguelen Archipelago. In a linear Nb/Y vs Zr/Y plot (inset of Fig. 13), 1140 units plot on a linear trend. Trace element and isotope systematics indicate that mixing of melts is the most likely explanation because mixing of solid mantle components is not likely to generate such strong correlations.
|
A mixing model for the 1140 basaltic units involving depleted MORB-like and plume-derived components is presented in Table 3. Inputs to the model are isotopic ratios and abundance of Sr, Nd and Pb in the end-members (SEIR MORB and Kerguelen plume magmas); outputs are the mixing proportions of the end-members required to match the observed isotopic ratios in each lava unit. Each isotopic system was fitted independently. For end-member components, two sets of data were used. For SEIR MORB, the two sets of isotopic data [SEIR-A and SEIR-B, representing respectively an Indian SEIR MORB composition calculated from literature data (Michard et al., 1986
; Dosso et al., 1988
; Johnson et al., 2000
; Mahoney et al., 2002
) and a depleted composition plotting near the depleted end of the SEIR MORB field, in Table 3] are within the field of N-MORB found by Mahoney et al. (2002)
. For the Kerguelen plume end-member, we used the average value of the Crozier basalt section in the Kerguelen Archipelago and the Crozier sample with the most radiogenic Pb (D. Weis, unpublished data, 1998). The modeling results based on SrNd isotopic ratios indicate that the relative proportion of SEIR MORB component varies from 99 to 69% in the order: Unit 1 > 6 > 5 > 2 > 3. Similar mixing calculations in the HfNd isotopic system indicate 92% of SEIR for Unit 1 and 65% in Unit 3 (Mattielli et al., 2000
). For Pb isotopic systematics, the results are more variable (from 90 to 63% for 206Pb/204Pb, from 95 to 67% for 207Pb/204Pb, and from 94 to 82% for 208Pb/204Pb), using the most radiogenic Pb end-member for the Kerguelen plume. In detail, the mixing model does not exactly fit all the data for each unit, and the relative proportion of the SEIR and Kerguelen plume components is not the same for each isotope ratio (Table 3). Nevertheless, in all cases, a large proportion (>60%) of the SEIR component is required.
Our mixing model results have implications for the genesis of volcanic sections in the northern part of the Kerguelen Archipelago, where some of the oldest basalts (2928 Ma) of tholeiitictransitional composition show considerable isotopic variation, extending from the enriched Kerguelen plume composition to isotopic compositions that are comparable with those of Units 2 and 3 of Site 1140 (Yang et al., 1998
; Doucet et al., 2002
). Mixing calculations, using the same end-members as for Site 1140 basalts (an average SEIR MORB and the average of the Crozier basalt section), indicate up to 6070% of a depleted asthenospheric component for the most depleted signatures in Kerguelen Archipelago basalts, the group D lavas (87Sr/86Sr
0·7040 and 143Nd/144Nd
0·5129). For these archipelago basalts, this component could represent depleted material that was incorporated into the plume during its ascent to the surface with the additional complexity that the Group D lavas have distinctive Sr/Nd and Ba/Th ratios and Eu anomalies, indicative of the involvement of a plagioclase-rich component (Yang et al., 1998
). For the Tourmente section, which is slightly younger (
26 Ma), and presents the transition from transitional to alkalic volcanism on the archipelago, the calculated values drop to
55% for Sr and Nd and
7080% for Pb. This section is nearly homogeneous in isotopic composition and if the basalts of this section were formed by mixing, the mixing process most probably involved solids (mantle components) rather than melts and was extremely efficient (Frey et al., 2002a
). We emphasize that the results of the mixing calculations are dependent on the end-member compositions as well as on the assumption that the end-members have homogeneous compositions. Recent SEIR MORB are fairly heterogeneous and the plume composition has probably not been homogeneous over the
115 Myr history of the Kerguelen plume. In addition, these modeled mixing proportions are strongly dependent on the Sr, Nd, Hf and Pb concentrations of the end-members and, especially in the case of the SEIR magmas, one could expect a large range of variations depending on where the magmas formed, i.e. close to or far from a hotspot or the AmsterdamSt. Paul Platform.
Mechanism for mixing plume-derived magmas and MORB-like magmas at Site 1140
Mixing between plume-derived and MORB magmas is typically manifested by along-axis geochemical gradients that reflect flow of a plume component to the spreading center axis (e.g. Schilling et al., 1985
). There are both geophysical (Small, 1995
; Yale & Phipps-Morgan, 1998
) and geochemical indications (Dosso et al., 1988
; Storey et al., 1989
; Johnson et al., 2000
; Mahoney et al., 2002
) that the Kerguelen plume has locally affected MORB erupted along the SEIR axis. At 34 Ma, Site 1140 lavas erupted
50 km away from the ridge axis (Fig. 14). The location of the Kerguelen plume at 34 Ma is uncertain, but it is inferred that the plume was coincident with the axis of the newly formed SEIR at
40 Ma (e.g. Frey et al., 2000a
). If this assumption is correct, the plume was probably within 200 km of Site 1140 at 34 Ma. It is possible that plume-derived melts at Site 1140 mixed with melts formed at the periphery of the melting regime that created oceanic crust at the ridge axis (e.g. Spiegelman, 1996
).
|
Recent geophysical studies and modeling indicate that there may be a significant Kerguelen hotspot-to-ridge flow followed by along-ridge flow (Yale & Phipps Morgan, 1998
). The Kerguelen hotspot is at present located beneath the stationary Antarctic Plate,
800 km to the SW of the AmsterdamSt. Paul Platform, which jumped closer to the hotspot than adjacent sections of the SEIR, which are
1200 km away from the Kerguelen hotspot (Small, 1995
). In the past, the Indian Plate was moving northwards at
70 mm/yr, but is now moving more slowly, at
35 mm/yr. Satellite gravity data reveal that there are transitions in the structure of the spreading center where it is influenced by the Kerguelen and the Amsterdam hotspots and the AustralianAntarctic Discordance (AAD) (Small, 1995
; Small et al., 1999
). Changes in axial morphology and segmentation may be controlled by an upper-mantle thermal gradient, possibly a flux of asthenosphere from the hotspots to the AAD (Small et al., 1999
). The rapid northward movement of the Indian Plate may also have allowed the preservation of preferential flow paths.
A possible mechanism for explaining the SEIRKerguelen plume mixing process observed in Site 1140 basalts and in the oldest flood basalts on the Kerguelen Archipelago (Yang et al., 1998
; Doucet et al., 2002
is that the proximity of a spreading center to a plume raises the temperature of the surrounding upper mantle and its viscosity, which could lead to more entrainment by the ascending plume (Hauri et al., 1994
). Once established, this mechanism may have allowed for the existence of preferential pathways between the SEIR and the plume lasting a few million years. Significant hotspotridge interactions have been documented along the SEIR. He isotopic compositions of dredged basalts along the SEIR indicate that there is either widespread contamination of the asthenosphere by the Kerguelen hotspot or localized input of some Kerguelen plume material along ridge segments to the SE of the AmsterdamSt. Paul Platform (Graham et al., 1999
). Clearly, the Kerguelen hotspot has even today a significant role on the SEIR, especially around the ASP segment.
No role for a continental component in creating the northernmost Northern Kerguelen Plateau
The older Cretaceous Kerguelen Plateau (Sites 738, 747, 749, 750 and 1137) is characterized by significantly lower 206Pb/204Pb (<18·2) than the NKP (Fig. 12b and c). For several of these sites, the low 206Pb/204Pb and higher 207Pb/204Pb have been interpreted as reflecting contamination by a continental crust component (Mahoney et al., 1995
; Frey et al., 2000a
, 2002b
; Weis et al., 2001a; Ingle et al., 2002
). In contrast to the Central and Southern Kerguelen Plateau and Elan Bank, where evidence for a continental crust component has been found, there is no geochemical evidence for a continental component in Site 1140 basalts nor in the flood basalts that form the majority of the Kerguelen Archipelago (Weis et al., 1998a
, 2002
; Yang et al., 1998
; Frey et al., 2000b
, 2002a
; Mattielli et al., 2002
).
The temporal evolution of 143Nd/144Nd in lavas associated with the Kerguelen plume is complex, but reflects a change of tectonic setting for volcanism related to the Kerguelen plume (i.e. relative plumeridgecontinent proximity) (Fig. 15). From 115 to 95 Ma, six drill sites range widely in 143Nd/144Nd. Basalts from the three sites with the lowest ratios clearly contain a continental component, i.e. Sites 738 (Mahoney et al., 1995
), 747 (Frey et al., 2002b
) and 1137 (Weis et al., 2001a; Ingle et al., 2002
). At each of these sites, low 143Nd/144Nd is also accompanied by (Nb/La)PM < 1. Basalts from these sites are all >95 Ma and were emplaced soon after the break-up of Gondwana, when continents were still relatively close to each other. These isotopic and geochemical data indicate that during the early stages of Kerguelen Plateau formation, small amounts of partially melted continental fragments were incorporated in Kerguelen plume magmas. Such low ratios (143Nd/144Nd and (Nb/La)PM < 1) are not found in the younger basalts from the Ninetyeast Ridge, Northern Kerguelen Plateau and Kerguelen Archipelago. Within the last 34 Myr, 143Nd/144Nd has decreased in the NKP, from Site 1140 to the 2522 Ma archipelago flood basalts, indicating that the plume component has become dominant as the distance between the plume and the SEIR has increased (Fig. 14). Although our sampling of the submarine NKP is still extremely limited, data for Site 1140 add to a growing body of evidence from the archipelago lavas that there was no involvement of continental crust during the formation of the NKP.
|
| CONCLUSIONS |
|---|
|
|
|---|
Site 1140 in the NKP is the first drilled site on the entire Kerguelen Plateau where there is evidence for submarine eruption. Five tholeiitic basalt units were recovered with pillow lavas. These basalts are dated at 34 Ma, which indicates that the northern NKP is much younger than the Cretaceous Central and Southern Kerguelen Plateau (>95 Ma). All geochemical and isotopic characteristics of Site 1140 tholeiitic basalts are intermediate between those of SEIR MORB and those of basalts derived from the Kerguelen plume. At 34 Ma, the SEIR was only
50 km away from the NKP and its proximity to the plume-derived material allowed extensive interaction between melts derived from the Kerguelen plume and from those forming the SEIR. The proportion of the depleted component varies between
9990% in Unit 1 and 7663% in Unit 3. As in Kerguelen Archipelago basalts
270 km to the south, there is no evidence for a component derived from continental crust in Site 1140 basalts; thus the NKP appears oceanic in origin.
| ACKNOWLEDGEMENTS |
|---|
We thank M. Lo Cascio for sample preparation, S. Ingle and J. S. Scoates for help in editing the manuscript, P. Ila for supervision of the MIT Neutron Activation Analysis Facility, C. Maerschalk for chemical processing and for isotopic analyses of 20 basalts, and S. Miller and S. Huang for their help in obtaining ICP-MS analyses. This paper benefited greatly from group discussions at Thursday seminars at ULB as well as from the numerous interactions that led to this special issue of Journal of Petrology. We thank D. Pyle and K. Johnson for their reviews. This research was supported by Belgium ARC Grant 98/03-233 from the Communauté Française de Belgique and Fonds National de la Recherche Scientifique (FNRS) Grant 2.4579.99, USSP Grant 418927-BA184 and US NSF EAR Grant 9814313. The FNRS funds the Belgian membership in the ODP program and supported the first authors participation in ODP Leg 183.
| FOOTNOTES |
|---|
*Corresponding author. Present address: Earth and Ocean Sciences, University of British Columbia, Vancouver, B.C. V6T 1Z4, Canada.Telephone: 1-604-822.1697. Fax: 1-604-822.6088. E-mail: dweis{at}eos.ubc.ca
| REFERENCES |
|---|
|
|
|---|
Barling, J., Goldstein, S. L. & Nicholls, I. A. (1994). Geochemistry of Heard Island (Southern Indian Ocean): characterization of an enriched mantle component and implications for enrichment of the sub-Indian Ocean mantle. Journal of Petrology 35, 10171053.
Dosso, L., Bougault, H., Beuzart, P., Calvez, J. Y. & Joron, J. L. (1988). The geochemical structure of the South-East Indian ridge. Earth and Planetary Science Letters 88, 4759.
Doucet, S., Weis, D., Scoates, J. S., Nicolaysen, K., Frey, F. A. & Giret, A. (2002). The depleted mantle component in the Kerguelen Archipelago basalts: petrogenesis of tholeiitictransitional basalts from the Loranchet Peninsula. Journal of Petrology 43, 13411366.
Duncan, R. A. (2002). A time frame for construction of the Kerguelen Plateau and Broken Ridge. Journal of Petrology 43, 11091119.
Fisk, M., Giovannoni, S. & Thorseth, I. (1998). Alteration of oceanic volcanic glass: textural evidence of microbial activity. Science 281, 978980.
Fitton, J. G., Saunders, A. D., Norry, M. J., Hardarson, B. S. & Taylor, R. N. (1997). Thermal and chemical structure of the Iceland plume. Earth and Planetary Science Letters 153, 197208.[Web of Science]
Frey, F. A. & Weis, D. (1995). Geochemical constraints on the origin and evolution of the Ninetyeast Ridge: a 5000 km hotspot trace in the eastern Indian Ocean. Contributions to Mineralogy and Petrology 121, 1828.
Frey, F. A., Coffin, M. F., Wallace, P. J., Weis, D., et al. (2000a). Origin and evolution of a submarine large igneous province: the Kerguelen Plateau and Broken Ridge, southern Indian Ocean. Earth and Planetary Science Letters 176, 7389.
Frey, F. A., Weis, D., Yang, H.-J., Nicolaysen, K., Leyrit, H. & Giret, A. (2000b). Temporal geochemical trends in Kerguelen Archipelago basalts: evidence for decreasing magma supply from the Kerguelen Plume. Chemical Geology 164, 6180.
Frey, F. A., Nicolaysen, K., Kubit, B. K., Weis, D. & Giret, A. (2002a). Flood basalt from Mont Tourmente in the Central Kerguelen Archipelago: the change from transitional to alkalic basalt at
25 Ma. Journal of Petrology 43, 13671387.
Frey, F. A., Weis, D., Borisova, A. Yu. & Xu, G. (2002b). Involvement of continental crust in the formation of the Cretaceous Kerguelen Plateau: new perspectives from ODP Leg 120 sites. Journal of Petrology 43, 12071239.
Gautier, I., Weis, D., Mennessier, J.-P., Vidal, P., Giret, A. & Loubet, M. (1990). Petrology and geochemistry of Kerguelen basalts (South Indian Ocean): evolution of the mantle sources from ridge to an intraplate position. Earth and Planetary Science Letters 100, 5976.
Graham, D. W., Johnson, K. T. M., Douglas Priebe, L. & Lupton, J. E. (1999). Hotspotridge interaction along the Southeast Indian Ridge near the Amsterdam and St. Paul islands: helium isotope evidence. Earth and Planetary Science Letters 167, 297310.
Hassler, D. R. & Shimizu, N. (1998). Osmium isotopic evidence for ancient subcontinental lithospheric mantle beneath the Kerguelen Islands, Southern Indian Ocean. Science 280, 418421.
Hauri, E. H., Whitehead, J. A. & Hart, S. R. (1994). Fluid dynamics and geochemical aspects of entrainment in mantle plumes. Journal of Geophysical Research 99(B12), 2427524300.[Web of Science]
Huang, S. & Frey, F. (2002). Trace element abundances of Mauna Kea basalt from Phase 2 of the Hawaiian Scientific Drilling Project: petrogenetic implications of correlations with major element content and isotopic ratios. Geochemistry, Geophysics, Geosystems (submitted).
Ila, P. & Frey, F. A. (1984). Utilization of neutron activation analysis in the study of geologic materials. In: Harling, O. K., Clark, L. & von der Hardt, P. (eds) Use and Development of Low and Medium Flux Research Reactors. Atomkernenergie Kerntechnik 44(supplement), 710716.
Ila, P. & Frey, F. A. (2000). Trace element analysis of USGS standards AGV2, BCR2, BHVO2, DTS2 and GSP2 by INAA. Journal of Radioanalysis and Nuclear Chemistry 244, 599602.
Ingle, S., Weis, D., Scoates, J. & Frey, F. A. (2002). Relationship between the early Kerguelen plume and continental flood basalts of the paleo-Eastern Gondwana margins. Earth and Planetary Science Letters 197, 3550.
Johnson, K., Graham, D., Rubin, K., Nicolaysen, K., Scheirer, D., Forsyth, D., Baker, E. & Douglas-Priebe, L. (2000). Boomerang Seamount: the active expression of the AmsterdamSt. Paul hotspot, Southeast Indian Ridge. Earth and Planetary Science Letters 183, 245259.[Web of Science]
Kieffer, B., Arndt, N. T. & Weis, D. (2002). A bimodal alkalic shield volcano on Skiff Bank; its place in the evolution of the Kerguelen Plateau. Journal of Petrology 43, 12591286.
Langmuir, C. H., Voche, R. D., Jr, Hanson, G. N. & Hart, S. R. (1978). A general mixing equation with applications to Icelandic basalts. Earth and Planetary Science Letters 37, 380392.
Macdonald, G. A. & Katsura, T. (1964). Chemical composition of Hawaiian lavas. Journal of Petrology 5, 82133.
Mahoney, J., Storey, M., Duncan, R., Spencer, K. & Pringle, M. (1993). Geochemistry and age of the Ontong Java Plateau. In: Pringle, M., Sager, W. & Sliter, W. (eds) The Mesozoic Pacific: Geology, Tectonics, and Volcanism Geophysical Monograph, American Geophysical Union., 77, 233261.
Mahoney, J. J., Jones, W. B., Frey, F. A., Salters, V. J. M., Pyle, D. G. & Davies, H. L. (1995). Geochemical characteristics of lavas from Broken Ridge, the Naturaliste Plateau and southernmost Kerguelen plateau: Cretaceous plateau volcanism in the southeast Indian Ocean. Chemical Geology 120, 315345.[Web of Science]
Mahoney, J. J., Graham, D. W., Christie, D. M., Johnson, K. T. M., Hall, L. S. & VonderHaar, D. L. (2002). Between a hotspot and a cold spot: isotopic variation in the Southeast Indian Ridge asthenosphere, 86°E118°E. Journal of Petrology 43, 11551176.
Mattielli, N., Weis, D., Scoates, J. S., Shimizu, N., Mennessier, J.-P., Grégoire, M., Cottin, J.-Y. & Giret, A. (1999). Evolution of heterogeneous lithospheric mantle in a plume environment beneath the Kerguelen Archipelago. Journal of Petrology 40, 17211744.
Mattielli, N., Weis, D., Ingle, S., Blichert-Toft, J. & Barling, J. (2000). Hafnium isotopic constraints on the source of Kerguelen Plateau lavas (Leg 183 Sites 1140 and 1137). EOS Transactions, American Geophysical Union 81, S430.
Mattielli, N., Weis, D., Blichert-Toft, J. & Albarède, F. (2002). Hf isotope evidence for a Miocene change in the Kerguelen mantle plume composition. Journal of Petrology 43, 13271339.
Michard, A., Montigny, R. & Schlich, R. (1986). Geochemistry of the mantle below the Rodriguez triple junction and the South-East Indian Ridge. Earth and Planetary Science Letters 78, 104114.
Müller, R. D., Royer, J. Y. & Lawver, L. A. (1993). Revised plate motions relative to the hotspots from combined Atlantic and Indian Ocean hotspot tracks. Geology 21, 275278.
Nicolaysen, K., Frey, F. A., Hodges, K., Weis, D. & Giret, A. (2000). 40Ar/39Ar geochronology of flood basalts, from the Kerguelen Archipelago, southern Indian Ocean: implications for Cenozoic eruptive rates of the Kerguelen Plume. Earth and Planetary Science Letters 174, 313328.
Rhodes, J. M. (1996). Geochemical stratigraphy of lava flows samples by the Hawaii Scientific Drilling Project. Journal of Geophysical Research 101(B5), 1172911746.
Schilling, J. G., Thompson, G., Kingsley, R. & Humphris, S. E. (1985). Hotspotmigrating ridge interaction in the South Atlantic. Nature 313, 187191.
Shipboard Scientific Party (2000). Site 1140. In: Coffin, M. F., Frey, F. A., Wallace, P. J., et al. (eds) Proceedings of the Ocean Drilling Program, Initial Reports, 183 (CD-ROM). College Station, TX: Ocean Drilling Program, pp. 1122.
Small, C. (1995). Observations of ridgehotspot interactions in the Southern Ocean. Journal of Geophysical Research 100(B9), 1793117946.
Small, C., Cochran, J., Sempéré, J.-C. & Christie, D. (1999). The structure and segmentation of the Southeast Indian Ridge. Marine Geology 161, 112.[Web of Science]
Spiegelman, M. (1996). Geochemical consequences of melt transport in 2-D: the sensitivity of trace elements to mantle dynamics. Earth and Planetary Science Letters 139, 115132.[Web of Science]
Storey, M., Saunders, A. D., Tarney, J., Leat, P., Thirlwall, M. F., Thompson, R. N., Menzies, M. A. & Marriner, G. F. (1988). Geochemical evidence for plumemantle interactions beneath Kerguelen and Heard Islands, Indian Ocean. Nature 336, 371374.
Storey, M., Saunders, A. D., Tarney, J., Gibson, I. L., Norry, M. J., Thirlwall, M. F., Leat, P., Thompson, R. N. & Menzies, M. A. (1989). Contamination of Indian Ocean asthenosphere by the KerguelenHeard mantle plume. Nature 338, 574576.
Sun, S.-S. & McDonough, W. F. (1989). Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. In: Saunders, A. D. & Norry, M. J. (eds) Magmatism in the Ocean Basins. Geological Society, London, Special Publications 42, 313345.
Wallace, P. J. (2002). Volatiles in submarine basaltic glasses from the Northern Kerguelen Plateau (ODP Site 1140): implications for source region compositions, magmatic processes, and plateau subsidence. Journal of Petrology 43, 13111326.
Weis, D. & Frey, F. A. (1991). Isotope geochemistry of Ninetyeast Ridge basalts: Sr, Nd, and Pb evidence for the involvement of the Kerguelen hot spot. In: Weissel, J., Peirce, J., Taylor, E., Alt, J., et al. (eds) Proceedings of the Ocean Drilling Program, Scientific Results, 121. College Station, TX: Ocean Drilling Program, pp. 591610.
Weis, D., Demaiffe, D., Cauët, S. & Javoy, M. (1987). Sr, Nd, O and H isotopic ratios in Ascension Island lavas and plutonic inclusions: cogenetic origin. Earth and Planetary Science Letters 82, 255268.
Weis, D., Frey, F. A., Leyrit, H. & Gautier, I. (1993). Kerguelen Archipelago revisited: geochemical and isotopic study of the SE Province lavas. Earth and Planetary Science Letters 118, 101119.
Weis, D., Damasceno, D., Frey, F. A., Nicolaysen, K. & Giret, A. (1998a). Temporal isotopic variations in the Kerguelen Plume: evidence from the Kerguelen Archipelago. (Goldschmidt Conference, Toulouse 98.) Mineralogical Magazine 62A, 16431644.
Weis, D., Frey, F. A., Giret, A. & Cantagrel, J. M. (1998b). Geochemical characteristics of the youngest volcano (Mount Ross) in the Kerguelen Archipelago: inferences for magma flux and composition of the Kerguelen Plume. Journal of Petrology 39, 973994.
Weis, D., Ingle, S., Nicolaysen, K., Frey, F. A., Damasceno, D., Barling, J. & Shipboard Scientific Party (2001). Origin of continental components in Indian Ocean basalts: evidence from Elan Bank (Kerguelen Plateau, ODP Leg 183, Site 1137). Geology 29, 147150.
Weis, D., Scoates, J. S., Debaille, V. & Giret, A. (2001). St. Paul Island, Indian Ocean: geochemical and isotopic constraints on the formation of the AmsterdamSt. Paul Plateau on the Southeast Indian Ridge. EOS Transactions, American Geophysical Union 82, F1169.
Weis, D., Frey, F., Schlich, R., Schaming, M., Montigny, R., Damasceno, D., Mattielli, N., Nicolaysen, K. & Scoates, J. (2002). Trace of the Kerguelen mantle plume: evidence from seamounts between the Kerguelen Archipelago and Heard Island, Indian Ocean. Geochemistry, Geophysics, Geosystems (in press).
Yale, M. M. & Phipps Morgan, J. (1998). Asthenosphere flow model of hotspotridge interactions: a comparison of Iceland and Kerguelen. Earth and Planetary Science Letters 161, 4556.
Yang, Y. J., Frey, F. A., Weis, D., Giret, A., Pyle, D. & Michon, G. (1998). Petrogenesis of flood basalts forming the Northern Kerguelen Archipelago: implications for the Kerguelen Plume. Journal of Petrology 39, 711748.
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
D. Hanano, J. S. Scoates, and D. Weis Alteration mineralogy and the effect of acid-leaching on the Pb-isotope systematics of ocean-island basalts American Mineralogist, January 1, 2009; 94(1): 17 - 26. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. S. Scoates, D. Weis, M. Franssens, N. Mattielli, H. Annell, F. A. Frey, K. Nicolaysen, and A. Giret The Val Gabbro Plutonic Suite: A Sub-volcanic Intrusion Emplaced at the End of Flood Basalt Volcanism on the Kerguelen Archipelago J. Petrology, January 1, 2008; 49(1): 79 - 105. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Yamamoto, J. P. Morgan, and W. J. Morgan Global plume-fed asthenosphere flow--II: Application to the geochemical segmentation of mid-ocean ridges Geological Society of America Special Papers, January 1, 2007; 430(0): 189 - 208. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. A. FREY, D. WEIS, A. YU. BORISOVA, and G. XU Involvement of Continental Crust in the Formation of the Cretaceous Kerguelen Plateau: New Perspectives from ODP Leg 120 Sites J. Petrology, July 1, 2002; 43(7): 1207 - 1239. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. KIEFFER, N. T. ARNDT, and D. WEIS A Bimodal Alkalic Shield Volcano on Skiff Bank: its Place in the Evolution of the Kerguelen Plateau J. Petrology, July 1, 2002; 43(7): 1259 - 1286. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. J. WALLACE Volatiles in Submarine Basaltic Glasses from the Northern Kerguelen Plateau (ODP Site 1140): Implications for Source Region Compositions, Magmatic Processes, and Plateau Subsidence J. Petrology, July 1, 2002; 43(7): 1311 - 1326. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. DOUCET, D. WEIS, J. S. SCOATES, K. NICOLAYSEN, F. A. FREY, and A. GIRET The Depleted Mantle Component in Kerguelen Archipelago Basalts: Petrogenesis of Tholeiitic-Transitional Basalts From the Loranchet Peninsula J. Petrology, July 1, 2002; 43(7): 1341 - 1366. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. A. FREY, K. NICOLAYSEN, B. K. KUBIT, D. WEIS, and A. GIRET Flood Basalt from Mont Tourmente in the Central Kerguelen Archipelago: the Change from Transitional to Alkalic Basalt at ~25 Ma J. Petrology, July 1, 2002; 43(7): 1367 - 1387. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. DAMASCENO, J. S. SCOATES, D. WEIS, F. A. FREY, and A. GIRET Mineral Chemistry of Mildly Alkalic Basalts from the 25 Ma Mont Crozier Section, Kerguelen Archipelago: Constraints on Phenocryst Crystallization Environments J. Petrology, July 1, 2002; 43(7): 1389 - 1413. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||














; Frey et al., 2000b




