Coupling of Redox Conditions of Mantle Melting and ... - Springer Link

0 downloads 0 Views 2MB Size Report
melting of 1300°C, sulfide exists as sulfide liquid (SL). With allowance for this ..... gold concentrations in sulphide-bearing magma under oxi- dizing conditions ...
ISSN 0869-5911, Petrology, 2018, Vol. 26, No. 2, pp. 145–166. © Pleiades Publishing, Ltd., 2018. Original Russian Text © N.L. Mironov, M.V. Portnyagin, 2018, published in Petrologiya, 2018, Vol. 26, No. 2, pp. 140–162.

Coupling of Redox Conditions of Mantle Melting and Copper and Sulfur Contents in Primary Magmas of the Tolbachinsky Dol (Kamchatka) and Juan de Fuca Ridge (Pacific Ocean) N. L. Mironova, * and M. V. Portnyagina, b aVernadsky Institute of Geochemistry and Analytical Chemistry, Russian Academy of Sciences, ul. Kosygina 19, Moscow, 119991 Russia b GEOMAR Helmholtz Center for Ocean Research Kiel, Wischhofstrasse 1-3, 24148 Kiel, Germany *е-mail: [email protected]

Received October 22, 2016; in final form, January 23, 2017

Abstract⎯The compositions of parental melts of Tolbachinsky Dol (Kamchatka) basalts were estimated from the compositions of olivine-hosted (Fo90.5-83.1) primitive melt inclusions in the rocks of the Northern breakthrough of the Great Tolbachik Fissure Eruption (1975 A.C.) and of the late-Holocene cone “1004”. The parental melts contain 100–150 ppm Cu and 0.16–0.30 wt % S. These concentrations are much higher than those determined for the initial magmas of mid-ocean ridge basalts (MORB), for example of the Juan de Fuca ridge (Cu = 55–105 ppm, S=0.09–0.12 wt %). Modeling of mantle melting under variable redox conditions demonstrated that the high Cu and S contents in the Tolbachinsky Dol melts can be obtained by 6– 12% melting of DMM-like source under oxidized conditions (ΔQFM = +1.2 ± 0.1) and do not require a significant (>30–35% for S) subduction-related influx of these elements to the mantle source. The high contents of Cu and S in the Tolbachinsky Dol melts are largely explained by the increase of sulfide solubility in a silicate melt under oxidized conditions. In contrast, relatively reduced (ΔQFM ~ 0) conditions of MORB generation result in low contents of Cu and S in their initial magmas. The estimated ΔQFM values agree well with the data obtained using the Cr-spinel–olivine oxybarometer. The high oxygen potential of Tolbachinsky Dol primary magmas is inherited by more evolved magmas, thus favouring Cu enrichment up to 270 ppm during magma fractionation, approaching maximum copper contents in the global systematics of island-arc rocks. DOI: 10.1134/S0869591118020030

INTRODUCTION Redox regime usually described by the oxygen fugacity ( fO2 ) is one of the most important factors of the Earth’s evolution (e.g., Kadik, 2006; Ryabchikov and Kogarko, 2010). It determines the composition of mantle magmas (e.g., Green et al., 1987), their differentiation trends (e.g., Ariskin and Barmina, 2004; Portnyagin et al., 2012), and the behavior of many chemical elements, including chalcophile ore elements (for instance, Cu, Ag, and Au), in magmatic processes (Mungall, 2002; Borisov, 2005; Botcharnikov et al., 2011). The more oxidized nature of island-arc magmas and their sources as compared to oceanic magmas for a long time was regarded as reliably established fact (e.g., Ballhaus et al., 1991). However, this paradigm has been challenged in recent decades. Based on trace element (Sc, V, Zn, Сu) systematics and Fe isotopes, some researchers proposed that the oxygen fugacity during formation of island-arc magmas is close to that of the mid-ocean ridge basalts (Lee et al., 2005, 2012)

(MORB, ΔQFM ~ 0). These data also arised some doubts concerning the initial enrichment of the island-arc magmas in chalcophile elements, for instance, in Cu (Lee et al., 2012). Other authors found, however, several lines of additional evidence in support of the more oxidized state of primary islandarc magmas (Kelley and Cottrell, 2009; Evans et al., 2012). Variations in the redox conditions at ΔQFM > 0 during mantle melting can be traced using sulfur content in melts. For instance, according to experimental data for water-saturated basaltic melt at P = 0.2 GPa and T = 1050°С (Jugo et al., 2010), the total solubility of sulfur species (sum of S2– and S6+) sharply increases from 0.13 to 1.3 wt % in response to the increase of S6+ solubility. This effect is documented within a narrow ΔQFM range from 0 to +1.5 in the field of the possible sulfide stability. The S content in melt remains high (1.3–1.2 wt %) under more oxidized conditions (ΔQFM > 1.5) in the sulfate saturation field (Jugo et al., 2010). However, the information on the initial sul-

145

146

MIRONOV, PORTNYAGIN

fur content is inaccessible from the composition of subaerially erupted island-arc rocks (Wallace, 2005), while data on the sulfur content in melt inclusions need to be treated with caution, with allowance, for instance, of possible melt degassing before its trapping as inclusion (e.g., Sadofsky et al., 2008). The contribution of different sources (mantle vs. subduction component) to the sulfur content in magmas is also uncertain (Richards, 2015). All these complications seem to prevent researchers from direct application of data on sulfur content in natural melts to reconstruct the redox conditions of magmas and their sources (Lee et al., 2012). This work reports new data on the Сu and S contents in the parental melts of Tolbachinsky Dol at Kamchatka, which were obtained from study of melt inclusions in high-Mg olivine. Data on quenched MORB glasses from the Juan De Fuca Ridge, Pacific Ocean, were taken for comparison (Gale et al., 2013; Jenner and O’Neill, 2012). These data are used for modeling of peridotite melting under different redox conditions. Our results indicate that the mantle melting beneath Tolbachinsky Dol occured under oxidized conditions (ΔQFM ~ +1.1…+1.3), which were preserved during subsequent magma differentiation. The data serve in support of more oxidized conditions of generation of island-arc magmas as compared to those of MORB, but raise doubts as to significant fluid influx of Cu and S in their source. OBJECTS Volcanism of Tolbachinsky Dol Tolbachinsky Dol (TD) exists as an active volcanic zone during at least Holocene (~10 000 years) and produces fissure areal monogenic basaltic eruptions (The Great Tolbachik …, 1983). Owing to two recent large eruptions (GTFE in 1975–1976 and in 2012– 2013), TD is geologically well-studied region and continues to attract attention of researchers, including petrologists and geochemists (2012–13 Tolbachik eruption, 2015). It was shown (Portnyagin et al., 2015) that the observed rock diversity of TD can be explained by the open-system cyclic fractionation of the primary high-Mg moderate-K magma and mixing of its fractionation products. The primary magmas were likely derived from a mantle-wedge peridotite compositionally similar to the moderately depleted MORB source (Churikova et al., 2001; Portnyagin et al., 2007a), which was melted (up to ~10%) at temperature of 1250–1300°C and pressure ~2 GPa under the influence of slab-derived hydrous melts (Portnyagin et al., 2015). The rocks of TD are ascribed to the tholeiitic series (Portnyagin et al., 2007b, Fig. 3), with the “tholeiitic index” (THI) according to (Zimmer et al., 2010) of 1.04 (THI > 1 for tholeiitic and 8.0 wt %) of the JDF show the same range of Cu contents (52–105 ppm, on average 74 ppm; Fig. 2а, Table 2). Similar concentrations were also found in the melt inclusion glasses in olivine from the Garrett Fracture Zone in the East Pacific Rise (70–100 ppm, on average ~80 ppm; Fig. 2а; Beaudoin et al., 2007). Based on these data, the typical initial Cu contents in MORB are assumed to be close to the average copper content in the primitive JDF basalts (MgO > 8.0 wt %, Cu 74 ± 12 ppm; Fig. 2a, Table 2). It should be noted that these values are also close to the copper content in the primitive rocks of some other volcanoes of Kamchatka, for instance, Klyuchevskoy volcano (Cu 50–100 ppm, MgO > 10 wt %; Portnyagin and Manea, 2008; GEOROC, 2017). Unlike rocks, melt inclusions can retain information on the volatile content, in particular, S content, in primary magmas (Wallace, 2005; Sadofsky et al., 2008). The highest S content up to 0.33 wt % for TD was found in primitive MIs. Matrix glasses and bulk rocks have S contents close to the detection limit of microprobe analysis (~0.01 wt %), which is related to the substantial magma degassing during eruption (Fig. 2b). The primitive TD melts differ from MORB in much higher concentrations of S and Cu (Fig. 2). The high-Mg glasses (MgO > 8 wt %) of the JDF ridge contain from 0.09 to 0.12 wt % S (Jenner and O’Neill, 2012), whereas the primitive melts of TD contain 0.14–0.33 wt % S (Fig. 2b; Kamenetsky et al., 2017). Since sulfur variations in the TD melts can be related to the partial degassing of melts prior to their trapping by host mineral and to the possible presence of S-bearing daughter phases in the inclusions (Kamenetsky et al., 2002), the range of S content for primary TD magmas was estimated from the average and maximum sulfur contents in primitive MIs. The estimated S contents in primary melts of the Northern Breakthrough of the TD (sample NT) are 0.19–0.27 wt % (average-maximum), while those of the Cone “1004” (sample K01-30) are 0.21–0.30 wt % (Figs. 2b, 3). The upper estimate for the Northern Breakthrough was obtained using expanded data set of inclusions, including those not analyzed for copper content. The S content in the primary magmas of the JDF ridge was estimated to be ~0.11 wt % based on data on the most magnesium-rich glasses (MgO > 8 wt %) (Figs. 2b, Table 2). This value is well consistent with sulfur contents (0.08–0.11 wt %) in the most primitive melts trapped by high-Mg olivine phenocrysts (Fo88.6-90.7)

