Evolution of Deeper Basaltic and Shallower Andesitic Magmas during the AD 1469-1983 Eruptions of Miyake-Jima Volcano, Izu-Mariana Arc: ...

Page created by Ellen Olson
 
CONTINUE READING
JOURNAL OF PETROLOGY                   VOLUME 44        NUMBER 12        PAGES 2113±2138         2003      DOI: 10.1093/petrology/egg072

Evolution of Deeper Basaltic and Shallower
Andesitic Magmas during the AD 1469---1983
Eruptions of Miyake-Jima Volcano, Izu---
Mariana Arc: Inferences from Temporal
Variations of Mineral Compositions in

                                                                                                                                             Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
Crystal-Clots

MIZUHO AMMA-MIYASAKA1* AND MITSUHIRO NAKAGAWA2
1
HOKKAIDO BRANCH, GEOLOGICAL SURVEY OF JAPAN, AIST, N8W2, KITA-KU, SAPPORO 060-0808, JAPAN
2
DEPARTMENT OF EARTH AND PLANETARY SCIENCES, GRADUATE SCHOOL OF SCIENCE, HOKKAIDO
UNIVERSITY, SAPPORO 060-0810, JAPAN

                          RECEIVED SEPTEMBER 2, 2002; ACCEPTED MAY 29, 2003

Miyake-jima volcano has erupted at least 13 times during the        end-member magmas suggest that the basaltic magma has
period 1469---1983. To understand the historic magmatic pro-        differentiated gradually since 1469, and that its magmatic
cesses, we focus on the mineral assemblage and chemical com-        temperature has fallen from 1220 to 1180 C. Conversely, the
positions of crystal-clots in single samples from each of the       andesitic magma has changed in a complex fashion to become
eruptions. Most of the historic lavas consist of nearly aphyric     more mafic (the magmatic temperature rose from 1050 to
to weakly porphyritic basalt to andesite, but there also exist      1100 C). As a result of this study, it is estimated that the
megacryst-bearing rocks. The megacrysts are considered to be        basaltic magma after the 1983 eruption was the least mafic, and
xenocrysts from a deep-seated plutonic body. Many samples of        the andesitic magma the most mafic, of the historic eruptions.
each eruption contain two types of clots beside megacrysts,
termed here B-type and A-type. The B-type clots are composed
of olivine, clinopyroxene and plagioclase, whereas the A-type       KEY WORDS: andesite; basalt; crystal-clots; evolution of magma;

clots additionally contain magnetite and orthopyroxene. Compo-      Miyake-jima volcano; magma mixing
sitional relationships between these mafic minerals suggest that
the minerals in the same type of clots are in equilibrium.
Comparing the chemical compositions of the minerals in the          INTRODUCTION
two types of clots in each sample, they are derived from distinct   The present state of the magma plumbing system
magmas: the B-type clots from basaltic magma and the A-type         beneath an active volcano provides important clues
clots from andesitic magma. During the historic activity, the       about the nature of the volcanic activity and the poten-
magma plumbing system appears to have included two magma            tial for future eruptions. Temporal changes in the
storage systems: a deep-seated basaltic and a shallower andesitic   magma plumbing system can be evaluated through
one. In many cases, basaltic magma has injected into shallower      detailed petrological studies of stratigraphically well-
andesitic magma to form mixed magma; however, andesitic             constrained eruption sequences. Many studies have
magma has sometimes erupted alone without extensive injections      focused on the magma plumbing system and eruption
of basaltic magma. Temporal variations of mineral composi-          processes based on the eruptive sequence of a single
tions in the clots and estimated whole-rock compositions of the     eruption (e.g. Reagan et al., 1987; Wolfe et al., 1987;

*Corresponding author. Telephone: ‡81-11-709-1813.          Fax:    Journal of Petrology 44(12) # Oxford University Press 2003; all rights
‡81-11-709-1817. E-mail: m.miyasaka@aist.go.jp                      reserved
JOURNAL OF PETROLOGY                 VOLUME 44       NUMBER 12       DECEMBER 2003

Garcia et al., 1989, 1992, 1996, 2000; Wolf &                  focused on the historic 1940 and 1962 eruptions and
Eichelberger, 1997; Marianelli et al., 1999; Nakagawa          established that deeper basaltic magma and shallower
et al., 1999, 2002; Streck et al., 2002). A series of papers   andesitic magma existed beneath the volcano during
on the Pu'u O'o eruption of Kilauea (Garcia et al.,            this period. It is not clear, however, how long the
1989, 1992, 1996, 2000) has revealed changing mag-             magma plumbing system has existed and how it has
matic processes over a period of 17 years; Reagan et al.       evolved through time.
(1987) and Streck et al. (2002) discussed the role of             In this paper, we describe the evolution of the whole-
basalt replenishment in the long-lived (about 30 years)        rock and mineral chemistry of the historic lavas from
eruption of Arenal volcano, Costa Rica. Clearly, to            1469 to 1983. We focus on the mineral assemblage and
study the evolution of magma plumbing systems, it is           chemical compositions of crystal-clots in single samples
important to investigate as long a period of activity          to clarify whether minerals in each sample are in equi-
as possible [e.g. Wright & Fiske (1971) for Kilauea,           librium or not. We demonstrate the existence of two

                                                                                                                            Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
Hawaii; Borgia et al. (1988) for Arenal, Costa Rica;           magma storage systems since 1469, one basaltic and
Fichaut et al. (1989) for Mt. Pelee, Martinique; Belkin        one andesitic. We carefully investigate the evolution
et al. (1993), Villemant et al. (1993) and Tedesco et al.      of both magma types for the past 500 years and the
(1998) for Vesuvius, Italy; D'Antonio et al. (1999) for        interaction between these magmas during each erup-
Campi Flegrei caldera, Italy; Nakano & Yamamoto                tion episode. This provides important new insights into
(1991) for Izu---Oshima, Japan]. In most of these stu-         the evolution of the historic magma plumbing system,
dies, however, the evolution of the magma plumbing             and allows us to estimate the possible state of the system
systems was discussed mainly in terms of the temporal          just before the 2000 eruption.
change in whole-rock compositions. Although detailed
analysis of phenocryst minerals is useful to understand
magmatic processes (e.g. Nakamura, 1995; Umino &               GEOLOGY AND 1469---1983
Horio, 1998; Nakagawa et al., 1999, 2002; Streck et al.,
2002), a systematic mineralogical study of the eruptive
                                                               ERUPTIONS OF MIYAKE-JIMA
activity of an active volcano for a considerable period        VOLCANO
(for hundreds to thousands years), based on the evolu-         Miyake-jima volcano is a composite volcano with two
tion of both whole-rock and mineral chemistry, has not         nested calderas (Fig. 1), and is composed of tholeiitic
been carried out, except for that by Borgia et al. (1988).     basalt and andesite. Based on the eruption style and
   Miyake-jima volcano has erupted at least 13 times           petrological characteristics of lavas and scoria, Tsukui
