A differential equation for specific catchment area

Analysis of the behavior of specific catchment area in a stream tube leads to a simple nonlinear differential equation describing the rate of change of specific catchment area along a flow path. The differential equation can be integrated numerically along a flow path to calculate specific catchment area at any point on a digital elevation model without requiring the usual estimates of catchment area and width. The method is more computationally intensive than most grid-based methods for calculating specific catchment area, so its main application is as a reference against which conventional methods can be tested. This is the first method that provides a benchmark for more approximate methods in complex terrain with both convergent and divergent areas, not just on simple surfaces for which analytical solutions are known. Preliminary evaluation of the D8, M8, digital elevation model networks (DEMON), and D-infinity methods indicate that the D-infinity method is the best of those methods for estimating specific catchment area, but all methods overestimate in divergent terrain.