<p>A geochemical model, like phreeqc calculates the saturation index (SI) for minerals based on the composition of the solution, temperature, and possibly pressure. These are used to calculate the activities of the ions, based on the activity coefficients. These are then compared to the activities that would be present at saturation as defined by the solubility product (Ksp). The saturation index is &lt; 0 if the solution is undersaturated with respect the mineral, = 0 at saturation, and > 0 if supersaturated.</p>

<p>So when you say the model shows you are saturated with respect to gypsum, it indicates that by definition you are in equilibrium with solid phase gypsum. It does not indicate if you got there by dissolving gypsum or by precipitating gypsum. It is highly unlikely (but theoretically possible) that your SI is exactly 0 without one of those things happening.</p>

CC BY-SA 3.0
2017-11-17T06:13:22.100Z
Thank you so much for your answer. It is very helpful. When you say "It does not indicate if you got there by dissolving gypsum or by precipitating gypsum" don't you mean both are happening? Isn't that what maintains equilibrium (SI = 0) with the solid phase? From my understanding, chemical equilibrium with a mineral means both reactions (i.e. dissolution and precipitation) should be happening reversibly, so there is zero net precipitation or net dissolution.
chrisuw92
CC BY-SA 3.0
2017-11-18T21:25:40.077Z
I was talking about how the solution approaches equilibrium with the solid. You could have a solid containing gypsum and an undersaturated solution so that the gypsum dissolves until saturation is reached. Or you could, for example, have a solution containing calcium and not much sulfate in contact with sulfides that are oxidizing to provide sulfate causing gypsum to precipitate.
haresfur
CC BY-SA 3.0
2017-11-16T08:08:22.510Z
2017-11-16T08:08:22.510Z
haresfur
<p>I did some geochemical modelling and found groundwater samples I have are saturated with respect to the mineral gypsum. I need some help with interpretation, more information is given below. </p>

<p>The groundwater comes from an engineered aquifer that is comprised of sand. Prior to construction, the sand was in contact with water containing high sulfate and calcium concentrations (dissolved gypsum). There was also some pyrite oxidation occurring which undoubtedly lead to increased sulfate concentrations. </p>

<p>Does my geochemical modelling work give evidence that the groundwater dissolved $Ca$ and $SO_4$ concentrations are in equilibrium with solid phase gypsum? meaning the mineral is present even though it was completely dissolved initially?</p>

<p>thanks</p>

CC BY-SA 3.0
2017-11-15T11:05:34.810Z
What type of model do you use? The samples you have are real-world samples -- or from your model?
daniel.heydebreck
CC BY-SA 3.0
2017-11-15T11:07:40.293Z
they are real groundwater samples i collected. i used default database in phreeqc.
chrisuw92
<geochemistry><mineralogy><groundwater>
Gypsum saturation index
chrisuw92
# Gypsum saturation index

(Asked by: username_0 on 2017-11-15)

I did some geochemical modelling and found groundwater samples I have are saturated with respect to the mineral gypsum. I need some help with interpretation, more information is given below.

The groundwater comes from an engineered aquifer that is comprised of sand. Prior to construction, the sand was in contact with water containing high sulfate and calcium concentrations (dissolved gypsum). There was also some pyrite oxidation occurring which undoubtedly lead to increased sulfate concentrations.

Does my geochemical modelling work give evidence that the groundwater dissolved $Ca$ and $SO\_4$ concentrations are in equilibrium with solid phase gypsum? meaning the mineral is present even though it was completely dissolved initially?

thanks

**Tags:** geochemistry, mineralogy, groundwater

> **Comments:**
> **username_1**: What type of model do you use? The samples you have are real-world samples -- or from your model?  
> **username_0**: they are real groundwater samples i collected. i used default database in phreeqc.  

---

## 1 answers

### Answer by username_2 on 2017-11-16. Score: 1 | Accepted answer

A geochemical model, like phreeqc calculates the saturation index (SI) for minerals based on the composition of the solution, temperature, and possibly pressure. These are used to calculate the activities of the ions, based on the activity coefficients. These are then compared to the activities that would be present at saturation as defined by the solubility product (Ksp). The saturation index is < 0 if the solution is undersaturated with respect the mineral, = 0 at saturation, and > 0 if supersaturated.

So when you say the model shows you are saturated with respect to gypsum, it indicates that by definition you are in equilibrium with solid phase gypsum. It does not indicate if you got there by dissolving gypsum or by precipitating gypsum. It is highly unlikely (but theoretically possible) that your SI is exactly 0 without one of those things happening.

> **Comments:**
> **username_0**: Thank you so much for your answer. It is very helpful. When you say "It does not indicate if you got there by dissolving gypsum or by precipitating gypsum" don't you mean both are happening? Isn't that what maintains equilibrium (SI = 0) with the solid phase? From my understanding, chemical equilibrium with a mineral means both reactions (i.e. dissolution and precipitation) should be happening reversibly, so there is zero net precipitation or net dissolution.  
> **username_2**: I was talking about how the solution approaches equilibrium with the solid. You could have a solid containing gypsum and an undersaturated solution so that the gypsum dissolves until saturation is reached. Or you could, for example, have a solution containing calcium and not much sulfate in contact with sulfides that are oxidizing to provide sulfate causing gypsum to precipitate.  

---

