Electron microprobe data for tephra attributed to Glacier Peak,Washington

Reference samples of three prominent pumice units of Glacier Peak tephra collected east of the volcano within a distance of 100 km are similar petrographically to units described by earlier workers. Glass shards isolated from these samples were analyzed by electron microprobe to determine the content of Ca, Fe, and K. Resulting data, plus those published for two other references samples, provide a basis for attributing certain outlying tephra layers from 14 locations in eastern Washington, Idaho, Wyoming, and Montana to eruptions of Glacier Peak. Ten of the samples have properties of both Glacier Peak tephra and Mount St. Helens set J tephra, but proportions of Ca:Fe:K in glass shards indicate that 9 of the 10 outlying samples came from Glacier Peak, whereas one is assigned to Mount St. Helens set J. The remaining six outlying samples, all from southeastern Washington, contain cummingtonite phenocrysts and are chemically similar to some parts of Mount St. Helens tephra sets that are older than 12,000 BP.     

The distribution of divalent trace elements between sulfides,oxides, silicates and hydrothermal solutions: I. Thermodynamic basis

Experimental data for the standard Gibbs free energies of formation from the elements of a wide variety of metal sulfides and oxides, spinels, olivines and pyroxenes at 25°C and 1 bar define linear correlations, within about ±900 cal·mole^?1, with the corresponding conventional standard partial molal Gibbs free energies of formation of the aqueous M²⁺ cations of the form where a_MaZ and b_MaZ are empirically determined constants characteristic of the structure M_nZ. The only exceptions to correlations of this type are compounds of the heavy alkaline earths Ca, Sr and Ba, which appear to follow correlations with cation radius instead. The linear free energy correlations enable prediction of standard Gibbs free energies of formation of compositional end-members of a particular structure M_nZ provided that a_MaZ and b_MaZ are known accurately. When only the free energy of the Mg end-member is known, the standard Gibbs free energy of formation at 25°C and 1 bar of the Fe endmember, and hence a_MaZ and b_MaZ Can be predicted from the temperature independence of and estimated entropies and heat capacities for the Fe end-member. Using this approach, the free energies of ferrosilite, hedenbergite and annite at 25°C and 1 bar were predicted to within ±1000 cal·mole^?1 of the helgesonet al. (1978) values. Free energies of formation of talc (M₃Si₄O₁₀(OH)₂), clinchlore (M₅Al₂Si₃O₁₀(OH)₈), and tremolite (Ca₂M₅(Si₄O₁₁)₂(OH)₂)-type compounds where M is Mg, Mn, Zn, Fe, Co, or Ni were then predicted at 25°C and 1 bar.Calculation of the equilibrium distribution of Mg, Zn and Sr between galena and hydrothermal solution, and Zn, Mg, Fe and Mn between chlorite and hydrothermal solution demonstrates: (1) that the Sr contents of low temperature galenas (e.g. Mississippi Valley-type) should be negligible (reported analyses of Sr content and Sr isotopic composition of such galenas are probably attributable to fluid inclusions or carbonate inclusions); and (2), that the Zn contents of hydrothermal chlorites in a model of the midoceanic ridge hydrothermal systems are sensitive to temperature, to complexing in the aqueous phase, and to the overall Fe/Mg ratio of the chlorite.     

Origin of trondhjemitic biotite—quartz—oligoclase gneisses from the Venezuelan Guyana Shield

