Bedrock erosion and canyon formation during extreme floods have dramatically altered landscapes on Earth and Mars. Grand Coulee was carved by outburst floods from Pleistocene glacial Lake Missoula and is the largest canyon in the Channeled Scabland, a megaflood‐scoured landscape in the northwestern USA. Quantifying paleo‐discharge is required to understand how landscapes evolve in response to extreme events, but there are few constraints on the magnitude of the floods that incised Grand Coulee; hence, we used hydraulic modeling and geologic evidence to quantify paleo‐flood discharges during different phases of canyon incision. When upper Grand Coulee was incising by headward waterfall retreat, the paleo‐discharge was 2.6 × 106 m3 s−1, which produced shear stresses great enough to cause the waterfall to retreat via toppling of basalt columns. The largest possible flood through upper Grand Coulee, a Missoula flood which raised glacial Lake Columbia to a stage of 750 m, produced a modeled discharge of 7.6 × 106 m3 s−1. The discharges associated with waterfall retreat and drainage of glacial Lake Columbia are >80% and ∼50% lower, respectively, than the 14–17 × 106 m3 s−1 discharge predicted by assuming the present‐day topography was inundated to the elevation of high‐water marks. Due to bedrock incision, high‐water marks may overestimate paleo‐flow depth in canyons carved by floods, hence bedrock erosion should be considered when estimating paleo‐discharge in flood‐carved canyons. Our results indicate that outburst floods with discharges and flow depths much lower than those required to inundate high‐water marks are capable of carving deep canyons.

During the last deglaciation, dozens of glacial outburst floods—among the largest known floods on Earth—scoured the Channeled Scabland landscape of eastern Washington. Over this same period, deformation of the Earth’s crust in response to the growth and decay of ice sheets changed the topography by hundreds of meters. Here, we investigated whether glacial isostatic adjustment affected routing of the Missoula floods and incision of the Channeled Scabland from an impounded, glacial Lake Columbia. We used modern topography corrected for glacial isostatic adjustment as an input to flood models that solved the depth-averaged, shallow water equations and compared the results to erosion constraints. Results showed that floods could have traversed and eroded parts of two major tracts of the Channeled Scabland—Telford-Crab Creek and Cheney–Palouse—near 18 ka, whereas glacial isostatic adjustment limited flow into the Cheney–Palouse tract at 15.5 ka. Partitioning of flow between tracts was governed by tilting of the landscape, which affected the filling and overspill of glacial Lake Columbia directly upstream of the tracts. These results highlight the impact of glacial isostatic adjustment on megaflood routing and landscape evolution.

New findings about old puzzles occasion rethinking of the Grand Coulee, greatest of the scabland channels. Those puzzles begin with antecedents of current upper Grand Coulee. By a recent interpretation, the upper coulee exploited the former high-level valley of a preflood trunk stream that had drained to the southwest beside and across Coulee anticline or monocline. In any case, a constriction and sharp bend in nearby Columbia valley steered Missoula floods this direction. Completion of upper Grand Coulee by megaflood erosion captured flood drainage that would otherwise have continued to enlarge Moses Coulee. Upstream in the Sanpoil valley, deposits and shorelines of last-glacial Lake Columbia varied with the lake’s Grand Coulee outlet while also recording scores of Missoula floods. The Sanpoil evidence implies that upper Grand Coulee had approached its present intake depth early the last glaciation at latest, or more simply during a prior glaciation. An upper part of the Sanpoil section provides varve counts between the last tens of Missoula floods in a stratigraphic sequence that may now be linked to flood rhythmites of southern Washington by a set-S tephra from Mount St. Helens. On the floor of upper Grand Coulee itself, recently found striated rock and lodgement till confirm the long-held view, which Bretz and Flint had shared, that cutting of upper Grand Coulee preceded its last-glacial occupation by the Okanogan ice lobe. A dozen or more late Missoula floods registered as sand and silt in the lee of Steamboat Rock. Some of this field evidence about upper Grand Coulee may conflict with results of recent two-dimensional simulations for a maximum Lake Missoula. In these simulations only a barrier high above the present coulee intake enables floods to approach high-water marks near Wenatchee that predate stable blockage of Columbia valley by the Okanogan lobe. Above the walls of upper Grand Coulee, scabland limits provide high-water targets for two-dimensional simulations of watery floods. The recent models sharpen focus on water sources, prior coulee incision, and coulee’s occupation by the Okanogan ice lobe. Field reappraisal continues downstream from Grand Coulee on Ephrata fan. There, some of the floods exiting lower Grand Coulee had bulked up with fine sediment from glacial Lake Columbia, upper coulee till, and a lower coulee lake that the fan itself impounded. Floods thus of debris-flow consistency carried outsize boulders previously thought transported by watery floods. Below Ephrata fan, a backflooded reach of Columbia valley received Grand Coulee outflow of small, late Missoula floods. These late floods can—by varve counts in post-S-ash deposits of Sanpoil valley—be clocked now as a decade or less apart. Still farther downstream, Columbia River gorge choked the largest Missoula floods, passing peak discharge only one-third to one-half that released by the breached Lake Missoula ice dam.

In May 2012, a sediment‐laden flood along the Seti Khola (= river) caused 72 fatalities and widespread devastation for > 40 km in Pokhara, Nepal's second largest city. The flood was the terminal phase of a hazard cascade that likely began with a major rock‐slope collapse in the Annapurna Massif upstream, followed by intermittent ponding of meltwater and subsequent outburst flooding. Similar hazard cascades have been reported in other mountain belts, but peak discharges for these events have rarely been quantified. We use two hydrodynamic models to simulate the extent and geomorphic impacts of the 2012 flood and attempt to reconstruct the likely water discharge linked to even larger medieval sediment pulses. The latter are reported to have deposited several cubic kilometres of sediment in the Pokhara Valley. The process behind these sediment pulses is debated. We traced evidence of aggradation along the Seti Khola during field surveys and from RapidEye satellite images. We use two steady‐state flood models, HEC‐RAS and ANUGA, and high‐resolution topographic data, to constrain the initial flood discharge with the lowest mismatch between observed and predicted flood extents. We explore the physically plausible range of simplified flood scenarios, from meteorological (1000 m3 s−1) to cataclysmic outburst floods (600,000 m3 s−1). We find that the 2012 flood most likely had a peak discharge of 3700 m3 s−1 in the upper Seti Khola and attenuated to 500 m3 s−1 when arriving in Pokhara city. Simulations of larger outburst floods produce extensive backwater effects in tributary valleys that match with the locations of upstream‐dipping medieval‐age slackwater sediments in several tributaries of the Seti Khola. Our findings are consistent with the notion that the medieval sediment pulses were linked to outburst floods with peak discharges of >50,000 m3 s−1, though discharge may have been an order of magnitude higher.

No abstract