Journal of Petrology Advance Access originally published online on April 17, 2007
Journal of Petrology 2007 48(6):1155-1183; doi:10.1093/petrology/egm013
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Constraining Fluid and Sediment Contributions to Subduction-Related Magmatism in Indonesia: Ijen Volcanic Complex
1Department of Earth Sciences, Durham University, Durham DH1 3LE, UK
2Department of Earth Sciences, Bristol University, Bristol BS8 1RJ, UK
3Department Of Geology, Royal Holloway University Of London, Egham Tw20 0ex, UK
RECEIVED JULY 10, 2006; ACCEPTED MARCH 7, 2007
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTION AND SETTINGOVERVIEW OF THE IJEN...SAMPLE GROUPINGSANALYTICAL TECHNIQUESGEOCHEMISTRYDISCUSSIONSUMMARY AND CONCLUSIONSSUPPLEMENTARY DATAREFERENCES |
|---|
Ijen Volcanic Complex (IVC) in East Java, Indonesia is situated on thickened oceanic crust within the Quaternary volcanic front of the Sunda arc. The 20 km wide calderas complex contains around 22 post-caldera eruptive centres, positioned either around the caldera-rim (CR) or along a roughly NE–SW lineament inside the caldera (IC). The CR and IC lavas exhibit separate differentiation histories. Major element and trace element modelling shows that fractionation of olivine, clinopyroxene, Fe–Ti oxide ± plagioclase is important in the CR group, whereas plagioclase is the dominant fractionating mineral in the same assemblage for the IC group. Spatial controls on magmatic differentiation highlight important structural controls on magma storage and ascent at IVC. Mantle-like
18O values, restricted ranges in Sr, Nd and Hf isotope ratios, and a lack of correlation between isotope ratios and indices of differentiation in IVC lavas are consistent with magmatic evolution through fractional crystallization. Furthermore, the small ranges in isotopic ratios throughout the complex indicate little heterogeneity in the mantle. IVC lavas possess higher 176Hf/177Hf and 143Nd/144Nd isotope ratios than other volcanoes of Java, representing the least contaminated source so far analysed and, therefore, give the best estimate yet of the pre-subduction mantle wedge isotopic composition beneath Java. Trace element and radiogenic isotope data are consistent with a two-stage, three-component petrogenetic model for IVC, whereby an Indian-type mid-ocean ridge basalt (I-MORB)-like fertile mantle wedge is first infiltrated by a small percentage of fluid, sourced from the altered oceanic crust, prior to addition of <1% Indian Ocean sediment dominated by pelagic material. KEY WORDS: differentiation; geochemistry; source components; Sr, Nd, Hf and O isotopes; Sunda arc
|
INTRODUCTION AND SETTING |
|---|
|
TOPABSTRACTINTRODUCTION AND SETTINGOVERVIEW OF THE IJEN...SAMPLE GROUPINGSANALYTICAL TECHNIQUESGEOCHEMISTRYDISCUSSIONSUMMARY AND CONCLUSIONSSUPPLEMENTARY DATAREFERENCES |
|---|
Understanding magma genesis and evolution in subduction zone environments is crucial to understanding the formation of the continents and crustal recycling in the mantle. However, determining original source compositions of volcanic arc rocks can be challenging, as lavas rarely reach the Earth's surface without experiencing some processes that modify their composition; for example, crystal fractionation, magma mixing or crustal contamination (Davidson, 1996
; Thirlwall et al., 1996; Mandeville et al., 1996; Reubi & Nicholls, 2004; Davidson et al., 2005). In most petrogenetic models, island arc magmas originate from a mantle wedge that is modified by slab-derived components (Hawkesworth et al., 1991; McCulloch & Gamble, 1991; Pearce & Peate, 1995). The magnitude and nature of the contribution from the subducting slab is debated, although there is a consensus that a slab-derived fluid is involved, consisting of fluids from either the altered oceanic crust (Hawkesworth et al., 1997; Turner et al., 1997; Turner & Foden, 2001) or overlying sediments (Class et al., 2000), and other components, possibly partial melts, derived from subducted sediment (Edwards et al., 1993; Elliott et al., 1997; Turner & Foden, 2001; Vroon et al., 2001).
Ijen volcanic complex (IVC) is located on the eastern edge of Java within the Quaternary volcanic front of the Sunda arc (Fig. 1). The arc forms the western part of the Indonesian subduction zone system, which extends over 3000 km from the Andaman Islands north of Sumatra to Flores in the Banda Sea. It formed as a result of the northward subduction of the Indo-Australian Plate beneath the Eurasian Plate at a rate of around 6–7 cm/year (DeMets et al., 1990; Tregoning et al., 1994). Despite this apparently simple tectonic setting, understanding magma genesis and evolution at the Sunda arc is complicated by the variable nature of the arc crust, the changing age of the subducting oceanic crust, which increases eastwards from 
80 to 130 Ma (Plank & Langmuir, 1998), and variation in the type and amount of subducted sediment, which ranges from bimodal siliceous ooze and clays with minor turbidites beneath the western part of the arc to strongly clay-dominated material with secondary carbonate and minor siliceous ooze in the east (Plank & Langmuir, 1998). Along the length of the arc the composition and thickness of the overriding Eurasian plate is also thought to change, from continental in the west to oceanic in the east (Hamilton, 1979). Beneath Java, the crust is 20 km thick and has a seismic velocity structure intermediate between continental and oceanic (Whitford, 1975). The Java crust is said to consist of ophiolite slivers, mélange and older volcanic rocks (Hamilton, 1979). The Sunda Land (SE Asian continental part of the Sunda block–Eurasian plate with pre-Tertiary basement) boundary shown in Fig. 1 suggests that in western Java the island arc is built on continental material. In eastern Java beneath the IVC (outside the hypothesized Sunda Land boundary) the overriding plate comprises thickened oceanic crust. However, recent research on inherited zircons (Smyth, 2005) provides evidence for old continental basement beneath the Southern Mountains Arc (south of the active volcanic axis) in East Java.
|
Despite a number of petrogenetic studies of Javan volcanoes (Edwards et al., 1993
; Turner & Foden, 2001; Gertisser & Keller, 2003; and references therein), uncertainty still prevails over the nature of the subduction component and the mechanisms by which it is added to the mantle wedge. For example, Edwards et al. (1993) proposed a homogeneous slab contribution along the Sunda arc, whereas Turner & Foden (2001) identified along-arc heterogeneity in this component. The IVC is located on the eastern edge of Java within the Quaternary volcanic front of the Sunda arc (Fig. 1) and is therefore in an ideal location to investigate the mantle source composition beneath Java; the relatively thin crust on which the volcano sits minimizes the likelihood of high-level continental crust contamination (compared with studies in West Java) and potentially allows us to constrain better crustal contributions in the mantle source. Furthermore, apparent spatial controls on some major and trace element variations at IVC mean that the volcanic complex provides a useful opportunity to investigate the relationship between volcanic structure and shallow-level processes influencing the composition of arc lavas. The aims of this study are, therefore, two-fold: (1) to identify the impact, if any, of volcanic and subjacent structure upon magmatic differentiation; (2) to characterize the IVC magma source and place constraints on the nature of the slab component(s) and the transfer mechanism(s) to the mantle wedge. |
OVERVIEW OF THE IJEN VOLCANO COMPLEX |
|---|
|
TOPABSTRACTINTRODUCTION AND SETTINGOVERVIEW OF THE IJEN...SAMPLE GROUPINGSANALYTICAL TECHNIQUESGEOCHEMISTRYDISCUSSIONSUMMARY AND CONCLUSIONSSUPPLEMENTARY DATAREFERENCES |
|---|
The Ijen volcanic complex consists of several stratovolcanoes and cinder cones constructed within and around a 20 km wide caldera (Fig. 2). The complex takes its name from Kawah Ijen volcano, the only volcano of the complex that is currently active. Based on an episode of caldera collapse, the stratigraphic succession has been divided into pre-caldera, caldera and post-caldera groups (van Bemmelen, 1949
). Little is known about the age, structure or volcanic history of the pre-caldera cone (Old Ijen); it is thought to be either early (van Bemmelen, 1949) or late (Sitorus, 1990; K–Ar dating) Pleistocene. Sitorus (1990) considered the pre-caldera volcanic structure to be a large, single, composite volcano whereas Berlo (2001), suggested that the pre-caldera collapse volcano may have had multiple vents. Kemmerling (1921) has suggested the possibility that twin volcanoes made up Old Ijen, because of the oval shape of the subsequent caldera. The deposits of Old Ijen disconformably overlie Miocene limestone (van Bemmelen, 1949) and consist of pyroclastic flow deposits, a series of pumice airfall deposits and several lava flows. The caldera wall forms a prominent arcuate ridge in the north, but elsewhere in the complex the caldera rim is buried by the post-caldera volcanoes (Fig. 2). The collapse of the pre-caldera cone(s) is thought to have occurred in stages more than 50 000 years ago (Sitorus, 1990). The caldera group consists of two main units, comprising large volumes of caldera-forming ignimbrite and minor layers of pumice airfall and lahar deposits. Most of the eruptive products were deposited to the north (van Bemmelen, 1949). Shortly after caldera collapse a resurgent dome (Blau) formed inside the caldera. The post-caldera group is characterized by many types of eruption (phreatomagmatic, phreatic, strombolian and plinian) and comprises 58 lithological volcanic units, erupted from 22 separate vents (Sitorus, 1990). The post-caldera eruptive centres can be further divided on the basis of their geographical position within the complex: those located on the caldera rim (CR), dominantly composite cones, or intra-caldera (IC) centres situated inside the southern part of the caldera along a ENE–WSW lineation (Fig. 2). The latter are predominantly cinder cones with minor composite cones and a lava dome. The IC volcanoes are generally younger (with the exception of Blau) than the CR volcanoes (Sitorus, 1990; Berlo, 2001). Faulting is evident within the caldera complex in the north (Fig. 2) and large faults have affected the caldera-rim volcanoes of Merapi, Ringgih and Djampit. Present volcanic activity is limited to Ijen crater (Kawah Ijen), an acidic crater lake with a maximum diameter of 1·2 km and a depth of around 200–300 m. The last major eruption created phreatomagmatic air fall and pyroclastic flow deposits 2590 years ago. Historical eruptions, which were mainly phreatic, occurred in 1796, 1817, 1917, 1936, 1952, 1993, 1994, 1999 and 2000 (Volcanic Survey of Indonesia: http://www.vsi.esdm.go.id/volcanoes/ijen/history.html).
|
|
SAMPLE GROUPINGS |
|---|
|
TOPABSTRACTINTRODUCTION AND SETTINGOVERVIEW OF THE IJEN...SAMPLE GROUPINGSANALYTICAL TECHNIQUESGEOCHEMISTRYDISCUSSIONSUMMARY AND CONCLUSIONSSUPPLEMENTARY DATAREFERENCES |
|---|
This study builds on the established stratigraphy and preliminary geochemical work of Sitorus (1990
). Few samples were available from the pre-caldera and caldera groups, therefore emphasis is on the geochemistry of the post-caldera rocks. Based on temporal variations in geochemistry, particularly in SiO2, FeO, CaO, Ni and Sr contents, Sitorus (1990) divided the post-caldera samples into six groups, but admitted that the variations are not entirely consistent through time, either within or between groups. Therefore, in this study the majority of the samples have been simply divided into CR and IC groups, based on their spatial distribution rather than age (Fig. 2). The centres of Kawah Ijen and Blau are considered separately from the two main post-caldera groups. Kawah Ijen is located on the western flank of Merapi (a CR volcano) in line with the other intra-caldera centres and appears to have geochemical affinities to both groups. Blau, although easily identifiable as a volcanic centre inside the caldera, is located north of the linear trend of the other intra-caldera volcanoes and geochemical data show scatter between the two groups. Stratigraphic relationships, geomorphology (deep incised valleys) and K–Ar dating suggest that Blau is significantly older than the rest of the intra-caldera centres (Sitorus, 1990).
A total of 43 samples from the IVC were obtained from the collection of Sitorus (1990) as rock powders, with the exception of five samples, for which hand specimens were also obtained. Major element data, mineralogical data and petrographic descriptions of the IVC rocks are taken from Sitorus (1990) and presented in Electronic Appendices 1–3, respectively (which may be downloaded from the Journal of Petrology website at http://petrology.oxfordjournals.org).
|
ANALYTICAL TECHNIQUES |
|---|
|
TOPABSTRACTINTRODUCTION AND SETTINGOVERVIEW OF THE IJEN...SAMPLE GROUPINGSANALYTICAL TECHNIQUESGEOCHEMISTRYDISCUSSIONSUMMARY AND CONCLUSIONSSUPPLEMENTARY DATAREFERENCES |
|---|
Five bulk-rock and 43 powdered samples were available for trace element and isotopic study. The generally low range of loss on ignition values (0·39–1·25 wt %) for the post-caldera rocks with the exception of three samples (2·24, 2·26 and 4·39 wt % LOI) indicates that most of the samples are fresh. Four of the five caldera samples with reported major element data (Electronic Appendix 1) have high LOI values and were not considered suitable for trace element and isotopic analysis.
Trace element concentrations in 43 IVC rocks were determined using a PerkinElmer ELAN 6000 quadrupole inductively coupled plasma mass spectrometry (ICP-MS) system at Durham University. Full details of the analytical procedure and instrument operating conditions have been given by Ottley et al. (2003) and Handley (2006). Calibration of the ELAN was achieved via the use of in-house standards and international reference materials [W2, BHVO-1, AGV1, BE-N and BIR1, using the accepted standard values given by Potts et al. (1992)], together with procedural blanks (three per batch). Internal and external reproducibility of standard values are generally less than 3% relative standard deviation (Ottley et al., 2003). Comparison of Zr and Ti elemental concentrations in IVC rocks obtained by X-ray fluorescence (Sitorus, 1990) and ICP-MS (this study) show excellent agreement, with mean relative differences of 6% and 7%, respectively.
Sr, Nd and Hf isotope analysis was undertaken at the Arthur Holmes Isotope Geology Laboratory (AHIGL), Durham University. The sample dissolution procedure and separation of Hf and Nd from rock samples is based on that presented by Dowall et al. (2003). The chemical separation procedure for Sr in this study follows that outlined by Charlier et al. (2006). Sr, Nd and Hf fractions were measured for isotope ratios using the AHIGL ThermoElectron Neptune multi-collector ICP-MS system. Detailed instrument operating conditions (including cup configurations and interference corrections) have been given by Dowall et al. (2003) and Nowell et al. (2003). Samples were analysed on several separate occasions. Data quality was monitored over the separate sessions by frequent analysis of standard reference materials throughout each run. Measured values for the NBS 987, J&M and JMC 475 standards ± 2SD error over the period of study were 0·710262 ± 11 (n = 53), 0·511106 ± 10 (n = 74) and 0·282146 ± 4 (n = 51), respectively. Data are plotted relative to NBS 987, J&M and JMC 475 standard values of 0·71024 (Thirlwall, 1991), 0·511110 (Royse et al., 1998) and 0·282160 (Nowell et al., 1998), respectively. Total procedural blanks for Sr, Nd and Hf were determined by ICP-MS on the PerkinElmer ELAN 6000 quadrupole ICP-MS system at Durham University and were below 1·2 ng for Sr, 219 pg for Nd and 73 pg for Hf. These values are considered insignificant in relation to the quantity of Sr, Nd and Hf typically processed from IVC rocks (12 µg, 3 µg and 500 ng, respectively).
Oxygen isotope analyses of mineral separates were performed by laser-fluorination at Royal Holloway, University of London using the procedures outlined by Macpherson et al. (2000). Oxygen results are reported as per mil deviations relative to the standard mean ocean water (SMOW) standard. In-house standard values of SC olivine 2 and GMG II over the period of study were +5·24
± 0·04 (1
, n = 6) and +5·69 ± 0·08 (1, n = 15); within 0·01 of accepted values. Precision and accuracy within runs were similar to the overall standard results. Oxygen yields were greater than 98% for olivine, 93% for plagioclase and between 92 and 94% for clinopyroxene. Although the yields for plagioclase and clinopyroxene are slightly low, they were consistent over the three sessions. Replicate analyses of olivine, plagioclase and clinopyroxene were within 0·06, 0·09 and 0·04, respectively.
|
GEOCHEMISTRY |
|---|
|
TOPABSTRACTINTRODUCTION AND SETTINGOVERVIEW OF THE IJEN...SAMPLE GROUPINGSANALYTICAL TECHNIQUESGEOCHEMISTRYDISCUSSIONSUMMARY AND CONCLUSIONSSUPPLEMENTARY DATAREFERENCES |
|---|
Major element data
Major element data are reported in Electronic Appendix 1 and illustrated in Fig. 3. Although SiO2 is commonly used as the index of differentiation for typically evolved arc lavas (Davidson et al., 2005
), MgO was chosen here as (1) IVC includes a large number of relatively less evolved (SiO2 < 55 wt %) rocks, and (2) there appears to be some division in SiO2 contents between the CR and IC post-caldera groups at lower MgO content (Fig. 3a). Lavas from IVC display a broad range in SiO2 content from 48 to 63 wt % and relatively low MgO contents (< 5·8 wt %) typical of island arc volcanoes, which suggests that no true, primary basalts were erupted within the complex. MgO correlates positively with Fe2O3, TiO2 and CaO, whereas Na2O, SiO2 and K2O generally increase with decreasing MgO. Al2O3 and P2O5 data are considerably more scattered, especially at lower MgO contents (Fig. 3).
|
The IC samples form a coherent, linear array in the majority of variation diagrams, with a narrower range of SiO2 and MgO than for the CR volcanic rocks. IC and CR rocks generally appear to follow different trends on plots of SiO2, Al2O3, K2O and CaO against MgO (Fig. 3a–d), with the exception of four CR samples (KI 202, KI 136, KI 164 and KI 108) that follow the IC trend. The IC group has higher SiO2 and K2O contents than the CR group at similar MgO and displays strong negative correlations, whereas within the CR group there is a weaker and less obviously negative correlation between SiO2 and K2O and MgO (Fig. 3a and 3c). The more differentiated CR rocks (below
4 wt % MgO) have elevated CaO and Al2O3 contents when compared with the IC samples at the same degree of differentiation (Fig. 3b and 3d). Samples from Kawah Ijen and Blau appear to straddle the two trends in the majority of the variation diagrams (e.g. Fig. 3b–d).
Trace element geochemistry
Ni (<40 ppm) and Cr (<60 ppm) concentrations are low in IVC rocks (Table 1) and Ni is positively correlated with MgO (e.g. Fig. 4). There is a large decrease in Ni (and Cr) abundance from the most mafic rocks to those with intermediate MgO contents, whereas the more evolved lavas display more constant concentrations. Large ion lithophile element (LILE; Rb, Cs, Ba and including U and Th) concentrations (e.g. Fig. 4b and 4c) generally increase with decreasing MgO. Variation of Rb with MgO mirrors that of K2O (Fig. 3c), where IC group contents are higher than the relatively constant concentrations in the CR group. Sr concentrations, on the other hand, are more elevated in the CR group and generally increase with decreasing MgO, whereas Sr concentrations in the IC rocks are in the main lower at similar MgO (Fig. 4d); Al2O3 and CaO exhibit similar behaviour. The IC samples have higher abundances of high field strength elements (HFSE), such as Zr, Hf, Nb and Ta, than the CR rocks at a given MgO (Fig. 4e and 4f). The variation of these elements also mirrors that of K2O. Rare earth element (REE) variation with MgO is illustrated in Fig. 4g and 4h. Light REE (LREE, e.g. La) concentrations in the IVC rocks generally increase with decreasing MgO (Fig. 4g); however, heavy REE (HREE) variation trends are more difficult to discern; there is a wide variation in Yb at low MgO content (Fig. 4h).
|
|
Mid-ocean ridge basalt (MORB)-normalized trace element abundance patterns (Fig. 5) are typical of island arc volcanoes, with enrichment of the more mobile LILE and LREE relative to HFSE and HREE. The highly to moderately incompatible elements are all enriched relative to MORB, similar to those from other Sunda arc volcanoes (Turner & Foden, 2001
) and local sediments (Vroon et al., 1995). The IC group are generally more enriched than the CR group. This enrichment is emphasized in the inset diagram of Fig. 5, where the IC are normalized to the CR group. There is a general decrease in the degree of enrichment of the IC as element incompatibility decreases, with the exception of Ba and Sr, which are considerably less enriched relative to their neighbouring elements, and in the case of Sr, which is actually depleted relative to the CR group. There is greater enrichment of Zr and Hf in the IC group compared with the more compatible elements, whereas the Ti concentrations of the IC and CR rocks are identical.
|
Radiogenic isotopes
The Sr–Nd–Hf isotope compositions of the IVC samples are listed in Table 1 and plotted in Fig. 6. In contrast to their major and trace element compositions, the intra-caldera and caldera-rim groups cannot be distinguished from one another on the basis of radiogenic isotopes. The relatively large dataset displays a restricted range in 87Sr/86Sr (0·704169–0·704483), 143Nd/144Nd (0·512814–0·512895) and 176Hf/177Hf (0·283078–0·283133), and plots within the range of other Indonesian island arc volcanoes. The Sr and Nd isotope ratios of lavas from the complex are among the least radiogenic Sr and most radiogenic Nd isotopic values reported for Javanese volcanoes (Fig. 6a); only samples from Guntur (Edwards, 1990
) and Galunggung (Gerbe et al., 1992) have higher 143Nd/144Nd and lower 87Sr/86Sr ratios. The Sr isotope ratios of the IVC lavas are consistent with the previously recognized eastward decrease in 87Sr/86Sr from West Java to Bali (Whitford, 1975). In Nd–Hf isotope space the IVC lavas lie on the edge of the Java field close to the Indian Ocean MORB field (Fig. 6b). The pre-caldera group has relatively low 87Sr/86Sr ratios (0·704169–0·704238) and are indistinguishable from one another in 143Nd/144Nd (0·512892–0·512895). However, their range in 176Hf/177Hf (0·283081–0·283118) is comparable with the range of Hf isotope ratios seen in the volcanic complex as a whole (Fig. 6b inset). IC 143Nd/144Nd values are similar to those of the CR group with no obvious division between the two groups. The IC lavas display a more limited range in 176Hf/177Hf ratios, but plot in the centre of the larger spread of the CR group (Fig. 6b inset).
|
Oxygen isotope data
The
18O values of clinopyroxene, olivine and plagioclase are +5·38 to +5·58 (n = 5), +5·02 (n = 1) and +6·08 (n = 1), respectively (Table 2). The range of clinopyroxene 18O values is narrow and comparable with the average mantle 18O value of +5·57 ± 0·32 (Ionov et al., 1994; Mattey et al., 1994). IVC 18O clinopyroxene values are lower than those of plagioclase and higher than those of olivine; with coexisting mineral pair 
cpx–ol and plag–cpx values of 0·36 and 0·5, respectively, suggesting isotopic equilibrium at typical magmatic temperatures for andesite liquids (Macpherson & Mattey, 1998; Macpherson et al., 1998). The limited dataset displays no noticeable difference in the cpx 18O values of the different geographical groups.
|
The oxygen isotope ratios of clinopyroxene phenocrysts from IVC rocks are within the range reported for clinopyroxene from the primitive Galunggung lavas from West Java (+5·3 to +5·6
, Harmon & Gerbe, 1992), and clinopyroxene and olivine 18O values lie at the lower end of the Banda arc 18O range (+5·18 to +7·04 and +4·92 to +5·59, respectively; Vroon et al., 2001). The 18O value for the only plagioclase sample analysed from the IVC is within error of the 18O range of Galunggung lavas (+5·6 to +6·0; Harmon & Gerbe, 1992) and is lower than 18O values recorded in plagioclase phenocrysts from Merapi volcanic rocks (+6·5 to +7·00, Gertisser & Keller, 2003). The 18O values of the IVC rocks also are close to the values postulated for the upper mantle in other subduction zones (Smith et al., 1996; Thirlwall et al., 1996; Macpherson & Mattey, 1998; Macpherson et al., 1998, 2000; Eiler et al., 2000). |
DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTION AND SETTINGOVERVIEW OF THE IJEN...SAMPLE GROUPINGSANALYTICAL TECHNIQUESGEOCHEMISTRYDISCUSSIONSUMMARY AND CONCLUSIONSSUPPLEMENTARY DATAREFERENCES |
|---|
Magmatic differentiation processes
Differentiation processes can significantly change the composition of magmas as they rise through the lithosphere towards the surface. It is important to identify the impact of these processes on the composition of the magma to remove uncertainty in establishing the composition and magnitude of slab-derived contributions to the mantle source in subduction zones. Low MgO, Ni and Cr abundances in IVC volcanic rocks indicate that they are not primary mantle melts, and that magma compositions were, therefore, modified en route to the Earth's surface.
Correlations of various major and trace elements with indices of differentiation (i.e. MgO) in IVC lavas (Figs 3 and 4) suggest that the concentrations of some elements are controlled by shallow-level differentiation processes such as fractional crystallization, magma mixing or contamination. Major and trace element variations for most of the post-caldera samples can be described as one of two trends, which are exemplified in the CaO vs MgO diagram (Fig. 3d), where a low-Ca group, dominantly consisting of the IC samples, displays a positive correlation, and a high-Ca group, composed of most CR samples along with some Kawah Ijen and Blau samples, shows little change in CaO with decreasing MgO. Similar trends are apparent in plots of Al2O3 and Sr against MgO (Figs 3b and 4d). Higher modal plagioclase abundances (10%, Table 3), and a less pronounced negative Eu anomaly in high-Ca lavas (Fig. 7) suggest that the different geochemical trends exhibited by the low-Ca and high-Ca groups may be due to more extensive plagioclase fractionation in the former, or plagioclase accumulation in the latter. Possible evidence for open-system processes in the petroenesis of the IVC lavas includes: coexistence of phenocrysts (plagioclase and clinopyroxene) that display normal, oscillatory and reverse zoning (Sitorus, 1990); clinopyroxene overgrowths on orthopyroxene phenocrysts; and bimodal distributions in plagioclase core compositions (Sitorus, 1990, see Electronic Appendix 2; Berlo 2001). One group of plagioclase crystal cores cluster below An60, and the other generally above An80. However, bimodal plagioclase compositions are found within volcanic rocks from both the high- and low-Ca groups, indicating that although mixing may occur, it is not the dominant control on chemical variation, otherwise we would only expect to find bimodal plagioclase populations in high-Ca rocks. Also, higher modal plagioclase abundance (Table 3; Electronic Appendix 3) does not necessarily indicate accumulation; it may just be a result of a greater degree of crystallization or more plagioclase in the modal equilibrium assemblage. Variations in the fractionating mineral assemblage can produce contrasting differentiation trends on variation diagrams (e.g. Davidson, 1996). Crystallization vectors drawn in Fig. 3b and 3d illustrate the magmatic evolution predicted for primitive IVC magma during fractionation of different mineral assemblages. The different trends observed in the high- and low-Ca groups can result from differences in the phases that make up the fractionating assemblage and their proportions. Assemblages consisting of clinopyroxene and olivine ± plagioclase can replicate the high-Ca array (Fig. 3b and 3d), whereas the same assemblage but with more plagioclase fractionation relative to clinopyroxene and olivine is suggested for the low-Ca trend (see labelled arrows, Fig. 3b and 3d). Less plagioclase fractionation in the evolution of the high-Ca lavas, relative to the low-Ca lavas, is also consistent with the variations of Al2O3, Sr and Ba with respect to MgO, along with a smaller negative Eu anomaly.
|
|
Fractional crystallization
The XLFRAC, least-squares major element modelling technique (see Stormer & Nicholls, 1978
) was used to ascertain whether the high-Ca and low-Ca differentiation trends can be explained by differences in the amount of plagioclase fractionation, and also to determine whether compositional variation within and between geographical groups (e.g. IC, CR, Kawah Ijen) can be explained by fractional crystallization. Models 1 and 2 in Table 3 suggest that the most evolved post-caldera rocks of the high-Ca and low-Ca groups can be produced by fractionation of plagioclase, clinopyroxene, olivine and Fe–Ti oxide from the least evolved samples of their respective groups (
r2 = 0·2 and 0·19, respectively). Models of fractionation to other high-Ca and low-Ca group daughter compositions (models 3–6 and 7–9, respectively) yield excellent r2 values (<0·16). All of these models also suggest that significantly less plagioclase fractionation is required from the parent magma to reach high-Ca daughter compositions (10–12%) compared with those with low-Ca (20–37%). Therefore, we suggest that the different (high- and low-Ca) trends (e.g. Fig. 3b and 3d) are due to a greater amount of plagioclase fractionation in the IC rocks than in the CR rocks.
Least-squares analysis demonstrates that intra-group CR fractionation models yield low r2 values (models 2–9, Table 3). A good solution (r2 = 0·02) is also obtained for a model using a less evolved IC sample as parent to a more evolved IC rock (model 10). Low r2 values are obtained in model 11, between the two Kawah Ijen samples, which straddle the high-Ca and low-Ca trends on several major element diagrams. Some inter-group models can also generate very acceptable results (models 13–15), suggesting that members of the different geographical groups can be related to each other through fractional crystallization. However, least-squares modelling cannot generate low r2 in models where IC rocks represent the parent magma composition and CR is the daughter, unless crystals are accumulated rather than removed (models 16 and 17).
Utilizing the phase proportions and degree of crystallization predicted from the major element modelling (Table 3) it is possible to test the conclusions of least-squares analysis by modelling trace element concentrations using the Rayleigh fractionation equation Cl = CoF(D–1), where Cl and Co represent the concentration of an element in the daughter and parental liquids, respectively, F is the fraction of liquid remaining and D is the bulk distribution coefficient. The distribution coefficients used in modelling are given in Table 4. The results of selected trace element models are given in Table 5 and show excellent agreement between calculated and observed daughter concentrations. The strongest agreement (model 10) is between the intra-caldera samples KI 92 and KI 34 in the low-Ca group, where calculated values lie within 10% of the measured concentrations. In summary, trace element modelling validates models of fractional crystallization developed from major element data; this suggests that fractional crystallization is the main control on melt evolution at IVC. It also confirms that the separate trends of the low-Ca and high-Ca groups are likely to be due to variable influence of plagioclase in the fractionating mineral assemblage and are not a result of plagioclase accumulation in the high-Ca group.
|
|
Structural controls on magma ascent and storage
The location of the CR volcanoes around the caldera rim suggests that ring fractures may facilitate the movement of magma below these volcanoes. This has been proposed for other post-caldera vents located along radial fractures within volcanic complexes, such as Roccamonfina volcano in Italy (Giannetti, 2001
) and Chichontepec volcanic centre in El Salvador (Rotolo & Castorina, 1998). The linear orientation (roughly NE–SW, Berlo, 2001) of the intra-caldera volcanoes also suggests that magma below the central part of the caldera is utilizing lines of weakness and further highlights the structural control on the location of post-caldera volcanoes at IVC.
The different fractionation trends exhibited by the low-Ca (dominantly the intra-caldera rocks) and high-Ca (dominantly the caldera-rim rocks) suggest that spatial variations in chemistry within the volcanic complex might be linked to sub-volcanic structure. The contrast in geochemistry is proposed to result from differences in the amount of plagioclase fractionation from the respective magmas. Experimental evidence (Grove et al., 2003) suggests that plagioclase crystallization is suppressed by high water contents and high pressure in basaltic andesite magmas. Therefore, the more extensive plagioclase fractionation inferred for the IC group magmas could result from either shallower level (i.e. lower pressure) storage or lower water contents (or both) in the IC magmas relative to the CR magmas. The volatile content of magmas is also controlled by pressure; thus, the suppression of plagioclase crystallization in the CR group could be a result of magmatic differentiation at deeper levels in the crust below the caldera-rim volcanoes (Fig. 8). Magmas stored at greater depth, with higher volatile contents, might be expected to erupt more explosively than those stored at shallower levels, from which eruptions might be more passive. The presence of large stratovolcanoes on the caldera-rim and dominantly small cinder cones inside the caldera at IVC is consistent with the type of eruption, compatible with this interpretation. Therefore, at IVC it is proposed that volcanic structure exerts some control on the depths at which the rising magma can pond. Kawah Ijen erupts lavas belonging to both the low-Ca and high-Ca trends. Kawah Ijen lies at the intersection between the caldera rim and the northeastern edge of the intra-caldera lineation (Fig. 2) and could, therefore, tap both shallow and deeply ponded magmas (Fig. 8). There is no temporal division in the production of high-Ca and low-Ca lava types at Kawah Ijen, and no evidence for physical mingling or mixing (Berlo, 2001), implying separate magmatic pathways for magmas erupted from different depths (Fig. 8).
|
Role of crustal contamination
There is considerable evidence from other island arcs (Thirlwall & Graham, 1984
; Davidson et al. 1987; Ellam & Harmon, 1990; Davidson, 1996; Smith et al., 1996; Thirlwall et al., 1996; Macpherson et al., 1998) for the contamination of primary magmas by the arc crust. Contamination is also believed to be an important process, responsible for modifying isotope ratios, in the western Sunda arc (Gasparon et al., 1994; Gasparon & Varne, 1998) and at Sangeang Api volcano in the East Sunda arc (Turner et al., 2003). Therefore, prior to discussing magma source compositions it is important to assess the role of crustal contamination at IVC. The absence of correlations between IVC Sr, Nd and Hf isotope ratios and indices of differentiation (e.g. MgO, Fig. 9) and the restricted ranges of Sr, Nd and Hf isotope ratios in IVC samples (Table 1) are consistent with a negligible input of isotopically distinct crust during differentiation. Clinopyroxene 18O values in the IVC lavas are homogeneous (+5·38 to +5·58, Table 2) and low, lying within the range of mantle 18O values (+5·57 ± 0·32) reported by Mattey et al. (1994) and Ionov et al. (1994). Therefore, it is unlikely that the lavas have been contaminated by upper crustal materials, which typically possess high 18O values, such as those reported by Gertisser & Keller (2003) for the local upper crust (calcareous sediments) in Central Java (+20·5 and +18·9). Furthermore, only igneous cumulate xenoliths have been found within the Ijen complex lavas; continental-type crustal xenoliths have not been detected (Sitorus, 1990; Berlo, 2001). Monomineralic plagioclase xenoliths are present in lavas erupted from Anyar (IC) and Blau, and cumulate xenoliths containing plagioclase, olivine, clinopyroxene and Ti-magnetite have been found in Rante (CR) and Kawah Ijen lavas. A cumulate containing olivine, ortho- and clinopyroxene was also found in lava from Kukusan (Berlo, 2001).
|
Similar to conclusions reached in studies of other Sunda arc volcanoes (Gerbe et al., 1992
; Elburg et al., 2002; Gertisser & Keller, 2003), the evidence above suggests that any interaction of IVC magmas with the arc crust during differentiation has had a negligible impact on their geochemistry. Fractional crystallization, therefore, appears to be the dominant differentiation process controlling geochemical variations in IVC lavas. Higher Sr and lower Nd and Hf isotope ratios in IVC rocks relative to MORB point towards contamination of the IVC source by an isotopically distinct crustal component.
Magma source components
Most models of magma petrogenesis at island arcs involve three main source components: (1) the mantle wedge; (2) the subducting slab (oceanic crust and associated sediments); (3) the arc lithosphere. The majority of island arc magmas are thought to originate in the mantle wedge (Ringwood, 1974; Ellam & Hawkesworth, 1988; McCulloch & Gamble, 1991), which has been inferred by several workers to be similar to the source of MORB (Davidson, 1987; Woodhead et al., 1993; Gamble et al., 1996; Turner et al., 2003). It has been proposed that there might be contribution from an enriched mantle component [ocean island basalt (OIB) source] in some Sunda arc lavas (Wheller et al., 1987; Edwards et al., 1991, 1993; van Bergen et al., 1992), although a MORB-source-like mantle source has been advocated for the Sunda arc by others (e.g. White & Patchett, 1984; Turner & Foden, 2001; Elburg et al., 2002). Helium isotope values of olivine crystals in mantle xenoliths and island arc volcanic rocks also implicate a MORB-source-like mantle source in the western Sunda arc (Hilton & Craig, 1989; Gasparon et al., 1994).
Mantle source characteristics
The HREE and HFSE are thought to remain relatively immobile compared with other elements, such as LILE, during slab dehydration (Tatsumi et al., 1986; Kessel et al., 2005) and abundances are generally too low in oceanic and continental crust (Taylor & McLennan, 1985) to significantly alter ratios of these elements in magmas derived from the mantle. We can use these trace elements, therefore, to help ascertain the pre-subduction composition of the mantle wedge. HREE concentrations are 10–15 times chondrite values in IVC lavas and display relatively flat profiles (Fig. 10). This suggests that garnet is not an important residual mineral in the source region and that the IVC magmas are derived largely from a shallow mantle source, above the garnet–spinel transition for wet peridotite. The La/Lu ratios of CR (41) and IC (44) lavas are very similar, which is consistent with a similar source mineralogy and degree of partial melting for these two post-caldera groups. The large difference in La/Lu ratios within the three pre-caldera group samples (34–80), when compared with the fairly constant ratio of the post-caldera volcanics (42) might be due to variable degrees of partial melting and/or different source mineralogies for the pre-caldera rocks. If post-caldera rocks from the IC and CR volcanoes, Kawah Ijen and Blau share the same source then they should possess similar ratios of immobile trace elements generally assumed to be unmodified by subduction processes (Table 6). Zr/Nb [17–23 (except one CR sample at 14)] and Ta/Nb (0·07–0·13) ratios do not change significantly with differentiation (Fig. 11a and 11b) and are relatively homogeneous in all IVC eruptive rocks. Zr/Nb ratios are similar to those of MORB [normal and Indian type (N-MORB and I-MORB)] whereas Ta/Nb ratios are comparable with both MORB and OIB (Table 6; e.g. Fig. 11). The similarity between HFSE and HREE concentrations in primitive IVC basalts and I-MORB, rather than N-MORB, is also apparent in Fig. 12. Therefore, we suggest that the mantle wedge beneath the IVC is similar to the source of I-MORB. Several studies utilizing Pb isotopic and trace element data have reached a similar conclusion for magmatism in SE Asia throughout the Cenozoic (Taylor et al., 1994; Hickey-Vargas, 1998; Macpherson & Hall, 2001, 2002; Elburg et al., 2002; Macpherson et al., 2003). The HFSE and HREE concentrations and ratios of IVC lavas also indicate that the mantle wedge is not significantly depleted beneath East Java, in contrast to others arcs such as Izu–Bonin and Mariana (Fig. 12), where mantle sources are thought to have experienced melt extraction prior to their involvement in arc petrogenesis (Woodhead et al., 1993; Elliott et al., 1997; Taylor & Nesbitt, 1998).
|
|
|
|
Sr is fluid-mobile during slab dehydration, therefore, Sr isotope data are not well suited for identifying the isotopic composition of the precursor mantle wedge. Experimental work (Tatsumi et al., 1986
; Brenan et al., 1995; You et al., 1996) and studies of arc lavas (McCulloch & Gamble, 1991; Pearce & Peate, 1995; Münker et al., 2004) indicate that Hf and Nd are relatively immobile during slab dehydration, particularly with respect to the formation of aqueous fluids (Kessel et al., 2005). However, Woodhead et al. (2001) have cast doubt upon the status of Hf as a truly fluid-immobile element; they see evidence in the New Britain subduction system for the transport of Hf in aqueous fluids when Nd is immobile. Their observations are based on the contrast in Hf isotope ratios of the mantle wedge beneath New Britain and the subducting crust of the Woodlark Basin. As we advocate an I-MORB composition for both the mantle wedge and the down-going subducted crust, any Hf mobility should have negligible impact on the Hf isotope ratio of the IVC magma source. Therefore, combined Nd and Hf isotope data can better document the source (mantle wedge) characteristics of the IVC lavas than combinations of isotope ratios involving 87Sr/86Sr isotope data (e.g. Pearce et al., 1999). IVC lavas display a restricted range in 143Nd/144Nd and especially 176Hf/177Hf and lie within the Hf isotope range of I-MORB. Lavas from other Javan volcanoes diverge further from typical mantle (MORB) values towards more ‘crustal’ Hf–Nd isotopic compositions. Therefore, although the major element data indicate that the IVC lavas are moderately fractionated and are not primary magmas, the IVC source appears to be among the least affected by crustal inputs from either the slab or arc lithosphere (noting that the Hf isotope database is slightly more limited in comparison with Sr–Nd isotope data available). Volcanic rocks from the IVC and other volcanic centres on Java lie above the dividing line that separates MORB of Indian and Pacific provenance drawn in Hf–Nd isotope space by Pearce et al. (1999) (Fig. 6b), implying derivation from an I-MORB-like mantle source rather than a source similar to Pacific MORB, and corroborating conclusions from ratios of immobile trace elements (above). In summary, restricted variation in isotope ratios and immobile trace element ratios suggests that the mantle sources of the different groups of IVC magmatic rocks are not distinct from each other. The mantle wedge beneath East Java is relatively fertile compared with that of many other island arcs, with HFSE (especially Nb and Ta) and HREE concentrations similar to those in the source of I-MORB.
Identification of subduction zone components: trace element constraints
Enrichment of LILE and LREE relative to HFSE in arc volcanic rocks has long been attributed to slab involvement (Hawkesworth et al., 1979; Kay, 1980; White & Dupré, 1986; Davidson, 1987; McCulloch & Gamble, 1991), but there remains considerable uncertainty as to the source of this contribution (altered oceanic crust or sediments or both) and the transfer mechanism of this component to the sub-arc mantle (as a fluid or melt or both).
Th and the LREE are thought to be less mobile in aqueous fluids than the LILE (Tatsumi et al., 1986; Pearce et al., 1999) and, consequently, high Ce/Pb and Th/Yb ratios in arc lavas are taken to indicate a sedimentary contribution from the slab (Woodhead et al., 2001; Gertisser & Keller, 2003). The large contrast in Ce/Pb and Th/Yb ratios between MORB (N- and I-MORB) and Indian Ocean sediments near Java (Table 6; Fig. 13a), suggests that relatively low Ce/Pb (3–7) and high Th/Yb (1–4) ratios of the IVC could be attributed to a sediment component in the arc mantle source. However, the slightly higher Th/Yb ratios observed in lavas of the low-Ca group suggest that this ratio was also affected by fractional crystallization. A plot of Th/Yb vs MgO (Fig. 13b) reveals slight elevation and correlation of Th/Yb ratios in the IC group compared with the CR group. The CR group displays limited variation (Th/Yb 1·5) over the entire range of MgO, and we infer that this ratio is representative of the IVC source.
|
|
Elevations of LILE/LREE ratios in arc lavas compared with MORB are thought to indicate fluid addition to the mantle wedge from dehydration of the subducting slab (Ben Othman et al., 1989
; Elliott et al., 1997). The Ba/La (23–36) and Sr/Nd (16–44) ratios of the IVC lavas are slightly to significantly higher than those of Indian Ocean sediment and also much greater than I-MORB (Table 6; Fig. 13a), indicating that slab fluids may also be added to the IVC mantle source, possibly from the altered oceanic crust (AOC). The Ba/La ratios of the IVC volcanic rocks are not as high as those observed in the Kermadec, Mariana and New Britain arcs, where the dominant slab input is thought to be a hydrous fluid (Woodhead et al., 2001). However, the overall flux of hydrous fluid might be the same in all these arcs; the lower Ba/La in the IVC lavas could reflect the addition of La from a subducted sedimentary component.
Modelling slab input: constraints from radiogenic isotopes
Two-component model: I-MORB + bulk sediment. Figure 14 shows simple mixing between an I-MORB source and bulk sediment for a variety of Indian Ocean sediment end-members (Table 7). Mixing arrays involving three of the four end-member sediments (A–C; Mn nodule, nanno-ooze and average Java sediment, respectively) produce reasonable fits to the general range of Java volcanic rocks. Addition of 0·5–1% average bulk Java sediment (C; Plank & Langmuir, 1998) to I-MORB produces a curve that passes through the most isotopically primitive IVC volcanic rocks (Fig. 14), but that is highly oblique to the IVC array. The terrigenous–biogenic sediment end-member (D) has high Sr/Nd and produces a convex-up mixing trend in Sr–Nd isotope space when added to I-MORB mantle (Fig. 14). Bulk terrigenous–biogenic sediments are unlikely to dominate the sediment contribution to the mantle source of the IVC volcanic rocks. This is supported by sediment data collected from the Indian Ocean. Terrigenous sands and silts, such as turbidites, are present in the average composition of the sedimentary column subducted at the Java trench as described by Plank & Langmuir (1993), but are minor components compared with clay and siliceous ooze. Boreholes 211 and 261 of Deep Sea Drilling Project (DSDP) leg 22 (SEATAR, 1981: http://www.ngdc.noaa.gov/mgg/geology/seatar.html) located in the Indian Ocean crust south of Java are dominated by nanno-ooze, nanno-claystone, a small amount of volcanic ash and a 12 cm thick surface layer of small manganese nodules.
|
|
|
Two-stage, three-component model: I-MORB + AOC component + bulk sediment. Two-component mixing fails to provide a good fit to the isotopic variation of the IVC lavas, suggesting that an additional isotopically distinct component is required in the source. Furthermore, some trace element ratios (e.g. Ba/La) are too high for MORB-source mantle and sediment mixtures to explain, but are consistent with a fluid contribution derived from the subducted AOC. Contributions from AOC have been inferred for many Sunda arc lavas (Turner & Foden, 2001
) and there is considerable support for two separate contributions from the subducting slab to the mantle wedge in the petrogenesis of many other arc lavas (White & Patchett, 1984; Ellam & Hawkesworth, 1988; Thirlwall et al., 1996; Elliott et al., 1997; Hawkesworth et al., 1997; Turner & Hawkesworth 1997; Turner et al., 1997; Taylor & Nesbitt, 1998). Interaction of seawater with basalt produces significant chemical and isotopic changes (Hart et al., 1974; White & Patchett, 1984). Altered oceanic crust displays a wide range in 87Sr/86Sr ratios whereas 143Nd/144Nd isotope ratios are generally indistinguishable from those of fresh MORB (Fig. 15). Because of the lack of isotopic data for altered Indian Ocean MORB we have used data for Atlantic MORB (Staudigel et al., 1996) to consider the effect of adding this component to the mantle wedge. Regardless of how this component is added to the mantle wedge, either as a fluid (most commonly proposed; e.g. Ellam & Hawkesworth, 1988; Turner et al., 1997) or melt (particularly where the subducting crust is young; e.g. Defant & Drummond, 1990), the isotopic composition of the resultant mixture will be generally the same; mixing I-MORB source with AOC fluid or melt will create a mixture with higher Sr isotope ratios and similar Nd isotope ratios to I-MORB. However, the quantity of AOC component required will be dependent on the transfer mechanism as a result of higher Sr/Nd ratios in fluids produced from the AOC compared with melts. The trajectory of the IVC data crosscuts any conceivable two-component mixing curve produced between the mantle source and an AOC–sediment mixture; therefore a two-stage model is preferred. In the mixing models, we have chosen to add the AOC component as a fluid, as it is hard to envisage (in a two-stage process whereby the mantle source is first contaminated by AOC prior to bulk addition of sediment; see Elliott et al., 1997; Turner & Hawkesworth, 1997) melting of the AOC and addition of this to the mantle wedge without melting the sediments and incorporating them too, simultaneously, to the mantle wedge. Furthermore, the IVC lavas exhibit none of the geochemical traits thought to record slab melting (Defant & Drummond, 1990). Figure 15 shows the results obtained from models of two-stage source contamination, where an I-MORB source is first contaminated by fluid from the AOC and then bulk sediment [Mn nodule (curve A-2) and nanno-ooze (curve B-2)]. These types of sediments are consistent with the findings from DSDP leg 22 boreholes 211 and 261 mentioned above. To create the fluid-modified mantle (FMM) component in this model, 3% AOC fluid was added (see Table 7 for AOC fluid calculation details). Mixing proportions of AOC and I-MORB are highly dependent on the Sr and Nd contents and isotopic ratios of AOC, which are themselves extremely variable. The fraction of fluid required is smaller if more highly altered basalt compositions are used, and greater if partial melting or bulk mixing of AOC is considered. The three-component models produce excellent fits not only to the trend of IVC lavas but also a large proportion of other Java volcanic rocks. Only a small percentage of sediment (0·25% Mn nodule or 0·5–1% nanno-ooze) is required to generate the isotopic values seen in the IVC volcanic rocks.
Three-component mixing in Nd and Hf isotope space. Modelling of two-stage source contamination in Nd–Hf isotope space is hampered by the paucity of Hf isotope data reported for altered oceanic crust. Nd isotope ratios remain relatively constant during hydrothermal and low-temperature seawater alteration of MORB, therefore, White & Patchett (1984) predicted that Hf isotope ratios will also be unaffected by the alteration of oceanic crust. In these circumstances, contamination of the mantle source by a slab fluid from the AOC would be undetectable using Nd and Hf isotopes.
The relatively shallow array of IVC lavas projects back towards a mantle source with lower 176Hf/177Hf ratios than average I-MORB (Fig. 16). Mixing of an I-MORB source with a lower than average Hf isotope composition (IM1, IM2 and IM3; Table 7) with sediment A (Hf isotope data are unavailable for the nanno-ooze, B) produces mixing curves that correspond perfectly to the IVC data array, requiring <0·5% addition of sediment (all mixing curves in Fig. 16) and also match the percentage of sediment required to explain the IVC data in three-component Sr–Nd mixing (Fig. 15). Using low 176Hf/177Hf I-MORB source end-members instead of average I-MORB compositions does not affect the conclusions reached from three-component modelling in Sr–Nd isotope space (Fig. 16 inset); almost identical results are obtained (compare with curve A-2 in Fig. 15). This implies that, if a hydrous fluid is added to the mantle source from the altered oceanic crust, there is a minimal effect on the resulting Nd and Hf isotope ratios.
In conclusion, radiogenic isotope and trace element data for IVC lavas suggest that a three-component, two-stage source contamination model is applicable. In this model I-MORB mantle is metasomatized by slab fluids prior to addition of subducted sediment. Three-component models have also been suggested in geochemical studies of other island arc volcanoes (Ellam & Hawkesworth, 1988; Elliott et al., 1997; Turner et al., 1997; Taylor & Nesbitt, 1998). However, the order in which the two components from the slab (fluid and sediment) are introduced to the mantle wedge beneath IVC differs from that proposed in other models. The U-series isotopic compositions of Mariana lavas indicate that addition of a sedimentary component took place at least 350 kyr before fluid addition to the mantle wedge beneath the Mariana arc (Elliott et al., 1997). The authors propose two possible models for the subduction component addition: (1) shallow (<120 km depth) sedimentary melts produce enriched sub-arc mantle, which is then dragged down to deeper depths where dehydration of the altered oceanic crust then releases fluid into the enriched sub-arc mantle causing melting; (2) altered oceanic crust dehydrates at shallow depth to produce hydrated peridotite mantle, which is down-dragged and dehydrates at around 120 km depth causing extensive melting. This material then encounters sedimentary material rising from greater depth. Elliot et al. (1997) acknowledge the inherent problems in their most simplified model (1) of satisfying thermal models of subduction zones (Peacock, 1991; Johnson & Plank, 1999). The reverse order in the addition of subduction components proposed at IVC (i.e. slab dehydration and fluid addition to the mantle wedge prior to sediment addition) appears mechanically more viable. However, in this model, fluid sourced predominantly from the oceanic crust must pass through the overlying sediment and avoid any sedimentary imprint. This is possible, as Elliott et al. (1997) pointed out that the sedimentary pile on top of the oceanic crust is likely to dehydrate and lose its fluid-mobile elements at shallower depths before major dehydration of the basaltic crust. This earlier fluid addition to the wedge may avoid being down-dragged to the region of arc magma genesis (Fryer, 1992) and, therefore, not contribute to magma geochemistry. Turner & Hawkesworth (1997) also proposed a model in which partial melts of sediment are added to the mantle wedge at relatively shallow levels prior to fluid addition from the slab. The timescale envisaged for fluid addition is of the order of 30–50 kyr before eruption.
In comparison, at volcanoes of the Andonara–Pantar section in the East Sunda arc, U-series, Sr–Nd–Pb isotope and trace element data support a similar transfer order for the addition of fluid- and sediment-sourced slab components (Hoogewerff et al., 1997) as that proposed at IVC, although differing slightly in the exact nature of the slab components (sediment-derived hydrous fluid and siliceous melt). The hydrous fluid is proposed to dominate the input in the shallow part of the subduction zone, whereas the sedimentary melt dominates the flux at deeper levels. Evidence for relatively recent fluid addition is provided by 226Ra excesses in volcanic rocks closest to the arc front in the Andonara–Pantar section and also in other Sunda arc lavas (Turner & Foden, 2001; Turner et al., 2003).
Comparison of slab contributions
The amount of local sediment incorporated in the source region of IVC lavas (<0·5 Mn nodule or 0·5–1% nanno-ooze) is similar to or less than the amount of subducted sediment proposed in magma genesis at other arcs; for example, the Lesser Antilles (White & Dupré, 1986), Banda (Vroon, 1992), Izu–Bonin (Taylor & Nesbit, 1998) and Aleutian arcs (McCulloch & Perfit, 1981). Radiogenic isotopic compositions of other Sunda arc lavas are also accounted for by incorporation of similar amounts of subducted Indian Ocean sediment into the mantle wedge to that proposed beneath IVC: 1–2% (Merapi, Gertisser & Keller, 2003), <2% (general Sunda arc lavas, Turner & Foden, 2001), <5% (eastern Sunda arc, Stoltz et al., 1990), 2·5% (Sangeang Api, East Sunda arc, Turner et al., 2003). Several workers have proposed that the sedimentary component is added to mantle source in the form of a partial melt (Turner & Foden, 2001; Turner et al., 2003). However, in eastern Sunda arc lavas calculated melt fractions of pelagic and biogenic sediments at low degrees of melting (1–10%) have low Sr/Nd ratios that produce curves that do not fit the volcanic data (even up to 40% melting) (Stoltz et al., 1990). Sediment addition as a partial melt is a possibility beneath the IVC, but is not required to explain the data. However, considering this possibility we propose that the bulk sediment to mantle- or fluid-modified mantle mixing curves provide an upper estimate of the amount of sediment contribution in the IVC magma source.
|
SUMMARY AND CONCLUSIONS |
|---|
|
TOPABSTRACTINTRODUCTION AND SETTINGOVERVIEW OF THE IJEN...SAMPLE GROUPINGSANALYTICAL TECHNIQUESGEOCHEMISTRYDISCUSSIONSUMMARY AND CONCLUSIONSSUPPLEMENTARY DATAREFERENCES |
|---|
Differentiation processes can account for the spatial variations observed in the geochemistry of volcanic rocks at the Ijen Volcanic Complex. Major element and trace element modelling suggests that trends towards high CaO, Al2O3 and Sr contents with increasing MgO in the CR group (high-Ca group) and positive correlations of these elements in the IC group (low-Ca group) are a result of slight differences in fractionating mineral assemblages: olivine, clinopyroxene ± plagioclase fractionation dominating the CR group, and plagioclase > clinopyroxene and olivine in the IC group (Fig. 8). Spatial variations in lava geochemistry are controlled to some extent by sub-volcanic structure. It is proposed that magma storage at shallower depths is prevented beneath volcanoes located on the caldera-rim (compare with the intra-caldera magmas) because of the presence of ring fractures, which focus and facilitate migration of magma towards the surface and inhibit the formation of shallow-level storage chambers. As a result, CR magmas do not undergo substantial shallow-level degassing or extensive plagioclase fractionation, like that envisaged for the IC volcanoes (Fig. 8). The linear orientation of the IC volcanoes further indicates the importance of sub-volcanic structure for volcanism at IVC. Location of Kawah Ijen at the intersection between the caldera rim and the northeastern edge of the intra-caldera lineation, combined with the eruption of lavas with the characteristics of both differentiation trends, implies that Kawah Ijen may be sourced by both the shallow and deep magma reservoirs by separate pathways (Fig. 8).
Low, mantle-like 18O values and a lack of correlation between isotope ratios and indices of differentiation imply that contamination by continental ‘crustal’ material is insignificant at IVC, consistent with IVC being located on thickened oceanic crust. Assimilation of isotopically indistinct material (i.e. older volcanic crust) cannot be discounted at IVC.
Despite the variation observed in major and trace element data, the HFSE ratios, Nd, Hf isotope ratios and 18O values of the IVC lavas are remarkably constant, suggesting that the subduction-modified mantle wedge source was relatively homogeneous. The IVC lavas have higher Hf and Nd isotope ratios than those of other Sunda arc volcanoes and lie on the edge of the I-MORB field in Hf–Nd isotope space, suggesting that they reflect the composition of the uncontaminated mantle wedge in East Java.
IVC magma genesis can be explained by two-stage, three-component source mixing. In this model, a fertile I-MORB source mantle wedge is infiltrated first by a hydrous slab fluid and then by sediment. As a result of the orientation of the IVC array and the shape of the mixing curves produced by different types of sediment addition in Hf–Nd isotope space (and constraints on sediment type from Sr–Nd isotope modelling), the I-MORB source component proposed has a lower than average 176Hf/177Hf isotope composition. In the model, the I-MORB source is contaminated by a small percentage of fluid created by the dehydration of the subducting slab (dominated by the chemistry of altered oceanic crust); this generates a mixture with higher 87Sr/86Sr, and similar 143Nd/144Nd to I-MORB. At a later stage, less than 1% of Indian Ocean sediment is transferred from the slab into the fluid-modified mantle source. The type of subducted sediment proposed in magma genesis at IVC (nanno-ooze and Mn nodule) is consistent with those found as modern deposits on the sea bed south of East Java.
|
SUPPLEMENTARY DATA |
|---|
|
TOPABSTRACTINTRODUCTION AND SETTINGOVERVIEW OF THE IJEN...SAMPLE GROUPINGSANALYTICAL TECHNIQUESGEOCHEMISTRYDISCUSSIONSUMMARY AND CONCLUSIONSSUPPLEMENTARY DATAREFERENCES |
|---|
Supplementary data for this paper are available at Journal of Petrology online.
|
ACKNOWLEDGEMENTS |
|---|
We thank John Gamble and Joel Baker for provision of IVC samples. Invaluable assistance was provided by Chris Ottley during trace element analysis and by Geoff Nowell throughout analytical work in the Arthur Holmes Isotope Geology Laboratory (AHIGL) at Durham University. We are grateful to D. Peate, I. Smith and R. Gertisser for constructive reviews, which greatly improved the manuscript. We also wish to thank Robert Hall and Helen Smyth of Royal Holloway University of London, and Val Troll and Jane Chadwick of Trinity College Dublin for their valuable comments and contributions. This work was funded by an NERC studentship (NER/S/A/2001/06127) and supported by the South East Asia Research Group (SEARG) at Royal Holloway University of London.
*Corresponding author. Present address: GEMOC Key Centre, Department of Earth and Planetary Sciences, Macquarie University, Sydney, NSW 2109, Australia. Telephone: +61 2 9850 4405. Fax: +61 2 9850 8943. E-mail: hhandley{at}els.mq.edu.au
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTION AND SETTINGOVERVIEW OF THE IJEN...SAMPLE GROUPINGSANALYTICAL TECHNIQUESGEOCHEMISTRYDISCUSSIONSUMMARY AND CONCLUSIONSSUPPLEMENTARY DATAREFERENCES |
|---|
Ben Othman DB, White WM, Patchett J. The geochemistry of marine sediments, island arc magma genesis, and crust–mantle recycling. Earth and Planetary Science Letters (1989) 94:1–21.[CrossRef][Web of Science]
Berlo K. The magmatic evolution of the Ijen Caldera. In: Unpublished MSc. Thesis (2001) East Java, Indonesia: University of Utrecht.
Brenan J, Shaw HF, Phinney DL, Ryerson JF. Mineral–aqueous fluid partitioning of trace elements at 900°C and 2·0 GPa: constraints on the trace element geochemistry of mantle and deep crustal fluids. Geochimica et Cosmochimica Acta (1995) 59:3331–3350.[CrossRef][Web of Science]
Charlier BLA, Ginibre C, Morgan D, Nowell GM, Pearson DG, Davidson JP, Ottley CJ. Methods for the microsampling and high-precision analysis of strontium and rubidium isotopes at single crystal scale for petrological and geochronological applications. Chemical Geology (2006) 232:114–133.[CrossRef][Web of Science]
Chase CG. Plate kinematics: the Americas, East Africa, and the rest of the world. Earth and Planetary Science Letters (1978) 37:355–368.[CrossRef][Web of Science]
Chauvel C, Blichert-Toft J. A hafnium isotope and trace element perspective on melting of the depleted mantle. Earth and Planetary Science Letters (2001) 190:137–151.[CrossRef][Web of Science]
Class C, Miller DM, Goldstein SL, Langmuir CH. Distinguishing melt and fluid subduction and components in Umnak volcanics, Aleutian Arc. Geochemisty, Geophysics, Geosystems (2000) 1(6). paper number 20000601.
Davidson JP. Crustal contamination versus subduction zone enrichment: examples from the Lesser Antilles and implications for mantle source compositions of island arc volcanic rocks. Geochimica et Cosmochimica Acta (1987) 51:2185–2198.[CrossRef][Web of Science]
Davidson JP. Deciphering mantle and crustal signatures in subduction zone magmatism. Subduction: Top to Bottom. Geophysical Monograph, American Geophysical Union—Bebout GE, Scholl DW, Kirby SH, Platt JP, eds. (1996) 95:251–262.
Davidson JP, Dungan MA, Ferguson KM, Colucci MT. Crust–magma interactions and the evolution of arc magmas: the San Pedro–Pellado volcanic complex, southern Chilean Andes. Geology (1987) 15:443–446.
Davidson JP, Hora JM, Garrison JM, Dungan MA. Crustal forensics in arc magmas. Journal of Volcanology and Geothermal Research (2005) 140:157–170.[CrossRef][Web of Science]
Defant MJ, Drummond MS. Derivation of some modern arc magmas by melting of young subducted lithosphere. Nature (1990) 347:662–665.[CrossRef]
DeMets C, Gordon RG, Argus DF, Stein S. Current plate motions. Geophysical Journal International (1990) 101:425–478.[Web of Science]
Dowall DP, Nowell GM, Pearson DG. Chemical pre-concentration procedures for high-precision analysis of Hf–Nd–Sr isotopes in geological materials by plasma ionisation multi-collector mass spectrometry (PIMMS) techniques. In: Plasma Source Mass Spectrometry. Special Publication, Royal Society of Chemistry—Holland G, Tanner SD, eds. (2003) 321–337.
Edwards CMH. Petrogenesis of tholeiitic, calc-alkaline and alkaline volcanic rocks. In: PhD thesis (1990) Sunda arc, Indonesia: Royal Holloway, University of London.
Edwards CMH, Menzies MA, Thirlwall MF. Evidence from Muriah, Indonesia, for the interplay of supra-subduction zone and intraplate processes in the genesis of potassic alkaline magmas. Journal of Petrology (1991) 32:555–592.
Edwards CMH, Morris JD, Thirlwall MF. Separating mantle from slab signatures in arc lavas using B/Be and radiogenic isotope systematics. Nature (1993) 362:530–533.[CrossRef]
Eiler JM, Crawford AJ, Elliott TR, Farley KA, Valley JW, Stolper EM. Oxygen isotope geochemistry of oceanic-arc lavas. Journal of Petrology (2000) 41:229–256.
Elburg M, van Bergen M, Hoogewerff J, Foden J, Vroon P, Zulkarnain I, Nasution A. Geochemical trends across an arc–continent collision zone: magma sources and slab–wedge transfer processes below the Pantar Strait volcanoes, Indonesia. Geochimica et Cosmochimica Acta (2002) 66:2771–2789.[CrossRef][Web of Science]
Ellam RM, Harmon RS. Oxygen isotope constraints on the crustal contribution to the subduction-related magmatism of the Aeolian Islands, southern Italy. Journal of Volcanology and Geothermal Research (1990) 44:105–122.[CrossRef][Web of Science]
Ellam RM, Hawkesworth CJ. Elemental and isotopic variations in subduction related basalts: evidence for a three component model. Contributions to Mineralogy and Petrology (1988) 98:72–80.[CrossRef][Web of Science]
Elliott T, Plank T, Zindler A, White W, Bourdon B. Element transport from slab to volcanic front at the Mariana Arc. Journal of Geophysical Research (1997) 102:14991–15019.[CrossRef]
Fryer P. A synthesis of Leg 125 drilling of serpentine seamounts on the Mariana and Izu–Bonin forearcs. In: Proceedings of the Ocean Drilling Program, Scientific Results, 125—Fryer P, Pearce JA, Stokking LB, et al, eds. (1992) 593–614. College Station, TX: Ocean Drilling Program.
Gamble J, Woodhead J, Wright I, Smith I. Basalt and sediment geochemistry and magma petrogenesis in a transect from oceanic island arc to rifted continental margin arc: the Kermadec–Hikurangi Margin, SW Pacific. Journal of Petrology (1996) 37:1523–1546.
Gasparon M, Varne R. Crustal assimilation versus subducted sediment input in west Sunda arc volcanics: an evaluation. Mineralogy and Petrology (1998) 64:89–117.[CrossRef][Web of Science]
Gasparon M, Hilton DR, Varne R. Crustal contamination processes traced by helium isotopes: examples from the Sunda arc, Indonesia. Earth and Planetary Science Letters (1994) 126:15–22.[CrossRef][Web of Science]
Gerbe M-C, Gouraud A, Sigmarsson O, Harmon RS, Joron J-L, Provost A. Mineralogical and geochemical evolution of the 1982–1983 Galunggung eruption (Indonesia). Bulletin of Volcanology (1992) 54:284–298.[CrossRef][Web of Science]
Gertisser R, Keller J. Trace element and Sr, Nd, Pb and O isotope variations in medium-K and high-K volcanic rocks from Merapi Volcano, Central Java, Indonesia: evidence for the involvement of subducted sediments in Sunda Arc magma genesis. Journal of Petrology (2003) 44:457–489.
Giannetti B. Origin of the calderas and evolution of Roccamonfina volcano (Roman Region, Italy). Journal of Volcanology and Geothermal Research (2001) 106:301–319.[CrossRef][Web of Science]
Grove TL, Elkins-Tanton LT, Parman SW, Chatterjee N, Muntener O, Gaetani GA. Fractional crystallization and mantle-melting controls on calc-alkaline differentiation trends. Contributions to Mineralogy and Petrology (2003) 145:515–533.[CrossRef][Web of Science]
Hamilton WB. Tectonics of the Indonesian Region. US Geological Survey, Professional Papers (1979) 1078.
Handley HK. Geochemical and Sr–Nd–Hf–O isotopic constraints on volcanic petrogenesis at the Sunda arc, Indonesia. In: PhD thesis (2006) Durham University.
Harmon RS, Gerbe MC. The 1982–83 eruption at Galunggung volcano, Java (Indonesia): oxygen isotope geochemistry of a zoned magma chamber. Journal of Petrology (1992) 33:585–609.
Hart SR, Erlank AJ, Kable EJD. Sea floor basalt alteration: some chemical and Sr isotopic effects. Contributions to Mineralogy and Petrology (1974) 44:219–230.[CrossRef][Web of Science]
Hawkesworth CJ, O’Nions RK, Arculus RJ. Nd and Sr isotope geochemistry of island arc volcanics, Grenada, Lesser Antilles. Earth and Planetary Science Letters (1979) 45:237–248.[CrossRef][Web of Science]
Hawkesworth CJ, Hergt JM, Ellam RM, McDermott F. Element fluxes associated with subduction related magmatism. Philosophical Transactions of the Royal Society of London, Series A (1991) 335:393–405.[CrossRef]
Hawkesworth CJ, Turner SP, McDermott F, Peate DW, van Calsteren PJ. U–Th isotopes in arc magmas: implications for element transfer from the subducted slab. Science (1997) 276:551–555.[CrossRef][Web of Science][Medline]
Hickey-Vargas R. Origin of the Indian Ocean-type isotopic signature in basalts form the Philippine Sea plate spreading centres: an assessment of local versus large-scale processes. Journal of Geophysical Research (1998) 103:20963–20979.[CrossRef]
Hilton DR, Craig H. A helium isotope transect along the Indonesian archipelago. Nature (1989) 342:906–908.[CrossRef]
Hoffmann-Rothe A, Ritter O, Haak V. Magnetotelluric and geomagnetic modelling reveals zones of very high electrical conductivity in the upper crust of Central Java. Physics of the Earth and Planetary Interiors (2001) 124:131–151.[CrossRef][Web of Science]
Hoogewerff JA, van Bergen MJ, Vroon PZ, Hertogen J, Wordel R, Sneyers A, Nasution A, Varekamp JC, Moens HLE, Mouchel D. U-series, Sr–Nd–Pb isotope and trace element systematics across an active island arc–continent collision zone: implications for element transfer at the slab–wedge interface. Geochimica et Cosmochimica Acta (1997) 61:1057–1072.[CrossRef][Web of Science]
Ionov DA, Harmon RS, France-Lanord C, Greenwood B, Ashchepkov IV. Oxygen isotope composition of garnet and spinel peridotites in the continental mantle: evidence from the Vitim xenolith suite, southern Siberia. Geochimica et Cosmochimica Acta (1994) 58:1463–1470.[CrossRef][Web of Science]
Ito E, White WM, Gopel C. The O, Sr, Nd and Pb isotope geochemistry of MORB. Chemical Geology (1987) 62:157–176.[CrossRef][Web of Science]
Johnson MC, Plank T. Dehydration and melting experiments constrain the fate of subducted sediments. Geochemistry, Geophysics, Geosystems (1999) 1. paper number 1999GC000014.
Kay RW. Volcanic arc magmas: implications of a melting–mixing model for element recycling in the crust–upper mantle system. Journal of Geology (1980) 88:497–522.[Web of Science]
Kemmerling GLL. Het Idjen Hoogland de geologie en geomorphologie van den Idjen. Kolff G, et al, eds. (1921) Batavia,Weltevreden: de Koninklijke Natuurkundige Vereniging. 162.
Keppler H. Constraints from partitioning experiments on the composition of subduction zone fluids. Nature (1996) 380:237–240.[CrossRef]
Kessel R, Schmidt MW, Ulmer P, Pettke T. Trace element signature of subduction-zone fluids, melts and supercritical liquids at 120–180 km depth. Nature (2005) 437:724–727.[CrossRef][Medline]
Macpherson CG, Hall R. Tectonic setting of Eocene boninite magmatism in the Izu–Bonin–Mariana forearc. Earth and Planetary Science Letters (2001) 186:215–230.[CrossRef][Web of Science]
Macpherson CG, Hall R. Timing and tectonic controls in the evolving orogen of SE Asia and the western Pacific and some implications for ore generation. The Timing and Location of Major Ore Deposits in an Evolving Orogen—Blundell DJ, Neubauer F, von Quadt A, eds. (2002) 204. Geological Society, London: Special Publication. 49–67.
Macpherson CG, Mattey DP. Oxygen isotope variations in Lau Basin lavas. Chemical Geology (1998) 144:177–194.[CrossRef][Web of Science]
Macpherson CG, Gamble JA, Mattey DP. Oxygen isotope geochemistry of lavas from an oceanic to continental arc transition, Kermadec–Hikurangi margin, SW Pacific. Earth and Planetary Science Letters (1998) 160:609–621.[CrossRef][Web of Science]
Macpherson CG, Hilton DR, Mattey DP, Sinton JM. Evidence for an 18O-depleted mantle plume from contrasting 18O/16O ratios of back-arc lavas from the Manus Basin and Mariana Trough. Earth and Planetary Science Letters (2000) 176:171–183.[CrossRef][Web of Science]
Macpherson CG, Forde EJ, Hall R, Thirlwall MF. Geochemical evolution of magmatism in an arc–arc collision; the Halmahera and Sangihe Arcs, eastern Indonesia. Intra-oceanic Subduction Systems: Tectonic and Magmatic Processes—Larter RD, Leat PT, eds. (2003) 219. Geological Society, London: Special Publication. 207–220.
Mandeville CW, Carey S, Sigurdsson H. Magma mixing, fractional crystallisation and volatile degassing during the 1883 eruption of Krakatau volcano, Indonesia. Journal of Volcanology and Geothermal Research (1996) 74:243–274.[CrossRef][Web of Science]
Mattey D, Lowry D, Macpherson C. Oxygen isotope composition of mantle peridotite. Earth and Planetary Science Letters (1994) 128:231–241.[CrossRef][Web of Science]
McCulloch MT, Gamble JA. Geochemical and geodynamical constraints on subduction zone magmatism. Earth and Planetary Science Letters (1991) 102:358–374.[CrossRef][Web of Science]
McCulloch MT, Perfit MR. 143Nd/144Nd, 87Sr/86Sr and trace element constraints on the petrogenesis of Aleutian island arc magmas. Earth and Planetary Science Letters (1981) 56:167–179.[CrossRef][Web of Science]
Münker C, Worner G, Yogodzinski G, Churikova T. Behaviour of high field strength elements in subduction zones: constraints from Kamchatka–Aleutian arc lavas. Earth and Planetary Science Letters (2004) 224:275–293.[CrossRef][Web of Science]
Nowell GM, Kempton PD, Noble SR, Fitton JG, Saunders AD, Mahoney JJ, Taylor RN. High precision Hf isotope measurements of MORB and OIB by thermal ionisation mass spectrometry: insights into the depleted mantle. Chemical Geology (1998) 149:211–233.[CrossRef][Web of Science]
Holland G, Tanner SD. Long-term performance characteristics of a plasma ionisation multi-collector mass spectrometer (PIMMS): the ThermoFinnigan Neptune. In: Royal Society of Chemistry (2003) Spec. Pub. 307–320. Plasma Source Mass Spectrometry.
Ottley CJ, Pearson DG, Irvine GJ. A routine method for the dissolution of geological samples for the analysis of REE and trace elements via ICP-MS. In: Plasma Source Mass Spectrometry. Special Publication, Royal Society of Chemistry—Holland G, Tanner SD, eds. (2003) 221–230.
Patchett PJ. Hafnium isotope results from mid-ocean ridges and Kerguelen. Lithos (1983) 16:47–51.[CrossRef][Web of Science]
Patchett PJ, Tatsumoto M. Hafnium isotope variations in oceanic basalts. Geophysical Research Letters (1980) 7:1077–1080.[Web of Science]
Peacock SM. Numerical simulation of subduction zone pressure–temperature–time paths: constraints on fluid production and arc magmatism. Philosophical Transactions of the Royal Society of London, Series A (1991) 335:341–353.
Pearce JA, Peate DW. Tectonic implications of the composition of volcanic arc magmas. Annual Review of Earth and Planetary Sciences (1995) 24:251–285.[CrossRef]
Pearce JA, Kempton PD, Nowell GM, Noble SR. Hf–Nd element and isotope perspective on the nature and provenance of mantle and subduction components in western Pacific arc–basin systems. Journal of Petrology (1999) 40:1579–1611.[CrossRef][Web of Science]
Plank T, Langmuir CH. Tracing trace elements from sediment input to volcanic output at subduction zones. Nature (1993) 362:739–743.[CrossRef]
Plank T, Langmuir CH. The chemical composition of subducting sediment and its consequences for the crust and mantle. Chemical Geology (1998) 145:325–394.[CrossRef][Web of Science]
Potts PJ, Tindle AG, Webb PC. Geochemical Reference Material Compositions: Rocks, Minerals, Sediments, Soils, Carbonates, Refractories and Ores Used in Research and Industry (1992) Caithness: Whittles Publishing.
Price RC, Kennedy AK, Riggs-Sneeringer M, Frey FA. Geochemistry of basalts from the Indian Ocean triple junction: implications for the generation and evolution of Indian Ocean ridge basalts. Earth and Planetary Science Letters (1986) 78:379–396.[CrossRef][Web of Science]
Rehkämper M, Hofmann AW. Recycled ocean crust and sediment in Indian Ocean MORB. Earth. and Planetary Science Letters (1997) 147:93–106.[CrossRef]
Reubi O, Nicholls IA. Magmatic evolution at Batur volcanic field, Bali, Indonesia: petrological evidence for polybaric fractional crystallisation and implications for caldera-forming eruptions. Journal of Volcanology and Geothermal Research (2004) 138:345–369.[CrossRef][Web of Science]
Ringwood AE. The petrological evolution of island arc systems. Journal of the Geological Society, London (1974) 130:183–204.
Rotolo SG, Castorina F. Transition from mildly-tholeiitic to calc-alkaline suite: the case of Chichontepec volcanic centre. El Salvador, Central America. Journal of Volcanology and Geothermal Research (1998) 86:117–136.[CrossRef][Web of Science]
Royse K, Kempton PD, Darbyshire DPF. NERC Isotope Geosciences Laboratory Report Series (1998) 121. Procedure for the analysis for rubidium-strontium and samarium-neodymium isotopes at the NERC Isotope Geosciences Laboratory.
Salters VJM. The generation of mid-ocean ridge basalts from the Hf and Nd isotope perspective. Earth and Planetary Science Letters (1996) 141:109–123.[CrossRef][Web of Science]
Salters VJM, Hart SR. The mantle sources of ocean islands and arc basalts: the Hf isotope connection. Earth and Planetary Science Letters (1991) 104:364–380.[CrossRef][Web of Science]
Salters VJM, White WM. Hf isotope constraints on mantle evolution. Chemical Geology (1998) 145:447–460.[CrossRef][Web of Science]
Sitorus K. Volcanic stratigraphy and geochemistry of the Idjen Caldera Complex, East Java, Indonesia. In: MSc thesis (1990) New Zealand: University of Wellington.
Smith TE, Thirlwall MF, Macpherson C. Trace element and isotope geochemistry of the volcanic rocks of Bequia, Grenadine Islands, Lesser Antilles Arc; a study of subduction enrichment and intra-crustal contamination. Journal of Petrology (1996) 37:117–143.
Smyth H. Eocene to Miocene basin history and volcanic activity in East Java, Indonesia. In: PhD thesis (2005) Royal Holloway University of London.
Staudigel H, Davies GR, Hart SR, Marchant KM, Smith BM. Large scale isotopic Sr, Nd and O isotopic anatomy of altered oceanic crust: DSDP/ODP sites 417/418. Earth and Planetary Science Letters (1995) 130:169–185.[CrossRef][Web of Science]
Staudigel H, Plank T, White B, Schimincke HU. Geochemical fluxes during seafloor alteration of the basaltic upper oceanic crust: DSDP Sites 417 and 418. In: Subduction: Top to Bottom. Geophysical Monograph, American Geophysical Union 96—Bebout GE, Scholl DW, Kirby SH, Platt JP, eds. (1996) 19–36.
Stille P, Unruh DM, Tatsumoto M. Pb, Sr, Nd, and Hf isotopic constraints on the origin of Hawaiian basalts and evidence for a unique mantle source. Geochimica et Cosmochimica Acta (1986) 50:2303–2319.[CrossRef][Web of Science]
Stolz AJ, Varne R, Davies GR, Wheller GE, Foden JD. Magma source components in an arc–continent collision zone: the Flores–Lembata sector, Sunda arc, Indonesia. Contributions to Mineralogy and Petrology (1990) 105:585–601.[CrossRef][Web of Science]
Stormer JC Jr, Nicholls J. XLFRAC; a program for the interactive testing of magmatic differentiation models. Computers & Geosciences (1978) 4:143–159.[CrossRef][Web of Science]
Sun S, McDonough WF. Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. Magmatism in the Ocean Basins—Saunders AD, Norry MJ, eds. (1989) 42. Geological Society, London: Special Publications. 313–345.
Tatsumi Y, Hamilton DL, Nesbitt RW. Chemical characteristics of fluid phase released from a subducted lithosphere and origin of arc magmas: evidence from high-pressure experiments and natural rocks. Journal of Volcanology and Geothermal Research (1986) 29:293–309.[CrossRef][Web of Science]
Taylor SR, McLennan SM. The Continental Crust; its Composition and Evolution (1985) Oxford: Blackwell Scientific.
Taylor RN, Nesbitt RW. Isotopic characteristics of subduction fluids in an intra-oceanic setting, Izu–Bonin Arc, Japan. Earth and Planetary Science Letters (1998) 164:79–98.[CrossRef][Web of Science]
Taylor RN, Nesbitt RW, Vidal P, Harmon RS, Auvray B, Croudace IW. Mineralogy, chemistry, and genesis of the boninite series volcanics, Chichijima, Bonin Islands, Japan. Journal of Petrology (1994) 35:577–617.
Thirlwall MF. Long-term reproducibility of multicollector Sr and Nd isotope ratio analysis. Chemical Geology; Isotope Geoscience Section (1991) 94:85–104.[CrossRef]
Thirlwall MF, Graham AM. Evolution of high-Ca, high-Sr C-series basalts from Grenada, Lesser Antilles: the effects of intra-crustal contamination. Journal of the Geological Society, London (1984) 141:427–445.
Thirlwall MF, Graham AM, Arculus RJ, Harmon RS, Macpherson CG. Resolution of the effects of crustal assimilation, sediment subduction, and fluid transport in island arc magmas: Pb–Sr–Nd–O isotope geochemistry of Grenada, Lesser Antilles. Geochimica et Cosmochimica Acta (1996) 60:4785–4810.[CrossRef][Web of Science]
Tregoning P, Brunner FK, Bock Y, Puntodewo SSO, McCaffrey R, Genrich JF, Calais E, Rais J, Subarya C. First geodetic measurement of convergence across the Java Trench. Geophysical Research Letters (1994) 21:2135–2138.[CrossRef][Web of Science]
Turner S, Foden J. U, Th and Ra disequilibria, Sr, Nd and Pb isotope and trace element variations in Sunda arc lavas: predominance of a subducted sediment component. Contributions to Mineralogy and Petrology (2001) 142:43–57.[Web of Science]
Turner S, Hawkesworth C. Constraints on flux rates and mantle dynamics beneath island arcs from Tonga–Kermadec lava geochemistry. Nature (1997) 389:568–573.[CrossRef]
Turner S, Hawkesworth C, Rogers N, Bartlett J, Worthington T, Hergt J, Pearce J, Smith I. 238U–230Th disequilibria, magma petrogenesis and flux rates beneath the depleted Tonga–Kermadec island arc. Geochimica et Cosmochimica Acta (1997) 61:4855–4884.[CrossRef][Web of Science]
Turner S, Foden J, George R, Evans P, Varne R, Elburg M, Jenner G. Rates and processes of potassic magma evolution beneath Sangeang Api volcano, East Sunda Arc, Indonesia. Journal of Petrology (2003) 44:491–515.
van Bemmelen RW. The Geology of Indonesia. Vol 1A (1949) The Hague: Government Printing Office.
van Bergen MJ, Vroon PZ, Varekamp JC, Poorter RPE. The origin of the potassic rock suite from Batu Tara Volcano (East Sunda Arc, Indonesia). Lithos (1992) 28:261–282.[CrossRef][Web of Science]
Vervoort JD, Patchett PJ, Blichert-Toft J, Albarède F. Relationships between Lu–Hf and Sm–Nd isotopic systems in the global sedimentary system. Earth and Planetary Science Letters (1999) 168:79–99.[CrossRef][Web of Science]
Vroon PZ. Subduction of continental material in the Banda Arc, Eastern Indonesia: Sr–Nd-Pb isotope and trace-element evidence from volcanics and sediments. In: PhD thesis (1992) University of Utrecht.
Vroon PZ, van Bergen MJ, Klaver GJ, White WM. Strontium, neodymium, and lead isotopic and trace-element signatures of the East Indonesian sediments: provenance and implications for Banda Arc magma genesis. Geochimica et Cosmochimica Acta (1995) 59:2573–2598.[CrossRef][Web of Science]
Vroon PZ, Lowry D, van Bergen MJ, Boyce AJ, Mattey DP. Oxygen isotope systematics of the Banda arc: low 18O despite involvement of subducted continental material in magma genesis. Geochimica et Cosmochimica Acta (2001) 65:589–609.[CrossRef][Web of Science]
Wheller GE, Varne R, Foden JD, Abbott MJ. Geochemistry of Quaternary volcanism in the Sunda–Banda Arc, Indonesia, and three-component genesis of island arc basaltic magmas. Journal of Volcanology and Geothermal Research (1987) 32:137–160.[CrossRef][Web of Science]
White WM, Dupré B. Sediment subduction and magma genesis in the Lesser Antilles: isotopic and trace element constraints. Journal of Geophysical Research (1986) 91:5927–5941.
White WM, Patchett J. Hf–Nd–Sr isotopes and incompatible element abundances in island arcs: implications for magma origins and crustal–mantle evolution. Earth and Planetary Science Letters (1984) 67:167–185.[CrossRef][Web of Science]
White WM, Patchett J, Ben Othman D. Hf isotope ratios of marine sediments and Mn nodules; evidence for a mantle source of Hf in seawater. Earth and Planetary Science Letters (1986) 79:46–54.[CrossRef][Web of Science]
Whitford DJ. Strontium isotopic studies of the volcanic rocks of the Sunda arc, Indonesia, and their petrogenesis. Geochimica et Cosmochimica Acta (1975) 39:1287–1302.[CrossRef][Web of Science]
Whitford DJ, White WM, Jezek PA. Neodymium isotopic composition of Quaternary island arc lavas from Indonesia. Geochimica et Cosmochimica Acta (1981) 45:989–995.[CrossRef][Web of Science]
Woodhead JD, Eggins SM, Gamble J. High field strength and transition element systematics in island arc and back-arc basin basalts: evidence for multi-phase melt extraction and depleted mantle wedge. Earth and Planetary Science Letters (1993) 114:491–504.[CrossRef][Web of Science]
Woodhead JD, Hergt JM, Davidson JP, Eggins SM. Hafnium isotope evidence for ‘conservative’ element mobitiy during subduction processes. Earth and Planetary Science Letters (2001) 192:331–346.[CrossRef][Web of Science]
You C-F, Castillo PR, Gieskes JM, Chan LH, Spivack AJ. Trace element behaviour in hydrothermal experiments: implications for fluid processes at shallow depths in subduction zones. Earth and Planetary Science Letters (1996) 140:41–52.[CrossRef][Web of Science]
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  T. Kuritani, T. Yokoyama, and E. Nakamura Generation of Rear-arc Magmas Induced by Influx of Slab-derived Supercritical Liquids: Implications from Alkali Basalt Lavas from Rishiri Volcano, Kurile Arc J. Petrology, July 1, 2008; 49(7): 1319 - 1342. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||