My Account: Log In | Join | Renew
Search
Author
Title
Vol.
Issue
Year
1st Page

Abstract

 

This article in SSSAJ

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


    * Corresponding author(s): pvaughan@ussl.ars.usda.gov
 View
 Download
 Alerts
 Permissions
 Share

doi:10.2136/sssaj2002.4740

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

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

Abstract

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.