since AD 1469. During the period 1469---1983, magma            et al. (2001) divided the volcanic activity of the last
effused mainly from flank fissures. The most recent            10 000 years into four stages: 10 000---7000, 4000---2500,
eruption, which began in June 2000 (Nakada et al.,             2500 years BP to AD 1154, and since AD 1469. During
2001), however, is distinctive compared with other             the first stage, the main cone was constructed, and the
historic eruptions. Earthquake swarms occurred begin-          outer caldera filled. Erupted materials are porphyritic
ning the night of June 26, and a submarine eruption            basalts. The second stage began after a 3000 year
took place on the morning of June 27. A summit erup-           repose period, consisting of andesitic lavas and scoria
tion (ash plume) followed on July 8, and the summit            erupted from lateral and central vents. The third stage
area suddenly subsided. The collapse continued                 began with a large-scale eruption that formed the inner
`silently' until mid-August, resulting in the formation        caldera, which has subsequently been filled with nearly
of a new caldera with a volume of about 06 km3 . The          aphyric to weakly porphyritic basalts.
volcano is now (April, 2003) discharging a large quan-            The latest stage, since1469, began after a 300 year
tity of volcanic gases (SO2, etc.) and inhabitants of          repose period. The eruptions were characterized by
Miyake-jima have been evacuated since September                the effusion of lavas and scoria mainly from NE---SW-
2000. Several petrological studies have been carried           trending flank fissures (Fig. 1). Eruptions occurred
out on recent historic eruptions, except the 2000 erup-        every 50---70 years before 1811, and have become
tion, of Miyake-jima volcano (e.g. Iwasaki et al., 1982;       more frequent (every 20---70 years) since then
Fujii et al., 1984; Soya et al., 1984; Sato et al., 1996;      (Table 1). Most of the eruptions continued for a short
Amma-Miyasaka & Nakagawa, 1998, 2002). Sato et al.             duration (typically a day to a month) whereas the
(1996) dealt with most of the historic lavas since 1643.       1763 eruption was prolonged, lasting for 6 years.
They showed that whole-rock compositions became                The volume of magma erupted is estimated to be
more silicic from 1643 to 1874, and then more mafic            50066 km3 for each eruption. Although these values
from 1940 until 1983. They proposed that mafic                 may be underestimated, especially for the older
magma was recently injected into the magma plumb-              eruptives, there seems to be little correlation between
ing system. Amma-Miyasaka & Nakagawa (2002)                    the eruptive mass and the duration of each eruption.

                                                           2114
AMMA-MIYASAKA AND NAKAGAWA                           MAGMA EVOLUTION, MIYAKE-JIMA

                                                                                                                                                      Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
Fig. 1. Location of Miyake-jima volcano and distribution of lavas and pyroclastic deposits erupted during the period 1469---1983. These
eruptions were characterized by the effusion of lavas and scoria mainly from NE---SW-trending flank fissures. The 1643 lava flows are
composed of three lobes: main-N, main-S and central. Modified from Isshiki (1984).

                    Table 1: Summary of eruptive history of Miyake-jima volcano from 1469 to 1983

Age (date     AD)   Vent location   Duration      Volume* (km3 DRE)   Rocks                                 Total phenocryst (vol. %)   SiO2 (wt %)

1983                SW flank        15 h           0.007              nearly aphyric---weakly porphyritic   1.0---8.5                   52.0---54.4
1962                NE flank        30 h           0.006              nearly aphyric                        0.2---1.5                   53.1---55.5
1940                summit          25 days        0.015              nearly aphyric---weakly porphyritic   2.8---8.0                   54.2---54.7
                    NE flank        23 h                              nearly aphyric                        0.7---2.6                   54.1---56.4
                                                                      megacryst bearing                     1.4---27.3                  51.4---55.5
1874                N flank         4---5 days     0.01               nearly aphyric                        1.6---3.0                   54.5---55.3
                                                                      megacryst bearing                     3.7---15.7                  52.8---55.2
1835                W flank         10 days       50.001              nearly aphyric---weakly porphyritic   1.3---7.9                   52.1---55.5
1811?               NE flank            7 days    50.001              megacryst bearing                     3.1---10.7                  54.6---55.2
1763---1769         SW flank            6 years    0.066              nearly aphyric                        0.6---1.4                   53.0---54.5
1712                SW flank        14 days        0.001              nearly aphyric                        1.3---3.2                   50.6---52.1
1643                SW flank        21 days        0.012              nearly aphyric                        0.3---1.3                   50.8---53.7
1595                SE flank?       ?             50.001              nearly aphyric                        1.2---1.7                   51.8---52.6
1535                SE flank?       ?              0.003              nearly aphyric                        0.9---1.7                   51.8---52.6
1469                W flank?        ?              0.002              nearly aphyric                        1.4---2.7                   53.2---54.9

*Quoted from Tsukui & Suzuki (1998) except for the 1811 lavas.

                                                                      2115
JOURNAL OF PETROLOGY                 VOLUME 44       NUMBER 12       DECEMBER 2003

ANALYTICAL METHODS                                             1835 lavas, however, do not contain olivine pheno-
                                                               crysts, whereas the 1643 and 1712 lavas lack ortho-
For each of the 1469---1983 eruptions, we collected
                                                               pyroxene. In addition, in some of the 1940 and 1962
samples to represent all source vents as well as the
                                                               lavas both olivine and orthopyroxene phenocrysts are
complete eruptive sequence. Mineral core composi-
                                                               absent (Amma-Miyasaka & Nakagawa, 1998).
tions were determined by a single analysis from the
                                                                  Plagioclase is the dominant phenocryst phase and its
centre of the minerals for representative samples that
                                                               maximum size is 46 mm in length. The An content of
reflect the whole-rock variations of each eruption. The
                                                               the plagioclase phenocrysts is in the range 53---96 mol %
JEOL 733 and 8800 electron probe microanalysers
                                                               (Fig. 2). Histograms for the 1835, 1962 and 1983
at Hokkaido University were used for the mineral
                                                               samples show unimodal distributions with a peak of
analyses. Operating conditions were 15 kV accelerat-
                                                               An ˆ 80---90, whereas the other samples show bimodal
ing voltage and 20 nA beam current with a minimum
                                                               distributions. Although the peaks of the An contents

                                                                                                                           Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
spot size of 1 mm. Each element was counted for 30 s on
                                                               are variable in each eruption, most of them are in the
the peak and 20 s on the background. Corrections were
                                                               range of An ˆ 80---90 and An ˆ 60---75.
made according to the ZAF method. Whole-rock com-
                                                                  Clinopyroxene less than 24 mm in length is the
positions were determined by X-ray fluorescence
                                                               dominant mafic phenocryst. Most of the phenocrysts
(XRF), using a Philips PW-1404 system with a Rh
                                                               are augite, although pigeonite sometimes occurs in
tube at Hokkaido University. Glass beads were used
                                                               samples of the 1712 and 1940 lavas (Fig. 3). The Mg-
for major element analysis, and pressed pellets for trace
                                                               number of most of the clinopyroxene phenocrysts is in
element analysis. Major element compositions were
                                                               the range 60---79; however, one phenocryst in the 1962
determined for 230 samples, and trace element compo-
                                                               lava is more mafic (Mg-number 480). There is little
sitions for 198 samples from the 1469---1983 eruptive
                                                               difference in the compositions of the phenocrysts from
products. Whole-rock compositions of all the samples
                                                               samples of different eruption ages, although composi-
for which we have analysed mineral compositions are
                                                               tional variations in the 1835 and 1962 lavas seem to be
listed in Table 2 and representative mineral composi-
                                                               smaller compared with other eruptives (Fig. 3).
tions are given in the Appendix. The complete whole-
                                                                  Olivine phenocrysts are usually smaller than 20 mm
rock and mineral composition dataset is included in
                                                               in length. The Fo content of the olivine phenocrysts is
Electronic Appendices 1---5, which may be downloaded
                                                               in the range 52---78 mol % (Fig. 3). The variations in
from the Journal of Petrology web site at http://www.
                                                               the phenocryst composition in the samples older than
petrology.oupjournals.org/.
                                                               the 1763 eruption are greater than those in the samples
                                                               younger than the 1874 eruption.
PETROGRAPHY AND MINERAL                                           Orthopyroxene and magnetite phenocrysts are rare
                                                               and smaller than 16 mm and 04 mm in length, respect-
CHEMISTRY                                                      ively. The Mg-number of orthopyroxene phenocrysts
Although most of the 1469---1983 lavas are nearly              is in the range 52---75 (Fig. 3). Although most of them
aphyric to weakly porphyritic, with less than 9 vol. %         are in the range Mg-number ˆ 60---75, some of the
of phenocrysts, there also exist rocks characterized by        orthopyroxene phenocrysts in the 1962 and 1983
the presence of anorthite (up to 3 cm in length) and           lavas have Mg-number 560.
olivine megacrysts (Table 1). Phenocryst minerals are
plagioclase, olivine, clinopyroxene, orthopyroxene and
magnetite. The groundmass is composed of plagioclase,          Megacryst-bearing rocks
clinopyroxene, magnetite and brown glass, and ranges           Megacryst-bearing rocks occur only in the 1811, 1874
in texture from intersertal to hyalo-ophitic.                  and 1940 flank eruption sequences (Table 1), and are
                                                               orthopyroxene-bearing clinopyroxene---olivine basalt
                                                               to andesite (Figs 2 and 3). The total volume of pheno-
Megacryst-free rocks (nearly aphyric to                        crysts varies from 1 to 27 vol. %, and gradually
weakly porphyritic rocks)                                      increases from the 1811 (511 vol. %), to 1874
All of the dated lavas except for the 1811 lavas include       (516 vol. %) and to 1940 (527 vol. %) eruptions.
nearly aphyric rocks (Table 1), and their modal                Petrographic characteristics and mineral compositions
volumes of phenocrysts are less than 3 vol. %. Some of         of the megacrysts in the 1940 samples were described
the 1835 and 1940 summit and 1983 lavas are weakly             by Amma-Miyasaka & Nakagawa (2002). Although
porphyritic, containing 3---9 vol. % phenocrysts. We           the mineral assemblage and chemical compositions of
generally call these `megacryst-free rocks'. Mafic pheno-      the megacrysts in the 1874 samples are similar to those
crysts in the megacryst-free rocks usually consist of clino-   of the 1940 eruption, clinopyroxene megacrysts are also
pyroxene, orthopyroxene and olivine (Figs 2 and 3). The        found in the 1811 megacryst-bearing rocks in addition

                                                           2116
