| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Journal of Petrology Volume 42 Number 12 Pages 2303-2331 2001
© Oxford University Press 2001
Geochemical Evolution of Akagi Volcano, NE Japan: Implications for Interaction Between Island-arc Magma and Lower Crust, and Generation of Isotopically Various Magmas
PHEASANT MEMORIAL LABORATORY OF GEOCHEMISTRY AND COSMOCHEMISTRY, INSTITUTE FOR STUDY OF THE EARTH’S INTERIOR, OKAYAMA UNIVERSITY AT MISASA, TOTTORI, 682-0193, JAPAN
Received April 30, 2000; Revised typescript accepted June 10, 2001
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONRESULTS AND DISCUSSIONINTERACTION BETWEEN ISLAND-ARC...EXISTENCE OF A WATER-RICH...EVOLUTION OF MAGMA IN...TEMPORAL ISOTOPIC VARIATION OF...CONCLUSIONSREFERENCES |
|---|
Major and trace element, and Sr, Nd and Pb isotopic compositions were determined for whole-rock samples from the ‘isotopically anomalous’ Akagi volcano in the volcanic front of the NE Japan arc. Sr and Nd isotopic compositions of phenocrysts were also analyzed together with their major and trace element compositions. Compared with the other volcanoes from the volcanic front, the whole-rock isotope compositions of Akagi show highly enriched characteristics; 87Sr/86Sr = 0·7060–0·7088,
Nd = -0·40 to -8·6, and 208Pb/204Pb = 38·4–38·8. The rare earth element (REE) patterns are characterized by heavy REE (HREE) depletions with U-shaped patterns from middle REE (MREE) to HREE, suggesting that amphibole fractionation was induced by a reaction between clinopyroxene and H2O-rich magma in the lower crust. The integrated isotope and trace elements systematics, and tectonic structure beneath Akagi volcano, suggest that lower-crustal assimilation by the H2O-rich primary magma could have been affected by the double subduction of Philippine Sea and Pacific oceanic plates. This double subduction could have supplied larger amounts of water to the magma source region in the wedge mantle than in the case of a single subduction zone. Significant differences in isotopic compositions are observed between phenocrysts and the coexisting melts. Such isotopic disequilibrium may have resulted from magma mixing between an isotopically depleted aphyric and an enriched porphyritic magma in a shallow magma chamber. The geochemical characteristics of these end-member magmas were retained in the lower crust, despite differing extents of lower-crustal assimilation by the H2O-rich magmas. KEY WORDS: Akagi volcano; H2O-rich magma; isotopic disequilibrium; lower-crustal assimilation; magma mixing
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTS AND DISCUSSIONINTERACTION BETWEEN ISLAND-ARC...EXISTENCE OF A WATER-RICH...EVOLUTION OF MAGMA IN...TEMPORAL ISOTOPIC VARIATION OF...CONCLUSIONSREFERENCES |
|---|
Island-arc volcanism has been considered to result from partial melting of the wedge mantle induced by addition of slab-derived fluid-rich materials in subduction zones. This general model is supported by many geochemical studies employing major and trace element compositions, and radiogenic isotope systematics of arc volcanic ejecta (e.g. Nakamura et al., 1985
; Woodhead & Fraser, 1985; Tatsumi et al., 1986; Ishikawa & Nakamura, 1994; Miller et al., 1994; Ryan et al., 1995; Shibata & Nakamura, 1997). These previous works mainly focused on the source characteristics in relation to the evolution of the mantle wedge in subduction zones. The magmatic evolution of individual volcanoes, after partial melting of the wedge mantle, and in shallow magma chambers has received less attention. Isotope and trace element geochemistry has typically been restricted to basaltic volcanic rocks, not the felsic rocks, to avoid the complexities of processes such as crystal fractionation, magma mixing and crustal assimilation. It is, however, essential to characterize the shallower magmatic processes along with the source characteristics. Magmas erupted at the Earth’s surface possess original source-related characteristics considerably modified by post-melting physico-chemical processes that may obscure the source signature.
In this study, Akagi volcano in the NE Japan arc was investigated because of its distinctive isotopic composition relative to neighboring volcanoes from the volcanic front of the NE Japan arc. This arc is one of the ideal localities for discussing island-arc magmatism, because the structure of the subduction zone is well defined by abundant seismological data (e.g. Utsu, 1974; Yoshii, 1979; Zhao et al., 1990; Zhao, 1992). Moreover, systematic across-arc variations of chemical and isotopic compositions were correlated with the depth of the Wadati–Benioff Zone (Notsu, 1983; Nakamura et al., 1985; Sakuyama & Nesbitt, 1986; Shibata & Nakamura, 1997). These across-arc systematics have been explained by continuous dehydration of the subducting oceanic slab resulting in a continuous decrease in the slab component in the source region with increasing depth to the subducted slab (Shibata & Nakamura, 1997). Volcanoes in the northern part of the volcanic front show a relatively narrow range of 87Sr/86Sr ratio (0·7038–0·7045) (Notsu, 1983; Shibata & Nakamura, 1997). However, the Sr isotope composition in the southern part of the volcanic front increases to the south and has the highest 87Sr/86Sr ratio of 
0·7087 at Akagi volcano (Notsu, 1983; Kersting et al., 1996). According to Gust et al. (1997), such systematic isotopic change along the volcanic front in the NE Japan arc could be explained by the contribution of subcontinental lithospheric mantle with some crustal contamination based on Sr–Nd–Pb isotope systematics and trace element geochemistry. They also suggested that the primary magmas of the volcanoes located south of the Tanakura Tectonic Line might have been involved with Indian-Ocean-type mantle, although they pointed out the possibility of lower-crustal contamination in the petrogenesis of the felsic volcanic rocks. However, their data are limited for individual volcanoes (only one to four samples without petrological description), making it difficult to discriminate samples that might have been involved with crustal contamination. It is, therefore, risky to discuss their source materials in the mantle beneath individual volcanic centers based on previous studies with insufficient petrological description of samples or a systematic dataset for each volcano.
From such a point of view, we investigated Akagi volcano, measuring major and trace element compositions and multi-isotope compositions including Sr, Nd and Pb isotopes for whole-rock samples, to comprehensively understand the magmatic processes involved in the formation of the isotopically most anomalous volcano in the volcanic front of the NE Japan arc, and to assess the mechanisms of along-arc variation with respect to tectonic setting. In addition, the trace element, and Sr and Nd isotope compositions of phenocryst minerals were determined to further understand the shallow magma chamber processes and thereby elucidate the genesis of andesite magma in subduction zones.
Geological and petrological features of Akagi volcano
Akagi volcano (36°33'N, 139°12'E) is located in the southern part of the volcanic front of the NE Japan arc (Fig. 1). This area is characterized by the presence of many Quaternary volcanoes (e.g. Haruna, Hotaka, Nikko-Shirane volcano), and is one of the most volcanically active regions in Japan. As illustrated in Fig. 1, an unusual structure exists in the mantle wedge beneath the volcanic area including Akagi volcano; that is, a double subduction zone, in which the Philippine Sea plate is subducted into the wedge mantle above the Pacific plate (Ishida, 1991). The tip of the Philippine Sea plate is subducted to depth of 90 km just under Akagi volcano, and it probably adheres closely to the subducting Pacific plate, of which the interface depth is 110 km. Moreover, the depth of the Moho under Akagi volcano is 39 km, some 4 km deeper than is typical beneath the volcanic front of the NE Japan arc (Zhao et al., 1990; Zhao, 1992). These unusual tectonic conditions probably contribute to the formation of the distinctive geochemical characteristics of Akagi volcano.
|
Although Akagi volcano is known to be mostly Quaternary in age, the lack of systematic radiometric ages precludes a clear understanding of the commencement of volcanic activity. On the basis of stratigraphical relationships between Akagi volcano (Moriya, 1968) and the neighboring Komochi volcano (Iizuka, 1996), the activity of Akagi volcano is thought to have started before 350 ka, and lasted intermittently to 30 ka. The latter age is based on the only available radiometric age, determined by the 14C method by Koga (1981).
The volcanic history of Akagi volcano is divided into three main stages based on the eruption style: the older stratovolcano formation stage (O) characterized by the eruption of andesitic lava flow and ejecta of scoria, the younger stratovolcano formation stage (Y) characterized by pyroclastics without lava flows, and the central cone formation stage (Cc) (Moriya, 1968). The older stratovolcano formation stage (O) is further divided into three substages based on the volcanic stratigraphy: the early (Oe), middle (Om) and late (Ol) substages. The main volcanic stages and substages are geologically distinguishable. However, the eruptive sequence of volcanic ejecta within each stage is unclear based on stratigraphy. The present volcanic flank is mainly formed of pyroclastics of the Y stage and of secondary deposits. Lava flows can be observed in the caldera, which formed after the activity of the Y stage, and near the caldera wall.
Sample locality and description
To investigate the magmatic processes of Akagi volcano, volcanic ejecta were collected based on the volcano stratigraphy established by Moriya (1968). The sampling localities are presented in Fig. 1. Samples were picked from fresh parts of outcrops to minimize alteration effects, and examined under the petrographic microscope to confirm the absence of significant secondary alteration.
Petrographic observations indicate that the Akagi rock samples are highly porphyritic, with phenocrysts mainly consisting of plagioclase, hypersthene, augite and opaque minerals (see Table 1). Some of the felsic samples with SiO2 contents exceeding 58 wt % in the Ol substage and Y and Cc stages contain phenocrysts of hornblende instead of augite, even without quartz phenocrysts. Olivine phenocrysts are absent in the studied samples except for a sample AK1303, which contains an olivine pseudomorph surrounded by orthopyroxene. The mineral assemblage in the groundmass of the Oe and Om substages is plagioclase + pigeonite ± augite ± hypersthene. In contrast, samples from the Ol substage and the Y and Cc stage samples do not contain groundmass pigeonite.
|
Sample preparation
In the preparation of whole-rock samples for major, trace element and isotope analyses, chunks free from surface alteration were picked. These unaltered chunks were washed with Milli-Q water using an ultrasonic bath, and dried. Then, they were ground to fine powders with grain sizes under 200 mesh using a silicon nitride mortar.
To separate phenocrysts, rock powders obtained using a disk mill were sieved to collect materials with grain size ranging from 100 to 200 mesh, as the average grain size of phenocrysts of plagioclase, orthopyroxene, and clinopyroxene or hornblende is 500 µm. The sieved fraction was then processed to purify the minerals using conventional heavy liquid methods followed by an isodynamic separator. The mineral separation was finally accomplished by hand-picking, and the purity is regarded as being better than 99%.
Analytical methods
All geochemical analyses were carried out at the Pheasant Memorial Laboratory (PML), Institute for Study of the Earth’s Interior, Okayama University at Misasa. Major element compositions of whole rocks were determined by X-ray fluorescence (XRF) using a Philips PW-2400 system; details of the analytical procedure are described elsewhere (Takei et al., in preparation). Trace element analyses of the whole rocks were performed by inductively coupled plasma mass spectrometry (ICP-MS) using a Yokogawa PMS 2000 instrument with a flow-injection method developed by Makishima & Nakamura (1997). The analytical reproducibility for trace element analyses of andesitic samples was <5%, and typically 3% (Makishima & Nakamura, 1997).
The analytical procedures, including the chemical separations and mass spectrometry utilized in this study, are from Yoshikawa & Nakamura (1993), Makishima & Nakamura (1991a, 1991b) and Koide & Nakamura (1990) for Sr, Nd and Pb isotope analysis, respectively. Mass spectrometry was carried out on Finnigan MAT 261 and MAT 262 instruments equipped with five Faraday cups, applying a static-multi-collection mode. Normalizing factors used to correct isotopic fractionation of Sr and Nd are 86Sr/88Sr = 0·1194 and 146Nd/144Nd = 0·7219, respectively. The Pb isotope analyses were carried out with a nearly fixed filament temperature of 1200°C. Measured ratios of reference materials were 87Sr/86Sr = 0·710233 ± 8 (2
m) for 100 ng of NIST SRM 987, 143Nd/144Nd = 0·511818 ± 13 (2m) for 30 ng of La Jolla, and 206Pb/204Pb = 16·940 ± 15 (2m), 207Pb/204Pb = 15·494 ± 16 (2m) and 208Pb/204Pb = 36·709 ± 38 (2m) for 100 ng of NIST SRM 981.
Major element compositions of phenocrysts were determined by electron probe microanalysis (EMPA), using a JEOL JXA-8800R instrument, at Misasa following the techniques of Iizuka (1996). Trace element analyses of clinopyroxene, orthopyroxene and plagioclase phenocrysts were carried out employing a Cameca ims 5f ion microprobe at PML, Misasa, following procedures described by Nakamura & Kushiro (1997, 1998). Standard materials used for the calibration of trace element compositions were clinopyroxene from mantle xenoliths for pyroxene and hornblende analyses, and labradorite megacrysts for plagioclase analyses. These standards have been well characterized in terms of homogeneity and concentrations by both ion microprobe and ICP-MS. Clinopyroxene phenocrysts in the thin sections were sputtered with an O- primary beam of 10 nA intensity, resulting in 10 µm beam diameter, and orthopyroxene and plagioclase with 15–20 nA intensity, resulting in 15 µm beam diameter. Positive secondary ions were collected by ion counting using an energy offset of -60 V for Si, Rb, Sr, Y and Zr, and of -45 V for other elements from 4500 V acceleration with an energy bandpass of ±10 V. These operational conditions resulted in (1-2) x 105 c.p.s. for 30Si secondary ion in the analysis of phenocrysts. Analytical reproducibilities (RSD 1%, n = 10) for trace elements in the clinopyroxene standard were typically 5%, except for Ba (60%), and for the labradorite standard also typically 5%, except for Zr, Nb, Sm and heavy rare earth elements (HREE) (20–30%).
|
RESULTS AND DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTS AND DISCUSSIONINTERACTION BETWEEN ISLAND-ARC...EXISTENCE OF A WATER-RICH...EVOLUTION OF MAGMA IN...TEMPORAL ISOTOPIC VARIATION OF...CONCLUSIONSREFERENCES |
|---|
Major element compositions of whole-rock samples
Major and trace element compositions of the Akagi samples are given in Table 2. Selected major element oxides are plotted against SiO2 content in Fig. 2, classifying the samples into the five stages defined by Moriya (1968)
with reference to the Nasu volcanic zone, which forms the volcanic front of NE Japan (Kawano et al., 1961). As shown in Fig. 2, the range of SiO2 contents of the Akagi samples is between 53·4 and 71·4 wt %, consistent with those of previous studies, indicating that basaltic rocks have not been discovered at Akagi (Koga, 1984; Yamaguchi, 1990). TiO2, Fe2O3, MgO and CaO show negative trends against SiO2, and these trends are essentially consistent with those of the Nasu volcanic zone. Two different trends, the relatively steeper trend formed by the samples in the Oe and Om stages and that in the Ol, Y and Cc stages, can be recognized in the Fe2O3 and MgO diagrams, as well as in the Na2O diagram.
|
|
The Al2O3–SiO2 diagram shows a more complicated feature, with a positive trend for the earlier stages of Oe and Om, and a negative trend for the later stages of Ol, Y and Cc. The former positive trend may be explained by considerable plagioclase accumulation in a shallow magma chamber, as the deduced Al2O3 contents in the melts obtained by subtracting the phenocryst compositions from those of the whole rocks form a negative trend against SiO2 overlapping the samples from the later stages as shown in Fig. 3. The major element compositions, after subtraction of the phenocrysts in Adatara volcano (the neighboring volcano located in the volcanic front of NE Japan), are also shown based on the data given by Fujinawa (1988, 1990). The corrected Al2O3 concentrations for Akagi volcano indicate a similar negative trend to that of Adatara volcano; however, that of Akagi is still systematically higher than for Adatara at similar SiO2 contents. The relatively high content of Al2O3 in the Akagi samples is consistent with the results of Yamaguchi (1990). The differentiation trend of corrected CaO against SiO2 (Fig. 3), which is also influenced by plagioclase accumulation, shows similar negative trends with increasing SiO2 content. However, the corrected Na2O contents are scattered and the differentiation trends are less obvious.
|
Two distinct trends in the co-variation of FeO* (total Fe as FeO) and MgO of Akagi volcano can be divided into different rock series using the Miyashiro diagram (Fig. 4a; Miyashiro, 1974). In this diagram, the Oe and Om suites belong to the tholeiitic rock series, and the Ol, Y and Cc suites to the calc-alkaline rock series. On the other hand, although the Oe and Om suites are slightly less depleted in FeO* than the Ol, Y and Cc suites, the entire Akagi volcanic suite follows a calc-alkaline differentiation trend in the AFM diagram (Fig. 4b). These major element characteristics of Akagi volcano cannot be explained by a simple magmatic differentiation process, because there appear to be considerable gaps in the SiO2 and total alkali contents in the transition from the Om to the Ol suites (Figs
2 and 4a). It may be, therefore, necessary to invoke more complex processes, such as magma mixing and crustal contamination associated with crystal fractionation. Moreover, the difference in the eruption style, which is the basis for the discrimination of the stages, does not correspond to the characteristics of the major element compositions. It may thus be more suitable to discuss the geochemical evolution of Akagi volcano without distinctions between stages and substages. Here, we discriminate simply between samples from the five stages (Oe, Om, Ol, Y, Cc), which are distinguishable by geological discontinuities.
|
Trace element compositions
Trace element compositions of whole rocks are listed in Table 1, and are plotted in primitive mantle (PM)-normalized diagrams in Fig. 5. The contents of incompatible elements are enriched compared with most basaltic rocks of the volcanic front of the NE Japan arc (Shibata & Nakamura, 1997). However, the less incompatible elements such as the HREE are less enriched. These trace element characteristics cannot be explained by simple fractional crystallization from a common primary basaltic magma in the volcanic front using the phenocryst assemblage in the Akagi volcanic rocks. Such a fractionation process should increase not only the most incompatible trace elements but also the HREE. The behavior of REE in the Akagi samples is discussed below in more detail in a section describing their isotopic compositions. With progress within a volcanic stage, the contents of highly incompatible elements become gradually higher. However, no correlation is observed between trace elements and SiO2 contents in the Om and Ol stages. Akagi volcanic rocks show positive spikes of Sr and Pb, and remarkably negative spikes for Nb and Ta, consistent with the general features of trace elements in island-arc volcanic rocks (e.g. Wood et al., 1979; Perfit et al., 1979; Sakuyama & Nesbitt, 1986; Shibata & Nakamura, 1997; Woodhead et al., 1998).
|
Sr–Nd–Pb isotope systematics of Akagi volcano
The whole-rock Sr, Nd and Pb isotopic compositions are listed in Table 3 and are plotted in Fig. 6 along with the reference variations of volcanic rocks from the NE Japan arc, lower-crustal materials, oceanic sediments and MORBs. The 87Sr/86Sr ratios of the Akagi samples vary widely from 0·70603 to 0·70879, consistent with the data reported by Notsu et al. (1985), and are significantly higher than those of other typical volcanoes in the northern part of Tanakura Tectonic Line in the volcanic front of the NE Japan arc, which have an average 87Sr/86Sr ratio of 0·7045 (Notsu, 1983; Gust et al., 1997; Shibata & Nakamura, 1997). The Nd isotopic compositions of Akagi range from -0·4 to -8·6 in Nd, and are much lower than those of other volcanoes in the northern part of the volcanic front with Nd of 3–10 (Gust et al., 1997; Shibata & Nakamura, 1997). The 206Pb/204Pb, 207Pb/204Pb and 208Pb/204Pb ratios vary in the range of 18·21–18·44, 15·56–15·64 and 38·44–38·79, respectively (Fig. 6b and c). These isotope characteristics of Akagi, therefore, are regarded as ‘unusual’ in the NE Japan arc, confirming the original observations based on Sr isotopic compositions by Notsu et al. (1985).
|
|
In the Sr and Nd isotope diagram (Fig. 6a), the Akagi samples have a linear trend with an extremely isotopically enriched signature compared with those of other ‘typical’ volcanoes in the volcanic front of the NE Japan arc, which have simple isotopic across-arc variations (Shibata & Nakamura, 1997). These isotopic compositions are widely varied not only over the entire history of the volcano but also with each volcanic stage of Akagi volcano. An extrapolation of the Akagi trend to the depleted direction nearly intersects the isotopic trend of the NE Japan arc at a value typical of volcanoes on the front. As is shown in Fig. 6b and c, the Pb isotopic compositions of the Akagi samples form clusters that are clearly distinct from the isotopic trend of the NE Japan arc (Shibata & Nakamura, 1997). The 206Pb/204Pb ratios are lower, whereas the 207Pb/204Pb and 208Pb/204Pb ratios are higher than those of the ‘typical’ volcanoes on the front. In the Nd–206Pb/204Pb diagram (Fig. 6d), the Akagi data also define a linear trend clearly different from the tendency in the isotope variations of other primary magmas in the northern part of the NE Japan arc, and which extrapolates to the values characteristic of the volcanic front.
The above isotope systematics suggest that the Akagi volcanic rocks were formed by two-component mixing of an isotopically depleted and an enriched end-member. Intersections between the ‘typical’ trend of the NE Japan arc and the extrapolations of the Akagi trends in the Sr–Nd and Nd–Pb diagrams indicate that the depleted end-member is similar to the primary magma in the volcanic front of the NE Japan arc, which was defined by Shibata & Nakamura (1997). It is, however, still difficult to identify the source of the enriched end-member using these isotope systematics.
On the basis of the existence of an enriched component deduced from the isotope systematics and the tectonic setting of the Akagi volcano, the following candidates for the enriched end-member can be postulated: (1) oceanic sediments; (2) upper crust; (3) lower crust. Isotopic variations of oceanic sediments in the NW Pacific and lower-crustal materials are shown in Fig. 6. The isotopic composition of JG-1, the granitic standard of the Geological Survey of Japan, which outcrops near Akagi volcano and is considered to be the upper-crustal basement of Akagi, is shown in the diagrams. It should be noted that the JG-1 composition lies within the field for oceanic sediments (data from Nohda & Wasserberg, 1981; Koide & Nakamura, 1990), therefore the upper-crustal field is not shown in Fig. 6. The isotopic characteristics of lower-crust materials under the NE Japan arc have not been well understood. The lower crust under the arc might be considered as a fragment of continental material and/or a piled-up sequence of young metamorphic rocks originally derived from the arc magmas (Arculus & Johnson, 1981). In the latter case, the isotopic characteristics should be similar to those of the arc volcanic rocks, as a result of less time-integrated isotopic evolution. We prefer, therefore, that the isotopic characteristics of the lower-crustal materials are similar to ‘continental’ type lower crust in Fig. 6.
In the Sr–Nd diagram (Fig. 6a), the isotopic compositions of oceanic sediments and lower crust overlap, making it difficult to unambiguously identify an enriched end-member. However, the Pb–Pb and Nd–Pb diagrams clearly discriminate between lower crust and oceanic sediments and/or upper crust. On the basis of Fig. 6b–d, it appears that oceanic sediments and/or upper crust are not possible candidates for the enriched end-member. In practice, the isotopic compositions of sediments, which occur on the Philippine Sea plate, cannot be the enriched end-member for Akagi volcano. Consequently, it is most likely that the enriched end-member involved in magma formation of Akagi volcano is lower-crustal material, although Notsu et al. (1985) proposed a sediment component derived from the subducting Philippine Sea plate as the enriched end-member in Akagi volcano based on 87Sr/86Sr and 
18O data. The complete Sr, Nd and Pb isotopic dataset from the along-arc volcanoes by Kersting et al. (1996) and Gust et al. (1997) is also shown in Fig. 6. The isotopic trends defined by the volcanoes located to the south of the Tanakura Tectonic Line (STTL: Akagi, Nikkoshirane, Nantai, Takahara, Nasu) essentially have the same direction as those of the primary magmas in the volcanic front of the NE Japan arc (Fig. 6). Therefore, it is likely that the depleted end-member for the source of the volcanoes from the STTL also is similar to the MORB-type wedge mantle metasomatized by the fluid derived from dehydration of the subducted slab as proposed by Shibata & Nakamura (1997) for the primary magma source of the NE Japan arc. Moreover, the Nd–Pb diagram (Fig. 6d) for the STTL group clearly indicates the involvement of lower-crustal materials as an enriched end-member, as well as at Akagi volcano. It is, therefore, not necessary to introduce a unique mantle source such as Indian Ocean-type mantle beneath the volcanoes from the STTL to explain the isotope systematics of the STTL group based on the dataset given in this study and in previous studies (Kersting et al., 1996; Gust et al., 1997).
Major and trace element compositions of phenocrysts
Seven samples from the Ol stage, which is characterized by large variations in whole-rock major element and isotopic composition (SiO2 = 55·8–63·5 wt %, 87Sr/86Sr = 0·70683–0·70879), were selected to investigate detailed magmatic processes based on the chemical and isotopic compositions of phenocrysts.
Phenocrysts in Akagi volcano rocks are always zoned. The mg-numbers [defined as 100 x Mg/(Mg + Fe2+)] of cores in orthopyroxene and clinopyroxene are 75–65 and 80–73, respectively. Reverse zoning is commonly observed in the pyroxene of samples AK1102, AK0910, AK1010 and AK1108. The an-number in the cores of plagioclase phenocrysts does not have a bimodal signature, which has been regarded as evidence of magma mixing according to Sakuyama (1981). Hornblende phenocrysts occur in only the most siliceous sample in this stage. The hornblende phenocrysts have a zonal structure with decreasing mg-number and edenite component from core to rim, and a pseudomorph of hornblende phenocrysts, consisting of fine-grained plagioclase + orthopyroxene + clinopyroxene + opaque minerals, also occurs in five siliceous samples listed in Table 1.
Trace element compositions of orthopyroxene, clinopyroxene and plagioclase phenocrysts in the Ol stage are given in Table 4 and their primitive mantle-normalized patterns are shown in Figs
7 and 8. The trace element patterns of pyroxenes show clear negative Sr and Nb spikes. Experimental determinations of the Kd for Sr between pyroxenes and silicate melt show no Sr depletion compared with the adjacent elements Pr and Nd (Green, 1994). The Sr depletion, therefore, indicates that the pyroxenes crystallized from Sr-depleted melts, although Sr enrichment is observed in whole-rock analyses of the Ol samples (see Fig. 5). Such Sr enrichment could be expected as a result of the accumulation of plagioclase, which is characterized by a strong Sr positive spike (Fig. 8). The influence of plagioclase accumulation is expected based on the major element compositions of the whole-rock samples (see Fig. 3). The positive Sr spikes are removed by consideration of phenocryst subtraction; for example, Sr/Sr* value [Sr* is defined as (Pr + Nd)/2 on the primitive mantle-normalized value] of the melt part in AK-A decreases from 2·0 to 0·7.
|
|
|
The Nb depletions in the phenocrysts can be attributed to the whole-rock compositions that are unique to island-arc volcanic rocks (see also Fig. 5). The Zr depletion in clinopyroxene and the slightly positive Zr spike in orthopyroxene may be simply explained by the difference in Kd [i.e. KdHFSE/KdLILE,LREE > 1 for orthopyroxene, KdHFSE/KdLILE,LREE < 1 for clinopyroxene (where HFSE indicates high field strength elements, LILE indicates large ion lithophile elements and LREE are light REE); Green, 1994]. The Y depletion in the plagioclase phenocrysts may be a characteristic feature of plagioclase, because such depletion is always observed in the plagioclases used for our secondary ion mass spectrometry standards (also precisely characterized by ICP-MS). In samples AK1010 and AK1108, the incompatible elements in the orthopyroxene rims are significantly depleted compared with those of the core. This is consistent with reverse zoning (based on mg-number) in the pyroxene of these samples.
Sr and Nd isotope compositions of phenocrysts in the Ol stage
The isotopic compositions of Sr and Nd for phenocrysts in the Ol stage are given in Table 5 and plotted in a Sr–Nd isotope diagram in Fig. 9. In this diagram, the isotopic compositions of the melts are calculated using the Sr and Nd isotopic and elemental compositions of the whole rocks, and the modal abundance, densities, isotopic and elemental compositions of the phenocrysts. The results are given in Table 6. The isotopic compositions of the phenocrysts tend to have more enriched compositions, especially for Sr, exceeding the analytical uncertainties, than the calculated compositions of the melts. In other words, the phenocrysts are isotopically in disequilibrium with the surrounding melt.
|
|
|
The uncertainty of calculated melt isotopic compositions mainly depends on the measured modal abundance of the phenocrysts, especially plagioclase, which is a main sink for Sr. The statistical error of the plagioclase modal abundance is estimated to be 3%, based on the counting statistics. This uncertainty affects the 87Sr/86Sr ratio of melts up to ±0·0002 for AK-A and AK1012, but the uncertainties in the others are less than ±0·0001. In contrast, the calculated Nd isotopic compositions of the melts are not significantly affected by the uncertainty in the modal abundance. Heterogeneity in the sample also causes the error in the estimation of the Sr and Nd isotopic composition of the melt. However, the distributions of phenocrysts are homogeneous in the samples from the Ol stage. Therefore, the isotopically depleted characteristics of the calculated melts compared with the phenocrysts are apparently significant except in sample AK1010. In AK1010, the Sr isotopic composition of the melt is identical to that of the phenocrysts within the analytical uncertainties. However, the Nd isotopic composition of the plagioclase is significantly lower than that of the calculated melt.
|
INTERACTION BETWEEN ISLAND-ARC MAGMA AND LOWER CRUST |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTS AND DISCUSSIONINTERACTION BETWEEN ISLAND-ARC...EXISTENCE OF A WATER-RICH...EVOLUTION OF MAGMA IN...TEMPORAL ISOTOPIC VARIATION OF...CONCLUSIONSREFERENCES |
|---|
In Fig. 10, chondrite-normalized REE patterns of Akagi volcano samples are compared with those for Adatara volcano (Fujinawa, 1992
), which is a neighboring volcano in the volcanic front of the NE Japan arc and has similar petrological and mineralogical features to those of Akagi. This comparison shows that Akagi samples have (1) significantly lower HREE abundance, (2) significantly higher LREE/HREE ratios, and (3) U-shaped REE patterns between the middle REE (MREE) and HREE. The REE patterns from Adatara volcano were explained by crystal fractionation (Fujinawa, 1992). However, the REE characteristics of Akagi volcano cannot be formed by such a simple fractionation involving removal of minerals such as plagioclase, orthopyroxene, clinopyroxene and opaque minerals occurring as phenocrysts in the Akagi samples. Because the mineral–melt partition coefficients of the REE are much smaller than one, fractional crystallization in a shallow magma chamber would increase the total REE with LREE enrichment, as is observed in Adatara volcano. This mechanism, therefore, cannot explain characteristics (1) and (2) of Akagi’s REE compositions.
|
On the basis of a high-pressure melting experiment involving a basaltic composition, Forden & Green (1992) suggested that the liquidus phase is clinopyroxene under high-H2O conditions in the lower-crustal environment, and olivine is the liquidus phase at lower pressure conditions. When the temperature of the melts drops, clinopyroxene becomes unstable, and reacts with the melt to form amphibole, resulting in the evolved melt being more siliceous and decreasing the H2O content of the magma itself. The lower-crustal granulites commonly contain clinopyroxene; ascending H2O-rich magma from the mantle, therefore, might react with the clinopyroxene and fractionate amphibole as a residual phase. The MREE are preferentially partitioned into amphibole, with mineral–melt partition coefficients greater than one (Adam & Green, 1994; Sisson, 1994; Witt-Eickschen & Harte, 1994). Consequently, the evolved magma acquires a U-shaped REE pattern because of the buffering effect of residual amphibole. Partial melting experiments on lower-crustal materials under high-pressure conditions indicate that the resultant liquid has a felsic composition caused by the persistence of mafic minerals in the residue (e.g. Beard & Lofgren, 1991; Beard et al., 1993, 1994). Although amphibole fractionation processes in shallow magma chambers could also explain the U-shaped REE patterns, the occurrence of hornblende phenocrysts is rare and only very small amounts occur in felsic samples in the later volcanic stages (Table 1). This process is, therefore, unlikely to explain the overall U-shaped signature retained in the Akagi volcanic rocks.
If the interaction between H2O-rich magma and lower-crustal material is a major process controlling the REE patterns of Akagi volcano, a process of ‘assimilation–fractional crystallization’ (AFC, DePaolo, 1981) under lower-crustal conditions might be invoked. To examine the assimilation process in the lower crust, the AFC model (DePaolo, 1981) was applied to the isotope and trace element systematics in the Akagi volcano. In this model calculation, a primary magma deduced from ‘normal’ primary magma on the northern part of the volcanic front in the NE Japan arc, defined by Shibata & Nakamura (1997), is assumed as an ‘assimilator’. It is, however, difficult to determine the other parameters such as isotope and element ratios of assimilated material (hereafter denoted as ‘assimilant’) under Akagi volcano, and bulk partition coefficients between the melt and the assimilant. Thus, we fixed these parameters, and only evaluated the assimilation rate (r) as given by
 |
