April 2010 LIP of the Month

Short duration for the eruption of the Grande Ronde Basalt lavas of the Columbia River Basalts, USA.

Excerpts modified from paper in press in Lithos, entitled ‘New 40Ar/39Ar dating of the Grande Ronde lavas, Columbia River Basalts, USA: implications for duration of flood basalt eruption episodes’.

T.L. Barry1, S. Self1, S.P. Kelley1, S. Reidel2, P. Hooper1, M. Widdowson1

1Volcano Dynamics Group, CEPSAR, Open University, Walton Hall, Milton Keynes, MK7 6AA, UK.; t.l.barry@open.ac.uk

2 Washington State University, Tri-Cities Campus, Richland, WA 99354, USA.

Summary

Grande Ronde Basalt (GRB) lavas represent the most voluminous eruptive pulse of the Columbia River-Snake River-Yellowstone hotspot volcanism. With an estimated eruptive volume of 150,000 km3, GRB lavas form at least 66% of the total volume of the Columbia River Basalt Group (Fig. 1) [using the latest estimates of total CRB-related lava volume (Camp et al., 2003; Camp & Ross, 2004) compared with 90% of the volume of the CRBG as previously defined by Tolan et al. (1989)]. New 40Ar/39Ar dates for GRB lavas reveal they were emplaced within a maximum period of 0.42 ±0.18 My. A well-documented stratigraphy indicates at least 110 GRB flow fields (or individual eruptions), and on this basis suggests an average inter-eruption hiatus of less than 4,000 years. Isotopic age-dating cannot resolve time gaps between GRB eruptions, and it is difficult to otherwise form a picture of the durations of eruptions because of non-uniform weathering in the top of flow fields and a general paucity of sediments between GRB lavas. Where sediment has formed on top of the GRB, it varies in thickness from zero to 20-30 cm of silty to fine-sandy material, with occasional diatomaceous sediment. Individual GRB eruptions varied considerably in volume but many were greater than 1000 km3 in size. Most probably eruptive events were not equally spaced in time; some eruptions may have followed short periods of volcanic repose (perhaps 102 to 103 of yrs), whilst others could have been considerably longer (many 1000s to >104 yrs). Recent improvements in age-dating for other continental flood basalt (CFB) lava sequences have yielded estimates of total eruptive durations of less than 1 My for high-volume pulses of lava production. The GRB appears to be a similar example, where the main pulse occupied a brief period. Even allowing for moderate to long-duration pahoehoe flow field production, the amount of time the system spends in active lava-producing mode is small – less than c. 2.6% (based on eruption durations of approximately 10,000 yrs, as compared to the duration of the entire eruptive pulse of c. 420,000 yrs). A review of available 40Ar/39Ar data for the major voluminous phases of the Columbia River Basalt Group suggests that activity of the Steens Basalt-Imnaha Basalt-GRB may have, at times, been simultaneous, with obvious implications for climatic effects. Resolving intervals between successive eruptions during CFB province construction, and durations of main eruptive pulses, remains vital to determining the environmental impact of these huge eruptions.


Figure 1. Geographic extent of the Grande Ronde Basalt (grey field) within the wider limits of the whole Columbia River Basalt Group (including Steens Basalts), western USA. Sample location sites are: Vantage (V), Frenchmans Spring Coulee (FC), Wildhorse Monument (WM), Burr Canyon (BC), Palouse Falls Park (PF), Weissenfels Ridge (WR), Joseph Creek (JC), Maloney Creek (MC). Inset map: location of detailed map and the GRB within the USA.

Background

The CRBG is the Earth’s youngest continental flood basalt province with volcanism beginning around 17 Ma. Previous age determinations placed the eruption of the whole province between 17 and 6 Ma (Baksi, 1989; Tolan et al., 1989), with the voluminous GRB lavas estimated to have erupted between ~ 16.0 and ~ 15.0 Ma ago (Long and Duncan 1983; Hooper et al., 2002). The GRB lavas overlie the readily identifiable phenocryst-rich Imnaha Basalt lavas in the south of the province, and immediately underlie the Vantage Sandstone Member which forms part of the Ellensburg Formation (Swanson et al., 1979), that is distributed throughout much of the western CRBG on the Columbia Plateau. In many places throughout the province, the lavas of the Wanapum Basalt lie above the GRB, with or without the Vantage sediments between.