Table 2: Representative whole-rock compositions during the period 1469---1983

                                                                                                                                                                    Downloaded
                                                                                                                                                                     AMMA-MIYASAKA
       Age:          1983      1983       1962      1962      1962       1962      1962       1962      1940      1940      1940      1940      1940      1940
       Sample no.:   1522      2203       1802      1810      122508     2233      122005     1901      2003      2008      1803      1820*     122002    1815

                                                                                                                                                                               from https://academic.oup.com/petrology/article/44/12/2113/1460172
       wt %
       SiO2            52.69     52.87     53.24      52.77     54.62      52.78     53.23      51.67     53.39     53.99     55.67    50.46     53.46      55.60
       TiO2             1.31      1.39      1.36       1.38      1.37       1.40      1.39       1.40      1.32      1.32      1.31     0.97      1.39       1.27
       Al2O3           14.80     15.39     14.46      14.46     14.53      14.65     14.60      14.49     14.78     14.90     14.87    16.26     14.56      14.72

                                                                                                                                                                                             AND NAKAGAWA
       Fe2O3           14.21     14.44     14.19      14.55     13.67      14.74     14.25      14.92     13.90     13.71     13.16    12.49     14.21      12.93
       MnO              0.20      0.20       0.20      0.20      0.20       0.20      0.20       0.20      0.20      0.20      0.20      0.19      0.20      0.20
       MgO              4.36      4.00       3.51      3.77      3.52       3.98      3.98       3.97      3.86      3.91      3.54      6.30      3.91      3.31
       CaO              9.13      9.23       8.01      8.34      8.04       8.64      8.64       8.65      8.52      8.47      7.99     10.10      8.46      7.70
       Na2O             2.79      2.77       3.11      3.01      3.06       2.84      2.80       2.88      2.91      2.93      2.97      2.18     2.85       3.25
       K2O              0.58      0.53       0.66      0.62      0.64       0.56      0.55       0.57      0.60      0.59      0.62      0.41     0.61       0.71
       P2O5             0.16      0.15       0.17      0.17      0.17       0.15      0.16       0.16      0.16      0.15      0.15      0.11     0.15       0.18
2117

       Total          100.23    100.97      98.91     99.27     99.80      99.94     99.80      98.91     99.64    100.17    100.48    99.47     99.80      99.87

                                                                                                                                                                                                                      MAGMA EVOLUTION, by
       FeO*/MgO         2.93      3.25       3.64      3.47      3.49       3.33      3.22       3.38      3.24      3.15      3.34     1.78      3.27       3.51
       ppm
       Sc              49        47         48        49        46         52        51         52        46        47        45        36        51        41
       V              386       437        357       389       301        439       418        434       365       360       313       312       411       289
       Cr              21        10          7         7         9          9         9          7        11        12        10        33         9        11
       Ni              12        10          8         8         7         10         9         10        10         8         8        48         9         7
       Cu             119       147        134       142        98        142       142        138       132       116        82       109       139        98
       Rb               8         7          8         8         8          8         8          7         8         8         8         6         8         9
       Zn             117       120        122       123       116        124       119        121       120       117       117        98       125       120

                                                                                                                                                                                                                                                  MIYAKE-JIMA
       Sr             235       236        235       235       234        234       233        234       236       234       235       235       233       234

                                                                                                                                                                                                                                                   guest on 14 January 2021
       Y               34        35         37        36        38         35        35         35        37        35        37        25        35        40
       Zr              72        71         81        77        83         71        70         71        74        74        82        50        74        87
       Ba             190       187        213       203       231        191       186        196       203       204       224       136       201       237
       Pb               4         4          3         4         5          4         5          4         4         3         4         5         8         5
       Ga              16        17         16        16        16         16        16         16        16        16        15        14        15        16
Table 2: continued

       Age:          1940      1874       1874       1874       1874       1874        1874       1835       1835       1811?       1811?     1763      1763      1763
       Sample no.:   1818*     1828*      1841       1842       11008*     11207*      31910*     1918       1919       1806*       1808*     1603      1609      2021

       wt %
       SiO2            50.79     52.14      53.11      54.93      52.72      52.81       52.97      52.38      53.26       54.10      53.20     50.12     51.50     52.97
       TiO2             0.96      1.18       1.35       1.36       1.21       1.24        1.22       1.25       1.29        1.13       1.17      1.29      1.35      1.36
       Al2O3           16.55     16.13      14.56      14.89      15.72      15.62       15.58      14.95      14.92       15.27      14.88     14.90     14.68     14.82
       Fe2O3           12.31     13.08      14.03      13.62      13.22      13.35       13.08      13.85      13.48       12.90      13.38     15.00     14.82     14.35
       MnO              0.19      0.20       0.20       0.20       0.19       0.20        0.23       0.20       0.20        0.19       0.19      0.20      0.20      0.20
       MgO              6.24      4.47       3.60       3.67       4.27       3.96        4.10       4.47       3.45        3.92       3.94      4.88      3.97      3.99
                                                                                                                                                                            JOURNAL OF PETROLOGY

       CaO             10.19      9.57       8.18       8.25       9.35       9.15        9.17       9.21       8.15        8.39       8.26      9.96      8.74      8.53
       Na2O             2.20      2.59       3.06       3.01       2.60       2.76        2.74       2.61       3.06        2.92       2.99      2.45      2.86      2.93
       K2O              0.39      0.50       0.64       0.64       0.53       0.55        0.57       0.57       0.66        0.64       0.65      0.48      0.59      0.59
       P2O5             0.11      0.14       0.17       0.18       0.14       0.15        0.15       0.14       0.18        0.14       0.15      0.12      0.15      0.15

2118
       Total           99.93    100.00      98.90     100.75      99.96      99.78       99.80      99.63      98.65       99.60      98.81     99.40     98.86     99.89
                                                                                                                                                                            VOLUME 44

       FeO*/MgO         1.77      2.63       3.50       3.34       2.78       3.03        2.87       2.79       3.51        2.96       3.05      2.76      3.36      3.23
       ppm
       Sc              41        43         46         48         45         42          41         48         42          42         40        53        49        47
       V              321       332        340        334        349        338         347        368        320         304        309       452       417       393
       Cr              36        20          9         10         20         17          12         28         10          17         15        32        10         7
                                                                                                                                                                            NUMBER 12

       Ni              44        18          7         10         17         11          12         13          8          15         13        15        10         8
       Cu             115       124        130        129        132        130         137         95         97         140        140       134       137       124
       Rb               6         7          8          8          7          8           8          8          9           8          8         6         8         8
       Zn             101       105        122        123        110        110         106        112        116         113        114       112       119       123
       Sr             233       241        238        237        237        244         243        231        242         228        227       230       234       237
       Y               25        31         38         38         33         33          34         34         37          35         34        31        35        36
       Zr              50        63         81         79         65         68          73         72         83          76         77        60        72        74
                                                                                                                                                                            DECEMBER 2003

       Ba             145       182        218        219        179        191         208        186        224         207        212       165       189       213
       Pb               5         5          4          4          1         11           3          5          5           4          4         4         5         6
       Ga              14        17         16         16         15         16          18         15         16          16         16        16        17        17

                                 Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
Age:          1763      1712         1712       1643      1643       1643      1595       1595      1595       1535      1535       1535      1469      1469      1469
       Sample no.:   2023      1504         1936       1620      1627       2222      2110       2113      2501       2116      2117       2502      1921      2227      2315

       wt %
       SiO2            51.75     49.36        51.24      49.64     51.08      51.81     51.24      50.68     51.08      50.64     51.38      51.63     51.78     53.46     54.17
       TiO2             1.36      1.34         1.31       1.31      1.38       1.19      1.14       1.15      1.15       1.14      1.15       1.16      1.17      1.18      1.18
       Al2O3           14.44     14.55        14.93      14.71     14.45      15.12     15.34      15.33     15.39      15.23     15.35      15.47     14.95     15.04     15.09
       Fe2O3           14.73     15.54        14.73      15.30     15.33      13.76     13.68      13.88     13.74      13.86     13.68      13.70     13.61     13.33     13.01
       MnO              0.20      0.21         0.21       0.21      0.20       0.20      0.20       0.20      0.20       0.20      0.20       0.20      0.20      0.20      0.19
       MgO              3.85      4.94         4.98       4.89      4.31       4.80      5.13       5.05      5.09       5.03      5.08       5.14      4.47      4.40      4.18
       CaO              8.48     10.12        10.06      10.05      9.18       9.61     10.02      10.03      9.99      10.06     10.02      10.06      9.25      9.09      8.72
       Na2O             2.93      2.48         2.40       2.51      2.67       2.52      2.40       2.41      2.37       2.43      2.37       2.35      2.59      2.59      2.68
       K2O              0.60      0.49         0.45       0.46      0.53       0.52      0.44       0.45      0.45       0.46      0.45       0.44      0.59      0.59      0.62
       P2O5             0.16      0.13         0.13       0.12      0.14       0.13      0.12       0.12      0.12       0.12      0.12       0.12      0.14      0.14      0.15
       Total           98.50     99.16       100.44      99.20     99.27      99.66     99.71      99.30     99.58      99.17     99.80     100.27     98.75    100.02     99.99
                                                                                                                                                                                   AMMA-MIYASAKA AND NAKAGAWA

       FeO*/MgO         3.44      2.83         2.66       2.81      3.20       2.58      2.40       2.47      2.43       2.48      2.42       2.40      2.74      2.72      2.80

