Thermodynamics, aqueous speciation and equilibrium

Why is a dissolved total concentration not equivalent to the chemical activity controlling mineral stability?

Synthetic field investigation. Fictional limestone spring samples span temperature, dissolved carbon and salinity; students compare simplified activity models with field constraints. The setting, sampling values and numerical observations are invented for instruction; they do not describe a monitored site or an operational geochemical assessment. Compare at least two hypotheses before drawing conclusions, and make every material balance dimensionally explicit.

By the end of this unit, students should be able to name the relevant reservoirs, express governing relations with units, calculate at least one interpretable quantity, check a competing explanation, and specify what extra observation would increase confidence. The worked calculations are tightly simplified and do not constitute a validated geochemical speciation model, a professional environmental investigation or a public-health threshold determination.

Mechanism and why it matters

The chemical potential includes standard-state and logarithmic activity terms; equilibrium occurs when the reaction Gibbs energy is zero. Identify the named chemical phase and its exchange boundaries. Distinguish measured species totals from inferred reactions; one pattern can have several mechanisms.

Sampling and measurement protocol

Write the reaction quotient with activities and compare it with a temperature-specific equilibrium constant. Record phase, sampling support, date, concentration or ratio units, calibration and detection limit. Define which quantities are observed and which depend on an assumed model.

Competing explanations and identifiability

An equilibrium constant copied from 25 °C need not describe a deep hot aquifer. Specify an independently testable observation capable of distinguishing this explanation from an alternative source, reaction or laboratory bias.

Worked numerical example — Reaction quotient ratio

Problem. A classroom reaction has Q=0.008 and K=0.02; calculate Q/K.

Mathematical setup: 0.008/0.02. After checking the listed mass, time and amount units, the calculated answer is 0.4 dimensionless. Write down the inputs separately from the chosen constant or reference so an independent reviewer can recompute the result. If the expression uses a logarithm or a ratio, confirm that its argument is dimensionless; if a molar amount is involved, convert grams to moles before applying reaction coefficients. Because the values are invented for this lesson, the arithmetic demonstrates the physical scale and not a measurement-derived regional conclusion.

Sensitivity check. Compare the sign and magnitude with a limiting case and perturb the main input by ten percent. Report whether the relationship is linear or nonlinear, and avoid overstating numerical precision.

Targeted investigation and practical interpretation

Compute the ratio Q/K and reason about the direction of the affinity. Create a six-row synthetic table with inputs, predicted outputs and one replicate. Compare a rival process and document a falsifiable difference.

Advanced reasoning: Thermodynamic scope

State the standard state and whether an input is an activity, molality or analytical total. Constants vary with temperature; ionic association and complexation change free-ion proportions. The simple classroom formula does not perform full mass-action speciation, so a physically credible field interpretation should demonstrate charge closure and consistent carbonate or redox boundaries.

Knowledge check and research questions

1. Recalculate the numerical example using consistent units.
2. Which competing physical process could reproduce the observation?
3. Which new measurement would discriminate the mechanisms?
4. Which assumed parameter creates the greatest model sensitivity?

Instructor benchmark. 0.4 dimensionless, conditional on the stated model and units.

Mechanism and why it matters

Ions interact electrostatically; a chemical analysis reports analytical totals while thermodynamic expressions require species activities. Identify the named chemical phase and its exchange boundaries. Distinguish measured species totals from inferred reactions; one pattern can have several mechanisms.

Sampling and measurement protocol

Estimate ionic strength with ion charges, state molality conventions and label any ideal-solution approximation. Record phase, sampling support, date, concentration or ratio units, calibration and detection limit. Define which quantities are observed and which depend on an assumed model.

Competing explanations and identifiability

Highly saline brines can invalidate a dilute-solution activity approximation even if every analytical ion total is precise. Specify an independently testable observation capable of distinguishing this explanation from an alternative source, reaction or laboratory bias.

Worked numerical example — Ionic strength

Problem. For idealized 0.015 mol/kg NaCl dissociated fully, I = 1/2 sum(mz²).

Mathematical setup: 0.5*(0.015+0.015). After checking the listed mass, time and amount units, the calculated answer is 0.015 mol/kg. Write down the inputs separately from the chosen constant or reference so an independent reviewer can recompute the result. If the expression uses a logarithm or a ratio, confirm that its argument is dimensionless; if a molar amount is involved, convert grams to moles before applying reaction coefficients. Because the values are invented for this lesson, the arithmetic demonstrates the physical scale and not a measurement-derived regional conclusion.

Sensitivity check. Compare the sign and magnitude with a limiting case and perturb the main input by ten percent. Report whether the relationship is linear or nonlinear, and avoid overstating numerical precision.

Targeted investigation and practical interpretation