FORMATION

Member

Previous

K/Ar

ages (Ma)

Previous

40Ar/39Ar

ages (Ma)

This study:

sample

number

Material

analysed

Steps on plateau/

Total no. of steps

MSWD

This study:

40Ar/39Ar

ages (Ma)

WANAPUM

Frenchman Springs

Basalt of Sand Hollow

 

 

 

 

PF7

PF7 (rpt)

 

 

gl matrix

gl matrix

 

 

7/7

6/6

 

 

1.40

1.50

 

 

14.70 ±0.21

15.12 ±0.38

Eckler Mountain

 

CRB05-081

gl matrix

9/16

0.33

15.70 ±0.34

Vantage sediments

 

 

 

 

 

 

GRANDE RONDE

N2 group

7 - 10 members

 

 

15.0 ±0.4 (F)

16.0 ±1.4 (G)

 

  15.05 ±0.60 (A)

  15.36 ±0.60 (A)

  15.56 ±0.40 (A)

  15.79 ±0.40 (B)

  15.87 ±0.40 (A)

 

WH1a-1

WH1b-2

WH1b-1

 

CRB05-106

 

gl matrix

alt. matrix

seg. ves matrix

alt. matrix

 

6/8

6/6

6/8

 

11/14

 

1.80

0.95

1.80

 

0.85

 

15.48 ±0.22

15.55 ±0.32

15.71 ±0.26

 

15.94 ±0.20

R2 group

4 - 5 members

15.8 ±0.6 (F)

  14.95 ±0.60 (A)

  16.07 ±0.28 (C)

  16.32 ±0.20 (C)

CRB05-033

gl matrix

6/9

1.20

15.46 ±0.21

N1 group

2 - 3 members

15.7 ±0.6 (F)

15.9 ±0.6 (F)

 

 

 

 

 

 

 

R1 group

4 - 5 members