Biotite—quartz—oligoclase gneisses constitute the dominant lithology in a 3-b.y.-old metamorphic assemblage in the Venezuelan Guyana Shield. The assemblage includes basaltic amphibolites as well as granitic gneisses, iron formation, and other metasedimentary lithologies.Consideration of major- and trace-element compositions indicates that both the biotite gneisses and basaltic amphibolites are meta-igneous. The amphibolites have oceanic tholeiite compositions.Two groups of biotite gneisses can be differentiated by chemical criteria. Both groups have major and trace-element compositions which allow their derivation as partial melts of tholeiite compositions at mantle depth, as has been suggested for similar rocks in other areas. However, the compositional correspondence of the gneisses with low variance liquidus loci in low-P_T synthetic systems, and their relatively oxidized character, appear more compatible with an origin by partial melting of graywackes at crustal levels. For both groups of gneisses, primary melts with relatively low $\text{Na}{}_2\text{O}\text{K}{}_2\text{O}$ ratios can be postulated which are credibly derived by melting of graywackes; reasonable fractionation processes can be hypothesized to explain the compositional variations within each group.Although stratigraphic relations in the Venezuelan Guyana Shield are uncertain, geologic relations in and near the map area allow, as valid working hypotheses, stratigraphic sequences which parallel those which have been recognized in better-known Archean terranes.     

Some chemographic relationships in n-component systems

A petrologic problem of fundamental importance is to determine whether 2 or more mineral assemblages can be related to one another by continuous or discontinuous facies changes, or whether their bulk compositions occupy non-overlapping regions of composition space. A general method is developed by which 2 regions of n-dimensional space whose vertices are defined by the phases present are tested for compositional overlap. This is accomplished by generating mass balance equations of the type:

$\text{∑}\text{i = 1}\text{m}\text{a}{}_{\mathrm{i}}\text{A}{}_{\mathrm{i}}\text{=}\text{∑}\text{j = 1}\text{k}\text{b}{}_{\mathrm{j}}\text{B}{}_{\mathrm{j}}$

where A_i is the ith phase in one region and B_j is the th phase in the other. If any such equation satisfies the requirement that the sign of each a_i is the same, and that the sign of each b_j is the opposite for all i, j such that: k + m = n + 1 then the 2 regions overlap in phase space.By eliminating all overlapping assemblages in a given set, the bivariant fields bounded by univariant equilibria in n-dimensional systems are completely specified. All bulk compositions are considered within the space defined by the phases that participate in the bounding reactions. An extension of the method generates in sequence all bivariant fields and associated reactions about any invariant point. A further extension is applied to multi-system analysis.     

The effect of ionic interactions on the oxidation of metals in natural waters

The effect of ionic interactions of the major components of natural waters on the oxidation of Cu(I) and Fe(II) has been examined. The various ion pairs of these metals have been shown to have different rates of oxidation. For Fe(II), the chloride and sulfate ion pairs are not easily oxidized. The measured decrease in the rate constant at a fixed pH in chloride and sulfate solutions agrees very well with the values predicted. The effect of pH (6 to 8) on the oxidation of Fe(II) in water and seawater have been shown to follow the rate equation

$\text{-d in [Fe(II)]/dt = k}{}_1\text{β}{}_1\text{α}{}_{\mathrm{Fe}}\text{/[H}{}^{\mathrm{+}}\text{] + k}{}_2\text{β}{}_2\text{α}{}_{\mathrm{Fe}}\text{/[H}{}^{\mathrm{+}}\text{]}{}^2$

where k₁ and k₂ are the pseudo first order rate constants, β₁ and β₂ are the hydrolysis constants for Fe(OH)⁺ and Fe(OH)⁰. The value of α_FE is the fraction of free Fe²⁺. The value of k₁ (2.0 ±0.5 min^?1) in water and seawater are similar within experimental error. The value of k₂ (1.2 × 10⁵ min^?1) in seawater is 28% of its value in water in reasonable agreement with predictions using an ion pairing model.For the oxidation of Cu(I) a rate equation of the form

$\text{?d ln [Cu(I)]/dt = k}{}_0\text{α}{}_{\mathrm{Cu}}\text{+ k}{}_1\text{β}{}_1\text{α}{}_{\mathrm{Cu}}\text{[Cl]}$

was found where k₀ (14.1 sec^?1) and k₁ (3.9 sec^?1) are the pseudo first order rate constants for the oxidation of Cu⁺ and CuCl⁰, β₁ is the formation constant for CuCl⁰ and α_Cu is the fraction of free Cu⁺. Thus, unlike the results for Fe(II), Cu(I) chloride complexes have measurable rates of oxidation.     

