ABSTRACT Emplacement models for voluminous sheet flows of the Columbia River flood basalts vary significantly in style and duration, with the latter ranging from as little as one week to decades and even centuries. Testing the efficacy of such models requires detailed field studies and close examination of each stratigraphic unit. The Steens Basalt, the oldest formation of the Columbia River flood basalts, differs from the later formations in that it is composed of stacked successions of thin, commonly inflated flow lobes combined into thicker compound flows, or flow fields. These flow lobes are of limited geographic extent, with relatively high emplacement rates, but they are otherwise similar to modern examples. Evidence for flow inflation in the much larger sheet flows of the Grande Ronde Basalt, Wanapum Basalt, and Saddle Mountains Basalt is also apparent, but with more variable rates of emplacement. For example, the Asotin and Umatilla Members (Saddle Mountains Basalt) and Sentinel Bluffs Member flows (Grande Ronde Basalt) erupted distinct compositions along their linear vent systems, but over 200 km west of their vents, these flows are no longer distinct. Instead, they exist as compositional zones of a single, moderately mixed lava flow. Such flows must have been emplaced rapidly, in perhaps weeks to months, while others have been shown to erupt over much longer time periods. We conclude that emplacement rates may be quite variable throughout the Columbia River flood basalt province, with thin flow units of Steens Basalt erupting continuously and rapidly, and larger inflated sheet flows erupting over variable time spans, some from a few weeks to months, and others over a duration of years.

Abstract The Middle Miocene Columbia River Basalt Group (CRBG) is the youngest and smallest continental flood basalt province on Earth, covering over 210,000 km 2 of Oregon, Washington, and Idaho and having a volume of 210,000 km 3 . A well-established regional stratigraphic framework built upon seven formations, and using physical and compositional characteristics of the flows, has allowed the areal extent and volume of the individual flows and groups of flows to be calculated and correlated with their respective dikes and vents. CRBG flows can be subdivided into either compound flows or sheet flows, and are marked by a set of well-defined physical features that originated during their emplacement and solidification. This field trip focuses on the Lewiston Basin, in southeastern Washington, western Idaho, and northeastern Oregon, which contains the Chief Joseph dike swarm, where classic features of both flows and dikes can be easily observed, as well as tectonic features typical of those found elsewhere in the flood basalt province.

Abstract The Moscow-Pullman basin, located on the eastern margin of the Columbia River flood basalt province, consists of a subsurface mosaic of interlayered Miocene sediments and lava flows of the Imnaha, Grande Ronde, Wanapum, and Saddle Mountains Basalts of the Columbia River Basalt Group. This sequence is ~1800 ft (550 m) thick in the east around Moscow, Idaho, and exceeds 2300 ft (700 m) in the west at Pullman, Washington. Most flows entered from the west into a topographic low, partially surrounded by steep mountainous terrain. These flows caused a rapid rise in base level and deposition of immature sediments. This field guide focuses on the upper Grande Ronde Basalt, Wanapum Basalt, and sediments of the Latah Formation. Late Grande Ronde flows terminated midway into the basin to begin the formation of a topographic high that now separates a thick sediment wedge of the Vantage Member to the east of the high from a thin layer to the west. Disrupted by lava flows, streams were pushed from a west-flowing direction to a north-northwest orientation and drained the basin through a gap between steptoes toward Palouse, Washington. Emplacement of the Roza flow of the Wanapum Basalt against the western side of the topographic high was instrumental in this process, plugging west-flowing drainages and increasing deposition of Vantage sediments east of the high. The overlying basalt of Lolo covered both the Roza flow and Vantage sediments, blocking all drainages, and was in turn covered by sediments interlayered with local Saddle Mountains Basalt flows. Reestablishment of west-flowing drainages has been slow. The uppermost Grande Ronde, the Vantage, and the Wanapum contain what is known as the upper aquifer. The water supply is controlled, in part, by thickness, composition, and distribution of the Vantage sediments. A buried channel of the Vantage likely connects the upper aquifer to Palouse, Washington, outside the basin. This field guide locates outcrops; relates them to stratigraphic well data; outlines paleogeographic basin evolution from late Grande Ronde to the present time; and notes structures, basin margin differences, and features that influence upper aquifer water supply.