(# Birch Creek)

 

#14.49 ±0.80 (D)

#15.20 ±0.60 (D)

#16.00 ±0.20 (D)

 16.41 ±0.22 (C)

CRB05-030

CJC-7

gl matrix

gl matrix

2/8

11/13

0.62

1.05

15.53 ±0.20

16.25 ±0.27

IMNAHA

(# Upper Pole Creek)

 

 

 

 

 

15.7 ±0.6 (F)

17.2 ±0.4 (H)

 

# 15.30 ±5.60 (D)

  15.50 ±0.40 (D)

  15.60 ±0.60 (D)

# 15.60 ±1.20 (D)

  16.53 ±0.09 (E)

  16.73 ±0.14 (E)

# 16.70 ±1.60 (D)

# 16.90 ±0.60 (D)

  17.60 ±0.60 (A)

  17.67 ±0.32 (C)

MLC-2

MLC-2 (rpt)

CJC-3(1)

CJC-1

CJC-1 (rpt)

CJC-5

wr matrix

wr matrix

wr matrix

wr matrix

wr matrix

wr matrix

3/5

8/12

6/9

3/8

10/14

4/10

1.90

0.47

0.76

2.00

0.63

0.18

15.29 ±0.22

15.70 ±0.49

15.89 ±0.60

16.08 ±0.67

16.06 ±0.29

16.11 ±0.54

Table 1: Summary of published and new age determinations for the Imnaha Basalt, GRB and the base of the Wanapum Basalt. All new data and published data are reported with 2s errors. Formations are presented in stratigraphic order, with the oldest at the base. Directly underlying Imnaha Basalt is the Steens Basalt, and directly overlying the Wanapum Basalt is Saddle Mountains Basalt (not shown). Bracketed letters (A) to (H) refer to reference sources: (A) previously unpublished data; details available from S. Reidel; (B) Long and Duncan (1983); (C) Baksi and Farrar (1990); (D) Hooper et al. (2002); (E) Jarboe et al. (2006); (F) Watkins and Baksi (1974); (G) Snavely et al. (1973); and (H) McKee et al. (1977). All previously published Ar  data has been corrected to a Fish Canyon Tuff accepted age of 28.201 Ma (Kuiper et al., 2008) and therefore will appear differently from published text. # indicates samples that are from Malheur Gorge succession, eastern Oregon – Birch Creek Member and Upper Pole Creek Member are lateral equivalent of GR Basalt – R1 and Imnaha Basalt, respectively (Hooper et al., 2002). Data for other newly analysed samples that did not produce a plateau age are given along with the full data set in an online Electronic Appendix accompanying Barry et al. (2010). Abbreviations for material are: gl = glassy; wr = wholerock; alt = altered; seg.ves. = segregation vesicle.

 

Discussion
Duration of the Grande Ronde Basalt
On the basis of our new Ar data, we can assess individual ages for the top and bottom of the GRB lavas and calculate a period in which the GRB erupted. Furthermore, the interval in which the GRB volcanism occurred can be examined by its relationship with the Imnaha Basalt below and the Wanapum Basalt above.

A calculated weighted mean average age for all the previous data (with no exclusions) for the timing of the Imnaha Basalt is 16.61 ±0.28 Ma (MSWD=11.9; Fig. 2). [Weighted mean average ages are calculated using ISOPLOT and include standard error propagation]. The same calculation for our Imnaha Basalt data, excluding the poor analysis of MLC-2 (Table 1), gives a result of 15.99 ±0.20 Ma (MSWD=0.5; Fig. 2). This latter date is outside of error of the weighted mean average age for all published data for Imnaha Basalt, but is our best estimate for the top of the lava pile. The wide scatter in published ages for the Imnaha Basalt (15.30 ±5.60 to 17.67 ±0.32 Ma; Fig. 2) may be the result of samples representing different parts of the Imnaha lava pile and possible diachroneity within the Imnaha Basalt, as a result of complex volcanic architecture such as overlapping lava lobes of differing ages, rather than a single horizon. In addition, the variation in the results may reflect variable degrees of sample alteration or excess argon in some samples. Regardless of such complexities, we suggest that Imnaha Basalt volcanism in southern Washington State - where GRB was subsequently emplaced - finished at around 16.00 Ma. In the absence of any evidence for a significant hiatus between the Imnaha Basalt and the GRB, it seems likely that the climactic GRB volcanism began shortly thereafter. No field-based evidence for interstratification of Imnaha and GRB lavas has yet been found.

A calculation of all published data (Table 1) for N2 lavas at the top of the GRB succession reveals a weighted mean average age of 15.62 ±0.37 Ma (MSWD = 1.7; Fig. 2). The same calculation for our new data alone, excluding sample CRB05-106 for reasons given in the Results section, gives a date of 15.57 ±0.15 Ma (MSWD = 0.92; Fig. 2). Clearly, these results are well within error of each other and fit with data for the overlying Wanapum Basalt, e.g. Eckler Mountain Member (15.70 ±0.34 Ma) and basalt of Sand Hollow (15.12 ±0.38 Ma). Between the top of the GRB succession and the Wanapum Basalt are the Vantage sediments. It was once considered that the Vantage sediments were deposited in a 250, 000 year period, but it is likely that they formed in a much shorter interval, given constraints from the new argon data, the lack of any sediment in places, and the lack of evidence of erosion on the top surface of the GRB.

On the basis of this new evidence it would appear that the GRB lava succession started erupting shortly after 15.99 ±0.20 Ma and finished around 15.57 ±0.15 Ma. This suggests a total duration of 0.42 ±0.18 (1s) My, or approximately 420,000 years. Given uncertainties in the timing of the end of the Imnaha volcanism and that dates younger than 16 Ma are reported for the Imnaha (Table 1), the duration of the Grande Ronde volcanism could be considerably less than our estimate.


Figure 2. A comparison of available 40Ar/39Ar data for Grande Ronde Basalt, and Imnaha and Steens Basalts, plotted against the geomagnetic polarity timescale. All 40Ar/39Ar errors are plotted here as 1s, rather than the usual 2s used throughout the text, to aid clarity of the figure. Data from this study are plotted as solid black squares. Published data are plotted as solid white squares. * refers to silicic units that are interpreted to overlie or interfinger with Steens Basalt. Reference abbreviations: A-E are the same as for Table 1; plus additional data from (F) Henry et al., 2006; (G) Jarboe et al., 2008; (H) Brueseke et al., 2007; (I) Swisher et al., 1990; (J) Baksi et al., 1991. Geomagnetic timescale from Gradstein et al. (2004). For the GRB-N2 data and for Imnaha Basalt, the weighted mean average age of our data is shown as a heavy dashed line with 1s errors shown as the upper and lower limits of the light grey box (refer to text for full description). Similar calculations for published data are shown by the dotted line within the dark grey box.

Absolute age of the Grande Ronde Basalt
As described by Swanson et al. (1979), the GRB lava pile contains flows that have been demonstrated to have two normal and two reverse magnetizations (R1-N1-R2-N2). Two full cycles of polarity overturn within 420,000 years is more rapid than generally observed in the geomagnetic record, although similarly rapid overturns are recorded in detailed studies of the Oligocene/Miocene boundary (e.g., Billups et al., 2004).

Using the geomagnetic polarity timescale and Ar data for Steens Basalt and one Imnaha Basalt sample (Table 1), Jarboe et al. (2006) similarly suggest a rapid succession of eruption for these lava formations; they propose that the bulk of the Imnaha and GRB lavas erupted within a 0.75 My period. However, on the basis of Ar data and paleomagnetic results they suggest that the Steens Basalts erupted in geomagnetic chron C5Cr (Fig. 2; 17.235 to 16.721 Ma; Gradstein et al., 2004), and that Imnaha erupted within the chron C5Cn.3n (16.721 to 16.543 Ma; Gradstein et al., 2004). As a consequence of this and the geomagnetic timescale, Jarboe et al. (2008) neatly suggest that rapid reversals C5Cn.2r (16.543 to 16.472 Ma), C5Cn.2n (16.472 to 16.303 Ma), C5Cn.1r (16.303 to 16.268 Ma) and C5Cn.1n (16.268 to 15.974 Ma) are GRB R1-N1-R2-N2, respectively (Fig. 2). However, they do not have age data directly from the GRB succession to support this.

Chron C5Cn.3n ends at 16.543 Ma, and following Jarboe et al. would imply that Imnaha magmatism did not occur after this time. Yet there are many Ar dates suggesting younger ages than 16.50 Ma for Imnaha Basalt (Table 1; Fig. 2). Similarly with the GRB, chrons C5Cn.2r to C5Cn.1n are between approximately 16.50 and 16.00 Ma, yet almost all available Ar-age data for the GRB are younger than 16.00 Ma (Fig. 2). The spread of available age data for Steens Basalt, Imnaha Basalt and GRB make it difficult to be certain when a particular volcanic phase ended and suggests that activity may have even overlapped at times. However, the overriding evidence from all the available age data is that the GRB erupted after 16 Ma and therefore cannot be constrained to the C5Cn.2r to C5Cn.1n chrons, and that Imnaha Basalt eruption cannot be restricted to have occurred only in chron C5Cn.3n (Fig. 2). It is worth noting that very little paleomagnetic work has been carried out on the Imnaha Basalt, and it may yet prove to be more complex than initially thought (V.E. Camp, pers. comm., 2010). The age data even suggest that Imnaha Basalt volcanism could have been occurring for some time after GRB volcanism had begun. Potential evidence supporting the eruption of highly voluminous lavas such as Imnaha Basalt and GRB at around 16.0 Ma and 15.5 Ma, respectively, are dissolution events recorded in sediments off the coast of Africa. These dissolution events are thought to record climatic disturbances caused by changes in the oxygen and carbon dioxide composition of the atmosphere and would fit well with voluminous basaltic eruptions (Kender et al., 2009). In summary, with the lack of knowledge of the detailed stratigraphy and paleomagnetic signatures of the Imnaha lavas, it is unwise to speculate further on the exact absolute age of the GRB lava succession from their relationship with Imnaha Basalts.

Our age data, along with other published information suggest that GRB lavas most likely erupted within the younger geomagnetic chron C5Br (15.974 to 15.160 Ma; Gradstein et al., 2004; Fig. 2). We recognise that this poses a problem for the significance of the reversals within the GRB; there appears to be a lack of rapid magnetic reversals during C5Br (e.g. Gradstein et al., 2004). We have no solution to this dilemma at present, and can only suggest that either (a) the reversals measured in the GRB are short excursions within the geomagnetic cycle, or b) that rapid reversals C5Cn.2r to C5Cn.1n should be younger than 16 Ma. In support of younger geomagnetic chrons for Imnaha and GRB lavas than suggested by Jarboe et al. (2008) are 40Ar/39Ar dates for Steens Basalts that are younger than 16.60 Ma (Fig. 2).

Average estimates of eruption frequency
The stratigraphy of the 150,000 km3 GRB lavas has been well constrained and mapped in previous studies (e.g. Reidel et al., 1989). The GRB lavas built up from at least 110 individual eruptions, which over a 420,000 year time span, suggest an average periodicity of one eruption every ~4200 years. Assuming an eruption duration of 10 to 100 years per eruption, that would add up to between 1100 and 11,000 years total eruption time, within the 420,000 year time span. This would suggest that the total amount of time that the system was actively erupting, was no more than 2.6% of the total duration time of formation of the GRB, and quite probably much less than 1% of the time.

In more detail, the individual GRB eruptions have estimated eruption volumes between 90 and as much as 10,000 km3 (Reidel et al., 1989) – the larger estimates rival the outputs of even the greatest volume explosive rhyolitic super-eruption deposits such as at Yellowstone (c. 2500 km3; Mason et al., 2004). In reality, the eruptions would most likely have been irregularly distributed through time, possibly in pulses with some eruptions following another after a few 1000 years, while others may have erupted after, possibly, up to 10,000 years of quiescence. With eruptions lasting decades to perhaps centuries, the chances of distinguishing a hiatus that is intra-eruption versus inter-eruption becomes very difficult.

Conclusions

Our new 40Ar/39Ar data provide constraints for the duration of time in which the Grande Ronde Basalt, the climactic phase of the Columbia River Basalt Province, erupted. The Grande Ronde Basalt is bracketed by an age of 15.99 ±0.20 Ma for the top of the underlying Imnaha Basalt and by an age of 15.57 ±0.15 Ma for the end of Grande Ronde Basalt volcanism. A duration of approximately 420,000 years implies that, on average, there may have been an eruption every 4000 years. Based on published estimates for single eruption durations, we find that within a 420,000 year period less than 2.6%, and quite possibly less than 1% of that time, would have been volcanically active. A review of all the available 40Ar/39Ar data for the Steens Basalt-Imnaha Basalt-GRB suggests that activity of these petrochemically distinct groups may have, at times, been simultaneous, and if confirmed would have significant implications for potential environmental effects.

References

Baksi, A.K., 1989. Reevaluation of the timing and duration of extrusion of the Imnaha, Picture Gorge, and Grande Ronde Basalts, Columbia River basalt group. In: Reidel, S.P., and Hooper, P.R., eds: Geological Society of America Special Paper 239, 105-112.

Baksi, A.K. and Farrar, E., 1990. Evidence for errors in the geomagnetic polarity time-scale at 17-15 Ma: 40Ar/39Ar dating of basalt from the Pacific northwest, USA. Geophysical Research Letters 17, 1117-1120.

Baksi, A.K., Hall, C.H., York, D., 1991. Laser probe 40Ar/39Ar dating studies on sub-milligram whole-rock basalt samples: the age of the Steens Mountain geomagnetic polarity transition (revisited). Earth and Planetary Science Letters 104, 292-298.

Barry, T.L., Self, S., Kelley, S.P., Reidel, S., Hooper, P., Widdowson, M. 2010,  New 40Ar/39Ar dating of the Grande Ronde lavas, Columbia River Basalts, USA: Implications for duration of flood basalt eruption episodes. In Press, Available online 27 March 2010, doi:10.1016/j.lithos.2010.03.014

Billups, K., Pälike, H., Channell, J.E.T., Zachos, J.C., Shackleton, N.J., 2004. Astronomic calibration of the late Oligocene through early Miocene geomagnetic polarity time scale. Earth and Planetary Science Letters 224, 33-44.

Brueseke, M.E., Heizler, M.T., Mertzman, S.A., 2007. Distribution and geochronology of Oregon Plateau (USA) flood basalt volcanism: The Steens Basalt revisited. Journal of Volcanology and Geothermal Research 161, 187-214.

Camp, V.E. and Ross, M.E., 2004. Mantle dynamics and genesis of mafic magmatism in the intermontane Pacific northwest. Journal of Geophysical Research 109, B08204, doi:10.1029/2003JB002838.

Camp, V. E., Ross, M.E., Hanson, W.E., 2003. Genesis of flood basalts and basin and range volcanic rocks from Steens mountain to the Malheur River Gorge, Oregon. The Geological Society of America Bulletin 115 (1), 105-128.

Gradstein, F., Ogg, J., Smith, A., 2004. A geologic time scale. 589 pp., Cambridge University Press, Cambridge, UK.

Henry, C. D., Castor, S. B., McIntosh, W. C., Heizler, M. T., Cuney, M., Chemillac, R., 2006. Timing of oldest Steens Basalt magmatism from precise dating of silicic volcanic rocks, McDermitt Caldera and Northwest Nevada volcanic field. Eos Transactions American Geophysical Union 87 (52), Fall Meet. Suppl., Abstract V44C-08.

Hooper, P.R., Binger, G.B., Lees, K.R., 2002. Ages of the Steens and Columbia River flood basalts and their relationship to extension-related calc-alkalic volcanism in eastern Oregon. Geological Society of America Bulletin 114, 43-50.

Jarboe, N.A., Coe, R.S., Renne, P.R., Glen, J.M., 2006. 40Ar/39Ar ages of the early Columbia River Basalt Group: determining the Steens Mountain geomagnetic polarity reversal (R0-N0) as the top of the C5Cr chron and the Imnaha normal (N0) as the C5Cn.3n chro.. EOS Transactions, American Geophysical Union 87(52), Fall meeting supplement, Abstract V51D – 1702.

Jarboe, N.A., Coe, R.S., Renne, P.R., Glen, J.M.G., Mankinen,, E.A., 2008. Quickly erupted volcanic sections of the Steens Basalt, Columbia River Basalt Group: secular variation, tectonic rotation, and the Steens Mountain reversal. Geochemistry Geophysics Geosystems 9 (11), Q11010, doi:10.1029/2008GC002067.

Kender, S., Peck, V.L., Jones, R.W., Kaminski, M.A., 2009. Middle Miocene oxygen minimum zone expansion offshore West Africa: evidence of global cooling precursor events. Geology 37, 699-702.

Kuiper, K. F., Deino, A., Hilgen, F.J., Krijgsman, W., Renne, P.R., Wijbrans, J.R., 2008. Synchronizing rock clocks of Earth history. Science 320, 500-504.

Long, P.E., Duncan, R.A., 1983. 40Ar/39Ar ages of the Columbia River Basalt from deep boreholes in south-central Washington [abs], EOS Transactions American Geophysical Union 64, 90.

Mason, B.G., Pyle, D.M., Oppenheimer, C., 2004. The size and frequency of the largest explosive eruptions on Earth. Bulletin of Volcanology 66, 735-748, doi: 10.1007/s00445-004-0355-9.

McKee, E.H., Swanson, D.A., Wright, T.L., 1977. Duration and volume of Columbia River basalt volcanism, Washington, Oregon, and Idaho. Geological Society of America Abstracts with Programs 9, 463.

Reidel, S.P., Tolan, T.L., Hooper, P.R., Beeson, M.H., Fecht, K.R., Bentley, R.D., Anderson, J.L., 1989. The Grande Ronde Basalt, Columbia River Basalt Group; stratigraphic descriptions and correlations in Washington, Oregon, and Idaho. Geological Society of America Special Paper 239, 21-53.

Snavely, P.D., Jr., MacLeod, N.S., Wagner, H.C., 1973. Miocene tholeiitic basalts of coastal Oregon and Washington and their relations to coeval basalts of the Columbia Plateau. Geological Society of America Bulletin 84, 387-424.

Swanson, D.A., Wright, T.L., Hooper, P.R., Bentley, R.D., 1979. Revisions in stratigraphic nomenclature of the Columbia River Basalt Group. U.S. Geological Survey Bulletin 1457-G, 59 pp.

Swisher, C.C., Ach, J.A., Hart, W.K., 1990. Laser fusion 40Ar/39Ar dating of the type Steens Mountain Basalt, southeastern Oregon and the age of the Steens geomagnetic polarity transition. EOS Transactions American Geophysical Union 71, 1296.

Tolan, T.L., Reidel, S.P., Beeson, M.H., Anderson, J.L., Fecht, K.R., Swanson, D.A., 1989. Revisions to the estimates of the areal extent and volume of the Columbia River Basalt Group. In: Reidel, S.P. and Hooper, P.R., (eds), Volcanism and tectonism in the Columbia River flood-basalt province, Geological Society of America Special Paper 239, 1-20.

Watkins, N.D., Baksi, A.K., 1974. Magnetostratigraphy and oroclinal folding of the Columbia River, Steens, and Owyhee Basalts in Oregon, Washington, and Idaho. American Journal of Science 274, 148-189.