Liquidus relationships in the system CaCO3Ca(OH)2CaS and the solubility of sulfur in carbonatite magmas

CaCO₃Ca(OH)₂CaS serves as a model system for sulfide solubility in carbonatite magmas. Experiments at 1 kbar delineate fields for primary crystallization of CaCO₃, Ca(OH)₂ and CaS. The three fields meet at a ternary eutectic at 652°C with liquid composition (wt%): CaCO₃ = 46.1%, Ca(OH)₂ = 51.9%, CaS = 2.0%. Two crystallization sequences are possible for liquids that precipitate calcite, depending upon whether the liquid is on the low-CaS side, or the high-CaS side of the line connecting CaCO₃ to the eutectic liquid. Low-CaS liquids precipitate no sulfide until the eutectic temperature is reached leading to sulfide enrichment. The higher-CaS liquids precipitate some sulfide above the eutectic temperature, but the sulfide content of the melt is not greatly depleted as the eutectic temperature is approached. Theoretical considerations indicate that sulfide solubility in carbonate melts will be directly proportional to and inversely proportional to ; it also is likely to be directly proportional to melt basicity, defined here by . A strong similarity exists in the processes which control sulfide solubility in carbonate and in silicate melts. By analogy with silicates, ferrous iron, which was absent in our experiments, may also exert an important influence on sulfide solubility in natural carbonatite magmas.     

Changes in the content and composition of pedogenic iron oxyhydroxides in a chronosequence of soils in southern California

Studies of the pedogenic iron oxyhydroxides in suites of latest Holocene to middle Pleistocene soils formed on fluvial deposits of the transverse ranges, southern California, indicate that the content and composition of iron oxyhydroxide change in a systematic manner. Analysis of total secondary free iron oxides (dithionite extractable, Fe₂O₃d) and ferrihydrite (oxalate extractable, Fe₂O₃o) shows that (1) a single-logarithmic model (Y = a + b log X) or double logarithmic model (log Y = a + b log X), where Y is the total mass of pedogenic Fe oxides (g/cm²-soil column) and X is soil age, describes the rate of increase in Fe₂O₃d with time; (2) the Fe₂O₃d content correlates linearly with soil reddening and clay content; (3) the $\text{Fe}{}_2\text{O}{}_3\text{o}\text{Fe}{}_2\text{O}{}_3\text{d}$ ratio, which indicates the degree of Fe oxide crystallinity, is moderately high to very high (0.22–0.58) in middle Holocene to latest Pleistocene soils and progressively decreases to less than 0.10 in older soils; (4) the value of the $\text{Fe}{}_2\text{O}{}_3\text{o}\text{Fe}{}_2\text{O}{}_3\text{d}$ ratio also appears to be infuenced by climate; and (5) temporal changes in Fe oxide content and mineralogy are accompanied by related, systematic changes in clay mineralogy and organic matter content. These relationships are attributed to a soil environment that must initially favor ferrihydrite precipitation and/or organic matter-Fe complexation. Subsequent transformation to hematite causes increasingly intense reddening and a concomitant decrease in the $\text{Fe}{}_2\text{O}{}_3\text{o}\text{Fe}{}_2\text{O}{}_3\text{d}$ ratio. The results demonstrate that iron oxide analysis is useful for numerical age studies of noncalcic soils and shows potential as an indicator of paleoclimates.     

Oxygen and sulfur isotope equilibria in the BaSO4HSO4−H2O system from 110 to 350°C and applications

Oxygen isotope exchange between BaSO₄ and H₂O from 110 to 350°C was studied using 1 m H₂SO₄-1 m NaCl and 1 m NaCl solutions to recrystallize the barite. The slow exchange rate (only 7% exchange after 1 yr at 110°C and 91% exchange after 22 days at 350°C in 1 m NaCl solution) prompted the use of the partial equilibrium technique. However, runs at 300 and 350°C were checked by complete exchange experiments. The temperature calibration curve for the isotope exchange is calculated giving most weight to the high temperature runs where the partial equilibrium technique can be tested. Oxygen isotope fractionation factors (α) in 1 m NaCl solution (110–350°C), assuming a value of 1.0407 for α_CO2H2O at 25°C, are:

