|
|
||||||||
Salinity Laboratory USDA-ARS, 450 W. Big Springs Rd., Riverside, CA 92507
* Corresponding author (pvaughan{at}ussl.ars.usda.gov)
Most existing models of soil solution speciation utilize either NewtonRaphson 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.
Abbreviations: CEC, cation-exchange capacity MRE, mean relative error RMS, root mean squares
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Journal of Natural Resources and Life Sciences Education |
Vadose Zone Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||