are used for the assimilant. The isotopic compositions of the assimilant are obtained by the extrapolation of the linear correlations in the multi-isotope systematics of Akagi volcano. The bulk partition coefficients between magma and the assimilant were obtained by using the modal abundance of residual phases consisting of plagioclase, amphibole, clinopyroxene and orthopyroxene, as determined by hydrous melting experiments at lower-crustal conditions by Beard & Lofgren (1991). The REE signature of Akagi volcanic rocks indicates that amphibole fractionation was involved in the evolution of magma under lower-crustal conditions as discussed in the previous sections. All parameters used in the AFC calculations are compiled in Table 7.
|
Representative results from the AFC model calculations are presented together with data for Akagi, Adatara and other volcanoes from NE Japan in Fig. 11a–d. There appear to be correlations between isotope and element ratios in the arc data, although the correlation is poor for Pb isotopes, probably as a result of the small variation of Pb isotope compositions and the small differences between those of the lower crust and the primary magma. These correlations are obviously different from those of basaltic rocks elsewhere in the NE Japan arc and at Adatara volcano, for which the source heterogeneity and crystal fractionation processes have caused variations in isotope and elemental ratios without crustal assimilation (Shibata & Nakamura 1997). As shown in Fig. 11, most of the Akagi data plot near the mixing lines with an assimilation rate of 0·9, although those of the earliest and final stages tend to plot above the mixing lines with assimilation rates <0·9. This observation indicates that the extent of crustal assimilation relative to the crystallization in the middle stage is larger than in the earliest and latest stages of Akagi volcano. In the middle stage, therefore, the enrichment of incompatible elements together with isotope variations may largely have resulted from significant lower-crustal assimilation by primary magma similar to ‘normal’ basaltic magmas in the northern part of the volcanic front in the NE Japan arc.
|
When the assimilation rate is nearly 1·0, the latent heat of crystallization and the heat of fusion for the assimilation process are thermally balanced, approximately. Such a circumstance is likely to be achieved under lower-crustal conditions, because the geotherm in the lower crust is such that temperatures are 600–700°C, making it feasible to reach the solidus temperature of the assimilant by intrusion of magma and with less heat loss. On the other hand, the latent heat of crystallization could be mostly consumed by raising the temperature of the wallrock to its solidus, when assimilation occurs under shallower crustal conditions. Consequently, this results in an assimilation rate smaller than that occurring under lower-crustal conditions. From the above discussion based on the AFC model calculations, it may be concluded that the assimilation occurred primarily in the lower-crustal region, resulting in the unique isotope characteristics of the magmas of Akagi volcano.
The seismologically determined Moho depth gradually increases from the northern part of NE Japan along the volcanic front, and it shows the maximum value beneath Akagi volcano (38 km, 1·3 GPa), which is 4 km deeper than the typical depth beneath the volcanic front of the NE Japan arc (Zhao et al., 1992). Hildreth & Moorbath (1988) suggested that chemical interaction between magma and crust is closely related to the crustal thickness in the Andean arc. The thicker crust might effectively increase the assimilation rate, because the higher temperatures in the bottom of the crust mean that it more easily attains its solidus temperature. Consequently, volcanoes such as Akagi located on thicker crust might experience a larger effect of assimilation of lower crust on their trace element and isotope compositions.
|
EXISTENCE OF A WATER-RICH MAGMA |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTS AND DISCUSSIONINTERACTION BETWEEN ISLAND-ARC...EXISTENCE OF A WATER-RICH...EVOLUTION OF MAGMA IN...TEMPORAL ISOTOPIC VARIATION OF...CONCLUSIONSREFERENCES |
|---|
The above discussion based on isotope and trace element systematics requires a primary magma sufficiently enriched in H2O to stabilize amphibole in the lower crust. It is, therefore, essential to understand the mechanism whereby an H2O-rich source region for the primary magma in the mantle wedge beneath Akagi volcano could form.
On the basis of petrological observations, quartz phenocrysts generally occur in felsic volcanic rocks without hydrous phenocrysts in the volcanic front of the NE Japan arc (Sakuyama, 1979). However, quartz does not appear as a phenocryst in Akagi samples, even in the felsic samples. Sakuyama (1979, 1983a, 1983b) determined the H2O content of magmas based on the crystallization sequence of phenocrysts. Following his definition, the Akagi volcanic rocks belong to Type II–III, which allows magmas to contain 3–4 wt % of H2O at 0·5 GPa. This suggests that the H2O content of the Akagi magmas is considerably higher than that of the Type I volcanoes, typical of the volcanic front in NE Japan, which contain <3 wt % of H2O.
Yamaguchi (1990) inferred that the groundmass in the Akagi volcanic rocks has a significantly higher Al2O3 content than that of the adjacent Hotaka volcano. This is because the plagioclase stability field becomes narrower relative to the diopside stability field with increasing H2O content in the magma. This observation is also consistent with the hypothesis that the primary magma of Akagi is enriched in H2O relative to other volcanoes in the volcanic front in NE Japan.
The peculiar petrological and geochemical characteristics of Akagi volcano discussed above might be explained by the development of the unusual tectonic setting beneath Akagi volcanic area, where the Philippine Sea plate overlaps onto the Pacific plate. Such a tectonic setting might lead to a significantly larger amount of water in the magma source region than that expected from subduction of a single Pacific plate.
|
EVOLUTION OF MAGMA IN THE Ol STAGE |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTS AND DISCUSSIONINTERACTION BETWEEN ISLAND-ARC...EXISTENCE OF A WATER-RICH...EVOLUTION OF MAGMA IN...TEMPORAL ISOTOPIC VARIATION OF...CONCLUSIONSREFERENCES |
|---|
Recently, magma mixing at Akagi has been discussed, on the basis of the major element composition of phenocrysts and their textures, in the pumice from the Y stage, by Horio & Umino (1995)
and Umino & Horio (1998). Those workers have suggested that high-T, less fractionated magma was periodically injected into a low-T, silicic, mushy magma in a shallow magma chamber. However, previous studies have not addressed the origin of these two different magmas because of the limited geochemical information available, especially isotope compositions, for Akagi volcano.
On the basis of the major element compositions of pyroxenes, the temperatures of the magmas in the Ol stage were determined by using the QUILF program of Andersen et al. (1993). Average compositions of the rims of clinopyroxene and orthopyroxene were used for the temperature estimation. Ferric–ferrous ratios were calculated to satisfy the stoichiometry of the pyroxenes. The samples investigated from the Ol stage can be divided into high-T (1000–1200°C) and low-T (<1000°C) groups, and tend to have a negative correlation between calculated temperature and Sr isotopic composition (Table 8 and Fig. 12). This indicates that such a variation in magma temperatures is not a consequence of the simple cooling process in the shallow magma chamber.
|
|
In the Ol samples, melts are isotopically more depleted than the coexisting phenocrysts (Fig. 9) as discussed in the previous section. Such isotopic disequilibrium observed could be explained only by a mechanical mixing between the isotopically enriched magma equilibrated with the phenocrysts and the relatively depleted aphyric magma. Mass-balance calculations, based on the trace element composition and modal abundance of the phenocrysts, indicate that the Sr/Nd ratio in the glass of the phenocryst-rich rock is lower than that in the aphyric-primary magma as a result of the equilibration of phenocrysts, especially plagioclase. If the aphyric-primary magma has a more depleted isotopic signature than the phenocryst-rich magma, the mixing curve forms a convex hyperbola in the Sr and Nd isotope diagram (Fig. 13a). That is, the Sr isotopic compositions of glasses are characteristically more depleted than those of the coexisting phenocrysts. However, the difference in Nd isotopic composition between them is smaller than the difference in Sr isotopic composition. These observations, together with model mixing, may suggest that the involvement of depleted aphyric magma mixed with phenocryst-rich magma was relatively small, and that isotopic re-equilibration among the phenocrysts and the mixed melt has not been achieved.
|
|
To examine the above argument further, a model calculation was carried out using the samples with remarkable isotope disequilibrium in Sr between melt and phenocrysts (AK1012 and AK1102) representing the high-T and the low-T group, respectively. The results are presented in Fig. 13b, based on the parameters given in Table 9. In this model, the isotope and trace element compositions of the primary aphyric magma are assumed to be those of the primary basaltic magma (a depleted end-member, Pm) in the volcanic front of NE Japan (Shibata & Nakamura, 1997
). The isotopic compositions of two end-members (Fm1 and Fm2) are represented by the values of the most enriched phenocrysts in AK1012 and AK1102, assuming that the phenocrysts in the enriched end-members were isotopically equilibrated with their host melts, when the phenocrysts crystallized. The Sr and Nd abundances in Fm1 and Fm2 were obtained by the subtraction of phenocryst compositions from the whole rock of the most mafic sample AK0910 in the Ol stage. In this model, clinopyroxene and opaque phenocrysts are not considered. Because the Kd values of Sr and Nd for clinopyroxene are almost identical (Green, 1994), and because of the extremely low concentrations of these elements in the opaque minerals, very little change in Sr/Nd is produced. The amounts of phenocrysts subtracted are assumed to be 30 wt % of plagioclase for Fm1, and 40 wt % of plagioclase and 10 wt % of hornblende for Fm2. After subtraction of the phenocrysts, the resulting major element compositions are similar to the whole-rock composition of the examined samples, AK1012 and AK1102, except for Mg and Fe (e.g. calculated CSiO2Fm1 = 55·8, CSiO2Fm2 = 59·2 wt %, and the whole rocks 56·4 and 59·5 wt %, respectively). In Fig. 13b, the calculated isotopic compositions of melts in AK1012 and AK1102 are plotted near the fractions of 0·25 (Mm1) and 0·12 (Mm2) of the Pm component on the mixing curves 1 and 2, respectively. The Sr isotopic compositions show significantly more depleted signatures than those of the whole rocks with less variation in Nd isotopic composition. This is consistent with the observation that the Sr isotopic compositions of most of the glasses are characteristically more depleted than those of the whole rocks and their constituent phenocryst minerals. As shown in Fig. 9, the differences in both Sr and Nd isotopic compositions between glasses and phenocrysts are variable from sample to sample. This may be attributed to the extent of mixing of the depleted aphyric magma with the enriched porphyritic magma and/or to the residence time after mixing in the magma chamber, or to the evolution of the primary magma by reaction with lower-crustal materials.
Although reverse zoning of pyroxenes is believed to be evidence of magma mixing (i.e. Sakuyama 1979, 1981), such a signature was not observed in the sample belonging to the high-T group, AK1012. If the depleted primary magma (Pm), derived from a strongly metasomatized wedge mantle as a result of the double subduction, is less affected by assimilation processes in the lower crust, the mode of phenocrysts could be smaller than for anhydrous magma under the same temperature conditions. This is because the liquidus temperature falls with increasing H2O contents in melt (e.g. 1100°C at 3 wt % H2O; Baker & Eggler, 1983). It may, therefore, be reasonable that a depleted primary magma enriched in water can cool without crystallization to a temperature similar to that of the isotopically enriched end-member (Fm1) of the high-T group. This consequently reduces the temperature difference between the aphyric-depleted and the enriched magma, making it difficult to form reverse zoning in pyroxenes. On the other hand, the low-T group represented by AK1102 commonly contains pyroxenes with reverse zoning and hornblendes surrounded by pyroxene reaction rims. Under low-pressure conditions, the hornblende stability field is sensitive to magmatic temperature, and it dissociates to form crystal aggregates of pyroxene and melts above 900–1000°C (Boettcher, 1977). Thus, the occurrence of both reaction rims in hornblende and reverse zoning of pyroxenes also indicates that the temperature of the enriched end-member (Fm2) was significantly lower than that for Pm and Fm1. Furthermore, the temperature of Fm2 was raised by mixing with Pm. It is, however, difficult to produce Fm2 directly from Fm1 by simple cooling processes associated with fractional crystallization in a shallow magma chamber, because of the significant difference in isotopic composition between them (Fig. 12). It is more likely that the high-T and low-T magmas behaved independently in the shallow region through the Ol stage. The temperature and isotopic differences between the high-T and low-T magmas could have depended on the difference in the degree of assimilation in the lower crust, before their migration into the shallow magma chamber. On the other hand, it may be possible to produce a high-T magma from a low-T magma, when a single magma batch such as low-temperature Fm2 is repeatedly injected by a ‘hotter’ primary magma (Pm) in the shallow magma chamber, resulting in rising magmatic temperature and increasing isotopic depletion. However, a relatively large discontinuity in isotopic composition and magmatic temperature exists between the high-T and low-T groups (Fig. 12), indicating that two enriched magma reservoirs have been involved in the formation of the Ol stage magmas. The relatively small variations of isotopic compositions and temperature in each group may have been caused by more complicated processes such as repeated injections of relatively small amounts of new hot magma into the enriched reservoirs in the shallow magma chamber.
It may be concluded that the isotopic variations and disequilibrium between glass and coexisting phenocrysts in the Ol stage were caused by two-component mixing of relatively depleted aphyric magma and enriched phenocryst-rich magma in the shallow magma chamber. These two magmas evolved independently in the lower crust before their mixing. This model is, therefore, different from that for the generation of calc-alkaline series magmas proposed by Sakuyama (1979, 1981), in which the original and its fractionated daughter magmas are mixed, so that isotopic variations are not expected. Our model is rather similar to the models of Eichelberger (1978) and Takahashi (1986) based on petrological arguments, proposing a mixing of magma derived from the mantle and a partially fused lower crust.
|
TEMPORAL ISOTOPIC VARIATION OF AKAGI VOLCANO |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTS AND DISCUSSIONINTERACTION BETWEEN ISLAND-ARC...EXISTENCE OF A WATER-RICH...EVOLUTION OF MAGMA IN...TEMPORAL ISOTOPIC VARIATION OF...CONCLUSIONSREFERENCES |
|---|
The isotope signatures changed temporally during the activity of Akagi volcano. The early stage volcanic rocks are isotopically most depleted and the final stage becomes depleted again, whereas the middle stage has more enriched signatures (Fig. 14). These temporal variations are probably attributable to the evolution of the interface (wall rock) between lower crust and the magma derived from the mantle wedge beneath Akagi. When magma starts to intrude into the lower crust, the margins of the magma conduit could be chilled as a result of the significant temperature difference between the basaltic magma and the lower crust. This could have led to less interaction between the lower crust and the magma or less assimilation by the magma resulting in a relatively depleted isotopic signature at the early stage of Akagi volcano. Repeated injection of magma should raise the temperature of the wall rock near the magma conduit and exceed the solidus temperatures of the surrounding wall-rock materials in the lower crust. Thus, the volcanic rocks of the middle stage could have acquired the relatively enriched isotopic signature by lower-crustal assimilation relative to magmas of the early stage (Fig. 14). In addition to the thermal evolution of the wall rock in the lower crust, the interface between the magma and the lower crust may have compositionally evolved, as a result of the solidification of magma and the formation of amphibole by the reaction between a water-rich primary magma and pyroxene in the lower crust (see discussion in previous section). Such an evolved interface could act as an insulator that would prevent interaction between the lower crust and the magma. Therefore, the isotope signature in the final stage would become similar to that of the early stage, with depleted isotope characteristics relative to the middle stages (Fig. 14).
|
|
CONCLUSIONS |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTS AND DISCUSSIONINTERACTION BETWEEN ISLAND-ARC...EXISTENCE OF A WATER-RICH...EVOLUTION OF MAGMA IN...TEMPORAL ISOTOPIC VARIATION OF...CONCLUSIONSREFERENCES |
|---|
The evolution of magma to form Akagi volcano is schematically summarized in Fig. 15.
|
The major element compositions and the phenocryst assemblages of Akagi volcano indicate that water content in the magma is significantly higher than that of other volcanoes in the volcanic front in the NE Japan arc. This may be attributed to the unique geometry of subduction beneath Akagi volcano. The Philippine Sea plate subducts into the mantle wedge above the downgoing Pacific plate underneath Akagi. Such double subduction of the oceanic plates should provide a larger amount of water as a result of the dehydration of hydrous minerals in the slabs, to be contributed to the source region of magmas for Akagi volcano, than in the case of single subduction of the Pacific plate, which produced other volcanoes in the volcanic front in NE Japan.
On the other hand, multi-isotope and trace element systematics of Akagi volcano clearly require interaction between the water-rich magma and the lower crust along with amphibole fractionation. Such assimilation–fractional crystallization in the lower crust may largely have been induced by a fall in the solidus temperature of the lower crust as a result of the higher water content in the primary magma of this volcano relative to that of other volcanoes in the volcanic front. The AFC model calculation suggests that the latent heat of crystallization and the heat of fusion for the wall rock were nearly balanced in the lower crust, implying that assimilation proceeded easily because of the high temperature in the lower crust and the high water content in the magma.
Isotope disequilibrium was observed between phenocrysts formed in the shallow magma chamber and the coexisting melts. Furthermore, the apparent magma temperatures correlate with Sr isotope compositions but to a lesser extent with Nd isotope compositions. These observations suggest that at least two recognizable and distinct magmas, which had been separately evolving in the lower crust, existed in the shallow magma chamber in a single stage. Then, the magmas were mixed independently with injections of isotopically depleted aphyric magmas. Although the other stages were not examined for phenocryst geochemistry in this study, it may be also possible, on the contrary, that the depleted magma is mixed by the injection of enriched magma. Such a difference should be attributable to the residence time of magma, and the evolution of the interface between the magma and the wall rock in the lower crust. Therefore, we conclude that the evolution of Akagi volcano has been strongly controlled by the assimilation process in the lower crust induced by the water-rich magma compositions.
As mentioned repeatedly in this paper, a primary factor in enhancing assimilation in the lower crust is the higher water content in the primary magma than that in magmas producing the surrounding volcanoes in the volcanic front. If this is the case, the along-arc Sr isotope variation, which decreases on either side of the ‘unique’ Akagi volcano in the volcanic front, may indicate a decrease in the extent of assimilation in the lower crust as a result of a decrease in the water content of the primary magmas. This may be due to the location of Akagi volcano, which stands above the top of a ‘tongue’ of the Philippine Sea plate. The position of this feature is shown by the isodepths to the subducting plate in Fig. 1. This would further increase the water content in the source region of the primary magma by the dehydration of the Philippine Sea plate, in addition to the water from the Pacific plate. If this is true, it may be possible to more clearly deduce the extent of the Philippine Sea plate in the mantle wedge of the Pacific plate and the role of water in the generation of island-arc magmas by extending detailed trace element and multi-isotope systematics to the volcanoes in central Japan as has been performed for Akagi volcano in this study.
|
ACKNOWLEDGEMENTS |
|---|
We thank T. Shibata, M. Yoshikawa and A. Makishima for their technical help in isotopic analyses and ICP-MS analysis. Thanks are also due to our many colleagues with whom we have discussed various issues during the course of this study. We are deeply indebted to Professor I. Kushiro for his encouragement. Professor I. Moriya is also sincerely thanked for geological information on the study area. Dr G. E. Bebout is also acknowledged for improving the manuscript. This research is supported by JSPS Fellowships for Japanese Junior Scientists (K.K.), and by the Ministry of Education, Science, Sports and Culture of the Japanese Government (Monbukagakusho) and the Japanese Society for the Promotion of Science (JSPS) (E.N.).
|
FOOTNOTES |
|---|
*Corresponding author. E-mail: katsura{at}pheasant.misasa.okayama-u.ac.jp
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONRESULTS AND DISCUSSIONINTERACTION BETWEEN ISLAND-ARC...EXISTENCE OF A WATER-RICH...EVOLUTION OF MAGMA IN...TEMPORAL ISOTOPIC VARIATION OF...CONCLUSIONSREFERENCES |
|---|
Adam, J. & Green, T. H. (1994). The effects of pressure and temperature on the partitioning of Ti, Sr and REE between amphibole, clinopyroxene and basanitic melts. Chemical Geology 117, 219–233.
Andersen, D. J., Lindsley, D. H. & Davidson, P. M. (1993). QUILF: a Pascal program to assess equilibria among Fe–Mg–Mn–Ti oxides, pyroxenes, olivine, and quartz. Computers and Geosciences 19, 1333–1350.
Arculus, R. J. & Johnson, R. W. (1981). Island-arc magma sources: a geochemical assessment of the roles of slab-derived components and crustal contamination. Geochemical Journal 15, 109–133.[Web of Science]
Baker, D. R. & Eggler, D. H. (1983). Fractionation paths of Atka (Aleutians) high-alumina basalts: constraints from phase relations. Journal of Volcanology and Geothermal Research 18, 387–404.[Web of Science]
Beard, J. S. & Lofgren, G. E. (1991). Dehydration melting and water-saturated melting of basaltic and andesitic greenstones and amphibolites at 1, 3, and 6·9 kb. Journal of Petrology
32, 365–401.
Beard, J. S., Abitz, R. J. & Lofgren, G. E. (1993). Experimental melting of crustal xenolith from Kilbourne Hole, New Mexico and implications for magma contamination and genesis. Contributions to Mineralogy and Petrology 115, 88–102.
Beard, J. S., Lofgren, G. E., Sinha, A. K. & Tollo, R. P. (1994). Partial melting of apatite bearing charnockite, granulite, and diorite: melt compositions, restite mineralogy, and petrologic implications. Journal of Geophysical Research 99, 21591–21603.
Ben Othman, D., White, W. M. & Patchett, J. (1989). The geochemistry of marine sediments, island arc magma genesis and crust–mantle recycling. Earth and Planetary Science Letters 94, 1–21.[Web of Science]
Bindeman, I., Davis, A. M. & Drake, M. J. (1998). Ion probe study of plagioclase–basalt partition coefficients at natural concentration levels of trace elements. Geochimica et Cosmochimica Acta 62, 1175–1193.
Boettcher, A. L. (1977). The role of amphiboles and water in circum-Pacific volcanism. In: Manghnani, M. H. & Akimoto, S. (eds) High-Pressure Research: Application in Geophysics . New York: Academic Press, pp. 107–125.
Cohen, R. S., Evensen, N. W., Hamilton, N. W. & O’Nions, R. K. (1980). The geochemistry of marine sediments, island arc magma genesis and crust–mantle recycling. Earth and Planetary Science Letters 94, 1–21.
Cousens, B. L., Allan, J. F. & Gorton, M. P. (1994). Subduction-modified pelagic sediments as the enriched component in back-arc basalts from the Japan Sea: Ocean Drilling Program Site 797 and 794. Contributions to Mineralogy and Petrology 117, 421–434.
DePaolo, D. J. (1981). Trace element and isotopic effects of combined wallrock assimilation and fractional crystallization. Earth and Planetary Science Letters 53, 189–202.[Web of Science]
DePaolo, D. J. & Wasserburg, G. J. (1979). Petrogenetic mixing models and Nd–Sr isotopic patterns. Geochimica et Cosmochimica Acta 43, 615–627.
Downes, H., Kempton, P. D., Briot, D., Harmon, R. S. & Leyreloup, A. F. (1991). Pb and O isotope systematics in granulite facies xenoliths, French Massif Central: implication for crustal process. Earth and Planetary Science Letters 102, 342–357.
Dupre, B. & Allègre, C. J. (1980). Sr–Nd–Pb isotopic correlation and the chemistry of the North Atlantic mantle. Nature 286, 17–22.
Eichelberger, J. C. (1978). Andesite volcanism and crustal evolution. Nature 275, 21–27.
Forden, J. D. & Green, D. H. (1992). Possible role of amphibole in the origin of andesite: some experimental and natural evidence. Contributions to Mineralogy and Petrology 109, 470–493.
Fujinawa, A. (1988). Tholeiitic and calc-alkaline magma series at Adatara volcano, Northeast Japan: 1. Geochemical constraints on their origin. Lithos , 22, 135–158.
Fujinawa, A. (1990). Tholeiitic and calc-alkaline magma series at Adatara volcano, Northeast Japan: 2. Mineralogy and phase relation. Lithos 24, 217–236.
Fujinawa, A. (1992). Distinctive REE patterns for tholeiitic and calc-alkaline magma series co-occurring at Adatara volcano, Northeast Japan. Geochemical Journal 26, 395–409.
Gill, J. B. (1981). Orogenic Andesites and Plate Tectonics . Berlin: Springer, 390 pp.
Green, T. H. (1994). Experimental studies of trace-element partitioning applicable to igneous petrogenesis—Sedona 16 years later. Chemical Geology 117, 1–36.[Web of Science]
Gust, D. A., Arculus, R. J. & Kersting, A. B. (1997). Aspects of magma sources and processes in the Honshu arc. Canadian Mineralogist 35, 347–365.
Hildreth, W. & Moorbath, S. (1988). Crustal contribution to arc magmatism in the Andes of Central Chile. Contributions to Mineralogy and Petrology 98, 455–489.[Web of Science]
Horio, A. & Umino, S. (1995). Mushy magma chamber beneath island arc volcanoes—evidence from Yunokuchi Pumice, Akagi volcano (in Japanese). Bulletin of Volcanological Society of Japan 40, 375–393.
Iizuka, Y. (1996). Experimental study on the slab–mantle interactions in subduction zones and its implications for the material recycling through the subduction zones. Ph.D. thesis, Okayama University, 232 pp.
Ishida, M. (1992). Geometry and relative motion of the Philippine Sea plate and Pacific plate beneath the Kanto–Tokai District, Japan. Journal of Geophysical Research 97, 489–513.
Ishikawa, T. & Nakamura, E. (1994). Origin of the slab component in arc lavas from across-arc variation of B and Pb isotopes. Nature 370, 205–208.
Kawano, Y., Yagi, K. & Aoki, K. (1961). Petrography and petrochemistry of the volcanic rocks of Quaternary volcanoes of northeastern Japan. Scientific Reports of Tohoku University, Series III 7, 1–46.
Kempton, P. D., Harmon, R. S., Hawkesworth, C. J. & Moorbath, S. (1990). Petrology and geochemistry of lower crustal granulites from the Geronimo Volcanic Field, southeastern Arizona. Geochimica et Cosmochimica Acta 54, 3401–3426.
Kersting, A. B., Arculus, R. J. & Gust, D. A. (1996). Lithospheric contribution to arc magmatism: isotope variations along strike in volcanoes of Honshu, Japan. Science 272, 1464–1468.[Abstract]
Koga, S. (1981). 14C age of pyroclastic flow-like deposits distributed at the northeastern foot of Akagi volcano, Gunma Prefecture (in Japanese). Bulletin of Volcanological Society of Japan 26, 239–240.
Koga, S. (1984). Geology and petrology of Akagi volcano, Gunma Prefecture, Japan. Scientific Reports, Institute of Geoscience, University of Tsukuba, Section B 5, 1–67.
Koide, Y. & Nakamura, E. (1990). Lead isotope analysis of standard rock samples. Mass Spectroscopy 38, 241–252.
Le Roex, A. P., Dick, H. J. B., Erlank, A. J., Reid, A. M., Frey, F. A. & Hart, S. R. (1983). Geochemistry, mineralogy and petrogenesis of lavas erupted along the southwest Indian Ridge between the Bouvet triple junction and 11 degree east. Journal of Petrology 24, 267–318.
Makishima, A. & Nakamura, E. (1991a). Precise measurement of cerium isotope composition in rock samples. Chemical Geology 94, 1–11.
Makishima, A. & Nakamura, E. (1991b). Calibration of Faraday cup efficiency in a multicollector mass spectrometer. Chemical Geology 94, 105–110.[Web of Science]
Makishima, A. & Nakamura, E. (1997). Suppression of matrix effect in ICP-MS by high power operation of ICP: application to precise determination of Rb, Sr, Y, Cs, Ba, REE, Pb, Th and U at ng g-1 level in a few milligram silicate sample. Geostandards Newsletter 21, 307–319.
McDonough, W. F. & Sun, S.-s. (1995). The composition of the Earth. Chemical Geology 120, 223–253.[Web of Science]
Miller, M. M., Goldstein, S. L. & Langmuir, C. H. (1994). Cerium/lead and lead isotope ratios in arc magmas and the enrichment of lead in the continents. Nature 368, 514–519.
Miyashiro, A. (1974). Volcanic rock series in island arcs and active continental margins. American Journal of Science 274, 321–355.[Abstract]
Moriya, I. (1968). Geomorphology and Geology of Akagi Volcano (in Japanese). Maebashi Forestry Bureau, pp. 1–65.
Nakamura, E. & Kushiro, I. (1997). Diffusivities of rare earth elements and Ba in magmatic silicate melts at high pressure. Proceedings of the Japan Academy 73B, 44–47.
Nakamura, E. & Kushiro, I. (1998). Trace element diffusion in jadeite and diopside melts at high pressures and its geochemical implication. Geochimica et Cosmochimica Acta 62, 3151–3160.
Nakamura, E., Campbell, I. H. & Sun, S.-s. (1985). The influence of subduction processes on the geochemistry of Japanese alkaline basalt. Nature 316, 55–58.
Nohda, S. & Wasserburg, G. J. (1981). Nd and Sr isotopic study of volcanic rocks from Japan. Earth and Planetary Science Letters 52, 264–272.
Notsu, K. (1983). Strontium isotope composition in volcanic rocks from the Northeast Japan Arc. Journal of Volcanology and Geothermal Research 18, 531–548.
Notsu, K., Kita, I. & Yamaguchi, T. (1985). Mantle contamination under Akagi volcano, Japan, as inferred from combined Sr–O isotope relationships. Geophysical Research Letters 12, 365–368.
O’Nions, R. K., Carter, S. R., Cohen, R. S., Evensen, N. W. & Hamilton, P. J. (1978). Pb, Nd and Sr isotopes in oceanic ferromanganese deposits and ocean floor basalts. Nature 273, 435–438.
Perfit, M. R., Gust, A. E., Bence, A. E., Arculus, R. J. & Taylor, S. R. (1979). Chemical characteristics of island-arc basalt: implications for mantle sources. Chemical Geology 30, 227–256.
Rudnick, R. L. & Fountain, D. M. (1995). Nature and composition of the continental crust: a lower crustal perspective. Reviews of Geophysics 33, 267–309.[Web of Science]
Rudnick, R. L. & Goldstein, S. L. (1990). The Pb isotopic compositions of lower crustal xenoliths and the evolution of lower crustal Pb. Earth and Planetary Science Letters 98, 192–207.[Web of Science]
Rudnick, R. L., McDonough, W. F., McCulloch, M. T. & Taylor, S. R. (1986). Lower crustal xenoliths from Queenland, Australia: evidence for deep crustal assimilation and fractionation of continental basalts. Geochimica et Cosmochimica Acta 50, 1099–1115.
Ryan, J. G., Morris, J., Tera, F., Leeman, W. P. & Tsvetkov, A. (1995). Cross-arc geochemical variations in the Kurile arc as a function of slab depth. Science
270, 625–627.
Sakuyama, M. (1979). Evidence of magma mixing: petrological study of Shirouma–Oike calc-alkaline andesite volcano, Japan. Journal of Volcanology and Geothermal Research 5, 179–208.
Sakuyama, M. (1981). Petrological study of the Myoko and Kurohime volcanoes, Japan: crystallization sequence and evidence for magma mixing. Journal of Petrology
22, 553–583.
Sakuyama, M. (1983a). Phenocryst assemblages and H2O contents in circum-Pacific arc magmas. In: Hilde, T. W. C. & Uyeda, S. (eds) Geodynamics of the Western Pacific–Indonesian Region. American Geophysical Union, Geodynamics Series 11, 143–158.
Sakuyama, M. (1983b). Petrology of arc volcanic rocks and their origin by mantle diapir. Journal of Volcanology and Geothermal Research 18, 297–320.
Sakuyama, M. & Nesbitt, R. W. (1986). Geochemistry of the Quaternary volcanic rocks of the Northeast Japan arc. Journal of Volcanology and Geothermal Research 29, 413–450.
Shibata, T. & Nakamura, E. (1991). Interaction between subducted oceanic slab and wedge mantle inferred from across-arc variations of Pb, Sr and Nd isotopic compositions in Northeastern Japan. Proceedings of the Japan Academy 67B, 115–120.
Shibata, T. & Nakamura, E. (1997). Across-arc variations of isotope and trace element compositions from Quaternary basaltic volcanic rocks in northeastern Japan: implications for interaction between subducted oceanic slab and mantle wedge. Journal of Geophysical Research 102, 8051–8064.
Shimoda, G., Tatsumi, Y., Nohda, S., Ishizaka, K. & Jahn, B. M. (1998). Setouchi high-Mg andesites revisited: geochemical evidence for melting of subducting sediments. Earth and Planetary Science Letters 160, 479–492.
Sisson, T. W. (1994). Hornblende–melt partitioning measured by ion microprobe. Chemical Geology 117, 331–344.
Stolz, A. J. & Davies, G. R. (1989). Metasomatised lower crustal and upper mantle xenoliths from north Queensland: chemical and isotopic evidence bearing on the composition and source of the fluid phase. Geochimica et Cosmochimica Acta 53, 649–660.
Takahashi, E. (1986). Genesis of calc-alkali andesite magma in a hydrous mantle–crust boundary: petrology of lherzolite xenoliths from the Ichinomegata crater, Oga peninsula, Northeast Japan, Part II. Journal of Volcanology and Geothermal Research 29, 355–395.
Tatsumi, Y., Hamilton, D. L. & Nesbitt, R. W. (1986). 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 29, 293–309.[Web of Science]
Umino, S. & Horio, A. (1998). Multistage magma mixing revealed in phenocryst zoning of the Yunokuchi Pumice, Akagi volcano, Japan. Journal of Petrology 39, 101–124.
Utsu, T. (1974). Epicenter distribution around Japan (in Japanese). Kagaku 44, 739–746.
White, W. M. & Hofmann, A. W. (1982). Sr and Nd isotope geochemistry of oceanic basalts and mantle evolution. Nature 296, 821–825.
White, W. M. & Patchett, J. (1984). Hf–Nd–Sr isotopes and incompatible element abundances in island arcs: implications for magma origins and crust–mantle evolution. Earth and Planetary Science Letters 67, 167–185.
White, W. M. & Dupre, B. (1986). Sediment subduction and magma genesis in the Lesser Antilles: isotopic and trace element constraints. Journal of Geophysical Research 91, 5927–5941.
Witt-Eickschen, G. & Harte, B. (1994). Distribution of trace elements between amphibole and clinopyroxene from mantle peridotites of the Eifel (western Germany): an ion-microprobe study. Chemical Geology 117, 235–250.
Wood, D. A., Joron, J. L. & Treuil, M. (1979). A re-appraisal of use of trace elements to classify and discriminate between magma series erupted in different tectonic setting. Earth and Planetary Science Letters 45, 326–336.
Woodhead, J. D. & Fraser, G. D. (1985). Pb, Sr and 10Be isotopic studies of volcanic rocks for magma genesis and crustal recycling in the Western Pacific. Geochimica et Cosmochimica Acta 49, 1925–1930.
Woodhead, J. D., Eggins, S. M. & Johnson, R. W. (1998). Magma genesis in the New Britain island arc: further insights into melting and mass transfer processes. Journal of Petrology 39, 1641–1668.
Yamaguchi, T. (1990). Comparative petrology of the Hotaka and Akagi volcanoes. Journal of the Japanese Association of Mineralogists, Petrologists and Economic Geologists 85, 229–248.
Yoshii, T. (1979). A detailed cross-section of the deep seismic zone beneath northeastern Honshu, Japan. Tectonophysics 55, 349–360.
Yoshikawa, M. & Nakamura, E. (1993). Precise isotopic determination of trace amounts of Sr in magnesium-rich samples. Journal of the Japanese Association of Mineralogists, Petrologists and Economic Geologists 88, 548–561.
Zartman, R. E. & Doe, B. R. (1981). Plumbotectonics—the model. Tectonophysics 75, 135–162.[Web of Science]
Zartman, R. E. & Haines, S. M. (1988). The plumbotectonic model for Pb isotopic systematics among major terrestrial reservoirs—a case for bi-directional transport. Geochimica et Cosmochimica Acta 52, 1327–1399.
Zhao, D. (1992). Seismic velocity structure of the crust beneath the Japan islands. Tectonophysics 212, 289.
Zhao, D., Horiuchi, S. & Hasegawa, A. (1990) 3-D seismic velocity structure of the crust and the uppermost mantle in the northeastern Japan Arc. Tectonophysics 181, 1–15.
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  J.-I. KIMURA and T. YOSHIDA Contributions of Slab Fluid, Mantle Wedge and Crust to the Origin of Quaternary Lavas in the NE Japan Arc J. Petrology, November 1, 2006; 47(11): 2185 - 2232. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J. DUFEK and G. W. BERGANTZ Lower Crustal Magma Genesis and Preservation: a Stochastic Framework for the Evaluation of Basalt-Crust Interaction J. Petrology, November 1, 2005; 46(11): 2167 - 2195. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
P. D. CLIFT, G. D. LAYNE, Y. M. R. NAJMAN, A. KOPF, N. SHIMIZU, and J. HUNT Temporal Evolution of Boron Flux in the NE Japan and Izu Arcs Measured by Ion Microprobe from the Forearc Tephra Record J. Petrology, July 1, 2003; 44(7): 1211 - 1236. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 Y. Tamura Some geochemical constraints on hot fingers in the mantle wedge: evidence from NE Japan Geological Society, London, Special Publications, January 1, 2003; 219(1): 221 - 237. [Abstract] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||