Abstract. Estimation of erosion rate is an important component of landscape
evolution studies, particularly in settings where transience or spatial
variability in uplift or erosion generates diverse landform morphologies.
While bedrock rivers are often used to constrain the timing and magnitude of changes in baselevel lowering, hilltop curvature (or convexity), CHT, provides an additional opportunity to map variations in erosion rate given that average slope angle becomes insensitive to erosion rate owing to threshold slope processes. CHT measurement techniques applied in prior studies (e.g., polynomial functions), however, tend to be computationally
expensive when they rely on high-resolution topographic data such as lidar,
limiting the spatial extent of hillslope geomorphic studies to small study
regions. Alternative techniques such as spectral tools like continuous
wavelet transforms present an opportunity to rapidly document trends in
hilltop convexity across expansive areas. Here, we demonstrate how
continuous wavelet transforms (CWTs) can be used to calculate the Laplacian
of elevation, which we utilize to estimate erosion rate in three catchments
of the Oregon Coast Range that exhibit varying slope angle, slope length,
and hilltop convexity, implying differential erosion. We observe that
CHT values calculated with the CWT are similar to those obtained from
2D polynomial functions. Consistent with recent studies, we find that
erosion rates estimated with CHT from both CWTs and 2D polynomial
functions are consistent with erosion rates constrained with cosmogenic
radionuclides from stream sediments. Importantly, our CWT approach
calculates curvature at least 103 times more quickly than 2D
polynomials. This efficiency advantage of the CWT increases with domain
size. As such, continuous wavelet transforms provide a compelling approach
to rapidly quantify regional variations in erosion rate as well as
lithology, structure, and hillslope sediment transport processes, which are
encoded in hillslope morphology. Finally, we test the accuracy of CWT and 2D
polynomial techniques by constructing a series of synthetic hillslopes
generated by a theoretical nonlinear transport model that exhibit a range of
erosion rates and topographic noise characteristics. Notably, we find that
neither CWTs nor 2D polynomials reproduce the theoretically prescribed
CHT value for hillslopes experiencing moderate to fast erosion rates,
even when no topographic noise is added. Rather, CHT is systematically
underestimated, producing a power law relationship between erosion rate and
CHT that can be attributed to the increasing prominence of planar
hillslopes that narrow the zone of hilltop convexity as erosion rate
increases. As such, we recommend careful consideration of measurement length
scale when applying CHT to estimate erosion rate in moderate to
fast-eroding landscapes, where curvature measurement techniques may be prone to systematic underestimation.