Introduction
============

Yucca Mountain, NV, is the designated site for a permanent geologic repository for high-level nuclear waste in the USA. The current waste package design consists of a double-walled container with an inner barrier of stainless steel, an outer barrier of highly corrosion resistant nickel--chromium--molybdenum alloy, and a titanium alloy drip-shield that covers the containers. Corrosion resistance and long-term integrity of the metal containers and shields are important for the safe disposal of the waste. Characterization of the compositional evolution of waters that affect waste package corrosion is necessary. If the site is licensed, the waste packages will be placed in tunnels several hundred meters below the ground surface and above the groundwater table in partially saturated volcanic tuff. Once the waste packages are in place, the repository will heat up due to the thermal energy of the nuclear waste. Although the waste packages will be above the groundwater table, pore water present in rock formations within (Topopah Spring Tuff) and above (Paintbrush Tuff) the repository may come into contact with the metal containers and shields. Additionally, brines may form from the deliquescence of salts found in dusts deposited on the containers.\[[@B1]\] In this study we focus on seepage brines formed by the evaporation of pore water at elevated temperature.

One method of evaluating the evolution of a brine is the chemical divide theory, which has been used to describe saline lake geochemistry.\[[@B2]-[@B5]\] The chemical divide theory generally describes the chemical evolution of dilute waters upon evaporation in terms of their equivalent calcium, sulfate, and bicarbonate ratios and is shown in Fig. [1(a)](#F1){ref-type="fig"}. The chemical evolution of evaporating water is controlled by the high solubility of salt minerals relative to the moderate solubility of calcium sulfate and low solubility of calcium carbonate minerals. A bicarbonate alkaline brine (Na--K--**CO**~**3**~--Cl--SO~4~--NO~3~) forms from dilute waters with dissolved calcium concentrations that are less than dissolved carbonate (Ca\<HCO~3~+CO~3~, equivalent%). A sulfate brine with near neutral *p*H (Na--K--Mg--Cl--**SO**~**4**~--NO~3~) forms from dilute waters with dissolved calcium concentrations that are greater than the dissolved carbonate, but less than the combined dissolved sulfate and carbonate concentrations (Ca\< SO~4~+HCO~3~, equivalent%). A calcium chloride brine with near neutral *p*H (Na--K--**Ca**--Mg--Cl--NO~3~) forms from dilute waters with a dissolved calcium concentration that is greater than the combined dissolved sulfate and carbonate concentrations (Ca\>SO~4~+HCO~3~, equivalent%). The measured compositions of Yucca Mountain pore water vary, but can be generally classified as waters that should evolve toward sulfate and sodium bicarbonate type brines, with a few calcium chloride brines as they evaporate \[Fig. [1(b)](#F1){ref-type="fig"}\].\[[@B6]-[@B8]\]

In Fig. [1](#F1){ref-type="fig"}, the simple ratios of calcium, sulfate, and carbonate illustrate the dominant carbonate and sulfate chemical divides that occur as waters evaporate.\[[@B2],[@B9],[@B10]\] However, Fig. [1](#F1){ref-type="fig"} does not show important chemical divides for magnesium, silica, or fluoride, nor does it show the relative amount of these salts to other major ions such as sodium, chloride, and nitrate. These are important chemical parameters, because they may mitigate or enhance the corrosiveness of brines at the waste package surfaces. Temperature will also impact the evolution of the pore water. Previous modeling of Yucca Mountain pore water initially within the sulfate brine field at 25°C evolved toward calcium chloride brines at 95°C, and not toward a sulfate brine as predicted by Fig. [1(a)](#F1){ref-type="fig"}.\[[@B7]\]

Thermodynamic equilibrium codes use a range of seepage water chemistry, temperature, and relative humidity to model the chemical environment on Yucca Mountain waste package surfaces.\[[@B1]\] This requires Pitzer parameters that can account for nonideal solutions at higher ionic strengths and elevated temperatures. Unfortunately, only a few relevant experiments are available to validate this modeling approach. The available data consist of evaporation of two bicarbonate type waters modeled after Yucca Mountain J-13 well water and a calcium chloride type water modeled after unsaturated pore water at Yucca Mountain and seawater.\[[@B1],[@B5]\]

In this paper we focus on the brine chemistry formed by the evaporation of a synthetic Yucca Mountain sulfate type pore water\[[@B8]\] at 95°C over a concentration range of 1× to \~3500×. This study provides additional benchmark data needed to both understand the brines and salts that form upon evaporation of various waters and to validate the EQ3/6 geochemical code and a high temperature Pitzer parameter thermodynamic database currently used by the Yucca Mountain Program to model aqueous chemical systems.

Methods
=======

A. Experimental methods
-----------------------

A synthetic pore water solution, based on HD-PERM-1 pore water composition,\[[@B8]\] was evaporated about 3400 times at 95°C in a series of experiments to iteratively concentrate the solution in manageable quantities of approximately 2 l and to monitor precipitation. The chemical composition of each successive leg was based on the brine composition toward the end of the previous leg (Table [I](#T1){ref-type="table"}). Experiments FEC 9, 12, 13, and 14 are referred to as legs 1, 2, 3, and 4 throughout this paper. The solutions were prepared at room temperature using analytical grade salts. Differences between the ending and starting compositions for each leg reflect the difficulty in exactly synthesizing the solutions and changes that occur when the solution is prepared at 25°C and then heated to the experimental temperature of 95°C. At the beginning of each new leg, the prepared solution did not include the undissolved solids that were present at the end of the previous leg. During the preparation of the starting solution for legs 3 and 4, an amorphous magnesium silicate precipitate formed in the 25°C solutions. This precipitate was not removed from the starting solutions and only partially dissolved when the solution was heated up to the experimental run temperature (95°C).

Evaporation was conducted in a vented, halar-lined vessel heated to 95°C in a fluidized sand bath furnace which provided optimal heat transfer for this method. The solution was stirred constantly and HEPA (high efficiency particulate air) filtered air was streamed over the solution to help control the evaporation rate. The solution vapor was refluxed to prevent evaporative water loss as it was heated to run temperature. Once the solution was at 95°C, the evolving water vapor was condensed into a separate container to monitor the extent of evaporation. Samples of the evaporating solution were periodically extracted and filtered at 95°C and analyzed to determine the water chemistry. The 95°C filtered samples withdrawn for cation and anion analysis were immediately diluted by directly injecting the sample into a known quantity of room temperature de-ionized water to prevent precipitation on cooling. Undiluted samples withdrawn for total dissolved inorganic carbon analysis, ∑CO~2~(aq), were immediately stored by filling gas tight vials to prevent equilibration with air at room temperature. Separate undiluted samples for solution *p*H were stored in a closed container and *p*H was measured as soon as they cooled to room temperature. In the last two samples of leg 4, precipitates formed in the *p*H and carbon samples as they cooled. Samples of the condensed water vapor were also periodically extracted and analyzed to monitor gas volatility. After the last sample was taken for each leg, the evaporation was continued to dryness. The solid precipitate was collected at the end of each leg of the experiment, dried in an oven at 40°C to facilitate sample preparation, and analyzed by powder x-ray diffraction (XRD).

B. Analytical methods
---------------------

Sample *p*H was measured at room temperature with a combination electrode, which is reliable in solutions with an ionic strength less than 0.1 molal.\[[@B11]\] The sample cooled to room temperature and *p*H was measured within a half hour of sampling. The *p*H measurements for legs 3 and 4, where the ionic strength of the solution was greater than 0.1 molal, represent uncorrected values and have not been corrected for ionic strength. Total dissolved carbon, ∑CO~2~(aq), was measured with an infrared carbon analyzer and had a detection limit of 1 ppm. Dissolved calcium, magnesium, silica, and sodium were measured with an inductively coupled plasma-atomic emission spectrometer, dissolved potassium was measured using an atomic absorption spectrophotometer, and fluoride, chloride, nitrate, and sulfate anions were determined using ion chromatography. Reproducibility of these techniques is typically better than ±2%. The mineralogical composition was determined by powdered XRD using a Cu *K*~*α*~source from 10° to 90° 2*θ*at 0.02° per step. The XRD instrument was calibrated using NIST traceable silicon (\# 640 c) and mica (\# 675) standards for high angle and low angle peaks, respectively. XRD cannot detect amorphous solids or minerals that are present at \<2 wt %. Mineral identification was based on the presence of the three most intense peaks in the XRD pattern for a given mineral. In some cases where the most intense peaks overlapped with other mineral peaks, identification was based on the presence of lower intensity diagnostic peaks.

C. Thermodynamic modeling calculations
--------------------------------------

Solution compositions were modeled using the EQ3/6 geochemical code, and a high temperature Pitzer ion-interaction model that is further described in Table [II](#T2){ref-type="table"}.\[[@B1],[@B12]-[@B14]\] The high temperature Pitzer ion-interaction model approximates nonideal behavior of solutions at elevated ionic strength and temperature. The predictive models were generated to mirror the experimental design and analysis, in which synthetic pore water was evaporated over a discrete range for each leg, with a cumulative evaporation up to 3500× for the overall experiment. The concentration factor, CF, can then be defined as

![](1467-4866-6-31-i9.gif)

Where H~2~O~(*i*)~is the initial mass of H~2~O solvent and H~2~O~(*n*)~is the mass of H~2~O solvent remaining after the *n*th step in the evaporation process.

The evaporation model consisted of three steps. In the first step, the measured composition of the first sample at 95°C (Table [I](#T1){ref-type="table"}, FEC\#-1) was speciated, suppressing all mineral precipitation in the calculation. At this point in the experiment, the solution was simply brought up to temperature, refluxing any water vapor to prevent evaporation. During the second step, the speciated water was evaporated by stepwise removal of solvent water at a fixed rate of 0.25 mol H~2~O reactant per mol of solute at 95°C. In the evaporation step, all minerals were allowed to precipitate with the exception of quartz and dolomite because of known slow kinetics; glaserite, hydromagnesite, and magnesite because the available thermodynamic data are questionable; and cristobalite because it forms above 1470°C and is not relevant to this experiment. In the speciation and evaporation steps, oxygen and carbon dioxide fugacities were fixed at 21% and 0.033%, respectively, to simulate atmospheric experimental conditions. A CO~2~(g)-sink was added to the model to remove excess buildup of CO~2~(g) from the reaction surface. Finally, in the third step, the predicted *p*H values at 95°C were recalculated to 25°C to compare with measured *p*H at room temperature. This was achieved by performing a further calculation that reduced the temperature of the reaction from 95 to 25°C, while fixing ∑CO~2~(aq) to the predicted 95°C value and suppressing mineral formation. For all calculations, electrical balance was achieved by automatically changing the sodium concentration with a convergence tolerance of 0.1 ppb. Charge balancing was necessary due to analytical errors generated in the experimental analysis and also potentially incomplete analysis. The calculated activity of water was interpreted as a function of relative humidity, and predicted solution composition, *p*H, and mineral composition were compared with experimental ion analysis and XRPD results with respect to the overall concentration factor. The results are shown graphically in Figs. [3](#F3){ref-type="fig"} and [4](#F4){ref-type="fig"}, and discussed in the following sections.

Results
=======

Evaporation of the dilute water yielded a sulfate brine as predicted by chemical divide theory based on its initial Ca:SO~4~:HCO~3~ratio (Fig. [2](#F2){ref-type="fig"}). Although this water is classified as a sulfate brine, sulfate concentrations are minor compared to the concentration of sodium and chloride, which dominate the solution chemistry. At the conclusion of leg 4, sodium and chloride were 45 and 40 mol %, respectively, while calcium and sulfate were relatively minor constituents at 0.1 and 0.8 mol %, respectively (Table [III](#T3){ref-type="table"}). The initial experimental solution contained trace amounts of fluoride that were rapidly removed from solution presumably as highly insoluble fluorite (CaF~2~) (Table [I](#T1){ref-type="table"}). Fluoride was not included in the model since it was not detected in the first sample analysis at 95°C (leg 1, FEC9-1). The minerals identified by XRD in the precipitates are halite (NaCl), bassanite (2CaSO~4~·H~2~O), anhydrite (CaSO~4~), niter (KNO~3~), and nitratine (NaNO~3~), and are listed for each leg in Table [IV](#T4){ref-type="table"}.

We see no evidence of volatility for HC1, HNO~3~, and HF gases in these experiments. Concentrations of fluoride, chloride, nitrate, and sulfate in the condensed vapor were all below the detection limits. This is in contrast to evaporation of a concentrated calcium chloride type water at around 140°C (based on the 1000× solution results from Rosenberg *et al*.\[[@B5]\]) where significant volatilization of HCl(g) was measured by acidic condensates at 90% evaporation at \~75 000×.\[[@B15]\] While our evaporation is less than the aforementioned research, our results indicate that gas volatility is not a major concern for the evaporation and concentration of sulfate waters at 95°C and \~3 400×.

Figures [3](#F3){ref-type="fig"} and [4](#F4){ref-type="fig"} show the experimental and predicted solution composition and the predicted mineral precipitation as a function of overall concentration factor. There is excellent agreement between the model predictions of potassium and nitrate concentrations with those measured by experiment. Solution data show conservative concentration of both potassium and nitrate in each evaporation leg, indicating no mineral precipitation \[Figs. [3(a)](#F3){ref-type="fig"} and [4](#F4){ref-type="fig"}\]. This is supported by the XRD data where only a small amount of niter was identified in the last leg after the solution had completely evaporated when precipitation of all mineral phases is expected (Table [IV](#T4){ref-type="table"}). Calculations required the suppression of pentasalt (gorgeyite, K~2~Ca~5~(SO~4~)~6~·H~2~O) precipitation at a concentration factor above 1000× in leg 4 \[Fig. [3(a)](#F3){ref-type="fig"}\] to achieve agreement with experimental solution composition and solid characterization. Pentasalt was not detected by XRD analysis.

For sodium and chloride behavior, we observe good agreement between experimental composition and model prediction \[Figs. [3(b)](#F3){ref-type="fig"} and [4](#F4){ref-type="fig"}\]. Both ions concentrate in solution with increasing evaporation until the solution is saturated with respect to halite. Halite is the dominant salt in samples taken to dryness and identified by XRD (Table [IV](#T4){ref-type="table"}). There is a discrepancy between the experimental sodium concentration and that predicted by the model in leg 1, due to the charge balance correction using sodium ions in the model. This observation is explained by the overprediction of positively charged magnesium, and the subsequent decrease in calculated sodium concentration to neutralize the charge discrepancy. The charge balance correction resulted in a sodium concentration that differed from the original (measured) concentration by 21% in the first leg, 13% in the second leg, 1% in the third leg, and 5% in the fourth leg. Since sodium saturation did not occur until the formation of halite late in fourth leg, we believe this sodium correction did not significantly affect the quality of our calculations.

There is also good agreement between experimental compositions and model prediction for calcium and sulfate \[Figs. [3(c)](#F3){ref-type="fig"} and [4](#F4){ref-type="fig"}\]. At a concentration factor of roughly 10 × in leg 2, both calcium and sulfate begin to precipitate as can be seen in the decrease in their slopes. Although calcium sulfate precipitation continues over the duration of the experiment, dissolved calcium decreases as the sulfate increases with continued evaporation at a concentration factor of about 1000×. This is consistent with the chemical divide theory and the initial composition of the water which contained SO~4~: Ca\> 1. This behavior was best modeled by using bassanite as the solubility limiting phase (using gypsum or anhydrite resulted in an overprediction and underprediction, respectively, of the calcium and sulfate concentrations). This finding is in partial agreement with the experiments, which identified both anhydrite and bassanite (Table [IV](#T4){ref-type="table"}). This difference may be an artifact of the experimental protocol because calcium sulfate hydration states can be readily altered by changes in temperature and humidity, such as those found in the drying and preparation of the precipitate prior to XKD analysis. The model overestimates calcium and sulfate by about a factor of 2 in concentrated brines (1000× concentration factor), and appears to increase with continued evaporation. It is possible that the overprediction in calcium and sulfate concentrations observed using bassanite and gypsum solubility controls, and the underprediction observed with anhydrite controls indicates a metastable mixture of these solubility limiting phases.

Comparison between both magnesium and silica experimental concentrations and model predictions during the evaporation show reasonable agreement in legs 2 and 3 and only fair agreement in legs 1 and 4 \[Fig. [3(f)](#F3){ref-type="fig"}\]. Experiment and prediction both show that magnesium and silica are removed from solution as solid precipitates. The model predicts that sepiolite, amorphous silica, and brucite are the solubility controls (Fig. [4](#F4){ref-type="fig"}). These phases were not observed in the experiment possibly because the amount was too small to be detected or because they were amorphous. It is possible that a noncrystalline magnesium silicate phase precipitated, similar to the solid that formed at 25°C (Legs 3 and 4).\[[@B16]\] The Pitzer database contained only two magnesium silicate minerals, talc (Mg~3~Si~4~O~10~(OH)~2~) and sepiolite (Mg~4~Si~6~O~15~(OH)~2~·6H~2~O). A better model fit was achieved when sepiolite was allowed to precipitate and talc formation was suppressed.

There is good agreement between the experimental and predicted total dissolved carbonate in the final leg \[shown in Fig. [3(d)](#F3){ref-type="fig"} as ∑CO~2~(aq)\], assuming equilibrium with atmospheric CO~2~(g) at temperature. In more dilute legs 1 to 3, total dissolved carbonate was not detected after the first few samples for each leg. Predicted carbonate concentrations for legs 1, 2, and 3 are all below the detection limit and are consistent with the experiment. The failure of the model to capture the initial dissolved carbonate concentrations suggests that ∑CO~2~(aq) in starting solution synthesized at 25°C had not degassed to the lower equilibrium amount at 95°C. Model predictions show that carbonate concentrations decrease throughout the evaporation process as carbonate is lost to the atmosphere as gaseous CO~2~in conjunction with a decreasing *p*H. No carbonate minerals were predicted to form until leg 4, where a very small amount of calcite precipitates, which was too small to be detected by XRD.

Model prediction of *p*H in the first two legs is in reasonable agreement with experimental data showing that the *p*H values decrease during the evaporation in each leg \[Fig. [3(e)](#F3){ref-type="fig"}\]. However, in legs 3 and 4 where the ionic strength exceeds 0.1 molal, measured *p*H values are uncorrected and are as much as 2 *p*H units lower than the predicted *p*H. The measured values were not corrected for ionic strength effects in these complex solutions. Rai and Felmy\[[@B17]\] report that measured *p*H will be lower than the actual *p*H by 0.14 units in 1 molal NaCl and by 0.97 0.97 *p*H in 6 molal NaCl due to the ionic strength effects on the liquid junction potential of a commercially available 3 M KC1 combination electrode similar to that used in these experiments. Even larger discrepancies between measured and real *p*H values are seen in more complex systems containing mixtures of mono and divalent ions at high ionic strength. The difference in measured *p*H at the end of one leg and the start of the next leg is an artifact of the experimental protocol. The waters were synthesized at room temperature, equilibrated with atmospheric CO~2~(g) and have higher dissolved carbon and higher *p*H than they possess at 95°C.

Accurate prediction of *p*H is very important because the solubility of many solid phases, such as sepiolite, calcite, and amorphous silica are strongly influenced by solution *p*H. Therefore, it is possible that the discrepancy between experimental and predicted values of magnesium, silica, and calcium reflect an underprediction of solution *p*H, although the *p*H in our models is below that which would affect amorphous silica. Unfortunately, in our study, *p*H is one of the most difficult parameters to predict because values at elevated temperature must be extrapolated to 25°C to compare with the measured values. The measured and predicted *p*H are then subject to change due to possible mineral precipitation and equilibration with atmospheric CO~2~(g) at room temperature. Furthermore, measured *p*H values in concentrated solutions are unconnected and do not represent real H^+^activity. We have minimized the contribution of *p*H uncertainty to the observed discrepancy by constraining predicted *p*H by fitting the ∑CO~2~(aq) concentration in leg 4, where we observed measurable concentrations. This yields a solution in equilibrium with atmospheric CO~2~and an initial *p*H of 7.5 at 95°C. This is consistent with the experiment, because filtered laboratory air was continually passed over the solution as the waters evaporated, and because carbonate samples were stored in gas-tight vials eliminating exchange with atmospheric CO~2~at room temperature. We also assume that the undiluted, sealed samples taken for *p*H measurement did not re-equilibrate with atmospheric CO~2~as they cooled from 95°C to room temperature. The good agreement between measured and predicted *p*H at lower ionic strengths in legs 1 and 2 supports these modeling constraints.

Discussion
==========

A. Chemical divides
-------------------

The important chemical divides that control the composition of brines formed from dilute sulfate type waters are halite, bassanite (or other calcium sulfates), magnesium silicate, amorphous silica, and possibly fluorite and brucite based on experimental results and model predictions. The early removal of fluoride from the starting solution is an important geochemical control for eliminating the evolution of a potentially corrosive fluoride containing brine. The precipitation of calcium as a calcium carbonate is not a major chemical divide for this solution. At 95°C and atmospheric CO~2~(g), carbonate is partitioned into the gas phase rather than the precipitation of calcium carbonate. The solution is below calcite solubility for most of the experiment. The precipitation of calcium as bassanite and magnesium as a magnesium silicate are important geochemical controls that limit the calcium and magnesium content in these brines. Additionally, the very low fluoride solubility limits chloride to be the most corrosive agent of Yucca Mountain sulfate type pore waters. High nitrate to chloride ratios of brines are known to limit susceptibility to localized corrosion of corrosion resistant materials such as the candidate waste package material.\[[@B18],[@B19]\] The brine contained a nitrate to chloride ratio of 0.2:1 at 99.97% evaporation. This ratio will increase with increasing evaporation because the chloride will be controlled by halite solubility, and nitrate will continue to concentrate until the solution reaches saturation with respect to nitratine (NaNO~3~) and/or niter (KNO~3~).

Although the Ca:SO~4~:HCO~3~ternary diagrams do not capture all of the important chemical divides that affect the composition of Yucca Mountain pore waters, they can be used to categorize the types of brines that will form from the wide range of Yucca Mountain pore waters. Evaporation of dilute sulfate (this study), bicarbonate and calcium chloride\[[@B5]\] type Yucca Mountain waters evolve toward their respective sulfate, carbonate, and calcium chloride brines indicated by their initial Ca:SO~4~:HCO~3~ratios (Fig. [2](#F2){ref-type="fig"}).

B. High temperature Pitzer model
--------------------------------

The comparison of experimental results and model predictions of the evaporation of synthetic Topopah Spring Tuff pore water at 95°C indicates that the current high-temperature Pitzer database used by the Yucca Mountain program adequately describes the chemical evolution of brines at elevated temperature for most species (Fig. [3](#F3){ref-type="fig"}). In this section we discuss our results in light of the high temperature Pitzer ion interaction database that includes solubility products for solids. Table [II](#T2){ref-type="table"} lists temperature-dependent Pitzer interaction parameters for binary and ternary reactions and Table [V](#T5){ref-type="table"} details the solubility products relevant to this work.

The Yucca Mountain Project high temperature Pitzer ion interaction database is the most comprehensive database available to account for the nonideal behavior of highly concentrated electrolytes over a wide range of temperature (0--140°C). The database was founded on the original variable-temperature Pitzer parameters\[[@B20],[@B21]\] supplemented by parameter data from several other sources.\[[@B13],[@B22]-[@B32]\] It also includes thermodynamic parameters converted from non-standard Pitzer equations from the published literature.\[[@B14]\]Temperature-independent parameters based on 25°C data are used for several parameters where temperature dependent data are lacking.\[[@B13]\] The database contains temperature-dependent ion interaction parameters for most ion groups relevant to our experimental system at 95°C. Exceptions include 25°C models for potassium nitrate interactions and some calcium and magnesium ion interactions;\[[@B13]\] and 20--90°C models for CO~2~(aq) ion interactions.\[[@B26]\] Substantial database and model validation has been performed,\[[@B1]\] however, it is acknowledged that the results of any model are only as good as the input parameters and database used, and typically, improvements will always be made to bridge the gaps between experimental observations and model predictions.

### 1. Na^+^, Cl^-^, K^+^, and N![](1467-4866-6-31-i5.gif)

The excellent agreement between model prediction and experiment for sodium, chloride, potassium, and nitrate suggests that the high temperature Pitzer ion interaction database and halite solubility product adequately describes the nonideal solution chemistry at high ionic strength and elevated temperature for these elements, as well as the halite solubility product. The agreement is expected for sodium and chloride because interaction parameters for binary and ternary ion groups are well defined as a function of temperature.\[[@B14],[@B21],[@B25],[@B27],[@B33]-[@B35]\]

The excellent agreement between prediction and experiment for potassium and nitrate suggests that Pitzer parameters for the K^+^--N![](1467-4866-6-31-i5.gif) ion interaction will yield accurate predictions of pore water concentrations in evaporating solutions despite being measured at 25°C and applied to 95°C systems. The reason for this is that high sodium to potassium, and chloride to nitrate mole ratios seen in our experiments may effectively mask any mismatch due to constant 25°C parameters. Clearly, this is true for undersaturated brines with respect to KNO~3~(niter) because Na^+^--N![](1467-4866-6-31-i5.gif) and K^+^--C1^-^interactions will be more important than K^+^--N![](1467-4866-6-31-i5.gif) interactions in determining solution properties. However, when brines are saturated with niter (KNO~3~), then K^+^--N![](1467-4866-6-31-i5.gif) interactions are more important. The extent that the constant temperature K^+^--N![](1467-4866-6-31-i5.gif) interaction parameter does not accurately predict behavior was shown in KNO~3~--NaNO~3~deliquescence experiments at 90°C. In solutions saturated with respect to niter/nitratine, they showed an underprediction of niter solubility by as much as 50%.\[[@B35]\]

### 2. Ca^2+^and S![](1467-4866-6-31-i10.gif)

The discrepancy between experimental results and model prediction for dissolved calcium and sulfate concentrations observed in the high ionic strength solutions of leg 4 may be due to either the temperature-dependent Ca^2+^--S![](1467-4866-6-31-i10.gif) Pitzer parameters or an uncertainty in the physical nature of the calcium sulfate precipitate formed during the course of the experiment. The parameters in the Yucca Mountain database are derived from a refit of the Moller,\[[@B20]\] and Greenberg and Moller\[[@B21]\] models, but without an explicit temperature dependent CaSO~4~(aq) ion pair. To examine the ability of the high temperature Pitzer parameters to accurately predict calcium sulfate solubility, we compare model calculations to literature anhydrite and bassanite solubility in solutions of varying sodium chloride concentrations.\[[@B36],[@B37]\] The results are presented in Fig. [5](#F5){ref-type="fig"} and show agreement within 0.01 molal between predicted and measured anhydrite and bassanite solubility, which is similar to the agreement between measured and predicted calcium in our evaporation experiment, but it is as much as seven times lower than the measured and predicted sulfate in our evaporation experiments.

The exact nature of the calcium sulfate solid phase formed in the experiment and controlling the solution concentrations of both calcium and sulfate may also contribute to the discrepancy between evaporation models and experiment. Gypsum changes rapidly to bassanite in contact with water at 97.5°C\[[@B38]\] while bassanite is considered metastable and in turn is converted slowly to anhydrite at similar temperatures.\[[@B39]\] The gypsum--anhydrite transition temperature of 42°C is known to increase in the presence of additional salts in solution. Recall that both bassanite and anhydrite were identified by XRD at the end of the experimental evaporation. Clearly, our experimental conditions are such that gypsum--bassanite--anhydrite transitions may occur both during the evaporation experiment at 95°C and during the drying of the final precipitates at 40°C prior to XRD analysis. Our experimental calcium concentrations during the evaporation in leg 4 generally fall between those predicted by the presence of bassanite and anhydrite \[Fig. [5(b)](#F5){ref-type="fig"}\]. The CaSO~4~solubility comparison shown in Fig. [5(a)](#F5){ref-type="fig"} is validated only to 6 molal, prior to the precipitation of halite. A theoretical solid solution composition with end members resembling both bassanite and anhydrite may be more appropriate in describing the nature of our evaporation solids.

### 3. ∑CO~2~(aq) and *p*H

The prediction of solution *p*H is important because *p*H greatly affects the solubility of many carbonate and silicate phases. The Pitzer database contains parameters for several ion interactions that were derived from highly acidic or alkaline solutions (0 to \> 6 molal) that are not representative of the *p*H of natural solutions or brines (4 \<*p*H \<12). In natural systems, the *p*H is strongly tied to the partial pressure of CO~2~(g), as well as the dissociation of CO~2~(aq) to HC![](1467-4866-6-31-i5.gif) and C![](1467-4866-6-31-i11.gif). He and Morse\[[@B26]\] derived Pitzer parameters for ion interactions involving CO~2~(aq), HC![](1467-4866-6-31-i5.gif) and C![](1467-4866-6-31-i11.gif) in strong electrolytes containing sodium, potassium, calcium, magnesium, chloride, and sulfate from calcite solubility experiments from 0 to 90°C. While the parameters based on He and Morse\[[@B26]\] are incorporated into the Yucca Mountain Pitzer database, the Pitzer parameter temperature function embedded in the EQ3/6 code is not always consistent with published Pitzer data (e.g., too few or too many Pitzer parameters). Consequently, Pitzer data are required to be refit to be internally consistent with both the EQ3/6 code. As stated earlier, we measured only uncorrected *p*H values in our evaporation study and most measured carbonate concentrations fell below the analytical limit of detection.

To evaluate our ability to model H~2~O--CO~2~systems, we utilized the equilibria between calcite, dissolved carbonate, and *p*H. We compare model calculations with known literature calcite solubility data approached from undersaturation in 1 to 5 molal sodium chloride brines (with lesser concentrations of potassium, calcium, and sulfate) from 25 to 90°C.\[[@B40],[@B41]\] We do not compare similar calcite seeded experiments approached from supersaturation because He and Morse\[[@B26]\] observed higher solubility, which they attributed to cation substitution in the precipitated calcite. The results are shown in Fig. [6](#F6){ref-type="fig"}. At higher temperatures approaching those of our evaporation study, model calculations match calcite solubility to within 0.003--0.004 molal. Thus model calculations of *p*H and ∑CO~2~(aq) are reasonable reflections of the experimental equilibrium.

### 4. Mg and Si

The factor of 2 discrepancy between model and experiment for magnesium and silica solubility translates to absolute concentrations on the order of 0.01 molal for magnesium and only 10^-5^molal for silica. This agreement suggests that the use of sepiolite and brucite (Mg(OH)~2~) to model magnesium and silica solubility, and the high temperature Pitzer parameters are adequate for the evaporation of waters in a repository environment. Since precipitates formed in the synthesis of the initial solution were not carried over to evaporation studies, formation (and loss through precipitation) of magnesium and silica solids could result in a discrepancy between model and experimental results. Uncertainty in the sepiolite solubility product is well defined at temperatures less than 100°C, with log *K*~sp~= 24 (± 0.3) at 95°C (Mg~4~Si~6~O~15~(OH)~2~·4H~2~O + 8H^+^= 4 Mg^2+^+ 6SiO~2~(aq) + 11H~2~O).\[[@B42]\] This translates to less than 10% uncertainty. Sepiolite has been observed as a likely control for the removal of magnesium and/or silica in other systems such as at Saline Valley, CA that have a brine evolution similar to this experimental pore water.\[[@B43]\] It is doubtful that an amorphous magnesium silicate forms, because the model slightly overpredicts magnesium and silica concentrations. The observed match between model and experiment also suggests that the suppression of talc is valid because its formation is thought to be kinetically unfavorable. However, additional studies are needed to identify the magnesium silicate phase that controls solubility.

The absolute overprediction of silica solubility is well within the uncertainty of Pitzer parameters in the high temperature database. We evaluate the Mg^2+^--SiO~2~(aq) Pitzer parameters by comparing model calculations of SiO~2~(am) solubility data in magnesium and sodium containing electrolytes at 100°C in Fig. [7](#F7){ref-type="fig"}\[[@B44],[@B45]\] Mg^2+^--SiO~2~(aq) Pitzer parameters were initially derived from amorphous silica solubility measured in a range of electrolyte solutions \[[@B32],[@B44]-[@B46]\] in which parameter accuracy is limited by low silica concentrations and uncertainty in the chemical potential for amorphous silica.\[[@B32]\] Model calculations in systems in which amorphous silica controls solubility (sepiolite is suppressed) slightly overpredict silica solubility by 0.001 molal on average. If sepiolite and amorphous silica are allowed to control silica solubility, then the scatter in the data noticeably increases to an average of 0.003 molal. These values are more than two orders of magnitude higher than the overprediction observed in our evaporation study.

V. Conclusion
=============

The assessment of water chemistry that may contact the waste containers and drip shields at the Yucca Mountain high-level radioactive waste repository in Nevada requires the calculation of a wide range of water composition over 10 000 year time periods in an environment in which temperature and relative humidity will change as the repository heats up and cools down. Understanding and verifying the chemical divides that control brine composition as seepage water evaporates is important to placing boundaries on the corrosiveness of the chemical environment to waste package materials. The model data generated for this solution by EQ3/ 6 geochemical code and high temperature Pitzer database indicate that they can provide satisfactory predictions as compared to the experimental data in most cases. The water evolved toward a complex sulfate type brine that contained about 45 mol % Na, 40 mol % Cl, 9 mol % NO~3~, 5 mol % K, and less than 1 mol % each of SO~4~, Ca, Mg, ∑CO~2~(aq), F, and Si at a concentration factor of about 3500×. Minerals predicted to form include halite, anhydrite, bassanite, niter, and nitratine in addition to fluoride, carbonate, sulfate, and magnesium silicate precipitates. This work is of importance to the continued validation of the Pitzer ion-interaction parameter database and geochemical modeling used by the Yucca Mountain Project, and supports its adequate capabilities in predicting brine evolution in complex aqueous systems at elevated temperatures.

Acknowledgements
================

The authors would like to thank Brian Viani and Tom Wolery at LLNL for their guidance and discussion regarding x-ray diffraction data and thermodynamic modeling, respectively, and Kirk Staggs for his invaluable skills in making these experiments possible. This work was performed under the auspices of the U.S. Department of Energy by the University of California Lawrence Livermore National Laboratory under Contract No. W-7405-Eng-48.

Figures and Tables
==================

![(Color) (a) Calcium carbonate and calcium sulfate chemical divides for evaporation of dilute waters. (b) Yucca Mountain, NV pore waters as measured.](1467-4866-6-31-1){#F1}

![(Color) Chemical evolution of dilute calcium chloride (synthetic Topopah Spring tuff porewater -- Ref. 5), Na-bicarbonate (synthetic J-13 groundwater -- Ref. 5), and sulfate (synthetic Topopah Spring tuff porewater, this study) waters upon evaporation.](1467-4866-6-31-2){#F2}

![(Color) Evaporation of dilute sulfate water based on a Topopah Spring tuff porewater chemistry. Comparison of experimental and model solution concentrations vs concentration factor. Symbols indicate experimental data and lines indicate model data.](1467-4866-6-31-3){#F3}

![(Color) Evaporation of dilute sulfate water based on a Topopah Spring tuff porewater chemistry. Predicted mineral amounts vs concentration factor.](1467-4866-6-31-4){#F4}

![(Color) (a) Comparison of literature and calculated anhydrite and bassanite solubility data. Literature values were only obtained up to 6 molal sodium, therefore our calcium sulfate validation is limited to this upper value, (b) Comparison of modeled calcium and sulfate concentrations using bassanite, anhydrite, and gypsum as the solubility controlling mineral phase and leg 4 calcium and sulfate concentrations prior to modeled halite formation.](1467-4866-6-31-5){#F5}

![(Color) Comparison of model data to literature calcite solubility data approached from undersaturation in synthetic brines (Ref. 26). Lines represent LLNL model and symbols represent (Ref. 40) and Lervson (Ref. 41) experimental data.](1467-4866-6-31-6){#F6}

![(Color) Comparison of silica solubility data in MgCl~2~and NaCl electrolytes at 100°C (Ref. 46) with predicted model calculations using the high temperature Pitzer ion interaction data base. Note that the smaller number of data plotted in the calculations, where both sepiolite and amorphous silica are allowed to control solubility, reflects the MgCl~2~electrolytes only.](1467-4866-6-31-7){#F7}

###### 

Starting and final compositions for the evaporation of a synthetic \"sulfate type\" Topopah Spring Tuff pore water.

           Leg l FEC 9 molal^a^starting   Leg l FEC 9 molal, 95°C starting   Leg l FEC 9 molal, 95°C final   Leg 2 FEC 12 molal, 95°C starting   Leg 2 FEC 12 molal, 95°C final   Leg 3^b^FEC 13 molal, 95°C starting   Leg 3 FEC 13 molal, 95°C final   Leg 4^b^FEC 14 molal, 95°C starting   Leg 4 FEC 14 molal, 95°C final
  -------- ------------------------------ ---------------------------------- ------------------------------- ----------------------------------- -------------------------------- ------------------------------------- -------------------------------- ------------------------------------- --------------------------------
  PH^a^    9.65                           na^c^                              7.15                            7.60                                7.32                             7.69                                  na^c^                            8.08                                  6.97
  F        5.94 × 10^-5^                  nd^d^                              nd^d^                           nd^d^                               nd^d^                            nd^d^                                 nd^d^                            nd^d^                                 nd^d^
  Ca       1.57 × 10^-3^                  1.25 × 10^-3^                      4.36 × 10^-3^                   4.69 × 10^-3^                       2.25 × 10^-2^                    2.17 × 10^-2^                         4.21 × 10^-2^                    4.06 × 10^-2^                         1.03 × 10^-2^
  Mg       1.01 × 10^-3^                  6.48 × 10^-4^                      5.58 × 10^-4^                   6.05 × 10^-4^                       4.52 × 10^-3^                    1.06 × 10^-3^                         6.43 × 10^-3^                    1.88 × 10^-3^                         9.02 × 10^-3^
  Na       2.84 × 10^-3^                  3.04 × 10^-3^                      9.23 × 10^-3^                   1.10 × 10^-2^                       1.05 × 10^-1^                    1.11 × 10^-1^                         1.15                             1.12                                  7.42
  SiO~2~   1.39 × 10^-3^                  9.21 × 10^-4^                      1.15 × 10^-3^                   1.38 × 10^-3^                       7.91 × 10^-3^                    2.24 × 10^-3^                         6.32 × 10^-3^                    3.09 × 10^-4^                         1.89 × 10^-5^
  K        1.76 × 10^-4^                  1.63 × 10^-4^                      5.15 × 10^-4^                   5.82 × 10^-4^                       5.66 × 10^-3^                    6.54 × 10^-3^                         6.99 × 10^-2^                    7.89 × 10^-2^                         8.09 × 10^-1^
  Cl       3.26 × 10^-3^                  3.31 × 10^-3^                      1.00 × 10^-2^                   9.80 × 10^-3^                       9.80 × 10^-2^                    1.05 × 10^-1^                         9.93 × 10^-1^                    1.09                                  6.59
  NO~3~    3.53 × 10^-4^                  3.94 × 10^-4^                      1.15 × 10^-3^                   1.32 × 10^-3^                       1.12 × 10^-2^                    1.18 × 10^-2^                         1.29 × 10^-1^                    1.32 × 10^-1^                         1.42
  SO~4~    1.23 × 10^-3^                  1.85 × 10^-3^                      4.54 × 10^-3^                   4.64 × 10^-3^                       2.35 × 10^-2^                    2.24 × 10^-2^                         5.85 × 10^-2^                    5.80 × 10^-2^                         1.36 × 10^-1^
  HCO~3~   7.93 × 10^-4^                  3.11 × 10^-4^                      nd^d^                           1.59 × 10^-4^                       nd^d^                            1.92 × 10^-4^                         nd^d^                            1.95 × 10^-4^                         2.66 × 10^-4^

^a^Measured at room temperature.

^b^Magnesium and silica precipitated from the solution at 25°C. Initial gravimetric concentrations are Mg = 4.5 × 10^-3^mol/kg-solution and Si = 7.8 × 10^-3^mol/kg-solution for Leg3 and are Mg = 6.0 × 10^-3^mol/kg-solution and Si = 5.9 × 10^-3^mol/kg-solution for Leg4.

^c^Not analyzed.

^d^Not detected, detection limits: F = 0.25 ppm, HCO~3~= 1 ppm.

###### 

Temperature dependence of Pitzer interaction parameters. 25°C-centric equation used to derive parameters: χ~(*T*)~= *a*~1~+ *a*~2~\[{1/*T*}-{l/298.15}\] + *a*~3~In{*T*/298.15} + *a*~4~{*T*- 298.15}.

                              Pitzer interaction parameters                                                                            
  ------------- ------------- ------------------------------- ----------------- ----------------- ----------------- ------------------ -----------------------------
  Ca^++^        Cl^-^                                         4.46 × 10^-1^     2.21 × 10^2^      1.61 × 10^-11^    2.28 × 10^-4^      Sterner *et al*. \[34\]
                                                              -1.66 × 10        -8.83 × 10^3^     7.10 × 10^-11^    -2.49 × 10^-2^     
                                                              0.00              0.00              0.00              0.00               
                                                              -1.73 × 10^-2^    -1.30 × 10^1^     -3.87 × 10^-13^   -3.15 × 10^-5^     
  Ca^++^        HC                                            4.00 × 10^-1^     0.00              0.00              0.00               Pitzer \[13\]
                                                              2.98              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
  Ca^++^        HS                                            2.15 × 10^-1^     0.00              0.00              0.00               Pitzer \[13\]
                                                              2.53              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
  Ca^++^        N                                             1.48 × 10^-1^     -4.88 × 10^1^     -7.47 × 10^-2^    -1.70 × 10^-4^     Oakes and Felmy \[28\]
                                                              2.44              -2.24 × 10^4^     -9.93 × 10^1^     1.19 × 10^-1^      
                                                              0.00              0.00              0.00              0.00               
                                                              -4.12 × 10^-3^    -1.79 × 10^1^     -1.19 × 10^-1^    1.87 × 10^-4^      
  Ca^++^        S                                             1.50 × 10^-1^     0.00              0.00              0.00               Greenberg and Moller \[21\]
                                                              3.00              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
  H^+^          Cl^-^                                         1.77 × 10^-1^     -3.35 × 10^1^     -2.62 × 10^-1^    1.26 × 10^-4^      Holmes *et al*. \[27\]
                                                              2.93 × 10^-1^     3.40 × 10^3^      1.98 × 10^1^      -2.79 × 10^-2^     
                                                              0.00              0.00              0.00              0.00               
                                                              3.62 × 10^-4^     -2.91 × 10^-11^   0.00              -3.04 × 10^-5^     
  H^+^          HS                                            2.09 × 10^-1^     1.05 × 10^3^      5.96              -8.78 × 10^-3^     Holmes and Mesmer \[24\]
                                                              4.41 × 10^-1^     2.96 × 10^2^      2.37              -4.63 × 10^-3^     
                                                              0.00              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
  H^+^          N                                             1.26 × 10^-1^     5.60 × 10^2^      4.92              -9.95 × 10^-3^     Felmy *et al*. \[32\]
                                                              2.88 × 10^-1^     2.71 × 10^-7^     1.59 × 10^-9^     1.34 × 10^-3^      
                                                              0.00              0.00              0.00              0.00               
                                                              -5.60 × 10^-3^    -6.58             -5.77 × 10^-2^    1.10 × 10^-4^      
  H^+^          S                                             9.86 × 10^-2^     -6.70 × 10^3^     -4.17 × 10^1^     6.37 × 10^-2^      Pabalan and Pitzer \[47\]
                                                              0.00              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              5.93 × 10^-2^     -3.03 × 10^3^     -1.61 × 10^1^     2.02 × 10^-2^      
  K^+^          Cl^-^                                         4.78 × 10^-2^     -3.43 × 10^2^     -1.38             1.34 × 10^-3^      Greenberg and Moller \[21\]
                                                              2.16 × 10^-1^     -5.76 × 10^2^     -2.88             4.64 × 10^-3^      
                                                              0.00              0.00              0.00              0.00               
                                                              -7.49 × 10^-4^    3.65 × 10^1^      1.48 × 10^-1^     -1.47 × 10^-4^     
  K^+^          C                                             1.29 × 10^-1^     0.00              0.00              0.00               Pitzer \[13\]
                                                              1.43              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              4.99 × 10^-4^     0.00              0.00              0.00               
  K^+^          F^-^                                          8.09 × 10^-2^     0.00              0.00              0.00               Pitzer \[13\]
                                                              2.02 × 10^-1^     0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              9.30 × 10^-4^     0.00              0.00              0.00               
  K^+^          HC                                            -1.07 × 10^-2^    0.00              0.00              0.00               Pitzer \[13\]
                                                              4.78 × 10^-2^     0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
  K^+^          HS                                            -3.00 × 10^-4^    0.00              0.00              0.00               Pitzer \[13\]
                                                              1.74 × 10^-1^     0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
  K^+^          N                                             -8.16 × 10^-2^    0.00              0.00              0.00               Pitzer \[13\]
                                                              4.94 × 10^-2^     0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              6.60 × 10^-3^     0.00              0.00              0.00               
  K^+^          OH^-^                                         1.30 × 10^-1^     0.00              0.00              0.00               Pitzer \[13\]
                                                              3.20 × 10^-1^     0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              4.10 × 10^-3^     0.00              0.00              0.00               
  K^+^          S                                             5.55 × 10^-2^     -1.42 × 10^3^     -6.75             8.27 × 10^-3^      Greenberg and Moller \[21\]
                                                              7.96 × 10^-1^     2.07 × 10^3^      2.33 × 10^-10^    2.36 × 10^-2^      
                                                              0.00              0.00              0.00              0.00               
                                                              -1.88 × 10^-2^    0.00              9.09 × 10^-13^    -2.66 × 10^-15^    
  Mg^++^        Cl^-^                                         3.51 × 10^-1^     -6.56 × 10^1^     -5.25 × 10^-1^    4.47 × 10^-4^      Pabalan and Pitzer \[47\]
                                                              1.65              -2.87 × 10^3^     -2.30 × 10^1^     4.95 × 10^-2^      
                                                              0.00              0.00              0.00              0.00               
                                                              6.53 × 10^-3^     -2.67 × 10^1^     -2.14 × 10^-1^    3.11 × 10^-4^      
  Mg^++^        HC                                            3.30 × 10^-2^     0.00              0.00              0.00               Pitzer \[13\]
                                                              8.50 × 10^-1^     0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
  Mg^++^        HS                                            4.75 × 10^-1^     0.00              0.00              0.00               Pitzer \[13\]
                                                              1.73              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
  Mg^++^        N                                             3.67 × 10^-1^     0.00              0.00              0.00               Pitzer \[13\]
                                                              1.58              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              -2.06 × 10^-2^    0.00              0.00              0.00               
  Mg^++^        S                                             2.23 × 10^-1^     -5.69 × 10^3^     -3.28 × 10^1^     4.73 × 10^-2^      Pabalan and Pitzer \[47\]
                                                              3.38              -2.32 × 10^4^     -1.39 × 10^2^     2.18 × 10^-1^      
                                                              -3.53 × 10^1^     2.17 × 10^6^      1.40 × 10^4^      -2.28 × 10^1^      
                                                              2.44 × 10^-2^     1.89 × 10^3^      1.08 × 10^1^      -1.56 × 10^-2^     
  MgOH^+^       Cl^-^                                         -1.00 × 10^-1^    0.00              0.00              0.00               Pitzer \[13\]
                                                              1.66              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
  Na^+^         Cl^-^                                         7.46 × 10^-2^     -4.71 × 10^2^     -1.85             1.66 × 10^-3^      Greenberg and Moller \[21\]
                                                              2.75 × 10^-1^     -5.21 × 10^2^     -2.88             4.71 × 10^-3^      
                                                              0.00              0.00              0.00              0.00               
                                                              1.54 × 10^-3^     4.81 × 10^1^      1.75 × 10^-1^     -1.56 × 10^-4^     
  Na^+^         C                                             3.62 × 10^-2^     1.11 × 10^3^      1.12 × 10^1^      -2.33 × 10^-2^     He and Morse \[26\]
                                                              1.5               4.41 × 10^3^      4.46 × 10^1^      -9.99 × 10^-2^     
                                                              0.00              0.00              0.00              0.00               
                                                              5.20 × 10^-3^     0.00              0.00              8.88 × 10^-16^     
  Na^+^         F^-^                                          2.15 × 10^-2^     0.00              0.00              0.00               Pitzer \[13\]
                                                              2.11 × 10^-1^     0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
  Na^+^         HC                                            2.80 × 10^-2^     6.83 × 10^2^      6.90              -1 .45 × 1 0^-2^   He and Morse \[26\]
                                                              4.40 × 10^-2^     1.13 × 10^3^      1.14 × 10^1^      -2.45 × 10^-2^     
                                                              0.00              0.00              0.00              0.00               
                                                              0.00              0.00              0.00              0.00               
  Na^+^         HS                                            7.34 × 10^-2^     5.26 × 10^1^      4.21 × 10^-^1     -8.21 × 10^-4^     Holmes and Mesmer \[24\]
                                                              3.00 × 10^-1^     4.70 × 10^3^      2.68 × 10^1^      -3.74 × 10^-2^     
                                                              0.00              0.00              0.00              0.00               
                                                              -4.62 × 10^-3^    -5.82 × 10^-11^   -2.27 × 10^-13^   6.82 × 10^-6^      
  Na^+^         N                                             3.57 × 10^-3^     -7.03 × 10^2^     -3.35             3.98 × 10^-3^      Archer \[33\]
                                                              2.32 × 10^-1^     -2.73 × 10^3^     -1.30 × 10^1^     2.07 × 10^-2^      
                                                              0.00              0.00              0.00              0.00               
                                                              -4.15 × 10^-5^    6.48 × 10^1^      3.18 × 10^-1^     -3.84 × 10^-4^     
  Na^+^         OH^-^                                         8.83 × 10^-2^     -1.20 × 10^3^     -6.11             7.43 × 10^-3^      He and Morse \[26\]
                                                              2.44 × 10^-1^     1.63 × 10^3^      9.48              -1.16 × 10^-2^     
                                                              0.00              0.00              0.00              0.00               
                                                              4.00 × 10^-3^     8.82 × 10^1^      4.07 × 10^-1^     -4.76 × 10^-4^     
  Na^+^         S                                             1.21 × 10^-2^     -2.19 × 10^3^     -1.01 × 10^1^     1.20 × 10^-2^      Greenberg and Moller \[21\]
                                                              1.12              -1.27 × 10^4^     -6.72 × 10^1^     8.91 × 10^-2^      
                                                              0.00              0.00              0.00              0.00               
                                                              6.57 × 10^-3^     3.39 × 10^2^      1.48              -1.64 × 10^-3^     
  Ca^++^        H^+^          ^*s*^*θ*~*MM\'*~                9.20 × 10^-2^     0.00              0.00              0.00               Pitzer \[13\]
  Ca^++^        K^+^          ^*s*^*θ*~*MM\'*~                1.16 × 10^-1^     9.31 × 10^-10^    0.00              1.42 × 10^-14^     Greenberg and Moller \[21\]
  Ca^++^        Na^+^         ^*s*^*θ*~*MM\'*~                5.00 × 10^-2^     1.86 × 10^-9^     7.28 × 10^-12^    0.00               Greenberg and Moller \[21\]
  Ca^++^        Mg^++^        ^*s*^*θ*~*MM\'*~                7.00 × 10^-3^     0.00              0.00              0.00               Pitzer \[13\]
  H^+^          K^+^          ^*s*^*θ*~*MM\'*~                5.00 × 10^-3^     0.00              0.00              0.00               Pitzer \[13\]
  H^+^          Mg^++^        ^*s*^*θ*~*MM\'*~                1.00 × 10^-^1     0.00              0.00              0.00               Pitzer \[13\]
  H^+^          Na^+^         ^*s*^*θ*~*MM\'*~                3.60 × 10^-2^     0.00              0.00              0.00               Pitzer \[13\]
  K^+^          Mg^++^        ^*s*^*θ*~*MM\'*~                0.00              0.00              0.00              0.00               Pabalan and Pitzer \[47\]
  K^+^          Na^+^         ^*s*^*θ*~*MM\'*~                -3.20 × 10^-3^    1.40 × 10^1^      9.09 × 10^-13^    -2.66 × 10^-15^    Greenberg and Moller \[21\]
  Mg^++^        Na^+^         ^*s*^*θ*~*MM\'*~                7.00 × 10^-2^     0.00              0.00              0.00               Pitzer \[13\]
  Cl^-^         C             ^*s*^*θ*~*XX\'*~                -2.00 × 10^-2^    0.00              0.00              0.00               Pitzer \[13\]
  Cl^-^         HC            ^*s*^*θ*~*XX\'*~                3.00 × 10^-2^     0.00              0.00              0.00               Pitzer \[13\]
  Cl^-^         HS            ^*s*^*θ*~*XX\'*~                -6.00 × 10^-3^    0.00              0.00              0.00               Pitzer \[13\]
  Cl^-^         N             ^*s*^*θ*~*XX\'*~                1.60 × 10^-2^     0.00              0.00              0.00               Pitzer \[13\]
  Cl^-^         OH^-^         ^*s*^*θ*~*XX\'*~                -5.00 × 10^-2^    0.00              0.00              0.00               Pabalan and Pitzer \[47\]
  Cl^-^         S             ^*s*^*θ*~*XX\'*~                7.03 × 10^-2^     1.83 × 10^2^      1.33              -2.33 × 10^-3^     Greenberg and Moller \[21\]
  C             HC            ^*s*^*θ*~*XX\'*~                -4.00 × 10^-2^    0.00              0.00              0.00               Pitzer \[13\]
  C             OH-           ^*s*^*θ*~*XX\'*~                1.00 × 10^-1^     0.00              0.00              0.00               Pitzer \[13\]
  C             S             ^*s*^*θ*~*XX\'*~                2.00 × 10^-2^     0.00              0.00              0.00               Pitzer \[13\]
  HC            S             ^*s*^*θ*~*XX\'*~                1.00 × 10^-2^     0.00              0.00              0.00               Pitzer \[13\]
  HS            S             ^*s*^*θ*~*XX\'*~                -1.17 × 10^-1^    -3.09 × 10^3^     -1.40 × 10^1^     1.54 × 10^-2^      Holmes and Mesmer \[24\]
  OH^-^         S             ^*s*^*θ*~*XX\'*~                -1.30 × 10^-2^    0.00              0.00              0.00               Pabalan and Pitzer \[47\]
  CO~2~(aq)     Ca^++^        ^*s*^*λ*~*NM*~                  1.98 × 10^-1^     -6.51 × 10^4^     -4.03 × 10^2^     -6.18 × 10^-1^     He and Morse \[26\]
  CO~2~(aq)     H^+^          ^*s*^*λ*~*NM*~                  0.00              0.00              0.00              0.00               He and Morse \[26\]
  CO~2~(aq)     K^+^          ^*s*^*λ*~*NM*~                  4.58 × 10^-2^     1.37 × 10^4^      8.36 × 10^1^      -1.28 × 10^-1^     He and Morse \[26\]
  CO~2~(aq)     Mg^++^        ^*s*^*λ*~*NM*~                  1.95 × 10^-1^     -3.96 × 10^3^     -3.17 × 10^1^     5.70 × 10^-2^      He and Morse \[26\]
  CO~2~(aq)     Na^+^         ^*s*^*λ*~*NM*~                  7.75 × 10^-2^     -1.92 × 10^4^     -1.17 × 10^2^     1.75 × 10^-1^      He and Morse \[26\]
  SiO~2~(aq)    H^+^          ^*s*^*λ*~*NM*~                  0.00              0.00              0.00              0.00               Felmy *et al*. \[32\]
  SiO~2~(aq)    Mg^++^        ^*s*^*λ*~*NM*~                  8.21 × 10^-2^     -3. 83 × 10^-9^   -2.33 × 10^-11^   -1.81 × 10^-3^     Felmy *et al*. \[32\]
  SiO~2~(aq)    Na^+^         ^*s*^*λ*~*NM*~                  -7.89 × 10^-2^    -1.44 × 10^-7^    -8.64 × 10^-10^   1.19 × 10^-4^      Felmy *et al*. \[32\]
  CO~2~(aq)     Cl^-^         ^*s*^*λ*~*NX*~                  2.02 × 10^-2^     5.08 × 10^3^      3.01 × 10^1^      -4.47 × 10^-2^     He and Morse \[26\]
  CO~2~(aq)     HS            ^*s*^*λ*~*NX*~                  -3.00 × 10^-3^    0.00              0.00              0.00               He and Morse \[26\]
  CO~2~(aq)     S             ^*s*^*λ*~*NX*~                  1.39 × 10^-1^     4.99 × 10^4^      .09 × 10^2^       -4.55 × 10^-1^     He and Morse \[26\]
  SiO~2~(aq)    Cl^-^         ^*s*^*λ*~*NX*~                  1.42 × 10^-1^     1.67 × 10^2^      1.51              -3.22 × 10^-3^     Felmy *et al*. \[32\]
  SiO~2~(aq)    N             ^*s*^*λ*~*NX*~                  1.34 × 10^-1^     6.30 × 10^3^      3.55 × 10^1^      -5.02 × 10^-2^     Felmy *et al*. \[32\]
  SiO~2~(aq)    S             ^*s*^*λ*~*NX*~                  7.76 × 10^-2^     1.56 × 10^-7^     9.33 × 10^-10^    5.53 × 10^-4^      Felmy *et al*. \[32\]
  Ca^++^K^+^    Cl^-^         *ψ*~*MM\'X*~                    -4.32 × 10^-2^    -2.71 × 10^1^     -3.64 × 10^-12^   3.55 × 10^-15^     Greenberg and Moller \[21\]
  Ca^++^Na^+^   Cl^-^         *ψ*~*MM\'X*~                    -3.00 × 10^-3^    0.00              1.14 × 10^-13^    -6.66 × 10^-16^    Greenberg and Moller \[21\]
  Ca^++^Na^+^   S             *ψ*~*MM\'X*~                    -1.20 × 10^-2^    0.00              4.55 × 10^-13^    -2.66 × 10^-15^    Greenberg and Moller \[21\]
  K^+^Mg^++^    Cl^-^         *ψ*~*MM\'X*~                    -2.20 × 10^-2^    -1.43 × 10^1^     0.00              0.00               Pabalan and Pitzer \[47\]
  K^+^Na^+^     Cl^-^         *ψ*~*MM\'X*~                    -3.69 × 10^-3^    -5.10             -3.41 × 10^-13^   6.66 × 10^-16^     Greenberg and Moller \[21\]
  K^+^Na^+^     S             *ψ*~*MM\'X*~                    7.32 × 10^-3^     -7.16 × 10^1^     -3.94 × 10^-1^    6.07 × 10^-4^      Greenberg and Moller \[21\]
  Mg^++^Na^+^   Cl^-^         *ψ*~*MM\'X*~                    -1.20 × 10^-2^    -9.51             0.00              0.00               Pabalan and Pitzen \[47\]
  Ca^++^Cl^-^   S             *ψ*~*MXX\'*~                    -1.80 × 10^-2^    -6.98 × 10^-10^   -2.73 × 10^-12^   1.78 × 10^-15^     Greenberg and Moller \[21\]
  K^+^Cl^-^     S             *ψ*~*MXX\'*~                    -1.62 × 10^-3^    3.76 × 10^1^      2.90 × 10^-12^    2.85 × 10^-4^      Greenberg and Moller \[21\]
  H^+^HS        S             *ψ*~*MXX\'*~                    0.00              0.00              0.00              0.00               Holmes and Mesmer \[24\]
  Mg^++^Cl^-^   S             *ψ*~*MXX\'*~                    -7.96 × 10^-3^    3.26 × 10^1^      0.00              0.00               Pabalan and Pitzer \[47\]
  Na^+^Cl^-^    OH^-^         *ψ*~*MXX\'*~                    -6.01 × 10^-3^    -9.93             0.00              0.00               Pabalan and Pitzer \[47\]
  Na^+^Cl^-^    S             *ψ*~*MXX\'*~                    -9.09 × 10^-3^    -7.86 × 10^1^     -5.52 × 10^-1^    9.46 × 10^-4^      Greenberg and Moller \[21\]
  Na^+^HS       S             *ψ*~*MXX\'*~                    1.44 × 10^-2^     2.58 × 10^2^      1.16              -1.26 × 10^-3^     Holmes and Mesmer \[24\]
  Na^+^OH^-^    S             *ψ*~*MXX\'*~                    -9.10 × 10^-3^    -1.17 × 10^1^     0.00              0.00               Pabalan and Pitzer \[47\]
  CO~2~(aq)     Ca^++^Cl^-^   *ξ*~*NMX*~                      -1.61 × 10^-2^    6.25 × 10^3^      3.90 × 10^1^      -6.04 × 10^-2^     He and Morse \[26\]
  CO~2~(aq)     H^+^Cl^-^     *ξ*~*NMX*~                      -4.65 × 10^-3^    -1.31 × 10^3^     -7.26             9.96 × 10^-3^      He and Morse \[26\]
  CO~2~(aq)     H^+^S         *ξ*~*NMX*~                      0.00              0.00              0.00              0.00               He and Morse \[26\]
  CO~2~(aq)     K^+^Cl^-^     *ξ*~*NMX*~                      -1.27 × 10^-2^    -9.33 × 10^3^     -5.65 × 10^1^     8.56 × 10^-2^      He and Morse \[26\]
  CO~2~(aq)     K^+^S         *ξ*~*NMX*~                      -4.10 × 10^-4^    -1.12 × 10^5^     -6.84 × 10^2^     1.04               He and Morse \[26\]
  CO~2~(aq)     Mg^++^Cl^-^   *ξ*~*NMX*~                      -1.53 × 10^-2^    -3.32 × 10^3^     -1.97 × 10^1^     2.94 × 10^-2^      He and Morse \[26\]
  CO~2~(aq)     Mg^++^S       *ξ*~*NMX*~                      -9.28 × 10^-2^    -6.09 × 10^4^     -3.64 × 10^2^     5.44 × 10^-1^      He and Morse \[26\]
  CO~2~(aq)     Na^+^Cl^-^    *ξ*~*NMX*~                      -5.50 × 10^-4^    -3.97 × 10^3^     -2.44 × 10^1^     3.73 × 10^-2^      He and Morse \[26\]
  CO~2~(aq)     Na^+^S        *ξ*~*NMX*~                      -3.73 × 10^-2^    -8.84 × 10^3^     -5.48 × 10^1^     8.49 × 10^-2^      He and Morse \[26\]
  SiO~2~(aq)    H^+^N         *ξ*~*NMX*~                      -3.30 × 10^-3^    0.00              0.00              0.00               Felmy *et al*. \[32\]
  SiO~2~(aq)    Mg^++^Cl^-^   *ξ*~*NMX*~                      -5.15 × 10^-2^    1.50 × 10^-8^     8.99 × 10^-11^    5.94 × 10^-4^      Felmy *et al*. \[32\]
  SiO~2~(aq)    Na^+^Cl^-^    *ξ*~*NMX*~                      -8.48 × 10^-15^   -1.84 × 10^-8^    -1.11 × 10^-10^   -2.00 × 10^-4^     Felmy *et al*. \[32\]

###### 

Concentration (molal) of a synthetic \"sulfate type\" Topopah Spring Tuff pore water as it evaporated at 95°C.

          Sample     Concentration factor   *p*H^a^   HC              Ca              Mg              Si              Na              K               Cl              NO~3~           SO~4~
  ------- ---------- ---------------------- --------- --------------- --------------- --------------- --------------- --------------- --------------- --------------- --------------- ---------------
  Leg l   FEC9-1     1.00                   na^*b*^   3.11 × 10^-4^   1.25 × 10^-3^   6.48 × 10^-4^   9.21 × 10^-4^   3.04 × 10^-3^   1.63 × 10^-4^   3.31 × 10^-3^   3.94 × 10^-4^   1.85 × 10^-3^
          FEC9-2     1.07                   na^*b*^   2.46 × 10^-4^   1.64 × 10^-3^   3.11 × 10^-4^   6.27 × 10^-4^   3.11 × 10^-3^   1.73 × 10^-4^   3.56 × 10^-3^   3.97 × 10^-4^   1.95 × 10^-3^
          FEC9-3     1.23                   na^*b*^   1.16 × 10^-4^   1.91 × 10^-3^   2.56 × 10^-4^   6.07 × 10^-4^   3.68 × 10^-3^   1.99 × 10^-4^   4.06 × 10^-3^   4.48 × 10^-4^   2.23 × 10^-3^
          FEC9-4     1.55                   7.94      nd^*c*^         2.31 × 10^-3^   2.95 × 10^-4^   6.97 × 10^-4^   4.42 × 10^-3^   2.46 × 10^-4^   5.09 × 10^-3^   5.27 × 10^-4^   2.80 × 10^-3^
          FEC9-5     2.08                   7.18      nd^*c*^         3.11 × 10^-3^   3.86 × 10^-4^   8.86 × 10^-4^   5.93 × 10^-3^   3.37 × 10^-4^   6.97 × 10^-3^   6.94 × 10^-4^   3.74 × 10^-3^
          FEC9-6     3.15                   7.57      nd^*c*^         4.36 × 10^-3^   5.58 × 10^-4^   1.15 × 10^-3^   9.23 × 10^-3^   5.15 × 10^-4^   1.00 × 10^-2^   1.15 × 10^-3^   4.54 ×10^-3^
                                                                                                                                                                                      
  Leg 2   FEC12-1    3.15                   8.6       1.59 × 10^-4^   4.69 × 10^-3^   6.05 × 10^-4^   1.38 × 10^-3^   1.10 × 10^-2^   5.82 × 10^-4^   9.80 × 10^-3^   1.32 × 10^-3^   4.64 × 10^-3^
          FEC12-2    3.44                   8.29      nd^*c*^         5.12 × 10^-3^   6.05 × 10^-4^   1.46 × 10^-3^   1.20 × 10^-2^   5.65 × 10^-4^   1.08 × 10^-2^   1.40 × 10^-3^   5.12 × 10^-3^
          FEC12-3    3.63                   7.91      nd^*c*^         5.43 × 10^-3^   6.17 × 10^-4^   1.51 × 10^-3^   1.28 × 10^-2^   6.70 × 10^-4^   1.14 × 10^-2^   1.45 × 10^-3^   5.42 × 10^-3^
          FEC12-4    4.23                   7.87      nd^*c*^         6.53 × 10^-3^   7.08 × 10^-4^   1.77 × 10^-3^   1.53 × 10^-2^   7.70 × 10^-4^   1.37 × 10^-2^   1.69 × 10^-3^   6.46 × 10^-3^
          FEC12-5    4.85                   7.59      nd^*c*^         7.49 × 10^-3^   8.14 × 10^-4^   2.05 × 10^-3^   1.74 × 10^-2^   8.91 × 10^-4^   1.62 × 10^-2^   1.92 × 10^-3^   7.65 × 10^-3^
          FEC12-6    5.93                   7.19      nd^*c*^         9.25 × 10^-3^   9.95 × 10^-4^   2.53 × 10^-3^   2.16 × 10^-2^   1.15 × 10^-3^   1.91 × 10^-2^   2.29 × 10^-3^   9.00 × 10^-3^
          FEC12-7    7.65                   7.61      nd^*c*^         1.17 × 10^-2^   1.26 × 10^-3^   3.18 × 10^-3^   2.72 × 10^-2^   1.42 × 10^-3^   2.49 × 10^-2^   2.92 × 10^-3^   1.17 × 10^-2^
          FEC12-8    11.79                  7.47      nd^*c*^         1.77 × 10^-2^   1.92 × 10^-3^   4.79 × 10^-3^   4.14 × 10^-2^   2.28 × 10^-3^   3.79 × 10^-2^   4.43 × 10^-3^   1.77 × 10^-2^
          FEC12-9    28.42                  7.32      nd^*c*^         2.25 × 10^-2^   4.52 × 10^-3^   7.91 × 10^-3^   1.05 × 10^-1^   5.66 × 10^-3^   9.80 × 10^-2^   1.12 × 10^-2^   2.35 × 10^-2^
                                                                                                                                                                                      
  Leg 3   FEC13-1    28.42                  7.69      1.92 × 10^-4^   2.17 × 10^-2^   1.06 × 10^-3^   2.24 × 10^-3^   1.11 × 10^-1^   6.54 × 10^-3^   1.05 × 10^-1^   1.18 × 10^-2^   2.24 × 10^-2^
          FEC13-2    38.63                  7.07      nd^*c*^         2.06 × 10^-2^   8.99 × 10^-4^   2.87 × 10^-3^   1.37 × 10^-1^   8.90 × 10^-3^   1.29 × 10^-1^   1.64 × 10^-2^   2.28 × 10^-2^
          FEC13-3    47.95                  6.99      nd^*c*^         2.04 × 10^-2^   9.80 × 10^-4^   3.38 × 10^-3^   1.57 × 10^-1^   1.27 × 10^-2^   1.52 × 10^-1^   2.00 × 10^-2^   2.42 × 10^-2^
          FEC13-4    64.87                  6.4       nd^*c*^         2.71 × 10^-2^   1.54 × 10^-3^   4.24 × 10^-3^   2.37 × 10^-1^   1.34 × 10^-2^   2.00 × 10^-1^   2.89 × 10^-2^   2.80 × 10^-2^
          FEC13-5    129.02                 6.21      nd^*c*^         2.91 × 10^-2^   2.62 × 10^-3^   6.41 × 10^-3^   4.55 × 10^-1^   2.55 × 10^-2^   3.66 × 10^-1^   4.59 × 10^-2^   3.24 × 10^-2^
          FEC13-6    187.64                 6.34      nd^*c*^         3.38 × 10^-2^   3.52 × 10^-3^   6.60 × 10^-3^   6.18 × 10^-1^   4.13 × 10^-2^   5.46 × 10^-1^   7.36 × 10^-2^   4.27 × 10^-2^
          FEC13-7    306.40                 na^*b*^   nd^*c*^         4.21 × 10^-2^   6.43 × 10^-3^   6.32 × 10^-3^   1.15 × 10^0^    6.99 × 10^-2^   9.93 × 10^-1^   1.29 × 10^-1^   5.85 × 10^-2^
                                                                                                                                                                                      
  Leg 4   FEC14-1    306.40                 8.076     1.95 × 10^-4^   4.06 × 10^-2^   1.88 × 10^-3^   3.09 × 10^-4^   1.12 × 10^0^    7.89 × 10^-2^   1.09 × 10^0^    1.32 × 10^-1^   5.80 × 10^-2^
          FEC14-2    361.93                 7.82      1.74 × 10^-4^   4.27 × 10^-2^   1.95 × 10^-3^   1.98 × 10^-4^   1.35 × 10^0^    9.22 × 10^-2^   1.28 × 10^0^    1.54 × 10^-1^   6.28 × 10^-2^
          FEC14-3    409.63                 7.476     1.67 × 10^-4^   4.38 × 10^-2^   2.10 × 10^-3^   1.53 × 10^-4^   1.52 × 10^0^    1.05 × 10^-1^   1.46 × 10^0^    1.76 × 10^-1^   6.67 × 10^-2^
          FEC14-4    465.88                 7.899     1.52 × 10^-4^   4.51 × 10^-2^   2.27 × 10^-3^   1.25 × 10^-4^   1.76 × 10^0^    1.21 × 10^-1^   1.68 × 10^0^    2.00 × 10^-1^   7.09 × 10^-2^
          FEC14-5    632.48                 7.772     1.39 × 10^-4^   4.19 × 10^-2^   2.81 × 10^-3^   1.05 × 10^-4^   2.37 × 10^0^    1.62 × 10^-1^   2.23 × 10^0^    2.73 × 10^-1^   7.24 × 10^-2^
          FEC14-6    820.18                 7.946     1.69 × 10^-4^   3.89 × 10^-2^   3.41 × 10^-3^   8.62 × 10^-5^   3.09 × 10^0^    2.13 × 10^-1^   2.90 × 10^0^    3.51 × 10^-1^   7.70 × 10^-2^
          FEC14-7    1227.81                7.376     1.53 × 10^-4^   2.82 × 10^-2^   4.83 × 10^-3^   6.16 × 10^-5^   4.66 × 10^0^    3.15 × 10^-1^   4.40 × 10^0^    5.26 × 10^-1^   8.41 × 10^-2^
          FEC14-8    1734.09                6.87      1.30 × 10^-4^   1.94 × 10^-2^   6.40 × 10^-3^   4.65 × 10^-5^   6.60 × 10^0^    4.40 × 10^-1^   6.20 × 10^0^    7.44 × 10^-1^   9.97 × 10^-2^
          FEC14-9    2601.94                6.971     2.94 × 10^-4^   1.23 × 10^-2^   8.48 × 10^-3^   3.12 × 10^-5^   7.20 × 10^0^    6.44 × 10^-1^   6.61 × 10^0^    1.09 × 10^0^    1.28 × 10^-1^
          FEC14-10   3389.64                7.322     2.66 × 10^-4^   1.03 × 10^-2^   9.02 × 10^-3^   1.89 × 10^-5^   7.42 × 10^0^    8.09 × 10^-1^   6.59 × 10^0^    1.42 × 10^0^    1.36 × 10^-1^

^a^Measured at room temperature.

^b^Not analyzed.

^c^Not detected.

Detection limits: F = 0.25 ppm, HCO~3~= 1 ppm.

###### 

Results of x-ray diffraction analysis of precipitates formed from complete evaporation.

                               Leg 1 (Exp. FEC 9)   Leg 2 (Exp. FEC 12)   Leg 3 (Exp. FEC 13)   Leg 4 (Exp. FEC 14)
  ---------------------------- -------------------- --------------------- --------------------- ---------------------
  Halite (NaCl)                X                    X                     X                     X
  Anhydrite (CaSO~4~)          X                    X                     X                     X
  Bassanite (2CaSO~4~·H~2~O)                        X                                           X
  Niter (KNO~3~)                                                                                X
  Nitratine (NaNO~3~)                                                                           X

###### 

Log *K*temperature grid.

                                                            log *K*temperature grid (°C)                                    
  ------------------------ -------------------------------- ------------------------------ ---------- ---------- ---------- -----------------------------
  Anhydrite                CaSO~4~                          -4.1258                        -4.3193    -4.7307    -5.3507    Greenberg and Moller \[21\]
  Aragonite                CaCO~3~                          2.3715                         1.9931     1.4762     0.9179     Wolery and Jarek \[12\]
  Arcanite                 K~2~SO~4~                        -2.1629                        -1.7916    -1.4973    -1.397     Greenberg and Moller \[21\]
  Artinite                 Mg~2~CO~3~(OH)~2~:3H~2~O         21.7371                        19.656     17.2642    15.1018    Wolery and Jarek \[12\]
  Brucite                  Mg(OH)~2~                        18.0898                        16.298     14.2674    12.4514    Wolery and Jarek \[12\]
  Burkeite                 Na~6~CO~3~(SO~4~)~2~             No_data                        9.5671     No_data    No_data    Harvie *et al*. \[48\]
  Calcite                  CaCO~3~                          2.2257                         1.8487     1.333      0.7743     Wolery and Jarek \[12\]
  Fluorite                 CaF~2~                           -10.3098                       -10.0371   -9.9067    -9.967     Wolery and Jarek \[12\]
  Glaserite                NaK~3~(SO~4~)~2~                 No_data                        -3.8027    No_data    No_data    Harvie *et al*. \[48\]
  Glauberite               Na~2~Ca(SO~4~)~2~                -4.7768                        -5.1827    -5.7677    -6.5632    Greenberg and Moller \[21\]
  Gypsum                   CaSO~4~:2H~2~O                   -4.5436                        -4.5199    -4.6536    -4.9448    Greenberg and Moller \[21\]
  Halite                   NaCl                             1.5012                         1.5857     1.6084     1.569      Greenberg and Moller \[21\]
  Hemihydrate              2CaSO~4~:H~2~O                   -3.4301                        -3.7773    -4.2698    -4.9372    Greenberg and Moller \[21\]
  Huntite                  CaMg~3~(CO~3~)~4~                12.9081                        10.301     7.0008     3.6895     Wolery and Jarek \[12\]
  Kieserite                MgSO~4~:H~2~O                    0.72                           -0.0239    -1.0716    -2.2058    Pabalan and Pitzer \[47\]
  Labile_Salt              Na~2~Ca~5~(SO~4~)~6~:3H~2~O      18.2726                        8.4195     6.7359     7.3639     Greenberg and Moller \[21\]
  Magnesite                MgCO~3~                          2.9734                         2.2936     1.4383     0.5875     Wolery and Jarek \[12\]
  Natrite                  Na~2~CO~3~                       11.6028                        10.984     10.2336    9.5237     Barin and Platzki \[49\]
  Niter                    KNO~3~                           -0.8418                        -0.2344    0.3772     0.8541     Barin and Platzki \[49\]
  Pentasalt                K~2~Ca~5~(SO~4~)~6~:H~2~O        -23.9466                       -24.5549   -26.5842   -29.0537   Greenberg and Moller \[21\]
  Quartz                   SiO~2~                           -4.1605                        -3.7501    -3.3553    -3.0132    Wolery and Jarek \[12\]
  Sellaite                 MgF~2~                           -9.2699                        -9.3939    -9.7091    -10.1577   Barin and Platzki \[49\]
  Sepiolite                Mg~4~Si~6~O~15~(OH)~2~:6H~2~O    32.3876                        30.4439    27.171     23.8968    Wolery and Jarek \[12\]
  Amorphous Silica         SiO~2~                           -3.124                         -2.7136    -2.4067    -2.1843    Wolery and Jarek \[12\]
  Nitratine (Soda Niter)   NaNO~3~                          0.7192                         1.1009     1.4544     1.6902     Barin and Platzki \[49\]
  Sylvite                  KC1                              0.6012                         0.9148     1.1871     1.3418     Greenberg and Moller \[21\]
  Syngenite                K~2~Ca(SO~4~)2:H~2~O             -7.5463                        -7.2618    -7.152     -7.2929    Greenberg and Moller \[21\]
  Talc                     Mg~3~Si~4~O~10~(OH)~2~           25.0261                        22.1646    18.6126    15.3725    Wolery and Jarek \[12\]
  Thenardite               Na~2~SO~4~                       -0.1329                        -0.2547    -0.4272    -0.6877    Greenberg and Moller \[21\]
  Foshagite                Ca~4~Si~3~O~9~(OH)~2~:0.5H~2~O   72.8341                        66.691     59.4664    52.8906    Wolery and Jarek \[12\]
  Gyrolite                 Ca~2~Si~3~O~7~(OH)~2~:1.5H~2~O   25.3869                        23.6797    21.5505    19.6199    Wolery and Jarek \[12\]
  Hillebrandite            Ca~2~SiO~3~(OH)~2~:0.17H~2~O     40.4592                        37.0757    33.142     29.565     Wolery and Jarek \[12\]
  Okenite                  CaSi~2~O~4~(OH)~2~:H~2~O         11.4698                        10.8948    10.147     9.5024     Wolery and Jarek \[12\]
  Plombierite              Ca~5~Si~6~H~11~O~22.5~           69.8789                        65.3842    59.9716    55.1602    Wolery and Jarek \[12\]
  Riversideite             Ca~5~H~2~(SiO~3~)~6~:2H~2~O      76.6335                        70.6194    63.3144    56.6105    Wolery and Jarek \[12\]
  Tobermorite              5CaO:6SiO~2~:5.5H~2~O            72.5004                        67.1518    60.6578    54.7498    Wolery and Jarek \[12\]
  Xonotlite                Ca~6~Si~6~O~17~(OH)~2~           102.0247                       93.3664    82.9922    73.4523    Wolery and Jarek \[12\]