150

MIRONOV, PORTNYAGIN 300 (а) 250

Rocks and melts of TD

Cu, ppm

200 Kamchatka rocks 150

100

50 MORB (JDF) 0

2

4

6

8

10

12

14

0.35 (b) 0.30

TD melts

S, wt %

0.25 0.20 0.15 MORB (JDF) 0.10 0.05 Primary melts 0

2

4

Tolbachinsky Dol (TD): Primitive inclusions MORB:

6 8 MgO, wt %

10

12

14

Rocks (Portnyagin et al., 2015)

NT Northern Breakthrough (sample NT) Hybrid Cone “1004” (sample K01-30) K01-30 melts Sample NT n.a. for Cu Groundmass (sample NT) JDF (Jenner, O’Neill, 2012) JDF (Gale et al., 2013) Primary melts:

NT

MI in olivine K01-30

FZ Garrett TF Siqueiros MORB (JDF)

Fig. 2. Variations of Cu, S, and MgO contents in melt inclusions, matrix glasses, and rocks of Tolbachinsky Dol (TD) and Pacific MORB basalts of the Juan de Fuca (JDF) Ridge. Compositions of studied melt inclusions in olivine and matrix glasses from TD rocks are shown here. Calculated average compositions of initial, primary melts (recalculated to equilibrium with olivine Fo91) are shown by large symbols (see also their composition in Table 2). Additional set of primitive melt inclusions from NT sample analyzed only for sulfur content (not analyzed, n.a. for Cu) is shown with smaller (yellow) circles. The compositions of oceanic glasses of the JDF ridge basalts (>40° N) (Jenner, O’Neill, 2012; Gale et al., 2013), melt inclusions in olivine from basalts of the Garrett Fracture Zone (FZ), East Pacific Rise (Beaudoin et al., 2007), melt inclusions in olivine from the Siqueiros Transform Fault (TF), East Pacific Rise (Danyushevsky et al., 2003), compositions of the TD rocks (Portnyagin et al., 2015), and volcanic rocks of Kamchatka (GEOROC, 2017) are also shown for comparison. PETROLOGY

Vol. 26

No. 2

2018

COUPLING OF REDOX CONDITIONS OF MANTLE MELTING AND COPPER...

151

Table 1. Composition of melt inclusions in olivine and groundmass glasses in the studied samples of Tolbachinsky Dol (TD) Sample

NT

NT

NT

NT

NT

NT

NT

NT

NT

NT

Inclusion

1

2

3a

3b

4

5a

6

7

8

10

Experiment, Tquench, °C

no

no

no

no

no

no

no

no

no

no

Inclusion size, а, μm

77

85

79

85

60

189

45

89

117

94

Inclusion size, b

74

113

64

113

94

83

145

57

115

85

Averaged diameter, d

75

99

72

99

77

136

95

73

116

90

Bubble, diameter, μm

19

36

17

26

25

32

19

26

38

26

87.5

88.5

83.1

83.3

87.9

90.1

85.4

87.9

88.9

90.0

2

1

2

2

2

1

2

1

1

1

Host olivine, Fo Melt group (1, 2, 3)

Composition of inclusion glasses, wt % SiO2 TiO2 Al2O3 FeO MnO MgO CaO Na2O K2O P2O5 S Cl H2O Total CaO/Al2O3 CaO/K2O Cu, ppm

49.53 47.32 50.10 50.34 49.64 46.90 49.85 47.02 1.26 1.08 1.37 1.35 1.26 1.10 1.30 1.08 16.28 15.87 16.28 16.62 16.30 15.91 16.01 15.48 7.09 7.10 8.32 8.06 7.00 6.41 8.85 8.16 0.08 0.14 0.14 0.10 0.15 0.12 0.10 0.16 5.90 6.60 5.46 5.46 5.75 6.77 5.72 6.53 11.80 14.04 9.52 9.38 11.93 14.60 9.57 13.65 3.06 2.78 3.22 3.41 3.06 2.82 3.15 2.87 1.34 0.98 1.62 1.63 1.32 0.91 1.50 0.93 0.32 0.22 0.39 0.45 0.26 0.21 0.38 0.17 0.176 0.239 0.089 0.088 0.166 0.254 0.081 0.223 0.092 0.129 0.059 0.049 0.075 0.129 0.059 0.125 n.a. n.a. n.a. n.a. n.a. n.a. n.a. n.a. 96.91 96.49 96.56 96.93 96.90 96.12 96.57 96.38 0.73 0.88 0.58 0.56 0.73 0.92 0.60 0.88 8.84 14.38 5.88 5.75 9.07 15.97 6.38 14.64 182 143 201 206 168 145 209 163 Inclusion composition calculated to equilibrium with host olivine, wt %

47.96 47.18 1.01 1.09 15.04 16.22 7.20 6.45 0.08 0.14 6.60 5.83 13.85 14.56 2.65 2.76 0.94 0.86 0.16 0.20 0.202 0.280 0.116 0.131 n.a. n.a. 95.80 95.69 0.92 0.90 14.76 16.92 116 126