2119
       ppm
       Sc              47        54           52         53        53         44        46         48        47         45        50         45        45        45        43
       V             404       462           455        446       452        362       386        380       383        388       388        373      357       352       340
       Cr               7        31           32         32         8         37        42         42        43         42        45         42        34        30        30
       Ni               9       15            15         14        11         16        15         17        16         16        17         17       14         12       13
       Cu            147       157           163        160       141        118       140        122       126        127       133        128      123       117       103
       Rb               8         6            7          6         7          7         7          6         6          6         6          6         8         8         8
       Zn            121       114           114        112       122        111       106        108       106        109       110        104      109       111       110
       Sr            236       229           229        230       233        229       232        230       227        230       228        229      227       227       226
       Y              35        30            30         31        33         32        29         29        29         30        29         29       34         34       36
       Zr             74        60            59         59        67         65        59         58        58         57        58         58       74         75       80
       Ba            200       164           160        168       178        177       170        161       161        164       160        168      198       203       214
       Pb               4         4            4          4         3          3         5          3         4          3         2          3        5          3        4
       Ga             17        16            16         16        16         16        16         15        15         15        16         16       16        16        15
                                                                                                                                                                                   MAGMA EVOLUTION, MIYAKE-JIMA

       *Megacryst-bearing rocks.

                                      Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
JOURNAL OF PETROLOGY                      VOLUME 44         NUMBER 12          DECEMBER 2003

                                                                                                                                             Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
Fig. 2. Histograms of An contents of plagioclase phenocryst cores in megacryst-free and megacryst-bearing rocks erupted during the period
1469---1983. The megacryst-bearing rocks have more Ca-rich plagioclase compared with the megacryst-free rocks. The 1811 lavas consist only
of the megacryst-bearing rocks.

to olivine and plagioclase. The chemical compositions of              butions of plagioclase with An 490 and An ˆ 80---90
these megacrysts are more Mg rich or Ca rich compared                 are recognized in all the samples, phenocrysts with An
with other phenocrysts in the megacryst-free rocks.                   590 are rare in the 1811 samples.
  Plagioclase up to 30 mm in length forms the domi-                     Olivine is the dominant mafic phenocryst. The
nant phenocryst phase. The compositions of plagio-                    compositions in the 1811 and 1940 samples range from
clase phenocrysts are in the range An ˆ 53---97,                      Fo ˆ 70 to Fo ˆ 86 (Fig. 3). The 1874 lavas, however,
nearly the same as those in the megacryst-free rocks                  also include Fo-poor phenocrysts (Fo 570). Olivine
(Fig. 2). Plagioclase megacrysts (longer than 5 mm) in                megacrysts (longer than 2 mm) have similar composi-
all the rocks have An 490. Although bimodal distri-                   tions (Fo 480) among the three eruptions.

                                                                  2120
AMMA-MIYASAKA AND NAKAGAWA                          MAGMA EVOLUTION, MIYAKE-JIMA

                                                                                                                                        Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
Fig. 3. Core compositions of olivine and pyroxene phenocrysts in megacryst-free and megacryst-bearing rocks erupted during the period
1469---1983. Continuous lines with numerals are isotherms of Lindsley (1983). All data are plotted in terms of quadrilateral pyroxene
components diopside---enstatite---hedenburgite---ferrosilite.

  Clinopyroxene less than 40 mm in length is also an                to andesite, whereas some of the megacryst-bearing
abundant mafic phenocryst. All of the phenocrysts                   rocks plot in the low-K field. Although lavas from the
are augite (Fig. 3), and their Mg-numbers are in                    1469 to the 1983 eruptions show variations in
the range 60---83. The compositional range in the                   SiO2 (506---564 wt %) and K2O (035---082 wt %),
rocks is nearly the same as that for the megacryst-                 there are no systematic changes in their whole-rock
free rocks. The 1811 megacryst-bearing rocks show a                 composition with time (Fig. 5). Comparing the differ-
compositionally bimodal distribution of the clino-                  entiated megacryst-free rocks in each eruption,
pyroxene phenocrysts (Fig. 3). The phenocrysts with                 SiO2 and K2O contents seem to increase from 1469
high Mg-number often coexist with plagioclase and                   to 1712, then decrease to 1940 and increase again to
olivine megacrysts.                                                 1983. On the other hand, the chemical compositions of
  Orthopyroxene and magnetite phenocrysts are rare,                 the differentiated megacryst-bearing rocks in each
and smaller than 12 mm and 03 mm in length, respec-               eruption are nearly constant, whereas the entire com-
tively. The compositional range of orthopyroxene phe-               positional ranges of SiO2 and K2O have gradually
nocrysts is narrow (Mg-number ˆ 62---73) compared                   become wider from 1811 to 1940. The compositional
with the megacryst-free rocks (Fig. 3).                             variations of both megacryst-free and megacryst-
                                                                    bearing rocks are shown in representative Harker
                                                                    diagrams in Fig. 6.
WHOLE-ROCK GEOCHEMISTRY                                               The SiO2 contents of the megacryst-free rocks range
All of the rocks during the period 1469---1983 are clas-            from 506 to 564 wt % (Figs 4---6). Among these rocks,
sified as tholeiites in a SiO2---FeO*/MgO diagram                   at least three distinct trends (I, II and III) can be
(Fig. 4). Most of the lavas consist of medium-K basalt              recognized, especially in SiO2---TiO2, SiO2---FeO*/

                                                               2121
JOURNAL OF PETROLOGY                     VOLUME 44        NUMBER 12         DECEMBER 2003

                                                                                                                                         Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
                                                                     Fig. 5. Temporal variations of whole-rock SiO2 and K2O contents
                                                                     for megacryst-free and megacryst-bearing rocks erupted during the
Fig. 4. (a) SiO2---FeO*/MgO and (b) SiO2---K2O diagrams for          period 1469---1983.
megacryst-free and megacryst-bearing rocks erupted during the per-
iod 1469---1983. Dividing line in the SiO2---FeO*/MgO diagram is
from Miyashiro (1974) and that in the SiO2---K2O diagram from Gill
(1981). All analyses are normalized to 100 wt % volatile-free with   FeO*/MgO and higher Cr (e.g. No. 1522 of 1983 and
total iron (FeO*) calculated as FeO.                                 No. 1918 of 1835 in Table 2), compared with other
                                                                     1763---1983 lavas.
                                                                       The SiO2 content of the megacryst-bearing rocks
MgO and SiO2---Cr diagrams, and these seem to corres-                ranges from 514 to 555 wt % (Figs 4---6). With an
pond to eruption age (Fig. 7). Trend I consists mainly               increase of SiO2, the Al2O3, Cr and Ni contents
of the 1469, 1535 and 1595 lavas. Two samples (e.g.                  decrease, whereas the TiO2 contents and the FeO*/
No. 2222 in Table 2) from one of the lava lobes of the               MgO ratio increase. The compositional variation of
1643 eruption (main-N, Fig. 1) are also classified as                the megacryst-bearing rocks in each eruption forms
part of this trend. The SiO2 content of the Trend I                  distinct trends in most of the variation diagrams
rocks ranges from 518 to 549 wt %. This trend is the               (Fig. 6). Although the 1811 lavas consist of SiO2 ˆ
most mafic (poor in SiO2, TiO2, FeO*/MgO and rich                    546---552 wt % andesites, the 1874 and 1940 lavas
in Cr) of the three trends. Trend II consists mainly of              have a wider compositional range of SiO2 ˆ 528---
lavas from the 1643 main-S lava lobe (Fig. 1; No. 1620               552 wt % and SiO2 ˆ 514---555 wt %, respectively.
in Table 2) and 1712 lavas. Compared with the rocks                  The 1811 lavas are rich in TiO2 and have higher
of Trend I, those of Trend II are in the range of SiO2               FeO*/MgO, whereas the 1874 lavas are poor in TiO2
506---532 wt %, and are characterized by lower                     and FeO*/MgO among the megacryst-bearing rocks.
SiO2 and Cr, and by higher TiO2 and FeO*/MgO.                        Comparing the whole-rock compositions of the
Trend III is composed mainly of the 1763---1983 lavas,               megacryst-bearing rocks with the megacryst-free sam-
and also contains the 1643 lavas from the central lava               ples, the megacryst-bearing rocks are rich in Al2O3 and
lobe (Fig. 1; No. 1627 in Table 2). The rocks of this                Ni, and poor in TiO2, FeO*/MgO and V. Further-
trend have SiO2 ˆ 520---564 wt %, and are the most                 more, the felsic ends of the trends of variations of
differentiated (lowest Cr, and highest SiO2, TiO2 and                the megacryst-bearing rocks intersect with those of the
FeO*/MgO). There also exist some samples with lower                  megacryst-free samples at around SiO2 ˆ 55---56 wt %.

                                                                 2122
