New Insights into Andesite Genesis: the Role of

JOURNAL OF PETROLOGY
VOLUME 49
NUMBER 11
PAGES 1971^2008
2008
doi:10.1093/petrology/egn055
New Insights into Andesite Genesis: the Role of
Mantle-derived Calc-alkalic and Crust-derived
Tholeiitic Melts in Magma Differentiation
beneath Zao Volcano, NE Japan
Y. TATSUMI1*, T. TAKAHASHI1, Y. HIRAHARA1, Q. CHANG1,
T. MIYAZAKI1, J.-I. KIMURA1, M. BAN2 AND A. SAKAYORI3
1
INSTITUTE FOR RESEARCH ON EARTH EVOLUTION (IFREE), JAPAN AGENCY FOR MARINE^EARTH SCIENCE AND
TECHNOLOGY (JAMSTEC), YOKOSUKA 237-0061, JAPAN
2
DEPARTMENT OF EARTH AND ENVIRONMENTAL SCIENCES, YAMAGATA UNIVERSITY, YAMAGATA 990-8560, JAPAN
3
DEPARTMENT OF EARTH SCIENCES, KANAZAWA UNIVERSITY, KANAZAWA 920-1192, JAPAN
RECEIVED APRIL 5, 2008; ACCEPTED OCTOBER 10, 2008
ADVANCE ACCESS PUBLICATION DECEMBER 4, 2008
Two distinctive differentiation trends, tholeiitic and calc-alkalic,
are recognized in Zao volcano, which is located immediately behind
the volcanic front of the NE Japan arc. The genetic relationship
between these two magma series is critical for a better understanding
of andesite genesis, because they often coexist in close spatial and
temporal proximity in arc volcanoes. Petrographic features indicative
of ‘disequilibrium’, such as reversely zoned pyroxene phenocrysts,
the wide and bimodal compositional distribution in Ca/(Ca þ Na)
of plagioclase phenocrysts, honeycomb textures and dusty zones that
these plagioclase phenocrysts often exhibit, and the presence of olivine^pyroxene pairs with different Mg/Fe, are observed exclusively
in calc-alkalic rocks. In tholeiitic rocks the Sr isotopic ratios of
plagioclase phenocrysts, determined by both micromilling combined
with thermal ionization mass spectrometry, and laser-ablation
inductively coupled plasma mass spectrometry techniques, are constant at 07042^07044. On the other hand, those in calc-alkalic
rocks (07033^07042) show more complex characteristics, which
can be best understood if at least three end-member components,
a calc-alkalic basaltic melt, a tholeiitic basaltic melt and a tholeiitic
felsic melt, contribute to the production of mixed calc-alkalic
magmas. The 87Sr/86Sr and trace element compositions of the leastdifferentiated basalt magmas, which are inferred from the composition of the calcic plagioclase [Ca/(Ca þ Na) 409], suggest that
two types of basaltic magma, calc-alkalic and tholeiitic, exist
beneath the volcano. The tholeiitic basalt magma has a higher
*Corresponding author. E-mail: [email protected]
87
Sr/86Sr than the calc-alkalic magma (07042 vs 07038) and a
characteristic trace element signature consistent with the presence
of plagioclase and amphibole as melting residues. This suggests that
the tholeiitic magmas are produced via anatexis of amphibolitic crust
caused by underplating and/or intrusion of mantle-derived calcalkalic basalt magmas into the sub-Zao crust. The mantle-derived
calc-alkalic basalt magma mixes with crust-derived tholeiitic melts
to form calc-alkalic andesite magmas. The hypothesis proposed here
requires revision (or even abandonment) of the general consensus
that calc-alkalic magmas have greater contributions of a crustal
component than tholeiitic magmas.
KEY WORDS:
andesite; calc-alkalic: crust; mantle; tholeiitic
I N T RO D U C T I O N
How andesite is generated has long been a central question
of igneous petrology. The reason for this is two-fold. First,
andesite erupts in more than 80% of arc volcanoes, typifies
subduction zone magmatism that creates over 20% of current terrestrial magmatic products, and is the dominant
volcanic rock in mature continental arcs. Second, the continental crust, the most differentiated end-member among
components within the solid Earth, is overall andesitic or
ß The Author 2008. Published by Oxford University Press. All
rights reserved. For Permissions, please e-mail: journals.permissions@
oxfordjournals.org
JOURNAL OF PETROLOGY
VOLUME 49
intermediate in composition (e.g. Rudnick, 1995; Taylor,
1995; Rudnick & Gao, 2003). Knowledge of andesite genesis should therefore provide key constraints on the origin
of continental crust and differentiation processes during
the evolution of the solid Earth. The following five models
are currently favored for the production of andesitic
magmas (sensu lato).
(1) Crystallization differentiation of mantle-derived
basaltic magmas either in shallow crustal reservoirs or in
the deep crust close to the Moho (Sisson & Grove, 1993;
Mu«ntener et al., 2001; Annen & Sparks, 2002; Pichavant
et al., 2002; Prouteau & Scaillet, 2003).
(2) Anatexis of mafic lower crust by intrusion or underplating of mantle-derived basaltic magma (Takahashi,
1986; Smith & Leeman, 1987; Petford & Atherton, 1996;
Kimura et al., 2002; Annen et al., 2006; Tatsumi et al., 2008).
(3) Open-system differentiation, such as mixing
between felsic and mafic magmas, and crustal assimilation
(Eichelberger, 1975; DePaolo, 1981; Sakuyama, 1981;
Hildreth & Moorbath, 1988; Clynne, 1999; Tatsumi &
Kogiso, 2003; Dungan & Davidson, 2004).
(4) Generation of melts in equilibrium with mantle
peridotite under hydrous conditions as a result of either
direct fluxing of slab-derived fluids or slab-melting and
subsequent melt^mantle interaction (Kay, 1978; Tatsumi,
1981; Crawford et al., 1989; Pearce et al., 1992; Yogodzinski
et al., 1994; Kelemen, 1995; Blatter & Carmichael, 2001;
Tatsumi & Hanyu, 2003; Parman & Grove, 2004).
(5) Production of andesitic to more felsic melts by
dehydration melting of the subducted oceanic crust
(Kay, 1978; Martin, 1986; Stern & Kilian, 1986; Defant &
Drummond, 1990).
Two distinctive differentiation trends, tholeiitic and
calc-alkalic, are recognized in the sub-alkalic volcanic
rocks, denoting the presence or absence of relative iron
enrichment during magmatic differentiation (Wager &
Deer, 1939; Nockolds & Allen, 1953; Kuno, 1959; Irvine &
Baragar, 1971; Miyashiro, 1974). It has been well documented that tholeiitic rocks are dominant in juvenile oceanic
arcs, whereas calc-alkalic rocks are the major magmatic
products in mature continental arcs with thicker crust
(e.g. Miyashiro, 1974; Ewart, 1982). In several arc^trench
systems, however, tholeiitic and calc-alkalic andesites to
dacites have been observed to coexist in close temporal
and spatial proximity (Kuno, 1959; Sakuyama, 1981; Grove
& Baker, 1984; Fujinawa, 1988; Brophy, 1990; Hunter &
Blake, 1995; Hunter, 1998). Resolution of the genetic relationship between these two types of andesitic magmas
should, therefore, provide a better understanding of andesite genesis and arc crust evolution.
The primary aim of this paper is to investigate the
mechanism that produces tholeiitic and calc-alkalic
magmas at a single volcano, by combining petrographical
and geochemical data, including micro-analyses of the
NUMBER 11
NOVEMBER 2008
isotopic and trace element compositions of plagioclase
phenocrysts.
O V E RV I E W O F T H O L E I I T I C A N D
CA LC -A L K A L IC A N DESI T E
Definition and chemical characteristics
Bowen (1928) proposed that silica content increases and
iron content decreases with differentiation of sub-alkalic
magmas. On the other hand, Fenner (1929) emphasized
that some suites show iron enrichment during magmatic
differentiation. These two differentiation trends have been
referred to as calc-alkalic and tholeiitic, respectively
(Wager & Deer, 1939). The identification of these two
trends is commonly based on a ternary plot of Na2O þ
K2O, FeO (total iron as FeO), and MgO (Irvine &
Baragar, 1971). For a more quantitative distinction of the
two magma series, FeO/MgO vs SiO2 variation plots
(Miyashiro, 1974) are often used; calc-alkalic and tholeiitic
rock series show steeper and gentler slopes, respectively,
than the straight line: SiO2 (wt %) ¼ 64 FeO/
MgO þ 428. Unfortunately, however, Miyashiro’s discriminant line is frequently misused and applied as a simple
compositional discriminant outside this range rather than
its intended use as a ‘trend slope’ comparison.
The terms medium-K and calc-alkalic series are often
used interchangeably (Fig. 1). This misusage may be due to
poor understanding of the following two observations.
First, concentrations of incompatible elements such as K,
Rb, and Nb in lavas at a constant SiO2 content increase
with distance from the volcanic front or the height of a
volcano above the slab surface (Fig. 1), which is generally
known as the K^h relationship (Dickinson, 1975). Second,
calc-alkalic series rocks are generally more enriched
in incompatible elements than tholeiitic series rocks
(e.g. Masuda & Aoki, 1979; Kimura & Yoshida, 2006).
A schematic illustration showing the across-arc variation
in magma series from tholeiitic, via calc-alkalic, to highK or alkalic with distance from the volcanic front (Fig. 1a)
is then often cited (e.g. Hess, 1989; Wilson, 1989), although
calc-alkalic rocks do exist along the volcanic front and tholeiitic rocks belonging to the medium-K series often erupt
at volcanoes behind the volcanic front as shown in Fig. 1b
(Kuno, 1960; Kawano et al., 1961; Gill, 1981; Yoshida & Aoki,
1984; Tatsumi & Eggins, 1995; Tatsumi & Kogiso, 2003).
The tholeiitic and calc-alkalic series should be defined by
the presence and absence of iron enrichment, respectively,
whereas low-, medium-, and high-K series should be
defined on the basis of K2O concentrations.
Yoder & Tilley (1962) emphasized the presence of at
least two types of basalt magmas, one with normative
hypersthene differentiating to silica-saturated and the
other with normative nepheline differentiating to silicaundersaturated liquids, and coined the terms olivine
1972
TATSUMI et al.
(a)
CALC-ALKALIC VS THOLEIITIC SERIES
An awkward scheme
(b)
4
An unfailing scheme
4
Shoshonitic
(Alkalic)
Shoshonitic
(Alkalic)
Hig
h-K
3
K2O (wt.%)
2
lic
-alka
Calc
2
Med
Low-K
Low
60
65
0
50
70
55
SiO2 (wt.%)
Volcanic
Front
Forearc
Volcanic Arc
Medium-K
(CA/TH)
High-K(CA/TH)
Increasing K, H2O
Trench
Calc-alkalic
High-K
Backarc
Backarc
Forearc
Volcanic Arc
Mantle Wedge
s
tho
ere
ere
ph
i
s
tho
i
dL
te
uc
bd
Su
70
Mantle Wedge
ph
dL
65
SiO2 (wt.%)
Low-K tholeiitic
Increasing K, H2O
60
Trench
55
K
ium-
1
iitic)
-K (thole
0
50
h-K
Volcanic
Front
1
Hig
Low-K (CA/TH)
K2O (wt.%)
3
te
uc
bd
Su
Fig. 1. Two schemes for the classification and spatial distribution of arc magmas. (a) The interchangeable usage of the terms medium-K and
calc-alkalic series leads to an awkward understanding of across-arc variation; that is, magma types change from low-K tholeiitic via calc-alkalic
to high-K series. (b) Low-, medium-, and high-K series should be defined by K2O contents, whereas calc-alkalic and tholeiitic series should be
identified based on the absence and presence of iron enrichment, respectively. Concentrations of incompatible elements such as K increase
towards the back-arc side of a volcanic arc. Calc-alkalic rocks (CA) do occur in the low-K zone along the volcanic arc and tholeiitic rocks
(TH) often coexist with calc-alkalic rocks in volcanoes behind the volcanic front.
tholeiite and alkali basalt for these magmas, respectively.
Tholeiites can be clearly defined based on normative compositions; however, it should be noted that Yoder & Tilley’s
tholeiites are identical to sub-alkalic magmas and thus
include both tholeiitic and calc-alkalic rock series defined
based on the presence or absence of iron enrichment
during differentiation of sub-alkalic magmas.
Calc-alkalic and tholeiitic series can be well identified
for rocks with intermediate compositions; in other words,
it may be difficult to classify mafic (basaltic) or felsic
(rhyolitic) rocks into these two rock series. We tentatively
use the term calc-alkalic or tholeiitic for basalts and rhyolites, if they, together with andesites, form calc-alkalic or
tholeiitic trends, respectively.
Petrographical characteristics
The calc-alkalic vs tholeiitic series should be defined exclusively on the basis of differentiation trends. However, the
identification of such chemical trends is, in some cases,
difficult because of a lack of sufficient data for defining the
1973
JOURNAL OF PETROLOGY
VOLUME 49
chemical trend. As a result, this can be supplemented with
petrographical characteristics to identify these two magma
series. On the basis of groundmass mineralogy Kuno
(1950, 1959, 1968) divided sub-alkalic volcanic rocks into
two series: hypersthenic and pigeonitic. These are distinguished by the presence or absence of orthopyroxene in
the groundmass, and are synonymous with the calc-alkalic
and tholeiitic series, respectively. Although it may not be
easy to identify fine-grained groundmass pyroxenes under
the microscope, orthopyroxene phenocrysts with a reaction rim of clinopyroxene, which occur solely in the pigeonitic rock series, can be recognized (e.g. Kawano et al.,
1961). Furthermore, phenocrysts of hornblende and biotite
are limited to the hypersthenic rock series (Kuno, 1950).
It has been well established that Kuno’s scheme is valid for
Quaternary arc volcanoes along the trench-side volcanic
chain in the NE Japan and Izu^Bonin arcs (Kuno, 1950;
Kawano et al., 1961; Wada, 1981, 1985; Fujinawa, 1988, 1990),
where tholeiitic rocks are broadly equivalent to the low-K
series of Gill (1981).
Sakuyama (1981) examined the petrographical characteristics of volcanic rocks from Myoko^Kurohime volcanoes, Central Japan, where both calc-alkalic and tholeiitic
magmas have erupted from single vents, and divided these
rocks into two types, N-type and R-type, on the basis of the
absence and presence of reversely zoned mafic phenocrysts,
respectively. The following ‘disequilibrium’ petrographical
features characterize the R-type volcanic rocks:
(1) the presence of reversely zoned pyroxene phenocrysts
with a lower Mg-number [¼100 Mg/(Mg þ Fe)] core
surrounded by a higher Mg-number rim;
(2) the presence of groundmass pyroxenes with higher
Mg-number than the phenocryst cores;
(3) bimodal distribution in the core compositions of
plagioclase phenocrysts;
(4) disequilibrium phenocryst assemblages such as
Mg-rich olivine and quartz;
(5) patchy groundmass with different colors and/or
amount of mafic minerals.
These disequilibrium features are not observed in
N-type rocks. It was emphasized by Sakutyama (1981)
that the N- and R-type rocks are broadly equivalent to
Kuno’s pigeonitic and hypersthenic rock series, and hence
more generally to tholeiitic and calc-alkalic series,
respectively.
The consistency between Sakuyama’s petrographical
classification and bulk-rock chemical characteristics (i.e.
N- vs R-type and tholeiitic vs calc-alkalic, respectively)
has been well established at least for Quaternary NE
Japan arc volcanoes (Wada, 1981; Sakuyama, 1983;
Fujinawa, 1988, 1990). However, applied elsewhere there
are exceptions. Particular calc-alkalic andesites (high-Mg
andesites) from the Setouchi volcanic belt, SW Japan,
characterized by unusually high MgO contents or high
NUMBER 11
NOVEMBER 2008
Mg-number and hence believed to represent leastdifferentiated mantle-derived magmas, are petrographically classified as N-type rocks (Tatsumi, 2006). On the
other hand, high-Mg andesite from Mt. Shasta, USA,
which is considered as representative of primitive andesite
(Baker et al., 1994; Grove et al., 2002), has now been identified as R-type andesite, and hence is likely to be the product of magma mixing (Streck et al., 2007).
Occurrence
It has long been known that some volcanic arcs are characterized by either tholeiitic or calc-alkalic magmatism
(Jakes & Gill, 1970; Plank & Langmuir, 1988). Calc-alkalic
rocks are clearly dominant in continental arcs rather than
oceanic arcs (Miyashiro, 1974; Ewart, 1982), and calcalkalic/tholeiitic volume ratios tend to increase with
increasing age or arc maturity (Baker, 1973) and crustal
thickness (Gill, 1981). It should be stressed, however, that
tholeiitic and calc-alkalic magmas do coexist in some
single volcanic systems; e.g. Mt. Shasta, USA (Baker et al.,
1994), Chichontepec, El Salvador (Bau & Knittel, 1993),
Aso in SW Japan (Hunter, 1998), and Myoko^Kurohime in
Central Japan (Sakuyama, 1981). Furthermore, along the
volcanic front of the NE Japan arc about one-third of
Quaternary volcanoes erupt both tholeiitic and calcalkalic rocks (Kawano et al.,1961). Examining the geochemistry of the two coexisting magma series in the above
volcanoes reveals that the calc-alkalic rocks are generally
more enriched in both compatible elements such as Mg, Ni
and Cr and incompatible elements such as Rb, K, Th and
U (e.g. Masuda & Aoki,1979; Kimura & Yoshida, 2006).
The misunderstanding caused by the interchangeable
usage of the terms medium-K and calc-alkalic series
(Fig. 1a) may further mislead us to a ‘likely’ conclusion
that calc-alkalic magmas are more hydrous than tholeiitic
magmas, because the H2O content in mantle-derived primary magmas may increase together with incompatible
elements towards the back-arc side within a volcanic arc
(Sakuyama, 1979; Tatsumi et al., 1983). The experimental
constraints that calc-alkalic trends can be reproduced
under hydrous conditions but tholeiitic trends under
H2O-poor conditions (e.g. Grove & Baker, 1984; Sisson &
Grove, 1993; Hamada & Fujii, 2008) may further reinforce
this idea. However, it should be stressed here that there are
no data available that demonstrate the difference in H2O
contents in these magma series occurring in a single
volcano.
G E O L O G I C A L B AC KG RO U N D
O F Z AO VO L C A N O
NE Japan arc
The NE Japan arc is formed by the subduction of the
Pacific plate beneath the North American and Eurasian
1974
TATSUMI et al.
44°N
Eurasia
Plate
CALC-ALKALIC VS THOLEIITIC SERIES
(3) the number of volcanoes and the eruptive volume are
greater in the trench-side volcanic chain;
(4) volcanic rocks in the back-arc side chain tend to be
more enriched in incompatible elements than those in the
trench-side chain.
These characteristics are well documented in the central
to northern NE Japan arc, but not in the southern arc
where the Philippine Sea plate is being subducted into
the mantle wedge above the descending Pacific plate and
the Izu^Bonin^Mariana arc is colliding with the Japanese
islands (Fig. 2). This complicated tectonic setting may cause
the unusual characteristics of the magmatism at the
southern end of the NE Japan arc.
From north to south the pre-Tertiary basement rocks of
the NE Japan arc consist of Cretaceous to Jurassic metamorphic rocks, Cretaceous sedimentary rocks and granitoids, and subduction zone and metamorphic complexes
of Ordovician to Cretaceous age. Granitoids in these basements show spatial variations in isotopic composition, with
more enriched isotopic signatures southwards (Kagami,
2005). Kersting et al. (1996) and Kimura & Yoshida (2006)
further demonstrated the correlation between the Sr^Nd
isotopic compositions of volcanic front lavas and their
underlying basement rocks.
km
200
m
0k
15
N American
Plate
J a p
a n
A r c
m
rc
0k eA
10 Kuril
60 km
40 km
20
km
ana
in-Mari
Izu-Bon
n
Pacific
Plate
Zao volcano
Arc
Philippine Sea
Plate
32°N
136°E
o
Ja
e-
pa
I
II
Zon
36°N
N
Zo
ne
-II
Zao Volcano
Tr
enc
h
rt
h e
a s
t
Z o n
e - I
40°N
140°E
144°E
Fig. 2. Distribution of Quaternary volcanoes in the NE Japan arc
and the adjacent Kurile and Izu^Bonin^Mariana arcs. The NE
Japan volcanic arc can be subdivided into three zones (filled circles,
Zone I; half-filled circles, Zone 2; open circles, Zone 3) based on
Sr^Nd isotopic characteristics and basement geology (Kimura &
Yoshida, 2006). Zao volcano straddles the volcanic front, which lies
110 km above the top of the subducting Pacific Plate (dashed lines,
depth contours to the slab surface; numbers, depths). Small filled
circles indicate volcanoes in the Kurile and Izu^Bonin^Mariana arcs.
plates at the Japan Trench (Fig. 2) at a rate of 10 cm/year.
The volcanic front of this arc, which is the trenchward
boundary of the volcanic arc, runs parallel to the trench
axis. The NE Japan arc exhibits the following tectonic and
magmatic characteristics, which are common to most arc^
trench system (Tatsumi & Eggins, 1995):
(1) dual volcanic chains are present, one trench-side and
one on the back-arc side, known as the Nasu and Chokai
chains, respectively;
(2) these chains are 110 km and 150^170 km above the
top of the subducted oceanic lithosphere, respectively, at
least in the central portion of the arc (Fig. 2);
A total of 55 volcanoes are distributed in the NE Japan
arc. Zao volcano is one such Quaternary volcano, and is
situated immediately behind the volcanic front in the
Nasu trench-side chain (Fig. 2). The following summary
of the geological history of Zao volcano is based on a
synthesis of previously published studies (Oba & Konda,
1989; Takaoka et al., 1989; Sakayori, 1992).
The basement around Zao consists of Cretaceous
granitoids of the Abukuma Terrane and Tertiary volcanic
rocks. Magmatic activity, which commenced at 10 Ma
and has formed a volcanic edifice with a current volume
of 25 km3, can be separated into four major stages
(Fig. 3): Stage 1 (10^06 Ma): formation of a small volcano
consisting of basaltic to basaltic andesite pyroclastic
rocks and dykes; Stage 2 (03 Ma): formation of a stratovolcano composed of andesitic to dacitic lavas and pyroclastic deposits; Stage 3 (03^01Ma): eruption of basaltic
andesite to basalt lava flows and pyroclastic deposits
from two vents near the summit; Stage 4 (501Ma):
caldera-forming eruption (522 km in diameter), after
which a pyroclastic cone composed of basaltic andesites
built up within the caldera. The rocks of Stage 1 are
classified as low-K tholeiitic series, whereas those of
Stages 2^4 belong to the medium-K, calc-alkalic series.
This study is based on 39 volcanic samples collected at
14 sites from Stages 1 and 3 of the evolution of Zao
volcano (Fig. 3).
1975
JOURNAL OF PETROLOGY
ZA3011, ZA3012,
ZA3013, ZA3014
ZA3041, ZA3042
VOLUME 49
NUMBER 11
NOVEMBER 2008
ZA3021, ZA3022,
ZA3023,ZA3024, ZA3031, ZA3032
ZA3025
ZA1041,
ZA1021 ZA1042,
ZA1043,
ZA1044 ZA3081,
ZA3082
2 km
Stage 1 (<0.1Ma)
Stage 2 (0.3-0.1Ma)
Stage 3 (~0.3Ma)
ZA3051
Stage 4 (1.0-0.6Ma)
ZA3071,
ZA3061, ZA3062
ZA3072
ZA3063, ZA3064
ZA1031, ZA1032, ZA1033
ZA3091, ZA3101,ZA3102-1 ZA1011,
ZA3093 ZA3102-2, ZA3103 ZA1012, ZA1013
Basement &
other volcanics
Crater &
caldera
Fig. 3. Simplified geological map of Zao volcano after Sakayori (1992), showing sampling localities.
A N A LY T I C A L M E T H O D S
Bulk-rock analyses
Rock samples for trace element and Sr^Nd^Pb isotope
analyses were crushed to coarse chips (505 mm3) and
fresh pieces were hand picked. To avoid surface contamination, the rock chips were washed with ethanol and then
leached with 10M HCl at room temperature for 1h.
These chips were rinsed three times with Milli-Q water
an then dried in an oven. The chips were pre-pulverized
using a tungsten carbide mortar. At this step, any remaining altered portions were removed. Finally, the coarse
grain samples were pulverized in a tungsten carbide vibration mill.
Major and trace element (Ni, Cu, Zn, Rb, Sr, Y, Zr, Nb,
Ba, Pb, and Th) compositions were measured on fused
glass beads and pressed powder pellets, respectively, using
RIGAKUÕ Simultix 12 and RIX3000 X-ray fluorescence
(XRF) spectrometers. Analytical procedures, precision
(mostly 51%), and accuracy (mostly 51%) have been
fully described by Tani et al. (2006).
Concentrations of rare earth and 17 other trace elements
(Sc, Co, Ni, Cu, Rb, Sr, Y, Zr, Nb, Cs, Ba, Hf, Ta, Tl, Pb,
Th, and U) were analyzed using by inductively coupled
plasma mass spectrometry (ICP-MS) using an Agilent
7500ce system fitted with PFA sample introducing and a
Pt-inject torch system. The ICP-MS system was operated
in no collision gas and multi-tune acquisition mode.
This combination allowed a wide range of elements to be
precisely determined using pulse counting detection with
the hydrofluoric acid containing sample solution being
delivered directly into the plasma. Sample dissolution,
preparation and measurement were described by Chang
et al. (2003) and Nakamura & Chang (2007). Analytical
accuracy and precision for ICP-MS analyses, estimated
from repeated measurements of international reference
rocks (JB-1a of GSJ, BCR-2 and BIR-1 of USGS) were
mostly better than 5% and 3%, respectively.
The analytical procedure used for chemical separation
and mass spectrometry for Sr, Nd and Pb isotope determinations was outlined by Miyazaki et al. (2007). Total procedural blanks for Sr, Nd and Pb were less than 15 pg, 3 pg
and 6 pg, respectively. Mass spectrometry was performed
on a Thermo-FinniganÕ Triton TI equipped with nine
Faraday cups, using a static multi-collection mode. Normalizing factors used to correct for isotopic fractionation
in Sr, Nd and Pb isotope analyses were 86Sr/88Sr ¼ 01194,
146
Nd/144Nd ¼ 07219, and 0148% per atomic mass unit,
respectively. Measured isotopic ratios for standard materials were 87Sr/86Sr ¼ 0710262 14 (2s) for NIST 987
(n ¼ 40), 143Nd/144Nd ¼ 0512098 13 (2s), for JNdi-1
(n ¼ 31), and 208Pb/204Pb ¼ 36717 7 (2s), 207Pb/204Pb ¼
15497 2 (2s) and 206Pb/204Pb ¼16940 2 (2s) for
NIST 981 (n ¼ 23).
Micro-analyses
Mineral compositions of major and minor elements were
analyzed by electron-probe micro-analysis (EPMA) using
a JEOL JXA-8800 instrument, following the method
described by Shukuno (2003). For all elements the excitation potential, specimen current, and analytical time of
peak and background were: on olivine 20 kV, 25 nA, 20 s,
and 10 s (except for Mn, Ca, and Ni, for which 25 kV,
20 nA, 100 s and 50 s were used); on spinel 15 kV, 12 nA,
1976
TATSUMI et al.
CALC-ALKALIC VS THOLEIITIC SERIES
and 20 s; on pyroxene and plagioclase 15 kV, 15 nA,
and 20 s. ZAF correction procedures were employed.
Trace element [Rb, Sr, Y, Zr, Nb, Ba, rare earth elements
(REE), Hf, Ta, Pb, Th, and U] micro-analyses of plagioclase were performed by laser-ablation (LA)-ICP-MS on
the same thick sections that were used for the Sr isotope
analyses (see below). The 193 nm excimer laser system
aerosol line was connected to a VG ElementalÕ PQ3
quadrupole-type ICP-MS system (Kimura et al., 2000).
We used the 43Ca peak for internal standardization,
which corrected any variation in the ablated sample
volume. CaO contents were determined by EPMA prior
to the LA-ICP-MS analysis. NIST 612 synthetic glass was
used as a standard, adopting the reference values of
Pearce et al. (1996). The element concentrations measured
in NIST 612 were within 85% of the reference values.
The overall homogeneity of the glass for most elements
was better than 10% (2 SD), when 20 mm laser spots were
used (Kimura et al., 2000). The LA conditions were set at
200 mJ laser source energy, 50 mm crater size, and repetition rate at 5 Hz, yielding signals of 2000 c.p.s./ppm at
115
In. Ar gas blank was measured for 60 s for blank subtraction prior to sample analysis. Three spots on the NIST
612 standard were analyzed and averaged to minimize
errors caused by heterogeneity. Typical analytical time for
a single crater was 60 s, generating craters about 30 mm
deep. Accuracy of the results can be affected by the
matrix effect between the synthetic glass standard and the
silicate minerals. However, accuracy is generally better
than 15% when the aerosol loading to the plasma is controlled by 50% (Kimura et al., 2000). Such accuracy is
sufficient for most geochemical and petrological purposes.
Sr isotopes were analyzed in the plagioclases using two
techniques; a combined micromilling followed by chemical
separation and thermal ionization mass spectrometry
(MM-TIMS), and in situ laser ablation multi-collector
inductively coupled plasma source mass spectrometry
(LA-MC-ICP-MS).
The MM-TIMS techniques we used were similar to
those of Davidson & Tepley (1997), Tepley et al. (1999) and
Charlier et al. (2006), and were fully described byTakahashi
et al. (2005). The rock sample was cut into a wafer with
a thickness of 10 mm, which was then bonded onto a
glass slide and polished. Micromilling of the sample was
carried out using a New WaveTM MicroMillTM. The diameter at the tip of the drill used for sampling is 027 mm.
After micromilling, the collected sample powder was
decomposed with HF, HCl and HNO3. Sr selective extraction resin (Sr Resin from EICHROM Technologies Inc.)
was used for chemical separation of Sr (Horwitz et al.,
1992). Resin was charged into a modified pipette tip
column with quartz wool filter. Bedded resin volume was
005 ml. Total procedural blanks for Sr in the MM-TIMS
procedure were less than 10 pg. Methods used for Sr
isotope MM-TIMS analyses are identical to those used
for the bulk-rock analyses. Repeated analyses of NIST
987 (10 ng Sr) gave 87Sr/86Sr ¼ 0710260 22 (2s, n ¼ 8).
Additionally, NIST 610 glass was used for estimating the
analytical precision. Analysis of NIST 610 by the bulkrock analysis procedure and the MM-TIMS techniques
gave 87Sr/86Sr ¼ 0709681 07 (2s, n ¼ 4) and 87Sr/86Sr ¼
0709696 29 (2s, n ¼ 5), respectively. These values are
almost identical to the value (87Sr/86Sr ¼ 0709699 18)
reported for NIST 610 (Woodhead & Hergt, 2001).
In situ Sr isotope micro-analyses were also performed
by LA-MC-ICP-MS using a VG ElementalÕ Plasma 54
MC-ICP-MS system equipped with a dry and solution
aerosol dual sample intake system. The dual intake system
consists of an in-house 193 nm ArF excimer laser ablation
(Kimura et al., 2000) with He carrier gas and an AridusÕ
desolvating nebulizer using Ar (þ trace N2) carrier gas.
The two carrier gas lines were mixed in a 50 cm3
volume TeflonÕ mixing chamber prior to the ICP torch.
While the laser aerosol alone was analyzed, the Milli-QÕ
deionized water was taken up by the AridusÕ solution line,
for which blanks were negligible. Standard solutions were
introduced into the ICP torch while LA was unfired with
the LA carrier gas on, for determination of instrumental
mass bias and isobaric overlap correction factors (see
below). We hereafter call this the ‘dual intake system’.
An on-peak background method was applied to subtract
the blank from the Kr impurity (82Kr, 83Kr, 84Kr, and
86
Kr) in the Ar plasma gas (Woodhead et al., 2005). This
is necessary rather than peak stripping using mass bias
factor by monitoring one of the Kr peaks, because mass
bias for Kr cannot be determined because of complex
interferences during the sample analyses (Vroon et al.
2008). Background signals were acquired for the first 15 s
on the gas blank and the signal from the sample or standard was then collected, typically for 5 min, depending on
the sample.
The 87Rb interference on 87Sr was corrected for by
monitoring 85Rb using an empirical overlap correction
factor, determined from a Rb-doped NIST 987 Sr isotope
standard solution (20 ppb Rb in 100 ppb Sr, Rb/Sr ¼ 02)
from the AridusÕ aerosol line. The mass bias factor of Rb
is not identical to that of Sr (Hirata, 1996) and Rb has only
two stable isotopes, which prevents internal mass bias
correction. The mass bias factor for Rb may be different
between solution and LA modes (Vroon et al., 2008).
However, our simultaneous dual intake system cancels
out the mass bias difference caused by the different ICP
operating parameters in discrete LA aerosol or solution
aerosol introduction. Therefore, our empirical Rb overlap
correction method is an alternative to correction factor
determination using a synthetic glass standard (Davidson
et al., 2001). The difference between LA and AridusÕ solution analyses with our dual intake system is in the presence
1977
JOURNAL OF PETROLOGY
VOLUME 49
or absence of matrix elements for the LA and the AridusÕ
aerosols, respectively. The non-spectral matrix effect
(Barling & Weis, 2008) should be present and is noted
below. The rate of 87Rb overlap on 87M during plagioclase
analyses was typically less than 1% (Rb/Sr 5 000025)
and occasionally exceeded 10% (Rb/Sr 4 0025) with
plagioclase containing glass inclusions. However, we did
not see any problems caused by the change in the Rb overlap correction factor on the 87Sr/86Sr ratios of plagioclase
within analytical precision (87Sr/86Sr ¼ 000005).
The CaAr and Ca dimmer molecular ion interferences
on masses 82M, 84M, 86M, and 88M were subtracted by
monitoring 82M using the correction method proposed
by Woodhead et al. (2005). The effect of the correction
was monitored by the most interference-sensitive isotope
ratio, 84Sr/86Sr, the value of which was typically 00565
(Woodhead et al., 2005). We have confirmed the validity of
the correction method by analyzing a Ca-doped NIST 987
standard solution. The reference isotope ratio of NIST 987
was reproducible within the analytical precision with
a standard solution containing 50 ppm Ca and 100 ppb Sr
(Ca/Sr ¼ 500), in which Ca/Sr was more than that in the
natural plagioclase (50^200) and comparable with
natural carbonates (500).
Other interferences from doubly charged ions such
as Yb2þ are not considered. Plagioclase crystals have low
element abundances of REE (see trace element compositions of the plagioclase crystals below). The effect of the
doubly charged heavy ions is negligible (Vroon et al., 2008).
Mass bias, which was not considered by internal correction using 86Sr/88Sr ¼ 01194, caused by day-to-day basis
changes in interface cone or operating parameters (gas
flows, sampling depths, etc., Woodhead et al., 2001), was
further corrected for using the beta factor (Pachett et al.,
1981; Iizuka & Hirata, 2005) determined by analyses of a
solution of NIST 987 with 87Sr/86Sr ¼ 070125 from the
AridusÕ line. The dual sample intake system was again
advantageous for performing this complex correction
procedure, as it allows immediate switching between the
two introduction lines without changing the condition of
the plasma. Absence and presence of matrix elements may
cause change in Sr mass bias between the LA and the standard solution aerosols. However, increase of the matrix element by the addition of 50 ppm Ca to the NIST 987
solution did not make any detectable change in 87Sr/86Sr,
confirming previous reports (Woodhead et al., 2005;
Vroon et al., 2008). The NIST 987 standard bracketing
method was used to perform the beta correction after a
series of 10 analyses on the plagioclase.
A source energy of 200 mJ was used for the laser, with
a repetition rate of 20 Hz and a diameter of 200 mm.
This yielded a stable signal of 2^3 V on 88Sr for plagioclase with 300 ppm Sr over 5 min by our LA system,
with crater penetration to about 400^500 mm in the
NUMBER 11
NOVEMBER 2008
thick section. All the plagioclase analyses were performed
in single spot mode rather than raster mode. The 87Sr/86Sr
ratios of individual plagioclase spots were measured for
5 s each. The downhole Sr isotopic zoning was carefully
rejected by observing the time resolving profiles. Integrated data from the homogeneous section of the profile
were used for further statistical treatment, including average and two standard error calculations. Downhole
element fractionation has been reported between U and
Pb in a zircon crystal (Horn et al., 2000). On the other
hand, such fractionation was not detectable during Hf
isotope analysis of zircon crystal including the overlap
correction factors of Lu and Yb on 176Hf (Woodhead
et al., 2004). The downhole elemental fractionation may
occur between Rb and Sr, likewise between Pb and U.
However, we did not see any detectable change in the Sr
isotopic ratio in homogeneous plagioclase crystals with
sporadic melt inclusions at different depth levels. In fact,
overall precision in the single spot analyses of the homogeneous crystals was better than 000005 (2 SE), irrespective of the presence or absence of melt inclusions that
contain Rb. This indicates that isotopic fractionation
during downhole laser ablation is unlikely or at least not
detectable with our analytical precision, for both Sr and
Rb, which is similar to the case of Pb, Hf, Yb, and Lu isotopes in zircon (Horn et al., 2000; Woodhead et al., 2004).
The accuracy of the LA-ICP-MS plagioclase method
was tested by comparing MM-TIMS data from the
same homogeneous plagioclase (J.-I. Kimura, unpublished
data, and the present study); the results agreed within
00001, which is almost the same level as those reported
elsewhere for LA-MC-ICP-MS Sr isotope analyses (Vroon
et al., 2008). This accuracy is acceptable for the purpose
of this study.
R E S U LT S
Bulk-rock chemistry
Major and trace element concentrations and Sr^Nd^Pb
isotopic compositions of Zao volcanic rocks are reported
in Table 1.
The volcanic rocks define two distinct chemical trends;
tholeiitic and calc-alkalic (Fig. 4). These are generally
regarded as equivalent to Kuno’s pigeonitic and hypersthenic rock series, and Sakuyama’s N-type and R-type rocks,
respectively. Tholeiitic magmatism occurred solely during
Stage 1, whereas the rocks of Stages 2^4 are classified as
calc-alkalic. It should be stressed that the Zao volcanic
rocks form chemical trends similar to those that typify the
tholeiitic and calc-alkalic rocks of the trench-side Nasu volcanic chain (Fig. 4), suggesting that the genetic relation
between the two magma series at Zao may be applied to
the NE Japan arc in general. Among the major elements,
the calc-alkalic rocks tend to be more depleted in Fe and
more enriched in Mg and K than the tholeiitic rocks.
1978
TATSUMI et al.
CALC-ALKALIC VS THOLEIITIC SERIES
Table 1: Major and trace element, isotopic and modal compositions of Zao volcanic rocks
Rock type:
Tholeiitic
Sample:
ZA1013
ZA1011
ZA1033
ZA1032
ZA1021
ZA1042
ZA1012
ZA1043
ZA1044
ZA1031
ZA1041
5531
SiO2 (wt %)
5198
5202
5215
5241
5261
5274
5275
5416
5472
5512
TiO2
080
078
095
093
081
080
078
083
100
094
104
Al2O3
1798
1774
1814
1792
1786
1781
1790
1802
1718
1705
1719
Fe2O3
1095
1109
1066
1111
1111
1026
1034
1087
962
1066
1010
MnO
019
018
019
019
020
018
018
017
019
018
019
MgO
549
597
475
494
563
582
547
465
408
411
394
779
CaO
973
1028
946
946
985
961
1024
924
862
891
Na2O
197
204
239
239
217
228
212
241
261
267
256
K2O
028
026
026
032
037
028
025
044
041
048
043
P2O5
009
011
016
016
016
014
009
016
020
018
011
Total
9960
10004
9956
9983
9992
10000
10065
9970
9967
9974
9951
Ni (ppm)
20
22
15
15
25
31
17
19
8
9
8
Cu
45
38
36
39
41
41
44
39
37
35
31
Zn
87
89
99
96
89
87
79
89
103
97
102
Rb
64
57
50
62
70
41
50
100
90
11
97
Sr
248
249
274
274
263
258
255
269
263
271
257
Y
26
18
18
20
21
20
16
19
20
20
28
Zr
44
46
54
56
55
55
43
65
69
72
73
Nb
Ba
Pb
11
105
40
Th
17
104
21
97
20
120
34
19
49
15
07
06
24
135
23
134
32
32
16
129
40
24
156
29
161
32
167
43
46
58
20
13
21
30
163
42
08
Sc (ppm)
362
375
335
329
319
345
342
Co
396
461
377
388
407
335
343
Ni
208
250
177
184
288
Cu
405
364
337
355
380
Rb
Sr
610
247
550
248
460
272
578
274
857
344
672
859
263
258
853
289
929
253
Y
264
187
186
207
221
199
278
Zr
459
482
565
578
563
714
745
Nb
147
166
225
229
225
309
Cs
0429
0361
0255
0354
0344
0348
Ba
La
Ce
Pr
Nd
102
663
141
229
112
987
447
111
106
482
120
117
563
136
127
160
649
638
135
158
156
173
196
205
214
760
85
942
972
993
319
0284
167
104
193
286
132
Sm
331
226
251
269
268
279
345
Eu
111
0825
0968
102
0982
104
126
Gd
433
295
313
341
346
335
439
Tb
074
051
0536
0581
0576
0573
0719
Dy
472
335
345
375
369
367
456
Ho
100
0721
0737
081
0794
0779
0988
Er
295
217
220
240
236
234
293
Tm
0416
031
0317
0342
0335
0336
041
Yb
270
206
212
229
220
229
272
(continued)
1979
JOURNAL OF PETROLOGY
VOLUME 49
NUMBER 11
NOVEMBER 2008
Table 1: Continued
Rock type:
Tholeiitic
Sample:
ZA1013
ZA1011
ZA1033
ZA1032
ZA1021
ZA1042
ZA1012
ZA1043
ZA1044
ZA1031
ZA1041
Lu
0414
0321
0329
0352
0342
0348
0425
Hf
144
146
167
173
166
207
215
Ta
0081
0086
0106
0109
011
0154
0155
Tl
0063
0048
0031
0082
0049
0046
0124
Pb
353
320
301
326
329
425
385
U
0649
0555
0537
0598
0664
0943
0973
Th
0171
015
0139
0158
0171
0245
87
0704189
0704227
0704361
0704351
0704332
Sr/86Sr
2s
143
Nd/144Nd
2s
9
8
0512797
9
10
0512793
10
7
0512760
9
9
0512770
10
0704327
9
0512777
10
0704176
8
0512779
9
0704236
10
0512818
10
0219
0704465
9
0512784
0704439
8
0512724
10
10
0704465
9
0512758
10
0512768
10
206
Pb/204Pb
18472
18472
18454
18452
18453
18462
18473
18451
18459
18465
18465
207
Pb/204Pb
15582
15582
15582
15575
15577
15588
15583
15576
15580
15586
15582
208
Pb/204Pb
38511
38506
38503
38482
38490
38526
38513
38486
38509
38528
10
17
10
09
07
08
01
olivine
02
—
—
38515
—
clinopyroxene
15
39
22
21
43
11
32
15
14
24
orthopyroxene
59
59
09
15
30
39
24
13
62
27
05
plagioclase
281
313
264
237
265
223
279
191
313
195
125
quartz
—
—
—
—
—
—
—
—
opaque
—
—
groundmass
643
589
Rock type:
Calc-alkalic
Sample:
ZA3093
ZA3091
SiO2 (wt %)
01
694
ZA3032
—
—
—
—
—
—
710
662
717
ZA3031
ZA3064
ZA3102-1
01
05
655
769
ZA3103
ZA3101
ZA3011
—
603
ZA3102-2
ZA3012
33
02
02
751
835
ZA3013
ZA3014
5245
5258
5562
5590
5643
5677
5679
5688
5689
5699
5709
5718
TiO2
085
086
082
082
085
084
079
083
081
084
081
082
5739
082
Al2O3
1739
1738
1705
1706
1698
1726
1740
1728
1663
1729
1702
1670
1670
Fe2O3
982
980
872
865
879
862
830
865
815
861
810
846
844
MnO
016
016
015
015
016
015
014
014
014
015
013
015
015
MgO
698
673
500
486
505
439
425
435
414
443
400
413
411
CaO
888
874
833
832
770
706
732
687
736
705
769
695
700
Na2O
247
251
253
260
271
258
262
259
259
258
270
253
253
K 2O
088
091
112
115
114
124
122
126
124
124
117
121
121
P 2O5
018
018
013
014
016
015
015
016
013
015
012
013
013
Total
10008
9985
9946
9965
9996
9905
9897
9901
9807
9932
9883
9824
9849
Ni (ppm)
58
27
31
21
21
21
Cu
36
42
25
38
38
36
Zn
70
68
66
77
77
72
Rb
27
36
32
37
37
36
Sr
391
317
340
306
306
277
Y
20
23
23
26
26
23
Zr
81
101
100
103
103
108
Nb
Ba
22
323
26
325
30
389
30
363
30
363
24
350
(continued)
1980
TATSUMI et al.
CALC-ALKALIC VS THOLEIITIC SERIES
Table 1: Continued
Rock type:
Calc-alkalic
Sample:
ZA3093
ZA3091
ZA3032
ZA3031
ZA3064
ZA3102-1
ZA3103
ZA3101
ZA3011
ZA3102-2
ZA3012
ZA3013
Pb
49
59
60
74
74
Th
26
44
32
40
40
Sc (ppm)
271
273
274
282
257
Co
422
454
343
302
410
Ni
682
322
368
224
240
Cu
352
390
241
344
322
Rb
266
349
309
362
Sr
396
Y
214
Zr
851
312
229
106
333
239
104
258
263
320
321
186
113
105
La
329
957
Ce
950
224
Pr
225
296
Nd
295
133
132
374
120
251
261
328
232
110
276
0884
378
109
145
338
268
273
0519
315
34
105
Nb
Ba
175
299
Cs
339
960
228
370
299
164
132
Sm
342
339
367
41
344
Eu
107
0951
108
115
0974
Gd
386
389
412
466
387
Tb
0631
0659
0687
0763
0658
Dy
391
414
429
472
417
Ho
0817
0875
0901
0995
0884
Er
242
261
270
296
267
Tm
0345
038
039
0424
039
Yb
230
254
261
283
262
Lu
0354
0386
0405
0435
0402
Hf
249
320
310
313
339
Ta
0161
0176
0211
0206
0193
Tl
0081
0113
0168
0243
0205
Pb
385
595
540
766
U
258
348
340
338
Th
0641
0863
0834
0831
091
87
0703596
0703786
0703733
0703815
0703850
Sr/86Sr
2s
8
143
Nd/144Nd
9
0512882
2s
0512864
10
Pb/204Pb
10
8
0512853
10
164
370
9
10
0512840
0512851
10
15
206
18389
18415
18410
18416
18420
207
15551
15558
15555
15561
15565
208
38357
38398
38394
38414
Pb/204Pb
Pb/204Pb
olivine
47
23
03
01
03
clinopyroxene
45
69
85
69
76
orthopyroxene
31
40
23
16
plagioclase
196
290
232
quartz
—
—
—
opaque
groundmass
—
ZA3014
38424
03
03
—
63
90
81
63
63
70
71
20
24
35
21
21
24
23
24
22
297
227
221
230
209
303
221
261
270
234
—
—
—
—
—
—
—
06
03
—
—
58
06
04
01
06
02
11
03
07
08
07
09
16
08
19
13
07
675
576
646
614
667
678
631
676
597
678
624
622
679
(continued)
1981
JOURNAL OF PETROLOGY
VOLUME 49
NUMBER 11
NOVEMBER 2008
Table 1: Continued
Rock type:
Calc-alkalic
Sample:
ZA3072
ZA3041
SiO2 (wt %)
ZA3071 ZA3042 ZA3051
ZA3061 ZA3062 ZA3025 ZA3063
ZA3021 ZA3081
ZA3024 ZA3023 ZA3022
ZA3082
6273
5741
5750
5752
5755
5861
5922
5966
5998
6060
6133
6161
6174
6187
6196
TiO2
081
082
082
083
081
077
076
079
076
073
071
073
072
074
070
Al2O3
1698
1656
1692
1668
1633
1680
1708
1605
1626
1602
1609
1560
1559
1546
1563
Fe2O3
839
879
849
879
868
794
765
849
757
748
720
748
753
747
698
MnO
015
015
015
015
014
014
014
014
013
013
013
013
013
013
013
MgO
469
437
475
434
413
393
365
364
359
326
334
321
309
316
319
CaO
707
753
702
751
713
634
609
665
610
617
609
606
570
575
571
Na2O
276
267
274
267
268
277
287
271
289
282
300
285
284
285
305
K 2O
123
124
124
124
129
120
113
149
145
156
151
173
175
179
163
P 2O5
015
013
016
013
012
013
013
013
013
011
012
011
012
012
012
Total
9963
9974
9981
9987
9991
9925
9915
10006
9946
9962
9980
9964
9932
9942
9987
Ni (ppm)
30
22
19
21
19
14
Cu
30
21
21
15
36
27
Zn
66
65
61
60
58
57
Rb
38
37
38
44
46
58
Sr
330
278
252
273
271
236
Y
23
24
24
29
26
29
Zr
102
110
116
122
124
146
Nb
Ba
28
372
24
338
25
31
350
418
32
415
30
485
Pb
63
84
97
75
67
Th
33
46
48
48
47
66
Sc (ppm)
269
259
260
223
217
243
Co
308
391
336
282
282
374
Ni
349
255
234
235
209
173
Cu
284
213
204
144
329
275
Rb
361
359
371
423
442
Sr
Y
Zr
318
243
104
273
240
114
246
262
247
293
120
120
262
265
117
68
558
231
288
122
Nb
316
276
290
358
365
375
Cs
106
0846
0916
125
121
155
Ba
382
La
109
Ce
244
Pr
Nd
321
141
341
960
224
293
130
346
424
970
225
294
447
125
125
279
270
283
396
130
427
134
172
359
154
368
161
Sm
360
343
344
422
382
406
Eu
105
0942
0914
111
101
098
Gd
407
396
402
484
433
463
Tb
0679
067
0681
0801
0729
0784
Dy
426
426
434
499
454
492
Ho
090
0899
0925
106
0968
105
Er
270
270
277
313
290
314
(continued)
1982
TATSUMI et al.
CALC-ALKALIC VS THOLEIITIC SERIES
Table 1: Continued
Rock type:
Calc-alkalic
Sample:
ZA3072
ZA3041
ZA3071 ZA3042 ZA3051
ZA3061 ZA3062 ZA3025 ZA3063
ZA3021 ZA3081
ZA3024 ZA3023 ZA3022
Tm
0393
0393
0405
0457
0421
0457
Yb
262
264
271
301
283
303
Lu
0405
0407
0419
0466
0441
0465
Hf
311
348
362
364
359
377
Ta
0211
0184
0196
0242
0244
0259
Tl
0149
0155
0155
0218
0177
0125
Pb
597
691
860
631
674
645
U
350
375
391
434
442
536
Th
0847
0926
0958
105
103
128
87
0703792
0703938
0704101
0703935
0703733
0703982
Sr/86Sr
2s
143
Nd/144Nd
2s
Pb/204Pb
9
0512838
12
10
9
0512847
7
0512837
11
8
0512820
10
8
0512798
9
0512834
10
11
206
18436
18409
18424
18437
18423
18422
207
15585
15555
15564
15580
15561
15563
208
38498
38388
38424
38485
38422
Pb/204Pb
Pb/204Pb
olivine
09
—
02
clinopyroxene
55
101
67
orthopyroxene
—
76
00
71
—
73
01
—
59
64
00
65
—
56
03
38
ZA3082
38422
—
53
—
83
—
71
02
61
25
10
34
25
21
17
15
29
11
27
29
24
08
16
12
plagioclase
209
253
215
237
263
245
226
239
236
307
250
303
297
290
288
quartz
—
—
—
—
21
opaque
groundmass
03
05
02
00
06
04
17
16
11
07
09
09
01
11
10
16
04
14
17
19
10
09
24
13
14
695
627
672
651
635
646
690
652
671
585
666
594
572
599
602
Total iron as Fe2O3.
Ni–Th by XRF; Sc–Th by ICP-MS.
As a result they belong to the medium-K series, whereas
the tholeiitic rocks belong to the low-K series (Fig. 4).
The normal mid-ocean ridge basalt (N-MORB) normalized incompatible trace element patterns of the Zao
volcanic rocks are shown in Fig. 5. Although elements
with higher incompatibility during mantle melting tend to
be more enriched, the high field strength elements such as
Nb, Ta and Zr do not show such enrichment. This selective
enrichment of particular elements results in strongly spiked
patterns (Fig. 5), which have also been observed in other
arc lavas (e.g. Tatsumi & Eggins, 1995). Calc-alkalic rocks
are characterized by higher concentrations of incompatible
trace elements than tholeiitic rocks.
The Sr^Nd^Pb isotope compositions of the Zao volcanic
rocks are also well within the range of the trench-side Nasu
lavas (Fig. 6). It has been well established that volcanic rocks
along the NE Japan volcanic front exhibit more enriched
Sr^Nd isotopic signatures towards the south (Notsu, 1983;
Kimura & Yoshida, 2006). Zao volcano is situated in the
southernmost part of Zone I of Kimura & Yoshida (Fig. 2)
and has Sr^Nd isotopic characteristics typical of this zone
(Fig. 6). The Pb isotopic ratios of Zao, and more generally
of Quaternary volcanic rocks from the NE Japan arc, form
a broad trend towards the compositions of Pacific sediments
from the MORB field (Fig. 6).
Calc-alkalic and tholeiitic rocks from Zao volcano exhibit systematic differences in both major/trace element and
isotopic characteristics. To emphasize the difference in
element concentrations, element abundances are compared
in Fig. 7 by normalizing the calc-alkalic compositions to
the tholeiitic composition at 55 wt % SiO2. It should be
stressed that the calc-alkalic rocks are more enriched in
both highly incompatible elements (U, Rb, Th, and K)
and compatible elements (Ni and Mg). Isotopically, the
calc-alkalic rocks exhibit more depleted characteristics in
terms of Sr, Nd and Pb isotopes (Fig. 6).
Petrography
Modal proportions of phenocrysts and representative mineral compositions are given in Tables 1^7.
1983
JOURNAL OF PETROLOGY
VOLUME 49
NUMBER 11
NOVEMBER 2008
SiO2 (wt.%)
5
45
50
55
60
SiO2 (wt.%)
65
70
45
12
Zao Tholeiitic
55
60
65
70
11
Zao Calc-alkalic
iitic lic
a
ole
Th c-alk
l
a
C
3
10
FeO* (wt.%)
FeO*/MgO
4
50
2
9
8
7
Nasu Chain TH
Nasu Chain CA
6
0
5
1.4
10
1.2
8
MgO (wt.%)
TiO2 (wt.%)
1
1.0
0.8
6
4
0.6
2
0.4
0
22
2.5
2.0
K2O (wt.%)
Al2O3 (wt.%)
20
18
16
14
45
1.5
-K
dium
1.0
Me
K
Low-
0.5
0
50
55
60
65
70
45
SiO2 (wt.%)
50
55
60
65
70
SiO2 (wt.%)
Fig. 4. Major element vs SiO2 variation diagrams for volcanic rocks of Zao and the Nasu trench-side volcanic chain. Zao volcano is composed
of two magma series, tholeiitic (TH) and calc-alkalic (CA), which form chemical trends that are typical of those of the magma series of the
trench-side Nasu volcanic chain. Calc-alkalic and tholeiitic rocks in Zao volcano belong to the medium- and low-K series, respectively.
Subsequent petrographic description of the Zao volcanic
rocks (below) highlights the differences between the calcalkalic and tholeiitic rocks.
Olivine phenocrysts (usually51vol. %) occur as euhedral
to subhedral crystals unrimmed or rarely mantled by orthopyroxene in basaltic andesites to andesites from both the
calc-alkalic and tholeiitic series. Although there is little difference in the Mg-number of the core compositions between
the two magma series (65^85), NiO contents are higher in
the calc-alkalic rocks at a constant Mg-number (Fig. 8).
Orthopyroxene phenocrysts are ubiquitous in the Zao
volcanic rocks. Those in rocks from the tholeiitic series
1984
TATSUMI et al.
CALC-ALKALIC VS THOLEIITIC SERIES
Plagioclase is the most abundant phenocryst phase and
generally makes up 10^30 vol. % of the analyzed samples.
The plagioclase phenocrysts in the tholeiitic rocks have a
limited compositional range, whereas in the calc-alkalic
rocks they exhibit a broader range of compositions with
a bimodal distribution (Fig. 10). Notably, the plagioclase
phenocrysts in the calc-alkalic rocks commonly have
honeycomb textures and dusty zones. Representative textures of plagioclase phenocrysts are shown in Fig. 11,
together with compositional profiles across the crystal for
Ca/(Ca þ Na) ratio and Sr concentration.
N-MORB Normalized
100
10
Calc-alkalic
Plagioclase microanalyses
1 Tholeiitic
Cs Rb Th Nb K Ce Pr Nd Hf Eu Tb Y Er Yb
Tl Ba U Ta La Pb Sr Zr Sm Gd Dy Ho Tm Lu
Incompatibility during mantle melting
Chondrite Normalized
100
50
Calc-alkalic
Tholeiitic
10
La Ce Pr Nd
Sm Eu Gd Tb Dy Ho Er Tm Yb Lu
Fig. 5. Incompatible trace element and REE characteristics of
the Zao volcanic rocks. (a) MORB-normalized; (b) chondritenormalized. Normalizing values are from Sun & McDonough (1989).
exhibit little or normal zoning with a higher Mg-number
core mantled by a lower Mg-number rim (Fig. 9). In
contrast, calc-alkalic rocks contain both normally and
reversely zoned orthopyroxene phenocrysts (Fig. 9), with
the cores of normally zoned orthopyroxene phenocrysts
having higher Mg-number (70^80) than those of the reversely zoned phenocrysts (60). Orthopyroxene phenocrysts
with reactions rims of clinopyroxene (pigeonite) occur
only in the tholeiitic rocks. Only in the calc-alkalic series
is orthopyroxene present as a groundmass phase.
Clinopyroxene (augite) phenocrysts occur in all Zao
volcanic rocks, although they are less abundant than
orthopyroxene. Augite phenocrysts have a rather limited
compositional range, both in tholeiitic and calc-alkalic
rocks, (Mg-number 70; Fig. 9). Reversely zoned clinopyroxene phenocrysts are found only in calc-alkalic rocks
(Fig. 9), whereas pigeonite exists only in the groundmass
of tholeiitic rocks.
Representative trace element concentrations and Sr isotopic compositions of plagioclase phenocrysts are listed in
Table 5. The correlations of selected trace element abundances and 87Sr/86Sr with Ca/(Ca þ Na) are shown in
Figs 12 and 13.
Plagioclase phenocrysts in calc-alkalic rocks tend to be
more enriched in Ba and Sr than those in tholeiitic rocks
(Fig. 12), consistent with the chemical characteristics of
the bulk-rocks (Fig. 7). Plagioclases from the two magma
series show little difference in REE concentrations,
whereas Yconcentrations in Ca-rich plagioclase from calcalkalic rocks are lower than those from tholeiitic rocks
(Fig. 12).
Sr isotopic ratios in plagioclases from tholeiitic rocks are
constant at 07042^07044 and show little correlation with
the anorthite content (Fig. 13), whereas those in the calcalkalic rocks show more complex characteristics. The core
compositions of plagioclase phenocrysts in calc-alkalic
rocks show a broad bimodal distribution in terms of
anorthite content (Fig. 10). Importantly, the high-An plagioclase from the most mafic calc-alkalic basaltic andesites
exhibit the lowest 87Sr/86Sr of 07034 (Fig. 13).
DISCUSSION
Closed vs open processes in the tholeiitic
series
Tholeiitic rocks from Zao volcano belong to the N-type
rocks of Sakuyama (1981) and Kuno’s pigeonitic rock
series, and exhibit little evidence for disequilibrium:
(1) phenocryst phases such as pyroxenes and plagioclase
are normally zoned (Figs 9 and 10); (2) olivine is in equilibrium with pyroxenes in terms of Fe^Mg partitioning
(Fig. 14) at temperatures of 10508C (Fig. 15). These observations, together with the relatively constant 87Sr/86Sr for
plagioclase phenocryst cores (Fig. 13), suggest that the Zao
tholeiitic magmas differentiate mainly in a closed system,
via fractional crystallization; this is a mechanism that has
been accepted as a general process of differentiation in arc
tholeiites (e.g. Sakuyama, 1981; Fujinawa, 1988, 1990;
Tatsumi & Kogiso, 2003; George et al., 2004; Villiger et al.,
2007). It should be stressed here that differentiation by
1985
JOURNAL OF PETROLOGY
VOLUME 49
NUMBER 11
NOVEMBER 2008
0.5129
15.60
38.5
0.5128
0.5134
15.55
39.5
38.3
18.4
18.5
18.4
15.7
18.5
39.0
15.6
Zao tholeiitic
0.5126
Zao calc-alkalic
15.5
NE Japan
38.5
Pb/204Pb
208
0.5130
Pb/204Pb
133
15.8
207
Nd/144Nd
0.5127
0.7035 0.7040 0.7045
38.0
Zone-I
0.5122
0.702
0.704
0.706
87
Sr/86Sr
0.708
Zone-II
17.8
18.2
Zone-III
206
18.6
204
Pb/
Pb
15.4
19.0
17.8
18.2
206
18.6
37.5
19.0
204
Pb/
Pb
Pacific MORB
Pacific sediments
Fig. 6. Sr^Nd^Pb isotope compositions of Zao volcanic rocks and Quaternary volcanoes in the trench-side Nasu chain. Data for NE Japan,
MORB and sediments are from Kimura & Yoshida (2006).
6
3
More enriched
in calc-alkalic
2
Comparable
More depleted in calc-alkalic
Calc-alkalic/Tholeiitic
4
1
the tholeiitic magmas. Assuming Abukuma granitic rocks
(Gr in Table 8, data from Kamei et al., 2003), which form
the basement of the Zao volcano, to be a possible upper
crustal contaminant, the observed variation in 87Sr/86Sr
for the tholeiitic rocks can be explained by simple bulk
contamination by upper crust (Fig. 16b).
Element Concentration
at 55% SiO2
5
Magma mixing in the calc-alkalic series
0.8
0.6
U Th K Pb Zr Pr Mg Sm Tb P Tm Na Al Eu Cu Ti Fe
Rb Ni Ba Ce Sr La Nd Yb Gd Lu Dy Ho Er Y Ca Nb Zn
Fig. 7. Geochemical differences between calc-alkalic and tholeiitic
rocks in Zao volcano. Element concentrations are normalized to the
tholeiitic composition at 55 wt % SiO2. Although for some of these
elements there is little difference between calc-alkalic and tholeiitic
rocks, calc-alkalic rocks are more enriched in both compatible
elements such as Ni and Mg and highly incompatible elements.
crystallization is not the only process that can occur
during closed-system differentiation. Partial melting, or
anatexis, for example, of basaltic lower crust could cause
differentiation from mafic to felsic compositions, which
will be discussed further below.
However, closed-system differentiation cannot solely
explain the geochemical characteristics of the Zao tholeiitic rocks. For example, plagioclase rims tend to show
higher 87Sr/86Sr (up to 07046) than plagioclase cores
(Fig. 13). Furthermore, a positive correlation between
SiO2 content and 87Sr/86Sr is observed for bulk-rock compositions (Fig. 16b). These observations suggest that an
open-system process, such as shallow-level crustal contamination, plays a role in controlling the final compositions of
‘Disequilibrium’ petrographic features in the calc-alkalic
rocks of Zao volcano indicate that they belong to the
R-type rocks of Sakuyama (1981). Such features include:
the occurrence of reversely zoned pyroxene phenocrysts
(Fig. 9); a wide compositional range in Ca/(Ca þ Na) of
plagioclase phenocrysts, which also show a bimodal compositional distribution (Fig. 10); the occurrence of plagioclase phenocrysts with honeycomb textures and dusty
zones; and the presence of disequilibrium olivine^pyroxene
pairs within a single specimen (Fig. 14).
Disequilibrium in the calc-alkalic rocks is further indicated by temperature estimates for pyroxene crystallization. Figure 15 summarizes the temperature estimates
from contiguous crystals of clinopyroxene and orthopyroxene using a two-pyroxene geothermometer (Wells,
1977). This thermometry cannot be applied to the pyroxenes in the groundmass of the tholeiitic rocks in Zao volcano, because they do not contain orthopyroxene; they
belong to Kuno’s pigeonitic rock series. Figure 15 indicates
that in calc-alkalic rocks the rim^rim pair and the groundmass pair of pyroxenes exhibit crystallization temperatures
higher than those obtained for the core^core pair.
One possible process that can provide a comprehensive
explanation of these disequilibrium textures would be
mixing of magmas having different compositions and temperatures (e.g. Eichelberger, 1975; Sakuyama, 1981, 1983;
1986
TATSUMI et al.
CALC-ALKALIC VS THOLEIITIC SERIES
Table 2: Representative compositions of olivine phenocrysts
Rock type:
Tholeiitic
Calc-alkalic
Sample:
ZA1011
Grain:
ol1-2
ol2-3
ol2-9
ol3-12
ol-12
ol-6
ol007
ol018
ol041
ol054
5ol-8
10ol-3
Position:
Core
Core
Core
Core
Core
Core
Core
Core
Core
Core
Core
Core
ZA1013
ZA1044
ZA3064
ZA3081
ZA3093
SiO2
4032
3888
3892
3882
3975
3786
4038
4034
3911
3786
3972
3864
FeO
1633
2115
2053
2645
1951
2915
1563
1589
197
3037
1803
2414
MnO
025
033
033
04
03
MgO
4504
4041
4085
3686
4241
CaO
012
013
013
017
013
017
012
012
013
01
015
011
NiO
009
003
004
004
004
004
017
014
011
003
010
005
10215
10093
10080
10274
10214
10206
10189
10179
10138
10234
10139
10190
Total
044
344
026
027
031
054
031
054
4533
4503
4202
3344
4308
3842
Si
0997
0996
0996
1000
0997
0997
0998
0999
0991
0999
0998
Fe
0338
0453
0439
0570
0409
0642
0323
0329
0417
0670
0379
0995
0520
Mn
0005
0007
0007
0009
0006
0010
0005
0006
0007
0012
0007
0012
Mg
1659
1542
1557
1415
1585
1349
1669
1661
1587
1315
1612
1473
Ca
0003
0003
0004
0005
0003
0005
0003
0003
0003
0003
0004
0003
Ni
0002
0001
0001
0001
0001
0001
0003
0003
0002
0001
0002
0001
Total
3004
3002
3004
3000
3001
3004
3001
3001
3007
3000
3002
3004
Mg/(Mg þ Fe)
0831
0773
0780
0713
0795
0678
0838
0835
0792
0662
0810
0739
Numbers of ions are calculated on the basis of four oxygens.
Wada, 1985; Fujinawa, 1988, 1990; Hunter & Blake, 1995;
Clynne, 1999; Streck et al., 2007). To better constrain the
magma mixing process, especially to decode the chemical
characteristics of the end-member components for the
mixed calc-alkalic magmas, the Sr isotopic compositions
recorded in plagioclase phenocrysts (Fig. 13) are considered in combination with the chemical composition of the
magma. Element concentrations for magmas that crystallize plagioclase in the Zao volcanic rocks are obtained
using element partitioning data (Bindeman et al., 1998;
Bindeman & Davis, 2000). The panels down the righthand side of Fig. 13 show the relationship between
87
Sr/86Sr and Sr concentration in the melt in equilibrium
with the plagioclase and allow the identification of
four melt components that mix to form the calc-alkalic
rocks. These are: a low 87Sr/86Sr (07034), Sr-rich
(650^700 ppm) melt (I in Fig. 13), probably having a
mafic/basaltic composition; and three components
(II, III, and IV in Fig. 13) with similar rather high
87
Sr/86Sr (07040), but with different Sr concentrations
(500^400, 300, 200^150 ppm, respectively), that are probably intermediate in composition. The above four melt
components may crystallize plagioclase with An490,
An80^90, An70, and An50^60, respectively (Figs 13 and 16a).
The next problem to address is how these four components are generated. Figure 16a summarizes the
geochemical characteristics of the four calc-alkalic melt
components (I^IV in Fig. 13) and the tholeiitic melts in
Zao volcano inferred from the core compositions of plagioclase phenocrysts. To understand the characteristics of
the Zao calc-alkalic melts, at least three end-member components may be required, L1, L2, and L3 in Fig. 16a and
Table 8. L1 and L2 are a calc-alkalic and a tholeiitic
basaltic melt, believed to be the host melts for Ca-rich
(An490) plagioclase found in the calc-alkalic and tholeiitic
basaltic andesites at Zao volcano, respectively. The SiO2
contents of these basaltic melts are assumed to be
50 wt %; Sr concentrations and 87Sr/86Sr are deduced
from plagioclase core compositions (Table 8). L3 is a tholeiitic felsic melt that is a liquid differentiated from L2 and
is assumed to contain 70 wt % SiO2. The 87Sr/86Sr of L3
is simply assumed to be identical to that of L2; the Sr concentration in L2 is 5200 ppm, and has been assumed to
be 100 ppm. The reason for this is that the identifiable
and the most differentiated tholeiitic melt in Zao volcano
contains 200 ppm Sr and crystallizes An60 plagioclase
(Fig. 13), which is more calcic than the plagioclase present
in the calc-alkalic intermediate component (An50^60;
Fig. 13). This suggests that a more differentiated tholeiitic
felsic melt with a Sr concentration lower than 200 ppm
is required as the felsic end-member for the calcalkalic rocks.
1987
JOURNAL OF PETROLOGY
VOLUME 49
NUMBER 11
NOVEMBER 2008
Table 3: Representative compositions of clinopyroxene
Rock type:
Tholeiitic
Calc-alkalic
Sample:
ZA1011
Grain:
px021 px021 pxg06 px042 px042 pxg58 px038 px038 pxg09 px014 px014 px024 px024 pxg030 px101 px101 px002 px002 pxg02
Position:
Core
ZA1013
Rim
GM
Core
ZA1044
Rim
GM
Core
ZA3064
Rim
GM
Core
ZA3093
Rim
Core
Rim
GM
Core
Rim
Core
Rim
GM
SiO2
5270
5294
5328
5285
5290
5268
5229
5261
5191
5313
5102
5287
5322
5403
5263
5277
5308
TiO2
042
034
036
033
040
031
042
035
044
038
063
035
068
044
027
031
067
057
070
Al2O3
275
207
157
174
186
140
195
212
168
251
282
097
349
102
088
090
182
211
166
Cr2O3
025
008
001
007
004
000
008
005
004
006
004
000
007
000
000
000
000
004
000
FeO
821
1038
1964
1176
1150 1929
992
1148
2111
777
881
1170
1028
2170
1067
2020
1252
1010
1056
5276 5240
MnO
023
026
045
033
047
015
021
032
020
025
044
030
066
037
057
035
031
039
MgO
1615
1649
1895
1548
1515 2111
1574
1570
1668
1618
1501
1380
1568
1830
1444
2311
1448
1473
1594
CaO
2050
1811
657
1834
1865
463
1943
1792
827
2091
2050
2080
1914
554
2052
193
1841
2055
1853
Na2O
021
023
007
023
028
008
022
021
017
024
026
025
024
014
020
007
023
024
032
Total
030
10142 10090 10090 10113 10094 9969 10081 10072 10100 10086 10023 10144 10090 10067 10057 10112 10111 10142 10118
Si
1922
1946
1973
1953
1952 1956
1949
1948
1960
1928
1923
1971
1887
1979
1979
1976
1952
1941
1951
Ti
0118
0090
0068
0076
0081 0062
0085
0092
0074
0108
0123
0043
0152
0045
0039
0039
0080
0092
0072
Al
0011
0009
0010
0009
0011 0009
0012
0010
0012
0010
0017
0010
0019
0012
0008
0009
0019
0016
0019
Cr
0007
0002
0000
0002
0001 0000
0002
0001
0001
0002
0001
0000
0002
0000
0000
0000
0000
0001
0000
Fe
0250
0319
0608
0363
0356 0602
0306
0355
0662
0238
0273
0363
0318
0679
0332
0618
0388
0311
0325
Mn
0007
0008
0014
0010
0009 0015
0005
0007
0010
0006
0008
0014
0009
0021
0012
0018
0011
0010
0012
Mg
0878
0903
1046
0852
0835 1174
0864
0865
0931
0884
0829
0763
0864
1021
0800
1259
0800
0807
0873
Ca
0801
0713
0261
0726
0739 0185
0767
0710
0332
0821
0814
0826
0758
0222
0818
0076
0731
0810
0730
Na
0015
0016
0005
0016
0020 0005
0015
0015
0012
0017
0019
0018
0017
0010
0014
0005
0017
0017
0023
4005
Total
4009
4006
3985
4007
4004 4008
4005
4003
3994
4014
4007
4008
4026
3989
4002
4000
3998
4005
Mg/(Mg þ Fe)
0778
0739
0632
0701
0701 0661
0739
0709
0585
0788
0752
0678
0731
0600
0707
0671
0673
0722
0729
Ca/(Ca þ Mg þ Fe)
0415
0368
0136
0374
0383 0094
0396
0368
0172
0423
0425
0423
0391
0116
0419
0039
0381
0420
0379
Mg/(Ca þ Mg þ Fe)
0455
0467
0546
0439
0433 0599
0446
0448
0484
0455
0433
0391
0445
0531
0410
0645
0417
0419
0453
Fe/(Ca þ Mg þ Fe)
0130
0165
0317
0187
0184 0307
0158
0184
0344
0122
0142
0186
0164
0353
0170
0316
0202
0161
0169
Total iron as FeO.
GM, groundmass. Numbers of ions are calculated on the basis of six oxygens.
Mixing between L1 and L2 (15% contribution of L1) produces a basaltic melt, L4, that has chemical characteristics
consistent with a melt component inferred from An80^90
plagioclase in the calc-alkalic rocks (Fig. 16a and Table 8).
A felsic melt component, L5, can be produced by 1:9
mixing of L1 and L3 (Fig. 16a and Table 8). Inferred intermediate calc-alkalic melt components can then be interpreted as mixing products between a felsic component
(L5), and a mafic component (L4) as shown in Fig. 16a
and Table 8.
Liquid mixing and cryptic mixing in
calc-alkalic magmas
Isotopic and elemental compositions can be used to
successfully identify the melt components that mixed to
produce the variety of calc-alkalic rocks in Zao volcano.
It should be stressed here that calc-alkalic andesites from
Zao volcano, whose end-member components are produced
by mixing between a calc-alkalic basaltic melt (L1 in
Fig. 16a), and tholeiitic basaltic and felsic melts (L2 and L3,
respectively), contain plagioclase crystallizing from mixed
end-member components (L4 and L5), but not from L1, L2,
nor L3 (see calc-alkalic andesite in Fig. 13). This observation
implies that mixing of melts or liquids, not magmas containing plagioclase crystals, and subsequent crystallization
of the mixed melts is the likely process that formed the
calc-alkalic end-member components. On the other hand,
the existing calc-alkalic andesites are produced by mixing
of these end-member ‘magmas’, which contain plagioclase
and other phenocryst phases.
The SiO2 contents and 87Sr/86Sr of the inferred mixed
magmas formed by the above processes can be calculated
1988
TATSUMI et al.
CALC-ALKALIC VS THOLEIITIC SERIES
Table 4: Representative compositions of orthopyroxene
Rock type:
Tholeiitic
Calc-alkalic
Sample:
ZA1011
Grain:
px032 px032 px002 px002 px009 px009 px004 px004 px004 px013 px013 px013 pxg017 px003 px003 px012 px012 pxg08 pxg12
Position:
Core
ZA1013
Rim
Core
ZA1044
Rim
Core
ZA3064
Rim
ZA3081
Core Mantle Rim
Core
Mantle Rim
GM
Core
Rim
Core
Rim
GM
GM
SiO2
5524
5511
5424
5466
5445
5512
5275
5456
5475
5470
5346
5333
5332
5416
5293
TiO2
022
025
024
018
019
036
021
012
033
020
024
022
025
020
020
015
016
037
038
Al2O3
136
133
143
121
133
086
062
082
073
165
092
076
131
062
066
130
090
254
271
Cr2O3
002
003
001
005
003
005
000
007
001
014
001
000
000
001
001
005
007
008
002
FeO
1600
1760
1842
1865
1779
2195 2370
1457
2193
1455
1646
1917
1863
2407
2440
1472 1608
1566
1940
MnO
043
043
042
045
024
MgO
2618
2520
2483
2437
2481
5325 5266
079
039
066
037
043
043
048
078
080
1955 2011
034
2749
1924
2779
2643
2453
2459
2092
2033
5568 5411
037
005
006
2751 2653
033
2579
2298
CaO
188
193
188
200
180
453
134
171
438
157
173
166
171
142
135
194
170
216
177
Na2O
002
004
001
002
004
004
006
004
007
003
004
005
004
001
004
003
002
004
004
Total
10135 10192 10148 10159 10068 10093 9949 10033 10010 10086 10101 10152 10047 10136 10111 10171 9994 10085 10029
Si
1973
1972
1959
1973
1973
1981 1991
1977
1983
1947
1968
1979
1955
1981
1987
1969 1964
1942
1940
Ti
0057
0056
0061
0051
0057
0038 0028
0035
0032
0069
0039
0032
0056
0027
0029
0054 0038
0107
0117
Al
0006
0007
0007
0005
0005
0010 0006
0003
0009
0005
0007
0006
0007
0006
0006
0004 0004
0010
0010
Cr
0001
0001
0000
0001
0001
0001 0000
0002
0000
0004
0000
0000
0000
0000
0000
0001 0002
0002
0001
Fe
0478
0526
0556
0563
0539
0683 0749
0437
0689
0434
0495
0580
0569
0747
0761
0435 0488
0469
0595
Mn
0013
0013
0013
0014
0007
0011 0025
0012
0021
0011
0013
0013
0015
0025
0025
0010 0011
0001
0002
Mg
1393
1344
1336
1311
1340
1084 1133
1469
1077
1478
1416
1323
1340
1158
1129
1450 1435
1378
1255
Ca
0072
0074
0073
0077
0070
0180 0054
0066
0176
0060
0067
0064
0067
0056
0054
0073 0066
0083
0070
Na
0002
0003
0001
0001
0003
0003 0005
0003
0005
0002
0003
0004
0003
0001
0003
0002 0001
0003
0003
Total
3995
3996
4006
3996
3995
3991 3991
4004
3992
4010
4008
4001
4012
4001
3994
3998 4009
3995
3993
Mg/(Mg þ Fe)
0745
0718
0706
0700
0713
0614 0602
0771
0610
0773
0741
0695
0702
0608
0598
0769 0746
0746
0679
Ca/(Ca þ Mg þ Fe)
0037
0038
0037
0039
0036
0092 0028
0033
0091
0030
0034
0033
0034
0029
0028
0037 0033
0043
0036
Mg/(Ca þ Mg þ Fe)
0717
0691
0680
0672
0688
0557 0585
0745
0555
0749
0716
0673
0678
0591
0581
0741 0721
0714
0654
Fe/(Ca þ Mg þ Fe)
0246
0271
0283
0289
0277
0351 0387
0222
0355
0220
0250
0295
0288
0381
0391
0222 0245
0243
0310
Total iron as FeO.
GM, groundmass. Numbers of ions are calculated on the basis of six oxygens.
using the compositions assumed for L1, L2, and L3 (Table 8
and Fig. 16b), and are plotted along the mixing curve
between L4 and L5 in Fig. 16b. The existing calc-alkalic
andesites have 87Sr/86Sr lower than the inferred mixed
magma (Fig. 16b). One possible mechanism to explain
this apparent dilemma would be that an above-liquidus
liquid, not a sub-liquidus magma, with low 87Sr/86Sr, such
as L1 (Fig. 16b), contributes to the existing calc-alkalic
andesite formation. Zao calc-alkalic rocks, however,
show no petrographic signs of such mixing. We thus propose to call this process ‘cryptic’ mixing of basaltic liquid
or melt. A calc-alkalic magma that experiences cryptic
mixing should not crystallize plagioclase phenocrysts.
If this is the case, then such cryptic mixing must have
taken place immediately before and could have caused
the eruption.
One petrographical observation that characterizes calcalkalic or Sakuyama’s R-type andesites is the presence of
pyroxenes in the groundmass that record higher temperatures than the phenocryst pyroxenes (e.g. Sakuyama, 1981,
1983; Wada, 1985; Fujinawa, 1988, 1990); this is the case for
Zao calc-alkalic rocks (Fig. 15). This observation can be
also explained by the cryptic mixing of a higher-T basaltic
melt with a lower-T magma.
Genesis of tholeiitic vs calc-alkalic basaltic
magmas: different fluid contributions
Petrological and geochemical observations suggest that at
least two basaltic magmas are simultaneously present in
the magma plumbing system of Zao volcano. One is a
magma that differentiates to form the tholeiitic series, and
the other a mafic end-member magma that mixes to form
1989
Rock type:
Tholeiitic
Sample:
ZA1011
Calc-alkalic
ZA1032
ZA3093
ZA3062: Andesite I
Grain:
bp002
bp002
ap001
ap001
cp001
cp001
cp003
cp003
bp001
bp001
ap002
ap002
cp001
cp001
cp001
cp001
Position:
Core
Rim
Core
Rim
Core
Rim
Core
Rim
Core
Rim
Core
Rim
Core
Rim
Core
Rim
4726
4505
4477
4438
4949
5355
5544
4626
4643
5040
5272
4652
4795
5467
000
001
000
000
002
005
004
005
000
001
001
007
000
004
001
5702
000
Al2O3
3469
3320
3418
3411
3423
3168
2920
2674
3389
3307
3040
2884
3281
3251
2876
2678
FeO
054
065
044
048
001
000
002
003
055
053
082
111
093
070
047
068
MnO
015
051
087
012
001
000
002
000
000
051
000
006
000
020
000
044
MgO
007
008
007
006
005
011
005
014
009
008
008
007
002
003
005
006
CaO
1861
1759
1862
1898
1903
1591
1267
1100
1794
1791
1481
1264
1807
1664
1238
1053
Na2O
060
144
067
064
061
258
426
437
080
108
307
405
105
177
449
513
K2O
002
001
000
002
001
003
009
023
004
005
011
025
002
007
017
028
Total
9949
10075
9990
9918
9835
9985
9990
9800
9957
9967
9970
9981
9942
9991
10100
10092
2082
2166
2091
2090
2084
2265
2425
2540
2139
2152
2314
2407
2161
2207
2452
2551
Ti
0000
0000
0000
0000
0001
0002
0001
0002
0000
0000
0000
0002
0000
0001
0000
0000
Al
1900
1793
1869
1876
1894
1709
1559
1443
1846
1806
1644
1551
1796
1763
1520
1411
Fe
0021
0025
0017
0019
0000
0000
0001
0001
0021
0020
0032
0042
0036
0027
0017
0026
Mn
0006
0020
0034
0005
0000
0000
0001
0000
0000
0020
0000
0002
0000
0008
0000
0017
Mg
0005
0006
0005
0004
0003
0008
0003
0010
0006
0005
0005
0005
0002
0002
0003
0004
Ca
0926
0863
0926
0949
0957
0780
0614
0540
0888
0889
0728
0618
0899
0820
0594
0505
Na
0054
0128
0061
0058
0056
0229
0374
0388
0071
0097
0273
0358
0095
0158
0390
0445
K
0001
0001
0000
0001
0000
0002
0005
0013
0003
0003
0006
0014
0001
0004
0010
0016
Total
4995
5002
5003
5002
4995
4995
4983
4937
4974
4992
5002
4999
4990
4990
4986
4975
O
8000
8000
8000
8000
8000
8000
8000
8000
8000
8000
8000
8000
8000
8000
8000
8000
Ca/(Ca þ Na)
0945
0871
0939
0942
0945
0773
0622
0582
0926
0902
0727
0633
0905
0839
0604
0532
Rb (ppm)
n.d.
0156
n.d.
0011
n.d.
0803
1420
2002
0031
0026
0127
4656
0015
3263
9061
1939
NOVEMBER 2008
Si
NUMBER 11
4481
TiO2
VOLUME 49
1990
SiO2
JOURNAL OF PETROLOGY
Table 5: Representative compositions of plagioclase phenocrysts
Sr
306
390
337
350
364
422
482
475
682
700
463
422
487
501
599
601
Y
0051
0145
0108
0161
0037
4903
0966
0598
0105
0105
0381
3575
0159
11593
3274
1402
Zr
0005
0259
0178
0262
0003
7914
2464
2025
0001
0006
0037
16299
0026
27001
21002
3996
Nb
n.d.
0009
0007
0004
n.d.
0559
0020
0180
n.d.
n.d.
0007
0730
0701
0782
Ba
504
487
585
436
3259
3548
La
0063
0308
0089
0129
0064
4262
1596
1457
0463
0552
1961
2688
0180
7970
6233
6598
Ce
0167
0534
0192
0267
0159
5788
3734
3405
0965
0993
3657
6018
0426
18418
12603
12409
1565
3610
10018
9081
7755
11493
n.d.
1223
17666
0236
25525
33913
0062
0020
0036
0018
1075
0380
0330
0092
0115
0378
0768
0054
1960
1301
1173
0247
n.d.
0162
0126
4548
1523
1229
0428
0447
1609
3274
0179
8068
4851
4223
Sm
0006
0056
0034
0052
0026
1009
0251
0309
0086
0102
0207
0792
0038
1994
0860
0580
Eu
0097
0191
0106
0108
0098
0584
1004
0853
0227
0204
0516
0802
0132
0710
1688
1978
Gd
0025
0047
0020
0025
0010
0883
0289
0217
0091
0110
0269
0732
0065
2131
1148
0914
Tb
0001
0008
0006
0006
0000
0132
0034
0022
0006
0010
0025
0114
0007
0367
0121
0071
Dy
0017
0038
0011
0036
0001
0823
0208
0127
0031
0039
0104
0733
0033
2353
0668
0405
Ho
n.d.
0003
0004
0003
0003
0161
0038
0024
0008
0009
0014
0124
0010
0436
0117
0057
Er
0006
0022
0006
0015
0008
0435
0120
0057
0003
0010
0029
0362
0023
1354
0352
0179
Tm
0002
0003
0001
0000
n.d.
0662
0009
0009
n.d.
0001
0003
0046
0000
0177
0054
0015
Yb
n.d.
0016
0004
0017
0008
0412
0114
0052
n.d.
0006
0016
0334
0017
1360
0358
0107
Lu
0001
0002
n.d.
n.d.
n.d.
0058
0014
0005
0003
n.d.
0001
0045
0002
0195
0059
0018
Hf
0002
0012
0011
n.d.
n.d.
0252
0060
0052
0002
n.d.
n.d.
0447
n.d.
0724
0582
0130
Ta
0001
0007
n.d.
n.d.
n.d.
0016
0007
0006
0001
n.d.
0001
0032
0003
0055
0057
0010
Pb
0113
0359
0161
0138
0112
0972
7406
6329
0291
0189
1533
1887
0407
4285
12205
6797
Th
n.d.
n.d.
0023
0004
n.d.
0119
0045
0039
n.d.
n.d.
0003
0584
n.d.
1044
0962
0169
U
0002
0003
0018
0001
0002
0053
0026
0033
0001
0001
0003
0150
n.d.
0658
0263
0073
0704204
0704486
0704374
0704469
0704549
070456
0703556
0703463
0704125
070398
0704165
0704193
0704094
87
Sr/86Sr
67
87
Sr/86Sr
0704119
15
85
89
0704202
40
0704328
25
0704260
15
71
0704396
42
61
82
70
0703418
33
72
80
0704033
26
75
59
0704205
20
54
57
0703971
17
(continued)
CALC-ALKALIC VS THOLEIITIC SERIES
0022
0102
TATSUMI et al.
1991
Pr
Nd
Rock type:
Calc-alkalic
Sample:
ZA3023: Andesite II
ZA3024: Andesite II
Grain:
ap003
ap003
ap004
ap004
ap004
ap004
ap005
ap005
ap007a
ap007a
Position:
Core
Rim
Core
Rim
Core
Rim
Core
Rim
Core
Rim
5418
5469
5093
4535
5576
5540
5577
5424
003
004
001
002
002
000
002
000
003
5657
001
Al2O3
3233
2797
2766
2996
3416
2738
2762
2735
2821
2684
FeO
048
042
041
044
000
003
001
001
000
001
MnO
040
000
000
033
000
002
000
001
000
000
MgO
003
003
004
003
002
002
000
000
004
003
CaO
1648
1172
1099
1417
1817
1021
1051
1024
1157
963
Na2O
162
492
491
348
096
533
514
539
477
560
K2O
004
026
026
015
002
030
027
029
023
035
Total
9879
9954
9897
9951
9870
9905
9897
9906
9909
9904
2203
2466
2495
2339
2115
2530
2517
2531
2471
Ti
0001
0001
0000
0001
0001
0000
0001
0000
0001
2562
0000
Al
1772
1500
1487
1621
1877
1464
1479
1463
1515
1432
Fe
0019
0016
0016
0017
0000
0001
0000
0000
0000
0000
Mn
0016
0000
0000
0013
0000
0001
0000
0000
0000
0000
Mg
0002
0002
0003
0002
0001
0001
0000
0000
0003
0002
Ca
0821
0571
0537
0697
0908
0496
0511
0498
0564
0467
Na
0146
0434
0434
0310
0087
0469
0452
0474
0421
0491
K
0002
0015
0015
0008
0001
0018
0016
0016
0013
0020
Total
4982
5005
4987
5008
4990
4980
4976
4982
4988
4974
O
8000
8000
8000
8000
8000
8000
8000
8000
8000
8000
Ca/(Ca þ Na)
0849
0568
0553
0692
0913
0514
0531
0512
0572
0487
Rb (ppm)
0405
6069
0326
0798
n.d.
0402
0402
0620
0344
0337
NOVEMBER 2008
Si
NUMBER 11
4738
TiO2
VOLUME 49
1992
SiO2
JOURNAL OF PETROLOGY
Table 5: Continued
Sr
584
507
480
642
389
490
477
466
545
490
Y
1198
0982
0875
1158
0060
0901
0785
1003
0908
0972
Zr
0499
2576
0150
1539
0017
0030
0119
0421
0536
0035
Nb
0031
0112
0004
0016
0001
n.d.
0013
0008
0020
Ba
14344
16266
18977
26710
763
29171
23557
29742
25335
0008
34454
3637
4008
6771
0127
5274
4887
5173
5065
6204
7117
7184
11140
0220
8611
9540
9595
7909
9609
Pr
0838
0699
0729
1105
0025
0883
0885
0873
0896
0961
Nd
2929
2411
2602
3951
0154
3177
3044
3537
3542
3824
Sm
0429
0360
0423
0588
0031
0521
0432
0465
0559
0553
Eu
1196
0980
1296
1675
0113
1701
1309
1600
1777
1919
Gd
0622
0524
0575
0911
0026
0751
0700
0815
0829
1044
Tb
0044
0035
0037
0051
0000
0043
0037
0063
0054
0050
Dy
0205
0171
0187
0275
0017
0239
0262
0243
0322
0248
Ho
0032
0031
0024
0032
0001
0026
0024
0044
0041
0033
Er
0096
0078
0055
0080
0006
0043
0080
0088
0125
0087
Tm
0010
0010
0006
0008
n.d.
0005
0008
0007
0011
0009
Yb
0081
0066
0037
0038
0012
0043
0060
0061
0046
0045
Lu
0013
0008
0004
0007
n.d.
0003
0005
0006
0005
0003
Hf
0027
0086
0019
0068
n.d.
0046
0042
0021
0062
0048
Ta
0001
0007
0002
0007
0000
0003
0000
0002
n.d.
n.d.
Pb
3894
3673
3794
5263
0399
5116
6187
7741
6013
6863
Th
0008
0101
0003
0038
0008
0002
0017
0085
0032
0013
U
0003
0027
0002
0010
0001
0006
0005
0015
0004
0001
0704198 62
0704144 63
0704056 69
0704058 65
0704042 38
0704052 68
0704007 67
0704164 62
0704079 44
0704168 57
87
Sr/86Sr
87
Sr/86Sr
LA-MC-ICP-MS,
0704018 19
MM-TIMS
0704035 29
0703945 17
0703942 33
CALC-ALKALIC VS THOLEIITIC SERIES
4339
8671
TATSUMI et al.
1993
La
Ce
JOURNAL OF PETROLOGY
VOLUME 49
NUMBER 11
NOVEMBER 2008
Table 6: Representative compositions of Fe^Ti oxides
Rock type:
Tholeiitic
Calc-alkalic
Sample:
ZA1011
Grain:
OQ-007
OQ-007c2
OQ-7
OQ-445
OQ-4
OQ-10
OQG-25
OQ-21
OQG-10
OQ-32
OQ-55
OQG-95
OQ-61
OQG-79
Mineral:
MT
MT
MT
MT
MT
MT
MT
IL
IL
MT
MT
MT
IL
IL
Position:
GM
GM
GM
GM
Core
Core
GM
Core
GM
Core
Core
GM
Core
GM
ZA1044
ZA3081
ZA3093
SiO2
012
020
016
012
005
006
010
001
007
005
005
007
007
011
TiO2
1181
1185
1564
1341
1099
1496
1564
4777
4813
1237
1219
981
4760
4261
Al2O3
381
302
214
281
149
176
118
007
010
222
224
161
009
011
FeO
7670
7733
7559
7391
8137
7755
7757
4934
4882
7621
7737
7993
4622
5302
MnO
025
032
034
042
035
046
040
066
053
033
037
045
126
039
MgO
267
180
089
148
092
101
105
148
147
221
238
090
302
086
Cr2O3
005
009
033
374
008
004
008
003
000
111
096
068
029
005
Total
9540
9463
9510
9589
9525
9584
9601
9936
9911
9449
9556
9345
9855
9715
Si
0004
0007
0006
0004
0002
0002
0004
0000
0002
0002
0002
0003
0002
0003
Ti
0325
0333
0444
0374
0311
0421
0441
0899
0906
0348
0339
0282
0896
0802
Al
0165
0133
0095
0123
0066
0077
0052
0002
0003
0098
0098
0073
0003
0003
Fe3þ
1174
1184
0995
1010
1307
1075
1057
0202
0180
1169
1193
1337
0209
0345
Fe2þ
1176
1229
1389
1284
1250
1352
1373
0830
0842
1216
1198
1219
0758
0764
Mn
0008
0010
0011
0013
0011
0015
0013
0014
0011
0010
0012
0015
0027
0008
Mg
0146
0100
0050
0082
0051
0056
0059
0055
0055
0123
0131
0051
0113
0032
Cr
0002
0003
0010
0110
0002
0001
0002
0001
0000
0033
0028
0020
0006
0001
Total
3000
3000
3000
3000
3000
3000
3000
2004
1998
3000
3000
3000
2013
1959
Total iron as FeO.
GM, groundmass; MT, titanomagnetite; IL, ilmenite. Fe3þ and Fe2þ are calculated assuming Fe–Ti oxides stoichiometry.
Numbers of ions are calculated on the basis of four oxygens.
the calc-alkalic series rocks. Although basaltic andesites in
Zao volcano are differentiated and do not represent primitive magma compositions, it is reasonable to assume that
calcic plagioclase (An490) in both tholeiitic and calcalkalic mafic andesites crystallizes from different primitive
magmas, as such plagioclase phenocrysts in the two rock
series exhibit different concentrations of some trace elements, especially of Sr (Fig. 12), If so, then the geochemical
characteristics of such primitive magmas, one tholeiitic
and other calc-alkalic, can be inferred from the plagioclase
compositions and element partitioning between plagioclase and silicate melt (Table 9 and Fig. 17).
Calc-alkalic rocks tend to be more enriched in incompatible trace elements than tholeiitic rocks, which is also
the case for the inferred primitive magma compositions
(Figs 7 and 17). This observation led Masuda & Aoki
(1979) to the conclusion that both calc-alkalic and tholeiitic
primary magmas tap a common mantle source but the
former is produced by lower degrees of melting. However,
this study on volcanic rocks from Zao volcano indicates
that the tholeiitic primitive magma has a higher 87Sr/86Sr
than the calc-alkalic magma, suggesting a more radiogenic
mantle source for the tholeiitic magmas. Such isotopic
heterogeneity within the mantle wedge could be caused
by larger contributions from isotopically enriched, slabderived fluids to the tholeiitic magma source than to the
calc-alkalic source.
To test this hypothesis quantitatively, geochemical modeling of element transport by slab-derived fluid addition to
the mantle wedge and subsequent incremental fractional
melting was conducted (Table 10 and Fig. 18). The model
compositions, including H2O contents of subducted altered
oceanic crust and terrigenous sediments, are from Tatsumi
& Hanyu (2003). Experimental data for the mobility of
Sr during dehydration of amphibolite (Kogiso et al., 1997)
and sediment (Aizawa et al., 1999) were used to estimate
the fluid compositions. Sr concentrations (357 and
664 ppm) and 87Sr/86Sr (070420 and 070343) in tholeiitic
and calc-alkalic primary magmas, respectively, are based
on the average compositions of melts in equilibrium
with An490 plagioclase phenocrysts (Fig. 13). The results
(Table 10 and Fig. 18) indicate that the geochemical
1994
TATSUMI et al.
CALC-ALKALIC VS THOLEIITIC SERIES
0.20
Table 7: Representative compositions of spinel inclusions
Rock type:
Tholeiitic
Calc-alkalic
Sample:
ZA1011
ZA3093
Grain:
OQ-10c1
OQ-10c2
OQ-7
OQ-445
Position:
Core
Core
Core
Core
Host mineral:
Cpx
Cpx
Ol
Ol
NiO (wt.%)
0.15
TH 1011
TH 1013
TH 1044
CA 3064
CA 3081
CA 3093
0.10
0.05
SiO2
006
008
010
010
TiO2
147
154
076
048
Al2O3
1961
1941
2821
3611
FeO
4334
4238
4383
3257
MnO
028
028
039
028
MgO
956
969
754
1236
Cr2O3
2447
2480
1785
1864
Total
9879
9816
9868
10053
0
60
65
70
75
80
85
100×Mg/(Mg+Fe)
Fig. 8. Compositions of olivine phenocrysts in tholeiitic (TH) and
calc-alkalic (CA) rocks of Zao volcano. Olivines in calc-alkalic rocks
tend to be more enriched in NiO than those in tholeiitic rocks.
Si
0002
0002
0003
0003
Ti
0035
0037
0018
0010
Al
0735
0732
1036
1232
Fe3þ
0576
0563
0482
0315
Fe2þ
0576
0570
0660
0473
Mn
0008
0008
0010
0007
Mg
0453
0462
0350
0533
Cr
0615
0627
0440
0427
Total
3000
3000
3000
3000
Cr/(Cr þ Al)
0456
0461
0298
0257
Mg/(Mg þ Fe2þ)
0440
0447
0347
0530
Cr/(Al þ Cr þ Fe3þ)
0319
0326
0225
0216
Al/(Al þ Cr þ Fe3þ)
0382
0381
0529
0624
Fe3þ/(Al þ Cr þ Fe3þ)
0299
0293
0246
0160
Total iron as FeO.
Cpx, clinopyroxene; Ol, olivine. Fe3þ and Fe2þ are calculated assuming spinel stoichiometry. Numbers of ions are
calculated on the basis of four oxygens.
characteristics of the tholeiitic and calc-alkalic primary
magmas can be explained by different contributions from
slab-derived fluids (1% vs 02%) and different degrees of
melting (40% vs 7%). It is suggested that these processes
produce primary magmas containing 263^265 wt %
H2O, which are acceptable values for arc magmas
(e.g. Hauri et al., 2006).
Although the different contributions from enriched slabderived fluids, together with different degrees of melting,
could generate two types of primary magma in the
mantle wedge, this may not be a plausible process to produce the tholeiitic and calc-alkalic magmas of Zao volcano.
This process does not explain why magma mixing plays
a major role in the differentiation of calc-alkalic magmas
but not in the tholeiitic magmas; the two types of
mantle-derived basalt magmas should have an equal
chance of magma mixing.
Genesis of tholeiitic vs calc-alkalic basaltic
magmas: crust vs mantle melts
When mantle-derived basaltic magmas are underplated
and/or intruded into the arc crust they transfer heat into
the overlying and surrounding crust, which can lead to
partial melting of the wall-rocks (e.g. Hildreth, 1981; Raia
& Spera, 1997; Annen & Sparks, 2002). However, herein
lies a problem: whether or not basaltic magmas emplaced
at the base of the lower crust could transfer enough heat to
continue to cause crustal anatexis; heat transfer from the
basaltic magma to the crust would cause a rapid temperature drop in the magma, which would lead to the magma
being unable to further melt the crust (e.g. Marsh, 1989;
Petford & Gallagher, 2001). This problem may be overcome
if the temperature of the pre-existing crust is high enough;
that is, the basaltic magma intrudes where crustal temperatures are near the basalt solidus or the crust is in a
partially molten state (Couch et al., 2001; Tatsumi et al.,
2006). Numerical simulations of heat transfer (Annen &
Sparks, 2002; Annen et al., 2006) further suggest a model
in which mantle-derived basalts emplaced as a succession
of sills into the lower crust generate a deep crustal ‘hot
zone’ where differentiated melts are produced from two
distinct sources: crystallization of mantle-derived magma,
and melting of crustal rocks.
A ‘hot zone’ is likely to occur beneath the NE Japan
arc. This is believed to be the case for the following
two reasons. First, the temperature of the NE Japanese
mantle-derived magma is much higher than the solidus
temperature of the lower crust. High-pressure melting
experiments (Tatsumi et al., 1983) suggest that primary
1995
JOURNAL OF PETROLOGY
VOLUME 49
NUMBER 11
Tholeiitic
70
60
60
50
90
50
90
70
70
60
60
50
50
90
60
70
80
100×Mg/ (Mg+Fe)
core
100×Mg/ (Mg+Fe)
rim
50
90
80
90
3093
clinopyroxene
60
70
3093
orthopyroxene
80
90
90
80
70
70
60
60
50
90
50
90
80
3064
clinopyroxene
3064
orthopyroxene
80
70
70
60
60
50
90
50
90
80
3062
clinopyroxene
3062
orthopyroxene
80
70
70
60
60
50
90
50
90
3081
clinopyroxene
80
70
reversely
zoned
60
normally
zoned
50
50
90
60
70
80
100×Mg/ (Mg+Fe)
core
3081
orthopyroxene
80
70
60
100×Mg/ (Mg+Fe)
rim
60
70
80
100×Mg/ (Mg+Fe)
core
70
100×Mg/ (Mg+Fe)
rim
80
100×Mg/ (Mg+Fe)
rim
1044
orthopyroxene
1044
clinopyroxene
100×Mg/ (Mg+Fe)
rim
80
100×Mg/ (Mg+Fe)
rim
1013
orthopyroxene
80
60
100×Mg/ (Mg+Fe)
rim
1013
clinopyroxene
100×Mg/ (Mg+Fe)
rim
50
90
70
80
50
90
100×Mg/ (Mg+Fe)
rim
100×Mg/ (Mg+Fe)
rim
80
60
50
90
100×Mg/ (Mg+Fe)
rim
90
90
70
reversely
zoned
normally
zoned
80
Calc-alkalic
100×Mg/ (Mg+Fe)
100×Mg/ (Mg+Fe)
core
core
100×Mg/ (Mg+Fe)
rim
60
1011
orthopyroxene
1011
clinopyroxene
80
70
100×Mg/ (Mg+Fe)
core
90
60
70
80
100×Mg/ (Mg+Fe)
rim
100×Mg/ (Mg+Fe)
rim
50
90
100×Mg/ (Mg+Fe)
core
60
70
80
NOVEMBER 2008
50
60
70
80
90
100×Mg/ (Mg+Fe)
core
Fig. 9. Compositions of pyroxene phenocrysts in tholeiitic and calc-alkalic rocks from Zao volcano. Tholeiitic rocks contains normally zoned
pyroxenes, whereas calc-alkalic rocks are characterized by the occurrence of both normally and reversely zoned orthopyroxene.
magmas beneath the volcanic front of this arc equilibrate
with the mantle at 14008C and 10 GPa, a pressure equivalent to the depth immediately beneath the Moho underlying the NE Japan arc. Second, low-frequency tremors
and micro-earthquakes, which may be caused by deformation associated with magma intrusion, are observed
at depths of 30^50 km only beneath the Quaternary volcanoes of the volcanic front of the NE Japan arc (Obara,
2002; Katsumata & Kamaya, 2003).
The petrographic and geochemical data for the Zao volcanic rocks presented here suggest the presence of two distinct basaltic magmas beneath a single volcano; one a high
87
Sr/86Sr tholeiitic magma and the other a low 87Sr/86Sr
calc-alkalic magma. If the formation of a ‘hot zone’,
resulting from basaltic underplating and subsequent
generation of both crust-derived and mantle-derived
magmas, is accepted, then it is reasonable to suggest that
the isotopically enriched tholeiitic and depleted calcalkalic magmas may be crust- and mantle-derived,
respectively.
Contributions from melting of crustal and mantle materials to generate tholeiitic and calc-alkalic magmas, respectively, are now examined on the basis of their geochemical
characteristics. Relative abundances of trace elements
between tholeiitic and calc-alkalic primitive magmas,
which are based on solid^melt partitioning between
Ca-rich plagioclase and silicate melts (Table 9), are plotted
as a function of the ionic radii of elements in Fig. 19.
Systematic patterns of enrichment and depletion of certain
elements can be observed for the two magma series.
1996
TATSUMI et al.
CALC-ALKALIC VS THOLEIITIC SERIES
Tholeiitic
Calc-allkalic
100xCa/(Ca+Na)
Freq.
15
50
60
70
80
100xCa/(Ca+Na)
90 100
40
20
1011
15
Freq.
40
20
10
5
20
20
15
15
10
10
5
5
0
20
0
20
1013
15
10
10
5
5
20
20
15
15
10
10
5
5
0
0
20
20
15
1044
15
10
10
5
5
20
20
15
15
10
10
5
5
0
40
50
60
70
80
60
70
80
90 100
3093
10
5
15
50
3064
3062C
0
20
90 100
100xCa/(Ca+Na)
3081
15
10
5
core
20
rim &
groundmass
15
10
5
0
40
50
60
70
80
90 100
100xCa/(Ca+Na)
Fig. 10. Compositions of plagioclase in Zao volcanic rocks. The plagioclase phenocrysts in the tholeiitic rocks have a much narrower range of
core compositions than those in the calc-alkalic rocks. Furthermore, the cores of plagioclase phenocrysts in calc-alkalic rocks exhibit a bimodal
compositional distribution.
1997
JOURNAL OF PETROLOGY
VOLUME 49
NUMBER 11
NOVEMBER 2008
A
100
0.704272±68
100×Ca/ (Ca+Na)
100×Ca/ (Ca+Na)
0.704296±66
0.704325±81
0.704262±16
B
800
80
600
TH-ZA1032bp005
60
Sr (ppm)
400
40
1000µm
0
1000
3000 µm
2000
0.703418±33
0.703690±96
0.703556±70
0.703463±72
A
B
80
600
CA basaltic andesite I
-ZA3093bp001
60
40
1000µm
800
0
Sr (ppm)
100×Ca/ (Ca+Na)
A
100
B
Sr (ppm)
A
B
400
4000 µm
2000
A
B
100
800
0.704035±29
0.704042±38
80
600
CA andesite II
-ZA3024ap004
Sr (ppm)
100×Ca/ (Ca+Na)
B
B
60
400
0.704052±68
40
1000µm
A
A
100
0
400
800
1600 µm
1200
A
B
800
CA basaltic andesite III
0.703914±28
A
80
60
400
1000µm
40
100
0.704105±75
0.704058±65
0.704056±69
B
0
500
1000
1500
A
2000 µm
B
800
CA andesite IV
-ZA3023ap004
80
600
Sr (ppm)
0.7041
19±15
100×Ca/ (Ca+Na)
A
600
Sr (ppm)
0.704113±116
-ZA3093bp002
100×Ca/ (Ca+Na)
B
0.704047±76
60
400
40
1000µm
0
500
1000
1500
2000 µm
Fig. 11. Scanning electron micrographs for representative plagioclase phenocrysts and compositional profiles of 100 Ca/Ca þ Na) (continuous line) and Sr concentration (dashed line) across the plagioclase crystals. The numbers (I^IV) labeled for plagioclase in calc-alkalic rocks
correspond to those given based on 87Sr/86Sr and Sr concentration (see Fig. 13). Large open circles and italic values on the images represent
sampling spots and 87Sr/86Sr by MM-TIMS, and large closed circles and regular-font values those by LA-MC-ICP-MS. Small grey circles are
sampling spots for LA-ICP-MS trace element analyses.
1998
TATSUMI et al.
CALC-ALKALIC VS THOLEIITIC SERIES
100×Ca/ (Ca+Na)
40
400
50
60
70
80
100×Ca/ (Ca+Na)
90
100
40
TH
60
70
80
90
100
50
60
70
80
90
100
700
CA
300
600
Sr
Ba
50
200
500
400
100
300
0
2.5
20
2.0
15
Eu
Ce
1.5
10
1.0
0.5
0
0
16
4
12
3
8
2
Y
Pb
5
1
4
0
40
50
60
70
80
90
0
40
100
100×Ca/ (Ca+Na)
100×Ca/ (Ca+Na)
Fig. 12. Trace element concentrations in plagioclase phenocrysts in Zao volcanic rocks as a function of Ca/(Ca þ Na). Plagioclase phenocrysts
in calc-alkalic rocks tend to be more enriched in Ba and Sr and more depleted in Y than those in tholeiitic rocks.
This pattern and a consideration of the crystal structure
control on trace element partitioning between melts and
solid phases (Matsui et al.,1977) may suggest that the tholeiitic basalt magma is more depleted in elements that are likely
to be partitioned into plagioclase and possibly amphibole
than the calc-alkalic basalt magma (Fig. 19). This can be
understood as the result of buffering of particular elements
by residual phases during partial melting; for example, Sr
and Eu by plagioclase, and K and Ba by amphibole. If so,
then the melting residue of the tholeiitic basalt magma is
distinct from that of the calc-alkalic basalt magmas in the
presence of plagioclase (and amphibole). These phases are
commonly observed in the melting residues of hydrous
mafic compositions such as amphibolite (e.g. Beard &
1999
JOURNAL OF PETROLOGY
VOLUME 49
NUMBER 11
NOVEMBER 2008
Sr (ppm) in Melt
100×Ca/(Ca+Na)
40
50
60
70
0.7048
90
0
100
0.7048
CA bulk TH bulk
Tholeiitic
0.7040
100 200 300 400 500 600 700 800
Tholeiitic
0.7044
87Sr/86Sr
0.7044
87Sr/86Sr
80
in equilibrium with plagioclase
MM-TIMS
0.7036
0.7040
Core TIMS
Rim TIMS
Core ICPMS
Rim ICPMS
0.7036
LA-ICP-MS
0.7032
0.7032
0.7048
Calc-alkalic
0.7048
Calc-alkalic
basaltic andesite
basaltic andesite
0.7044
87Sr/86Sr
87Sr/86Sr
0.7044
II
0.7040
0.7040
II
III
IV
IV
0.7036
0.7036
I
III
I
0.7032
0.7032
0.7048
0.7048
Calc-alkalic
Calc-alkalic
andesite
andesite
87Sr/86Sr
87Sr/86Sr
0.7044
II
0.7040
0.7044
0.7040
II
III
0.7032
IV
IV
0.7036
0.7036
I
40
50
60
70
80
90
III
I
0.7032
100
100×Ca/(Ca+Na)
0
100 200 300 400 500 600 700 800
Sr (ppm) in Melt
in equilibrium with plagioclase
Fig. 13. 87Sr/86Sr of plagioclase phenocrysts in tholeiitic and calc-alkalic rocks from Zao volcano as functions of the anorthite content and Sr
concentration in melts inferred from Sr partitioning between plagioclase and silicate melts (Bindeman et al., 1998; Bindeman & Davis, 2000).
Plagioclase in calc-alkalic rocks tends to have lower 87Sr/86Sr than in tholeiitic rocks; more importantly, calcic plagioclase in the most mafic calcalkalic andesites crystallizes from magmas with characteristically low 87Sr/86Sr (07034) and high Sr up to 700 ppm (I). At least four components can be identified (I^IV), based on plagioclase core compositions, for production of mixed calc-alkalic rocks. Accuracy of MM-TIMS
analyses (2s) and precision of LA-ICP-MS analyses are shown by bars.
Lofgren, 1991; Beard et al., 1993; Patin‹o Douce & Beard,
1995), but are unlikely to coexist with primary arc basalt
magmas in a peridotite system (e.g. Tatsumi et al., 1983).
Therefore, it is possible that the geochemical characteristics
of the tholeiitic and calc-alkalic basaltic magma series in
Zao volcano can be understood if they are produced by partial melting of amphibolitic lower crust and peridotitic
upper mantle, respectively.
2000
TATSUMI et al.
CALC-ALKALIC VS THOLEIITIC SERIES
100×Mg/(Mg+Fe) in olivine
90
50
60
70
80
100×Mg/(Mg+Fe) in olivine
90
50
90
Tholeiitic
70
80
90
Calc-alkalic
101 1
1044
70
1013
60
3064
80
100×Mg/(Mg+Fe)
in clinopyroxene
100×Mg/(Mg+Fe)
in clinopyroxene
80
60
3081
70
60
3093
50
90
50
90
Tholeiitic
Calc-alkalic
100×Mg/(Mg+Fe)
in orthopyroxene
100×Mg/(Mg+Fe)
in orthopyroxene
80
70
80
70
60
60
50 50
60
70
80
50
50
90
100×Mg/(Mg+Fe) in olivine
60
70
80
90
100×Mg/(Mg+Fe) in olivine
Fig. 14. Relationship between the core composition of olivine and pyroxene phenocrysts in Zao volcanic rocks. The compositional range for
each sample (numbered) is shown by a box. Broadly identical Mg/(Mg þ Fe) values in olivine and pyroxenes for tholeiitic rocks suggest that
these phases are in equilibrium. The pyroxenes with low Mg/(Mg þ Fe) that tend to be reversely zoned (Fig. 9) are not in equilibrium with
olivine with higher Mg/(Mg þ Fe).
1200
Temperature (°C)
1150
Higher-T Rim
in CA
Higher-T Rim &
Groundmass
in CA
1100
1050
1000
950
900
850
50
TH core
CA core
CA rim
CA gm
52
54
56
58
60
62
64
SiO2 (wt.%) in bulk rock
Fig. 15. Temperature estimates based on two-pyroxene geothermometry (Wells, 1977). Calc-alkalic rocks are characterized by the
occurrence of pyroxene phenocrysts with cores that show lower
temperatures than their rims and pyroxenes in the groundmass. CA,
calc-alkalic; TH, tholeiitic.
Subduction zone tholeiitic magmas have been considered to form from mantle-derived basaltic magmas via differential crystallization (e.g. Wada, 1981; Sakuyama, 1983;
Fujinawa, 1988, 1990; Tatsumi & Kogiso, 2003) for the
following reasons: (1) tholeiitic rocks show little evidence
of disequilibrium textures; (2) they exhibit systematic
changes in both phenocryst compositions and assemblages; (3) the tholeiitic trend can be explained by the fractionation of phenocryst phases. On the other hand, these
petrographic and compositional characteristics are also
consistent with inverse differential crystallization of a
parental basaltic magma (i.e. partial melting of a basaltic
source) having affected the tholeiitic rocks. As a basaltic
parental magma for the Zao tholeiitic rocks is generated
by melting of mafic lower crust caused by heat transfer
from an underplating mantle-derived, calc-alkalic basaltic
magma, it is reasonable to suggest that differentiated
tholeiitic melts are also created via crustal anatexis
rather than crystallization of a mantle-derived basaltic
magma.
2001
JOURNAL OF PETROLOGY
VOLUME 49
Sr (ppm) in melt
in equilibrium with plagioclase
(a)
0.7046
0
100
200
300
400
500
87Sr/ 86Sr
0.7042
An50-60
in CA-A
L2
0.7040
0.7036
Magma type Plagioclase SiO2
0.7034
L1
0.7032
Tholeiitic basalt
Tholeiitic rhyolite
Calc-alkalic basalt
Calc-alkalic mafic endmember
(b)
87
Sr/86Sr
L1
calc-alkalic
An490
50
664
070342
L2
tholeiitic
An490
50
357
070425
L3
tholeiitic
An550
70
100
070425
Gr
upper crust
70
300
070512
L4 (¼ 015L1 þ 085L3) calc-alkalic
An90–80
50
403
070401
L5 (¼ 01L1 þ 09L2)
An50
68
156
070386
calc-alkalic
Calc-alkalic felsic endmember
Upper crust
Tholeiitic
Calc-alkalic
0.7055
n
atio
in
tam
on
0.7050
Cru
0.7045
L3
L2
0.7040
0.7035
0.7030
45
Gr
lC
sta
87Sr/ 86Sr
Sr
(wt %) (ppm)
An>90
in CA-BA
An~70
in CAA & BA
L5
800
An80-90
in CA-A
L4
0.7038
700
NOVEMBER 2008
Table 8: Characteristics of end-member and mixed
components for Zao magmas
An>90 in TH
0.7044
L3
600
NUMBER 11
L4
L5
L1
Cryptic Mixing
50
55
60
65
70
75
SiO2 (wt.%)
Fig. 16. Compositions of end-member components that contribute
to the production of Zao magmas. (a) Variations in 87Sr/86Sr and Sr
concentrations for the melt components inferred from the plagioclase
compositions (see Fig. 13). Three principal end-member components
(L1, L2, and L3) are required to explain the calc-alkalic melt components, L4 and L5 , which further mix to produce calc-alkalic melts that
are able to crystallize An-poor plagioclase. TH, tholeiitic; CA, calcalkalic; BA, basaltic andesite; A, andesite. (b) Relationship between
87
Sr/86Sr and SiO2 contents inferred for the end-member components
of Zao magmas. The compositional characteristics in tholeiitic rocks
can be explained by contamination from basement granitic rock
(Gr), whereas those in calc-alkalic rocks suggest a contribution from
basaltic melt, L1, which forms an end-member component for cryptic
mixing (see text).
Magma plumbing system beneath
Zao volcano
The model for the generation of the two types of magmas,
tholeiitic and calc-alkalic, beneath Zao volcano described
below is shown schematically in Fig. 20.
A mantle-derived basaltic magma (L1), which finally
equilibrates with the upper mantle immediately below the
Moho (10 GPa; Tatsumi et al., 1983), underplates and
transfers heat to the base of the crust, causing both its
own crystallization and partial melting of the lower crust.
Although the temperature and melt fraction gradually
decrease upwards within the partially molten ‘hot zone’
formed at the base of the lower crust, the hot zone is simplified to consist of two sub-zones with higher and lower
melt fractions that generate basaltic and felsic melts.
These are melts L2 and L3, respectively, inferred from the
isotopic compositions of plagioclase phenocrysts in the tholeiitic rocks (Figs 16 and 20). The boundary between the
two sub-zones may be defined by the breakdown and melting of amphibole, which causes an abrupt increase in melt
fraction (e.g. Foden & Green, 1992; Annen et al., 2006).
If so, then the temperature at this boundary would be
10758C (Mu«ntener et al., 2001). This temperature estimate
is consistent with the following petrographic and experimental constraints; first, the highest temperature estimate
for tholeiitic magmas based on a two-pyroxene geothermometer is 10758C (Fig. 15); second, experiments at
03 GPa on a basalt (SiO2 49 wt %; Al2O3 18 wt %; FeO
11wt %; MgO 7 wt %) in the presence of 05 wt %
H2O yield a partial melt with a composition similar to
that of the tholeiitic andesitic basalts (SiO2 52 wt %) at
11008C by 40^50% of partial melting (Tatsumi &
Suzuki, in preparation).
A variable contribution of L2 vs L3 through liquid^
liquid mixing yields melts with mafic to felsic compositions, which crystallize in shallow-level magma reservoirs
to form the tholeiitic rocks (Fig. 20). Alternatively, differential crystallization of L2 may also contribute to the production of tholeiitic magmas with intermediate compositions.
In addition to these processes, contamination of the
magmas by upper crustal granitic rocks plays a role in
controlling the Sr isotopic compositions of the tholeiitic
magmas (Fig. 20), which is suggested by (1) the observation
that plagioclase rims tend to have higher 87Sr/86Sr
than cores in tholeiitic rocks (Fig. 13) and (2) bulk compositions exhibiting higher 87Sr/86Sr with increasing SiO2
(Fig. 16).
2002
TATSUMI et al.
CALC-ALKALIC VS THOLEIITIC SERIES
Table 9: Trace element compositioins of calcic plagicolase and inferred primitive melt
Tholeiitic
Plagioclase
Rb
00779
K
76
Ba
578
Sr
357
Calc-alkalic
1s
0123
50
140
19
Melt
D
001
1s
902
Plagioclase
145
004–005
1649
1065
007–005
109
23
118–094
358
25
1s
00872
0106
250
001
86
376
827
005
77
691
Melt
D
24
1s
101
4780
1711
007–006
601
138
109–100
665
38
La
0102
0040
014–013
0779
0298
0552
0110
014
Ce
0215
0080
008–007
291
105
107
033
008
Pr
00274
00118
007–006
0428
0181
0118
0049
007
163
068
Nd
0124
0051
009–008
157
063
0461
0271
009
523
307
Sm
00301
00162
006–005
0534
0283
00971
00939
007–006
152
147
Eu
0118
0110
008–007
167
154
0234
0059
008
297
077
Pb
0206
0369
050–045
0783
0326
0141
049–048
Y
0152
0149
002–001
0302
0714
002–001
0442
122
123
3899
130
0674
198
0781
40
0297
467
Distribution coefficients (D) are after Bindeman et al. (1998) and Bindeman & Davis (2000).
N-MORB Normalized
100
Primitive CA
CA
Primitive TH
TH
Table 10: Modeling for primary magma generation via
different fluid contributions
87
Sr/86Sr
Sr
10
Fluid
Melt
H2O
fraction (f)
fraction (F)
(wt %)
Altered oceanic crust
amphibolite
1
mobilityy
fluidz
500
070450
041
13667
Sediment
terriginous
0.1
Rb Ba K
mobilityy
La Ce Pb Pr Sr Nd Sm Eu Y
Fig. 17. N-MORB-normalized trace element patterns for average
calc-alkalic (CA) and tholeiitic (TH) rocks from Zao volcano and
those inferred for primitive magma compositions. Calc-alkalic rocks
tend to be more enriched in incompatible trace elements than
tholeiitic rocks, which is also the case for inferred primitive magma
compositions.
250
fluidz
2000
Slab fluidx
11333
Original mantle
070600
012
070455
25
070250
Tholeiitic
357
070420
00105
04
263
Calc-alkalic
664
070343
000182
0068
265
The generation of calc-alkalic magmas is distinct from
tholeiitic magmas in that it involves a mantle-derived
basaltic component, L1, either as a liquid or magma
(Fig. 20). This contributes to the generation of mafic (L4)
and felsic (L5) end-member magmas for andesites via
mixing with crust-derived, tholeiitic basalt (L2) and felsic
(L3) melts. Furthermore, the calc-alkalic primitive liquid
L1 plays a role in the production of all calc-alkalic
Tatsumi & Hanyu (2003).
yKogiso et al. (1997) and Aizawa et al. (1999).
zH2O in amphibolite and sediment is assumed to be 15 wt %.
xTwenty per cent sediment fluid contribution.
magmas via cryptic mixing (Fig. 16). On the other hand,
a magma with the composition of L1, containing Ca-rich
and low 87Sr/86Sr (07034) plagioclase phenocrysts, can
be identified in the most mafic basaltic andesite in Zao
volcano (Figs 16 and 20).
2003
JOURNAL OF PETROLOGY
VOLUME 49
0.7050
0.7025
0.7020
1
Primitive Melt
Bulk Rock
amph
Slab-derived
fluid
5.0
0.005
0.003
0.002
Y
amph
Rb
plag
f =0.001
Original
mantle
wedge
10
102
103
Tholeiitic/Calc-alkalic
0.7030
0.02
0.01
F=0.4
0.2
0.1
0.05
87Sr/86Sr
0.7035
NOVEMBER 2008
10.0
0.7045
0.7040
NUMBER 11
Tholeiitic
Calc-alkalic
104
105
106
Sr (ppm)
Fig. 18. Compositions of magmas produced by different contributions
from slab-derived fluid (f) and different degrees of partial melting (F).
The geochemical characteristics of tholeiitic and calc-alkalic primary
magmas in Zao volcano can be explained by different contributions
from slab-derived fluids (1% vs 02%) and different degrees of melting (40% vs 7%).
Pb
3+
1+
2+
Sm
1.0
Eu
Nd
Pr
Ce
K
0.5
La Sr
Ba
4+
CONC LUSIONS
It is generally accepted that the differentiation of calcalkalic magmas involves variable contributions from crustal components; for example, via wall-rock assimilation, or
mixing with crust-derived felsic magma (Eichelberger,
1975; DePaolo, 1981; Sakuyama, 1981; Hildreth &
Moorbath, 1988; Clynne, 1999; Dungan & Davidson, 2004;
Tatsumi & Kogiso, 2003), whereas tholeiitic magmas show
more pristine mantle signatures. The hypothesis presented
here, which proposes a crustal origin for the tholeiitic
magmas and a mantle origin for the calc-alkalic basaltic
magmas, requires that these models be revised and even
in some cases discarded.
One aspect that we need to re-examine concerns the
geochemical characteristics of the ‘mantle-derived’ basalt
magmas that are used to understand the contribution
of slab-derived components to arc magma generation. To
minimize the effect of shallow-level crustal contamination
and to assess the magma source characteristics, basaltic
rocks and/or tholeiitic rocks tend to be examined (e.g.
Notsu, 1983; Sakuyama & Nesbitt, 1986; Shibata &
Nakamura, 1997; Kimura & Yoshida, 2006). We suggest
that the tholeiitic basalt magmas in Zao volcano are
derived from melting of mafic lower crust via underplating
of calc-alkalic, mantle-derived basalt magmas and subsequent crustal anatexis. If so, then the tholeiitic basalt,
although it is relatively primitive, does influence the geochemical signatures of Quaternary mantle-derived arc
magmas, contributing to the production of calc-alkalic
magmas. However, only by examining the phenocryst
phases that crystallized from the least differentiated
0.1
60
80
100
120
140
160
Ionic Radius (pm)
Fig. 19. Relative abundances of elements between tholeiitic and
calc-alkalic primitive magmas in Zao volcano as a function of
ionic radius. The tholeiitic magma is more depleted in elements that
are likely to be partitioned into plagioclase (plag) and possibly
amphibole (amph) than the calc-alkalic basalt magma, suggesting
the presence of these phases in the melting residue of the tholeiitic
magma.
mantle-derived magmas can the effect of shallow-level
magma mixing processes on the mantle signatures be
hinted at. Therefore, analysis and examination of the compositions of minerals that crystallize from the primitive
calc-alkalic basalt magma could provide the only chance
to fully understand the geochemical characteristics of a
mantle-derived magma, and hence the source mantle and
slab-derived components.
AC K N O W L E D G E M E N T S
We thank Bogdan Vaglarov for analytical assistance, Miki
Fukuda for preparing the manuscript and figures, and
Richard Price, Bruce Charlier, Gene Yogodzinski, Alex
Nichols, and the editor John Gamble for their critical and
constructive comments on the manuscript. This work is
partially supported by Grant-in-Aid for Creative
Scientific Research (19GS0211).
2004
TATSUMI et al.
CALC-ALKALIC VS THOLEIITIC SERIES
Tholeiitic
Calc-alkalic
Basaltic Andesitic
Andesitic Basaltic
Mixing
Crystallization
Contamination
L5
Partially Molten
Hot Zone
L4
Crystallizing
Magma Reservoirs
Lower-T Partial Melt L3
+ Restite
Higher-T Partial Melt L2
+ Restite
Underplating & Crystallizing
Basaltic Magma (L1)
that causes crustal anatexis
Upper Mantle
Contribution of melt
Contribution of melt+crystal (magma)
Fig. 20. A schematic model for the magma plumbing system and the magma differentiation process in Zao volcano. The compositional
characteristics of melts (L1 to L5) contributing to tholeiitic and calc-alkalic magmas have been given in Figs 13 and 16a and Table 10.
R EF ER ENC ES
Aizawa, Y., Tatsumi, Y. & Yamada, H. (1999). Element transport
during dehydration of subducting sediments: implication for arc
and ocean island magmatism. Island Arc 8, 38^46.
Annen, C. & Sparks, R. S. J. (2002). Effects of repetitive emplacement
of basaltic intrusions on thermal evolution and melt generation in
the crust. Earth and Planetary Science Letters 203, 937^955.
Annen, C., Blundy, J. D. & Sparks, R. S. J. (2006). The genesis of intermediate and silicic magmas in deep crustal hot zones. Journal of
Petrology 47, 505^539.
Baker, M. B., Grove, T. L. & Price, R. (1994). Primitive basalts and
andesites from the Mt. Shasta region, N. California: products of
varying melt fraction and water content. Contributions to Mineralogy
and Petrology 118, 111^129.
Baker, W. E. (1973). The role of humic acids from Tasmanian podzolic
soils in mineral degradation and metal mobilization. Geochimica et
Cosmochimica Acta 37, 269^281.
Barling, J. & Weis, D. (2008). Influence of non-spectral matrix effects
on the accuracy of Pb isotope ratio measurement by MC-ICP-MS:
implications for the external normalization method of instrumental
mass bias correction. Journal of Analytical Atomic Spectrometry 23,
1017^1025.
Bau, M. & Knittel, U. (1993). Significance of slab-derived partial
melts and aqueous fluids for the genesis of tholeiitic and calcalkaline island-arc basalts; evidence from Mt. Arayat, Philippines.
Chemical Geology 105, 233^251.
Beard, J. S. & Lofgren, G. E. (1991). Dehydration melting and watersaturated melting of basaltic and andesitic greenstones and
amphibolites at 1, 3 and 69 kb. Journal of Petrology 32, 365^401.
Beard, J. S., Abitz, R. J. & Lofgren, G. E. (1993). Experimental melting of crustal xenoliths from Kilbourne Hole, New Mexico and
implications for the contamination and genesis of magmas.
Contributions to Mineralogy and Petrology 115, 88^102.
Bindeman, I. N. & Davis, A. M. (2000). Trace element partitioning
between plagioclase and melt; investigation of dopant influence on
partition behavior. Geochimica et Cosmochimica Acta 64, 2863^2878.
Bindeman, I. N., Davis, A. M. & Drake, M. J. (1998). Ion microprobe
study of plagioclase^basalt partition experiments at natural
concentration levels of trace elements. Geochimica et Cosmochimica
Acta 62, 1175^1193.
Blatter, D. L. & Carmichael, I. S. E. (2001). Hydrous phase equilibria
of a Mexican high-silica andesite: A candidate for a mantle origin?
Geochimica et Cosmochimica Acta 65, 4043^4065.
Bowen, N. L. (1928). The Evolution of the Igneous Rocks. Princeton, NJ:
Princeton University Press.
Brophy, J. G. (1990). Andesites from northeastern Kanga Island,
Aleutiansçimplication for calc-alkaline fractionation mechanisms
and magma chamber development. Contributions to Mineralogy and
Petrology 104, 568^581.
Chang, Q., Shibata, T., Shinotsuka, K., Yoshikawa, M. &
Tatsumi, Y. (2003). Precise determination of trace elements in
geological standard rocks using inductively coupled plasma mass
spectrometry (ICP-MS). Frontier Research on Earth Evolution 1,
357^362.
Charlier, B. L. A., Ginibre, C., Morgan, D., Nowell, G. M.,
Pearson, D. G., Davidson, J. P. & Ottley, C. J. (2006). Methods for
the microsampling and high-precision analysis of strontium and
rubidium isotope at single crystal scale for petrological and geochronological applications. Chemical Geology 232, 114^133.
2005
JOURNAL OF PETROLOGY
VOLUME 49
Clynne, M. A. (1999). A complex magma mixing origin for rocks
erupted in 1915, Lassen Peak, California. Journal of Petrology 40,
105^132.
Couch, S., Sparks, R. S. J. & Carroll, M. R. (2001). Mineral disequilibrium in lavas explained by convective self-mixing in open magma
chambers. Nature 411, 1037^1039.
Crawford, A. J., Falloon, T. J. & Green, D. H. (1989). Classification,
petrogenesis and tectonic setting of boninites. In: Crawford, A. J.
(ed.) Boninites and Related Rocks. London: Unwin Hyman, pp. 1^49.
Davidson, J. P. & Tepley, F. J. (1997). Recharge in volcanic systems:
evidence from isotope profiles of phenocrysts. Science 275, 826^829.
Davidson, J. P., Tepley, F. J., Palacz, Z. & Meffan-Main, S. (2001).
Magma recharge, contamination and residence times revealed by
in situ laser ablation isotopic analysis of feldspar in volcanic rocks.
Earth and Planetary Science Letters 184, 427^442.
Defant, M. J. & Drummond, M. S. (1990). Derivation of some modern
arc magmas by melting of young subducted lithosphere. Nature 347,
662^665.
DePaolo, D. J. (1981). Trace element and isotopic effects of combined
wallrock assimilation and fractional crystallization. Earth and
Planetary Science Letters 53, 189^202.
Dickinson, W. R. (1975). Potash^depth (K^h) relations in continental
margin and intra-oceanic magmatic arcs. Geology 3, 53^56.
Dungan, M. A. & Davidson, J. (2004). Partial assimilative recycling of
the mafic plutonic roots of arc volcanoes: An example from the
Chilean Andes. Geology 32, 773^776.
Eichelberger, J. C. (1975). Origin of andesite and dacite: evidence of
mixing at Glass Mountain in California and at other CircumPacific volcanoes. Geological Society of America Bulletin 86, 1381^1391.
Ewart, A. (1982). The mineralogy and petrology of Tertiary^Recent
orogenic volcanic rocks with special reference to the andesitic^
basaltic compositional range. In: Thorpe, R. S. (ed.) Andesites.
Chichester: John Wiley, pp. 26^87.
Fenner, C. N. (1929). The crystallization of basalts. American Journal of
Science 18, 225^253.
Foden, 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, 479^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 relations. Lithos 24, 217^236.
George, R., Turner, S., Hawkesworth, C., Bacon, C. R., Nye, C.,
Stelling, P. & Dreher, S. (2004). Chemical versus temporal controls
on the evolution of tholeiitic and calc-alkaline magmas at two
volcanoes in the Alaska^Aleutian arc. Journal of Petrology 45,
203^219.
Gill, J. B. (1981). Orogenic Andesites and Plate Tectonics. Berlin: SpringerVerlag, 390 pp.
Grove, T. L. & Baker, M. B. (1984). Phase equilibrium controls on the
tholeiitic versus calc-alkaline differentiation trends. Journal of
Geophysical Research 89, 3253^3274.
Grove, T. L., Parman, S. W., Bowring, S. A., Price, R. C. & Baker, M.
B. (2002). The role of an H2O-rich fluid component in the generation of primitive basaltic andesites and andesites from the
Mt. Shasta region, N California. Contributions to Mineralogy and
Petrology 142, 375^396.
Hamada, M. & Fujii, T. (2008). Experimental constraints on the
effects of pressure and H2O on the fractional crystallization of
high-Mg island arc basalt. Contributions to Mineralogy and Petrology
155, 767^790.
NUMBER 11
NOVEMBER 2008
Hauri, E. H., Shaw, A. M., Wang, J., Dixon, J. E., King, P. L. &
Mandeville, C. (2006). Matrix effects in hydrogen isotope analysis
of silicate glasses by SIMS. Chemical Geology 235, 352^365.
Hess, P. C. (1989). Origin of Igneous Rocks. Cambridge, MA: Harvard
University Press.
Hildreth, W. (1981). Gradients in silicic magma chambers: implications
for lithospheric magmatism. Journal of Geophysical Research 86,
10153^10192.
Hildreth, W. & Moorbath, S. (1988). Crustal contributions to arc
magmatism in the Andes of central Chile. Contributions to
Mineralogy and Petrology 98, 455^489.
Hirata, T. (1996). Lead isotopic analysis using multiple collectorinductively coupled plasma mass spectrometry coupled with
modified external correction method for mass discrimination
effect. Analyst 121, 1407^1411.
Horn, I., Rudnick, R. L. & McDonough, W. F. (2000). Precise elemental and isotopic ratio determination by simultaneous solution
nebulization and laser ablation-ICP-MS: application to U^Pb
geochronology. Chemical Geology 164, 281^301.
Horwitz, E. P., Chiarizia, R. & Deitz, M. L. (1992). A novel strontium-selective extraction chromatographic resin. Solvent Extraction
and Ion Exchange 10, 313^336.
Hunter, A. G. (1998). Intracrustal controls on the coexistence of
tholeiitic and calc-alkaline magma series at Aso Volcano, SW
Japan. Journal of Petrology 39, 1255^1284.
Hunter, A. G. & Blake, S. (1995). Petrogenetic evolution of a transitional tholeiitic^calc-alkaline series: Towada volcano, Japan.
Journal of Petrology 36, 1579^1605.
Iizuka, T. & Hirata, T. (2005). Improvements of precision and accuracy in in situ Hf isotope microanalysis of zircon using the laser
ablation-MC-ICPMS technique. Chemical Geology 220, 121^137.
Irvine, T. N. & Baragar, W. R. A. (1971). A guide to the chemical classification of the common volcanic rocks. Canadian Journal of Earth
Sciences 8, 523^548.
Jakes, P. & Gill, J. (1970). Rare earth elements and the island arc
tholeiitic series. Earth and Planetary Science Letters 9, 17^28.
Kagami, H. (2005). Formative periods and source materials of
Cretaceous^Paleogene granitoids from Honshu Arc. Journal of the
Geological Society of Japan 111, 441^457.
Kamei, A., Takagi, T. & Kubo, K. (2003). Geology and petrography of
the Abukuma Granites in the Hiyama District, Fukushima Prefecture, NE Japan. Bulletin of the Geological Survey of Japan 54, 395^409.
Katsumata, A. & Kamaya, N. (2003). Low-frequency continuous
tremor around the Moho discontinuity away from volcanoes in
the southwest Japan. Journal of Geophysical Research 30(1), 1020, doi:
10.1029/2002GL015981.
Kawano, Y., Yagi, K. & Aoki, K. (1961). Petrography and petrochemistry of the volcanic rocks of Quaternary volcanoes of northeastern
Japan. Science Reports, Tohoku University, Series 3: Mineralogy, Petrology,
and Economic Geology 7, 1^46.
Kay, R. W. (1978). Aleutian magnesian andesites: melts from subducted Pacific Ocean crust. Journal of Volcanology and Geothermal
Research 4, 117^132.
Kelemen, P. B. (1995). Genesis of high Mg-number andesites and the
continental crust. Contributions to Mineralogy and Petrology 120, 1^19.
Kersting, A. B., Arculus, R. J. & Gust, D. A. (1996). Lithospheric
contributions to arc magmatism: isotope variations along strike in
volcanoes of Honshu, Japan. Science 272, 1464^1468.
Kimura, J.-I. & Yoshida, T. (2006). Contributions of slab fluid, mantle
wedge and crust to the origin of Quaternary lavas in the NE Japan
arc. Journal of Petrology 47, 2185^2232.
Kimura, J.-I., Danhara, T. & Iwano, H. (2000). A preliminary report
on trace element distributions in zircon and apatite crystals using
2006
TATSUMI et al.
CALC-ALKALIC VS THOLEIITIC SERIES
excimer laser ablation inductively coupled plasma mass spectrometry (ExLA-ICP-MS). FissionTrack Newsletter 13, 11^20.
Kimura, J.-I., Yoshida, T. & Iizumi, S. (2002). Origin of low-K
andesite lavas at Nekoma volcano, NE Honshu arc, Japan: geochemical constraints for lower-crustal melts. Journal of Petrology 43,
631^661.
Kogiso, T., Tatsumi, Y. & Nakano, S. (1997). Trace element transport
during dehydration processes in the subducted oceanic crust: 1,
Experiments and implications for the origin of ocean island
basalts. Earth and Planetary Science Letters 148, 193^205.
Kuno, H. (1950). Petrology of Hakone volcano and the adjacent areas,
Japan. Geological Society of America Bulletin 61, 957^1014.
Kuno, H. (1959). Origin of Cenozoic petrographic provinces of Japan
and surrounding areas. Bulletin of Volcanology 20, 37^76.
Kuno, H. (1960). High-alumina basalt. Journal of Petrology 1, 121^145.
Kuno, H. (1968). Differentiation of basalt magmas. In: Hess, H. H. &
Poldervaart, A. (eds) BasaltsçThe Poldervaart Treatise on Rocks of
Basaltic Composition,Volume 2. New York: Interscience, pp. 623^688.
Marsh, B. D. (1989). Magma chambers. Annual Review of Earth and
Planetary Sciences 17, 439^474.
Martin, E. (1986). Effect of steeper Archean geothermal gradient on
geochemistry of subduction-zone magmas. Geology 14, 753^756.
Masuda, Y. & Aoki, K. (1979). Trace element variations in the volcanic
rocks from the Nasu Zone, Northeast Japan. Earth and Planetary
Science Letters 44, 139^149.
Matsui, Y., Onuma, N., Nagasawa, H., Higuchi, H. & Banno, S.
(1977). Crystal structure control in trace element partition between
crystal and magma. Bulletin de la Socie¤ te¤ Franc aise de Mine¤ ralogie et de
Cristallographie 100, 315^324.
Miyashiro, A. (1974). Volcanic rock series in island arcs and active
continental margins. AmericanJournal of Science 274, 321^355.
Miyazaki, T., Takahashi, T., Hirahara, Y., Chang, Q., Vaglarov, B. S.,
Suzuki, K. & Tatsumi, Y. (2007). Low blank and precise Sr^Nd^Pb
analysis and in-situ Sr analysis of minerals in volcanic rocks. Frontier
Research on Earth Evolution 3, 72^76.
Mu«ntener, O., Kelemen, P. B. & Grove, T. L. (2001). The role of H2O
during crystallization of primitive arc magmas under uppermost
mantle conditions and genesis of igneous pyroxenites: an experimental study. Contributions to Mineralogy and Petrology 141, 643^658.
Nakamura, K. & Chang, Q. (2007). Precise determination of ultralow (sub-ng g^1) level rare earth elements in ultramafic rocks by
quadrupole ICP-MS. Geostandards and Geoanalytical Research 31,
185^197.
Nockolds, S. R. & Allen, R. (1953). The geochemistry of some igneous
rock series [Part 1]. Geochimica et Cosmochimica Acta 4, 105^142.
Notsu, K. (1983). Strontium isotope composition in volcanic rocks
from the Northeast Japan arc. Journal of Volcanology and Geothemal
Research 18, 531^548.
Oba, Y. & Konda, T. (1989). Geology and petrology of central Zao
Volcano, Yamagata Prefecture. Bulletin of the Yamagata University,
Natural Science 12, 199^210.
Obara, K. (2002). Nonvolcanic deep tremor associated with subduction in Southwest Japan. Science 296, 1679^1681.
Pachett, P. J., Kuovo, O., Hedge, C. E. & Tatsumoto, M. (1981).
Evolution of continental crust and mantle heterogeneity; evidence from Hf isotopes. Contributions to Mineralogy and Petrology 78,
279^297.
Parman, S. W. & Grove, T. L. (2004). Harzburgite melting with and
without H2O: Experimental data and predictive modeling. Journal
of Geophysical Research 109, paper number B02201.
Patin‹o Douce, A. E. & Beard, J. S. (1995). Dehydration-melting of
biotite gneiss and quartz amphibolite from 3 to 15 kbar. Journal of
Petrology 36, 707^738.
Pearce, J. A, Van der Laan, S. R., Arculus, R. J., Murton, B. J.,
Ishii, T., Peate, D. W. & Parkinson, I. J. (1992). Boninite and harzburgite from Leg 125. (Bonin^Mariana forearc): a case study of
magma genesis during the initial stages of subduction. In:
Fryer, P., Pearce, J. A. & Stokking, L. B. et al. (eds) Proceedings of
the Ocean Drilling Program, Scientific Results, 125. College Station, TX:
Ocean Drilling Program, pp. 623^659.
Pearce, N. J. G., Perkins, W. T., Westgate, J. A., Gorton, M. G.,
Jackson, S. E., Near, C. R. & Chenery, S. P. (1996). A compilation
of new and published major and trace element data for NIST
SRM610 and NIST SRM612 glass reference materials. Geostandard
Newsletters 21, 115^144.
Petford, N. & Atherton, M. (1996). Na-rich partial melts from newly
underplated basaltic crust: the Cordillera Blanca batholith, Peru.
Journal of Petrology 37, 1491^1521.
Petford, N. & Gallagher, K. (2001). Partial melting of mafic (amphibolitic) lower crust by periodic influx of basaltic magma. Earth and
Planetary Science Letters 193, 483^499.
Pichavant, M., Martel, C., Bourdier, J. L. & Scaillet, B. (2002).
Physical conditions, structure, and dynamics of a zoned magma
chamber: Mont Pele¤e (Martinique, Lesser Antilles Arc). Journal of
Geophysical Research 107, 2093.
Plank, T. & Langmuir, C. H. (1988). An evaluation of the global variations in the major element chemistry of arc basalts. Earth and
Planetary Science Letters 90, 349^370.
Prouteau, G. & Scaillet, B. (2003). Experimental constraints on
the origin of the 1991 Pinatubo dacite. Journal of Petrology 44,
2203^2241.
Raia, F. & Spera, F. J. (1997). Simulations of crustal anatexis; implications for the growth and differentiation of continental crust.
Journal of Geophysical Research 102, 22629^22648.
Rudnick, R. L. (1995). Making continental crust. Nature 378, 571^578.
Rudnick, R. L. & Gao, S. (2003). The composition of the continental
crust. In: Rudnick, R. L. (ed.) The Crust. Oxford: Pergamon,
pp. 1^64.
Sakayori, A. (1992). Geology and petrology of Zao Volcano. Journal of
Mineralogy, Petrology and Economic Geology 87, 433^444.
Sakuyama, M. (1979). Lateral variation of H2O contents in
Quaternary magmas of north-eastern Japan. Earth and Planetary
Science Letters 43, 103^111.
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. (1983). Phenocryst assemblages and H2O content in
Circum-Pacific arc magmas. In: Hilde, T. W. C. & Uyeda, S. (eds)
Geodynamics of the Western Pacific^Indonesian Region. American
Geophysical Union, Geodynamics Series11, 143^158.
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. (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.
Shukuno, H. (2003). Quantitative analysis of rock-forming minerals
and volcanic glasses by electron probe microanalyzer. Frontier
Research on Earth Evolution 1, 129^136.
Sisson, T. W. & Grove, T. L. (1993). Experimental investigations of
the role of H2O in calc-alkaline differentiation and subduction
zone magmatism. Contributions to Mineralogy and Petrology 113,
143^166.
2007
JOURNAL OF PETROLOGY
VOLUME 49
Smith, D. R. & Leeman, W. P. (1987). Petrogenesis of Mount
St. Helens dacitic magmas. Journal of Geophysical Research 92,
10313^10334.
Stern, C. R. & Kikian, R. (1986). Role of the subducted slab, mantle
wedge and continental crust in the generation of adakites from the
Andean Austral volcanic zone. Contributions to Mineralogy and
Petrology 123, 263^281.
Streck, M. J., Leeman, W. P. & Chesley, J. (2007). High-magnesian
andesite from Mount Shasta; a product of magma mixing and contamination, not a primitive mantle melt. Geology 35, 351^354.
Sun, S. S. & McDonough, W. F. (1989). Chemical and isotopic
systematics of oceanic basalts: implications for mantle composition
and processes. In: Saunders, A. D. & Norry, M. J. (eds) Magmatism
in the Ocean Basins. Geological Society, London, Special Publications 42,
313^345.
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.
Takahashi, T., Yoshikawa, M., Shibata, T., Tatsumi, Y. & Shimizu, N.
(2005). Sr isotopic micro analyses of plagioclase in andesites from
Zao Volcano, NE Japan. Frontier Research on Earth Evolution 2, Sam-18.
Takaoka, N., Konno, K., Oba, Y. & Konda, T. (1989). K^Ar dating of
lavas from Zao Volcano, north-eastern Japan. Journal of the Geological
Society of Japan 95, 157^170.
Tani, K., Kawabata, H., Chang, Q., Sato, K. & Tatsumi, Y. (2006).
Quantitative analyses of silicate rock major and trace elements
by X-ray fluorescence spectrometer: Evaluation of analytical
precision and sample preparation. Frontier Research on Earth
Evolution 2, 8.
Tatsumi, Y. (1981). Melting experiments on a high-magnesian andesite.
Earth and Planetary Science Letters 54, 357^365.
Tatsumi, Y. (2006). High-Mg andesites in the Setouchi volcanic belt,
southwestern Japan: analogy to Archean magmatism and continental crust formation? Annual Review of Earth and Planetary Sciences 34,
467^499.
Tatsumi, Y. & Eggins, S. (1995). Subduction Zone Magmatism. Boston,
MA: Blackwell Science.
Tatsumi, Y. & Hanyu, T. (2003). Geochemical modeling of dehydration and partial melting of subducting lithosphere: Towards a comprehensive understanding of high-Mg andesite formation in the
Setouchi volcanic belt, SW Japan. Geochemistry, Geophysics, Geosystems
4, 1081, doi:10.1029/2003GC000530.
Tatsumi, Y. & Kogiso, T. (2003). The subduction factory: its role in
the evolution of the Earth’s crust and mantle. In: Larter, R. D. &
Leat, P. T. (eds). Intra-Oceanic Subduction Systems:Tectonic and Magmetic
Processes. Geological Society, London, Special Publications 219, 55^80.
Tatsumi, Y., Sakuyama, M., Fukuyama, H. & Kushiro, I. (1983).
Generation of arc basalt magmas and thermal structure of the
mantle wedge in subduction zones. Journal of Geophysical Research B
88, 5815^5825.
Tatsumi, Y., Kawabata, H., Sato, K., Miyazaki, T., Chang, Q.,
Takahashi, T., Tani, K., Shibata, T. & Yoshikawa, M. (2006). The
petrology and geochemistry of Oto-Zan composite lava flow on
Shodo-Shima Island, SW Japan: remelting of a solidified high-Mg
andesite magma. Journal of Petrology 47, 595^629.
NUMBER 11
NOVEMBER 2008
Tatsumi, Y., Shukuno, H., Tani, K., Takahashi, N., Kodaira, S. &
Kogiso, T. (2008). Structure and growth of the Izu^Bonin^
Mariana arc crust: II. The role of crust^mantle transformation
and the transparent Moho in arc crust evolution. Journal of
Geophysical Research 113, B02203, doi:10.1029/2007JB005121.
Taylor, S. R. (1995). Cosmochemistry; potassium tells a tale. Nature
376, 20^21.
Tepley, F. J., Davidson, J. P. & Clynne, M. A. (1999). Magmatic interactions as recorded in plagioclase phenocrysts of Chaos Crags,
Lassen Volcanic Center, California. Journal of Petrology 40, 787^806.
Villiger, S., Ulmer, P. & Mu«ntener, O. (2007). Equilibrium and fractional crystallization experiments at 07 GPa; the effect of pressure
on phase relations and liquid compositions of tholeiitic magmas.
Journal of Petrology 48, 159^184.
Vroon, P. Z., van der Wagt, B., Koornneef, J. M. & Davies, G. R.
(2008). Problems in obtaining precise and accurate Sr isotope analysis from geological materials using laser ablation MC-ICPMS.
Analytical and Bioanalytical Chemistry 390, 465^476, doi:10.1007/
s00216-007-1742-9.
Wada, K. (1981). Contrasted petrological relations between tholeiitic
and calc-alkalic series from Funagata Volcano, Northeast Japan.
Journal of the Japanese Association of Mineralogists, Petrologists and
Economic Geologists 76, 215^231.
Wada, K. (1985). Magma mixing process of calc-alkalic andesites from
Funagata Volcano. Journal of the Japanese Association of Mineralogists,
Petrologists and Economic Geologists 80, 467^483.
Wager, L. R. & Deer, W. A. (1939). Geological investigations in East
Greenland. Part III. The petrology of the Skaergaard intrusion,
Kangerdlugssuaq, East Greenland. Meddelelser om Gronland 105, 352.
Wells, P. R. A. (1977). Pyroxene thermometry in simple and complex
systems. Contributions to Mineralogy and Petrology 62, 129^139.
Wilson, M. (1989). Igneous Petrogenesis. London: Unwin Hyman.
Woodhead, J. & Hergt, J. (2001). Strontium, neodymium and lead isotope analyses of NIST glass certified reference materials: SRM 610,
612, 614. Geostandards Newsletter 25, 261^266.
Woodhead, J., Hergt, J., Shelley, M., Eggins, S. & Kemp, R. (2004).
Zircon Hf-isotope analysis with an excimer laser, depth profiling,
ablation of complex geometries, and concomitant age estimation.
Chemical Geology 209, 121^135.
Woodhead, J., Swearer, S., Hergt, J. & Maas, R. (2005). In situ
Sr-isotope analysis of carbonates by LA-MC-ICP-MS: interference
corrections, high spatial resolution and an example from otolith
studies. Journal of Analytical Atomic Spectrometry 20, 22^27.
Woodhead, J. D., Hergt, J. M., Davidson, J. P. & Eggins, S. M. (2001).
Hafnium isotope evidence for ‘conservative’ element mobility
during subduction zone processes. Earth and Planetary Science Letters
192, 331^346.
Yoder, H. S. & Tilley, C. E. (2003). Origin of basalt magmas; an
experimental study of natural and synthetic rock systems. Journal
of Petrology 3, 342^529.
Yogodzinski, G. M., Volynets, O. N., Koloskov, A.V., Seliverstov, N. I.
& Matvenkov, V. V. (1994). Magnesian andesites and the subduction
component in a strongly calc-alkaline series at Piip Volcano, far
western Aleutians. Journal of Petrology 35, 163^204.
Yoshida, T. & Aoki, K.-i. (1984). Geochemistry of major and trace elements in the Quaternary volcanic rocks from northeast Honshu,
Japan. Science Reports, Tohoku University, Series 3: Mineralogy, Petrology,
and Economic Geology 16, 1^34.
2008