Projections of Landscape Evolution on a 10,000 Year Timescale With Assessment and Partitioning of Uncertainty Sources

Long‐term erosion can threaten infrastructure and buried waste, with consequences for management of natural systems. We develop erosion projections over 10 ky for a 5 km2 watershed in New York, USA. Because there is no single landscape evolution model appropriate for the study site, we assess uncertainty in projections associated with model structure by considering a set of alternative models, each with a slightly different governing equation. In addition to model structure uncertainty, we consider the following uncertainty sources: selection of a final model set; each model's parameter values estimated through calibration; simulation boundary conditions such as the future incision of downstream rivers and future climate; and initial conditions (e.g., site topography which may undergo near‐term anthropogenic modification). We use an analysis‐of‐variance approach to assess and partition uncertainty in projected erosion into the variance attributable to each source. Our results suggest one sixth of the watershed will experience erosion exceeding 5 m in the next 10 ky. Uncertainty in projected erosion increases with time, and the projection uncertainty attributable to each source manifests in a distinct spatial pattern. Model structure uncertainty is relatively low, which reflects our ability to constrain parameter values and reduce the model set through calibration to the recent geologic past. Beyond site‐specific findings, our work demonstrates what information prediction‐under‐uncertainty studies can provide about geomorphic systems. Our results represent the first application of a comprehensive multi‐model uncertainty analysis for long‐term erosion forecasting. Plain Language Summary Erosion of ground material is a hazard to buildings and other infrastructure, and can pose an environmental risk when it occurs in areas such as radioactive waste repositories and post‐industrial sites. We make projections for erosion over the next 10,000 years and assess uncertainty sources at a 5 km2 watershed in New York state. Natural systems, like the study site, are not as well understood as engineered systems. However, information from studies like this one can provide useful insights to guide management decisions. Parts of the watershed experienced up to 50 m of erosion over the past 13,000 years. The type of model used to simulate land surface evolution over thousands of years is a Landscape Evolution Model (LEM). We use a set of alternative LEMs identified by prior