A numerical scheme for solving a stream nutrient model (Q1113776)

Water quality models have often neglected to incorporate the feedback effects of nutrient uptake and release by riparian vegetation. When this increasingle important effect is incorporated, the numerical solution methods used for standard models are no longer applicable. Such a model comprises a coupled system of first-order nonlinear partial differential equations, the characteristics of which are different and (in one case) depend on the state variables. This paper provides a method of numerical solution which can be used by practitioners in this and other applications. An illustrative example is given.