Continental flood basalt provinces are the subaerial expression of large igneous province volcanism. The emplacement of a continental flood basalt is an exceptional volcanic event in the geological history of our planet with the potential to directly impact Earth's atmosphere and environment. Large igneous province volcanism appears to have occurred episodically every 10–30 m.y. through most of Earth history. Most continental flood basalt provinces appear to have formed within 1–3 m.y., and within this period, one or more pulses of great magma production and lava eruption took place. These pulses may have lasted from 1 m.y. to as little as a few hundred thousand years. Within these pulses, tens to hundreds of volumetrically large eruptions took place, each producing 10 3 –10 4 km 3 of predominantly p3hoehoe lava and releasing unprecedented amounts of volcanic gases and ash into the atmosphere. The majority of magmatic gas species released had the potential to alter climate and/or atmospheric composition, in particular during violent explosive phases at the eruptive vents when volcanic gases were lofted into the stratosphere. Aside from the direct release of magmatic gases, magma-sediment interactions featured in some continental flood basalt provinces could have released additional carbon, sulfur, and halogen-bearing species into the atmosphere. Despite their potential importance, given the different nature of the country rock associated with each continental flood basalt province, it is difficult to make generalizations about these emissions from one province to another. The coincidence of continental flood basalt volcanism with periods of major biotic change is well substantiated, but the actual mechanisms by which the volcanic gases might have perturbed the environment to this extent are currently not well understood, and have been little studied by means of atmospheric modeling. We summarize current, albeit rudimentary, knowledge of continental flood basalt eruption source and emplacement characteristics to define a set of eruption source parameters in terms of magmatic gases that could be used as inputs for Earth system modeling studies. We identify our limited knowledge of the number and length of non-eruptive phases (hiatuses) during continental flood basalt volcanism as a key unknown parameter critical for better constraining the severity and duration of any potential environmental effects caused by continental flood basalt eruptions.

The middle Miocene Columbia River Basalt Group is the youngest and smallest continental flood basalt province on Earth, covering over 210,000 km 2 of mainly Oregon, Washington, and Idaho, with an estimated basalt volume of ~210,000 km 3 . A well-established regional stratigraphic framework built upon six formations contains numerous flows and groups of flows that can be readily distinguished by their physical and compositional characteristics, thus producing mappable units, the areal extent and volume of which can be calculated and correlated with their respective feeder dikes. The distinct physical features that help to define these units originated during their emplacement and solidification, as the result of variations in cooling rates, degassing, thermal contraction, and interaction with their paleoenvironment. Columbia River Basalt Group flows can be subdivided into two basic flow geometries. Sheet flows dominate the basalt pile, but the earliest flows comprising the Steens Basalt and some of the Saddle Mountains Basalt flows are compound flows with elongated bodies composed of numerous, local, discontinuous, and relatively thin lobes of basalt lava. The internal physical characteristics of the voluminous sheet flows are recognizable throughout their extent, thus allowing mechanistic models to be developed for their emplacement. The emplacement and distribution of individual Columbia River Basalt Group flows resulted from the interplay among the regional structure, contemporaneous deformation, eruption rate, preexisting topography, and the development of paleodrainage systems. These processes and their associated erosional and structural features also influenced the nature of late Neogene sedimentation during and after the Columbia River Basalt Group eruptions. In this paper, we describe and revise the stratigraphic framework of the province, provide current estimates on the areal extent and volume of the flows, and summarize their physical features and compositional characteristics.

