This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Unit U U U U U V V V W W W X1SG 2 2 2 2 2 2 2 2 2 2 2 3Sample P2150A P2150B P2150C P2151 P2152 P2161A P2162 P2162A P2155A P2155B P2156 P2157
Age BP 2850 2850 2850 2850 2850 2800 2800 2800 2750 2750 2750 2150
Lu 0.548 0.535 0.538 0.480 0.512 0.490 0.507 0.552
Hf 5.82 5.66 5.83 4.97 5.05 5.17 5.11 6.33
Pb 22.7 22.4 22.7 15.5 16.7 16.6 16.6 22.9
Th 11.25 10.98 10.99 9.61 10.27 10.24 10.42 11.20
U 2.50 2.44 2.44 2.18 2.32 2.27 2.31 2.44
Rb/Sr 0.74 0.72 0.73 0.79 0.83 0.82 0.89 0.60
Eu/Eu* 0.69 0.68 0.70 0.66 0.65 0.66 0.65 0.68
Oxide abundances normalised to 100 % on a volatile free basis, with original analytical totals and LOI (loss on ignition) values given. Eruption ages are given in years before present (BP) and sourced from Wilson (1993) and Hogg et al. (2012). Eu/Eu* calculated as Eu / √(Sm*Gd), normalised to chondrite. Major elements analysed by X-Ray Flourescence (XRF) and trace elements by solution Inductively Coupled Plasma Mass Spectrometry (ICPMS). See Barker et al. (2015) for full description of analytical methods and standards.
Table DR1. Major and trace element compositions of whole rock (WR) pumices and lava from the final SG2 eruptions (U, V and W) and the SG3 eruptions (X,Y and Z)
Barker et al. 2016, Data Repository
GSA Data Repository 2016100Rapid priming, accumulation, and recharge of magma driving recent eruptions at a hyperactive caldera volcano
Total 99.6 94.1-100.9 100.1 94.6-97.9 98.5 93.1-98.7 100.6
n 23 21 23 15 14 16 18
Glass trace elements (ppm)Li 41.26 38.86 38.59 40.55
Sc 10.2 10.5 10.4 11.2
V 0 0 1 1
Cr 0.1 0.5 0.3 0.3
Ni bdl 0.0 0.1 bdl
Cu 2.9 2.6 4.3 3.5
Zn 63 64 63 68
Ga 15.3 16.3 15.7 16.3
Rb 107.4 112.3 112.4 104.8
Sr 105 104 112 151
Y 34.0 36.9 35.9 36.1
Zr 215 224 219 244
Nb 7.9 2.1 8.7 2.8
Cs 5.26 5.52 5.39 5.10
Ba 637 698 681 713
La 25.9 28.1 27.5 27.3
Ce 56.6 60.2 58.9 59.2
Pr 6.59 7.25 7.22 7.15
Nd 26.2 28.4 28.2 27.0
Sm 5.66 6.10 5.93 5.99
Eu 1.11 1.19 1.21 1.29
Gd 5.59 6.39 6.16 6.28
Tb 0.89 0.97 0.95 0.93
Dy 5.59 6.19 6.16 6.01
Ho 1.21 1.31 1.27 1.25
Er 3.61 3.89 3.77 3.82
Tm 0.534 0.557 0.554 0.559
Yb 3.66 3.95 3.73 3.81
Lu 0.549 0.588 0.580 0.588
Hf 5.95 6.18 6.15 6.28
Pb 21.1 22.1 22.1 22.0
Th 10.90 11.52 11.21 10.79
U 2.43 2.63 2.67 2.60
Rb/Sr 1.02 1.08 1.01 0.69
Eu/Eu* 0.60 0.58 0.61 0.65
Temp 816 ˚C 813 ˚C 813 ˚C 849 ˚C
log η 5.17 5.17 5.32 4.72
Table DR2. Major and trace element compositions of groundmass (GM) glass and major element compositional range of melt inclusions (MI) from the final SG2 eruptions (U, V and W) and the SG3 eruptions (X,Y and Z) and subunits
Oxide abundances normalised to 100 % on a volatile free basis, with original analytical totals given. Major elements measured by Electron Probe Micro-Analysis (EPMA), trace elements measured by solution-ICPMS for pure groundmass glass separates only. Temperatures are averages from Barker et al. (2015) using the orthopyroxene-liquid thermometer of Putirka (2008). Viscocity (log η) of melt (in Pa s ) calculated using the model of Giordano et al. (2008) assuming 4.5 wt. % H2O content. Other details as in Table DR1. See Barker et al. (2015) for analytical methods and standards.
Barker et al. 2016, Data Repository
Unit Y1 Y3 Y7 ZSample P2119 P2296 P2287 P2173A
Type MI glass GM glass MI glass GM glass MI glass GM glass MI glass GM glass
Glass major elements (wt. %)SiO2 75.34-76.93 74.80 ±0.62 74.18-76.17 75.03 ±0.55 74.71-77.00 75.24 ±0.30 75.17-76.83 76.80 ±0.63
Total 94.0-99.7 100.4 93.9-96.9 100.13 94.2-100.2 98 94.5-96.4 98.7
n 12 8 14 11 21 10 28 17
Glass trace elements (ppm)Li 44.84 42.56 36.66 42.45
Sc 10.6 10.8 10.6 10.8
V 1 1 1 2
Cr 0.3 0.4 1.2 0.2
Ni 0.1 0.1 0.1 0.0
Cu 5.2 2.4 2.0 1.9
Zn 71 69 68 68
Ga 16.3 16.2 15.6 15.6
Rb 105.2 101.1 93.9 104.6
Sr 148 152 145 131
Y 35.0 34.3 32.6 35.1
Zr 238 231 217 238
Nb 3.9 5.5 4.4 5.3
Cs 5.77 5.64 5.29 5.70
Ba 661 644 614 651
La 27.4 26.7 25.2 27.1
Ce 59.5 57.8 54.4 58.5
Pr 6.91 6.74 6.48 6.76
Nd 27.4 26.7 25.3 27.1
Sm 5.96 5.86 5.45 5.80
Eu 1.25 1.34 1.23 1.23
Gd 6.13 6.06 5.66 6.10
Tb 0.98 0.94 0.86 0.93
Dy 6.07 5.96 5.46 5.96
Ho 1.27 1.23 1.14 1.25
Er 3.82 3.73 3.42 3.76
Tm 0.585 0.550 0.517 0.565
Yb 3.75 3.71 3.38 3.72
Lu 0.591 0.571 0.513 0.577
Hf 6.59 6.30 5.74 6.42
Pb 21.9 21.7 19.4 21.5
Th 11.12 11.05 9.91 11.07
U 2.62 2.52 2.28 2.52
Rb/Sr 0.71 0.67 0.65 0.79
Eu/Eu* 0.63 0.69 0.68 0.63
Temp 851 ˚C 852 ˚C 853 ˚C 853 ˚C
Viscocity 4.67 4.73 4.76 4.78
Table DR2. continued
Barker et al. 2016, Data Repository
Unit U 1 mm
plain light BSE (mirror image)
0
2
4
6
8
10
12
14
16
18
36 4038 42 44 46 48 50 52 54 56 58 60 62 64 66
Key
unzoned
outer rim
interior
core
no. o
f ana
lyse
s
Mg-number
A
B
normal zoning
Fig. DR1. (A) Orthopyroxene textures from Unit U (SG2) in plain light photograph (left) and corresponding Back-Scattered Electron (BSE) images (right). Dark shades in BSE zoning are relatively low FeO or high MgO, light shades are high FeO or low MgO (B) Histogram showing orthopyroxene compositions for zones within single crystals from A, analysed by EPMA. Mg-number is calculated as Mg= 100*Mg/Mg+Fetotal (modified from Barker et al., 2015)
Barker et al. 2016, Data Repository
Unit V1 mm
0
2
4
6
8
10
12
36 4038 42 44 46 48 50 52 54 56 58 60 62 64 66
Key
unzoned
outer rim
interior
core
no. o
f ana
lyse
s
Mg-number
A
B
Fig. DR2. (A) Orthopyroxene textures from Unit V (SG2) in plain light photograph (left) and corresponding BSE images (right). (B) Histogram showing orthopyroxene compositions for zones within single crystals from A, analysed by EPMA. Other details as in Fig. DR1.
normal zoning
Barker et al. 2016, Data Repository
Unit W1 mm
0
2
4
6
8
10
12
36 4038 42 44 46 48 50 52 54 56 58 60 62 64 66
Key
unzoned
outer rim
interior
core
no. o
f ana
lyse
s
Mg-number
A
B
Fig. DR3. (A) Orthopyroxene textures from Unit W (SG2) in plain light photograph (left) and corresponding BSE images (right). (B) Histogram showing orthopyroxene compositions for zones within single crystals from A, analysed by EPMA. Other details as in Fig. DR1.
normal zoning
Barker et al. 2016, Data Repository
Unit X 1 mm
0
2
4
6
8
10
14
36 4038 42 44 46 48 50 52 54 56 58 60 62 64 66
Key
unzoned
outer rim
interior
core
no. o
f ana
lyse
s
Mg-number
B
Fig. DR4. (A) Orthopyroxene textures from Unit X (SG3) in plain light photograph (left) and corresponding BSE images (right). (B) Histogram showing orthopyroxene compositions for zones within single crystals from A, analysed by EPMA. Cross-hatched area refers to dominant compositional mode for SG2 orthopyroxenes from Figs DR1-DR3. Other details as in Fig. DR1.
normal zoning
A
reverse zoning
SG2 inheritedcores
12
Barker et al. 2016, Data Repository
Unit Y1 1 mm
0
5
10
15
20
25
30
36 4038 42 44 46 48 50 52 54 56 58 60 62 64 66
Key
unzoned
outer rim
interior
core
no. o
f ana
lyse
s
Mg-number
B
Fig. DR5. (A) Orthopyroxene textures from Unit Y1 (SG3) in plain light photograph (left) and corresponding BSE images (right). (B) Histogram showing orthopyroxene compositions for zones within single crystals from A, analysed by EPMA. Cross-hatched area refers to dominant compositional mode for SG2 orthopyroxenes from Figs DR1-DR3. Other details as in Fig. DR1.
normal zoning
A
reverse zoning
SG2 inheritedcores
Barker et al. 2016, Data Repository
Unit Y3 1 mm
0
2
4
6
8
10
12
36 4038 42 44 46 48 50 52 54 56 58 60 62 64 66
Key
unzoned
outer rim
interior
core
no. o
f ana
lyse
s
Mg-number
B
Fig. DR6. (A) Orthopyroxene textures from Unit Y3 (SG3) in plain light photograph (left) and corresponding BSE images (right). (B) Histogram showing orthopyroxene compositions for zones within single crystals from A, analysed by EPMA. Cross-hatched area refers to dominant compositional mode for SG2 orthopyroxenes from Figs DR1-DR3. Other details as in Fig. DR1.
normal zoning
A
reverse zoning
SG2 inheritedcores
Barker et al. 2016, Data Repository
Unit Y7 1 mm
0
2
4
6
8
36 4038 42 44 46 48 50 52 54 56 58 60 62 64 66
Key
unzoned
outer rim
interior
core
no. o
f ana
lyse
s
Mg-number
B
Fig. DR7. (A) Orthopyroxene textures from Unit Y7 (SG3) in plain light photograph (left) and corresponding BSE images (right). (B) Histogram showing orthopyroxene compositions for zones within single crystals from A, analysed by EPMA. Cross-hatched area refers to dominant compositional mode for SG2 orthopyroxenes from Figs DR1-DR3. Other details as in Fig. DR1.
normal zoning
A
reverse zoning
SG2 inheritedcores
10
Barker et al. 2016, Data Repository
eruption Z 1 mm
0
2
6
10
14
36 4038 42 44 46 48 50 52 54 56 58 60 62 64 66
Key
unzoned
outer rim
interior
core
no. o
f ana
lyse
s
Mg-number
B
Fig. DR8. (A) Orthopyroxene textures from eruption Z (SG3) in plain light photograph (left) and corresponding BSE images (right). (B) Histogram showing orthopyroxene compositions for zones within single crystals from A, analysed by EPMA. Cross-hatched area refers to dominant compositional mode for SG2 orthopyroxenes from Figs DR1-DR3. Other details as in Fig. DR1.
normal zoning
A
reverse zoning
SG2 inheritedcores
18
4
8
12
16
Barker et al. 2016, Data Repository
Diffusion modelling in orthopyroxene A strong negative linear relationship (R2=0.95) has been observed between back-scattered
electron (BSE) image greyscale values and the Mg/(Mg + Fe) content of orthopyroxene
(Allan et al., 2013; Cooper, 2014; Chamberlain et al., 2014), and therefore the zoning
observed in BSE images is inferred to be an accurate representation of the Fe-Mg content.
This linear relationship allows compositional gradients in Fe and Mg concentrations to be
investigated at a much higher spatial resolution than is possible from spot analyses alone
(Morgan et al., 2004; Martin et al., 2008; Saunders et al., 2012; Allan et al., 2013).
Diffusional geochronometry is used here to model the evolution of compositional profiles
within zoned orthopyroxene crystals to determine the time elapsed since compositional
variations were introduced. Typically, diffusion modelling assumes that the compositional
boundaries initially had step-wise gradients, which over time at magmatic temperatures were
modified by element diffusion to form sigmoidal shaped concentration gradients until
quenching on eruption (e.g. Zellmer et al., 1999; Costa et al., 2003; Morgan et al., 2004,
2006; Costa and Dungan, 2005; Wark et al., 2007; Saunders et al., 2012; Allan et al., 2013).
The initial compositional gradient is modelled forwards in time until it matches the observed
profile and can be regarded as representing a maximum time, as the exact initial condition is
not known.
In this study we use the methods of Allan et al. (2013) to calculate timescales of Fe-
Mg diffusion in orthopyroxene. Image J (http://rsb.info.nih.gov/ij/) was used to extract
spatially resolved profiles of BSE intensity across crystal zonation boundaries to quantify
Mg/(Mg + ΣFe) profiles from rotated images. Fe-Mg profiles were obtained along the
crystallographic a- or b-axis to avoid anisotropy effects, as growth effects have been observed
along the c-axis of orthopyroxene (e.g. Allan et al. 2013; Figure DR9). Due to ƒO2
dependence of DFe-Mg in orthopyroxene (Ganguly and Tazzoli, 1994), the formula of Ganguly
and Tazzoli (1994) modified as in Allan et al. (2013) is used to calculate DFe-Mg:
where XFe is the molar proportion of the Fe end member (ferrosilite), T is temperature in
Kelvins, and ƒO2 is oxygen fugacity. Finite-difference software was used to generate a
database of simulated diffusion profiles, which obey composition-dependent diffusion under
1-D (linear) diffusion geometry, as detailed in Allan et al. (2013). This study adopts average
Barker et al. 2016, Data Repository
Figure DR9. Element maps of a representative orthopyroxene crystal from Unit Y showing the contrasts between fast diffusing (Mg) and slow diffusing (Ca+Al) elements and the relative preservation of initially sharp boundaries in a strongly zoned crystal. Light colours in WDS maps represent relatively high and dark colours are relatively low concentrations. Note the difference in zonation between the a-axis and c-axis directions, where the zoning parallel to the c-axis is smeared out and is kinematically controlled by a rapid growth regime, whereas sharp boundaries in slow-diffusing elements are largely preserved parallel to the a-axis with the slow-diffusing elements (Allan et al., 2013).
temperatures and ƒO2 values calculated using mineral-mineral and mineral-melt equilibria for
each eruptive unit from Barker et al. (2014). Oxygen fugacities were estimated using the Fe-
Ti oxide equilibrium models of Ghiorso and Evans (2008) and Sauerzapf et al. (2008). For
uncertainty calculations on single model-age determinations, conservative uncertainties of
±30 °C and ±0.3 ΔNNO log units for oxygen fugacity were used to generate inferred
maximum and minimum timescales. Statistical analyses were also conducted across each
profile, with the uncertainty on the greyscale profile essentially representing random thermal
noise in the BSE detector. The plateaux at either end of a sigmoidal diffusion profile should
be flat and profiles which had variation in the plateaux that exceeded the calculated 2
standard-error value, based on the number of averaged pixels were rejected. The relatively
large width of the diffusion profiles investigated in this study (Figs DR10-14) means that
common problems reported in other studies of convolution and pixel size were not considered
to be of significant effect (Morgan et al., 2004; Cooper, 2014). A comprehensive summary of
the measured and modelled Fe-Mg diffusion profiles across the crystal boundaries are
presented in Figs DR10-14.
Barker et al. 2016, Data Repository
horizontal distance (µm)0 10 20 30 40
Mg
#
45
46
47
48
50
horizontal distance (µm)
Mg
#
42
44
46
50X1_OPX7_reverse
0 20 40 60 80
horizontal distance (µm)
43
45
47
51X1_OPX8_Profile1_reverse
0 10 20 30 40
Mg
#
horizontal distance (µm)
Max: 20 years
Mg
#
48
52
56
58
Min: 4 years
0 20 40 60 80
X1_OPX2_normal
49
48
49
Av: 9 years
X1_OPX5_reverse
Fig. DR10. Caption over page
Max: 30 yearsMin: 6 years
Av: 14 years
Max: 42 yearsMin: 9 years
Av: 19 years
Max: 57 yearsMin: 12 years
Av: 21 years
Barker et al. 2016, Data Repository
horizontal distance (µm)
Mg
#49
51
53
57X1_OPX12_normal
0 15 30 45 60
horizontal distance (µm)
50
51
52
54X1_OPX13_Normal
0 10 20 30 40
Mg
#
X1_OPX8_Profile2_reverse
55
53
horizontal distance (µm)
43
45
47
51
0 10 20 30 40
Mg
#
49
horizontal distance (µm)
42
44
46
50
0 15 30 45 60
Mg
#
48
X1_OPX11_reverse
Fig. DR10. BSE images (left) and corresponding Fe-Mg diffusion models (right) of zoned orthopyroxene from Unit X. Yellow boxes in BSE images represent areas where diffusion modeling was undertaken. Red curves represent the modeled profile of an initially sharp compositional boundary and diamond symbols represent greyscale-calibrated Mg # (following Allan et al., 2013). Average (Av) diffusion model timescales given for parameter estimates of 860 °C, 1.5 MPa and log ƒO2 of 0.2 ∆NNO from Barker et al. (2015). Maximum (Max) and minimum (Min) timescales use uncer-tainties ±30 ºC or ± 0.3 log units ∆NNO.
Max: 47 yearsMin: 10 years
Av: 21 years
Max: 56 yearsMin: 12 years
Av: 25 years
Max: 120 yearsMin: 26 years
Av: 55 years
Max: 123 yearsMin: 26 years
Av: 56 years
Barker et al. 2016, Data Repository
Max: 355 yearsMin: 76 years
Av: 161 years
horizontal distance (µm)0 10 20 30 40
Mg
#
47
48
49
50
52
horizontal distance (µm)
Mg
#
44
46
48
50Y1_OPX3_reverse_inner
51
0 10 20 30 40
horizontal distance (µm)
45
47
49
51Y1_OPX5_reverse_inner
0 20 40 60 80
Mg
#
horizontal distance (µm)
Mg
#
44
47
49
51
0 20 40 60 80
Y1_OPX1_reverse_inner
Y1_OPX2_reverse_inner
Fig. DR11. Caption over page
Max: 21 yearsMin: 5 years
Av: 10 years
Max: 88 yearsMin: 19 years
Av: 40 years
Max: 173 yearsMin: 37 years
Av: 78 years
Barker et al. 2016, Data Repository
Max: 31 yearsMin: 7 years
Av: 14 years
Max: 142 yearsMin: 30 years
Av: 64 years
Y1_OPX6_normal
horizontal distance (µm)
Mg
#
46
47
48
50
0 10 20 30 40
horizontal distance (µm)
43
45
49
51
0 10 20 30 40
Mg
#
horizontal distance (µm)
45
47
49
51
0 20 40 60 80
Mg
#
horizontal distance (µm)
48
52
56
60
0 20 40 60 80
Mg
#
49
47
A
B
A B
Y1_OPX7_reverse_inner
Y1_OPX10_reverse_inner
Fig. DR11. Caption over page
Max: 31 yearsMin: 7 years
Av: 14 years
Max: 22 yearsMin: 5 years
Av: 10 years
Max: 130 yearsMin: 28 years
Av: 59 years
Max: 82 yearsMin: 18 years
Av: 36 years
Y1_OPX13_reverse_inners_A+B
Barker et al. 2016, Data Repository
Fig. DR11. BSE images (left) and corresponding Fe-Mg diffusion models (right) of zoned orthopyroxene from Unit Y1. Average (Av) diffusion model timescales given for parameter estimates of 860 °C, 1.5 MPa and log ƒO2 of 0.2 ∆NNO from Barker et al. (2015). Maximum (Max) and minimum (Min) timescales use uncertainties ±30 ºC or ± 0.3 log units ∆NNO. Other details as in Fig. DR10
Y1_OPX15_reverse_inner_1
Y1_OPX15_reverse_inner_2
horizontal distance (µm)
43
45
47
51
0 10 20 30 40
Mg
#
49
horizontal distance (µm)
43
45
47
51
0 10 20 30 40
Mg
#
49
Max: 33 yearsMin: 7 years
Av: 14 years
Max: 40 yearsMin: 8 years
Av: 18 years
Barker et al. 2016, Data Repository
Max: 18 yearsMin: 4 years
Av: 8 years
Y3_OPX3_reverse_inner
Y3_OPX4_reverse_inner
horizontal distance (µm)
Mg
#
44
46
48
52
0 10 20 30 40
Y3_OPX18_reverse_inner
50
horizontal distance (µm)
Mg
#
44
46
48
52
0 15 30 45 60
50
horizontal distance (µm)
Mg
#
43
45
47
51
0 30 60 90 120
49
horizontal distance (µm)
Mg
#
43
45
47
51
0 20 40 60 80
49
Y3_OPX15_reverse_inner
Fig. DR12. Caption over page
Max: 80 yearsMin: 17 years
Av: 36 years
Max: 87 yearsMin: 19 years
Av: 39 years
Max: 119 yearsMin: 25 years
Av: 54 years
Barker et al. 2016, Data Repository
Max: 91 yearsMin: 19 years
Av: 40 years
Y3_OPX9_normal_inner
Y3_OPX20_reverse_inner_A+B
horizontal distance (µm)
Mg
#
44
46
48
50
0 15 30 45 60
horizontal distance (µm)
Mg
#48
50
52
58
0 20 40 60 80
54
horizontal distance (µm)
Mg
#
43
45
47
51
0 20 40 60 80
49
A B
A
B
56
Y3_OPX5_reverse_inner
Max: 120 yearsMin: 26 years
Av: 55 years
Max: 157 yearsMin: 34 years
Av: 71 years
Max: 397 yearsMin: 85 years
Av: 180 years
Fig. DR12. BSE images (left) and corresponding Fe-Mg diffusion models (right) of zoned orthopyroxene from Unit Y3. Average (Av) diffusion model timescales given for parameter estimates of 860 °C, 1.5 MPa and log ƒO2 of 0.2 ∆NNO from Barker et al. (2015). Maximum (Max) and minimum (Min) timescales use uncertainties ±30 ºC or ± 0.3 log units ∆NNO. Other details as in Fig. DR10.
Max: 75 yearsMin: 16 years
Av: 34 years
Y7_OPX2_reverse_inner
Y7_OPX14_normal
horizontal distance (µm)
Mg
#
47
48
49
51
0 15 30 45 60
Y7_OPX1_reverse_inner_A
50
horizontal distance (µm)
Mg
#
47
48
49
51
0 15 30 45 60
50
horizontal distance (µm)
Mg
#
45
47
49
51
0 20 40 60 80
horizontal distance (µm)
Mg
#
48
52
56
60
0 10 20 30 40
Y7_OPX1_reverse_inner_B
Fig. DR13. Caption over page
Max: 76 yearsMin: 16 years
Av: 34 years
Max: 38 yearsMin: 8 years
Av: 17 years
Max: 57 yearsMin: 12 years
Av: 26 years
Barker et al. 2016, Data Repository
Max: 59 yearsMin: 13 years
Av: 27 years
Y7_OPX11_reverse_inner
Y7_OPX11_reverse_outer
horizontal distance (µm)
Mg
#
44
46
48
52
0 20 40 60 80
Y7_OPX13_reverse
50
horizontal distance (µm)
Mg
#
48
52
56
68
0 15 30 45 60
60
horizontal distance (µm)
Mg
#
48
50
52
56
0 10 20 30 40
64
horizontal distance (µm)
Mg
#
43
45
47
51
0 15 30 45 60
49
54
Y7_OPX12_normal
Fig. DR13. Caption over page
Max: 75 yearsMin: 16 years
Av: 34 years
Max: 59 yearsMin: 13 years
Av: 27 years
Max: 52 yearsMin: 11 years
Av: 24 years
Barker et al. 2016, Data Repository
Max: 101 yearsMin: 22 years
Av: 46 years
Max: 38 yearsMin: 8 years
Av: 17 years
Max: 75 yearsMin: 16 years
Av: 34 years
Max: 108 yearsMin: 23 years
Av: 49 years
Y7_OPX12_reverse_outer_2
Y7_OPX8_reverse_inner
Y7_OPX16_reverse
horizontal distance (µm)
Mg
#
50
51
52
53
0 10 20 30 40
horizontal distance (µm)
Mg
#
48
50
52
56
0 10 20 30 40
54
horizontal distance (µm)
Mg
#
45
47
49
51
0 15 30 45 60
horizontal distance (µm)
Mg
#
50
51
52
53
0 10 20 30 40
Y7_OPX12_reverse_outer_1
Fig. DR13. Caption over page
Barker et al. 2016, Data Repository
Fig. DR13. BSE images (left) and corresponding Fe-Mg diffusion models (right) of zoned orthopyroxene from Unit Y7. Average (Av) diffusion model timescales given for parameter estimates of 860 °C, 1.5 MPa and log ƒO2 of 0.2 ∆NNO from Barker et al. (2015). Maximum (Max) and minimum (Min) timescales use uncertainties ±30 ºC or ± 0.3 log units ∆NNO. Other details as in Fig. DR10.
Max: 125 yearsMin: 27 years
Av: 57 years
Y7_OPX5_reverse
horizontal distance (µm)
Mg
#
46
48
50
52
0 15 30 45 60
Barker et al. 2016, Data Repository
horizontal distance (µm)0 20 40 60 80
Mg
#
50
51
52
53
55
horizontal distance (µm)
Mg
#
46
48
50
52Z_OPX14_reverse
54
0 15 30 45 60
Z_OPX15_reverse
horizontal distance (µm)
Mg
#
48
49
50
53
0 10 20 30 40
Z_OPX3_reverse_outer
52
51
AB
AB
horizontal distance (µm)
Mg
#
46
48
50
52
0 15 30 45 60
Z_OPX13_normal_outer
Fig. DR14. Caption over page
Max: 8 yearsMin: 2 years
Av: 4 years
Max: 122 yearsMin: 26 years
Av: 58 yearsMax: 129 yearsMin: 27 years
Av: 59 years
Max: 46 yearsMin: 10 years
Av: 21 years
Max: 80 yearsMin: 17 years
Av: 36 years
Barker et al. 2016, Data Repository
Fig. DR14. BSE images (left) and corresponding Fe-Mg diffusion models (right) of zoned orthopyroxene from eruption Z. Average (Av) diffusion model timescales given for parameter estimates of 860 °C, 1.5 MPa and log ƒO2 of 0.0 ∆NNO from Barker et al. (2015). Maximum (Max) and minimum (Min) timescales use uncertainties ±30 ºC or ± 0.3 log units ∆NNO. Other details as in Fig. DR10.
Max: 9 yearsMin: 2 years
Av: 4 years
horizontal distance (µm)
Mg
#
48
49
50
53
0 10 20 30 40
Z_OPX19_reverse_outer
52
51
Barker et al. 2016, Data Repository
REFERENCES CITED IN SUPPLEMENTARY MATERIAL
Allan, A.S.R., Morgan, D.J., Wilson, C.J.N., and Millet, M-A., 2013, From mush to eruption in
centuries: assembly of the super-sized Oruanui magma body: Contributions to Mineralogy and
Petrology, v. 166, p. 143-164, doi:10.1007/s00410-013-0869-2.