Geology topics

Katherine R. Barnhart

Publications and source records attributed to Katherine R. Barnhart.

At least 19 recordsLinked to original sources

USGS2026

Regional models for postfire debris-flow likelihood and rainfall thresholds across the western United States

The U.S. Geological Survey (USGS) uses an empirical model developed with logistic regression (the ‘M1’ model) to rapidly assess debris-flow likelihood and to identify quantitative rainfall thresholds for debris flows after wildfire in the western United States. The M1 model was calibrated to a debris-flow inventory from southern California (United States) and has been applied throughout the western United States. Limited spatial coverage in the calibration dataset has motivated evaluation of M1 model accuracy outside the calibration region (e.g., the Sierra Nevada or the eastern Cascade Range, United States). Previous test cases showed that M1 overpredicts debris-flow likelihood and underpredicts rainfall thresholds for some locations (e.g., Arizona, northern California, Colorado, New Mexico, United States). We sought to improve the regional applicability of a debris-flow likelihood model by expanding the debris-flow inventory used for calibration, testing multiple potential models and generating an updated model framework. The updated inventory includes 3788 observations from 67 burned areas paired with short duration rainfall ratios. The updated model framework consists of a modified model structure and sets of coefficients calibrated separately to the entire updated inventory and to subsets of the inventory that intersect three Environmental Protection Agency (EPA) Level 2 ecoregions (Mediterranean California, Upper Gila Mountains and Western Cordillera). Comparisons of predictions from the updated models with observed rainfall and debris-flow activity show that the updated models outperform the M1 model by ~15%–60% and improve the uniformity of predictive performance across the western United States. The updated models also reduce false positive rates relative to M1 and generate rainfall thresholds that are better aligned with relative differences in regional climatology and debris-flow activity.

Arizona, California, Colorado, Idaho, Montana, Nev↗

USGS2026

An improved empirical model for predicting postfire debris-flow volume in the western United States

Reliable estimates of debris-flow volume can be used to help predict the magnitude of debris-flow hazards following wildfire in the western United States. In this study, we compiled and used a database of 227 postfire debris-flow volumes that were collected across the western United States to develop a multiple linear regression model for predicting postfire debris-flow volume. We explored 36 predictor variables related to rainfall, terrain, and fire characteristics, and selected the model with the combination of variables that yielded the most accurate predictions of debris-flow volume. We evaluated model performance against the entire volume database, as well as against four subsets of volume data from southern California, the Intermountain West, the Southwest, and regions with limited volume data, such as northern California and Washington. We also compared model performance against 3 existing postfire debris-flow volume models that were developed for use in southern California, the Intermountain West, and the Southwest. We demonstrate that the new volume model performs as well as the regional models in the regions for which they were developed and outperforms existing models when applied to volumes from data-limited regions in the western United States. These results indicate that the debris-flow volume model introduced in this study can be used to improve postfire hazard assessments across the western United States, especially outside of southern California.

Arizona, California, Colorado, New Mexico, Utah, W↗

USGS2026

A 481 m-high landslide-tsunami in a cruise ship-frequented Alaska fjord

Early in the morning of 10 August 2025, a >64 × 10 6 –cubic meter landslide struck Tracy Arm fjord in Alaska. The landslide was preconditioned by glacial retreat caused by climate change. The resulting 481-meter runup megatsunami followed an initial 100-meter-high breaking wave traveling at >70 meters per second. The landslide was preceded by several days of microseismicity, which increased in rate and magnitude until ~1 hour before failure. The landslide produced globally observed long-period seismic waves equivalent in size to a moment magnitude 5.4 earthquake. A long-period (~66 second) global seismic signal, produced by a landslide-induced seiche trapped within the fjord, persisted for up to 36 hours, the second time a days-long seiche had thus been observed. With fjord regions increasingly visited by cruise ships, and climate change making similar events more likely, this unanticipated, near-miss event highlights the growing risk from landslides and tsunamis in coastal environments.

Alaska↗

USGS2026

Channel morphology and large wood control postfire debris-flow erosion and deposition