.These data, when corrected for ion hydration effects in solution (Truesdell, 1974), give the fractionation factors in pure water:

.In the 1 m H₂SO₄-1 m NaCl runs, sulfur isotope fractionation between HSO^?₄ and BaSO₄ is less than the detection limit of 0.4%. A barite-sulfide geothermometer is obtained by combining HSO^?₄H₂S and sulfide-H₂S calibration data.Barite in the Derbyshire ore field, U.K., appears to have precipitated in isotopic equilibrium with water and sulfur in the ore fluid at temperatures less than 150°C. At the Tui Mine, New Zealand, the barite-water geothermometer indicates temperatures of late stage mineralization in the range 100–200°C. A temperature of 350 ± 20°C is obtained from the barite-pyrite geothermometer at the Yauricocha copper deposit, Peru, and oxygen isotope analyses of the barite are consistent with a magmatic origin for the ore fluids.     

Density calculations for silicate liquids. I. Revised method for aluminosilicate compositions

Equations are developed for calculating the density of aluminosilicate liquids as a function of composition and temperature. The mean molar volume at reference temperature T_r, is given by $\text{V}{}_{\mathrm{r}}\text{= ∑X}{}_{\mathrm{i}}\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{i}}\text{+ X}{}_{\mathrm{A}}\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$, where the summation is taken over all oxide components except A1₂O₃, X stands for mole fraction, terms are constants derived independently from an analysis of volume-composition relations in alumina-free silicate liquids, and $\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is the composition-dependent apparent partial molar volume of Al₂O₃. The thermal expansion coefficient of aluminosilicate liquids is given by $\text{α = ∑X}{}_{\mathrm{i}}\text{}\text{?}\text{ga}{}_{\mathrm{i}}{}^{\mathrm{o}}\text{+ X}{}_{\mathrm{A}}\text{}\text{?}\text{ga}{}_{\mathrm{A}}{}^{\mathrm{o}}$, where $\text{}\text{?}\text{ga}{}_{\mathrm{i}}{}^{\mathrm{o}}$ terms are constants independent of temperature and composition, and $\text{}\text{?}\text{ga}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is a composition-dependent term representing the effect of Al₂O₃ on the thermal expansion. Parameters necessary to calculate the volume of silicate liquids at any temperature T according to V(T) = V_rexp[α(T-T_r)], where T_r = 1400°C have been evaluated by least-square analysis of selected density measurements in aluminosilicate melts. Mean molar volumes of aluminosilicate liquids calculated according to the model equation conform to experimentally measured volumes with a root mean square difference of 0.28 $\text{cc}\text{mole}$ and an average absolute difference of 0.90% for 248 experimental observations. The compositional dependence of $\text{V}\text{?}{}^{\mathrm{o}}{}_{\mathrm{A}}$ is discussed in terms of several possible interpretations of the structural role of Al³⁺ in aluminosilicate melts.     

Electron microprobe analysis of glass shards from tephra assigned to set W,Mount St. Helens,Washington

We have extended the fallout areas for each of two members of tephra-set W, erupted from Mount St. Helens about 1500 ad, by several hundred kilometers beyond the limits mapped in 1975. We traced one member (We) east into Idaho, and the other (Wn) northeast into British Columbia. After using stratigraphic and petrographic observations to assign more than 100 tephra samples to set W, we found 26 of these, selected for chemical analysis, to be closely similar in content of Ca, Fe, and K in glass shards. But improved homogeneity was evident when the 26 sampling localities for tephra W were segregated geographically, east vs. northeast of the volcano. When Ca:Fe:K proportions were plotted on a ternary diagram, there was no overlap of the plotting areas for these two groups of tephra W samples. Without such data, tephra layers We and Wn are currently separable only from stratigraphic and geographic information. Partial glass analysis is also an aid, along with stratigraphic position and petrographic characteristics, in distinguishing tephra W from associated tephra layers. These include tephra layers T and Yn from Mount St. Helens, as well as older tephra layers from Mount Mazama and Glacier Peak.     