AMMA-MIYASAKA AND NAKAGAWA                            MAGMA EVOLUTION, MIYAKE-JIMA

                                                                                                                                           Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
Fig. 6. Selected SiO2 variation diagrams for major elements (TiO2, Al2O3), FeO*/MgO and trace elements (V, Cr and Ni) for megacryst-free
and megacryst-bearing rocks erupted during the period 1469---1983.

MINERAL ASSEMBLAGE AND                                               are, however, nearly aphyric and have less than 3 vol. %
                                                                     phenocrysts. This suggests that plagioclase accumula-
CHEMICAL COMPOSITIONS OF
                                                                     tion is a minor source of compositional variation in
CRYSTAL-CLOTS                                                        these rocks. In addition, SiO2 and K2O do not system-
It is widely accepted that two types of magmas                       atically change during the period 1469---1983 (Fig. 5).
have erupted from the tholeiitic volcanoes of the                    If differentiation had proceeded in a closed magma
Izu---Mariana arc: `plagioclase-controlled' (plagioclase-            storage system, these values should have increased
accumulated) and `differentiated' magmas (Nakano &                   with time. These observations suggest that simple dif-
Yamamoto, 1991; Nakano et al., 1991; Tsukui &                        ferentiation does not play an important role during the
Hoshino, 2002). This is consistent with the evidence                 period. On the other hand, phenocryst minerals show
that there exist no disequilibrium phenocryst assem-                 wide and polymodal compositional distributions (Figs 2
blages, such as olivine and quartz, or olivine and horn-             and 3), suggesting that all of the phenocrysts could
blende, in all of the rocks formed during the 1469---1983            not crystallize simultaneously from a single magma.
eruptions of Miyake-jima volcano. Most of these rocks                To investigate whether these rocks are equilibrium

                                                                 2123
JOURNAL OF PETROLOGY                     VOLUME 44        NUMBER 12       DECEMBER 2003

                                                                      crystallization products or not, we examined the
                                                                      mineral assemblage and chemical compositions of
                                                                      crystal-clots found in each sample, as described by
                                                                      Amma-Miyasaka & Nakagawa (2002), because it can be
                                                                      considered that minerals in the same crystal-clot have
                                                                      crystallized simultaneously from the same magma.

                                                                      Types of crystal-clots
                                                                      On the basis of the mineral assemblage and the chemi-
                                                                      cal compositions of the minerals within the crystal-clots
                                                                      in the 1469---1983 rocks, crystal-clots can be divided

                                                                                                                                  Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
                                                                      into four types (Figs 8 and 9): megacryst-type (M-type);
                                                                      basaltic-type (B-type); andesitic-type (A-type); basal-
                                                                      tic andesitic-type (AB-type). The M-type crystal-clots
                                                                      are recognized in megacryst-bearing rocks, and mostly
                                                                      consist of Ca-rich plagioclase and Mg-rich olivine with
                                                                      occasional Mg-rich clinopyroxene (e.g. in the 1811
                                                                      lavas). They are the least differentiated of the four
                                                                      types of clots: An ˆ 88---97, Fo ˆ 77---85 and Mg-
                                                                      number ˆ 77---84. The B-type crystal-clots consist of
                                                                      plagioclase, clinopyroxene and olivine, characterized
                                                                      by the absence of orthopyroxene and magnetite.
                                                                      Mineral core compositions are An ˆ 72---94, Mg-
                                                                      number ˆ 68---79 and Fo ˆ 62---78, respectively, more
                                                                      evolved than those of the M-type clots. In contrast, the
                                                                      A-type crystal-clots are characterized by the presence
                                                                      of magnetite (Usp ˆ 24---46) and orthopyroxene (Mg-
                                                                      number ˆ 59---73), as well as plagioclase, clino-
                                                                      pyroxene and olivine. Core compositions of plagio-
                                                                      clase, clinopyroxene and olivine in these clots are An
                                                                      ˆ 55---91, Mg-number ˆ 60---77 and Fo ˆ 60---73,
                                                                      slightly more evolved compared with those in the B-
                                                                      type clots. In addition to these clots, a fourth type of
                                                                      crystal-clot (AB-type)exists in the 1535 and 1595
                                                                      lavas. The crystal-clots contain Mg-rich orthopyrox-
                                                                      ene (Mg-number ˆ 70---76; Fig 3 and 9) without mag-
                                                                      netite. The core compositions of coexisting plagioclase
                                                                      and clinopyroxene are An ˆ 73---85 and Mg-number ˆ
                                                                      70---75, respectively. The mineral assemblage and
                                                                      chemical compositions are intermediate between the
                                                                      A-type and B-type clots.

                                                                      Compositional relationship of minerals
                                                                      within and among the clots
                                                                      The compositional relationship between olivine and
                                                                      pyroxenes can be investigated on the basis of the
                                                                      Fe---Mg distribution in these minerals. If these minerals
Fig. 7. Variation of TiO2, FeO*/MgO, and Cr vs SiO2 for               coexist in equilibrium, the Mg-number of olivine
megacryst-free rocks ( symbol in Fig. 6) erupted during the period   should be nearly the same as or slightly lower than
1469---1983. The megacryst-free rocks can be divided into three
trends. The 1643 lavas flowed down separately into three lobes:
                                                                      that of clinopyroxene and orthopyroxene (Obata et al.,
main N, main S and central lobes, which respectively correspond to    1974; Brey & Kohler, 1990). We have determined
Trend I, II and III.                                                  the average compositions and compositional ranges of

                                                                  2124
AMMA-MIYASAKA AND NAKAGAWA                            MAGMA EVOLUTION, MIYAKE-JIMA

                                                                      equilibrium from the same magma, and that different
                                                                      types of crystal-clots have crystallized from distinct
                                                                      magma compositions.
                                                                        The Fo content of olivine phenocrysts is strongly
                                                                      dependent on the FeO/MgO of the magma and
                                                                      KD(olivine/liquid) is constant over a wide range of P, T,
                                                                      fO2and H2O (Roeder & Emslie, 1970). Based on the Fo
                                                                      content of olivine in each type of clot, M-type clots are
                                                                      formed in the most mafic magma, and the A-type clots
                                                                      in the most differentiated magma. The presence or
                                                                      absence of magnetite and orthopyroxene in these clots
                                                                      suggest that the magma producing the M-type and B-

                                                                                                                                   Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
                                                                      type clots is basaltic, and that producing the AB-type
                                                                      and A-type clots is probably andesitic (e.g. Gill, 1981).

                                                                      DISCUSSION
                                                                      Origin of M-type clots and megacrysts
                                                                      Amma-Miyasaka & Nakagawa (2002) investigated the
                                                                      origin of M-type clots and megacrysts in the 1940 lava.
                                                                      We defined the M-type clots and megacrysts by their
                                                                      size [plagioclase of L 403 mm (where L is the radius of
                                                                      the equivalent circle) and olivine of L 402 mm] and
                                                                      the least differentiated composition of minerals (plagio-
                                                                      clase of An ˆ 88---97, olivine of Fo ˆ 77---85 and
                                                                      clinopyroxene of Mg-number ˆ 77---84). Examining
                                                                      these crystals in detail, we pointed out the following
                                                                      petrographical features: (1) spherical olivine in plagio-
                                                                      clase; (2) wide, homogeneous cores of plagioclase;
                                                                      (3) kink-banding of olivine megacrysts. These features
                                                                      cannot be explained by normal crystallization
                                                                      processes. The spherical olivine and wide, homoge-
                                                                      neous cores of plagioclase can be formed by long-term
                                                                      diffusion, and the kink bands can be formed under
                                                                      conditions that can transmit strain. Such a condition
                                                                      could be achieved within a plutonic body. This is
                                                                      supported by the evidence that the mineral assem-
                                                                      blage, chemical compositions and crystal size distribu-
                                                                      tions of both the megacrysts and M-type clots are
                                                                      identical to those of plutonic xenoliths (smaller than
                                                                      5 cm in diameter, consisting only of An ˆ 93---97
                                                                      anorthite and Fo ˆ 83---85 olivine crystals) found in