SiO2 TiO2 Al2O3 FeOtot MnO MgO CaO

49.6 1.15 14.87 9.5 0.07 9.54 10.78

47.84 0.99 14.6 9.5 0.13 10.16 12.92

51.22 1.36 16.11 9.5 0.14 6.89 9.42

51.13 1.32 16.26 9.5 0.1 6.95 9.18

49.56 1.13 14.68 9.5 0.14 9.89 10.75

47.2 0.97 13.99 9.5 0.11 11.76 12.84

50.92 1.27 15.61 9.5 0.1 8.16 9.33

48.01 1.03 14.73 9.5 0.15 9.61 12.99

48.68 0.93 13.8 9.5 0.07 10.7 12.71

47.51 0.94 14.04 9.5 0.12 11.79 12.6

Na2O

2.79

2.56

3.19

3.34

2.76

2.48

3.07

2.73

2.43

2.39

K2O

1.22

0.9

1.6

1.59

1.19

0.8

1.46

0.88

0.86

0.74

P2O5

0.29

0.2

0.39

0.44

0.23

0.18

0.37

0.16

0.15

0.17

Total Tcalc (1 atm, dry), °С

99.8

99.8

1254

1258

99.8 1195

99.8 1200

99.8 1262

99.8 1298

99.8 1230

99.8 1245

99.8 1271

99.8 1296

Corr. coeff.

0.91

0.92

0.99

0.98

0.90

0.88

0.97

0.94

0.92

0.86

Cl

0.084

0.119

0.058

0.047

0.068

0.113

0.057

0.118

0.106

0.113

S

0.160

0.221

0.088

0.086

0.150

0.223

0.079

0.211

0.185

0.241

Cu

166

132

198

201

152

127

204

154

106

108

S and Cu contents in a primary melt (Fo91) S

0.20

Cu

120 PETROLOGY

Vol. 26

No. 2

0.21 121

2018

0.19 137

0.17 97

0.23 103

152

MIRONOV, PORTNYAGIN

Table 1. (Contd.) Sample

NT

NT

NT

NT

NT

NT

NT

NT

NT

NT

Inclusion

11

12

13a

13b

14

15

16

17

18

19

Experiment, Tquench, °C

no

no

no

no

no

no

no

no

no

no

Inclusion size, а, μm

160

123

43

72

117

104

189

72

117

113

Inclusion size, b

208

75

72

40

49

60

79

40

121

34

Averaged diameter, d

184

99

58

56

83

82

134

56

119

74

Bubble diameter, μm

75

28

11

15

21

25

42

13

28

26

88.0

88.3

85.0

84.6

86.9

84.5

87.9

83.1

89.0

88.9

1

1

2

2

1

2

1

2

1

1

47.97

46.87

48.71

46.98

47.37

Host olivine, Fo Melt group (1, 2, 3) SiO2 TiO2 Al2O3 FeO MnO MgO CaO Na2O K2O P2O5 S Cl H2O Total CaO/Al2O3 CaO/K2O Cu, ppm

1.08 15.14 7.54 0.16 6.71 13.64 2.84 0.95 0.21 0.145 0.118 n.d. 96.49 0.90 14.40 154

Composition of inclusion glasses, wt % 49.31 49.65 47.25 50.12 47.25

1.08 1.32 1.32 1.05 1.39 1.11 15.92 16.30 16.21 15.71 16.58 15.72 7.47 8.65 8.91 8.19 7.78 7.18 0.11 0.15 0.18 0.14 0.12 0.11 6.49 5.53 5.62 6.57 5.54 6.37 13.65 10.75 9.47 13.08 9.30 13.45 2.75 3.23 3.28 2.87 3.36 2.79 1.03 1.41 1.51 0.95 1.59 0.99 0.24 0.30 0.37 0.18 0.41 0.20 0.238 0.133 0.080 0.232 0.086 0.211 0.139 0.081 0.061 0.123 0.066 0.130 n.d. n.d. n.d. n.d. n.d. n.d. 95.97 97.15 96.64 96.34 96.33 95.49 0.86 0.66 0.58 0.83 0.56 0.86 13.23 7.65 6.29 13.78 5.87 13.61 151 206 209 156 142 154

1.22 1.12 15.81 15.71 9.09 6.91 0.18 0.13 5.86 6.41 10.98 13.77 3.12 2.81 1.27 0.94 0.22 0.19 0.169 0.217 0.097 0.122 n.d. n.d. 96.72 95.30 0.69 0.88 8.65 14.64 145 148

1.05 15.64 7.42 0.13 6.51 13.89 2.89 0.92 0.17 0.246 0.124 n.d. 96.35 0.89 15.01 140

Composition of inclusions calculated to equilibrium with host olivine, wt % SiO2

48.63

47.81

TiO2

1.01

1.01

Al2O3

14.2

14.9

50.1

50.8

48.42

51.02

48.32

50.21

47.93

47.97

1.28

1.3

1.02

1.35

1.04

1.24

1.03

0.96

15.81

15.98

15.26

16.06

14.77

16.05

14.39

14.34

FeOtot

9.5

9.5

9.5

9.5

9.5

9.5

9.5

9.5

9.5

9.5

MnO MgO CaO Na2O

0.15 9.79 12.8

0.1 9.97 12.77

0.15 7.78 10.43

0.18 7.63 9.33

0.14 8.89 12.7

0.12 7.57 9.01

0.1 9.7 12.64

0.18 6.82 11.15

0.12 10.64 12.62

0.12 10.55 12.73

2.66

2.57

3.13

3.23

2.79

3.25

2.62

3.17

2.57

2.65

K2O

0.89

0.96

1.37

1.49

0.92

1.54

0.93

1.29

0.86

0.84

P2O5

0.2

0.22

0.29

0.36

0.17

0.4

0.19

0.22

0.17

0.16

Total Tcalc (1 atm, dry), °С Corr. coeff. Cl S Cu S Cu

99.8 1249 0.94 0.111 0.136 144 0.12 129

99.8 1253

99.9 1212

99.8 1217

99.8 1226

99.8 1218

99.8 1247

0.93 0.98 0.99 0.97 0.97 0.94 0.130 0.079 0.060 0.119 0.064 0.122 0.222 0.130 0.079 0.225 0.084 0.198 140 201 207 151 138 144 S and Cu contents in a primary melt (Fo91) 0.20 126

0.20 132

99.8 1177

99.8 1271

1.02 0.91 0.098 0.111 0.172 0.198 148 135

0.18 129 PETROLOGY

0.18 124 Vol. 26

99.8 1269 0.91 0.113 0.224 127 0.21 117 No. 2

2018

COUPLING OF REDOX CONDITIONS OF MANTLE MELTING AND COPPER...

153

Table 1. (Contd.) Sample

NT

NT

NT

NT

NT

Inclusion

20

4gm*

7gm

3gm*

15gm

2

3

4

5

6

Experiment, Tquench, °C

no

no

no

no

no

1180

1190

1250

1185

1167

8.40

10.00

14.45

11.15

7.10

Experimental time, min

K01-30 K01-30 K01-30 K01-30 K01-30

Inclusion size, а, μm

47

110

130

135

115

142

Inclusion size, b

94

110

110

130

130

112

Averaged diameter, d

71

110

120

133

123

127

Bubble, diameter, μm

26

30

35

50

42

37

Host olivine, Fo

90.3

73.9

1

3