Calculate ionic strength for a simple electrolyte mixture. Create a six-row synthetic table with inputs, predicted outputs and one replicate. Compare a rival process and document a falsifiable difference.

Advanced reasoning: Thermodynamic scope

State the standard state and whether an input is an activity, molality or analytical total. Constants vary with temperature; ionic association and complexation change free-ion proportions. The simple classroom formula does not perform full mass-action speciation, so a physically credible field interpretation should demonstrate charge closure and consistent carbonate or redox boundaries.

Knowledge check and research questions

1. Recalculate the numerical example using consistent units.
2. Which competing physical process could reproduce the observation?
3. Which new measurement would discriminate the mechanisms?
4. Which assumed parameter creates the greatest model sensitivity?

Instructor benchmark. 0.015 mol/kg, conditional on the stated model and units.

Mechanism and why it matters

Proton activity controls protonation distributions; total inorganic carbon can be partitioned among aqueous CO2, bicarbonate and carbonate. Identify the named chemical phase and its exchange boundaries. Distinguish measured species totals from inferred reactions; one pattern can have several mechanisms.

Sampling and measurement protocol

Track pH, total inorganic carbon, alkalinity and temperature as separate observations. Record phase, sampling support, date, concentration or ratio units, calibration and detection limit. Define which quantities are observed and which depend on an assumed model.

Competing explanations and identifiability

A laboratory pH shift from CO2 exchange may masquerade as a geological water–rock reaction. Specify an independently testable observation capable of distinguishing this explanation from an alternative source, reaction or laboratory bias.

Worked numerical example — pH from activity

Problem. If hydrogen-ion activity is 10^-6.6, calculate pH.

Mathematical setup: -math.log10(10(-6.6)). After checking the listed mass, time and amount units, the calculated answer is 6.6 pH units. Write down the inputs separately from the chosen constant or reference so an independent reviewer can recompute the result. If the expression uses a logarithm or a ratio, confirm that its argument is dimensionless; if a molar amount is involved, convert grams to moles before applying reaction coefficients. Because the values are invented for this lesson, the arithmetic demonstrates the physical scale and not a measurement-derived regional conclusion.

Sensitivity check. Compare the sign and magnitude with a limiting case and perturb the main input by ten percent. Report whether the relationship is linear or nonlinear, and avoid overstating numerical precision.

Targeted investigation and practical interpretation

Compute pH from specified hydrogen-ion activity. Create a six-row synthetic table with inputs, predicted outputs and one replicate. Compare a rival process and document a falsifiable difference.

Advanced reasoning: Thermodynamic scope

State the standard state and whether an input is an activity, molality or analytical total. Constants vary with temperature; ionic association and complexation change free-ion proportions. The simple classroom formula does not perform full mass-action speciation, so a physically credible field interpretation should demonstrate charge closure and consistent carbonate or redox boundaries.

Knowledge check and research questions

1. Recalculate the numerical example using consistent units.
2. Which competing physical process could reproduce the observation?
3. Which new measurement would discriminate the mechanisms?
4. Which assumed parameter creates the greatest model sensitivity?

Instructor benchmark. 6.6 pH units, conditional on the stated model and units.

Mechanism and why it matters

A saturation index is a logarithmic activity-product ratio relative to a mineral equilibrium constant; sign alone is not kinetic proof. Identify the named chemical phase and its exchange boundaries. Distinguish measured species totals from inferred reactions; one pattern can have several mechanisms.

Sampling and measurement protocol

Estimate an ion activity product from measured species and compare SI across plausible uncertainty bands. Record phase, sampling support, date, concentration or ratio units, calibration and detection limit. Define which quantities are observed and which depend on an assumed model.

Competing explanations and identifiability

A supersaturated solution may persist because nucleation barriers, inhibitors and reaction times limit precipitation. Specify an independently testable observation capable of distinguishing this explanation from an alternative source, reaction or laboratory bias.

Worked numerical example — Saturation index

Problem. If IAP = 4.0e-9 and K = 1.0e-8, calculate log10(IAP/K).

Mathematical setup: math.log10(4.0e-9/1.0e-8). After checking the listed mass, time and amount units, the calculated answer is -0.39794 SI units. Write down the inputs separately from the chosen constant or reference so an independent reviewer can recompute the result. If the expression uses a logarithm or a ratio, confirm that its argument is dimensionless; if a molar amount is involved, convert grams to moles before applying reaction coefficients. Because the values are invented for this lesson, the arithmetic demonstrates the physical scale and not a measurement-derived regional conclusion.

Sensitivity check. Compare the sign and magnitude with a limiting case and perturb the main input by ten percent. Report whether the relationship is linear or nonlinear, and avoid overstating numerical precision.

