Missoula Floods Research Articles

The timing of Missoula floods: Implications for the age of Grand Coulee (eastern Washington, USA)

The Channeled Scabland of eastern Washington (USA) was formed by outburst floods from glacial Lake Missoula. Despite chronological advances, the timing of erosion in the main flood channels is unresolved. In particular, it is still uncertain whether upper Grand Coulee, the largest canyon in the Channeled Scabland, was incised during or prior to the last glaciation. We report 10Be exposure ages from erratics in upper Grand Coulee, glacial Lake Columbia, and surrounding flood routes. Flood-transported boulders on the high-elevation east rim of Grand Coulee date to ca. 17−15 ka. Ages from boulders on the floor of Grand Coulee indicate later flooding at ca. 14 ka, which post-dated canyon incision and occurred after inundation of the Telford-Crab Creek scabland at ca. 15−14.5 ka. Prior hydraulic modeling and dating suggest the entrance to Grand Coulee was blocked by rock and that canyon incision was incomplete at ca. 17 ka; hence, we interpret the 17−15 ka exposure ages on the east rim to coincide with flow over a retreating cataract during canyon incision. Our results indicate incision of Grand Coulee was completed between 17 ka and 14 ka. The short duration of canyon incision suggests that glacial Lake Missoula generated some of the most erosive outburst floods in Earth’s history.

Okanogan lobe tunnel channels and subglacial floods into Moses Coulee, Channeled Scabland, northwestern United States

Abstract Outburst floods from glacial Lake Missoula largely explain erosion of the Channeled Scabland, a system of overfit, basaltic channels in Washington, northwestern United States. However, it is challenging to explain Missoula flood routing into Moses Coulee due to its topographic isolation from flood routes. To clarify flood pathways into Moses Coulee, we mapped channels that delineate a radial-anastomosing network connecting to Moses Coulee. Channels consist of coulee-like features eroded mainly in basalt. Channels climb adverse slopes and cross divides, and the network is draped with eskers and recessional moraines. These channel geometries and glacial landform associations suggest a tunnel channel network. Large channel dimensions and network anastomosis indicate formation in subglacial floods. The network connects to Moses Coulee with sufficient cross-sectional area to convey megafloods, hinting that subglacial floods may have been a significant source of Moses Coulee meltwater, in addition to possible diverted Missoula floods.

New Insights from Legacy Seismic Data regarding Basalt Elevations and Variability on the Hanford Site

Abstract Migration of groundwater contaminants in the Gable Gap area of the Hanford Site in southeastern Washington State is strongly influenced by the distribution and permeability of basalts that lie beneath an unconfined aquifer. Locally, folding and faulting of the Columbia River Basalt associated with the Yakima fold and thrust belt followed by erosion due to the Lake Missoula floods resulted in a complex basalt surface that represents either an impermeable lower boundary to the unconfined aquifer system or localized regions of increased permeability that potentially promote communication between the unconfined aquifer system and deeper, confined aquifer systems. Paleo-channels carved into the basalt by floodwaters are thought to provide preferential flow paths for groundwater contaminants. In 2011, a seismic landstreamer campaign was carried out to image the basalt surface and produced pre-stack depth migrated p-wave reflection images. The reflection images identified two large troughs that may represent paleo-channels and several areas of possible faulting. Here, the streamer data are re-analyzed using refraction travel-time and Rayleigh wave dispersion analyses to obtain images of compressional and shear wave velocities within the suprabasalt sediment sections and the upper basalt surface. The combined interpretation of reflection and seismic velocity images shows complexity in the basalt velocity and elevation, which varies by 50 m or more within the study area. These results, along with other ongoing geophysical investigations, will be used to inform the site geologic model and potentially guide placement of future boreholes needed to quantify vertical flow between the confined and unconfined aquifers.

Glacial isostatic adjustment directed incision of the Channeled Scabland by Ice Age megafloods

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.

Pleistocene Megaflood Discharge in Grand Coulee, Channeled Scabland, USA

AbstractBedrock 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.

Glacial processes and landforms