Melt group (1, 2, 3)

78.9

77.8

77.7

88.7

87.6

90.5

87.9

87.1

3

3

3

1

1

1

1

2

Composition of inclusion glasses, wt % SiO2 TiO2 Al2O3 FeO MnO MgO CaO Na2O K2O P2O5 S Cl H2O Total CaO/Al2O3 CaO/K2O Cu, ppm

47.38 51.83 51.05 52.34 51.45 46.84 46.80 49.78 1.10 1.98 1.39 1.62 1.48 0.87 0.90 0.81 16.55 14.83 16.76 16.04 16.73 13.49 13.23 11.31 5.90 10.78 8.53 9.77 9.01 8.53 9.40 8.12 0.04 0.22 0.11 0.20 0.18 0.15 0.18 0.17 6.46 4.31 4.74 5.04 4.61 9.63 11.06 13.51 13.86 8.95 9.64 9.66 9.57 13.38 11.23 11.78 2.80 3.57 3.34 3.37 3.56 2.34 2.20 1.73 1.04 2.48 1.69 1.88 1.76 0.78 0.64 0.65 0.20 0.59 0.39 0.44 0.39 0.17 0.13 0.13 0.232 0.017 0.011 0.008 0.022 0.261 0.197 0.169 0.121 0.076 0.063 0.056 0.064 0.140 0.127 0.107 n.d. n.d. n.d. n.d. n.d. 1.8 2.8 3.0 95.67 99.64 97.72 100.42 98.82 98.39 98.93 101.26 0.84 0.60 0.58 0.60 0.57 0.99 0.85 1.04 13.30 3.61 5.70 5.14 5.44 17.06 17.61 18.03 153 272 258 219 227 129 129 108 Inclusion compositions calculated to equilibrium with host olivine, wt %

SiO2 TiO2 Al2O3 FeOtot MnO MgO CaO Na2O K2O P2O5

Corr. coeff. Cl S Cu

47.54 51.83 51.05 52.34 51.45 48.3 49.27 50.33 50.32 51.44 0.95 1.98 1.39 1.62 1.48 0.88 0.99 0.81 0.96 1.45 14.23 14.83 16.76 16.04 16.73 13.65 14.54 11.35 13.82 13.51 9.5 10.8 8.53 9.77 9.01 9.54 9.52 9.50 9.51 9.54 0.03 0.22 0.11 0.20 0.18 0.15 0.2 0.17 0.18 0.12 12.18 4.31 4.74 5.04 4.61 10.46 9.7 13.28 10.18 9.32 11.92 8.95 9.64 9.66 9.57 13.54 12.34 11.83 11.67 9.44 2.41 3.57 3.34 3.37 3.56 2.37 2.42 1.74 2.33 3.35 0.89 2.48 1.69 1.88 1.76 0.79 0.70 0.65 0.67 1.30 0.17 0.59 0.39 0.44 0.39 0.17 0.14 0.13 0.18 0.36 99.8 99.5 97.6 100.4 98.7 99.85 99.82 99.79 99.82 99.83 1309 1115 1123 1131 1118 1261 1241 1320 1255 1264 0.85 1.26* 1.10* 1.01 1.10 0.99 1.06 0.96 0.103 0.076 0.063 0.056 0.064 0.141 0.140 0.106 0.127 0.080 0.198 0.017 0.011 0.008 0.022 0.262 0.216 0.168 0.242 0.215 131 342 258 241 227 130 142 107 123 130

S Cu

0.19 125

Total Tcalc (1 atm, dry), °С

48.45 51.79 0.90 1.51 13.01 14.09 9.52 8.20 0.17 0.12 10.72 8.23 10.98 9.84 2.19 3.49 0.63 1.36 0.17 0.38 0.227 0.224 0.119 0.084 2.8 1.7 99.89 101.06 0.84 0.70 17.45 7.24 116 136

S and Cu contents in a primary melt (Fo91)

PETROLOGY

Vol. 26

0.24 118

No. 2

2018

0.19 125

0.16 103

0.21 109

154

MIRONOV, PORTNYAGIN

Table 1. (Contd.) Sample Inclusion

K01-30 K01-30 K01-30 K01-30 K01-30 K01-30 K01-30 K01-30 K01-30 K01-30 7

8-1

8-2

9

10

11

12

14

18-1

19-1

Experiment, Tquench, °C

1235

1192

1192

1178

1161

1185

1166

1230

1207

1186

Experimental time, min

10.50

13.45

13.45

10.00

8.00

13.15

9.30

9.10

14.30

7.15

Inclusion size, а, μm

170

88

110

140

125

162

100

165

90

88

Inclusion size, b

120

125

70

120

160

132

132

140

55

45

Averaged diameter, d

145

107

90

130

143

147

116

153

73

66

Bubble, diameter, μm

48

35

22

40

52

50

44

52

20

20

88.7

88.5

88.5

89.0

88.9

88.6

89.1

88.9

89.3

89.9

1

2

2

1

1

1

1

1

1

1

Host olivine, Fo Melt group (1, 2, 3)

Composition of inclusion glasses, wt % SiO2 TiO2 Al2O3 FeO MnO MgO CaO Na2O K2O P2O5 S Cl H2O Total CaO/Al2O3 CaO/K2O Cu, ppm

46.58 0.87 11.54 9.31 0.15 13.04 11.61 1.88 0.64 0.13 0.199 0.115 2.6 98.69 1.01 18.11 130

48.78 1.13 12.66 8.18 0.14 9.97 9.79 2.92 1.31 0.38 0.185 0.087 3.0 98.55 0.77 7.50 137

49.11 48.96 1.32 0.86 12.63 12.69 8.24 8.64 0.16 0.16 10.07 10.21 11.00 13.08 2.61 2.04 1.31 0.72 0.61 0.15 0.183 0.190 0.077 0.103 2.5 2.6 99.89 100.35 0.87 1.03 8.39 154

18.20 116

47.33 0.94 13.92 8.25 0.16 9.50 12.17 2.71 0.61 0.12 0.328 0.166 2.6 98.83 0.87 19.84 143

47.45 0.92 12.61 8.99 0.19 11.14 12.29 2.22 0.65 0.13 0.210 0.114 2.6 99.55 0.97 18.83 109

47.60 0.96 13.86 8.37 0.14 9.55 12.69 2.38 0.70 0.19 0.269 0.124 2.1 98.92 0.92 18.15 151

47.78 0.82 11.86 9.00 0.14 12.78 11.78 1.92 0.62 0.18 0.243 0.120 2.6 99.86 0.99 19.06 122

49.24 0.90 12.28 8.43 0.14 10.66 12.98 2.20 0.68 0.12 0.293 0.163 n.d. 98.11 1.06 19.23 126

48.10 0.81 12.67 7.71 0.17 11.18 13.43 1.87 0.63 0.10 0.190 0.130 n.d. 97.01 1.06 21.16 131

Inclusion composition calculated to equilibrium with host olivine, wt % SiO2

49.26

50.69

50.05

49.72

48.85

49.16

48.71

49.64

49.88

TiO2

0.98

1.16

1.32

0.85

0.94

0.97

0.95

0.89

0.89

0.8

Al2O3

13.06

12.96

12.65

12.59

13.94

13.29

13.68

12.87

12.16

12.5

FeOtot MnO MgO CaO Na2O

9.52 0.17 10.69 13.14 2.13

9.53 0.14 10.57 10.02 2.99

9.50 0.16 10.56 11.01 2.61