Fig. 8. Photomicrographs illustrating the various types of crystal-   some of the 1940 lavas.
clots (M, B, A and AB) during the period 1469---1983. The M-type         M-type crystal-clots have also been recognized in the
photograph is through crossed nicols and the other photographs are    megacryst-bearing rocks of the 1811 and 1874 erup-
in plane-polarized light.
                                                                      tions. The clots in the rocks of the 1874 eruption consist
                                                                      only of plagioclase and olivine, as in the case of the
olivine and pyroxenes in each type of clot of each                    1940 megacryst-bearing lavas. On the other hand, the
eruption age (Fig. 10). The Mg-number of clinopyrox-                  clots in the 1811 lavas also contain clinopyroxene.
ene is slightly higher compared with that of olivine                  There is no significant difference in the petrographic
in each type of clot. Moreover, the Mg-number of                      characteristics and chemical composition of plagioclase
orthopyroxene is nearly the same as that of olivine in                and olivine megacrysts among the three lavas of known
the A-type clots. This suggests that the olivine and                  age. Thus, we consider that the megacrysts in the 1811
pyroxene in each type of clot have crystallized in                    and 1874 lavas have the same origin as those in the

                                                                  2125
JOURNAL OF PETROLOGY                      VOLUME 44          NUMBER 12          DECEMBER 2003

                                                                                                                                                Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
Fig. 9. Mineral assemblage and mineral core compositions (plagioclase, clinopyroxene and olivine) in four types of crystal-clots (M, B, A and
AB) in the rocks erupted during the period 1469---1983 (see text).

                                                                        from N---NE flank fissures. Therefore, it may be sug-
                                                                        gested that these lavas have captured both the mega-
                                                                        crysts and M-type clots as xenocrysts from the same
                                                                        body located beneath the N---NE sector of the volcano.

                                                                        Evidence for magma mixing
                                                                        Both A-type and B-type crystal-clots coexist in each
                                                                        sample of the 1469---1983 lavas, except for 1835, 1962
                                                                        and 1983. Histograms of minerals occurring in A-, B-
                                                                        and AB-type crystal-clots in each eruption are shown
                                                                        in Fig. 11. Plagioclase in these lavas shows composi-
                                                                        tionally bimodal distribution. The peaks of the An
                                                                        contents are variable among lavas of different eruption
                                                                        ages. Plagioclase phenocrysts in the B-type clots are,
                                                                        however, always more Ca rich than those in the A-type
                                                                        clots in each eruption. Histograms of olivine and clino-
                                                                        pyroxene do not clearly show compositionally bimodal
                                                                        distributions; however, these minerals are always more
                                                                        Mg rich in the B-type clots. This coexistence of two
                                                                        types of crystal-clots with distinct chemical composi-
                                                                        tions is considered to be a disequilibrium feature.
                                                                           The relationship between the crystal-clots and their
                                                                        host rocks can be investigated on the basis of Fe---Mg
Fig. 10. Mg-number of (a) clinopyroxene core and (b) orthopyrox-
                                                                        partitioning (Fig. 12). If minerals crystallize in equili-
ene core vs Mg-number of olivine core in the A-, B- and M-type clots    brium with the host magma, Fe---Mg distribution
in the rocks erupted during the period 1469---1983. The AB-type clots   between the mineral-cores and the magma must be
lack olivine phenocrysts. Dashed line in (a) is after Obata et al.      constrained by plausible partition coefficients, shown
(1974).
                                                                        by the lines in Fig. 12. In the case of the rocks with both
                                                                        A-type and B-type clots, most of the mafic minerals in
1940 lavas. Although the whole-rock chemistry of the                    the clots do not plot near equilibrium. Most olivine and
megacryst-bearing rocks of each eruption seems to                       clinopyroxene in the B-type clots are more Mg rich
define distinct trends, this may reflect the difference                 than possible equilibrium, whereas those in the A-
in the ratio of component minerals in the plutonic                      type clots are more Fe rich. Orthopyroxene, contained
rocks. All of the megacryst-bearing rocks erupted                       only in the A-type clots, is also more Fe rich than that

                                                                   2126
AMMA-MIYASAKA AND NAKAGAWA                             MAGMA EVOLUTION, MIYAKE-JIMA

                                                                                                                                               Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
Fig. 11. Histograms of core compositions of plagioclase, clinopyroxene, orthopyroxene and olivine in A-, B- and AB-type crystal-clots in the
rocks erupted during the period 1469---1983. The M-type clots are considered to be xenocrysts and are excluded. It should be noted that the
rocks containing only B-type clots have never erupted during the period 1469---1983.

required to be in equilibrium with the whole-rock                      (An ˆ 60---90). The compositional variation of clino-
FeO*/MgO ratio (Fig. 12). This suggests that the                       pyroxene in the 1983 lavas is also relatively extensive,
mafic minerals in the B-type clots should have crystal-                compared with that in the 1962 and 1835 lavas. The
lized from a more mafic (lower FeO*/MgO) magma                         compositional relationships between minerals and
than the host magma, and those of the A-type clots from                whole rocks also show that most of the mafic minerals
a more differentiated (higher FeO*/MgO) magma.                         in these rocks are in equilibrium with their host magma;
   Our data suggest that most of the 1469---1983 lavas                 however, Fe-rich clinopyroxene in the 1983 sample
contain two types of phenocryst that originated from                   (No. 1522 in Table 2) and olivine in the 1962 sample
distinct magmas, and indicate that these lavas could be                (No. 122508 in Table 2) might not be in equilibrium
produced by magma mixing of two end-member mag-                        (Fig. 12). The lack of B-type clots indicates that these
mas, a differentiated A-type (andesitic) and a mafic                   lavas may have been produced without mixing with the
B-type (basaltic) magma. The 1535 and 1595 lavas                       mafic B-type magma. Although it is also possible that
that have another type of crystal-clot (AB-type,                       nearly aphyric basaltic magma mixed with the A-type
Fig. 11), in addition to the A-type and B-type, might                  magma, the proportion of the basaltic magma would be
be the results of mixing between three magmas. In                      minor, because the mineral compositions are nearly in
summary, at least two distinct magmas have existed                     equilibrium with the host magma (Fig. 12). Consider-
beneath the Miyake-jima volcano since 1469, and have                   ing the wide compositional variations of the phenocryst
usually mixed during eruption.                                         minerals (plagioclase and clinopyroxene) especially
                                                                       in the 1983 lavas, however, we could not exclude
                                                                       the possibility of mixing between distinct andesitic
Rocks without obvious evidence                                         magmas. In conclusion, we suggest that lavas without
for magma mixing                                                       the B-type clots are not mixing products between the
The 1983, 1962 and 1835 lavas contain only the A-type                  A-type and B-type magmas, and that only pure A-type
crystal-clots (Fig. 11). Although histograms of phe-                   magma (or magmas) has erupted during the 1835,
nocryst compositions in these lavas do not show                        1962 and 1983 eruptions. Furthermore, it is noteworthy
clear bimodal distributions, the compositional range                   that rocks containing only the B-type clots have never
of the plagioclase phenocrysts is relatively wide                      erupted during the period 1469---1983.

                                                                   2127
JOURNAL OF PETROLOGY                       VOLUME 44           NUMBER 12          DECEMBER 2003

                                                                                                                                                     Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
Fig. 12. Fe---Mg partitioning between whole rock and cores of A-type and B-type olivine and pyroxenes from 1469---1983 megacryst-free rocks.
The megacryst-bearing rocks are excluded, because the FeO*/MgO of these rocks may be affected by the presence of xenocryst minerals. Lines
represent partition coefficients of KD [(XFeO/XMgO)mineral/(XFeO/XMgO)liquid] of Beattie (1993) for olivine and orthopyroxene, and Baker &
Eggler (1987) for clinopyroxene. FeO contents of whole rocks are calculated by assuming Fe2 ‡ /(Fe2 ‡ ‡ Fe3 ‡ ) ˆ 09. We have also tried Fe2 ‡ /
(Fe2 ‡ ‡ Fe3 ‡ ) ˆ 081 by wet chemical analysis (Isshiki, 1960; Iwasaki et al., 1982), but 09 better explains the equilibration relationships in
non-mixing rocks.

                                                                          Except for the 1835, 1962 and 1983 lavas, the FeO*/