AbstractFrom 1965 to 2000 glacial geomorphology became increasingly specialized and developed significantly due to technological improvements, particularly in remote sensing, surveying and field-based glaciological process studies. The better understanding of basal thermal regimes in ice sheets and glaciers led to the development of concepts such as spatial and temporal migration of ice divides in dynamic ice sheets that could overprint subglacial landform assemblages, debris entrainment processes related to polythermal glacier systems, and glacier and ice sheet beds composed of cold and warm based mosaics. Process observations at the ice–bed interface led to the discovery of the third glacier flow mechanism, substrate deformation, which provided the impetus to reconstruct the genesis of subglacial bedforms such as drumlins and to evaluate the origins and potential flow law for till. Numerical evaluations of glacial erosion led to a better understanding of abrasion and quarrying as well as the erection of genetic models and erosion rates for larger-scale features such as U-shaped valleys and cirques. Linkages were made between debris transport pathways and moraine construction in supraglacial environments, with the role of glacier structure being linked to specific landforms, such as medial, lateral, hummocky and ice-cored moraines as well as rock glaciers. Our appreciation of the erosional and depositional impacts of glacifluvial systems was enhanced significantly with the advent of process observations on the hydrology of modern glaciers as well as the final vindication of J.H. Bretz and his proposed jökulhlaup origins of the Channelled Scablands and the Missoula Floods. In addition to the increasing numbers of studies at modern glacier snouts, the embracing of sedimentology by glacial geomorphologists was to result in significant developments in understanding the process-form regimes of subglacial, marginal and proglacial landforms, particularly the recognition of landform continua and hybrids. Advances resulting from this included the recognition of different modes of moraine and glacitectonic thrust mass development, lithofacies models of the varied glacifluvial depositional environments, and the initial expansion of work on the sediments and depocentres of glacimarine settings, the latter being the result of glacial research taking to submersibles and ice-strengthened ships for the first time. A similarly new frontier was the expansion of research on the increasingly higher resolution images returning from Mars, where extraterrestrial glaciations were recognized based on comparisons with Earth analogues. Holistic appreciations of glaciation signatures using landform assemblages were developed, initially as process-form models and later as glacial landsystems, providing an ever-expanding set of templates for reconstructing palaeoglaciology in the wide variety of topographic and environmental settings, which also acknowledge spatial and temporal change in glacier and ice sheet systems.

Chronology of Missoula Flood Deposits at the Coyote Canyon Mammoth Site, Washington State, USA

Late Pleistocene outburst megafloods, mostly from glacial Lake Missoula, hydraulically ponded behind downstream constrictions in the Columbia River in southeastern Washington State, USA. Optically stimulated luminescence (OSL) ages for flood deposits from the Coyote Canyon Mammoth Site, located in a high (315–320 m asl) distal portion of the transient lake, indicate that at least seven floods ponded high enough to inundate the area during the period 20.9 ± 2.6 ka to 16.3 ± 2.8 ka. This is consistent with a radiocarbon age of about 17.4 ± 0.2 ka cal BP from the middle of the flood sequence. OSL ages from loess deposits overlying a paleosol at the top of the flood sequence range from 14.0 ± 2.3 ka to 10.9 ± 2.0 ka, suggesting a hiatus of about 2.3 thousand years. These datasets are consistent with current understanding that multiple late Pleistocene megafloods occurred between 20 ka and 14 ka and that earlier floods produced higher flood stages than later ones. The lack of flood deposits in the Coyote Canyon area younger than 16 ka supports the hypothesis that younger megafloods did not hydraulically pond in the Pasco Basin above about 230 m asl.

Empirical evidence for cosmogenic 3He production by muons

