Data Fitting Techniques with Applications to Mineral Dissolution Kinetics

A common problem in chemical kinetics is the development of a rate law that describes the dependence of the reaction rate on the surrounding conditions such as concentrations of reacting species or temperature of the reacting media (see Chap. 1). The most direct approach to solving this problem is to measure reaction rates under systematically varied conditions and then to perform mathematical analyses on these data to determine the form of the rate law and to generate estimates of any unknown constants, or parameters, that make up the proposed rate law. Other chapters in this book provide information both for designing kinetics experiments and for selecting appropriate rate laws for a variety of geochemical reactions. In this chapter we describe the mathematical analyses – known collectively as curve fitting or regression analysis – that can be used to select a rate equation that matches a given data set, to generate estimates for any unknown parameters in the rate equation (e.g., rate constants or reaction orders), and to quantify the uncertainty associated with the estimated values for the parameters. As we traverse this entirely quantitative process we will attempt to describe the underlying, qualitative process of looking at kinetic data: plots to make, features of these plots to examine, and conceptual sketches to draw. In order to fit a curve to data one must first obtain the data. Over the past few decades, researchers have produced numerous high quality studies on the kinetics of mineral dissolution and the effects of experimental conditions such as pH, temperature, concentration of dissolved metal ions, mineral composition, etc. We have compiled this dissolution rate data for a number of rocks and minerals including apatite, basalt, biotite, hornblende, kaolinite, olivine, plagioclase, potassium feldspar, pyroxene, and quartz. We will use these data, which are available in Appendix I of

[1]  Lei Chou,et al.  Study of the weathering of albite at room temperature and pressure with a fluidized bed reactor , 1984 .

[2]  Oleg S. Pokrovsky,et al.  Kinetics and mechanism of forsterite dissolution at 25°C and pH from 1 to 12 , 2000 .

[3]  J. Ganor,et al.  The effect of pH and temperature on kaolinite dissolution rate under acidic conditions , 2002 .

[4]  M. Temkin The Kinetics of Some Industrial Heterogeneous Catalytic Reactions , 1980 .

[5]  C. nsnv,et al.  Dissolution rates of plagioclase at pH : 2 and 3 , 2007 .

[6]  Lei Chou,et al.  Steady-state kinetics and dissolution mechanisms of albite , 1985 .

[7]  P. Bennett,et al.  The dissolution of quartz in dilute aqueous solutions of organic acids at 25°C , 1988 .

[8]  Abraham Charnes,et al.  The Equivalence of Generalized Least Squares and Maximum Likelihood Estimates in the Exponential Family , 1976 .

[9]  Paul G Tratnyek,et al.  Central limit theorem for chemical kinetics in complex systems , 2005 .

[10]  W. E. Stewart,et al.  Discrimination and goodness of fit of multiresponse mechanistic models , 1998 .

[11]  M. Hochella,et al.  Sample preparation and aging effects on the dissolution rate and surface composition of diopside , 1989 .

[12]  G. V. Gibbs,et al.  Mechanisms of silica dissolution as inferred from the kinetic isotope effect , 1990 .

[13]  Sylvie Huet,et al.  Statistical tools for nonlinear regression : a practical guide with S-PLUS examples , 1997 .

[14]  D. Marquardt An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[15]  Don M. Miller,et al.  Reducing Transformation Bias in Curve Fitting , 1984 .

[16]  John E. Dennis,et al.  Algorithm 573: NL2SOL—An Adaptive Nonlinear Least-Squares Algorithm [E4] , 1981, TOMS.

[17]  Robert E. Davis,et al.  Statistics for the evaluation and comparison of models , 1985 .

[18]  J. Acker,et al.  The influence of pH on biotite dissolution and alteration kinetics at low temperature , 1992 .

[19]  J. Tester,et al.  Correlating quartz dissolution kinetics in pure water from 25 to 625°C , 1994 .

