ABCDEFGHIJKLMNOPQRSTUVWXYZ
1
Cumulate and residual wallrock major element calculation tool for Magma Chamber Simulator
2
Developed by Jussi Heinonen (jussi.s.heinonen@helsinki.fi)
3
4
Tabs
Instructions
5
Color coding:
1. Clear contents (do not use delete rows or columns command!) of "Input" and "RS" tabs
6
Input
2. Copy and paste values from MCS Input to the "Input" tab (starting from cell A1)
7
Calculations (DO NOT MODIFY!)
3. Copy and paste values from MCS RunSummary to the "RS" tab (starting from cell A1)
8
Output
4. Collect the most fitting form of output from the three output tabs
9
10
Input
Copy and paste values of "Input" tab from MCS output
11
RS
Copy and paste values of "RunSummary" tab from MCS output
Error
12
R
Recharge magma compositions
The negative seeds are calculated out from incremental cumulate (to avoid small negative oxide concentrations),
13
CALC
Main calculation tab
but are included in the bulk cumulate (to ensure consistency with MCS throughout the run)
14
Rows
Rows compatible with MCS-Traces Magma and Wallrock output
Negative seed increments may thus cause minor discrepancies in the incremental vs. bulk cumulate output within one step
15
CopyToRS
This can be copied and pasted to "RunSummary" tab of MCS output
This discrepancy should usually not exceed total sum error of 0.5%
16
CopyToTracesM
This can be copied and pasted to "Traces_X_Magma" tab of MCS-Traces output
MAX SEED ERROR =
0.00000%
17
CopyToTracesWR
This can be copied and pasted to "Traces_X_Wallrock" tab of MCS-Traces output
18
19
Equations
20
The central equation behind the incremental cumulate calculations is:
21
Xic = (Xp-(Xc(1-ms/(ml+ms+mf))))/(ms/(ml+ms+mf)), where
22
Xic = oxide composition of incremental cumulate in the current step
23
Xp = oxide composition of liquid (melt) + fluid in the previous step
24
Xc = oxide composition of liquid (melt) + fluid in the current step
25
ms = mass of incremental solids in the current step
26
ml = mass of liquid (melt) in the current step
27
mf = mass of cumulative fluid in the current step
28
29
The central equation behind the residual WR calculations is:
30
Xr = Xp-Xl(ml/(ml+mr)), where
31
Xr = composition of WR residual in the current step
32
Xp = composition of bulk WR in the previous step
33
Xl = composition of WR liquid in the current step
34
ml = mass of WR liquid added to M in the current step
35
mr = mass of residual WR in the current step
36
This is then further refined to also calculate out melt (below percolation threshold) and fluid from the residual WR and output solid WR residual composition.
37
38
All calculations are normalized to 100 wt.% for output.
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100