Cosmic ray muons penetrate deeply into rock where they interact with atoms to produce cosmogenic nuclides. Incorporation of the muon contribution to the production rates of cosmogenic nuclides increases the accuracy of exposure dates, burial ages, and erosion rates inferred from measured nuclide concentrations. In the absence of empirical evidence, it is generally assumed that muons do not produce 3He, a cosmogenic nuclide commonly used for exposure dating. Here we assess whether muons produce 3He by measuring He isotope concentrations in pyroxene and ilmenite from a ∼300 m deep drill core and other subsurface samples of the mid-Miocene Columbia River Basalt in Washington, USA. 3He concentrations in our samples exhibit an exponential decline with depth with an e-folding length of 32.4 m, which corresponds to an attenuation length for 3He production of 8780 g cm−2. The deeply penetrating exponential is diagnostic of 3He production by cosmic ray muons. Assuming no erosion, we constrain the minimum surface muonogenic production rate to be 0.23 atom g−1 pyroxene yr−1, whereas when incorporating erosion the production rate is 0.45 atom g−1 pyroxene yr−1. 3He concentrations in samples deeper than ∼100 meters are consistent with model-based estimates of depth-independent nucleogenic production from the capture by 6Li of neutrons produced by alpha particle reactions on light elements. Measurements in other subsurface samples indicate that muon-produced 3He is prevalent across the Columbia Plateau. The penetration depth of muonogenic 3He production is substantially deeper, and the ratio of muon- to spallation-produced 3He is substantially lower, than found for other cosmogenic nuclides. Our results provide the first definitive empirical evidence for 3He production by muons, which has several implications for quantifying the timing and rates of Earth surface change and interpreting He isotope ratios. Importantly, despite the low production rates, landforms in the Channeled Scablands, which were formed by incision of the Columbia River Basalt by the late-Pleistocene Missoula floods, have high concentrations of 3He inherited from post-Miocene muon exposure. Hence 3He production by muons must be considered, particularly when dating rapid erosional events in old bedrock. Our findings indicate samples with less than several tens of meters of shielding by overlying rock will contain cosmogenic 3He that elevates 3He/4He ratios. Hence caution should be used when using 3He/4He ratios from samples at shallower depths to infer mantle sources of basalt.

Snake River: A Navigation and Spawning Dilemma

The Snake River in northwestern United States is 1735 km long, the largest tributary of the Columbia River and is the 13th longest river in the United States. The Snake River drainage basin includes parts of six U.S. states. The Snake River Plain was created by a volcanic hotspot that lies beneath Yellow-stone National Park. The previous Ice Age carved out canyons, cliffs and waterfalls along the middle and lower Snake River. The Missoula Flood was to the north and Bonneville Flood to the south altered the Snake River and surrounding landscape. The Snake River has a drainage basin of 282,000 km2 in the states of Oregon, Washington, Utah, Nevada and Idaho. The Snake River drops from mountain elevations of 3000 m to its confluence with the Columbia River. The river is one of the most biologically diverse freshwater systems in the United States with trails designed to promote recreational tourism, increase use of the Snake River and create generations of people who care about the river and are willing to protect and provide environmental stewardship of the river watershed resources.

Erratics and Other Evidence of Late Wisconsin Missoula Outburst Floods in Lower Wenatchee and Adjacent Columbia Valleys, Washington

The Pleistocene Missoula floods through eastern and central Washington are by peak fl ow rate (discharge) the greatest freshwater cataclysms known on Earth. Newly explored features along the Wenatchee reach of Columbia valley give new evidence and revise earlier interpretations of size, frequency, and routing of megafloods.Crystalline-rock erratics derived far northeast lie scattered about the sandstone hills of lower Wenatchee valley and adjacent Columbia valley up to 495 m altitude, 320 m above Columbia River. They can only have been ice-rafted by flood(s) running down the Columbia. Before the late Wisconsin Okanogan lobe of Cordilleran ice blocked the Columbia, at least one monstrous Missoula flood poured down the valley past Wenatchee and backflooded Wenatchee valley.Rhythmically bedded sandy silt in Columbia valley between Trinidad and Wenatchee records repeated silt-rich backfloods up the Columbia from Quincy basin—after Okanogan-lobe ice had blocked the Columbia upvalley. Rhythmically graded silt beds in Wenatchee valley at Dryden containing Columbia-derived dropstones record ten Missoula backfloods up the valley.Thick silt farther up Wenatchee valley between Peshastin and Leavenworth had been thought deposits of a long-lived lake, dammed supposedly by the Malaga landslide. But the heights and distribution of provable lake beds now make Moses Coulee bar the only viable dam—and only up to altitude 275 m. The silt above Dryden lying at 315–385 m altitude must also have been laid by Missoula floods.

Pre-Late Wisconsin Valley-Glacier Erratics between Leavenworth and Peshastin, Wenatchee Valley, Washington

