PACROFI VI - Electronic Program
CalcicBrine: a Microsoft Excel 5.0 Add-in for calculating salinities from microthermometric data in the system NaCl-CaCl2-H2O
The British Geological Survey
Sir Kingsley Dunham Centre
NG12 5GG, UK
Experimental data (Yanatieva, 1946; Vanko et al., 1988; Oakes et al., 1989) in the system NaCl-CaCl2-H2O have been regressed, using a stepwise procedure, to give the following equations:
W = co + aoT + a1R + a2T2 + a3R2 + a4RT2 + a5T3 + a6T4 + a7T5 + a8T5R ice field
W = co + aoT + a1R + a2R2 + a3R3 + a4TR hydrohalite field
W = co + aoT + a1R + a2T2 + a3R2 + a4R2T + a5R3 + a6TR halite field
R = co + aoT + a1T2 + a2T3 + a3T4 ice and halite cotectics
W = bulk salinity (weight percent)
R = NaCl/(NaCl+CaCl2) weight ratio
T = Tm ice/100, T m hyd/100 or T m halite/100
The coefficients for above equations are as follows:
|ice field||hydrohalite field||halite field||ice cotectic||halite cotectic|
An algorithm, written in Visual Basic for Applications, has been developed to calculate weight percent equivalent NaCl + CaCl2 and the NaCl/(NaCl+CaCl2) weight ratio from microthermometric data (Tm ice, Tm hydrohalite, Tm halite). The algorithm uses similar principles to those described in Bodnar et al. (1989) to solve the above equations for T, using paired input data of Tice-Thyd for compositions in the ice and hydrohalite fields or Tice-Thalite for compositions in the halite field.
Experimental bulk-salinity data are generally reproduced to within 3 percent for finalmelting in the icefield, 10 percent for bulk compositions in the hydrohalite-field and 6 percent for compositions in the halite-field. The regression equations tend to model the data of Oakes et al. (1990) and Vanko et al. (1989) better than the data Yanatieva (1946).
Limitations to the applicability algorithm are as follows:
The function will not calculate salinity data where input data do not meet criteria 1 to 4 or where microthermometric data indicate compositions in the antarcticite field. However, the function will calculate salinities where ice melts metastably in the presence of halite. This is due to the construction of the algorithm; it uses ice-melting temperatures to calculate Na to Ca ratios for final melting in all three fields of the system. Thus, to preserve flexibility, it does not discriminate between stable and metastable ice melting. The algorithm is currently being modified to include a metastable ice-hydrohalite cotectic (see discussion in Vanko et al. 1988 p. 2455).
- Final-ice melting (Tm ice)must be between -52 and 0.0oC
- Hydrohalite melting (Tm hydrohalite) must be between -52 and 0.1oC
- Final ice-melting and hydrohalite melting must not be above -21.2oC
- Halite (Tm halite) dissolution should be between 0.1 and 500oC
- For final-melting/dissolution in the halite field, ice-melting must be in the presence of hydrohalite
Copies of the Excel Add-in are available, by e-mail, from the author (firstname.lastname@example.org).
||System 7.0, 16MB RAM||Microsoft Excel 5.0|
|IBM PC or compatible||Windows 3.1 or `95, 16MB RAM||Microsoft Excel 5.0|
- Bodnar, R.J., Sterner, S.M., and Hall, D.L. (1989) Salty: a FORTRAN program to calculate compositions of fluid inclusions in the system NaCl-KCl-H2O. Computers and Geosciences, 15, p.19-41.
- Oakes, C.S., Bodnar, R.J. and Simonson, J.M. 1990. The system NaCl-CaCl2-H2O: I. The ice liquidus at 1 atm. total pressure. Geochimica et Cosmochimica Acta, 54, p.603-610.
- Vanko, D.A., Bodnar, R.J. and Sterner, S.M. 1988. Synthetic fluid inclusions: VIII. Vapour-saturated halite solubility in part of the system NaCl-CaCl2-H2O, with applications to fluid inclusions from oceanic hydrothermal systems. Geochimica et Cosmochimica Acta, 52, p.2451-2456.
- Yanatieva, O.K. 1946. Polythermal solubilities in the systems CaCl2-MgCl2-H2O and CaCl2-NaCl-H2O. Zhur. Priklad. Khim., 19, p.709-722 (in Russian).