Third-order geomorphometric variables (derivatives): definition, computation and utilization of changes of curvatures

Third-order geomorphometric variables (based on third derivatives of the altitudinal field) have been neglected in geomorphometry, but their application to the delimitation of surface objects will lead to their increasing significance in future. New techniques of computation, presented and evaluated here, facilitate their use. This paper summarizes recent knowledge concerning definition, computation and geomorphologic interpretation of these variables. Formulae defining various third-order variables are unified based on the physical definition of slope gradient. Methods for their computation are compared from the point of view of method error and error generated by digital elevation model (DEM) inaccuracy. For exact mathematical test surfaces, the most natural and simple variant of the method of central differences (CD2) shows a method error 2-3 times smaller than the other methods used recently in geomorphometry. However, success in coping with DEM inaccuracy depends (for a given grid mesh) on the number and weighting of points from which the derivative is computed. This was tested for surfaces with varying degrees of random error. Here least squares-based methods are the most effective for mixed derivatives (especially for finer grids and less accurate DEMs), while a variant of the CD method, that repeats numerical evaluation of first derivatives (CD1), is the most successful for derivatives in cardinal directions. The CD2 method is generally the most successful for coarser grids where the method error is dominant. Utilization of third-order variables is documented from examples of terrain feature (ridge, valley and edge) extraction and from a first statistical test of the hypothesis that real segments of the land surface have a tendency to a constant value of some morphometric variable. For detection of (sharp) ridges and valleys, it is shown that the rate of change of tangential curvature is inadequate: rate of change of normal curvature is also required. A basic confirmation of the constant-value tendency is provided.