The late Wisconsin Icicle Creek alpine glacier transported tonalite boulders from the Mount Stuart batholith to arcuate end moraines in Icicle valley and Wenatchee valley near Leavenworth. Some previous workers considered sparsely weathered Mount Stuart boulders lying outside these moraines and draped by silt as having been ice rafted in a late Wisconsin lake. But the boulders are more likely within drift or end moraines of pre-late Wisconsin glaciations correlative to lateral moraines on Boundary Butte. The silt—a younger deposit that covers the boulders—derives not from a physically ponded lake but apparently from brief Missoula flood(s) backflooding up the valley.

10Be dating of late Pleistocene megafloods and Cordilleran Ice Sheet retreat in the northwestern United States

During the late Pleistocene, multiple floods from drainage of glacial Lake Missoula further eroded a vast anastomosing network of bedrock channels, coulees, and cataracts, forming the Channeled Scabland of eastern Washington State (United States). However, the timing and exact pathways of these Missoula floods remain poorly constrained, thereby limiting our understanding of the evolution of this spectacular landscape. Here we report cosmogenic ^(10)Be ages that directly date flood and glacial features important to understanding the flood history, the evolution of the Channeled Scabland, and relationships to the Cordilleran Ice Sheet (CIS). One of the largest floods occurred at 18.2 ± 1.5 ka, flowing down the northwestern Columbia River valley prior to blockage of this route by advance of the Okanogan lobe of the CIS, which dammed glacial Lake Columbia and diverted later Missoula floods to more eastern routes through the Channeled Scabland. The Okanogan and Purcell Trench lobes of the CIS began to retreat from their maximum extent at ca. 15.5 ka, likely in response to onset of surface warming of the northeastern Pacific Ocean. Upper Grand Coulee fully opened as a flood route after 15.6 ± 1.3 ka, becoming the primary path for later Missoula floods until the last ones from glacial Lake Missoula at 14.7 ± 1.2 ka. The youngest dated flood(s) (14.0 ± 1.4 ka to 14.4 ± 1.3 ka) came down the northwestern Columbia River valley and were likely from glacial Lake Columbia, indicating that the lake persisted for a few centuries after the last Missoula flood.

Megafloods and Clovis cache at Wenatchee, Washington

Immense late Wisconsin floods from glacial Lake Missoula drowned the Wenatchee reach of Washington's Columbia valley by different routes. The earliest debacles, nearly 19,000 cal yr BP, raged 335 m deep down the Columbia and built high Pangborn bar at Wenatchee. As advancing ice blocked the northwest of Columbia valley, several giant floods descended Moses Coulee and backflooded up the Columbia past Wenatchee. Ice then blocked Moses Coulee, and Grand Coulee to Quincy basin became the westmost floodway. From Quincy basin many Missoula floods backflowed 50 km upvalley to Wenatchee 18,000 to 15,500 years ago. Receding ice dammed glacial Lake Columbia centuries more—till it burst about 15,000 years ago. After Glacier Peak ashfall about 13,600 years ago, smaller great flood(s) swept down the Columbia from glacial Lake Kootenay in British Columbia. The East Wenatchee cache of huge fluted Clovis points had been laid atop Pangborn bar after the Glacier Peak ashfall, then buried by loess. Clovis people came five and a half millennia after the early gigantic Missoula floods, two and a half millennia after the last small Missoula flood, and two millennia after the glacial Lake Columbia flood. People likely saw outburst flood(s) from glacial Lake Kootenay.

Record of glacial Lake Missoula floods in glacial Lake Columbia, Washington

During the last glaciation (marine oxygen isotope stage 2), outburst floods from glacial Lake Missoula deposited diagnostic sediments within glacial Lake Columbia. Two dominant outburst flood lithofacies are present within glacial Lake Columbia deposits: a flood expansion bar facies and a finer-grained hyperpycnite facies. We conclude that the flood sediments have a glacial Lake Missoula source because: (1) current indicators indicate westward flow through the lake, and upvalley flow followed by downvalley flow in tributary valleys; (2) no flood sediments are found north of a certain point; (3) there is a dominance of Belt–Purcell Supergroup clasts in a flood expansion bar; and (4) some of the finer-grained beds have a pink colour, reflective of glacial Lake Missoula lake-bottom sediments. A new radiocarbon age of 13,400 ± 100 14C BP on plant detritus found below 37 flood beds helps constrain the timing of outburst flooding from glacial Lake Missoula.

