About Us | Help Videos | Contact Us | Subscriptions

Soil Science Society of America Journal Abstract - DIVISION S-2—NOTES

Evaluation of numerical techniques applied to soil solution speciation including cation exchange


This article in SSSAJ

  1. Vol. 66 No. 2, p. 474-478
    Received: May 21, 2001

    * Corresponding author(s): pvaughan@ussl.ars.usda.gov
Request Permissions

  1. Peter J. Vaughan *
  1. Salinity Laboratory USDA-ARS, 450 W. Big Springs Rd., Riverside, CA 92507


Most existing models of soil solution speciation utilize either Newton–Raphson or Picard iteration to obtain a numerical solution to the nonlinear set of algebraic equations expressing mole balance and charge balance of free ions for major species as well as cation exchange. A computer program was written to test speed and accuracy of these and other methods. Picard iteration was fastest but produced a mean relative error (MRE) for mole and charge balance of 0.03 with no further convergence after a few iterations. A tensor (quadratic) method and two simplex methods converged to the correct result. The tensor method was preferable because the final result was more accurate; the rate of convergence was faster by 10 to 100 times, and convergence occurred for all compositions tested up to an ionic strength of 0.25. These results point out the value of testing various algorithms prior to implementation of speciation code applied to soils.

  Please view the pdf by using the Full Text (PDF) link under 'View' to the left.

Copyright © 2002. Soil Science SocietyPublished in Soil Sci. Soc. Am. J.66:474–478.