Runoff-generated debris flows are a known response to wildfire, and accurately predicting the volume of these debris flows is important for estimating the magnitude of downstream hazards. Prior data collection efforts have focused on debris-flow volume measurements at catchment outlets, but few studies have considered how erosion and deposition modulate the volume of sediment arriving at catchment outlets. This study takes advantage of a high-resolution dataset to document the factors that control the total debris-flow volume reaching the catchment outlet during a fatal postfire debris flow. Using pre- and post-event airborne lidar, satellite imagery and field mapping, we found that a postfire debris flow in the Black Hollow catchment in northern Colorado eroded 136,000 ± 30,000 m 3 and redeposited 27,000 ± 7,500 m 3 in the main channel. Most of the in-channel deposition (52% by volume) occurred where a confined channel reach transitioned to an unconfined channel reach downstream, allowing the flow to widen and deposit material. Wood jams played multiple roles in the debris-flow dynamics, both nucleating deposition (25% of the deposit volume was stored behind wood jams) and exacerbating erosion (50% of the total erosion occurred downstream from a wood dam break). The remaining deposition occurred due to spatial changes in channel slope as well as deposition observed at newly formed channel bars. Using these data in this study, we identified topographic and vegetation metrics that can be used (pre-event) to anticipate where deposition may occur in channels prior to a debris flow.

Colorado↗

USGS2026

Multi-site evaluation of a postfire debris-flow runout forecast method

Postfire debris flows pose a hazard to human life, property, and infrastructure when they travel from steep source areas to urbanized alluvial fans or other developed areas. Existing methods for rapid (<1 week) postfire debris-flow hazard assessment document the increase in the likelihood and size of debris flows as the magnitude of high-intensity rain necessary to initiate debris flows increases but do not indicate the extent of downstream debris-flow runout. Although many models for the simulation of debris-flow motion are available, there is no established approach for using these models to delineate locations susceptible to postfire debris-flow runout that (a) is feasible to use at the spatial scale of an entire fire; (b) is appropriate for runout onto unconfined areas; (c) reproduces observed relations between runout and rainfall intensity; and (d) characterizes inherent uncertainty in runout, even without spatiotemporally variable rainfall. We propose and evaluate a method for generating postfire debris-flow runout hazard maps that has all the above qualities. Selection of case studies prioritized events triggered by a range of rainfall intensities, locations within and outside of southern California, and observed runout onto unconfined topography. Qualitative and quantitative assessment of performance for four events indicate that simulation results broadly match observations albeit with some discrepancies at a scale larger than structure or land parcel level (approximately 20-m by 20-m). The method may be used to identify potentially hazardous areas immediately following a fire and to provide approximate runout forecasts when a storm is imminent.

Arizona, California↗

USGS2025

2024 Surprise Inlet landslides: Insights from a prototype landslide‐triggered tsunami monitoring system in Prince William Sound, Alaska

Alaska's coastal communities face growing landslide hazards owing to glacier retreat and extreme weather intensified by the warming climate, yet hazard monitoring remains challenging. As part of ongoing experimental monitoring in Prince William Sound, we detected three large landslides (0.5–2.3 M m 3 ) at Surprise Inlet on 20 September 2024, within the span of an hour. These events were identified in near real-time through seismic data and later confirmed using satellite imagery, tidal records, and infrasound. The landslides generated a modest tsunami, and a 4 cm wave was recorded by a tide gauge 18 km away, marking the first recorded landslide to reach water since monitoring began in this region in 2021. Here, we examine the detection and interpretation of these landslides using multiple data sources and modeling. We demonstrate the effectiveness of this regional seismic monitoring system and show how complementary instrumentation, where available, can enhance detection capabilities.

Alaska↗

USGS2025

A method to obtain remotely sensed grain size distributions from nonplanar granular deposits

Constraining the grain size distribution of granular deposits with complex surfaces is difficult with existing approaches. Field and laboratory techniques are time consuming and limited by the maximum grain size that laboratories can accommodate. In this study, we present a new method to identify the coarse fraction of the grain size distribution at a debris-flow fan deposit surveyed with terrestrial laser scanning (TLS) in Glenwood Canyon, Colorado, USA. This method is a novel grain segmentation algorithm developed for application to point cloud data of deposits with complex surfaces and angular grains ranging in size from centimeters to a meter. This approach combines an existing random forest machine learning method with a novel iterative clustering algorithm. We compared the grain size distribution from our algorithm with a Wolman pebble count conducted in the field, and found a root mean squared error of less than 2 cm from the 5th to 95th percentile of the grain size distribution of grains ranging from cobble to boulder sized (6.3–78 cm in our application). Finally, we compared our new algorithm with an existing open-source grain segregation algorithm, and our method outperformed the selected alternative when applied to the debris-flow deposit point cloud.