End-member magmas and their                                               MgO ratios of the whole rocks are intermediate
relationships                                                             between those of the A-type and B-type magmas in
Estimation of whole-rock compositions of                                  each eruption. This is consistent with our conclusion
end-member magmas                                                         that these rocks are mixing products between these
Using the Fe---Mg mineral---melt distribution coeffi-                     magmas. Whole-rock FeO*/MgO ratios of the 1962
cients for olivine/liquid and orthopyroxene/liquid                        and most of the 1835 and 1983 lavas are nearly the
(Beattie, 1993), we can estimate the FeO*/MgO ratios                      same as the estimated FeO*/MgO ratios of the A-type
of both the A-type and B-type magmas (Fig. 13).                           magmas, whereas the whole-rock ratios of the 1535,

                                                                      2128
AMMA-MIYASAKA AND NAKAGAWA                              MAGMA EVOLUTION, MIYAKE-JIMA

                                                                                                                                                 Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
Fig. 13. Estimated FeO*/MgO ratios of A-type and B-type magmas erupted during the period 1469---1983. We used KD(olivine/liquid) ˆ 0303,
KD(orthopyroxene/liquid) ˆ 0284 of Beattie (1993) and the representative core compositions of mafic minerals in each type of clot. Calculated
FeO*/MgO values of the A-type liquid are plotted at the SiO2-rich end and those of the B-type liquid at the SiO2-poor end.

1595, 1712 and most of the 1643 lavas are similar to                    27 km in radius and 100 m thick for the 1940 eruption,
those of the B-type magmas. FeO*/MgO ratios might                       assuming the depth of the magma to be 3 km. Hypo-
reflect the mixing ratio between the A-type and B-type                  centres of long-period earthquakes (Ueki et al., 1984)
magmas.                                                                 and geomagnetic change (Nakagawa et al., 1984) asso-
                                                                        ciated with the 1983 eruption suggested that magma
                                                                        was located at a depth of 2---3 km. In addition to these