Targeted investigation and practical interpretation

Calculate a saturation index and test sensitivity. Create a six-row synthetic table with inputs, predicted outputs and one replicate. Compare a rival process and document a falsifiable difference.

Advanced reasoning: Thermodynamic scope

State the standard state and whether an input is an activity, molality or analytical total. Constants vary with temperature; ionic association and complexation change free-ion proportions. The simple classroom formula does not perform full mass-action speciation, so a physically credible field interpretation should demonstrate charge closure and consistent carbonate or redox boundaries.

Knowledge check and research questions

1. Recalculate the numerical example using consistent units.
2. Which competing physical process could reproduce the observation?
3. Which new measurement would discriminate the mechanisms?
4. Which assumed parameter creates the greatest model sensitivity?

Instructor benchmark. -0.39794 SI units, conditional on the stated model and units.

Mechanism and why it matters

Oxidation and reduction couple electron transfer to element speciation; measured potential is not always a single equilibrium pe. Identify the named chemical phase and its exchange boundaries. Distinguish measured species totals from inferred reactions; one pattern can have several mechanisms.

Sampling and measurement protocol

Balance electrons and distinguish Eh electrodes from a fully constrained thermodynamic redox calculation. Record phase, sampling support, date, concentration or ratio units, calibration and detection limit. Define which quantities are observed and which depend on an assumed model.

Competing explanations and identifiability

Mixed electrode potentials or microbial activity can make a field Eh reading nonunique. Specify an independently testable observation capable of distinguishing this explanation from an alternative source, reaction or laboratory bias.

Worked numerical example — Idealized 25 C redox term

Problem. For a one-electron redox ratio [oxidized]/[reduced]=10, compute 0.05916log10(ratio).

Mathematical setup: 0.05916*math.log10(10). After checking the listed mass, time and amount units, the calculated answer is 0.05916 V. Write down the inputs separately from the chosen constant or reference so an independent reviewer can recompute the result. If the expression uses a logarithm or a ratio, confirm that its argument is dimensionless; if a molar amount is involved, convert grams to moles before applying reaction coefficients. Because the values are invented for this lesson, the arithmetic demonstrates the physical scale and not a measurement-derived regional conclusion.

Sensitivity check. Compare the sign and magnitude with a limiting case and perturb the main input by ten percent. Report whether the relationship is linear or nonlinear, and avoid overstating numerical precision.

Targeted investigation and practical interpretation

Evaluate a simplified redox ratio using a Nernst-style relation. Create a six-row synthetic table with inputs, predicted outputs and one replicate. Compare a rival process and document a falsifiable difference.

Advanced reasoning: Thermodynamic scope

State the standard state and whether an input is an activity, molality or analytical total. Constants vary with temperature; ionic association and complexation change free-ion proportions. The simple classroom formula does not perform full mass-action speciation, so a physically credible field interpretation should demonstrate charge closure and consistent carbonate or redox boundaries.

Knowledge check and research questions

1. Recalculate the numerical example using consistent units.
2. Which competing physical process could reproduce the observation?
3. Which new measurement would discriminate the mechanisms?
4. Which assumed parameter creates the greatest model sensitivity?

Instructor benchmark. 0.05916 V, conditional on the stated model and units.

('Carbonate spring chemistry', 'Prepare paired spring samples collected in sealed and aerated conditions. For each, invent pH, alkalinity and calcium with an explicit temperature. Compare a no-exchange hypothesis with atmospheric carbon dioxide equilibration. Recalculate an ion activity product under two assumptions for the activity coefficient and report the saturation result as a range. State that a positive index cannot alone measure precipitation rate.')

Design and controls

Before computing any numbers, list the intended mechanistic question, the sampling unit, plausible confounding influences and the criteria for rejecting a proposed explanation. Separate control treatments from replicates. Specify one independent check that was not used to choose model coefficients, since a model cannot be validated merely by predicting its own calibration inputs. Synthetic scenarios are safe for classroom inference but do not replace field permissions, contamination control or specialist laboratory work.

A complete synthetic data record

Use a header with the following fields: sample_id, phase, dry_mass_g or volume_L, chemical_species, measured_value, analytical_unit, uncertainty, model_assumption, replicate_id and observation_status. Every derived column should include a formula and an audit trail. Avoid blanks that could ambiguously mean zero; label missing or below-quantification values explicitly. The class should archive its original input table and a derived working table so updates to calculations cannot overwrite evidence.

Expected result and alternatives

Analyze the invented measurements with the simplest mass balance or equilibrium relationship appropriate to the question. Choose plausible parameter ranges before looking at the pattern. Next deliberately generate an alternative dataset under a different process that yields a similar headline gradient. Explain which additional variable, sampling depth, time series or reference material would distinguish them. This step addresses the practical reality that geological and chemical processes are often nonunique when seen through sparse observations.