Colorado↗

USGS2025

Detection of landslide-generated tsunami by shipborne GNSS precise point positioning

Precise point positioning (PPP) of ships using Global Navigation Satellite System (GNSS) data reveals the precise movements of marine vessels. This method may quantify anomalies in sea surface height with implications for oceanographic monitoring, exploration, and tsunami warning. The GNSS PPP data from the R/V Sikuliaq , a research ship of the University of Alaska Fairbanks, were processed to detect a small local tsunami generated by the Lowell Point landslide, which occurred near Seward, Alaska, on 8 May 2022 (UTC). The GNSS receiver aboard the R/V Sikuliaq recorded the waves generated by the landslide, with a maximum wave amplitude of 6 cm and wave periods between 40 and 50 s. These results are consistent with simulations of the landslide event.

Alaska↗

USGS2025

Uncertainty reduction for subaerial landslide-tsunami hazards

Subaerial rock slopes may generate a tsunami by rapidly moving into the water. Large uncertainty in landslide characteristics propagates into large uncertainty in tsunami hazard, making hazard assessment more difficult for land and emergency managers. Once a potentially tsunamigenic landslide is identified, it may not be clear which landslide characteristics contribute most significantly to uncertainty in the tsunami hazard. Our aim is to document the relative worth of different landslide characteristics (e.g., size, material properties) for reducing uncertainty in landslide-tsunami hazard assessments. Isolating the relative importance of specific landslide characteristics may inform prioritization of data collection and improve efficiency in understanding hazard. To accomplish this, we generated a set of 288 landslide-tsunami simulations in which we systematically varied the size and material properties of possible failure extents at the Barry Arm landslide complex in northwestern Prince William Sound, Alaska, USA. We find that for landslides smaller than the receiving waterbody, the landslide volume has the strongest effect on resulting wave characteristics and thus the highest leverage on reducing uncertainty in tsunami hazard. In contrast, for landslides substantially larger than the waterbody, the duration of rapid movement of the landslide has the highest leverage. Based on our results, we propose a classification scheme for subaerial landslides based on the relative size of the landslide and waterbody. Additionally, our results support the generation of a tsunami height transfer function between existing tide gages and a nearby coastal city. These results have direct implications for the practice of operational early warning.

Alaska↗

USGS2025

Cascading land surface hazards as a nexus in the Earth system

Earth’s surface is sculpted by numerous processes that move sediment, ranging from gradual and benign to abrupt and catastrophic. Although infrequent, high-magnitude sediment mobilization events can be hazardous to people and infrastructure, leaving topographic imprints on the landscape and remarkable narratives in the historical record. Hazardous events such as fires, storms, and earthquakes accelerate erosion and sediment transport, increasing landscape sensitivity to subsequent perturbations, thus forming a cascading hazard. Although the redistribution of sediment across Earth’s landscape can result in higher risks to vulnerable populations, cascading processes are commonly unaccounted for in hazard assessments. Cascading hazards can occur almost immediately after triggering events, such as coseismic landslides, or over months, years, or even decades after an initial perturbation, such as debris flows after wildfires or flooding in channels alluviated by volcanic debris. Sediment cascades span Earth’s surface, from mountaintops to river valleys, where erosion, deposition, and aggradation can lead to a myriad of hazardous processes, including decreased river conveyance capacity, which increases the likelihood of downstream flooding. An improved understanding of the magnitude, frequency, and persistence of cascading hazards is critical given the rapid changes in the frequency and severity of storms, fires, sea-level change, and cryospheric melting, as well as the expansion of high-population-density urban footprints in regions susceptible to solid Earth hazards. Understanding the full consequences and underlying physics of Earth’s cascading land surface hazards can help minimize future human and economic losses.

California↗

USGS2024

Constraining mean landslide occurrence rates for non-temporal landslide inventories using high-resolution elevation data