Evaluation of deep temperatures of hydrothermal systems by a new gas geothermometer

The chemical composition of gas mixtures emerging in thermal areas can be used to evaluate the deep thermal temperatures. Chemical analyses of the gas compositions for 34 thermal systems were considered and an empirical relationship developed between the relative concentrations of H₂S, H₂, CH₄ and CO₂ and the reservoir temperature. The evaluated temperatures can be expressed by: $\text{t°C =}\text{24775}\text{α + β + 36.05}\text{?273}$ where $\text{α = 2 log}\text{CH}{}_4\text{CO}{}_2\text{?log}\text{H}{}_2\text{CO}{}_2\text{?3 log}\text{H}{}_2\text{S}\text{CO}{}_2$ (concentrations in % by volume) and β = 7 logP_co2     

The role of CO2 in the chemical modification of deep continental crust

Compositional differences between granulite facies rocks and equivalent amphibolite facies rocks and the observation of CO₂-rich fluid inclusions in granulites, have led to the suggestion that CO₂ must play a role in modifying the composition of deep continental crust. How CO₂ effects this change has remained unclear. Using the thermodynamic properties of aqueous ions in a fluid of evolving $\text{CO}{}_2\text{H}{}_2\text{O}$ ratio, it is possible to model the incongruent dissolution of feldspars under conditions appropriate for granulite facies metamorphism. The results demonstrate that dissolution will be strongly enhanced at high $\text{CO}{}_2\text{H}{}_2\text{O}$ ratios, with ion solubilities being Na⁺ >K⁺ ? Ca⁺⁺. This enhancement is compatible with the reported compositional contrasts between granulite and amphibolite facies rock, but requires large fluid volumes.To test the dissolution model, a detailed field and petrologic study was conducted in a well exposed granulite facies terrane in West Greenland. Strong correlation between fluid composition and bulk rock chemistry can be documented; CO₂-rich regions contain rocks which consistently have low ratios, while H₂O-rich regions consistently have high ratios. Magnetite rims on sulfide grains are ubiquitous in high regions and are absent in high regions, and they provide evidence that CO₂ was introduced into the region. These correlations and observations are predictable from the properties of the dissolution process. These considerations, along with observations regarding graphite petrogenesis, provide strong arguments that the total fluid volume interacting with the rock during metamorphism was very large, in some cases equaling or exceeding total rock volume. Such large fluid volumes can lead to significant compositional modification of the crust, and will mask the original protolith chemistry. Such processes should lead to Ca- and Al-enriched, Na-, K-, S- and Si-depleted residues in the deep crust.     

Titan-aegirine from early Tertiary ash layers in northern Denmark

Microprobe analyses on a xenocrystic suite of salites, aegirine-augites, aegirines, titan-aegirines and acmites from a lower Tertiary ash layer in northern Denmark are presented. The sodic pyroxenes show an unusual titan-enrichment and up to 42 mol.% of the component $\text{NaTi}{}_{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}{}^{\mathrm{4+}}\text{M}{}_{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}{}^{\mathrm{2+}}\text{Si}{}_3\text{O}{}_6\text{(}\text{M = Fe}{}^{\mathrm{2+}}\text{,}\text{Mn or Mg}\text{)}$, is estimated. Optical absorption measurements show no evidence for Ti³⁺. The titan-aegirines were formed during late to post-magmatic crystallization in a system with a high Ti⁴⁺/Fe²⁺ ratio and were followed by acmite showing enrichment in jadeite. Comparison with experimentally investigated titan-aegirine indicates crystallization far below the Mn₂O₃Mn₃O₄f₀₂ buffer.     

Viscosity and configurational entropy of silicate melts