[20]  O. Pokrovsky,et al.  Experimental determination of the effect of dissolved CO2 on the dissolution kinetics of Mg and Ca silicates at 25 °C , 2005 .

[21]  K. Knauss,et al.  Diopside dissolution kinetics as a function of pH, CO2, temperature, and time , 1993 .

[22]  J. Donald Rimstidt,et al.  A high resolution study of forsterite dissolution rates , 2000 .

[23]  N. Draper,et al.  Applied Regression Analysis. , 1967 .

[24]  J. Ganor,et al.  The effect of pH on kaolinite dissolution rates and on activation energy , 1995 .

[25]  S. Brantley,et al.  Feldspar dissolution at 25°C and pH 3: Reaction stoichiometry and the effect of cations , 1995 .

[26]  Paul G Tratnyek,et al.  Applicability of single-site rate equations for reactions on inhomogeneous surfaces , 2004 .

[27]  P. Bennett,et al.  Quartz dissolution in organic-rich aqueous systems , 1991 .

[28]  W. Stumm Aquatic surface chemistry : chemical processes at the particle-water interface , 1987 .

[29]  Werner Stumm,et al.  Dissolution kinetics of kaolinite in acidic aqueous solutions at 25°C , 1992 .

[30]  J. Walther,et al.  A surface complex reaction model for the pH-dependence of corundum and kaolinite dissolution rates , 1988 .

[31]  L. Chou,et al.  Mechanism of kaolinite dissolution at room temperature and pressure Part II: kinetic study , 1999 .

[32]  C. Amrhein,et al.  Some factors affecting the dissolution kinetics of anorthite at 25°C , 1992 .

[33]  G. R. Holdren,et al.  Dissolution kinetics of experimentally shocked silicate minerals , 1989 .

[34]  K. Ishida Stochastic Approach to Nonequilibrium Thermodynamics of First-Order Chemical Reactions , 1966 .

[35]  C. Mörth,et al.  Dissolution of microcline and labradorite in a forest O horizon extract: the effect of naturally occurring organic acids , 2002 .

[36]  David S. Moore,et al.  Statistics: Concepts and Controversies , 1979 .

[37]  H. Hartley The Modified Gauss-Newton Method for the Fitting of Non-Linear Regression Functions by Least Squares , 1961 .

[38]  Rate control of weathering of silicate minerals at room temperature and pressure , 1988 .

[39]  H. Riedwyl Goodness of Fit , 1967 .

[40]  P. Frogner,et al.  Hornblende dissolution kinetics at 25°C , 1998 .

[41]  R. Ocampo-Pérez,et al.  Adsorption of Fluoride from Water Solution on Bone Char , 2007 .

[42]  S. Gíslason,et al.  Meteoric water-basalt interactions. I. A laboratory study , 1987 .

[43]  A. Putnis,et al.  The dissolution of apatite in the presence of aqueous metal cations at pH 2–7 , 1998 .

[44]  E. A. Boucher Stochastic approach to reaction and physico-chemical kinetics , 1974 .

[45]  D. Hudson.,et al.  Statistics for Physicists , 1972 .

[46]  F. Mackenzie,et al.  Experimental Study of Igneous and Sedimentary Apatite Dissolution: Control of pH, Distance from Equilibrium, and Temperature on Dissolution Rates , 2003 .

[47]  P. Green Iteratively reweighted least squares for maximum likelihood estimation , 1984 .

[48]  E. Busenberg,et al.  The dissolution kinetics of feldspars at 25°C and 1 atm CO2 partial pressure , 1976 .

[49]  J. Drever,et al.  Mineral dissolution rates in plot-scale field and laboratory experiments , 1993 .

[50]  G. R. Holdren,et al.  Reaction rate-surface area relationships during the early stages of weathering. II. Data on eight additional feldspars , 1987 .

[51]  G. Milliken,et al.  Estimable functions in the nonlinear model , 1976 .