9.52 0.16 11.11 12.98 2.02

9.53 0.16 10.76 12.19 2.71

9.53 0.2 10.55 12.96 2.34

9.53 0.14 11.04 12.53 2.35

9.53 0.15 11.03 12.78 2.08

9.52 0.14 11.4 12.86 2.18

9.51 0.17 12.05 13.25 1.84

K2O

0.72

1.34

1.31

0.71

0.61

0.69

0.69

0.67

0.67

0.62

P2O5 Total Tcalc (1 atm, dry), °С Corr. coeff. Cl S Cu S

49

0.15 0.39 0.61 0.15 0.12 0.14 0.19 0.2 0.12 0.1 99.82 99.79 99.78 99.81 99.81 99.83 99.81 99.84 99.82 99.84 1263 1287 1276 1273 1274 1262 1276 1271 1282 1292 1.12 1.03 1.00 0.99 0.99 1.06 0.99 1.08 0.99 0.98 0.129 0.089 0.077 0.101 0.165 0.121 0.122 0.130 0.162 0.127 0.223 0.190 0.183 0.187 0.326 0.222 0.266 0.263 0.290 0.185 146 140 153 114 142 116 149 132 125 128 0.20

Content of S and Cu in a primary melt (Fo91) 0.17 0.30 0.20 0.24

0.24

PETROLOGY

0.27

Vol. 26

0.17

No. 2

2018

COUPLING OF REDOX CONDITIONS OF MANTLE MELTING AND COPPER...

155

Table 1. (Contd.) Cu

133

Sample

K01-30 K01-30 K01-30 K01-30 K01-30 K01-30 K01-30 K01-30

Inclusion

105

130

105

137

121

23

27

28

31

32

33

35

39

Experiment, Tquench, °C

1180

1166

1213

1167

1183

1155

1206

1215

Experimental time, min

116

121

NT+K01-30 (average)

7.00

10.30

13.45

7.00

8.20

7.00

11.40

13.30

Inclusion size, а, μm

50

55

113

168

75

65

160

64

Inclusion size, b

75

75

83

100

85

70

50

90

Averaged diameter, d

63

65

98

134

80

68

105

77

Bubble, diameter, μm

19

20

31

48

23

20

30

23

89.1

90.4

88.9

88.9

89.4

90.3

88.8

88.6

88.9

85.7

77.1

1

1

1

1

1

1

1

1

1

2

3

48.13

47.91

46.13

47.52

49.74

51.67

Host olivine, Fo Melt group (1, 2, 3) SiO2

Composition of inclusion glasses, wt % 46.73 48.70 43.45 48.56 48.62

TiO2

1.09

1.21

0.81

0.83

1.03

1.15

0.94

0.87

0.98

1.31

1.62

Al2O3

14.81

14.68

12.72

13.18

13.63

15.14

13.09

12.88

14.11

15.48

16.09

FeO MnO MgO CaO Na2O

8.14 0.12 8.99 13.38

6.91 0.13 9.71 12.96

9.02 0.15 11.28 12.91

8.94 0.18 9.97 12.80

7.97 0.14 9.84 13.76

7.25 0.12 8.93 13.04

9.25 0.17 10.58 13.16

8.95 0.14 10.53 13.29

8.02 0.14 9.08 13.09

8.20 0.14 6.59 10.28

9.52 0.18 4.68 9.46

2.47

2.61

2.00

2.12

2.27

2.55

2.04

2.18

2.42

3.16

3.46

K2O

1.02

0.63

0.67

0.73

0.78

0.90

0.67

0.72

0.80

1.43

1.95

P2O5

0.25

0.19

0.16

0.20

0.20

0.21

0.13

0.18

0.17

0.37

0.45

S Cl H2O

0.243 0.129

0.261 0.161

0.214 0.140

0.235 0.133

0.208 0.152

0.189 0.110

0.237 0.128

0.225 0.153

0.228 0.129

0.138 0.074

0.014 0.065

n.d.

n.d.

n.d.

n.d.

n.d.

n.d.

n.d.

n.d.

2.5

2.4

Total CaO/Al2O3

98.80

97.37

96.82

98.02

93.45

98.17

99.04

96.29

97.5

97.5

0.90

0.88

1.01

0.97

1.01

0.86

1.01

1.03

CaO/K2O

13.16

20.71

19.21

17.53

17.67

14.56

19.67

18.53

Cu, ppm SiO2

169

136

129

119

155

141

135

129

Inclusion composition calculated to equilibrium with host olivine, wt% 48.06 48.02 48.49 49.47 46.28 48.16 49.12 48.08

TiO2

1.03

1.11

0.86

Al2O3

14.02

13.5

13.44

FeOtot

9.52

9.53

9.51

MnO MgO CaO Na2O

0.11 10.86 12.66

0.12 12.48 11.92

2.34

K2O

137 48.56

0.67

0.6

7.3

5.0

175 50.56

244 51.67

1.08

1.03

0.95

0.92

0.96

1.28

1.62

14.26

13.61

13.18

13.59

13.73

15.05

16.09

9.52

9.53

9.52

9.50

9.50

9.52

9.53

9.52

0.16 10.74 13.64

0.18 10.99 12.72

0.15 10.72 14.4

0.11 12.38 11.72

0.17 10.79 13.25

0.15 10.28 14.03

0.14 10.85 12.77

0.13 8.47 9.99

0.18 4.68 9.46

2.4

2.11

2.11

2.38

2.29

2.05

2.3

2.35

3.07

3.46

0.97

0.58

0.71

0.73

0.82

0.81

0.67

0.76

0.77

1.39

1.95

P2O5

0.24

0.17

0.17

0.2

0.21

0.19

0.13

0.19

0.17

0.36

0.45

Total

99.81

99.83

99.83

99.84

99.83

99.82

99.81

99.80

Tcalc (1 atm, dry), °С

1274

1312

1263

0.82

0.93 16.8

99.1

13.1

1271

1268

1310

1264

1254

Corr. coeff.

0.95

0.93

1.06

1.00

1.05

0.90

1.00

1.06

Cl

0.123

0.149

0.148

0.133

0.160

0.100

0.128

0.162

PETROLOGY

Vol. 26

No. 2

2018

99.8 1271 0.127

99.8 1233 0.072

99.1 1122 0.065

156

MIRONOV, PORTNYAGIN

Table 1. (Contd.) S

0.232

0.241

Cu

161

126

S Cu

0.21 148

0.23 121

0.226

0.235

0.218

0.171

0.237

136 119 164 128 135 S and Cu contents in a primary melt (Fo91) 0.21 125

0.21 109

0.20 153

0.16 122

0.22 123

0.238

0.223

136

133

0.22 124

0.204 122

0.135 170

0.014 267