Pre-and post-Missoula flood geomorphology of the Pre-Holocene ancestral Columbia River Valley in the Portland forearc basin, Oregon and Washington, USA

Geomorphic landscape development in the pre-Holocene ancestral Columbia River Valley (1–5km width) in the Portland forearc basin (~50km length) is established from depositional sequences, which pre-date and post-date the glacial Lake Missoula floods. The sequences are observed from selected borehole logs (150 in number) and intact terrace soil profiles (56 in number) in backhoe trenches. Four sequences are widespread, including (1) a vertically aggraded Pleistocene alluvial plain, (2) a steep sided valley that is incised (125–150m) into the Pleistocene gravel plain, (3) Missoula flood terraces (19–13ka) abandoned on the sides of the ancestral valley, and (4) Holocene flooding surfaces (11–8ka) buried at 70–30m depth in the axial Columbia River Valley. Weathering rims and cementation are used for relative dating of incised Pleistocene gravel units. Soil development on the abandoned Missoula flood terraces is directly related to terrace deposit lithology, including thin Bw horizons in gravel, irregular podzols in sand, and multiple Bw horizons in thicker loess-capping layers. Radiocarbon dating of sand and mud alluvium in the submerged axial valley ties Holocene flooding surfaces to a local sea level curve and establishes Holocene sedimentation rates of 1.5cmyear−1 during 11–9ka and 0.3cmyear−1 during 9–0ka. The sequences of Pleistocene gravel aggradation, river valley incision, cataclysmic Missoula flooding, and Holocene submergence yield complex geomorphic landscapes in the ancestral lower Columbia River Valley.

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 6 m 3 /s. However, constrictions in Columbia gorge limit peak discharges to 6 m 3 /s and greatly extend the duration of flooding. We compare these model results with field observations of scabland distribution and high-water indicators. Our model predictions of the locations of maximum scour (product of bed shear stress and average flow velocity) match the distribution of existing scablands. We compare model peak stages to high-water indicators from the Rathdrum-Spokane valley, Walulla Gap, and along Columbia gorge. Though peak stages from this less-than-maximal flood model attain or exceed peak-stage indicators along Rathdrum-Spokane valley and along Columbia gorge, simulated peak stages near Walulla Gap are 10–40 m below observed peak-stage indicators. Despite this discrepancy, our match to field observations in most of the region indicates that additional sources of water other than Glacial Lake Missoula are not required to explain the Missoula floods.

New stratigraphic markers in the late Pleistocene Palouse loess: novel fossil gastropods, absolute age constraints and non-aeolian facies

Four stratigraphic sections in the southern part of the Columbia Basin preserve a sequence of aeolian and non-aeolian sediments ranging in age from 9AE43 to >47AE0 14 C ka based on accelerator mass spectrometry radiocarbon dating of fossil molluscs, geochemistry of Cascade Mountain-sourced tephra and association with formally recognized pedostratigraphic units (the Washtucna and Old Maid Coulee soils). Study sections are interpreted as representing concurrent deposition of loess and distal Missoula Flood rhythmites in valleys tributary to main drainages backflooded during the Missoula Floods, and formation of carbonate and iron-rich soils. Sediments belong to the formally recognized L-1 and L-2 loess units established for the Palouse loess, which were deposited in the Columbia Basin subsequent to events of glacial outburst flooding. Sediments associated with the Mount Saint Helens set S and set C tephras in the study sections preserve a fauna of five species of gastropod mollusc which have not been reported previously from sediments of late Pleistocene age in the Palouse region. The fossils comprise two distinct faunules stratigraphically separated by the Mount Saint Helens So tephra. Accelerator mass spectrometry radiocarbon dating of the fossils collected above the tephra in two of the sections yielded ages of 12AE48 ± 0AE06 and 9AE43 ± 0AE05 14 C kyr. These ages suggest that independent determinations of the 13AE35 14 C kyr age of the So tephra in other areas where Missoula Flood sediments are preserved are probably accurate, and help to refine the age of the latest events in the most recent sequence of catastrophic glacial outburst flooding.

