July 2017 LIP of the Month

Mafic magmatism in the Siletz terrane, NW North America: Fragments of the youngest oceanic plateau?

Bethan A. Phillips1*, Andrew C. Kerr1*, Emily K. Mullen2, Dominique Weis2

1 School of Earth and Ocean Sciences, Cardiff University, Park Place, Cardiff, UK, CF10 3AT

2 Pacific Centre for Isotopic and Geochemical Research, Department of Earth, Ocean and Atmospheric Sciences, University of British Columbia, Vancouver, BC V6T 1Z4, Canada

*phillipsba@cf.ac.uk; kerra@cf.ac.uk


This article is an abridged version of the following Open Access paper:

Phillips, B.A., Kerr, A.C., Mullen, E.K., Weis, D. 2017. Oceanic mafic magmatism in the Siletz terrane, NW North America: Fragments of an Eocene oceanic plateau? Lithos, 274-275, 291-303.




  1. Introduction

The Siletz terrane, a mostly Eocene large igneous province (LIP) located in the present Cascadia forearc region of Oregon, Washington and Vancouver Island (Fig. 1), consists of accreted subaerial massive flows, submarine basaltic pillow lavas, and intrusive sheets (Snavely et al., 1968; Wells et al., 2014), with an estimated magmatic volume of 2.6 × 106 km3 (Trehu et al., 1994). Eruption and intrusion of the Siletz terrane magmas are thought to have taken place from ~56 to 49 Ma (Duncan, 1982; Haeussler et al., 2000; Hirsch and Babcock, 2009; Massey, 1986; Wells et al., 2014) during major plate reorganisation in the Pacific (Atwater, 1989; Engebretson, 1985). The province has been variously known as: Siletzia (Irving, 1979), the Coast Range Volcanic Province (Babcock et al., 1992), the Siletz–Crescent terrane (McCrory and Wilson, 2013) and the Siletz terrane (Fleming and Trehu, 1999; Snavely et al., 1968) used here.

There has been little consensus on the tectonic setting in which the Siletz terrane magmas were generated. Previously proposed models include initial melting of the Yellowstone mantle plume head (Pyle et al., 2009; Wells et al., 2014), accretion of a hotspot island chain (Duncan, 1982; McCrory and Wilson, 2013), Kula ridge subduction producing marginal rifting (Massey, 1986; Wells et al., 1984) and formation in a slab window-related setting (Babcock et al., 1992; Breitsprecher et al., 2003; Groome et al., 2003; Haeussler et al., 2003).

In this study samples have been collected and analysed from throughout the terrane and this article reports the results of a subset of the samples that have been analysed for radiogenic isotopes (Phillips et al. 2017).

Fig. 1. Map of the extent and exposed areas of the Siletz Terrane, adapted from McCrory and Wilson (2013). HR – Hurricane Ridge Fault.

  1. Geological history and tectonic setting

The ~ 650 km-long Siletz Terrane consists of three main volcanic formations: the Siletz River Volcanics in Oregon, the Crescent Formation in southwestern Washington and the Olympic Peninsula, and the Metchosin Igneous Complex on the southeastern coast of Vancouver Island (Fig. 1). The lower part of the terrane sequence contains mostly submarine pillow basalts, pillow breccia, gabbros and occasional sheeted dyke units (Duncan, 1982; Haeussler et al., 2000; Massey, 1986; Snavely et al., 1968; Wells et al., 2000). The upper part of the sequence contains both submarine pillow basalts and subaerial massive flows with oxidised flow tops along with rare felsic flows (Duncan, 1982; Haeussler et al., 2000; Massey, 1986; Snavely et al., 1968; Wells et al., 2000), that are interbedded with brecciated lapilli tuffs, laminated tuffs and volcaniclastic sediments. The thickest part of the terrane in central Oregon, reaches 27 ± 5 km (Trehu et al., 1994). Beneath northwestern Oregon and southwestern Washington, the terrane is approximately 20 km thick (Parsons et al., 1999), thinning to 10 km in northwestern Washington and southwestern Vancouver Island and down to 6 km offshore Vancouver Island (Hyndman et al., 1990).

The age of the Siletz terrane has generally been constrained using 40Ar-39Ar whole rock dating techniques, placing the maximum age range of the terrane between 56 and 49 Ma, with the majority of the terrane having erupted between 54 and 50 Ma (Brandon, 2014; Duncan, 1982; Pyle et al., 2009). In addition to this, limited U–Pb ages have also been determined. Haeussler et al. (2000) reported a U–Pb zircon age of 50.5 Ma for the Bremerton Igneous Complex in the Crescent Formation and U–Pb zircon ages of 52 and 54 Ma have been determined for the volcanics of the Metchosin Igneous Complex by Yorath et al. (1999)..

Accretion is thought to have occurred rapidly after, and probably even during, the eruption of the terrane, first in the south between 51 and 49 Ma (Wells et al., 2014) and subsequently in the north from 48 Ma to 45 Ma (McCrory and Wilson, 2013; Wells et al., 2014).

At the time of formation (Early Eocene) there were 4 (or 5) plates in the region; the Farallon, Kula, Pacific, North American (and Resurrection) plates (Fig. 2) (Haeussler et al., 2003; Seton et al., 2012). During this time a major adjustment in spreading direction occurred between the Pacific and Farallon plates from WSW–ENE to E– W (Atwater, 1989). The Pacific-Kula ridge also underwent a change in spreading direction from N–S to NW–SE (Seton et al., 2012) before an eventual cessation of spreading (Engebretson, 1985). Most early Eocene plate reconstruction models agree that the Farallon–Kula ridge was striking NE–SW at the time of formation of the Siletz terrane, intersecting the North American plate to the east, forming a triple junction with the Pacific plate in the west and a slab window at the adjacent margin of the North American continent (Atwater, 1989; Breitsprecher et al., 2003; Engebretson, 1985; Madsen et al., 2006; Seton et al., 2012).