Here shown are compositions of melt inclusions (MI) for two studied TD samples: sample of 1975 Northern breakthrough (NT) and cone “1004” (K01-30). Inclusions from sample NT are represented by naturally quenched MI. Inclusions from sample K01-30 were reheated and quenched, with shown heating time (exposure time from 1064°C to the quenching temperature, T = 1064°C was controlled by the moment of Au melting) in minutes and seconds (for instance, 7.10, i.e. 7 min and 10 s). Melt groups: (1) Primitive melts unaffected or slightly affected by mixing with differentiated subalkaline melts. The group was distinguished based on K2O < 1 wt % (for compositions corrected for equilibrium with host olivine) (Fig. 1a); (2) Hybrid melt group showing signs of mixing; (3) groundmass glass. MI compositions were corrected for equilibrium with host olivine with allowance for iron loss using an “Fe-loss” software (Danyushevsky et al., 2000); FeOinitial = 9.5 wt % (corresponds to the average content in high-Mg rocks of TD), Fe2+/Fe3+ in a melt was given as constant value equal 5.5 (FeO = 8.05, Fe2O3 = 1.63 wt %), estimated from study of syngenetic spinel inclusions in olivine (Table 4). Primary melts were recalculated to equilibrium with olivine Fo91 (only for group of primitive inclusions) using Petrolog3 software (Danyushevsky and Plechov, 2011) by reverse olivine fractionation at P = 10 kbar and fO2 = QFM+1. The Fo values for the groundmass were calculated using a Petrolog3 software at 1 atm and fO2 = NNO. Correction coefficients (Corr. coeff.) were used to correct sulfur and copper contents in melts calculated to equilibrium with their olivine. The coefficients were calculated from ratios: K2O in a melt (value corrected for equilibrium with olivine)/K2O in an inclusion glass (measured value). Correction coefficients for calculation of S and Cu contents in primary melts from equilibrium composition of melt inclusions (not listed in table) varied from 0.87 to 0.97, 0.92 on average. * In name (3gm* and 4gm*) is the correction for copper content in the groundmass glass due to trapping the Pl–Px microlites during laser analysis; uncorrected values for them are shown in Fig. 2a; n.a. for H2O means not analyzed. Oxides, S, and Cl are given in wt %, Fo, in mol %, Cu in ppm.

from the Siqueiros Transform Fault in the East Pacific Rise (Fig. 2b; Danyushevsky et al., 2003). DISCUSSION Model of the Cu and S Behavior during Mantle Melting Variations of Cu and S contents in the primary magmas of Tolbachinsky Dol as compared to MORB can be explained using mantle melting model of Lee et al. (2012). This model made it possible to predict the Cu and S contents in the melts derived by different degrees of peridotite melting at variable oxygen fugacity, Cu and S contents in a source, Cu sulfide–melt partition coefficient, pressure and temperature of melting. The model of Lee et al. (2012) is underlain by (1) model of Mavrogenes and O’Neill (1999) of S2– solubility in a melt depending on Т and Р, and (2) model of Jugo et al. (2010) describing the total sulfur content in melt as a function of oxygen fugacity and accounting for a sharp increase of S6+ solubility and, respectively, the total sulfur content within the range of QFM + 1 – QFM + 1.5. It is unclear, however, whether these models provide adequate description of sulfur solubility depending on P-T conditions and melt composition (for instance, depending on water or nickel content in melt) (Ariskin et al., 2013). The accuracy of the Jugo et al. (2010) model is also an open question. For instance, the dependence of S solubility on ΔQFM implied by Jugo et al. (2010) and, respectively, by Lee et al. (2012) suggests that log S6+/S2– iso-

pleths in the 1/T–log fO2 field are parallel to QFM equilibrium line, which is not likely the case (A.A. Borisov, personal communication). The pressure effect on the sulfur solubility is also uncertain (Jugo et al., 2010; Matjuschkin et al., 2016). In this work, we are based on the model (Lee et al., 2012), because no alternative mutually consistent model has been proposed yet that would take into account discovered shortcomings. Estimated absolute fO2 values are certainly model dependent, but the relative difference in the redox conditions of magma generation for TD and MORB is not. The difference in fO2 may increase if the higher pressure formation of the TD magmas (2 GPa) as compared to MORB magmas (1 GPa) is taken into account (Matjuschkin et al., 2016). The Cu content in the mantle sources of TD and MORB magmas was taken to be 30 ppm, which corresponds to the DMM composition (Salters and Stracke, 2004); the sulfur content—200 ppm (Lee et al., 2012). It should be noted that the estimated sulfur concentrations in the upper mantle varies from 119 ppm (Salters and Stracke, 2004), 146 ± 35 ppm (Saal et al., 2002) and up to 190 ± 90 ppm (Lorand, 1990). The sulfur content of 200 ppm taken in this work corresponds to the upper limit of the inferred contents. Such S content in a source allows us to closely reproduce the Cu and S contents in typical MORB and to compare our modeling results with those in Lee et al. (2012). The pressure of the primary magma generation for TD was estimated by Portnyagin et al. (2015) to PETROLOGY

Vol. 26

No. 2

2018

COUPLING OF REDOX CONDITIONS OF MANTLE MELTING AND COPPER...

157

Table 2. Average compositions of primitive and primary melts for the studied samples of Tolbachinsky Dol (TD) and Pacific basalts of the Juan de Fuca (JDF) ridge and estimated conditions of crystallization and formation Object (sample) TD (NT) TD (K01-30) MORB (JDF) Melts I st. dev. II I st. dev. II I st. dev. II (I—primitive, II—primary) Equilibrium olivine, Fo 88.6 1.0 91 89.1 0.7 91 85 89.4 Melt composition SiO2, wt % 48.07 0.52 47.34 48.91 0.95 48.2 49.55 0.78 48.77 TiO2 0.99 0.04 0.9 0.94 0.09 0.86 1.33 0.24 0.85 Al2O3 14.45 0.41 13.14 13.31 0.75 12.24 15.99 0.74 18.69 FeOtot 9.52 0.01 9.57 9.52 0.01 9.55 9.45 0.82 7.37 MnO 0.11 0.03 0.1 0.16 0.03 0.15 MgO 10.40 1.00 13.81 11.07 0.86 14.00 8.43 0.34 10.0 CaO 12.68 0.26 11.53 12.82 0.75 11.79 11.90 0.51 11.5 Na2O 2.57 0.12 2.34 2.22 0.22 2.04 2.55 0.27 2.43 K2O 0.87 0.06 0.79 0.71 0.09 0.65 0.13 0.05 0.08 P2O5 0.18 0.02 0.16 0.16 0.04 0.15 0.14 0.04 0.090 Total 99.83 0.01 99.7 99.8 0.01 99.6 99.46 0.40 99.78 CaO/Al2O3 0.88 0.03 0.88 0.97 0.07 0.96 0.75 0.05 0.62 Cl 0.115 0.007 0.105 0.134 0.020 0.123 0.010 0.002 S 0.207 0.026 0.188 0.232 0.039 0.212 0.106 0.010 0.106 (n = 19) S max 0.274 0.298 Cu, ppm 133 15 121 133 15 121 74 12 74 Number of analyses, n 12 20 44 P-T crystallization conditions P, GPa 1 1 0.5 T, °С (dry, 1 atm) 1206 T, °С (dry, 0.5 GPa) 1270 T, °С (dry, 1 GPa) 1314 25 1392 1324 20 1390 1295 1205 1283 1215 1281 T, °С (4 wt % H2O, 1 GPa) P-T-F formation conditions 2.0 1360

P, GPa TDPS,°С (Zhang and Hirschmann, 2016) H2O, wt % dTDPS,°С (Portnyagin et al., 2007a) T, °С (TDPS+dTDPS) Fmelt, % P, GPa (Lee et al., 2009) T, °C (Lee et al., 2009)

2

4 –75 1285 7 2.0 1332

1

2.0 1360

1.0 1240

4 –60 1300 10 1.8 1325

0 60 1300 8 1.2 1308

2