Constraining landslide occurrence rates can help to generate landslide hazard models that predict the spatial and temporal occurrence of landslides. However, most landslide inventories do not include any temporal data due to the difficulties of dating landslide deposits. Here we introduce a method for estimating the mean landslide occurrence rate of deep-seated rotational and translational slides derived solely from high-resolution (≤3 m) elevation data and globally available estimates of the diffusion coefficient for sediment flux. The method applies a linear diffusion model to the roughest landslide deposits until they reach a representative non-landslide roughness distribution. This estimates the time for a landslide deposit to be unrecognizable in high-resolution digital elevation data, which we term the mean lifetime of the landslide. Using the mean lifetime and number of landslides within an area of interest, we can estimate the mean occurrence rate of landslides over that domain. We validate this approach using a comprehensive temporal inventory of landslides in western Oregon created using age-roughness curves that are calibrated with high-resolution elevation data and radiocarbon data. We find good agreement between our diffusion method and the existing age-roughness-derived estimates, producing mean lifetimes of 4500 and 5200 years (4% difference), respectively. Hazard maps produced using the two methodologies generally agree, with the maximum differences in landslide probability reaching 0.1. Due to the relative abundance of high-resolution elevation data compared with age-dated landslides, our method could help constrain landslide occurrence rates in areas previously considered unfeasible.

Oregon↗

USGS2024

Probabilistic assessment of postfire debris-flow inundation in response to forecast rainfall

Communities downstream of burned steep lands face increases in debris-flow hazards due to fire effects on soil and vegetation. Rapid postfire hazard assessments have traditionally focused on quantifying spatial variations in debris-flow likelihood and volume in response to design rainstorms. However, a methodology that provides estimates of debris-flow inundation downstream of burned areas based on forecast rainfall would provide decision-makers with information that directly addresses the potential for downstream impacts. We introduce a framework that integrates a 24 h lead-time ensemble precipitation forecast with debris-flow likelihood, volume, and runout models to produce probabilistic maps of debris-flow inundation. We applied this framework to simulate debris-flow inundation associated with the 9 January 2018 debris-flow event in Montecito, California, USA. When the observed debris-flow volumes were used to drive the probabilistic forecast model, analysis of the simulated inundation probabilities demonstrates that the model is both reliable and sharp. In the fully predictive model, however, in which debris-flow likelihood and volume were computed from the atmospheric model ensemble's predictions of peak 15 min rainfall intensity, I 15 , the model generally under-forecasted the inundation area. The observed peak I 15 lies in the upper tail of the atmospheric model ensemble spread; thus a large fraction of ensemble members forecast lower I 15 than observed. Using these I 15 values as input to the inundation model resulted in lower-than-observed flow volumes which translated into under-forecasting of the inundation area. Even so, approximately 94 % of the observed inundated area was forecast to have an inundation probability greater than 1 %, demonstrating that the observed extent of inundation was generally captured within the range of outcomes predicted by the model. Sensitivity analyses indicate that debris-flow volume and two parameters associated with debris-flow mobility exert significant influence on inundation predictions, but reducing uncertainty in postfire debris-flow volume predictions will have the largest impact on reducing inundation outcome uncertainty. This study represents a first step toward a near-real-time hazard assessment product that includes probabilistic estimates of debris-flow inundation and provides guidance for future improvements to this and similar model frameworks by identifying key sources of uncertainty.

California↗

USGS2024

Evaluating post-wildfire debris-flow rainfall thresholds and volume models at the 2020 Grizzly Creek Fire in Glenwood Canyon, Colorado, USA

As wildfire increases in the western United States, so do postfire debris-flow hazards. The U.S. Geological Survey (USGS) has developed two separate models to estimate (1) rainfall intensity thresholds for postfire debris-flow initiation and (2) debris-flow volumes. However, the information necessary to test the accuracy of these models is seldom available. Here, we studied how well these models performed over a 2-year period in the 2020 Grizzly Creek Fire burn perimeter in Glenwood Canyon, Colorado, USA, through the development of a debris-flow response inventory. The study area had the advantage of a network of 11 rain gauges for rainfall intensity measurements and repeat lidar data for volume estimates. Our observations showed that 89 % of observed debris flows in the first year postfire were triggered by rainfall rates higher than the fire-wide rainfall threshold produced by the current USGS operational model (M1). No debris flows were observed in the second year postfire, despite eight rainstorms with intensities higher than the modeled rainfall threshold. We found that the operational model for debris-flow initiation rainfall thresholds works well in this region during the first year but may be too conservative in year 2 due to vegetation recovery and sediment depletion. However, rainfall thresholds in the second year can be improved by using updated remote sensing imagery to recalculate the debris-flow initiation likelihood with the M1 model. The current volume model overestimates for this region by a median value of 4.4 times. However, the offset between the model estimates and observations is approximately linear, and the volumes from the Grizzly Creek debris flows had a similar magnitude to historic postfire debris flows in the region. Consequently, the current volume model could be adjusted with a regional correction factor.