Fig. 2. Suggested plate configuration map at 53 Ma. Adapted from Haeussler et al. (2003), Seton et al. (2012) and Wells et al. (2014), where the stars represent the proposed positions of the Yellowstone hotspot; 1 (O'Neill et al. (2005), 2 (Doubrovine et al., 2012) and 3 (Müller et al., 1993).

The Kula-Resurrection ridge probably subducted by 50 Ma (Breitsprecher et al., 2003; Haeussler et al., 2003; Wells et al., 2014) and the Siletz terrane is suggested to have formed on, or adjacent to, the mid-ocean ridge located between the Resurrection plate and the Farallon plate (McCrory and Wilson, 2013). It has also been argued that at this time (and at least since 60 Ma) the Yellowstone hotspot was located offshore of the North American continent; however, its location remains unclear (McCrory and Wilson, 2013). O'Neill et al. (2005) proposed that at 55 Ma the Yellowstone hotspot was located on the Kula–Farallon or Farallon–Resurrection ridge. Alternatively, Doubrovine et al. (2012) suggested that the hotspot was located off-axis on the Farallon plate whereas Müller et al. (1993) argued that it was adjacent to the North American continent and had moved underneath it by 50 Ma.

  1. Geochemical results

4.1. Major and trace elements

The silica content of the analysed samples varies between 39.7 and 74.6 wt.% with the majority of samples having < 50 wt.% SiO2, while the Mg# varies between 18.5 and 70.8. The Cr abundance varies between 3.4 and 576 ppm and the Ni ranges between 3.1 and 1244 ppm. The LOIs of the samples are mostly between 0.5 and 2% but extend up to ~ 6%, which, along with the petrography, indicates that varying degrees of alteration have affected the rocks of the Siletz terrane. Under such conditions many elements, in particular the large-ion lithophile elements (i.e., K2O, Ba, Rb), may become mobile (e.g., Hastie et al., 2007; Pearce, 1996).

The samples have therefore been classified using the Zr/Ti vs. Nb/Y diagram (Pearce, 1996), and it it can be seen that the majority of samples in this study are tholeiitic basalts while a small number of samples, largely from the Siletz River Volcanics, trend towards alkaline compositions (Fig. 3).

Fig. 3. Zr/Ti–Nb/Y classification diagram adapted from Pearce (1996).

In Fig. 4 representative trace elements and ratios are plotted against Zr as it is relatively immobile and incompatible over the range of observed compositions (Cann, 1970). Many of the more immobile elements such as Nb, Th and La generally show a positive correlation with Zr (Fig. 3), and so are a more robust representation of the original magmatic variations among the rocks. Other, more mobile, trace elements such as Sr display a positive correlation in most of the samples but much more scatter is evident.

Fig. 4. Bivariate plots of trace element variations with Zr (ppm) for the Siletz terrane. Also shown in grey is the Ontong Java Plateau basalts compositional field (Fitton and Godard, 2004).

On a primitive mantle normalised diagram (Fig. 5), samples from the Siletz River Volcanics and Crescent South and some samples from Roseburg show enrichment in the most incompatible elements, while the remainder of samples show generally flat patterns. Chondrite-normalised REE patterns (Fig. 5) display flat to slightly depleted HREE and middle REE patterns (Fig. 6). The LREE in the Siletz samples show flat-to-enriched patterns, with several samples from the Siletz River Volcanics displaying the most enriched LREE patterns (La/Sm up to 9.8). There is little variation from east to west in the Dy/Yb ratio of the rocks of the Siletz terrane. The highest Dy/Yb values are found in the central part of the terrane (Fig. 6) and the lowest values are seen in the Metchosin Igneous Complex, the northernmost samples of the Crescent Terrane and from the Roseburg area in the southernmost part of the terrane. Some of these low Dy/Yb samples (from the Metchosin Igneous Complex) are also LREE depleted (Fig. 6).

Fig. 5. Primitive mantle and chondrite-normalised trace element diagrams (normalising values from Sun and McDonough, 1989) of samples also analysed for isotopic ratios. For comparison the range of compositions of Ontong Java Plateau basalts is also shown (Fitton and Godard, 2004).

Fig. 6. a) 206Pb/204Pb (initial) vs. longitude, showing an overall decrease from West to East. b) 176Hf/177Hf (initial) vs. Dy/Yb which displays a noticeable inverse correlation. c) Latitude vs. ?Nd indicating an overall increase in variability southward across the terrane. d) Latitude vs. Dy/Yb ratio, showing little variation overall with the very lowest values corresponding to the northernmost part of the terrane.

4.2. Sr-Nd-Hf-Pb radiogenic isotopes

