Simulations of cataclysmic outburst floods from Pleistocene Glacial Lake Missoula

Using a flow domain that we constructed from 30 m digital-elevation model data of western United States and Canada and a two-dimensional numerical model for shallow-water flow over rugged terrain, we simulated outburst floods from Pleistocene Glacial Lake Missoula. We modeled a large, but not the largest, flood, using initial lake elevation at 1250 m instead of 1285 m. Rupture of the ice dam, centered on modern Lake Pend Oreille, catastrophically floods eastern Washington and rapidly fills the broad Pasco, Yakima, and Umatilla Basins. Maximum flood stage is reached in Pasco and Yakima Basins 38 h after the dam break, whereas maximum flood stage in Umatilla Basin occurs 17 h later. Drainage of these basins through narrow Columbia gorge takes an additional 445 h. For this modeled flood, peak discharges in eastern Washington range from 10 to 20 × 10⁶ m³/s. However, constrictions in Columbia gorge limit peak discharges to