With the configurational entropy theory of relaxation processes of Adam and Gibbs (1965), one predicts that the viscosity depends on temperature according to log $\text{η = A}{}_{\mathrm{e}}\text{+}\text{B}{}_{\mathrm{e}}\text{TS}{}_{\mathrm{<mtext>conf</mtext>}}$, where S_conf is the configurational entropy of the liquid. Thermochemical calculations of S_conf performed for some mineral compositions show the importance of non-configurational contributions to the entropy differences between amorphous and crystalline phases. Except for the case of SiO₂, the available thermodynamic data indicate that the above equation for viscosity accounts quantitatively for the experimentally determined temperature dependence of the viscosity of silicate melts. The Adam and Gibbs theory also provides a simple rationale for the non linear variation of the logarithmic viscosity with composition in mixed alkali silicate liquids at low temperatures, the minimum of viscosity resulting from the contribution of the entropy of mixing to S_conf.     

Light hydrocarbons in Red Sea brines and sediments

Light hydrocarbon (C₁-C₃) concentrations in the water from four Red Sea brine basins (Atlantis II, Suakin, Nereus and Valdivia Deeps) and in sediment pore waters from two of these areas (Atlantis II and Suakin Deeps) are reported. The hydrocarbon gases in the Suakin Deep brine (T = ~ 25°C, $\text{Cl}{}^{\mathrm{?}}\text{= ~ 85‰}$, $\text{CH}{}_4\text{=}\text{~ 71}\text{1}$) are apparently of biogenic origin as evidenced by $\text{C}{}_1\text{(C}{}_2\text{+ C}{}_3$) ratios of ~ 1000. Methane concentrations (6–8 μl/l) in Suakin Deep sediments are nearly equal to those in the brine, suggesting sedimentary interstitial waters may be the source of the brine and associated methane.The Atlantis II Deep has two brine layers with significantly different light hydrocarbon concentrations indicating separate sources. The upper brine (T = ~ 50°C, $\text{Cl}{}^{\mathrm{?}}\text{= ~ 73‰}$, CH₄ = ~ 155 μl/l) gas seems to be of biogenic origin [$\text{C}{}_1\text{(C}{}_2\text{+ C}{}_3\text{) = ~1100}$], whereas the lower brine (T = ~ 61°C, $\text{Cl}{}^{\mathrm{?}}\text{= ~ 155‰}$, CH₄ = ~ 120μl/l) gas is apparently of thermogenic origin [$\text{C}{}_1\text{(C}{}_2\text{+ C}{}_3\text{) = ~ 50}$]. The thermogenic gas resulting from thermal cracking of organic matter in the sedimentary column apparently migrates into the basin with the brine, whereas the biogenic gas is produced in situ or at the seawater-brine interface. Methane concentrations in Atlantis II interstitial waters underlying the lower brine are about one half brine concentrations; this difference possibly reflects the known temporal variations of hydrothermal activity in the basin.     

Polyunsaturated fatty acids in a lacustrine sediment as a possible indicator of paleoclimate