The measured isotope ratios for the samples used in this study have been age corrected for in-situ decay between 54.5 and 50.5 Ma (Fig. 7. The Metchosin Igneous Complex samples display the most depleted and restricted isotope ratios of the entire Siletz terrane, and this is particularly evident in Hf–Nd isotopic space (Fig. 7d). The Crescent Formation samples show more variation in isotope ratios and this variation increases southward. The more-depleted Crescent Formation samples are still generally significantly more enriched than Pacific MORB compositions (Fig. 7). Samples CM230 and CM225 have elevated 87Sr/86Sr values compared to the other samples, which may be related to seawater alteration as all other isotope values for these samples are not anomalous. In particular Nd and Hf isotope systems are relatively resistant to alteration and so are likely to represent the primary composition of the rocks (White and Patchett, 1984; Nobre Silva et al., 2010).

Fig. 7. Initial isotope ratios of the Siletz terrane samples. Also plotted are age-corrected (53 Ma) isotope ratios for mantle components (Supplementary Material 7), PMM - Pacific MORB mantle (Chauvel and Blichert-Toft, 2001), DGM - Depleted Gorda MORB mantle (Davis et al., 2008; Salters et al., 2011), IM - Imnaha mantle component (source of the primary Columbia River Flood basalts and thought to represent the Yellowstone hotspot; (Wolff et al., 2008)), CS - Cobb seamount chain (a series of seamounts located near the East Pacific Rise; (Chadwick et al., 2014)), and KFB - Karmutsen flood basalts of the Wrangellia terrane, Vancouver Island (Greene et al., 2009b). Also shown are compositional fields of the Ontong Java Plateau (OJP) (Mahoney et al., 1993a, 1993b; Tejada et al., 2002, 2004), Caribbean Plateau (CP) (Hastie et al., 2008; Hauff et al., 2000; Kerr et al., 1997, 2002; Thompson et al., 2004), Kerguelen Plateau (KP) (Weis and Frey, 2002; Weis et al., 2002), the Hawaiian Islands (HI) (Garcia et al., 2010; Tanaka et al., 2002; Weis et al., 2011) and the Samoan Islands (SI) (Salters et al., 2011; Workman et al., 2004). Error bars for the Siletz terrane data are smaller than the symbols. Binary mixing curves between DGM and the Imnaha mantle component, EM2 and HIMU mantle components are also shown, where each increment represents 10% addition of the mixing component to the DGM.

The samples from northern Oregon and southern Washington (the Crescent South group of the Crescent formation, and the Siletz River Volcanics) show a similar amount of isotopic variation to one another, with the majority of samples from this area falling in the main cluster of data. The samples from the Roseburg area, in the south of the terrane, show the most isotopic variation. For example, Sample 073RB, which has the highest 206Pb/204Pb ratio of 19.668 and so plots closest to the HIMU mantle component in 206Pb/204Pb diagrams (Fig. 7), also has the highest La/Sm ratio of 9.8. Sample 095RB, also from the Roseburg area of the terrane, plots closest to EM2 and the Imnaha mantle component (IM) fields in all isotopic spaces (Fig. 7). The Imnaha component is thought to be the mantle plume source of the Columbia River Flood Basalts (CRFB) and the Snake–Yellowstone volcanics (Wolff et al., 2008).

There is little correlation between major and trace elements and isotope abundances with the exception of the more trace element-depleted samples (in comparison to the other more-enriched samples) (Fig. 5) also having the most depleted isotopic ratios of the terrane (Figs. 6 and 7). Samples with lower Pb and Sr isotopic ratios and higher εHf and εNd values also have the lowest La/Sm ratios (0.9–1.3) and the most-depleted incompatible element abundances. There is also a negative correlation between Dy/Yb and Hf isotopic ratios (Fig. 6). As this trend is distinguishable in both trace elements and radiogenic isotopes it indicates that the differences observed are a reflection of the composition of the source regions rather than varying degrees of partial melting alone. The overall variation in isotopes generally increases southward across the terrane (Fig. 6).

All samples analysed have higher Pb isotopic ratios and 87Sr/86Sr values and lower 176Hf/177Hf and 143Nd/144Nd than the Pacific MORB Mantle (PMM) (Fig. 7). The Siletz terrane samples also generally overlap the Karmutsen Flood Basalt (KFB) field, which are a series of predominantly tholeiitic basalts located on Vancouver Island and form part of the accreted Jurassic aged Wrangellia oceanic plateau (Greene et al., 2009b) (Fig. 7). In Pb-Pb isotopic space the Siletz terrane shares similar isotopic characteristics to the Caribbean Plateau (Fig. 7).

  1. Discussion

5.1. Mantle source composition

Although the samples analysed in this study span a range of compositions there is no correlation between the degree of fractionation undergone by each sample and the radiogenic isotope ratios. This indicates that the samples are unlikely to have been contaminated with lithospheric material during fractionation and are therefore likely to record the isotopic signature of the mantle source.

Binary mixing curves between MORB-source mantle and the Imnaha mantle component, EM2, and HIMU mantle components are shown in Fig. 7, with 10% increments. A 53 Ma age corrected depleted Gorda MORB mantle (DGM) composition (Davis et al., 2008; Salters et al., 2011) was used as the end-member to represent a possible depleted mantle component.

The positive εHf and εNd values of all the Siletz samples reflect long term incompatible element depletion of the source. In Pb-Pb, Pb-Nd and Pb-Hf isotopic space the Siletz data lie mostly between the DGM–HIMU and DGM–EM2 mixing curves (Fig. 7), while a dominance of EM2 is required in the other diagrams. The Imnaha component also lies adjacent to the DGM-EM2 mixing curve, which indicates that the 3 samples with the significant EM2 component share similar isotopic signatures to the Columbia River Flood Basalts and the Snake River– Yellowstone hotspot trace (Wolff et al., 2008). Carlson (1984) also recognised that this component was present in the source of the CRFB, but proposed that it represented subducted sediment rather than a mantle plume component. However, the trace element data rule out the possibility of subducted sediment in the analysed samples as there is no ‘arc signature’ present in the mantle-normalised diagrams (i.e., no marked negative Nb-Ta anomaly or Th spike) (Fig. 5). Therefore, the Siletz isotope data are generally distinct from that of the Imnaha component and require a different source than the Columbia River Flood Basalts and Snake River– Yellowstone hotspot track.

There is a significant EM2 signature present in the majority of the data, indicating a deep mantle derived component (Weaver, 1991; Workman et al., 2004); this EM2. The most enriched samples display higher Dy/Yb values, indicating that they have undergone smaller amounts of melting of a deeper source relative to the more depleted samples (Fig. 6). The distinctive isotopic ratios of some of the samples may indicate that the lower melt fractions are preferentially sampling isotopically distinct enriched pods or plums (Cousens, 1996) or streaks within a mantle plume (Kerr et al., 1995).

As noted above the data also extend partially towards the HIMU component in Pb-Pb, Pb-Nd and Pb-Hf isotopic space (Fig. 7), indicating that this component may be present in the source of the Siletz basalts. The involvement of a HIMU like signature is consistent with observations in other NE Pacific locations. The samples which plot closest to HIMU have a more alkaline signature (highlighted in Fig. 7c) which suggests a relationship between the HIMU component and less extensive (deeper) melting, as is also observed at the Cobb Seamounts (Chadwick et al., 2014).

Some of the samples analysed from the Crescent Formation and Metchosin Igneous Complex lie adjacent to, or at the PREMA component composition in the Pb-Pb, Pb-Hf and Pb-Sr isotopic diagrams (Fig. 7c, d, e, f). However, while Siletz samples do appear to trend towards PREMA in some isotope diagrams (Pb-Pb), the majority of the Siletz terrane data, in most isotopic spaces, do not overlap with this component (Fig. 7). Although the Siletz terrane isotopic data demonstrates some similarity to the trends observed in, for example, the Hawaiian Islands, most notably in Hf–Nd isotopic space, the Siletz terrane samples have distinctly more radiogenic Pb isotopic ratios along with lower 87Sr/86Sr (Fig. 7). Multiple components are therefore most likely required to explain the trends observed in the isotopic data (Fig. 7) i.e., a depleted component isotopically similar to DGM, along with an enriched signature that straddles between EM2 and HIMU.

5.2. Mantle source T and P conditions

We have also attempted to determine the pressure and temperature of the mantle source region of the basalts of the Siletz terrane. Due to the relatively low magnesium contents and abundant clinopyroxene fractionation, few of the samples are suitable for the calculations. The methods outlined in Lee et al. (2009), which are based on Si and Mg contents of the magma resulted in only 7 samples producing mantle temperature and pressures, while the PRIMELTS 3 modelling software of Herzberg and Asimow (2015) produced results for 30 samples, 7 of which are in agreement with the temperature and pressure values obtained using the modelling of Lee et al. (2009). The results for these 7 samples were relatively consistent across the two techniques, producing a range of initial mantle potential temperatures (Tp) between ~1400 and 1500 °C and pressures between 1.5 GPa and 2.4 GPa (i.e., melting at less than 100 km). Mid-ocean ridge basalts display Tp ranges between ~1280 and 1400 °C, while intra-plate magmas are produced at Tp above >1400 °C (Herzberg et al., 2007; Lee et al., 2009). The Tp conditions estimated for the Siletz terrane data therefore support a mantle plume related origin for the Siletz terrane. The PRIMELTS 3 modelling results also require an olivine addition of on average 39% which indicates that subsequent to melting of the source, the magma is likely to have undergone significant fractionation in a crustal magma chamber to generate the observed compositions.

The results of this modelling, in combination with the isotopic data discussed in the previous section, indicate that the melts forming the Siletz terrane were derived from a heterogeneous and partially enriched mantle source, with an elevated temperature. The volume of melt produced (~ 2.6 × 106 km3; Trehu et al., (1994)) is consistent with extensive melting in a mantle plume head (Fitton and Godard, 2004; Herzberg and Gazel, 2009). While the range in calculated pressures may indicate inaccuracies in the models they may also point to melting over a significant range of depths which may represent the interaction of a plume and depleted mantle (Fig. 7). The amount of melting modelled in PRIMELTS3 is also indicative of a hotter-than-ambient mantle source with melting percentages of ~27% (with the majority of samples producing results between 25 and 33%), which are comparable to the extent of melting calculated for the Ontong Java Plateau and the Caribbean Plateau (~30% melting) (Fitton and Godard, 2004; Herzberg and Gazel, 2009).

5.3. Tectonic implications

There have been many tectonic models proposed to explain the origin of the Siletz terrane and its resulting geochemical and physical characteristics. One model which may explain the characteristics of the Siletz terrane is its formation in a marginal basin setting (Brandon, 2014; Wells et al., 1984), yet this does not accurately explain the proposed geometry of the plates at the time of eruption (Wells et al., 2014), the enriched isotopic signatures, nor the lack of arc signatures in the Siletz terrane (Fig. 5).

Another suggested model of formation for the Siletz terrane proposes subduction initiation along a transform fault, which can explain both the linear geometry of the terrane and its suggested volume (Stern, 2004). In this model a transform fault parallel to the North American margin between the Farallon and the partially subducted Kula/Resurrection plate represents a zone of weakness along which subduction can propagate (Babcock et al., 1992; Denny, 2012; Haeussler et al., 2003; Stern, 2004; Wells et al., 1998). The composition of the Siletz basalts does not support an early-arc tectonic setting resulting from subduction initiation as the primary mechanism of formation of the Siletz terrane. Isotopically, the Siletz basalts are much more enriched than those related to subduction initiation, e.g., in the basalts of the Izu Bonin and Mariana areas (Arculus et al., 2015; Stern, 2004).

An alternative slab window model (Babcock et al., 1992) is also not supported by the initial temperatures calculated for the mantle source of the Siletz terrane, as the temperature of primary magmas derived from slab windows are thought to be comparable to mid-ocean ridges (1300–1400 °C) (Hastie and Kerr, 2010), compared to 1400–1500 °C calculated for the Siletz terrane. The suggested characteristics of the mantle source of the Siletz terrane are also more enriched than the DMM-like mantle that is thought to produce slab window related magmatism (Hastie and Kerr, 2010; Herzberg et al., 2007). Finally, the trace elements of Siletz rocks do not generally display typical subduction related characteristics, such as Nb-Ta depletions or Th enrichments (Figs. 3 and 4) (Pearce, 2008).

The lack of uniform resemblance to the Imnaha component in the majority of Siletz samples appears to rule out the Yellowstone hotspot mantle plume head as a source region (Pyle et al., 2009). However, the composition of the plume source is unlikely to have been consistent over time, or may be heterogeneous, and therefore variation in the amount of melting for example may result in the sampling of parts of the plume with markedly different compositions (Escuder-Viruete et al., 2007; Kerr et al., 1995). In addition to this, several plate reconstructions propose that the Yellowstone hotspot was present in this area at the time of the formation of the Siletz terrane (Doubrovine et al., 2012; Müller et al., 1993; O'Neill et al., 2005).

Although volcanic ash layers/tuffs are typically unusual in oceanic plateaus (Kerr, 2014), there are abundant volcanoclastic horizons throughout the terrane, especially in the upper subaerial flows. Progressive shallowing throughout the sequence from deep to shallow water to subaerial environments indicates continual dynamic uplift by the plume (as well as accumulation of the lava pile) while also explaining the abundance of ash and volcanoclastic material. This has also been observed at other well-characterised mantle plume-related oceanic plateaus, such as in the Western Cordillera of Colombia, in the Caribbean Plateau (Kerr, 2014) and also in limited occurrences in the Ontong Java Plateau where tuff layers and sedimentary reworking occur (Mahoney et al., 2001).

The Siletz basalts are most depleted in the north and northeast, becoming more enriched in the central and southern areas and towards the west. The basalts are represented by a range of moderately enriched compositions (E-MORB-like) with even the most depleted samples being more enriched than typical N-MORB. The depth of melting also appears to increase towards the centre of the terrane. Sample 073RB is the most western sample and additionally has the largest HIMU component.

Overall, the mantle characteristics of the Siletz terrane best represent mixing of depleted mantle with an additional enriched input from a mantle plume. The temperature of the primary magmas has been calculated at 1400–1500 °C, which is hotter than ambient mid-ocean-ridge related mantle. The enriched component and mantle plume source has been proposed to represent the Yellowstone hotspot (Fig. 8) (cf. Duncan, 1982; McCrory and Wilson, 2013; Wells et al., 2014). The depleted mantle source component recorded in the relatively depleted samples from the NE (the Metchosin Igneous Complex and Crescent Formation) is comparable to mid-ocean ridge source mantle (more extensive melting), most likely the Farallon–Kula/Resurrection ridge, and the off axis interaction of an additional hotter enriched mantle source region (Fig. 8). Alternatively, the depleted component may be sampling a depleted relatively refractory portion of the mantle plume (Kempton et al., 2000; Kerr et al., 1995).

Fig. 8. Schematic diagram of the tectonic setting in which the rocks of the Siletz terrane were formed.

5.4. Youngest oceanic plateau?

Oceanic plateaus represent vast areas of over-thickened predominantly basaltic oceanic crust (>5 × 105 km3) the majority of which erupted over a few million years (Kerr, 2014; Ernst, 2014). Elevated topography and greater crustal thickness in comparison with ‘normal’ oceanic crust lead to an increase in buoyancy in oceanic plateaus. Therefore, notably for plateaus that collide with subduction zones shortly after formation (<5 Ma), the probability of partial accretion to the adjacent upper continental plate margin and so preservation within the geologic record is greatly increased (Cloos, 1993). Accreted oceanic plateaus have had a significant role in the growth of the continental crust; however, secondary processes such as deformation and burial lead to difficulties in the recognition of their characteristics once accreted (Kerr et al., 2000). Despite this, many examples of crustal growth through oceanic plateau accretion have been recognised in the geological record, including sections of the Wrangellia terrane (Greene et al., 2009a, 2009b), which along with the Siletz terrane highlights the role of oceanic plateaus in the growth and evolution of the western margin of North America.

Several well-characterised oceanic plateaus have similar heterogeneous trace element compositions to the Siletz terrane, e.g., the Cretaceous Caribbean plateau (Kerr and Mahoney, 2007) and the Late Triassic (ca. 225–231 Ma) Karmusten flood basalts of Vancouver Island and Alaska (Greene et al., 2009a, 2009b). In contrast, the largest oceanic plateau, the Cretaceous Ontong Java Plateau, is more depleted and anomalously homogeneous in composition (Fitton and Godard, 2004) (Figs. 4 and 6). While geochemistry alone cannot be used to distinguish oceanic plateaus in the geological record, their Nb/La ratio is often ~1 (Kerr, 2014), and this is generally reflected in the Siletz terrane data, which clusters between ratios of 1 and 1.5. Samples also show similar flat-to-LREE-enriched patterns with no significant Nb-Ta depletions or Th enrichments (Fig. 5).

Features of the terrane which support an oceanic plateau origin include; moderately enriched isotopic compositions, flat to slightly enriched REE patterns and the vast size and magmatic volume of the terrane (estimated to be 2.6 × 106 km3; Trehu et al., 1994), the bulk of which erupted within 5–7 Ma. In addition to this, mantle melt modelling for the Siletz terrane (as discussed above) indicates elevated mantle source potential temperatures, and similar amounts of partial melting and subsequent fractionation to the Caribbean and Ontong Java Plateaus. The terrane does, however, thin to 10 km at Vancouver Island and 6 km offshore, which is more typical of normal oceanic crust, compared to the average ~20 km thickness estimated for the remainder of the terrane (Hyndman et al., 1990; Parsons et al., 1999; Trehu et al., 1994). Sheeted dykes present at Bremerton (Crescent Formation) and on Vancouver Island (Metchosin Igneous Complex), which while unusual in other oceanic plateaus (Kerr, 2014), are consistent with plume–ridge interaction. Finally, as discussed above, there are relatively depleted signatures present in some sequences from the north of the terrane, comparable to mid-oceanic ridge environments. Overall, the accumulated evidence appears to support the conclusion that the Siletz terrane is an accreted oceanic plateau, and hence represents the youngest oceanic plateau thus far characterised.

  1. Conclusions
  1. The Siletz accreted mafic terrane has a large magmatic volume (2.6 × 106 km3) (Trehu et al., 1994) which represents a sequence of rocks that progressively shallow from marine to terrestrial environments and erupted over a relatively short time period (56–49 Ma).
  2. The rocks of the terrane are geochemically comparable, both in trace element (generally flat to LREE enriched REE patterns) and radiogenic isotope composition, to several well characterised oceanic plateaus.
  3. The estimated initial mantle potential temperatures of 1400–1500°C and the amount of partial melting undergone, ~ 25–33%, is similar to other mantle plume derived provinces such as the Ontong Java Plateau and the Caribbean plateau (Fitton and Godard, 2004; Herzberg and Gazel, 2009).
  4. The mafic rocks of the Siletz terrane appear to have been derived from a heterogeneous and partially enriched mantle source with an above ambient temperature. This source composition is comprised of a relatively depleted signature and EM2and HIMU-like enrichments. The enriched components may therefore represent melting of a heterogeneous mantle plume, possibly the Yellowstone Hotspot, which likely interacted with a mid-oceanic ridge, the Kula–Farallon (or Farallon–Resurrection) ridge.
  5. Although individually, the geochemical signatures and physical characteristics of the Siletz terrane can be interpreted differently, when taken together, the evidence for the Siletz terrane representing an accreted oceanic plateau linked to a mantle plume is compelling.

Click to open/close ReferencesReferences

Arculus, R. J., Ishizuka, O., Bogus, K. A., Gurnis, M., Hickey-Vargas, R., Aljahdali, M. H., Bandini-Maeder, A. N., Barth, A. P., Brandl, P. A., and Drab, L., 2015, A record of spontaneous subduction initiation in the Izu-Bonin-Mariana arc: Nature Geoscience, 8, 728-733.

Arndt, N. T., Kerr, A. C., and Tarney, J., 1997, Dynamic melting in plume heads: The formation of Gorgona komatiites and basalts: Earth and Planetary Science Letters, 146, 289-301.

Atwater, T., 1989, Plate tectonic history of the northeast Pacific and western North America: The eastern Pacific Ocean and Hawaii: Boulder, Colorado, Geological Society of America, Geology of North America, N, p.21-72.

Babcock, R. S., Burmester, R. F., Engebretson, D. C., Warnock, A., and Clark, K. P., 1992, A rifted margin origin for the crescent basalts and related rocks in the Northern Coast Range Volcanic Province, Washington and British-Columbia: Journal of Geophysical Research-Solid Earth, 97, 6799-6821.

Brandon, M. T., 2014, New evidence for Backarc Basin interpretation for Eocene Coast -Range Terrane in Proceedings 2014 GSA Annual Meeting in Vancouver, British Columbia 2014.

Breitsprecher, K., Thorkelson, D., Groome, W., and Dostal, J., 2003, Geochemical confirmation of the Kula-Farallon slab window beneath the Pacific Northwest in Eocene time: Geology, 31, 351-354.

Chadwick, J., Keller, R., Kamenov, G., Yogodzinski, G., and Lupton, J., 2014, The Cobb hot spot: HIMU-DMM mixing and melting controlled by a progressively thinning lithospheric lid: Geochemistry, Geophysics, Geosystems, 15, 3107-3122.

Cloos, M., 1993, Lithospheric buoyancy and collisional orogenesis: subduction of oceanic plateaus, continental margins, island arcs, spreading ridges, and seamounts: Geological Society of America Bulletin, 105, 715-737.

Cousens, B. L., 1996, Depleted and enriched upper mantle sources for basaltic rocks from diverse tectonic environments in the northeast Pacific Ocean: The generation of oceanic alkaline vs. tholeiitic basalts: Earth Processes: Basu, A., and S. Hart (Eds) Reading the Isotopic Code, AGU Monograph, 96, 207-231.

Davis, A., Clague, D., Cousens, B., Keaten, R., and Paduan, J., 2008, Geochemistry of basalt from the North Gorda segment of the Gorda Ridge: Evolution toward ultraslow spreading ridge lavas due to decreasing magma supply: Geochemistry, Geophysics, Geosystems, 9. DOI: 10.1029/2007GC001775

Denny, A. C., 2012, Geochemistry of the Lower Crescent Formation, Washington, and Re evaluation of Proposed Formation Hypotheses. Masters thesis, Carleton College, Northfield, Minnesota (2012)

Doubrovine, P. V., Steinberger, B., and Torsvik, T. H., 2012, Absolute plate motions in a reference frame defined by moving hot spots in the Pacific, Atlantic, and Indian oceans: Journal of Geophysical Research: Solid Earth, 117, B9, 2156-2202

Duncan, R. A., 1982, A captured island chain in the Coast Range of Oregon and Washington: Journal of Geophysical Research, 87, NB13, 827-837.

Engebretson, D. C., 1985, Relative motions between oceanic and continental plates in the Pacific basin, Geological society of America. Special Publication 206, p. 1-60. http://dx.doi.org/10.1130/SPE206-p1

Ernst, R.E., 2014, Large Igneous Provincs. Cambridge University Press, 653 p.

Escuder-Viruete, J., Perez-Estaun, A., Contreras, F., Joubert, M., Weis, D., Ullrich, T. D., and Spadea, P., 2007, Plume mantle source heterogeneity through time: Insights from the Duarte Complex, Hispaniola, northeastern Caribbean: Journal of Geophysical Research-Solid Earth, 112, B4. doi/10.1029/2006JB004323

Fitton, J. G., and Godard, M., 2004, Origin and evolution of magmas on the Ontong Java Plateau, in Fitton, J. G., Mahoney, J. J., Wallace, P. J., and Saunders, A. D., eds., Origin and evolution of the Ontong Java Plateau, Geological Society of London, Special Publication, 229, 151-178.

Fleming, S. W., and Trehu, A. M., 1999, Crustal structure beneath the central Oregon convergent margin from potential-field modeling: Evidence for a buried basement ridge in local contact with a seaward dipping backstop: Journal of Geophysical Research: Solid Earth, 104, B9, 20431-20447.

Greene, A. R., Scoates, J. S., Weis, D., and Israel, S., 2009a, Geochemistry of Triassic flood basalts from the Yukon (Canada) segment of the accreted Wrangellia oceanic plateau: Lithos, 110, 1-19.

Greene, A. R., Scoates, J. S., Weis, D., Nixon, G. T., and Kieffer, B., 2009b, Melting history and magmatic evolution of Basalts and picrites from the accreted Wrangellia Oceanic Plateau, Vancouver Island, Canada: J Petrol, 50, 467-505.

Groome, W. G., Thorkelson, D. J., Friedman, R. M., Mortensen, J. K., Massey, N. W., Marshall, D. D., and Layer, P. W., 2003, Magmatic and tectonic history of the Leech River Complex, Vancouver Island, British Columbia: Evidence for ridge-trench intersection and accretion of the Crescent Terrane: Special Paper 371 - Geological Society of America, 327-354.

Haeussler, P., Clark, J., and Kenneth, P., 2000, Geologic map of the Wildcat Lake 7.5'Quadrangle Kitsap and Mason counties, Washington. Geological Survey Open File Report 00-356.

Haeussler, P. J., Bradley, D. C., Wells, R. E., and Miller, M. L., 2003, Life and death of the Resurrection plate: Evidence for its existence and subduction in the northeastern Pacific in Paleocene-Eocene time: Geological Society of America Bulletin, 115, 867-880.

Hastie, A. R., and Kerr, A. C., 2010, Mantle plume or slab window?: Physical and geochemical constraints on the origin of the Caribbean oceanic plateau: Earth-Science Reviews, 98, 3-4, 283-293.

Hastie, A. R., Kerr, A. C., Mitchell, S. F., and Millar, I. L., 2008, Geochemistry and petrogenesis of Cretaceous oceanic plateau lavas in eastern Jamaica: Lithos, 101, 3-4, 323-343.

Hastie, A. R., Kerr, A. C., Pearce, J. A., and Mitchell, S., 2007, Classification of altered volcanic island arc rocks using immobile trace elements: development of the Th–Co discrimination diagram: Journal of Petrology, 48, 12, 2341-2357.

Hauff, F., Hoernle, K., van den Bogaard, P., Alvarado, G. E., and Garbe-Schonberg, C. D., 2000, Age and geochemistry of basaltic complexes in western Costa Rica: Contributions to the geotectonic evolution of Central America: Geochemistry Geophysics Geosystems, 1. doi/10.1029/1999GC000020

Herzberg, C., and Asimow, P., 2015, PRIMELT3 MEGA. XLSM software for primary magma calculation: Peridotite primary magma MgO contents from the liquidus to the solidus: Geochemistry, Geophysics, Geosystems, 16, 2, 563-578.

Herzberg, C., and Gazel, E., 2009, Petrological evidence for secular cooling in mantle plumes. Nature, 458, 619-622.

Herzberg, C., Asimow, P. D., Arndt, N., Niu, Y., Lesher, C. M., Fitton, J. G., Cheadle, M. J., and Saunders, A. D., 2007, Temperatures in ambient mantle and plumes: Constraints from basalts, picrites, and komatiites: Geochemistry, Geophysics, Geosystems, 8, 2006GC001390.

Hirsch, D. M., and Babcock, R. S., 2009, Spatially heterogeneous burial and high-P/T metamorphism in the Crescent Formation, Olympic Peninsula, Washington: American Mineralogist, 94, 8-9, 1103-1110.

Hyndman, R., Yorath, C., Clowes, R., and Davis, E., 1990, The northern Cascadia subduction zone at Vancouver Island: Seismic structure and tectonic history: Canadian Journal of Earth Sciences, 27, 3, 313-329.

Irving, E., 1979, Paleopoles and paleolatitudes of North America and speculations about displaced terrains: Canadian Journal of Earth Sciences, 16, 3, 669-694.

Kempton, P. D., Fitton, J. G., Saunders, A. D., Nowell, G. M., Taylor, R. N., Hardarson, B. S., and Pearson, G., 2000, The Iceland plume in space and time: a Sr-Nd-Pb-Hf study of the North Atlantic rifted margin: Earth and Planetary Science Letters, 177, 3-4, 255-271.

Kerr, A. C., 2014, Oceanic plateaus. In: Holland, H.C. & Turekian K. (series eds.). Treatise on Geochemistry 2nd Edition, Volume 4, Rudnick, R. (ed.) The Crust. Chapter 18, 631-667. Elsevier.

Kerr, A. C., and Mahoney, J. J., 2007, Oceanic plateaus: Problematic plumes, potential paradigms: Chemical Geology, 241, 3-4, 332-353.

Kerr, A. C., Tarney, J., Kempton, P. D., Spadea, P., Nivia, A., Marriner, G. F., and Duncan, R. A., 2002, Pervasive mantle plume head heterogeneity: Evidence from the late Cretaceous Caribbean-Colombian oceanic plateau: Journal of Geophysical Research-Solid Earth, 107, B7, DOI 10.1029,

Kerr, A., Marriner, G., Tarney, J., Nivia, A., Saunders, A., Thirlwall, M., and Sinton, C., 1997, Cretaceous Basaltic Terranes in western Columbia: elemental, chronological and Sr–Nd isotopic constraints on petrogenesis: Journal of Petrology, 38, 6, 677-702.

Kerr, A. C., Saunders, A. D., Tarney, J., Berry, N. H., and Hards, V. L., 1995, Depleted mantle plume geochemical signatures: No paradox for plume theories: Geology, 23, 843-846.

Lee, C. T. A., Luffi, P., Plank, T., Dalton, H., and Leeman, W. P., 2009, Constraints on the depths and temperatures of basaltic magma generation on Earth and other terrestrial planets using new thermobarometers for mafic magmas: Earth and Planetery Science Letters, 279, 1-2, 20-33.

Madsen, J. K., Thorkelson, D. J., Friedman, R. M., and Marshall, D. D., 2006, Cenozoic to Recent plate configurations in the Pacific Basin: Ridge subduction and slab window magmatism in western North America: Geosphere, 2, 1, 11-34.

Mahoney, J. J., Fitton, J. G., Wallace, P. J., 2001, Proceedings of the Ocean Drilling Program, Initial Reports, College Station, Texas, Ocean Drilling Program, 1-75.

Mahoney, J. J., Storey, M., Duncan, R. A., Spencer, K. J., and Pringle, M., 1993a, Geochemistry and age of the Ontong Java Plateau, in Pringle, M. S., ed., The Mesozoic Pacific: Geology, Tectonics and Volcanism, Volume Geophysical Monograph 77, American Geophysical Union, 233-261.

Mahoney, J. J., Storey, M., Duncan, R. A., Spencer, K. J., and Pringle, M., 1993b, Geochemistry and geochronology of Leg 130 basement lavas: nature and origin of the Ontong Java Plateau, in Berger, W. H., Kroenke, L. W., and Mayer, L. A., eds., Proceedings of the Ocean Drilling Program, Scientific Results, Volume 130: College Station, TX, Ocean Drilling Program, Texas A&M University, 3-22.

Massey, N. W. D., 1986, Metchosin Igneous Complex, southern Vancouver Island - ophiolite stratigraphy developed in an emergent island setting: Geology, 14, 7, 602-605.

McCrory, P. A., and Wilson, D. S., 2013, A kinematic model for the formation of the Siletz-Crescent forearc terrane by capture of coherent fragments of the Farallon and Resurrection plates: Tectonics, 32, 3, 718-736.

Müller, R. D., Royer, J.-Y., and Lawver, L. A., 1993, Revised plate motions relative to the hotspots from combined Atlantic and Indian Ocean hotspot tracks: Geology, 21, 3, 275-278.

Nobre Silva, I. G., Weis, D., Barling, J., and Scoates, J. S., 2009, Leaching systematics and matrix elimination for the determination of high-precision Pb isotope compositions of ocean island basalts: Geochemistry, Geophysics, Geosystems, 10, 8. doi: 10.1029/2010GC003176

O’Neill, C., Muller, D., and Steinberger, B., 2005, On the uncertainties in hot spot reconstructions and the significance of moving hot spot reference frames - art. no. Q04003: Geochem Geophys Geosyst, 6, 4003. doi: 10.1029/2004GC000784

Parsons, T., Wells, R. E., Fisher, M. A., Flueh, E., and ten Brink, U. S., 1999, Three-dimensional velocity structure of Siletzia and other accreted terranes in the Cascadia forearc of Washington: Journal of Geophysical Research-Solid Earth, 104, B8, 18015-18039.

Pearce, J. A., 1996, A user's guide to basalt discrimination diagrams, in Wyman, D. A., ed., Trace Element Geochemistry of Volcanic Rocks: Applications for Massive Sulphide Exploration, Volume Geological Association of Canada, Short Course Notes, 12, 79-113.

Pearce, J. A., 2008, Geochemical fingerprinting of oceanic basalts with applications to ophiolite classification and the search for Archean oceanic crust: Lithos, 100, 1, 14-48.

Pyle, D. G., Duncan, R. A., Wells, R. E., Graham, D. W., Harrison, B., and Hanan, B., Siletzia: An oceanic large igneous province in the Pacific Northwest, in Proceedings 2009 Portland GSA Annual Meeting 2009.

Salters, V. J., Mallick, S., Hart, S. R., Langmuir, C. E., and Stracke, A., 2011, Domains of depleted mantle: New evidence from hafnium and neodymium isotopes: Geochemistry, Geophysics, Geosystems, 12, 8. doi: 10.1029/2011GC003617

Seton, M., Müller, R. D., Zahirovic, S., Gaina, C., Torsvik, T., Shephard, G., Talsma, A., Gurnis, M., Turner, M., Maus, S., and Chandler, M., 2012, Global continental and ocean basin reconstructions since 200Ma: Earth-Science Reviews, 113, 3-4, 212-270.

Snavely, P. D., MacLeod, N. S., and Wagner, H. C., 1968, Tholeiitic and alkalic basalts of the Eocene Siletz River volcanics, Oregon coast range: American Journal of Science, 266, 6, 454-481.

Stern, R. J., 2004, Subduction initiation: spontaneous and induced: Earth and Planetary Science Letters, 226, 3-4, 275-292.

Sun, S.-s., and McDonough, W. F., 1989, Chemical and isotope systematics of oceanic basalts: implications for mantle composition and processes, in Saunders, A. D., and Norry, M. J., eds., Magmatism in the Ocean Basins, Volume 42, Geological Society of London, Special Publication, 313-345.

Trehu, A. M., Asudeh, I., Brocher, T. M., Luetgert, J. H., Mooney, W. D., Nabelek, J. L., and Nakamura, Y., 1994, Crustal architecture of the Cascadia fore-arc: Science, 266, 5183, 237-243.

Tejada, M. L. G., Mahoney, J. J., Castillo, P. R., Ingle, S. P., Sheth, H. C., and Weis, D., 2004, Pin-pricking the elephant: evidence on the origin of the Ontong Java Plateau from Pb-Sr-Hf-Nd isotopic characteristics of ODP Leg 192 basalts., in Fitton, J. G., Mahoney, J. J., Wallace, P. J., and Saunders, A. D., eds., Origin and evolution of the Ontong Java Plateau, Geological Society of London, Special Publication, 229, 133-150.

Tejada, M. L. G., Mahoney, J. J., Neal, C. R., Duncan, R. A., and Petterson, M. G., 2002, Basement geochemistry and geochronology of central Malaita, Solomon islands, with implications for the origin and evolution of the Ontong Java Plateau: Journal of Petrology, 43, 3, 449-484.

Thompson, P. M. E., Kempton, P. D., White, R. V., Saunders, A. D., Kerr, A. C., Tarney, J., and Pringle, M. S., 2004, Elemental, Hf-Nd isotopic and geochronological constraints on an island arc sequence associated with the Cretaceous Caribbean plateau: Bonaire, Dutch Antilles: Lithos, 74, 1-2, 91-116.

Weaver, B. L., 1991, The origin of ocean island basalt end-member compositions: trace element and isotopic constraints: Earth and Planetary Science Letters, 104, 381-397.

Weis, D., Frey, F. A., Schlich, R., Schaming, M., Montigny, R., Damasceno, D., Mattielli, N., Nicolaysen, K. E., and Scoates, J. S., 2002, Trace of the Kerguelen mantle plume: Evidence from seamounts between the Kerguelen Archipelago and Heard Island, Indian Ocean: Geochemistry Geophysics Geosystems, 3, 6. doi: 10.1029/2001GC000251

Wells, R., Bukry, D., Friedman, R., Pyle, D., Duncan, R., Haeussler, P., and Wooden, J., 2014, Geologic history of Siletzia, a large igneous province in the Oregon and Washington Coast Range: Correlation to the geomagnetic polarity time scale and implications for a long-lived Yellowstone hotspot: Geosphere, 10, 4, 692-719.

Wells, R., Engebretson, D., Snavely, P., and Coe, R., 1984, Cenozoic plate motions and the volcano-tectonic evolution of western Oregon and Washington: Tectonics, 3, 2, 275-294.

Wells, R. E., Jayko, A. S., Niem, A., Black, G. L., Wiley, T. J., Baldwin, E., Molenaar, K., Wheeler, K., DuRoss, C., and Givler, R., 2000, Geologic map and database of the Roseburg 30 x 60 quadrangle, Douglas and Coos counties, Oregon, US Department of the Interior, US Geological Survey.

Wells, R. E., Weaver, C. S., and Blakely, R. J., 1998, Fore-arc migration in Cascadia and its neotectonic significance: Geology, 26, 8, 759-762.

White, W. M., and Patchett, P. J., 1984, Hf-Nd-Sr isotopes and incompatible element abundances in island arcs: Implications for magma origins and crust-mantle evolution Earth and Planetery Science Letters., 67, 167-185.

Wolff, J., Ramos, F., Hart, G., Patterson, J., and Brandon, A., 2008, Columbia River flood basalts from a centralized crustal magmatic system: Nature Geoscience, 1, 3, 177-180.

Workman, R. K., Hart, S. R., Jackson, M., Regelous, M., Farley, K. A., Blusztajn, J., Kurz, M., and Staudigel, H., 2004, Recycled metasomatized lithosphere as the origin of the enriched mantle II (EM2) end-member: Evidence from the Samoan volcanic chain: Geochemistry Geophysics Geosystems, 5.

Yorath, C. J., Sutherland Brown, A., and Massey, N. W. D., 1999, Lithoprobe, southern Vancouver Island, British Columbia: geology; With contributions by TDJ England, RM Bustin, Ottawa. Geological Survey of Canada, Bulletin 498, 1999; 145 pages (3 sheets), https://doi.org/10.4095/210350