The Spokane Flood debates: historical background and philosophical perspective

Abstract The 1920s–1930s debates over the origin of the ‘Channeled Scabland’ landscape of eastern Washington, northwestern USA, focused on the cataclysmic flooding hypothesis of J Harlen Bretz. During the summer of 1922, Bretz began leading field parties of advanced University of Chicago students into the region. In his first paper, published in the Bulletin of the Geological Society of America , Bretz took special care not to mention cataclysmic origins. However, in a subsequent paper in the Journal of Geology , to the editorial board of which he had recently been added, Bretz formally described his hypothesis that an immense late Pleistocene flood, which he named the ‘Spokane Flood’, had derived from the margins of the nearby Cordilleran Ice Sheet. This cataclysm neatly accounted for numerous interrelated aspects of the Channeled Scabland landscape and nearby regions. Nevertheless, the geological community largely resisted Bretz's hypothesis for decades, despite his enthusiastic and eloquent defence thereof. Resolution of the controversy came gradually, initially through the recognition by J. T. Pardee of a plausible source for the flooding: ice-dammed Pleistocene glacial Lake Missoula in northern Idaho and western Montana. Eventually, by the 1960s, the field evidence for cataclysmic flooding became overwhelming, and physical processes were found to be completely consistent with that evidence. The controversy is of philosophical interest in regard to its documentation of the attitudes of geologists toward hypotheses, which illustrate aspects of geological reasoning that are distinctive in degree from those of other sciences.

Cataclysmic Scabland flooding: Insights from a simple depth-averaged numerical model

This work presents numerical assessments of possible routing paths for the reported late Pleistocene Missoula floods, which involved Earth's largest known peak discharges of freshwater. For our numerical model, we adopt the diffusion wave approximation of the shallow-water equation with the empirical Manning coefficient. This simplification is required because (1) most previous studies are one-dimensional or analytical studies based on simple equations, so that our two-dimensional expansion is a natural step towards more sophisticated models; (2) a simple model can easily incorporate the effect of complicated topography; and (3) lower computational costs permit a broader exploration of parameters. Important and new insights obtained from this work include: (1) numerically, we confirm the presence of the Cordilleran Ice Sheet at the time of the Missoula floods; (2) floods with relatively small peak discharges cannot inundate some major reaches, including the Cheney Palouse Scabland Tract; (3) deposits indicative of multiple floods are mostly confined in the locations that can be inundated by relatively lower peak flood discharges; and (4) although Lake Missoula glacial lake failure scenario can reasonably reproduce extensive water coverage, the best fit result to the geological evidence of flooding is obtained when the total water volume is three times larger than that estimated for Glacial Lake Missoula. The above results do not negate previous works, but indicate the need for further detailed field investigation and more sophisticated modeling in order to pursue further understanding of the complex behavior of megaflooding in the Channeled Scabland.

Numerical simulations of large-scale cataclysmic floodwater: A simple depth-averaged model and an illustrative application

A simple numerical simulation code is developed to quantitatively discuss the behaviors of past cataclysmic floods, whose rigorous hydraulic descriptions would be too complex to be supported by typically available field data. The code aims to keep the basic equations simple and the required computational costs low, which allows greater exploration of parameter space. Examination of both the qualitative meaning of turbulent equations and the relative importance of terms in the shallow water equation have resulted in the adoption of the depth-averaged diffusion wave approximation in two dimensions, as well as Manning's empirical equation for simplification. On the other hand, topography is included as a boundary condition, since topography plays an important role in selecting flood routes. The model provides important advantages over theoretical analyses or one-dimensional simulations, including: (1) the reproduction of complicated water flow paths such as bifurcations and reconvergences; (2) the reconstruction of hydrological relationships among the water paths; and (3) the direct comparison among the calculated extent of the flood inundation (including floodwater depths) and geological and geomorphological observations. Preliminary reconstructions of the Missoula floods demonstrate the advantages of this new approach. The method of coupling field-based flood information with two-dimensional analysis of an entire flow path provides a basis for estimating the peak discharges and flow durations of cataclysmic megafloods on Earth with direct bearing on the understanding of regional paleohydrological histories of Mars.