Polyunsaturated fatty acids (C_18:2 and C_18:3ω3 were analyzed in the upper 20m layer of a 200 m long sediment core taken from Lake Biwa. Concentration maxima occur in layers at depths of 0.2, 1–5, 11–12, and 16m. The vertical changes in the $\text{(C}{}_{\mathrm{18:2}}\text{C(C}{}_{\mathrm{18:0}}$ ratio appear to correlate with paleoclimatic condition suggested from palynological evidence. On the basis of $\text{C}{}_{\mathrm{18:2}}\text{C}{}_{\mathrm{18:0}}$ ratios, it was suggested that it has been colder at 200, 1000–4000, 15,000 and 20,000 yr BP than at other times.     

Kiglapait Geochemistry III: Potassium and Rubidium

The regular geometry and completeness of the Kiglapait intrusion permit its bulk composition to be obtained by summation, and the composition of successive liquids to be obtained by subtraction. The summations for K and Rb give 1806 and 1.08 ppm, yielding Rfrsol|K/Rb= 1670 for the intrusion, taken as equal to the parent magma. R increases slightly from this initial value to 2000 at the end of crystallization where MgO approaches zero in the rocks. K and Rb are therefore closely coherent and their distribution coefficients can differ only by a small amount in the Kiglapait system.Apparent feldspar/liquid distribution coefficients (D^F/L) can be estimated from detailed plots of feldspar and liquid compositions against F_L. The Kiglapait data imply that these coefficients are linear 1:1 functions of plagioclase composition within experimental error, having values given by D_KF/L = 1.42? X_AnD_RbF/L = 1.13? X_An with minimum values of 0.75 and 0.49, respectively. The ratio $\text{R}{}^{\mathrm{F}}\text{R}{}^{\mathrm{L}}$ lies in the range of 1.53± 0.03 for the plagioclase composition range X_An= 0.34 to 0.67 showing that high-R rocks such as anorthosite crystallized from high-R liquids.The apparent feldspar distribution coefficients are much closer to 1.0 than common literature values. They can be reduced by assuming that the cumulate pile was continuously recharged by the circulating magma until an advanced stage of differentiation was reached, and assuming that alkalies were exchanged to the feldspars from the magma. When such an ‘aquifer recharge’ model is calibrated using olivine-liquid equilibria as a time marker for the liquid, the inferred minimum equilibrium values of the distribution coefficients are $\text{D}{}_{\mathrm{K}}\text{F}\text{L}$= 0.42, $\text{D}{}_{\mathrm{Rb}}\text{F}\text{L}$ = 0.25 at the base of the intrusion. Their variation is given by $\text{D}{}_{\mathrm{K}}\text{F}\text{L}\text{= 1.66?1.88X}{}_{\mathrm{An}}$, $\text{D}{}_{\mathrm{Rb}}\text{F}\text{L}\text{= 1.17?1.41X}{}_{\mathrm{An}}$, The equilibrium values are considered to be appropriate for deducing liquid compositions in plutonic bodies where alkali exchange can be shown or inferred to have been inhibited, such as in small intrusions. The apparent values are considered to be appropriate, even though they may be artificial, for large intrusions similar to the Kiglapait.The bulk K and Rb concentrations in the Kiglapait intrusion are consistent with a plagioclase-rich abyssal tholeiite magma. Clinopyroxene and olivine fractionation in the mantle may contribute to the production of such high-Rmagmas.     

Thermodynamics of brine-salt equilibria—I. The systems NaCl-KCl-MgCl2-CaCl2-H2O and NaCl-MgSO4-H2O at 25°C

A thermodynamic model for concentrated brines has been developed which is capable of predicting the solubilities of many of the common evaporite minerals in chloro-sulfate brines at 25°C and 1 atm. The model assumes that the behaviour of the mean stoichiometric ionic activity coefficient in mixtures of aqueous electrolytes can be described by the Scatchard deviation function and Harned's Rule. In solutions consisting of one salt and H₂O, the activity coefficient is described by the expression where $\text{a}\text{?}$ and B? salt specific parameters obtained from data regression. In a mixture of n electrolytes and H₂O, B? for the ith component is given by $\text{B}{}^{\mathrm{<mtext>i</mtext><mtext>?</mtext>}}{}_{\mathrm{i}}\text{=B}\text{i}\text{?}{}_{\mathrm{i}}\text{+σ α}{}_{\mathrm{ij}}\text{y}{}_{\mathrm{j}}$ where α_ij is a (constant) mixing parameter characterizing the interaction of the i and j components and y_j is the ionic strength fraction of the jth component. The activity of H₂O is obtained from a Gibbs-Duhem integration and does not require any additional parameters or assumptions. In this study, parameters have been obtained for the systems NaCl-KCl-MgCl₂-CaCl₂-H₂O and NaCl-MgSO₄-H₂O at 25°C and 1 atm. Computed solubility curves and solution compositions predicted for invariant points in these systems agree well with the experimental data. The model is flexible and easily extended to other systems and to higher temperatures.