The Columbia River flood basalt province, United States, is likely the most well-studied, radiometrically well-dated large igneous province on Earth. Compared with older, more-altered basalt in flood basalt provinces elsewhere, the Columbia River Basalt Group presents an opportunity for precise, accurate ages, and the opportunity to study relationships of volcanism with climatic excursions. We critically assess the available 40 Ar/ 39 Ar data for the Columbia River Basalt Group, along with K-Ar data, to establish an up-to-date picture of the timing of emplacement of the major formations that compose the lava stratigraphy. Combining robust Ar-Ar data with field constraints and paleomagnetic information leads to the following recommendations for the age of emplacement of the constituent formations: Steens Basalt, ca. 16.9 to ca. 16.6 Ma; Imnaha Basalt, ca. 16.7 to ca. 16 Ma; Grande Ronde Basalt, ca. 16 Ma to ca. 15.6 Ma; Wanapum Basalt, ca. 15.6 to ca. 15 Ma; and Saddle Mountains Basalt from ca. 15 Ma to ca. 6 Ma. The results underline the previously held observation that Columbia River Basalt activity was dominated by a brief, voluminous pulse of lava production during Grande Ronde Basalt emplacement. Under scrutiny, the data highlight areas of complexity and uncertainty in the timing of eruption phases, and demonstrate that even here in the youngest large igneous province, argon dating cannot resolve intervals and durations of eruptions.

The radiometric dating evidence for the timing and duration of volcanism for the Steens through Wanapum Basalt of the Columbia River Basalt Group is critically reviewed here. K-Ar dates generally underestimate the age of crystallization, though one important exception is detected, where excess argon led to dates that were too old. The 40 Ar/ 39 Ar results on whole-rock basalts from 1980 through 2010 are examined for statistical validity of plateau sections, as well as alteration state of the material dated. In most instances, listed ages are shown to be invalid. The 40 Ar/ 39 Ar total gas (fusion) ages are, in general, not accurate estimates of the time of formation of these rocks. The 40 Ar/ 39 Ar ages on plagioclase separates from basalts yield good estimates of the extrusion age of the lavas. New 40 Ar/ 39 Ar ages on whole-rock basalts are presented that are in good agreement with the plagioclase ages. Various forms of the geomagnetic polarity time scale for mid-Miocene time are examined, along with the ages of lavas and their magnetic polarity. The main sections of the Columbia River Basalt Group (Imnaha through Wanapum Basalt) were formed in ~0.5 m.y. between 16.3 and 15.8 Ma. Steens Basalt extrusion occurred about ~0.1 m.y. before the Imnaha Basalt and appears to have been a precursor to the more voluminous volcanism noted in the Columbia River Basalt Group.

We examined Grande Ronde Basalt lava flows from surface sections and boreholes throughout Washington, Oregon, and Idaho to determine chemical and physical properties that would allow the recognition and mapping of these flows on a regional scale. We estimate there are ~100 flows covering nearly 170,000 km 2 , with a total volume of ~150,400 km 3 , that were erupted over four polarity intervals (reverse 1, normal 1, reverse 2, and normal 2) in ~0.42 m.y. These flows are the largest known on Earth, with individual volumes ranging from ~100 km 3 to greater than 10,000 km 3 . Although all known Grande Ronde Basalt flows erupted in the eastern part of the Columbia River flood basalt province, the thickest and most complete sections (>3 km) occur in the central Columbia Basin. From the center of the basin, the number of flows decreases outward, resulting in a nearly complete stratigraphy in the interior and an abbreviated and variable stratigraphy along the margins. The areal extent of many flows suggests that the Chief Joseph dike swarm greatly expanded after Imnaha Basalt time, and now many dikes are buried beneath younger flows in the eastern part of the province. The Grande Ronde Basalt has a relatively uniform lithology with only a few distinctive flows. However, when compositions are combined with paleomagnetic polarity, lithology, and stratigraphic position, the Grande Ronde Basalt can be subdivided into at least 25 mappable units. Grande Ronde Basalt flows are siliceous, with typically SiO 2 >54 wt%, MgO contents ranging from ~2.5 to 6.5 wt%, and TiO 2 ranging from 1.6 to 2.8 wt%, with an enrichment in iron and incompatible elements relative to mid-ocean-ridge basalt. Although most Grande Ronde Basalt flows have homogeneous compositions, some are heterogeneous. Dikes that fed the heterogeneous flows show that the first composition erupted was not typical of the flow, but as the eruption progressed, the compositions gradually evolved to the bulk composition of flow. The average effusion rate was ~0.3 km 3 /yr, with basalt volume peaking during the R2 polarity with the eruption of the Wapshilla Ridge Member. Eruption and emplacement rates for the flows are controversial, but available data collected from the field suggest that many of the flows could have been emplaced in a few years to perhaps a decade.