Sensitivity and uncertainty

Build a one-at-a-time sensitivity table varying at least two inputs by 5%, 10% and 20%. Use a separate scenario for systematic calibration shift, because correlated bias does not vanish by averaging independent samples. Include sampling uncertainty when a large reservoir is represented by a few rock chips or bottles. Report the resulting range as a model-conditioned bound rather than an objectively measured confidence interval unless a statistical sampling model was justified.

Interpretation and reporting

The final deliverable is a short field-style memo with a figure, a unit-checked table, a reproducible calculation and two rival explanations. Distinguish measured variables from estimated parameters by typography or column label. Provide the domain where the simple model is defensible, including temperature, composition and phase assumptions. Avoid making health, engineering or management decisions from synthetic data; those require validated data and appropriate regulatory and professional guidance.

Further experimental design

Propose two follow-up observations and rank them by discrimination power, cost and likelihood of analytical bias. A new sample that merely repeats the same confounded measurement may reduce noise without reducing the ambiguity of the hypothesis. A strategically chosen cross-phase, isotope or seasonal measurement can be more informative. Set a stopping rule: what numerical discrepancy or missing diagnostic signature would be sufficient to reject the favored model?

Research notebook prompts

1. At checkpoint 1, identify which assumptions in investigation 2 are directly tested, which still require external constraints, and how one alternative reaction or transport pathway could change the conclusion without necessarily altering the measured primary tracer. Record the expected sign of the change and the observation that would distinguish it.
2. At checkpoint 2, identify which assumptions in investigation 2 are directly tested, which still require external constraints, and how one alternative reaction or transport pathway could change the conclusion without necessarily altering the measured primary tracer. Record the expected sign of the change and the observation that would distinguish it.
3. At checkpoint 3, identify which assumptions in investigation 2 are directly tested, which still require external constraints, and how one alternative reaction or transport pathway could change the conclusion without necessarily altering the measured primary tracer. Record the expected sign of the change and the observation that would distinguish it.
4. At checkpoint 4, identify which assumptions in investigation 2 are directly tested, which still require external constraints, and how one alternative reaction or transport pathway could change the conclusion without necessarily altering the measured primary tracer. Record the expected sign of the change and the observation that would distinguish it.
5. At checkpoint 5, identify which assumptions in investigation 2 are directly tested, which still require external constraints, and how one alternative reaction or transport pathway could change the conclusion without necessarily altering the measured primary tracer. Record the expected sign of the change and the observation that would distinguish it.
6. At checkpoint 6, identify which assumptions in investigation 2 are directly tested, which still require external constraints, and how one alternative reaction or transport pathway could change the conclusion without necessarily altering the measured primary tracer. Record the expected sign of the change and the observation that would distinguish it.
7. At checkpoint 7, identify which assumptions in investigation 2 are directly tested, which still require external constraints, and how one alternative reaction or transport pathway could change the conclusion without necessarily altering the measured primary tracer. Record the expected sign of the change and the observation that would distinguish it.
8. At checkpoint 8, identify which assumptions in investigation 2 are directly tested, which still require external constraints, and how one alternative reaction or transport pathway could change the conclusion without necessarily altering the measured primary tracer. Record the expected sign of the change and the observation that would distinguish it.

The correct output of a geochemical investigation is not just a fitted number. A defensible conclusion links an observation to a specific physical or chemical mechanism while making conditional assumptions explicit. In this unit, write a three-column comparison: observed or invented data, quantitative model result, and possible geological explanation. Any statement that relies on a presumed endmember, reference standard or equilibrium assumption belongs in the last two columns. For a stronger claim, propose observations that would separate mechanism from coincidence.

Independent short assessment

• Choose one of the five numerical cases and reconstruct it using a different sequence of unit conversions. Confirm that the dimensions still cancel.
• Choose the most important environmental, geological or laboratory confounder from the five lessons. Explain whether it creates random scatter, systematic bias or a genuine alternative Earth process.
• Choose two variables of different spatial scales and explain why plotting them on the same grid would not establish causation.
• Write a transparent uncertainty statement: data are fictional, reference assumptions are specified, sensitivity has been tested, and a professional application requires additional measurements.

Suggested group discussion

Compare a precise but unrepresentative measurement with a noisy but representative sampling design. The former may look more convincing in a graph yet produce a more biased inventory or flux. Discuss a case where adding more significant digits does not improve geological accuracy. Discuss whether a model with a near-perfect fit has proved its reaction pathway or merely fitted several free parameters. End by identifying one prediction not used in the original fit, and describe how a prospective measurement could potentially falsify the process model.