Table shows the average compositions of primitive melts (group 1 in Table 1) of TD and standard deviation from mean. For primitive MORB composition of the JDF ridge, we took the average composition of high-Mg basalts with MgO > 8 wt % according to (Jenner and O’Neill, 2012, except for data on seamount basalts) and (Gale et al., 2013). Primary melts were calculated for TD by reverse fractionation of olivine to equilibrium with Fo91, at 10 kbar, QFM + 1 using a Petrolog3 software (Danyushevsky and Plechov, 2011). For MORB compositions of JDF, the calculations were carried out by reverse fractionation of Ol (10%) and Pl (25%) up to MgO = 10 wt % at 5 kbar and QFM using the same program. Obtained composition is close to the primary melt TOR-2 (Sobolev et al., 1989). Crystallization temperature was estimated with allowance for pressure (+5 С/kbar) and influence of water content (according to the model of Almeev et al., 2007 T decrease for 4 wt % H2O is 110°C)). Temperature of formation (or melting) was calculated using formula TDPS + dTDPS, where TDPS is the temperature of dry peridotite solidus at P = 2 GPa and 1 GPa (Zhang and Hirschmann, 2016, Fig. 10), while dTDPS is the deviation of true melting temperature from TDPS (dTDPS, ~ –75–60°C for TD and +60°С for MORB). The dTDPS values were estimated according to (Portnyagin et al., 2007a, Fig. 7b and H2O–Fmelt relation in Supplementary there), using estimated degrees of melting for TD and JDF and water content in primary magmas of ~4 wt % and 0 wt %, respectively. In all three cases (NT, K01-30, JDF), melting temperatures are close to 1300°С. Melting degree, Fmelt, for TD was estimated using technique (Portnyagin et al., 2007а) on the basis of trace element contents in the average compositions of initial melts. Fmelt for JDF (MORB) was estimated from Fmelt– Na2O relation (Langmuir, 1992, Fig. 42 where Na2O ~ 2.1–2.7 wt % correspond to Fmelt ~ 10–6%, 8%, on average). The pressures and temperatures of the formation of the same primary melts according to model (Lee et al., 2009) are shown for comparison. PETROLOGY

Vol. 26

No. 2

2018

158

MIRONOV, PORTNYAGIN

Table 3. Composition of crystalline inclusions Cpx and Opx in olivine from Tolbachinsky Dol (TD), cone “1004” TD, cone “1004” (sample K01-30) Components 22

23

26

24

25a

25b

25c

27a

27b

Inclusion

Cpx

Cpx

Cpx

Opx

Opx

Opx

Opx

Opx

Opx

SiO2

53.91

53.10

53.48

55.89

56.80

56.77

55.32

57.57

56.45

TiO2

0.15

0.16

0.19

0.08

0.01

0.06

0.03

0.14

0.14

Al2O3

1.02

1.52

1.47

0.65

0.98

0.99

0.80

0.87

0.68

FeOtot MnO MgO CaO Na2O

3.45 0.07 18.82 21.70 0.24

3.60 0.10 17.72 22.04 0.25

3.67 0.08 18.14 22.16 0.25

6.77 0.12 33.30 1.90 0.01

6.43 0.16 33.42 1.79 0.00

6.27 0.18 33.47 1.80 0.00

6.71 0.14 34.81 1.69 0.00

6.66 0.19 33.40 2.03 0.01

6.84 0.18 33.90 1.97 0.03

K2O

0.01

0.01

0.00

0.00

0.01

0.01

0.01

0.01

0.01

0.69 0.07 100.1 90.7 1185 –2.2

0.79 0.01 99.3 89.8 1266 7.7

0.72 0.04 100.2 89.8 1245 5.2

0.29 0.14 99.2 89.8

0.41 0.09 100.1 90.3

0.45 0.05 100.1 90.5

0.43 0.10 100.0 90.2

0.43 0.02 101.3 89.9

0.31 0.05 100.6 89.8

Ol

Ol

Ol

Ol

Ol

Ol

38.70 10.20 0.19 48.63 0.19 0.27 0.04 98.2 89.5

40.09 9.77 0.17 48.94 0.16 0.36 0.05 99.5 89.9

39.98 9.80 0.19 48.92 0.15 0.31 0.04 99.4 89.9

40.08 9.76 0.19 49.16 0.16 0.34 0.04 99.7 90.0

Cr2O3 NiO Total Px, #Mg T, °С (Putirka et al., 2003) P, kbar (Putirka et al., 2003) Host mineral

Ol

Ol

SiO2

39.92 9.92 0.16 49.15 0.20 0.22 0.04 99.6 89.8

39.92 10.79 0.20 48.43 0.19 0.23 0.03 99.8 88.9

FeO MnO MgO CaO NiO Cr2O3 Total Fo

Ol 40.48 11.43 0.20 48.14 0.18 0.22 0.05 100.7 88.2

40.88 10.27 0.18 49.32 0.17 0.27 0.04 101.1 89.5

41.05 10.17 0.20 49.32 0.18 0.27 0.05 101.2 89.6

Pressure and temperature of crystallization were calculated using Cpx–melt model (Putirka et al., 2003); the average primitive composition for sample K01-30 (Table 2) estimated from melt inclusion data in olivine was used as a melt. For unclear reasons, one Cpx inclu+ sion (K01-30-22) yielded negative pressure and lowered temperature. Px, #Mg = Mg/(Mg + Fe2tot ).

be ~2 GPa. The formation pressure of the primary magmas of the JDF ridge was assumed to be ~1 GPa on the basis of previous estimates made for TOR-2 basalts (Dmitriev et al., 2006). Estimated melting temperature of a source of TD rocks and basalts of the JDF ridge was ~1300°C (Table 2). This value is well consistent with the liquidus temperature of crystallization of the initial magmas of the TD (1280–1285°C at 1 GPa and 4 wt % H2O) and of the JDF ridge (1270– 1295°C at 0.5–1 GPa) (Table 2). Similar pressures and temperatures of primary magma generation (2.0–1.8 GPa and 1330–1325°C for TD; ~1.2 GPa and 1310°С for JDF ridge) are also predicted by the model (Lee et al., 2009) (Table 2).

Calculated Cu and S contents in melts are dependent on the pressure and temperature variations owing to their effect on sulfur solubility in a melt. In particular, variations ±50°C or ±0.5 GPa have approximately the same effect as variations ±0.05–0.1 log units ΔQFM in calculations according to model (Lee et al., 2012). A temperature increase by 50°С leads to the increase of S solubility by ~0.03 wt % and is equivalent to a decrease of pressure by 0.5 GPa or increase of ΔQFM by 0.05–0.1 log units. Thus the existing uncertainties of the estimated temperature and pressure of magma formation do not affect the main conclusions of this work. Based on data (Zhang and Hirschmann, 2016), the temperature of sulfide liquidus at a pressure of 2 GPa PETROLOGY

Vol. 26

No. 2

2018

COUPLING OF REDOX CONDITIONS OF MANTLE MELTING AND COPPER...

is 1236°C. This denotes that at temperature of mantle melting of 1300°C, sulfide exists as sulfide liquid (SL). With allowance for this fact, we calculated the Cu partition coefficient between sulfide liquid (SL) and siliSL-SM cate melts (SM) (DCu ) using model (Li and Audetat, 2015). At a given temperature of 1300°C and ΔQFM = +1, FeO = 9.5 wt % in a silicate melt for TD and ΔQFM = 0 and FeO = 7.5–8 wt % for the JDF SL-SM Ridge, DCu were ~600 (TD) and ~700 (JDF), respectively. The calculated values are the least sensitive to variations of given ΔQFM (Li and Audetat, SL-SM 2015). The use of the lower DCu than in the work of SL-SM Lee et al. (2012) ( DCu = 800) results in the increase of the calculated copper contents in primary melts by 10–20 ppm.