The Columbia River system is one of the great river systems of North America, draining much of the Pacific Northwest, as well as parts of the western United States and British Columbia. The river system has had a long and complex history, slowly evolving over the past 17 m.y. The Columbia River and its tributaries have been shaped by flood basalt volcanism, Cascade volcanism, regional tectonism, and finally outburst floods from Glacial Lake Missoula. The most complex part of river development has been in the northern part, the Columbia Basin, where the Columbia River and its tributaries were controlled by a subsiding Columbia Basin with subtle anticlinal ridges and synclinal valleys superimposed on a flood basalt landscape. After negotiating this landscape, the course to the Pacific Ocean led through the Cascade Range via the Columbia Trans-Arc Lowland, an ancient crustal weakness zone that separates Washington and Oregon. The peak of flood basalt volcanism obliterated the river paths, but as flood basalt volcanism waned, the rivers were able to establish courses within the growing fold belt. As the folds grew larger, the major pathways of the rivers moved toward the center of the Columbia Basin where subsidence was greatest. The finishing touches to the river system, however, were added during the Pleistocene by the Missoula floods, which caused local repositioning of river channels.

Newly generated and previously published strontium isotope signatures of plagioclase phenocrysts in Columbia River Basalt Group lavas exhibit heterogeneity largely imposed by mantle-derived magmas assimilating variable crustal rocks. Steens basalts assimilated accreted terrane crust with 87 Sr/ 86 Sr ratios of <0.7040. In contrast, Imnaha, Grande Ronde, Wanapum, and Saddle Mountains basalts likely assimilated crust with more radiogenic Sr (>0.7040) including a cratonic component, perhaps as a result of residence in magma chambers partly located east of the 87 Sr/ 86 Sr = 0.7060 line. Strontium isotope ratios in plagioclase phenocrysts from early- erupted Imnaha basalts anticipate whole-rock signatures of later-erupted Grande Ronde basalts consistent with a geochemical continuum between the two formations, which is also seen in whole-rock trace element abundances, undermining the notion of an abrupt change in the magma sources generating Imnaha and Grande Ronde basalts.

The Grande Ronde Basalt lavas constitute ~63% of the Columbia River Basalt Group, a large igneous province in the NW United States. The lavas are aphyric or contain less than 5% phenocrysts of plagioclase, augite, pigeonite, and olivine (altered). Plagioclase hygrometry shows that the erupted lavas generally contained less than 0.3% dissolved H 2 O; however, the presence of rare disequilibrium An 96 plagioclase phenocrysts suggests that some magmas may have originally had 4.5 wt% dissolved H 2 O at depth, but they all degassed during ascent and eruption. The size of plagioclase phenocrysts suggests an average plagioclase phenocryst residence time in the magmas of 160 yr. Ignoring hiatuses between eruptions, we estimate that the ~110 flows of the Grande Ronde Basalt erupted over a cumulative time of 17,600 yr, with an average eruption rate of ~8.6 km 3 /yr. The average interval between eruptions is estimated to be 3658 yr. It is envisaged that a shallow intrusive dike-sill complex, rather than large kilometer-sized magma chamber(s), fed the Grande Ronde basalt lavas. We performed model simulations using the COMAGMAT software to retrace the pre-eruption histories of the Grande Ronde Basalt lavas. Based on such simulations and petrological reasoning, we find that the primary melts could have been generated from a spinel peridotite source at 1.5 GPa, perhaps under hydrous conditions. Extensive melting of lithospheric eclogite may have played an important role as well; however, this is not constrained by our simulations. All lavas were contaminated by the crust, and they were last processed (mixing, crystallization) during their short residence within shallow (6 km) intrusive rocks prior to eruption. Our petrologic conclusions lead us to present a petrotectonic model that supports the hypothesis that the Columbia River Basalt Group magma generation was greatly aided by a thinned lithosphere and H 2 O that may have come off the asthenospheric wedge.