[52]  C. Chatfield,et al.  Statistics for physicists , 1971 .

[53]  John E. Dennis,et al.  Numerical methods for unconstrained optimization and nonlinear equations , 1983, Prentice Hall series in computational mathematics.

[54]  Kenneth Levenberg A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .

[55]  T. Pačes The kinetics of base cation release due to chemical weathering , 1990 .

[56]  P. Brady,et al.  Controls on silicate dissolution rates in neutral and basic pH solutions at 25°C , 1989 .

[57]  J. Banfield,et al.  Effect of Microorganisms and Microbial Metabolites on Apatite Dissolution , 2002 .

[58]  P. Bloom,et al.  Rates and stoichiometry of hornblende dissolution over 115 days of laboratory weathering at pH 3.6 4.0 and 25 °C in 0.01 M lithium acetate , 1996 .

[59]  N. M. Rose Dissolution rates of prehnite, epidote, and albite , 1991 .

[60]  Narayanaswamy Balakrishnan,et al.  A General Purpose Approximate Goodness-of-Fit Test , 1995 .

[61]  H. Ahrens,et al.  Brownlee, K. A.: Statistical Theory and Methodology in Science and Engineering. John Wiley & Sons, New York 1965, 590 S., 70 Abb., Tafelanhang , 1968 .

[62]  G. R. Holdren,et al.  Dissolution rates of plagioclase at pH = 2 and 3 , 1991 .

[63]  D. Legates,et al.  Evaluating the use of “goodness‐of‐fit” Measures in hydrologic and hydroclimatic model validation , 1999 .

[64]  R. D'Agostino,et al.  Goodness-of-Fit-Techniques , 1987 .

[65]  D. Gillespie The chemical Langevin equation , 2000 .

[66]  W. House,et al.  Investigation of the pH dependence of the kinetics of quartz dissolution at 25 °C , 1992 .

[67]  John E. Dennis,et al.  An Adaptive Nonlinear Least-Squares Algorithm , 1977, TOMS.

[68]  S. Gíslason,et al.  The mechanism, rates and consequences of basaltic glass dissolution: I. An experimental study of the dissolution rates of basaltic glass as a function of aqueous Al, Si and oxalic acid concentration at 25°C and pH = 3 and 11 , 2001 .

[69]  P. Schweda,et al.  Kinetics of muscovite, phlogopite, and biotite dissolution and alteration at pH 1-4, room temperature , 1996 .

[70]  Paul G Tratnyek,et al.  Effects of carbonate species on the kinetics of dechlorination of 1,1,1-trichloroethane by zero-valent iron. , 2002, Environmental science & technology.

[71]  I. Langmuir THE ADSORPTION OF GASES ON PLANE SURFACES OF GLASS, MICA AND PLATINUM. , 1918 .

[72]  Patrick V. Brady,et al.  Kinetics of quartz dissolution at low temperatures , 1990 .

[73]  Eric H. Oelkers,et al.  An experimental study of forsterite dissolution rates as a function of temperature and aqueous Mg and Si concentrations , 2001 .

[74]  P. Dove,et al.  Dissolution rate of quartz in lead and sodium electrolyte solutions between 25 and 300°C: Effect of the nature of surface complexes and reaction affinity , 1994 .

[75]  Ralph G. Pearson,et al.  Kinetics and mechanism , 1961 .

[76]  K. Knauss,et al.  Dependence of albite dissolution kinetics on ph and time at 25°c and 70°c , 1986 .

[77]  N. Draper,et al.  Applied Regression Analysis , 1966 .

[78]  P. Bloom,et al.  The pH dependence of hornblende dissolution , 1999 .

[79]  R. Wogelius,et al.  Olivine dissolution kinetics at near-surface conditions , 1992 .

[80]  M. Hodson The influence of Fe-rich coatings on the dissolution of anorthite at pH 2.6 , 2003 .