Assessment of Conditions of Magma Formation on the Basis of Cu–S Systematics Our modeling makes it possible to compare the observed Cu and S contents in primary magmas with their model contents at the estimated temperature and pressure of magma formation, to quantify independently the degree and redox conditions of mantle melting beneath TD, and compare the estimated parameters with the conditions of MORB formation (Fig. 3). As demonstrated by the numerical modeling, the high observed contents of Cu (100–150 ppm, on average 120 ppm) and S (0.16–0.30 wt %, on average 0.21 wt %) in the initial TD melts correspond to the degree of mantle melting from ~6 to 12% and ΔQFM ranging from +1 to +1.3 (Fig. 3). Obtained degrees of mantle melting are close to those estimated independently from the average MIs trace-element compositions (5– 11%, Table 2) and for TD magmas in general (~6%, Portnyagin et al., 2015). However, it should be noted that independent estimates of degree of melting are required for each of the individual melt inclusions for a more accurate comparison of modeling results. The modeling results can be also used to estimate the possible sulfur influx in a mantle source and sulfur content in a source. Obtained data show that no sulfur influx is required to explain the average sulfur composition of primary magmas of the Northern Breakthrough. The maximum sulfur contents in these melts can be explained by 7% melting of MORB-type mantle source under oxidized conditions (ΔQFM=+1.25) (Fig. 3). The average composition of the initial magmas of Cone “1004” also does not require an additional sulfur influx and can be obtained by 9% melting at ΔQFM from +1.1 to +1.2. The maximum sulfur contents (0.30 wt %) in the primary magmas of Cone “1004” also can be obtained by melting of mantle with up to 200 ppm S, but at sufficiently low degrees of melting (~5–7%), which are lower than estimated for this sample from trace elements (9–11%, Table 2). PETROLOGY

Vol. 26

No. 2

2018

159

Alternatively, the observed sulfur and copper contents in MIs from Cone “1004” can be obtained at ΔQFM = +1.28 by higher degrees of melting Fmelt ~10–11% of a more S-rich source (300 ppm) than DMM (Fig. 3), which may imply up to 30–35% sulfur in primary melt derived from the subducted plate. The redox conditions of the crystallization of the studied primitive TD melts were also estimated using oxybarometer based on the dependence of composition of the Ol–Cr-Spl–(Opx) equilibrium assemblage on fO2 (Ballhaus et al., 1991; Nikolaev et al., 2016). An application of this oxybarometer for TD is well justified, because the primitive melts are in equilibrium with this assemblage comprising co-genetic crystalline Cr-spinel as well as orthopyroxene inclusions in highMg olivine phenocrysts (Tables 3, 4). The calculated ΔQFM values for olivine Fo88.2-90.7 according to the model (Ballhaus et al., 1991) (at P = 1 atm) vary within the range of +1.5…+2 (on average +1.8) regardless of olivine composition. There is no difference between two samples studied (Fig. 4, Table 4). The values of ΔQFM obtained using modified model (Nikolaev et al., 2016) are systematically by ~0.4 log units lower (ΔQFM +1.2…+1.6, on average +1.4) (Fig. 4, Table 4). The estimates obtained from modeling of Cu–S systematics (+1.1 < ΔQFM +1) and FeOtot = 9–11 wt % (Portnyagin et al., 2012, Fig. 10). The persistent or increasing oxidation during magma fractionation presumably also provides favorable conditions for postmagmatic abundant copper sulfate and chloride sulfide-free PETROLOGY

Vol. 26

No. 2

2018

COUPLING OF REDOX CONDITIONS OF MANTLE MELTING AND COPPER...

mineralization of exhalations, which are typical of TD (Pekov et al., 2014; I. Pekov, personal communication, 2016). Contribution in the Discussion of Conditions of the Island-Arc Magma Formation The Cu and S contents in primary melts are insensitive to the possible oxidation of magmas after their extraction from a mantle source. Oxidation can only retard the melt saturation in sulfide phase and prior to this moment, Cu and S will behave like elements strongly incompatible in silicate minerals (for instance, Ba, La, Nb). Therefore, our conclusion concerning the similarity of the estimated redox conditions during melting of island-arc mantle (Cu–S systematics) and crystallization of primitive melts (olivine–spinel oxybarometry) has an important petrogenetic significance. Firstly, this indicates that the redox state of magmas can be preserved during ascent of primary magmas from the mantle and their further crystallization at crustal conditions (P ~5–8 kbar, see above, and Table 3). Therefore olivine–spinel oxybarometry can be informative about the conditions of magma generation in the mantle (Fig. 4; Evans et al., 2012). This conclusion casts some doubts on the widely discussed hypothesis about significant oxidation of mantle magmas en route to the surface (Lee et al., 2012) Secondly, obtained data raise questions regarding the main conclusion of Lee et al. (2012) on similarity of the redox conditions during formation of midocean ridge basalts (MORB) and island-arc ( subduction-related) magmas based on the close initial content of only copper (50–100 ppm), without consideration of sulfur content data. In contrast to (Lee et al., 2012), our work involves data on sulfur contents in primitive melts, which in combination with Cu content unambiguously indicate the more oxidized conditions of the formation of Tolbachinsky island-arc magmas as compared to those of MORB. At the same time, the moderate oxidation of the island-arc mantle and formation of primary magmas in the presence of sulfide phase place significant restriction on the possible enrichment of these magmas in copper and other strongly chalcophile elements (Mungall, 2002). Future studies of primitive melt inclusions on regional scale would make it possible to analyze the Cu–S systematics for different island-arc complexes and to provide alternative estimate of redox conditions of the formation of the primary island-arc magmas. Data obtained in this work also demonstrate that the observed Cu and S contents in the initial TD magmas can be generated from MORB-like source without S enrichment, which has been suggested for suprasubduction zones (e.g., Richards, 2015 and references therein). The elevated sulfur content in the TD melts was likely caused by the increase of sulfide solubility at PETROLOGY

Vol. 26

No. 2

2018

163

a moderate increase of fO2 (Jugo et al., 2010). The missing or very limited subduction-related sulfur influx to the initial magmas of the Northern Breakthrough of Tolbachinsky Dol cast doubts concerning the key role of sulfate fluids in the oxidation of mantle sources of island-arc magmas (e.g., Jégo and Dasgupta, 2014). In this case, the more probable oxidizing agent in the mantle can be Fe2O3 dissolved in slab-derived melts (Mungall, 2002), as is also suggested by recent data on TD emphasizing a major role of slab-derived hydrous melts in the enrichment of its mantle source (Portnyagin et al., 2015). CONCLUSIONS (1) The study of melt inclusions in olivine (Fo83.1from the rocks of Tolbachinsky Dol (TD) (Kamchatka) showed that the most primitive melts have higher contents of Cu (100–150 ppm) and S (0.16– 0.30 wt %) as compared to typical MORB, for instance, from the Juan de Fuca (JDF) ridge (Pacific Ocean), which are characterized by the moderate contents of Cu (55–105 ppm) and S (0.09–0.12 wt %) in the parental melts. (2) Modeling algorithm of Lee et al. (2012) showed that the primary melts of TD and JDF ridge could be derived from a similar mantle source with S ~200 ppm and Cu ~30 ppm, but at different oxygen fugacity: +1.1