USGS

 

CALCULATION OF INDIVIDUAL ISOTOPE EQUILIBRIUM CONSTANTS FOR IMPLEMENTATION IN GEOCHEMICAL MODELS

U.S. Geological Survey
Water-Resources Investigations Report 02-4172

By Donald C. Thorstenson and David L. Parkhurst


The full report is available in pdf.  Links to the pdf ).

Abstract

Theory is derived from the work of Urey to calculate equilibrium constants commonly used in geochemical equilibrium and reaction-transport models for reactions of individual isotopic species. Urey showed that equilibrium constants of isotope exchange reactions for molecules that contain two or more atoms of the same element in equivalent positions are related to isotope fractionation factors by Equation relating fractionation factor to equilibrium constant, where is n the number of atoms exchanged. This relation is extended to include species containing multiple isotopes, for example and Carbon-13, oxygen-16, oxygen-18 carbon dioxide, and Hydrogen-1, hydrogen-2, oxygen-18 water to include the effects of nonideality. The equilibrium constants of the isotope exchange reactions provide a basis for calculating the individual isotope equilibrium constants for the geochemical modeling reactions. The temperature dependence of the individual isotope equilibrium constants can be calculated from the temperature dependence of the fractionation factors. Equilibrium constants are calculated for all species that can be formed from Hydrogen-1, hydrogen-2Carbon-12, carbon-13, carbon-14 Oxygen-16, oxygen-18 and selected species containing Tritium, in the molecules Carbon dioxide gas, carbon dioxide aqueous Bicarbonate ionCarbonate ion Calcite Water vapor Liquid water Hydronium ion Hydroxide ion and the Sodium ion, potassium ionCalcium ion, magnesium ion ion pairs with where the subscripts g, aq, l, and s refer to gas, aqueous, liquid, and solid, respectively. These equilibrium constants are used in the geochemical model PHREEQC to produce an equilibrium and reaction-transport model that includes these isotopic species. Methods are presented for calculation of the individual isotope equilibrium constants for the asymmetric bicarbonate ion. An example calculates the equilibrium of multiple isotopes among multiple species and phases.

CONTENTS

Abstract

Introduction

Isotope ratios in gaseous and aqueous CO2

Species and symmetry numbers

Equilbrium constants for intraspecies reactions among the CO2 species

Isotope ratios

Calculation of individual isotope equilibrium constants for the reaction CO2gas - CO2aqueous

CO18O

C18O2

13CO2

13CO18O

13C18O2

14CO2, 14CO18O, and 14C18O2

Symmetry assumptions in the isotope ratios

Individual isotope equilibrium constants for bicarbonate

Individual isotope equilibrium constants for ion pairs

Application to other isotopes

Nonideality

Temperature dependence

Fundamental individual isotope equilibrium constants

Comparison of methods

Example calculation using PHREEQC

Summary

References cited


 

The text and graphics are presented here in pdf format (print quality):

The full report is 912KB

 

If you have Adobe Acrobat Reader installed on your computer, you may view and/or print the PDF version of this report. If you do not have Acrobat Reader, you may download it here.

Users with visual disabilities can visit this site for conversion tools and information to help make PDF files accessible.




FirstGov button  Take Pride in America button