Isaac Scientific Publishing

Journal of Advances in Applied Mathematics

Mathematical Methods for Dissolved Phosphorus Concentrations

Download PDF (410.4 KB) PP. 143 - 148 Pub. Date: October 15, 2019

DOI: 10.22606/jaam.2019.44002


  • Kuiyuan Li*
    Department of Mathematics and Statistics, University of West Florida, Pensacola, FL 32514, USA
  • Josaphat Uvah
    Department of Mathematics and Statistics, University of West Florida, Pensacola, FL 32514, USA


Phosphorus is a nutrient contained in fertilizer that can run off from lawns and crops. Predicting the phosphorus concentrations is an important component in determining how phosphorus is cycled in the surrounding ecosystem. In this paper, we present the development of two mathematical methods. The first is a steady-state diffusion ordinary differential equation with given initial data. For this method, we estimate the unknown parameter values using least squares approximations for the data set without the boundary values. The second method is identical to the first except that the boundary values are imposed. We then distinguish our methods from existing approaches by deploying homotopy continuation to connect all time stages. With this approach, phosphorus concentrations can be estimated at all times and any depth. Using real-life data, we give an example to show that the methods are not only easy to use, but also provide estimates between any time stages and at any depth.


Steady-state diffusion, time stages, homotopy continuation.


[1] B. Berner, Early Diagenesis, A Theoretical Approach, Princeton University Press, New Jersey, (1980).

[2] B. Boudreau, Diagenetic Models and Their Implementation, New York, Springer-Verlag, (1997).

[3] B. Boudreau, A Method-of-Lines Code for Carbon and Nutrient Diagenesis in Aquatic Sediments, Computers and Geosciences (1996), 22, 479-496.

[4] Z. Ji, Hydrodynamics and Water Quality: Modeling Rivers, Lakes, and Estuaries, John Wiley & Sons, New Jersey, (2008).

[5] F. Meysman, et. al, Reactive Transport in Surface Sediments. I. Model Complexity and Software Quality, Computers and Geosciences (2003) 29, 91-300.

[6] K. Li, T. Y. Li, An Algorithm for Symmetric Tridiagonal Eigenproblems —— Divide and Conquer with Homotopy Continuation, SIAM J. Sci. Comput. Vol.14, No. 3, (1993), 735-751.

[7] K. Li, T. Y. Li, A Homotopy Algorithm for Symmetric Generalized Eigenproblem, Numerical Algorithms 4(1993), 167-195.

[8] N. Davila, K. Li, Mathematical Models for Dissolved Phosphorus Concentrations in the Gulf of Mexico, Tech. Report, Department of Mathematics and Statistics, University of West Florida, May 2008.