Colorado↗

USGS2024

Catchment coevolution and the geomorphic origins of variable source area hydrology

Features of landscape morphology—including slope, curvature, and drainage dissection—are important controls on runoff generation in upland landscapes. Over long timescales, runoff plays an essential role in shaping these same features through surface erosion. This feedback between erosion and runoff generation suggests that modeling long-term landscape evolution together with dynamic runoff generation could provide insight into hydrological function. Here we examine the emergence of variable source area runoff generation in a new coupled hydro-geomorphic model that accounts for water balance partitioning between surface flow, subsurface flow, and evapotranspiration as landscapes evolve over millions of years. We derive a minimal set of dimensionless numbers that provide insight into how hydrologic and geomorphic parameters together affect landscapes. Across the parameter space we investigated, model results collapsed to a single inverse relationship between the dimensionless relief and the ratio of catchment quickflow to discharge. Furthermore, we found an inverse relationship between the Hillslope number, which describes topographic relief relative to aquifer thickness, and the proportion of the landscape that was variably saturated. While the model generally produces fluvial topography visually similar to simpler landscape evolution models, certain parameter combinations produce wide valley bottom wetlands and non-dendritic, trellis-like drainage networks, which may reflect real conditions in some landscapes where aquifer gradients become decoupled from topography. With these results, we demonstrate the power of hydro-geomorphic models for generating new insights into hydrological processes, and also suggest that subsurface hydrology may be integral for modeling aspects of long-term landscape evolution.

Water Resources Research↗

USGS2024

Evaluation of debris-flow building damage forecasts

Reliable forecasts of building damage due to debris flows may provide situational awareness and guide land and emergency management decisions. Application of debris-flow runout models to generate such forecasts requires combining hazard intensity predictions with fragility functions that link hazard intensity with building damage. In this study, we evaluated the performance of building damage forecasts for the 9 January 2018 Montecito postfire debris-flow runout event, in which over 500 buildings were damaged. We constructed forecasts using either peak debris-flow depth or momentum flux as the hazard intensity measure and applied each approach using three debris-flow runout models (RAMMS, FLO-2D, and D-Claw). Generated forecasts were based on averaging multiple simulations that sampled a range of debris-flow volume and mobility, reflecting typical sources and magnitude of pre-event uncertainty. We found that only forecasts made with momentum flux and the D-Claw model could correctly predict the observed number of damaged buildings and the spatial patterns of building damage. However, the best forecast only predicted 50 % of the observed damaged buildings correctly and had coherent spatial patterns of incorrectly predicted building damage (i.e., false positives and false negatives). These results indicate that forecasts made at the building level reliably reflect the spatial pattern of damage but do not support interpretation at the individual building level. We found the event size strongly influences the number of damaged buildings and the spatial pattern of debris-flow depth and velocity. Consequently, future research on the link between precipitation and the volume of sediment mobilized may have the greatest effect on reducing uncertainty in building damage forecasts. Finally, because we found that both depth and velocity are needed to predict building damage, comparing debris-flow models against spatially distributed observations of building damage is a more stringent test for model fidelity than comparison against the extent of debris-flow runout.

Natural Hazards and Earth System Sciences↗

USGS2024

Unlearning Racism in Geoscience (URGE): Summary of U.S. Geological Survey URGE pod deliverables