Relationship of the end-member magmas                                   shallower `sources', the inflation and deflation accom-
Plots of incompatible-element ratios (K2O/Ba, Y/Ba                      panying the 1983 eruption suggested a magma reser-
and Y/Zr) in the 1469---1983 megacryst-free rocks                       voir about 1---2 km SW of the summit crater, at a depth
define linear trends that pass through the origin with                  estimated to be about 8 km (Tada & Nakamura,
R2 values of 089---093 (Fig. 14). This strongly suggests              1988). Sasai et al. (2001) also suggested that the pres-
that both the A-type and B-type magmas may have                         sure source during 1995---1999 lies beneath the south-
been derived from a single primary magma. This con-                     ern flank of the volcano at a depth of 8 km below sea
clusion is consistent with Sr isotopic data. Notsu et al.               level. Considering that the depths of short-period
(1983) and Notsu & Aramaki (1984) measured 87 Sr/                       earthquakes are 1---10 km below sea level (Minakami
86                                                                      et al., 1963; Miyazaki & Sawada, 1984; Ueki et al.,
   Sr ratios of the rocks from Miyake-jima volcano, and
confirmed that the isotopic ratio is concentrated in a                  1984), the magma storage systems appear to have
narrow range (070350---070369), suggesting that the                   existed at least at two depths (c. 2---3 km and 8 km).
source region of the magmas beneath Miyake-jima                         On the other hand, geobarometry did not successfully
volcano is isotopically uniform.                                        reveal the crystallization depth of the A-type and
                                                                        B-type magmas of the volcano. We used the pseudo-
                                                                        ternary normative diagram (cpx---ol---SiO2) of Walker
Depths of the end-member magmas                                         et al. (1979) and Baker & Eggler (1983) for lavas with
The depths of the magma storage systems have been                       only the A-type clots and those with both the A-type
estimated based on geophysical observations carried                     and B-type clots, and the clinopyroxene geobarometer
out since the 1940 eruption. Takahashi & Hirano                         (Nimis, 1995) for both types of clots. These calculations
(1941) estimated that the size of the chamber was                       show, however, no difference between the A-type and

                                                                    2129
JOURNAL OF PETROLOGY                      VOLUME 44        NUMBER 12      DECEMBER 2003

                                                                        (A-type) and a deeper basaltic (B-type), have existed
                                                                        during the period 1469---1983. The relative position of
                                                                        these magma storage systems may be consistent with
                                                                        the fact that the A-type magma (or magmas) has
                                                                        erupted alone during the 1835, 1962 and 1983 erup-
                                                                        tions, whereas rocks containing only the B-type clots
                                                                        have never erupted during the period 1469---1983. The
                                                                        absolute pressure of both magmas should be less than
                                                                        3 kbar; the shallower andesitic magma storage systems
                                                                        might be located at a depth of 2---3 km, and the deeper
                                                                        basaltic magma storage systems might be at 8 km.

                                                                                                                                   Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
                                                                        Magma storage systems and
                                                                        their evolution
                                                                        Our petrological analysis of the 1469---1983 lavas
                                                                        reveals that the magma plumbing system during this
                                                                        period has involved two distinct magma storage
                                                                        systems, filled with the andesitic (A-type) and basaltic
                                                                        (B-type) magmas. In many of the eruptions, the resul-
                                                                        tant lavas are the mixing products between these two
                                                                        magmas, although shallower andesitic magmas have
                                                                        sometimes erupted alone. Because we can identify
                                                                        phenocryst minerals that are derived from each end-
                                                                        member magma during mixing events, the evolution
                                                                        of both end-member magmas can be studied by using
                                                                        these phenocrysts. We can demonstrate the evolution
                                                                        of the chemical composition of these phenocrysts and of
                                                                        the deduced magmatic temperature during the past
                                                                        500 years.

                                                                        Evolution of the end-member magmas and their
                                                                        magmatic temperature
                                                                        The core composition of olivine in the B-type crystal-
                                                                        clots has gradually become Fe rich with time (Fig. 15).
                                                                        The average composition has decreased systematically
                                                                        from Fo ˆ 75 in the 1469 lavas to Fo ˆ 69 in the 1940
                                                                        lavas. Although the compositional range of clinopyr-
Fig. 14. Incompatible-element variation diagrams for the mega-          oxene in the B-type magma is wide, the average pyr-
cryst-free rocks erupted during the period 1469---1983. R2 values are   oxene composition in each eruption is nearly the same
for all points.                                                         in the range of Mg-number ˆ 72---74 throughout the
                                                                        period 1469---1983. According to Brey & Kohler
                                                                        (1990), the ratio of the Mg-number of clinopyroxene
B-type magmas. Whole-rock compositions are pro-                         to the Fo content of olivine increased from 1 to more
jected around 1 atm for both lavas on the normative                     than 1, with falling equilibrium temperature. There-
diagram, and equilibration pressures calculated based                   fore, a temporal change in both Mg-number of
on clinopyroxene compositions are 0---3 kbar for both                   clinopyroxene and Fo content of olivine in the B-type
the A-type and B-type clots.                                            clots would indicate that the temperature of the basal-
  The B-type magma appears to have entrained plu-                       tic magma has fallen systematically with time. The
tonic xenoliths before mixing with the A-type magma                     temperature of the basaltic magma calculated by
during the 1940 eruption (Amma-Miyasaka &                               Loucks (1996) actually fell from 1220 C in the 1469
Nakagawa, 2002). Thus, we suggest that basaltic                         eruption to 1180 C in 1940. Using the temporal
magma should exist at a deeper level beneath the                        variation of phenocryst minerals in the B-type clots,
volcano than andesitic magma. We propose that two                       the evolution of the whole-rock chemistry of the basal-
distinct magma storage systems, a shallower andesitic                   tic magma can also be evaluated. The FeO*/MgO

                                                                    2130
AMMA-MIYASAKA AND NAKAGAWA                           MAGMA EVOLUTION, MIYAKE-JIMA

                                                                                                                                         Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
Fig. 15. Temporal variations of core compositions of B-type clinopyroxene and olivine and magmatic temperature for the rocks erupted
during the period 1469---1983. Numbers in parentheses represent the eruption intervals of B-type magma. The B-type magma erupted every
48---69 years during the period 1469---1983. Large filled triangles are average core compositions of the minerals.

ratios of the basaltic end-member magmas during                     a complex fashion from 1070 C in the 1469 eruption
1469---1595 were the lowest (Fig. 13). Since then, the              to 1100 C in 1983 (Fig. 16). As in the case of the
ratio has increased, accompanied by falling magmatic                basaltic magma, temporal changes in the whole-rock
temperature.                                                        chemistry of the andesitic magma are also evident: the
  In contrast, the core composition of clinopyroxene                FeO*/MgO ratio has decreased from 44 in the 1469
and orthopyroxene in the A-type clots has not changed               eruption to 33 in 1983 (Fig. 13).
systematically as it has for the B-type clots (Fig. 16).              Thus, the two magmas, which have formed end-
The average core composition of clinopyroxene is Mg-                member components for magma mixing during the
number ˆ 67 in the 1469 lavas, which is the most Fe                 past 500 years, have evolved following separate paths
rich during the past 500 years. Although the composi-               (Fig. 17a). The basaltic magma has become more
tion has fluctuated since then, it has gradually become             differentiated with time while cooling, whereas the
more Mg rich with time to Mg-number ˆ 72 in the                     andesitic magma has become more mafic while heat-
1983 lavas. Except for the 1763 eruption, the average               ing. A temporal fall in the magmatic temperature of
Mg-number of orthopyroxene in the 1469 A-type clots                 the basaltic magma could be explained by fractional
reflects the most differentiated magma during the per-              crystallization of a primary basaltic magma, and it
iod 1469---1983, and has gradually become Mg rich up                suggests that the basaltic end-member magma has
to 1983. Although orthopyroxene in the A-type clots in              been continuously in existence at least for 500 years.
the 1763 lava is the most Fe rich among the 1469---1983             The complicated and reverse evolution of the andesitic
samples, it is rare in 1763. Using the average composi-             magma during the past 500 years cannot be explained
tions of both clinopyroxene and orthopyroxene in the                by a closed system as in the case of the basaltic magma.
A-type clots, magmatic temperatures were calculated                 In each eruption, the deeper basaltic magma was
according to the Wells (1977) geothermometer for each               injected into the shallower andesitic magma and
eruption. The temperature appears to have risen in                  erupted a mixed magma. It suggests that the shallower

                                                                2131
JOURNAL OF PETROLOGY                     VOLUME 44          NUMBER 12         DECEMBER 2003

                                                                                                                                            Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
Fig. 16. Temporal variations of core compositions of A-type pyroxenes and magmatic temperature for the rocks erupted during the period
1469---1983. Numbers in parentheses represent the eruption intervals of A-type magma. The A-type magma erupted every 21---69 years during
the period 1469---1983, and about every 20 years since 1940. Large filled circles are average core compositions of the minerals.

andesitic magma chamber was an open magma                             by fractional crystallization of the basaltic magma. We
storage system for injection of the basaltic magma. In                used the composition of the 1535 basaltic lava (No.
this case, if only a small quantity of the basaltic magma             2321 in Table 3) as the starting material because the
was injected, the shallower andesitic magma would                     FeO*/MgO values of the 1535 and 1595 lavas
become more differentiated mainly by the effect of                    (240---248) are the lowest among the 1469---1983
fractional crystallization. In contrast, when basaltic                megacryst-free rocks, and nearly the same as the
magma was injected on a large scale, the andesitic                    calculated FeO*/MgO of the B-type liquid (24---25,
magma could have become more mafic. Moreover,                         Fig. 13). According to MELTS calculations, the
the andesitic magma might be nearly replaced by the                   FeO*/MgO continuously increases as temperature
basaltic magma. In conclusion, the temporal fluctua-                  falls, in the range of SiO2 553 wt % (Fig. 18a). Liquid
tion of magmatic temperature and the changing                         lines of descent are partially dependent on pressure and
composition of the andesitic magma during the past                    water content. Liquid compositions jump to SiO2 455
500 years depend on the quantity of basaltic magma                    wt % when magnetite begins to crystallize. The FeO*/
injected.                                                             MgO continuously increases as temperature falls, and
                                                                      liquid lines of descent for SiO2 455 wt % have a wide
                                                                      range. The compositional range of the 1469---1983
Evaluation of our model by MELTS calculation                          megacryst-free rocks lies between the liquid lines of
Our explanation for the evolution of the basaltic                     descent for SiO2 553 wt % and SiO2 455 wt %.
magma can be tested using the MELTS calculation                          Comparing the FeO*/MgO ratios of the end-
(Ghiorso & Sack, 1995). The chemical evolution of the                 member magmas calculated from phenocryst compo-
andesitic magma might reflect mixing with injected                    sitions (44 for the 1469 andesitic magma and 24---32
basaltic magma. The andesitic magma existing before                   for the basaltic magma) with the results of MELTS
the 1469 eruption, however, might have been formed                    calculations, the FeO*/MgO of the basaltic magma

                                                                 2132
AMMA-MIYASAKA AND NAKAGAWA                      MAGMA EVOLUTION, MIYAKE-JIMA

                       Table 3: Whole-rock                    magma during the past 500 years, the composition of
                       compositions of end-member             basaltic magma after the 1983 eruption must be SiO2
                                                              452 wt %, FeO*/MgO 432 and Fo content of the
                       magmas of the 1469 eruption            B-type olivine 570 (Fig. 17b). On the other hand, it is
                                                              difficult to estimate the composition of the shallower
                       Basaltic                Andesitic      andesitic magma because it shows a more complicated
Sample no.:            2321                    estimated      temporal change compared with the basaltic magma.
                                                              It has, however, become more mafic from 1874 to 1983,
wt %                                                          probably because of the injections of the basaltic
SiO2                    5150                   6277         magma. After the 1983 eruption, the andesitic
TiO2                     117                    092         magma may have had SiO2 553 wt %, FeO*/MgO
                                                              533 with A-type clots (pyroxenes; Mg-number 470).

                                                                                                                              Downloaded from https://academic.oup.com/petrology/article/44/12/2113/1460172 by guest on 14 January 2021
Al2O3                   1553                   1296
FeO*                    1265                    833
                                                              In other words, the compositions of andesitic and
MnO                      020                    052
                                                              basaltic magmas are indistinguishable after the 1983
                                                              eruption. Before the 2000 eruption, we proposed that
MgO                      525                    189
                                                              the erupted magma would be a mixing product of these
CaO                     1027                    632
                                                              basaltic and andesitic magmas, or would be andesitic
Na2O                     235                    356
                                                              alone.
K2O                      045                    110
                                                                Geshi et al. (2002) have already discussed the mag-
P2O5                     012                    031
                                                              matic system during the 2000 eruptive activity mainly
H2O                      050                    130
                                                              based on the whole-rock chemistry of the eruptive
Total                  10000                  10000         rocks. However, we consider that detailed analysis of
FeO*/MgO                 241                    440         phenocryst minerals, as investigated in this paper, are
                                                              essential to reveal the magmatic system. We will
                                                              develop detailed petrological analysis of the 2000
                                                              eruption in another paper (Amma-Miyasaka et al.,
                                                              submitted).
can be explained by both 1 and 3 kbar pressures
(Fig. 18b). Geophysical constraints suggest that the
magmas have been stored at a depth of 8 km; there-
fore, we consider that the basaltic magma may have            Implications for studies of the magma
stagnated at around 3 kbar pressure. On the other             plumbing system and its evolution in
hand, the FeO*/MgO ratio of the andesitic magma               active volcanoes
can be explained only by the 1 kbar calculated results        Several studies have already focused on the structures
(Fig. 18b). This result is consistent with the geophy-        of magma plumbing systems and their evolution in
sical observations (the shallower andesitic magma             active volcanoes over hundreds to thousands of years:
storage system might be located at a depth of                 Fichaut et al. (1989), Belkin et al. (1993), Villemant et al.
2---3 km). It is suggested that the 1469 andesitic            (1993) and Sato et al. (1996). These workers argued
magma (Table 3) could have been formed by frac-               that the increase of SiO2 content or FeO*/MgO ratio
tional crystallization of basaltic magma at a shallow         with time is due to fractionation or to the change of the
level (about 1 kbar). It is evident that the 1469 ande-       mixing ratio. On the other hand, they also speculated
sitic magma has changed its composition over 500              that a sudden drop of SiO2 or FeO*/MgO may indi-
years in a complex fashion (Fig. 18b), and it is also         cate the injection of a new batch of mafic magma. If the
consistent with our model that the shallower andesitic        end-member magmas had not changed for consider-
magma chamber was an open magma storage system                able durations, these explanations might be correct.
for injection of the basaltic magma.                          Temporal change of the end-member magmas, how-
                                                              ever, may be common as in the case of Miyake-jima
                                                              volcano, in which case temporal variations of the
Possible state of the magma plumbing                          whole-rock chemistry cannot simply correspond to the
system just before the 2000 eruption                          degree of fractionation and/or to mixing ratio of
Our model for the evolution of the magma plumbing             the end-member magmas. It should be noted that
system and eruption processes should be able to esti-         the detailed analysis of eruptive rocks from each
mate the state of the system just before the 2000 erup-       eruption using the mineral assemblage and chemical
tion of Miyake-jima volcano. According to simple              compositions in crystal-clots can clarify the evolution of
temporal differentiation of the deep-seated basaltic          the magma system of active volcanoes.

                                                           2133
You can also read