We present a reappraisal of Columbia River basalt petrogenesis based on an internally consistent X-ray fluorescence and inductively coupled plasma–mass spectrometry data set for major and trace elements plus new and existing isotopic analyses of the Imnaha, Steens, Picture Gorge, and Grande Ronde Basalts. Source materials for the main-phase Columbia River Basalt Group are upwelling ocean-island basalt source–like mantle, depleted mantle variably fluxed by slab-derived fluids, Phanerozoic arc crust, and ancient North American cratonic crust. The mantle upwelling may be a deep-seated plume or material displaced and mobilized by fragmented sinking slabs. We endorse the conclusions of earlier workers that the Imnaha, Steens, and Picture Gorge Basalts represent different mixtures of upwelling mantle, depleted mantle, slab-derived fluids, and crust. Cratonic crust of the Idaho batholith has a minor role as a contaminant of Imnaha basaltic magma and a major role in the petrogenesis of the Grande Ronde basaltic andesites, which we model as contaminated Imnaha basalt. We emphasize the geochemical continuity of the Imnaha and Grande Ronde Basalts and propose that they derive from a single central crustal magma system (or chamber) that lasted from ca. 16.7 Ma to 16.0 Ma. The Imnaha–Grande Ronde magma system was centered beneath the location where the western Snake River Plain, Oregon-Idaho graben, and Chief Joseph dike swarm converge, and probably transgressed the cratonic boundary to the east. Magma from this system was injected into dikes and traveled hundreds of kilometers northward to erupt and feed the gigantic Grande Ronde lava flows. In contrast to previous studies, we question the idea that the dike swarm provides a geographic map of the magma source regions.

The Lewiston Structure is located in southeastern Washington and west-central Idaho and is a generally east-west–trending (~075°), asymmetric, noncylindrical anticline in the Columbia River Basalt Group that transfers displacement into the Limekiln fault system to the southeast and the Silcott fault system to the southwest. A serial cross-section analysis and three-dimensional (3-D) construction of this structure show how the fold varies along its trend and shed light on the deformational history of the Lewiston Basin. Construction of the fold’s 3-D form shows that the fold’s wavelength increases and amplitude decreases near its eastern and western boundaries. Balanced cross sections show ~5% shortening across the structure, which is consistent with the Yakima Fold Belt. An angular unconformity below the Grande Ronde Basalt N1 magnetostratigraphic unit, in addition to a variation of N1 unit thickness across the structure, suggests that the fold was forming before N1 time. Analysis of structural data using the Gauss method for heterogeneous fault-slip data indicates north-south (~350°) shortening prior to and after N1 emplacement. The presence of a reverse fault on the southern limb of the Lewiston Structure is controversial. This fault crops out to the east of the field area where Grande Ronde Basalt magnetostratigraphic unit R2 is thrust over Pliocene(?) gravels. However, better control on unit thicknesses and map contacts rules out an exposed reverse fault on the southern limb of the fold west of the Washington-Idaho border, suggesting the fault either dies out or becomes blind.