The U.S. Geological Survey (USGS) is in a unique position to be a leader in diversity, equity, inclusion, and accessibility in the Earth sciences. As one of the largest geoscience employers, the USGS wields significant community influence and has a responsibility to adopt and implement robust, unbiased policies so that the science it is charged to deliver is better connected to the diverse communities it serves. Meaningful and effective improvements in implementation of diversity, equity, inclusion, and accessibility principles made within the USGS will be seen across the geoscience community. Despite this opportunity, however, the community, including the USGS, has struggled to diversify the geoscience workforce, which does not reflect the Nation’s diversity. This disparity suggests that the implementation of past policies did not achieve their desired outcomes in sustainable ways. The persistent lack of diversity across the geosciences and the racial justice protests of 2020 motivated a group of geoscience scholars to launch a global, virtual initiative known as Unlearning Racism in Geoscience (URGE) in the winter of 2021. The 16-week facilitated URGE curriculum was designed to highlight existing literature, share expert opinions, and foster discussion within groups (pods) around the world to build action plans for affecting change within their home institutions ( https://urgeoscience.org/ ). Whereas most of the approximately 4,500 participants were university faculty and students or from professional societies and nonprofit organizations, 6 pods totaling over 120 participants from the USGS joined the program—the largest participation of any Federal science agency. Summarizing the recommendations of over 100 USGS employees across the USGS who participated in URGE, this Circular represents a grassroots plan for making the USGS workforce more diverse and inclusive. To identify how peer-reviewed best practices could be implemented at the USGS, participants read existing literature, reviewed USGS policies and procedures, explored new and existing demographic data, and listened to and discussed primary interviews with scholars, all facilitated by the URGE program. This Circular leverages the efforts and momentum of the USGS URGE participants to catalyze and support positive, systemic change throughout the organization in a holistic way that represents an achievable departure from past policies that fell short of success. These recommendations are starting points. Reflecting the overall lack of diversity in the USGS workforce, the USGS URGE participants were predominantly White, meaning recommendations in this document can and should continue to evolve to include the lived perspectives and input from the broader USGS community. Finally, this summary document only highlights some of the key findings of these working groups. Additional details that may be more suitable to specific USGS Centers and programs can be found in the primary documents developed by each USGS URGE pod; a list of all registered pods can be found at https://urgeoscience.org/pods/ .

Circular↗

USGS2024

Satellite interferometry landslide detection and preliminary tsunamigenic plausibility assessment in Prince William Sound, southcentral Alaska

Regional mapping of actively deforming landslides, including measurements of landslide velocity, is integral for hazard assessments in paraglacial environments. These inventories are also critical for describing the potential impacts that the warming effects of climate change have on slope instability in mountainous and cryospheric terrain. The objective of this study is to identify slow-moving landslides in the Prince William Sound region, southcentral Alaska, United States, which has had rapid deglaciation since the mid-1800s, and assess their tsunamigenic plausibility. We use an automated time series persistent scatterer interferometric synthetic aperture radar processing method with 7 years of Sentinel-1 data (2016–22) to identify 43 slow-moving slopes with average velocities ranging from approximately 0.2 to 21 millimeters per year. Landslide presence is confirmed using aerial imagery and previous landslide inventory records. We assess the tsunamigenic plausibility of the landslides using empirically derived estimates of landslide mobility based on modeled landslide volumes. Of the identified landslides, our preliminary analysis suggests that 11 have tsunamigenic potential if they were to fail rapidly and catastrophically. Although our estimate of tsunamigenic plausibility is preliminary and can be refined with additional observations and analyses, it can be used to prioritize ongoing and future hazard assessment, surveillance, and research efforts.

Alaska↗

USGS2023

Steady-state forms of channel profiles shaped by debris flow and fluvial processes

Debris flows regularly traverse bedrock channels that dissect steep landscapes, but our understanding of bedrock erosion by debris flows and their impact on steepland morphology is still rudimentary. Quantitative models of steep bedrock channel networks are based on geomorphic transport laws designed to represent erosion by water-dominated flows. To quantify the impact of debris flow erosion on steep channel network form, it is first necessary to develop methods to estimate spatial variations in bulk debris flow properties (e.g., flow depth, velocity) throughout the channel network that can be integrated into landscape evolution models. Here, we propose and evaluate two methods to estimate spatial variations in bulk debris flow properties along the length of a channel profile. We incorporate both methods into a model designed to simulate the evolution of longitudinal channel profiles that evolve in response to debris flow and fluvial processes. To explore this model framework, we propose a general family of debris flow erosion laws where erosion rate is a function of debris flow depth and channel slope. Model results indicate that erosion by debris flows can explain the occurrence of a scaling break in the slope–area curve at low-drainage areas and that upper-network channel morphology may be useful for inferring catchment-averaged erosion rates in quasi-steady landscapes. Validating specific forms of a debris flow incision law, however, would require more detailed model–data comparisons in specific landscapes where input parameters and channel morphometry can be better constrained. Results improve our ability to interpret topographic signals within steep channel networks and identify observational targets critical for constraining a debris flow incision law.

Earth Surface Dynamics↗