The Lunar Crater volcanic field (LCVF) in central Nevada (USA) is dominated by monogenetic mafic volcanoes spanning the late Miocene to Pleistocene. There are as many as 161 volcanoes (there is some uncertainty due to erosion and burial of older centers); the volumes of individual eruptions were typically ∼0.1 km3 and smaller. The volcanic field is underlain by a seismically slow asthenospheric domain that likely reflects compositional variability relative to surrounding material, such as relatively higher abundances of hydrous phases. Although we do not speculate about why the domain is in its current location, its presence likely explains the unusual location of the LCVF within the interior of the Basin and Range Province. Volcanism in the LCVF occurred in 4 magmatic episodes, based upon geochemistry and ages of 35 eruptive units: episode 1 between ca. 6 and 5 Ma, episode 2 from ca. 4.7 to 3 Ma, episode 3 between ca. 1.1 and 0.4 Ma, and episode 4, ca. 300 to 35 ka. Each successive episode shifted northward but partly overlapped the area of its predecessor. Compositions of the eruptive products include basalts, tephrites, basanites, and trachybasalts, with very minor volumes of trachyandesite and trachyte (episode 2 only). Geochemical and petrologic data indicate that magmas originated in asthenospheric mantle throughout the lifetime of the volcanic field, but that the products of the episodes were derived from unique source types and therefore reflect upper mantle compositional variability on spatial scales of tens of kilometers. All analyzed products of the volcanic field have characteristics consistent with small degrees of partial melting of ocean island basalt sources, with additional variability related to subduction-related enrichment processes in the mantle, including contributions from recycled ocean crust (HIMU source; high-µ, where µ = 238U/204Pb) and from hydrous fluids derived from subducted oceanic crust (enriched mantle, EM source). Geochemical evidence indicates subtle source heterogeneity at scales of hundreds of meters to kilometers within each episode-scale area of activity, and temporary ponding of magmas near the crust-mantle boundary. Episode 1 magmas may have assimilated Paleozoic carbonate rocks, but the other episodes had little if any chemical interaction with the crust. Thermodynamic modeling and the presence of amphibole support dissolved water contents to ∼5–7 wt% in some of the erupted magmas. The LCVF exhibits clustering in the form of overlapping and colocated monogenetic volcanoes that were separated by variable amounts of time to as much as several hundred thousand years, but without sustained crustal reservoirs between the episodes. The persistence of clusters through different episodes and their association with fault zones are consistent with shear-assisted mobilization of magmas ponded near the crust-mantle boundary, as crustal faults and underlying ductile deformation persist for hundreds of thousands of years or more (longer than individual episodes). Volcanoes were fed at depth by dikes that occur in en echelon sets and that preserve evidence of multiple pulses of magma. The dikes locally flared in the upper ∼10 m of the crust to form shallow conduits that fed eruptions. The most common volcanic landforms are scoria cones, agglomerate ramparts, and ‘a‘ā lava fields. Eruptive styles were dominantly Strombolian to Hawaiian; the latter produced tephra fallout blankets, along with effusive activity, although many lavas were likely clastogenic and associated with lava fountains. Eroded scoria cones reveal complex plumbing structures, including radial dikes that fed magma to bocas and lava flows on the cone flanks. Phreatomagmatic maar volcanoes compose a small percentage of the landform types. We are unable to identify any clear hydrologic or climatic drivers for the phreatomagmatic activity; this suggests that intrinsic factors such as magma flux played an important role. Eruptive styles and volumes appear to have been similar throughout the 6 m.y. history of the volcanic field and across all 4 magmatic episodes. The total volume and time-volume behavior of the LCVF cannot be precisely determined by surface observations due to erosion and burial by basin-fill sediments and subsequent eruptive products. However, previous estimates of a total volume of 100 km3 are likely too high by a factor of ∼5, suggesting an average long-term eruptive flux of ∼3–5 km3/m.y.

Intraplate volcanism has long been a topic of research focused on basalt petrogenesis and upper mantle composition and structure, due to relatively low degrees of crustal overprinting on geochemical signatures and, in many cases, abundant mantle xenoliths (e.g., Wörner et al., 1986; Farmer et al., 1989; Lum et al., 1989; Menzies, 1989; Luhr et al., 1995; Witt-Eickschen and Kramm, 1998; Gazel et al., 2012; Cousens et al., 2013; McGee and Smith, 2016). A growing body of literature in the past decade has also focused on physical processes of intraplate volcanoes, including magma ascent, relationships with crustal structure, and eruption (see reviews by Valentine and Gregg, 2008; Németh, 2010; Le Corvec et al., 2013a, 2013b; integrated case studies by Le Corvec et al., 2013c; Murcia et al., 2016). Recent work has been partly aimed at assessing the potential hazards associated with eruptions in monogenetic volcanic fields (e.g., Ort et al., 2008; Valentine and Perry, 2009; Houghton et al., 2006; Németh et al., 2012), but also has been driven by the increasing recognition of the complexity of these small volcanoes, which were traditionally thought to be simple and well understood (e.g., Houghton and Schmincke, 1986; Ort and Carrasco-Nuñez, 2009; Genareau et al., 2010; Martí et al., 2011; Bolós et al., 2012; van Otterloo et al., 2013; Pedrazzi et al., 2016).

Several fundamental questions arise in the study of intraplate, dominantly monogenetic, volcanic fields. (1) Are the mantle magma sources homogeneous or are they heterogeneous in time and space for a given field? (2) Are there relationships between magma sources, clusters, and preexisting crustal structures, and why? (3) What are the characteristics of the transition from crustal magma transport, mainly via dikes, to eruptions at the surface? (4) How do individual volcanoes grow and feed various eruptive processes? (5) Are there relationships between source characteristics and the scales and eruptive styles of individual volcanoes, and the time-volume behavior of the volcanic field? Here we report an integrated study of the late Miocene–Pleistocene Lunar Crater volcanic field (LCVF) in central Nevada (USA) that is aimed at exploring these general questions. We present new geologic maps of the field, build upon previous geochronological and petrogenetic studies including a reassessment of mantle sources through time and space, and summarize the time-space distribution of volcanism, relationships with crustal structures, hypabyssal plumbing of the volcanoes, volcano types, and eruptive styles. We interpret data related to these individual aspects of the volcanic field, and conclude by interpreting relationships between them where such can be defined. We also briefly summarize the evolution of the volcanic landscape over several million years, and then compare the LCVF with nearby volcanic fields.

Monogenetic volcanism is defined differently, so we first clarify our terminology. Consistent with Valentine and Connor (2015), we use the term monogenetic for volcanoes that have only one eruptive episode, after which they become extinct. An eruptive episode may span time scales of weeks to decades (see Wood, 1980; Luhr and Simkin, 1993). A monogenetic event encompasses any processes within a given eruptive episode. Our definition has no implications regarding the volume, number, or composition of magma batches that might feed a monogenetic event. Many eruptions, in the sense of individual explosions or periods of sustained activity of a given style (e.g., fountaining, effusion), typically occur within the context of a monogenetic eruptive episode; thus our definition does not imply a single eruptive process. We distinguish between volcanoes, also referred to as centers, and vents. An individual volcano can have multiple vents during a single monogenetic eruptive episode. An additional issue that sometimes arises in the literature is how to distinguish monogenetic volcanoes from polygenetic volcanoes that undergo multiple eruptive episodes separated by dormant periods (e.g., Németh, 2010). This distinction can be ambiguous because in a volcanic field there is always a chance that a new monogenetic volcano will form within the footprint of a previous one; as described herein, this is especially common where the volcanoes form dense, long-lived clusters. In our view, polygenetic volcanism typically requires a single sustained magma plumbing system in the crust that produces multiple eruptive episodes. Thus in the absence of petrologic evidence for sustained magma reservoirs shared by colocated centers of different ages, we assume that colocated volcanoes are simply overlapping monogenetic volcanoes. We understand that at some point the distinction between monogenetic and polygenetic volcanism is rather subjective, but this issue has no impact on the work presented here. We also distinguish eruptive episodes at individual volcanoes from magmatic episodes on the scale of the volcanic field. In the case of the LCVF we have identified four major magmatic episodes, each of which includes many monogenetic volcanoes. These large-scale episodes are each characterized by a specific range of volcano ages and by similar magma compositions.

The interior of the western United States is dotted with intraplate volcanic fields (Fig. 1) that individually contain between one and a few hundred Quaternary monogenetic volcanoes. Most of these fields are located near the edges of the actively extending Basin and Range Province (e.g., Smith and Luedke, 1984; Parsons, 1995), such as along the margin of the stable cratonic block that forms the Colorado Plateau, along the eastern front of the Sierra Nevada (California), or in the Rio Grande Rift (an eastern extension of Basin and Range).

The LCVF is relatively isolated from other major volcanic fields and is centrally located within the Basin and Range Province. Prior to the basaltic volcanism that is our focus, the area underwent voluminous ignimbrite-forming eruptions during the Oligocene to early Miocene, associated with the central Nevada caldera complex, and subsequent basaltic andesitic volcanism that continued through the mid-Miocene (Best et al., 2013, 2016). These older volcanic rocks are underlain by a thick sequence of Paleozoic carbonate and clastic sediments; the crystalline basement is Proterozoic in age (Menzies, 1989). Extensional tectonics, accompanied by formation of north- to north-northeast–trending Basin and Range topography, has been active since the Oligocene. Crustal thickness is ∼30–35 km, similar to most of the Basin and Range Province, but heat flow in the area of the LCVF is relatively low (<62 mW/m2, the so-called Eureka low; Thompson and Burke, 1974; Sass et al., 1994; Parsons, 1995). The lithosphere beneath the LCVF is estimated to be 60–75 km thick (Li et al., 2007; Levander and Miller, 2012). High-resolution P and S wave tomography of the western United States by the transportable array component of the EarthScope USArray (www.usarray.org/) provides an image of the upper mantle in the region (Schmandt and Lin, 2014). A profile through the youngest part of the LCVF (Fig. 2) indicates that the field is underlain by a domain of relatively slow upper mantle, suggesting the presence of some melt due to upwelling and/or volatile enrichment; the relatively low heat flow at the LCVF leads us to favor the latter. The slowest part of this domain correlates well with estimated melting depths for Pleistocene LCVF volcanoes (Cortés et al., 2015). Although we do not speculate about why the domain is in its current location, its presence likely explains the unusual location of the volcanic field within the interior of the Basin and Range Province.

Most of the LCVF is expressed on the surface as an ∼90-km-long, ∼5–15-km-wide belt of volcanoes exposed on and around two mountain ranges (Fig. 3, but see following discussion about the field boundaries). The older late Miocene–Pliocene part of the field is mainly in the Reveille Range and adjacent Kawich Valley (Fig. 4), and the younger Pliocene–Pleistocene part extends northeastward into the Pancake Range (Fig. 5). Note that the name Pancake Range is used for a belt of ignimbrite mesas and uplifted structural blocks, but in detail it includes low-lying areas such as the Lunar Lake basin, where many of the Pleistocene volcanoes are situated. Basalts (loosely referring to trachybasalts, tephrites, basanites, and basalts sensu stricto) and very minor volumes of relatively evolved compositions erupted beginning ca. 6 Ma (Foland and Bergman, 1992; Yogodzinski et al., 1996) and continued to the Marcath eruption, ca. 35 ka (Table 1), indicating that the LCVF is essentially an active field that is in a period of dormancy between monogenetic events. The field contains 142–161 monogenetic volcanoes (46–65 in the Reveille Range with uncertainty related to extensive surface modification and erosion there, and 96 in the Pancake Range), most of which consist of a pyroclastic construct (cone or rampart) and an attendant lava field. Three small centers in the Reveille Range are trachyandesite to trachyte domes and lavas (Yogodzinski et al., 1996), but most are highly eroded scoria cones and mesa-capping lava fields (e.g., Hintz and Valentine, 2012; Harp and Valentine, 2015). The Pancake Range volcanoes also include four maars, two of which are associated with pristine to partly eroded scoria cones or elongate agglomerate ramparts from the same monogenetic events (Valentine et al., 2011; Amin, 2013; Valentine and Cortés, 2013), and the poorly exposed remnant of one tuff cone.

Geochronology

The ages of lavas in the LCVF have been compiled in four studies. The first was by Kargel (1987); these K-Ar results were summarized along with isotopic data by Foland and Bergman (1992). Naumann et al. (1991) reported K-Ar ages for the northern part of the Reveille Range and, along with Yogodzinski et al. (1996), suggested that volcanism there occurred in two episodes (before 5 Ma, and ca. 4.7–3 Ma) that were also chemically and petrographically distinct. A third set of ages was presented by Dickson (1997) and by Stickney (2004), and the fourth set is our work. Previous determinations by Kargel (1987) are not included in our discussion because they were whole-rock K-Ar analyses that tended to return older ages (especially for the Pleistocene rocks) that are not consistent with the more recent work, although Kargel’s dates appear consistent with most stratigraphic relationships in terms of relative ages. For example, while Kargel (1987) correctly identified the Marcath eruption (described in the following) as the youngest in the volcanic field, his age is ca. 350 ka, which is inconsistent with the ca. 38–35 ka ages determined by independent methods (surface exposure and 40Ar/39Ar) and with the pristine morphology of the lavas and cone. The older ages of Kargel (1987) might be related to inclusions and megacrysts that are common in the lavas. The Stickney (2004) ages, which focus on the northern approximately one-third of the Pancake Range, are not included because we found inconsistencies between the reported values and stratigraphic relationships in the field.

Table 1 summarizes our ages along with those of Dickson (1997) and Naumann et al. (1991), a total of 35 eruptive units (lavas, feeder dikes, and plugs) as well as 3 duplicate analyses by the different authors. All but one of these determinations is based on groundmass plagioclase, using either K-Ar analyses for the late Miocene–Pliocene rocks, or 40Ar/39Ar analyses for the younger rocks. The exception is Lunar Crater maar, where limited ejecta exposures and phreatomagmatic eruption mechanisms resulted in very little juvenile material that could be used for 40Ar/39Ar analyses. Here we applied the cosmogenic 10Be surface exposure dating method, analyzing quartz phenocrysts from the Miocene rhyolitic ignimbrite that was ejected onto the tephra ring as ballistic blocks and exposed in the crater wall after having been buried by 3.82 Ma basalts (sample GVLC13–50; Table 1). The age of Lunar Crater maar is bracketed between ca. 72 ka (the minimum age from 10Be exposure dating [Supplemental Table1]) and 190 ka (the age of the youngest lava that the Lunar Crater maar ejecta appears to overlap; sample LC12–85; Table 1). This range of ages, especially toward the younger portion of the range, is consistent with a qualitative estimate based upon morphology and soil development (Valentine et al., 2011).

Volcanism in the LCVF has gradually shifted northward since initial activity in the Reveille Range, a direction that is oblique to the North American plate motion and therefore not a hotspot track. During any given time window volcanoes formed apparently randomly (in time) over a significant fraction of the total area (Fig. 6; Naumann et al., 1991; Foland and Bergman, 1992). We have identified four episodes of magmatism, two of which were originally discussed by Naumann et al. (1991) and Yogodzinski et al. (1996), based upon age, geochemistry, and petrographic similarities. Episode 1 (ca. 6–5.1 Ma) appears to have been restricted to the Reveille Range, probably extending 20–25 km farther south than indicated in Figure 6, where similar basaltic rocks occur but have not been dated (see Fig. 4). Where episode 1 is in contact with episode 2 rocks the contact is an unconformity, and episode 1 clasts are included in isolated fanglomerates around paleohighs that hosted some of the episode 2 centers, consistent with a significant time break between the episodes. Episode 2 (ca. 4.7–3 Ma) is prevalent in the northern Reveille Range and southern Pancake Range, probably extending ∼10 km farther north than indicated in Figure 6; small, faulted remnants of a vent structure and lava flows occur around the edge of the Lunar Lake basin (a half-graben within the Pancake Range; Fig. 5) that are petrographically and morphologically similar to episode 2 lavas (e.g., Broken Cone; Hintz and Valentine, 2012) but these were not analyzed. It is possible that additional episode 2 lavas are subsided and buried in Lunar Lake basin, and that these might fill the apparent ∼1.9 m.y. gap between episodes 2 and 3 (Fig. 6). Episodes 3 and 4 (1.1–0.4 Ma and 0.3–0.035 Ma, respectively) overlap each other spatially, as well as the northern extent of episode 2, in the Pancake Range.

Spatial Distribution of Volcanoes and Structural Relationships

Defining Boundaries

Definition of the boundaries of the LCVF is difficult due to subsidence of basins on either side of the Reveille and Pancake Ranges. It is likely that some basaltic centers occur in these basins but are buried by alluvial fill, and that the volcanic field is not as narrow as it appears from a surface geologic map, especially in its southern part. The presence of a small number of basaltic centers in the Quinn Canyon and Grant Ranges on the east side of Railroad Valley (Fig. 3) that have degrees of erosion and geochemistry similar to those of episode 1 and 2 centers in the Reveille Range (Emery, 2012) suggests that volcanism may have been continuous across that basin during late Miocene–Pliocene time. A line of three small, eroded, mostly buried cones occurs in the Kawich Valley to the west of the southern Reveille Range, along with a partly buried lava field (Fig. 4). These features suggest the possibility of additional buried volcanoes in that basin.

In the Pancake Range (Fig. 5) the eastern boundary of the volcanic field can be reasonably well defined and corresponds to the eastern, normally faulted margin of Lunar Lake basin (an escarpment referred to as the Wall; Ekren et al., 1972). The Big Sand Springs Valley that borders the northwestern part of the field might contain some buried basaltic centers. Drill core from a well located in Big Sand Springs Valley, ∼5 km northwest of the Marcath cone, intersected basaltic lava intercalated with alluvium at a depth of 43–53 m (Ekren et al., 1974). The northern boundary of the volcanic field is abrupt, with a stepped transition from dense clusters of cones to none at all.

Pancake Range Clustering and Alignments

Tadini et al. (2014) conducted a formal analysis of vent clustering and alignment for the Pancake Range, and the age data provide additional insight. Most volcanoes in the Pancake Range occur in clusters, while ∼20% of the volcanoes are isolated (>∼1 km from nearest neighbor). Some clusters are quite dense; one at the northern end of the field contains 12 separate cones within a 5 km2 area. Scoria cones commonly overlap within such clusters because their individual basal diameters are ∼400–1500 m. Two dated volcanoes within ∼500 m of each other in the northernmost cluster have ages of ca. 130 ka and 740–620 ka (Giggle Springs and Mizpah lavas, respectively; Table 1; Cortés et al., 2015; Rasoazanamparany et al., 2015). The youngest volcano in the field, Marcath, is on top of older cones, and the flow paths of its lavas were strongly influenced by cone-cluster topography. In the cluster that includes Easy Chair volcano, the field evidence suggests approximately colocated vents for eruptions ca. 840 ka (proto–Easy Chair; Table 1) and 140 ka. Thus the clusters repeatedly underwent monogenetic eruptions over periods of many hundreds of thousands of years. We do not interpret this to be evidence of a polygenetic volcano in the sense of having a sustained magma feeding and crustal reservoir system over tens or hundreds of thousands of years, because in both the cases specified, the two eruptions (Mizpah-Giggle Springs and proto–Easy Chair-Easy Chair) have significantly different source geochemistry (Rasoazanamparany et al., 2015). Rather, these are random intersections of event locations within clusters. Similar long-term reuse of clusters was documented in the San Francisco volcanic field (Arizona; Conway et al., 1997). It is likely that there are many other such situations in the LCVF that were not dated, or where later events completely buried the products of earlier ones. We also note that the ca. 140 ka Easy Chair eruption was quite complex, producing a series of agglomerate ramparts, two central scoria cones, and a maar crater, all apparently within the context of a single monogenetic event (Valentine and Cortés, 2013). The term monogenetic volcanism does not imply simple, isolated volcanoes.

The volcanic field as a whole is elongated ∼N30°E. Clusters within the Pancake Range are also preferentially elongated in this direction, as are many craters and nearly all eruptive fissures and ramparts (Tadini et al., 2014). Thus ∼N30°E is a preferred orientation for processes at a range of length scales. Tadini et al. (2014) summarized broad relationships between the orientations of volcanism and structures in the host rocks (mainly Oligocene–early Miocene volcanics) and Quaternary tectonics. The current maximum horizontal principal stress in the region, based mainly upon focal mechanisms of earthquakes on normal faults, is also oriented ∼N30°E (Heidbach et al., 2008). Quaternary faults in the region surrounding the LCVF have a similar, but broader, range of strikes; most are between north-south and N60°E. Older faults trend closer to north-south or north-northwest–south-southeast, reflecting earlier stress fields within the Basin and Range (older than ca. 10 Ma; Zoback et al., 1981), or have orientations more closely related to Oligocene–Miocene caldera structures.

Relationships of Preexisting Structures, Vents, and Feeder Dikes in the Pancake Range

Although burial by young volcanic products has obscured relationships between faults and vents, there are some locales in the Pancake Range where it is clear that vents are colocated with preexisting faults. These are mainly around the periphery of the field, where volcanoes are relatively isolated, so it is possible that the vent-fault relationships there do not hold in the more concentrated parts of the volcanic field (see Tadini et al., 2014). Near the southwestern (Fig. 7A) and southern edges of the Pancake Range scattered vents are located on normal faults associated with the hinge side of the Lunar Lake half-graben and with the southern escarpment of the Citadel Mountain structural block. Four vents are prominently located along the escarpment bounding the east side of Lunar Lake basin (the Wall; Figs. 7B, 7C). In some of these cases the vent-hosting faults have been active since the eruptions, and in others it is not clear (the latter are morphologically younger volcanoes, likely associated with episode 3). Previously workers have inferred the presence of buried faults from elongations and alignments of cones and ramparts (e.g., Scott and Trask, 1969; Snyder et al., 1972), but we do not make this inference without direct evidence for displacement before and/or after volcanic activity at a given location.

For a feeder dike to occupy a fault plane, the magma pressure in the dike must be sufficient to overcome the minimum compressive stress, σ3, plus the added components of stress due to the strike and dip of the fault plane (compared to being a vertical plane perpendicular to σ3 in the regional extensional setting; Delaney et al., 1986; Valentine and Krogh, 2006; Gaffney et al., 2007). In general, faults occupying vents in the Pancake Range, and by inference their feeder dikes, occur where the host faults have the north-northeast trend that is shared by vent alignments and the volcanic field elongation, rather than where faults have orientations that differ from the Quaternary principal stress orientations. The preference of feeder dikes for faults striking perpendicular to the Quaternary σ3 indicates that magmas did not have large overpressures. It is likely that the ascending dikes were captured by faults in the upper few hundred meters of the crust, where the component of overburden pressure acting on the fault planes was low enough for the low magma pressures to overcome it (see Valentine and Krogh, 2006; Gaffney et al., 2007). In this sense the faults did not represent pathways or weaknesses for magmas to follow from depth, as is often inferred when volcanoes are located on faults, but rather represent coincidences where a feeder dike was sourced beneath a favorably oriented fault (e.g., Valentine and Perry, 2006).

The northernmost part of the field (Fig. 7D) appears to have concentrated basaltic volcanism in a small graben-like domain, bounded by approximately north-south–trending blocks of Oligocene–early Miocene volcanic rocks, that is now partly filled with basalt. This ∼3-km-wide domain extends ∼5 km north-northeast beyond the limit of basaltic volcanism, forming a complex belt of normally faulted blocks of Oligocene–early Miocene ignimbrites (see northern edge of geologic map; Fig. 5). To the south-southeast the domain apparently plunges toward the larger half-graben of the Lunar Lake basin and disappears beneath basaltic lavas and sedimentary fill. Most of the volcanoes in the Lunar Lake basin occur along the south-southeast projection of the graben-like domain; however, some volcanoes might be subsided and buried beneath the basin, leaving us with only a partial picture of volcano distribution. The western horst of the graben-like domain in turn formed an ∼100-m-high escarpment bounding the edge of Big Sand Springs Valley. Basaltic volcanism has been intense along this escarpment and most of the episode 4 products have vented along it (e.g., Giggle Springs and Marcath volcanoes) or its projection beneath the Lunar Lake basin (e.g., Easy Chair volcano and Lunar Crater maar). Small normal faults have displaced episode 3 lavas to the west of the escarpment, indicating concurrent extensional deformation and volcanism.

Reveille Range

Although formal analysis of clustering and alignments has not been conducted for the Reveille Range volcanoes, it is clear that they share some common elements with those in the Pancake Range, especially the tendency for centers to occur in clusters (Fig. 4). Individual clusters in the northeastern part of the range underwent both episode 1 and 2 volcanism, with age differences of ∼1–2 m.y. (e.g., compare age of sample R8–1-17-LN, ca. 5.13 Ma, with that of nearby but stratigraphically higher sample R8–1-18-LN, ca. 3.99 Ma; Table 1). Vents, in the forms of highly eroded cone and agglomerate rampart remnants, are concentrated on either side of the range and some of these are clearly related to the range-bounding, north-northwest–trending fault along the western side (Fig. 4). The line of isolated vents in the Kawich Valley to the west has a similar orientation and might be associated with a buried, inactive fault. Other volcanoes might be colocated with preexisting structures; for example, Dark Peak, a cone remnant in the north-central part and along the crest of the Reveille Range, is centered within 200 m of an older caldera rim (Harp and Valentine, 2015).

Petrology and Geochemistry—Magma Sources and Ascent

Data Sources

Representative major and trace element analyses for the four eruptive episodes were compiled from Yogodzinski et al. (1996), Dickson (1997), Cortés et al. (2015), and Rasoazanamparany et al. (2015). Previously published Sr, Nd, and Pb isotope data for samples from episode 1 (Yogodzinski et al., 1996) and episodes 3 and 4 (Rasoazanamparany et al., 2015) were also compiled for comparison in this study. In addition to these previously published data, we performed new Sr, Nd, and Pb isotope analyses of subset of episode 2 basalt samples that were originally reported by Yogodzinski et al. (1996). Sample leaching, dissolution, chemical purifications, and mass spectrometry were performed following the procedures used previously for the episode 3 and 4 sample isotopic analyses (detailed in Rasoazanamparany et al., 2015), in order for these data sets to be directly comparable. Relatively minor differences between the new data and the previous results of Yogodzinski et al. (1996), likely due to differences of acid leached versus nonleached samples, respectively, indicate that the isotopic data sets for all four magmatic episodes can be meaningfully compared to one another.

Summary of Petrologic and Geochemical Data

Episodes 1 and 2. Episodes 1 and 2 were first described in studies of the Reveille Range by Naumann et al. (1991). Episode 1 (ca. 5.9–5.1 Ma; Table 1) products comprise primarily alkalic, hawaiitic basalts containing olivine and plagioclase phenocrysts and megacrysts of laboradorite, with whole-rock SiO2 and MgO contents ranging, respectively, from 45.5 to 49.4 wt% and 3.4 to 6.1 wt% (Yogodzinski et al., 1996). Basalts of episode 2 (ca. 4.6–3.0 Ma; Table 1) are also alkalic, but range from trachybasalt to basanite, and are characterized by megacrysts of plagioclase, clinopyroxene, and/or amphibole and phenocrysts of olivine, plagioclase, clinopyroxene, and Fe-Ti oxides. Most are nepheline-normative, and SiO2 and MgO both vary widely, from 41.8 to 47.8 wt% and 4.4 to 10.2 wt%, respectively (Figs. 8 and 9; Yogodzinski et al., 1996). Very small quantities of more evolved trachyandesite and trachytes erupted at the beginning of or just prior to episode 2. Gabbro inclusions as much as several centimeters in size also occur in many of these lavas; some localities contain abundant dunite and harzburgite inclusions (Yogodzinski et al., 1996).

Episodes 1 and 2, like the subsequent episodes 3 and 4, have trace element abundances that are similar to ocean island basalts (OIBs), with enrichments relative to primitive mantle in moderately to highly incompatible elements, including relative enrichments in the high field strength elements Nb and Ta, and in light rare earth elements (REEs) relative to heavy REEs (Fig. 10). These data are consistent with an asthenospheric, OIB-like mantle source region for both magmatic episodes. Episode 1 and 2 samples overlap in Zr/Nb (most in a narrow range from 6 to 8, with a few as high as 12 to 14), suggesting that they may share a geochemically similar source; in contrast, a large range in La/Yb ratios (Fig. 11) among samples from episode 1 (12–18) and especially episode 2 (10–25) could be indicative of more variable degrees of source melting (e.g., Cortés et al., 2015).

The Sr, Nd, and Pb isotope signatures of the episode 1 and 2 basalts likewise exhibit significant ranges. Most pronounced is the range in 87Sr/86Sr in episode 1 basalts (0.70427–0.70611) that correlates with a (more modest) variation in Nd isotopes (εnd +1 to +5; 143Nd/144Nd from 0.512680 to 0.512866), including basalts in the Reveille and Quinn Canyon Ranges (Yogodzinski et al., 1996; Emery, 2012). These isotope systematics, which result in a concave downward trend of the episode 1 basalts in Sr-Nd isotope space (Fig. 12), led Yogodzinski et al. (1996) to infer contamination by the thick carbonate section that is beneath much of the volcanic field. Foland and Bergman (1992) showed that Sr and Nd isotopes of basalts in the area are positively correlated with oxygen isotopes (δ18O) and also suggested that older Pliocene basalts in the Reveille Range were contaminated by continental crust.

New whole-rock Sr, Nd, and Pb isotope ratios for episode 2 basalts from the Reveille Range (Table 2) show a substantially more limited variation in 87Sr/86Sr (0.70337–0.70351), as well as 143Nd/144Nd (0.512810–0.512855) and Pb (206Pb/204Pb = 19.122–19.375; 207Pb/204Pb = 15.591–15.615; 208Pb/204Pb = 38.770–38.966) isotopic compositions. In comparison to the episode 1 basalts, episode 2 basalts are significantly less radiogenic in Sr, and comparable to or higher in 143Nd/144Nd than the most radiogenic of the episode 1 basalts (Fig. 12A). Episode 2 Pb isotope signatures plot on the same trend as episode 1 in 206Pb/204Pb versus 208Pb/204Pb, but most have higher 206Pb/204Pb (Fig. 13A). The 207Pb/204Pb variations are minimal, but on average, 207Pb/204Pb is slightly higher in episode 1 basalts (Fig. 13B); this has been attributed to the influence of crustal assimilation (Yogodzinski et al., 1996). The isotope data, as observed for the trace element data, are consistent with the two magmatic episodes sharing a similar but heterogeneous mantle source, the composition of which is more closely reflected by the episode 2 basalts that, unlike the episode 1 basalts, did not undergo significant crustal assimilation (Yogodzinski et al., 1996, and this study). We did not conduct new analyses on the small volumes of trachyandesite and trachyte that occur in the northeastern part of the Reveille Range. Their ages (Table 1) correspond to early episode 2. Yogodzinski et al. (1996) assessed these rocks and determined, based mainly upon isotopic evidence, that they were probably evolved from episode 1 basaltic magmas.

Episodes 3 and 4. Episodes 3 and 4 occur only in the Pancake Range. Episode 3 lavas (ca. 1.1–0.4 Ma; Table 1; referred to as the HiMUp group in Rasoazanamparany et al., 2015) (HIMU—high-µ, where µ = 238U/204Pb) range in composition from basalt to trachybasalts (Fig. 8), and are characterized by variable concentrations of plagioclase megacrysts, relatively sparse pyroxene megacrysts, and common gabbro inclusions, with variable proportions of phenocrysts of olivine, plagioclase, and pyroxene. Whole-rock SiO2 and MgO contents range, respectively, from 44.1 to 48.3 wt% and 7.0 to 10.7 wt%; some samples are close to primitive mantle melt compositions with Ni contents >120 ppm (and as much as 260 ppm) and Cr contents >220 ppm (and as much as 270 ppm) (Rasoazanamparany et al., 2015). Episode 4 (ca. 300 ka; Table 1; referred to as the GEM group in Rasoazanamparany et al., 2015) includes tephrites, basanites, and trachybasalts. They contain conspicuous euhedral megacrysts of plagioclase, pyroxene, olivine, and amphibole in varying proportions, all of which are interpreted as comagmatic with the host lavas, as well as a range of ultramafic inclusions (Bergman et al., 1981; Bergman, 1982; Foland and Bergman, 1992). Phenocrysts are olivine, plagioclase, and pyroxene, also in varying proportions. Whole-rock SiO2 and MgO contents are similar to those of episode 3, ranging, respectively, from 43.5 to 45.8 wt% and 7.2 to 10.7 wt%.; likewise, some episode 4 samples also approximate primitive mantle melt compositions with Ni contents reaching 360 ppm and Cr contents reaching 820 ppm (Rasoazanamparany et al., 2015).

Detailed whole-rock and mineral chemistry data from two episode 3 eruptive centers (Cortés et al., 2015) suggest that the magmas may have originated by 3%–5% partial melting of an olivine websterite-lherzolite source at depths of 60–130 km. The magmas partially crystallized at depths of 30–45 km, either while stalled near the base of the crust (30–35 km thick) or during ascent in that depth interval, after having risen directly from source depths (Spera and Fowler, 2009; Rasoazanamparany et al., 2015; Cortés et al., 2015). After ∼50% crystallization, the magmas then ascended rapidly to the surface, with no geochemical evidence for crustal contamination. In contrast, episode 4 magmas appear to have crystallized at two levels: first near the base of the crust (30–45 km depth), and second, at mid-crustal levels where, in the case of Marcath, amphibole megacrysts formed and indicate a dissolved water content of 5–7 wt% in magma at that depth (Cortés et al., 2015). Rapid ascent between (and possibly including) crystallization levels is supported by the relatively primitive nature of these magmas and the lack of detectable assimilation of crustal components (Cortés et al., 2015; Rasoazanamparany et al., 2015).

Trace element patterns of episode 3 and 4 lavas are similar to those of episodes 1 and 2, with OIB-like relative enrichments in Nb-Ta, and light REEs relative to heavy REEs. All sample groups largely overlap in their trace element abundance patterns (Fig. 10). However, in detail the episode 3 and 4 lavas differ from one another modestly in their trace element patterns, with episode 4 lavas having on average higher concentrations of highly incompatible trace elements, steeper REE patterns (higher La/Yb), and higher Zr/Hf (Rasoazanamparany et al., 2015), but slightly lower Zr/Nb ratios (Fig. 11) than episode 3 basalts.

Radiogenic isotope data are available for 34 lavas from episodes 3 and 4; a subset of those (19 and 13, respectively) was analyzed for Os and Hf isotopes (Rasoazanamparany et al., 2015). Episode 4 data include 20 samples from Marcath, Giggle Springs, and Easy Chair; episode 3 data comprise 14 samples from Hi Desert, Mizpah, proto–Easy Chair, and 8 unnamed cones. The radiogenic isotope data are consistent with the trace element data in that they likewise reflect OIB-like mantle source signatures. In Sr-Pb isotope space (Fig. 14), the lavas collectively produce a two-component mixing trend between a high 206Pb/204Pb, low 87Sr/86Sr (HIMU type) end member and a low 206Pb/204Pb, high 87Sr/86Sr (EM type) end member that compose the episode 3 and 4 samples, respectively. Based on the combined trace element and radiogenic isotope data, these end members have been interpreted to relate to OIB-type mantle sources, including an asthenospheric HIMU-type source containing a component of recycled, dehydrated oceanic crust, and an EM-type source that likely resulted from introduction of hydrous fluids from subducted oceanic crust that further metasomatized the HIMU-type mantle (Rasoazanamparany et al., 2015).

Synthesis of Petrology and Geochemistry in Magmatic Episodes 1–4

Comparing and contrasting the four magmatic episodes provide constraints on variations through time in magma sources and processes of melt generation, melt-lithosphere interaction, and magma fractionation. Among the four magmatic episodes, MgO contents exhibit a complex temporal relationship, but on average MgO contents increase with time (Fig. 9). Crustal assimilation is only significant in episode 1 (trachybasalts) and in the 4.3 Ma trachytes, as indicated most clearly by variable and highly radiogenic Sr isotope signatures; in contrast, basalts from episode 2, as well as those from episodes 3 and 4, show no evidence for significant crustal contamination (Yogodzinski et al., 1996; Rasoazanamparany et al., 2015). The apparent lack of crustal contamination in episode 2 trachybasalts suggests that episode 2 magmas were able to ascend from source to surface and undergo fractionation with minimal interaction with the crustal basement. The near primitive compositions of many episode 3 and 4 basanites and tephrites, some of which are mantle xenolith bearing, also suggest rapid source to surface transport.

In general, the magmas also become more silica undersaturated with time (Fig. 8), with episode 1 comprising solely basalts and trachybasalts and episode 4 solely basanites (episodes 2 and 3 contain both), consistent with an overall increase in depth of melting and decrease in degree of melting with time. Higher concentrations of incompatible trace elements at a given MgO content for episode 2 basanites relative to trachybasalts (Yogodzinski et al., 1996) are consistent with a lower degree of melting to produce the basanites. Trace element patterns are consistent with OIB-type sources for the magmas from all four episodes (Fig. 10), although subtle compositional differences between episodes are observed. In a Zr/Nb versus La/Yb plot (Fig. 11), episodes 1 and 2 are characterized by high Zr/Nb relative to episodes 3 and 4, and episodes 3 and 4 are largely distinct from one another in both La/Yb and Zr/Nb. Variations in the La/Yb ratio are indicative of variable degrees of melting of a garnet-bearing mantle source, whereas variations in Zr/Nb most likely reflect a heterogeneous source composition and indicate subtle differences in chemistry between the mantle sources tapped during magmatic episodes 1 and 2 compared to that of episode 3, which in turn differs from that of episode 4. Note that the Cortés et al. (2015) model was based solely on major and trace element data, and assumed a single, relatively dry source, thus inferring that variations in La/Yb are due to variable and relatively higher degrees of partial melting (a few percent). In contrast, in Rasoazanamparany et al. (2015), the isotopic variations were interpreted to require distinct and strongly metasomatized mantle sources, such that highly incompatible trace element ratios are influenced by source heterogeneity, resulting in substantially lower (<1%) degrees of melting.

Distinctions between mantle sources tapped during the respective magmatic episodes are further revealed by the highly variable radiogenic isotope signatures observed within and between episodes. In Sr-Pb isotope space (Fig. 14), episode 3 and 4 samples were shown previously to exhibit a negative linear correlation consistent with two-component mixing between a HIMU-type and enriched mantle, EM-type mantle sources, although a gap was apparent in the trend at intermediate isotopic compositions (Rasoazanamparany et al., 2015). The episode 2 basalts have intermediate 87Sr/86Sr and 206Pb/204Pb isotopic compositions that essentially fill this gap and generally plot along the same, albeit somewhat more scattered, trend (Fig. 14). Episode 1 trachybasalts extend the range of 206Pb/204Pb to lower values, but are significantly off the trend toward radiogenic 87Sr/86Sr. As a group, the episode 1–4 samples therefore appear at first glance to be consistent with two-component mixing between HIMU- and EM-type mantle sources, with episode 1 samples exhibiting widely variable and radiogenic Sr isotope signatures due to assimilation of carbonate-rich upper crust (Yogodzinski et al., 1996). Pb-Pb and Nd-Sr isotope systematics, however, indicate a more complex scenario requiring additional source components to explain the isotopic heterogeneity observed among the four magmatic episodes. In a 208Pb/204Pb - 206Pb/204Pb diagram (Fig. 13A), for example, the samples plot along two distinct trends that could reflect three-component mixing between a common HIMU-type source and two distinct low 206Pb/204Pb, low 208Pb/204Pb components; 207Pb/204Pb variations, although significantly more limited (Fig. 13B), are likewise suggestive of the contribution of at least three compositionally distinct sources. Although the Pb isotope signatures of episode 1 samples could be influenced by crustal assimilation (Yogodzinski et al., 1996), the strong correlation among episode 1 samples between Sr and Nd isotopes (Fig. 12), yet complete lack of correlation between Sr and Pb isotopes (Fig. 14), indicates that the Pb isotope variations of episode 1 basalts were partly controlled by crustal assimilation and partly by source heterogeneity.

Sr-Nd isotope systematics show clearly the impact of carbonate-rich crust assimilation on episode 1 samples, and further highlight the distinctions between episode 2, 3, and 4 basalts, each of which is in a distinct compositional field with steep trends relative to that of episode 1 (Fig. 12). Proterozoic lower crust within the Basin and Range is characterized by radiogenic Sr and unradiogenic Nd and Pb (Kempton et al., 1990), as well as positive Δ8/4, and thus could contribute to the isotopic variability among episode 2 samples (Δ8/4 and Δ7/4 refer to the deviation of a sample from the Northern Hemisphere Reference Line in 208Pb/204Pb and 207Pb/204Pb, respectively, at a given 206Pb/204Pb; Hart, 1984). However, the common occurrence of mantle xenoliths within the episode 2 basalts suggests a rapid rise from the mantle to the surface without ponding and storage in the crust, arguing against a significant role for crustal assimilation (e.g., Foland and Bergman, 1992; Yogodzinski et al., 1996). Furthermore, assimilation of Proterozoic lower crust has been ruled out for basalts of episodes 3 and 4, which exhibit behavior in Sr-Nd isotope space similar to that of episode 2 basalts, yet have Os isotope signatures and negative Δ8/4 signatures that are inconsistent with assimilation of local Proterozoic lower crust (Rasoazanamparany et al., 2015).

With the exception of the Sr isotope signatures in episode 1 trachybasalts, the observed trends in Sr-Nd-Pb isotope space are therefore attributed to mixing between at least three isotopically distinct mantle sources, including at least one HIMU-type source with radiogenic 206Pb/204Pb, high 143Nd/144Nd, and low 87Sr/86Sr, and two EM-type mantle sources with low 206Pb/204Pb, low 143Nd/144Nd, and, respectively, positive and negative Δ8/4 and slightly different corresponding 87Sr/86Sr signatures (Figs. 12 and 13). Whether two slightly different HIMU-type sources are implied by the data trends in Sr-Nd isotope space (Fig. 12A), one associated with episodes 2 and 3, and another with episode 4, is unknown. Nevertheless, the apparent mixing relationships in each magmatic episode between a HIMU-type and an EM-type mantle source are consistent with similar processes affecting the mantle beneath the broad LCVF region, producing heterogeneity on a small spatial scale in which compositionally distinct sources are tapped throughout time and space. Our previous modeling of the Pancake Range basalts, suggesting that subduction-related mantle enrichment processes are responsible for the observed compositional and isotopic variations (Rasoazanamparany et al., 2015), may thus be applied to the larger LCVF region. In this model, mixing between a HIMU mantle source containing ancient recycled oceanic crust and an EM source related by addition of slab-derived hydrous fluid to the HIMU source, produced the compositional variations in the episode 3 and 4 basalts of the Pancake Range. A similar process could produce the isotopic signatures of the second EM-type mantle that contributed to the episode 1 and 2 basalts, if the hydrous fluid in this case contained a minor component derived from pelagic sediment (low 206Pb/204Pb, positive Δ7/4 and Δ8/4). In this scenario, a regional HIMU-type mantle beneath the Pancake and Reveille Ranges, further metasomatized by compositionally variable subduction-related fluids from slab ± sediment, could explain the geochemical and isotopic signatures that we now identify within the LCVF as a whole.

Shallow Plumbing and Eruptive Products

The long lifetime of the volcanic field and the northward shifting of eruptive activity have resulted in a wide range of erosional exposure. The Reveille Range and walls of Lunar Crater maar expose some feeder dikes, the transition from dike to vent, and the internal plumbing of eroded scoria cone volcanoes. We use the term feeder dike system for dikes below, or laterally away from, the localized upward flaring that is common where vents formed, but in deeply eroded cases we cannot always clearly link such dikes with eruption centers. The plumbing information can be combined with the relatively pristine, younger volcanic surface features and deposits in the Pancake Range to develop an overall picture of the very shallow (uppermost ∼150 m) plumbing structures and their relationships to eruption processes.

Feeder Dikes

Reveille Range feeder dikes. Feeder dikes exposed in the Reveille Range typically range between 1 and 3 m wide. One example (Fig. 15) can be traced laterally for ∼700 m with an elevation change of 175 m on steep hillsides. In detail the lowest and southernmost exposure is a vertical dike striking N5°E that is ∼70 cm wide with local additional widening of as much as 110 cm, similar to higher altitude exposures where the dike is hosted by partly welded Oligocene ignimbrite (Figs. 15A, 15B). In all places where we observed this dike it exhibits vesicle bands that are symmetrically paired with respect to the dike’s center plane (Figs. 15A, 15B). We interpret this banding to result from multiple pulses of magma through the dike center, where each pulse would push the partly solidified remains of the preceding pulse to the side (see also Hintz and Valentine, 2012). In detail the dike is a series of en echelon segments that parallel the prominent jointing of the host ignimbrite and likely occupy joints; each segment can typically be traced a few tens of meters, often with minor bends in the trace, and lateral step distances between segments are <10 m. However, 240 m to the north-northeast of Figure 15A, and 45 m higher in altitude, the dike steps ∼125 m to the right, and a complex network of smaller (meters to tens of meters long) dikes and sills is present between the en echelon segments. The dike then can be traced clearly up steep slopes where it intrudes ignimbrite and, in the upper parts of the exposure, an eroded scoria cone of episode 1 (Fig. 15C). The dike cuts the scoria cone remnant and postdates it, and might be related to episode 2.

Other feeder dikes exposed in the Reveille Range have less vertical exposure, but are similar in that their total map lengths are several hundred meters, and they consist of en echelon segments that are typically several tens of meters long, and with local small sills and widened zones. One such feeder dike system that is clearly linked with the Dark Peak eruptive center was described in Harp and Valentine (2015); it can be traced ∼908 m along an arcuate trend that parallels a nearby Oligocene-age caldera wall structure. The dike varies in width between 1 and 3 m along its trace. Another dike system intrudes (likely Miocene) volcaniclastic rocks and can be traced 630 m (Fig. 16). It forms a set of en echelon segments that is as much as 100 m wide. Individual segments trend between ∼N10°W and N15°E, with dips between ∼75° and 90° both to the east and west. In detail the dike segmentation occurs on different scales; the larger scale forms 50–200-m-long segments, which are in turn broken into multiple segments of several meters to tens of meters in length (Fig. 16A). It is unclear whether this feeder dike system is related to episode 1 or 2, and it is not directly connected to any preserved eruptive centers.

Magma pressures in Reveille Range feeder dikes. For a dike in equilibrium with a brittle-elastic host medium, magma overpressure (difference between magma pressure, Pm, and horizontal compressive stress perpendicular to the dike plane, here assumed to be the regional σ3) can be estimated by (Pollard, 1987; Rubin, 1995):

Here w is the dike half-width and l is its half-length. The second term on the right side of Equation 1 is the host-rock stiffness, which is determined by the shear modulus (μ) and Poisson’s ratio (ν). Here we assume that the appropriate value of w is the average dike width measured on individual segments (i.e., not the width of the system of en echelon segments), and l is the length of a feeder dike system. We use these values rather than, for example, the local maximum dike widths, which are strongly influenced by host-rock variability and mechanical erosion of dike walls, and because we infer that the length of a dike system is related to the length of the feeder dike below the level at which it begins to break into small segments. For the feeder dike systems described here the typical value of w ≈ 1 m, and l ≈ 600–1000 m, or w/l = 1–2 × 10–3. Elastic properties of the host materials are not well constrained. For feeder dike systems within this range of w/l, hosted in moderately welded ignimbrites (e.g., Fig. 15), μ ≈ 4–7 GPa and ν ≈ 0.35 (reported for silicic tuffs at the Nevada Test Site; see Valentine and Krogh, 2006), or μ/(1 – ν) ≈ 6–11 GPa. These parameter ranges give magma overpressures of 6–22 MPa. Similar feeder dike systems hosted by clastic sedimentary rocks, such as the example in Figure 16, would have overpressures of ∼5–10 MPa, assuming μ/(1 – ν) ≈ 5 GPa (e.g., Delaney and Gartner, 1997). In addition to uncertainties in elastic properties or the host rocks, there are many other caveats with this approach, such as (1) its appropriateness for en echelon dike systems rather than distinct, coherent dikes; (2) the assumption of mechanical equilibrium between the dike-filling magma and its elastic host; (3) whether the host medium is truly elastic, or if it can deform inelastically; and (4) whether the final dike dimensions are reflective of the dimensions during dike emplacement (see also Daniels et al., 2012). The latter is especially problematic because paired symmetric banding in the dikes is suggestive of multiple pulses of magma that contributed to the final dike widths, and because relaxation of the host rock as the magma overpressure decays at the end of an emplacement event can modify its shape. Nevertheless, overpressures on the low end of the ranges stated above seem reasonable.

Reveille Range sills. Sills associated with feeder dikes appear to be localized features at the levels of exposure available; it is not clear if these would have affected eruption processes (e.g., with respect to gas segregation processes; Menand and Phillips, 2006), and it is unlikely that they were of sufficient scale to affect petrologic processes such as shallow fractionation or assimilation (e.g., as occurred at Parícutin volcano, Mexico; Erlund et al., 2010). The sills observed along the dike that fed the Dark Peak eruptive center each occur along only ∼20 m of dike length, and extend laterally as much as 15 m from the dike (Harp and Valentine, 2015). Sills exposed in cross section (Fig. 17) extend decimeters to ∼10 m into the host rocks, ranging in shape from irregular and stubby to thin and tabular. These are similar in scale to small sills observed in the upper ∼100 m of feeder dikes at eroded scoria cone volcanoes elsewhere in the region (e.g., Basalt Ridge; Keating et al., 2008) and likely represent the small-volume end of a spectrum of shallow sill sizes and geometries that can form in the upper few hundred meters of crust beneath monogenetic volcanoes (e.g., Valentine and Krogh, 2006; Németh and Martin, 2007; Erlund et al., 2010; Re et al., 2015; Richardson et al., 2015; Muirhead et al., 2016). Mechanisms for such sills include localized changes in the principal stresses along faults, proximity to the free ground surface which allows mechanical translation of roof rock in addition to elastic deformation, interaction of dikes with bedding plane weaknesses, and others. The complex network of small sills and dikes between (but not connecting, at least with current exposures) en echelon segments of one of the feeder dikes might be related to magma filling in preexisting joints in the host tuff that dilated due to stresses between the segments.

Pancake Range feeder dikes. The younger volcanoes in the Pancake Range do not have exposed feeder dikes, but some inferences on their lengths can be made based upon crater morphologies and the details of individual eruption reconstructions (potential relationships of these dikes with preexisting faults have been discussed). Many, possibly most, scoria cone-forming eruptions begin along a fissure that represents the intersection of the feeder dike with the surface, and as an eruptive event proceeds the magma discharge commonly focuses into one or a few point sources, around which cones build and bury the earlier fissure vents. Most individual scoria cones in the Pancake Range are approximately circular or slightly elongate in plan, with no coeruptive vents beyond the cone perimeters. This means that many initial fissure vents have been completely buried by the cones and their lengths cannot be directly measured, except to say that the cone diameters are the maximum lengths of the early fissures (see also Valentine and Perry, 2006). Many of the scoria cones cluster in elongate patterns, and it can be difficult to determine whether a line of cones is from the same eruptive event, and whether they were fed by the single, contemporaneous dike system. However, through surface mapping we identified six volcanoes that certainly erupted along identifiable fissure vents, forming agglomerate ramparts rather than cones (Table 3). These ranged from 200 to 2500 m long, but the longest ones were composed of three or more en echelon fissure segments (e.g., Easy Chair; Valentine and Cortés, 2013). Assuming that the fissure lengths were similar to the feeder dike lengths (realizing that the feeder dikes might be longer, and that only a portion reached the surface), and width:length ratios were similar to those exposed in the Reveille Range, these fissure lengths suggest dike widths of several decimeters to a few meters beneath Pancake Range volcanoes.

Feeder Dike to Vent Transition and Conduits

Observations. We define the transition from feeder dike to vent to occur where the dike begins to flare as it approaches the Earth’s surface. The flared portion is the vent (or conduit, without implication of plan form shape), and may contain any combination of coherent basalt, breccias, and pyroclastic material. Three locations in the LCVF expose volcanic edifices, which allow some inferences on eruption styles, and the transition from feeder dike to vent (Fig. 18). At Dark Peak (Fig. 4; Harp and Valentine, 2015), a remnant of an episode 1 scoria cone, the feeder dike transitions upward from a width of 1.4 m to an ∼18-m-wide body containing domains of agglomerate and coherent basalt in the upper 15 m of crust (Figs. 18A, 18B). The Dark Peak eruption is interpreted, based upon its deposits, to have involved Strombolian to Hawaiian explosive styles, along with effusion of extensive lavas.

Excavation and collapse of the Lunar Crater maar exposed a cross section of an agglomerate rampart that extended at least 300 m to the north, approximately perpendicular to the cross section (we reserve the term spatter for deposits composed nearly entirely of fluidal and/or ragged clasts and moderately to densely welded, preferring to use the more general term agglomerate for deposits dominated by bombs but possibly with some blocks, and any degree of welding; Fisher and Schmincke, 1984). The agglomerate rampart is part of episode 4 (sample GVLC13–53; Table 1). At its deepest exposed level, the feeder dike is ∼50 cm wide; the outer 15–20-cm-wide margins exhibit paired vesicle bands similar to those described here, but not as clearly developed, while the central portion of the dike consists of discontinuous centimeter- to decimeter-size vesicles and cavities with irregular shapes. These vesicles and cavities are lined with ropy, fluidal textured basalt. The feeder dike flares to ∼2 m width in the upper 2 m below the paleosurface, and extends a few meters into its own eruptive products in the form of densely welded agglomerate with relict clasts (Figs. 18C, 18D). The transition from coherent basalt in the dike to the densely welded agglomerate is gradational over decimeters. The densely welded agglomerate in turn grades outward over a distance of meters into moderately to poorly welded, crudely bedded agglomerate, which is overlain by a lava flow (likely clastogenic) in the outer slopes of the rampart. We interpret the agglomerate rampart to have formed by accumulation of coarse clasts ejected during low fountaining or spattering with fountain heights of a few tens of meters or less (e.g., Parcheta et al., 2012). Magma flux through the narrow feeder dike was probably relatively low, with slugs of gas rising to feed the fountains. Cavities in the central portion of the dike record the final slugs as they were locked in place by the solidifying magma.

An episode 2 vent in the Reveille Range records a similar, very rapid upward transition from a decimeters-wide feeder dike into the edifice that it built (Figs. 18E, 18F). Here the dike flares no more than a few decimeters as it approaches the paleosurface (top of volcaniclastic sediments; Fig. 18F). Above the paleosurface, in the volcanic edifice, it transitions into a 3–4-m-wide body of coherent basalt that probably records a very small lava pond in the crater of the agglomerate cone that accumulated around it. The eruptive processes associated with this feeder dike and vent were likely Strombolian to Hawaiian, similar to the preceding two examples.

Comparison with other dike-to-vent data. A growing body of data shows that the widening of feeder dikes and transitions to vents for Hawaiian and Strombolian eruptions is a very shallow process, limited to the upper ∼100 m of crust, and in most reported cases to the upper 20 m (Keating et al., 2008; Geshi et al., 2010; Geshi and Neri, 2014; Valentine, 2012; Friese et al., 2013; Harp and Valentine, 2015; Parcheta et al., 2015; Reynolds et al., 2016). Our observations of this transition in the LCVF are consistent with these published data. This geometry strongly contrasts with the shallow plumbing of volcanoes that likely had similar deep (>1–2 km) feeding systems, but that underwent dominant phreatomagmatic activity and formed large craters such as Lunar Crater maar, with vent structures (diatremes) that can extend downward to 2 km below the surface (White and Ross, 2011; Valentine and White, 2012; Brown and Valentine, 2013). The difference in vent structures between phreatomagmatic- and magmatic volatile–dominated volcanoes is due to strong subsurface magma-water explosions and resulting radially directed stress waves in the former, compared to more gradual upward acceleration of magma or gas slugs in the latter. Unfortunately there are no eroded maar diatremes in the LCVF to provide direct comparisons within the volcanic field.

Edifice Types and Their Internal Structures, and Maar Craters

Four types of volcanoes are well preserved in products of episodes 3 and 4 of the Pancake Range: simple central-vent scoria cones, ramparts, spatter cones and mounds, and maars. Maars are related to explosive phreatomagmatic activity, while the other three types were formed by eruptions that were driven primarily by magmatic volatiles, based upon the overall characteristics of the deposits (e.g., Valentine and Connor, 2015). Eroded cones and ramparts in the Reveille Range expose their internal plumbing structures. Two episode 2 trachytic and trachyandesitic lava domes form a fifth vent type.

Cones, ramparts, and mounds. Central-vent cones are circular to elliptical in plan, and range in basal diameter from ∼300 to 1200 m (most are in the range ∼600–1000 m), and heights of ∼40–200 m. Cones near the smaller end of the size range occur in the Lunar Lake basin and might be partly buried beneath alluvium and playa deposits. Surface processes have substantially modified all but the youngest cones, so their morphology is no longer strictly representative of the primary volcanic processes that produced them. The cones are composed of variably welded agglomerate consistent with Strombolian to Hawaiian eruption processes. Craters that are open on one side are common and often the result of cone fragments having been preferentially rafted away on top of lavas on one side of the cone (e.g., the downhill side if the cone is erupting on a slope or escarpment; e.g., Holm, 1987; Riggs and Duffield, 2008; Valentine et al., 2006, 2007), and/or of the effects of strong wind on pyroclast dispersal. Many older craters may be breached in one quadrant due to posteruptive erosion processes.

Fissure-fed eruptions form elongate agglomerate ramparts. These are typically several hundred meters long and tens of meters high (Table 3) and consist of en echelon ramparts except in the shortest examples. Early eruptive activity at the ca. 35 ka Marcath volcano (Figs. 5, 19A, and 19B) began along a fissure system (Table 3) and built elongate ramparts on either side of the fissures. Activity focused toward the northern end of the system, building a prominent cone that partly buried the earlier formed ramparts; the total volume of cone and rampart is ∼4 × 107 m3 (note that this corrects an earlier reported volume of 9 × 107 m3 in Johnson et al., 2014, that included a geometry error), which equates to a dense-rock equivalent (DRE) volume of ∼2 × 107 m3, assuming 50% porosity of the deposits. It is unclear whether a physically separate agglomerate rampart just to the north of the main Marcath cone is associated with the same eruption; the petrographic characteristics of the two are similar, but this is not diagnostic because most episode 4 volcanoes are petrographically similar. Easy Chair volcano (ca. 140 ka; Fig. 5) forms the longest set of agglomerate ramparts from a single eruptive event (Figs. 19C, 19D). Its eruption also began along a fissure system (Table 3), and then focused toward the center to build two overlapping scoria cones from which a small lava field emanated. Late explosive phreatomagmatic activity excavated the southern ramparts during maar crater formation (Valentine and Cortés, 2013). The ca. 130 ka Giggle Springs eruption formed low agglomerate ramparts along a 700-m-long fissure system (Table 3; Figs. 5, 19E, and 19F). This eruption did not focus toward a point source, but fed lavas that issued from both ends of the ramparts.

Additional landforms record weakly explosive to nearly purely effusive venting. Some spatter cones and mounds are associated with relatively young scoria cones and record weak, closing phases of eruptive activity. They are typically several meters high and can have steep sides with sharp peaks (e.g., Fig. 20) or more rounded shapes; the latter may have originated as more conical shapes but have partly collapsed and disaggregated. In one case (see Table 3; partly shown in Fig. 20) the spatter cones and/or mounds are arranged in an ∼500-m-long line that extends southwestward from the floor of a larger (200 m diameter) crater. One isolated vent, ∼2 km north-northeast of the Marcath cone, is an elongate spatter mound or rampart. It forms a very low profile (∼5 m high, ∼100 m wide) ridge ∼200 m long, trending north-northeast, and might have fed a relatively young lava that extends 2.5 km to the southwest and is partly covered by Marcath lava; the large number of overlapping centers in the area and partial blanketing by Marcath tephra and alluvium make it difficult to be certain of the lava’s origin. It is probable that there have been other such vents in the LCVF, but they were rapidly obscured by products from other centers due to their low topographic expression.

Internal plumbing of cones and ramparts. Old eroded centers demonstrate that cones and ramparts often have complex internal plumbing. Hintz and Valentine (2012) described the Kimana center, a remnant of a scoria cone (episode 1, Reveille Range; Fig. 4) that had a basal diameter of ∼1 km and a probable height of 250–300 m. The summit of Kimana is an ∼50-m-diameter area of moderately to densely welded lapilli that we infer records the main vent structure within the cone. Eight dikes and dike sets form a radial pattern around this inferred vent. Although their overall geometry is tabular with near vertical dips, in detail the dikes are composed of many segments that are meters to a few tens of meters long. Typical thicknesses are ∼1 m, but the shapes of segments can be quite irregular. This is probably related to the magma, which was variably degassed, crystallized, and therefore of high viscosity, being pushed into unconsolidated scoria and agglomerate (Mathieu et al., 2008; Delcamp et al., 2014). Many of the dikes exhibit paired vesicle bands and internal quench margins, and in some cases have central portions with larger (to ∼10 cm), irregularly shaped vesicles. Scoria cone remnant Dark Peak (episode 1, Reveille Range; Fig. 4), described in Harp and Valentine (2015), exposes two irregular dikes within its edifice, in addition to the transition from feeder dike to the main vent (Figs. 18A, 18B). One of these radiates from the main vent and, at its distal end where it neared the outer cone slopes during emplacement, produced bulbous intrusions several meters in size (small cryptodomes; see also Delcamp et al., 2014) and fed lava flows that vented near the cone’s base. Other eroded centers in and near the Reveille Range, although not studied in detail, exhibit similar complex internal plumbing with irregular radial dikes and bulbous intrusions.

The radial dikes are interpreted to record lateral intrusion of magma driven by magmastatic pressure when magma levels were relatively high in the cone conduits. Lateral intrusions were particularly facilitated by dense, degassed magma in the conduits, which resulted in larger magmastatic pressures for given magma column heights; in many cases, the intrusion was pulsatory, as recorded by paired vesicle banding, due to pressure fluctuations in the conduits (Hintz and Valentine, 2012; Harp and Valentine, 2015). Alternatively, radial dikes might have resulted from uplift or bulging of the volcanic edifices due to intrusions, a process that opens up radial fractures that could facilitate magma migration and dike formation (van Wyk de Vries et al., 2014). Some of the radial intrusions breached the outer flanks of the cones and fed lavas, similar to those observed in historical scoria cone eruptions (e.g., Parícutin; Luhr and Simkin, 1993). In some ways the internal structures of scoria cones are scaled-down versions of stratovolcano plumbing (e.g., Swanson, 1966; Nakamura, 1977; Poland et al., 2008).

Maars. The fourth type of vents in the LCVF is maars, inferred to have formed by phreatomagmatic explosions. Of the 96 volcanoes in the Pancake Range, 4 maars have been identified (none have been identified in the Reveille Range, but it is possible that any maars there were buried in the surrounding basins). The simplest of these is Lunar Crater maar, which is an ∼1100-m-diameter, 130-m-deep, nearly circular crater surrounded by a low-profile tephra ring (Figs. 5 and 21A; Valentine et al., 2011). The pre-maar landscape was a gently north-sloping surface underlain by episode 2 lavas, with episode 3 and 4 lavas, cones, and ramparts. These basaltic rocks are underlain by Oligocene silicic ignimbrites. The floor of the maar crater is ∼100 m below the paleosurface, thus exposing the lavas and uppermost ignimbrite in the crater walls as well as the aforementioned rampart (Figs. 18C, 18D), and a scoria cone and dikes in its southern wall (Fig. 21B). Tephra ring deposits are poorly exposed in the crater walls due to colluvium, and are covered by a veneer of eolian sediments beyond the crater rim. The limited exposure exhibits bedded and cross-bedded lapilli tuffs composed almost entirely of clasts derived from the pre-maar lavas; ignimbrite clasts compose as much as 25% of the deposits (Valentine et al., 2011). Blocks as large as 3.3 m were ejected as ballistic clasts and occur within, and litter the surface of, the proximal tephra ring (Valentine et al., 2011). The only clearly juvenile material observed was an occasional thin coating of basalt on ignimbrite clasts. Quartz from in situ exposures of the silicic ignimbrite in the crater wall and blocks on the tephra ring surface were analyzed for cosmogenic 10Be exposure dating, which constrains the minimum age of the maar (Table 1). Lunar Crater maar is partly filled by talus from the crater walls and alluvium brought in by an ephemeral stream that enters the southwest quadrant of the maar.

Another maar, informally referred to as Bea’s Crater, preserves a more complex eruptive history (Figs. 5, 21C, and 21D; Amin, 2013). It occurs in the southeastern part of a dense cone and rampart cluster that includes Marcath volcano, and as a result it disrupted at least three earlier edifices. The crater is actually composed of two overlapping craters (∼440 m and ∼1050 m diameter), the centers of which form a line that is parallel to the dominant north-northeast trend of ramparts and vent alignments in the area. One of the pre-maar volcanoes was an agglomerate rampart system that was likely elongated in the same north-northeast direction as nearby ramparts (such as the early deposits of the Marcath center). Most of this rampart system was located within the footprint of Bea’s Crater and was destroyed, but the outer flanks are exposed along the southeast quadrant of the crater wall. The western part of the crater disrupted a cone composed of moderately to densely welded agglomerate that was relatively competent and thus forms small cliffs along the crater wall. Tephra ring deposits are exposed along the southern rim of Bea’s Crater, but are covered by younger Marcath tephra in the northern part. The tephra ring deposits are massive to crudely bedded lapilli tuffs with lapilli, bombs, and blocks as large as 1.5 m, all supported in a coarse ash matrix. They are thickest (as thick as 15 m) along the southeastern part of the tephra ring. In most places 70%–80% of the blocks are angular to subangular dense basalt, while the remaining portion consists of angular to subrounded rhyolitic ignimbrite. The crater rim is interrupted by a cone in its northeastern sector. Variably welded agglomerate associated with the cone both underlies and overlies the tephra ring deposits, indicating that pyroclastic fountaining activity at the cone, driven by magmatic volatiles, was contemporaneous with phreatomagmatic maar-forming activity. The cone partly rebuilt itself after the phreatomagmatic activity, extending a new scoria cone rim into the maar (Amin, 2013), and possibly a lava flow, the top of which is exposed on the crater floor.

Easy Chair was mentioned here in the context of agglomerate ramparts (Figs. 5, 19C, and 19D); this ∼540-m-diameter maar that formed at the southern end of the ramparts represents one of the closing phases of activity there. Its tephra ring is draped onto the outer rampart slopes that formed earlier in the eruption, and thus it has a steeper surface than many tephra rings. The deposits are similar to those at Lunar Crater and Bea’s Crater, comprising crudely bedded lapilli tuffs with blocks derived from underlying rocks, in this case agglomerate associated with the earlier rampart-building phases, and pre–Easy Chair lavas. Clasts derived from older ignimbrites are rare (Valentine and Cortés, 2013). A thin veneer of scoria on top of the tephra ring deposits indicates minor Strombolian activity, probably vented from the central cones in the Easy Chair system, after phreatomagmatic explosions ended.

A final maar-like landform that is inferred to be related to phreatomagmatic or hydrothermal explosions is informally called Middle Maar (coordinates of crater center: 38°25′10″N, 116°01′03″W). This ∼550-m-diameter crater is on top of a preexisting cluster of scoria cones, such that its floor is ∼100 m above the surrounding valley floor (we refer to it as a maar because in detail its preexisting landscape was the scoria cone cluster). Material that was likely ejected during the crater-forming explosions is limited to one small, 50-m-wide drainage in the outer northeast slopes of the cone cluster. The deposits are crudely bedded and cross-bedded lapilli tuffs with isolated meter-sized blocks. Lapilli and blocks are composed of dense basalt and welded agglomerate, likely derived from the cone cluster through which the crater was excavated. Although deeper seated ignimbrite clasts were not observed in situ in the deposits, they are scattered in the colluvium at the site and probably weathered out of the lapilli tuff. The lack of other outcrops and eolian cover precludes a more detailed understanding of Middle Maar.

Discussion of maar occurrences. The four maars compose ∼4% of the total number of centers in the Pancake Range; this is typical for volcanic fields in the intraplate western United States (Valentine and Connor, 2015), and the question arises, were the maars related to unique hydrologic or climatic conditions that promoted subsurface magma-water interaction? The maars are all within a few kilometers of the lowest part of the Lunar Lake basin. Lunar Crater formed in a bedrock setting (bedrock here being older basaltic lavas and Oligocene ignimbrite) ∼30 m above the Lunar Lake playa, although its crater floor is lower than the playa surface (Valentine et al., 2011). Bea’s Crater is ∼5 km north of the playa and its floor is ∼80 m above it. Easy Chair, Bea’s Crater, and Middle Maar all formed on older clusters of basaltic cones and ramparts that must interfinger with basin sediments. The top of the Oligocene ignimbrites is 200–400 m below the surface (Snyder et al., 1972; Valentine and Cortés, 2013; Amin, 2013). Most of the volcanoes in Lunar Lake basin are expected to share this shallow subsurface geology, but most did not have significant phreatomagmatic activity. Thus it is difficult to speculate about any systematic hydrologic conditions that would have promoted subsurface phreatomagmatic explosions.

Climatic conditions that might have promoted such activity also are difficult to ascertain without better constraints on volcano ages with respect to potential lake stands in the past few hundred thousand years. The 140 ka age of Easy Chair is near the end of Marine Isotope Stage (MIS) 6, which was a relatively cool and wet time in the region, based upon regional paleoclimate records (Winograd et al., 1997; Jiménez-Moreno et al., 2007). One of these records, Devil’s Hole, is 220 km to the south, while the other, Bear Lake, is 550 km northeast of the volcanic field. To our knowledge, little paleoclimate research spanning the interval of interest has been published for the interior Great Basin area around the LCVF. The age of Lunar Crater maar, 190–72 ka, spans both MIS 5 (dry) and 6 (wet), although we infer that it more likely occurred during MIS 5, based upon qualitative geomorphic features (Valentine et al., 2011). It is important to note that in the cases of Easy Chair and Bea’s Crater, styles of activity alternated within the context of the individual monogenetic events. This indicates that variables other than local geohydrology or climate, such as magma flux during an event, probably played the major role in determining whether phreatomagmatic or volatile-driven eruptions occurred, given the presence of groundwater (e.g., Houghton et al., 1999; Gutmann, 2002; Valentine and Cortés, 2013).

Lava domes. Remnants of two lava domes of trachytic and trachyandesitic composition occur in the northeastern part of the Reveille Range, and have ages that correspond with early episode 2; they are described briefly here (including brief petrographic descriptions, which are not elsewhere discussed for these minor products of the LCVF). The largest example corresponds to sample R8–1-43-LN (Table 1). Here dense trachyte lava is underlain by massive lapilli tuff with scattered blocks of trachyte lava and older basaltic lithics, possibly recording early explosive activity. The lava is fine grained with dominant sanidine phenocrysts a few millimeters in size in a cryptocrystalline groundmass. The main body of the lava dome is an ∼700-m-diameter, irregularly shaped stack of 2 or more trachyte lavas totaling ∼100 m thick (centered at N38°07′45″N, 116°06′07″W). We infer that this body corresponds to the vent location for the lavas. Two major lobes or coulees extend 700 and 900 m to the east from this body and a third coulee extends 800 m to the north. The coulees are 70–80 m thick. We infer that viscous trachytic lavas extruded on the eastern flank of the Reveille Range, first forming a thick exogenous lava dome. The two eastern coulees flowed downslope toward the Railroad Valley. The northern coulee might have been diverted around one of the thick early lobes in order to flow northward. A smaller, ∼140 m diameter and ∼20 m thick, remnant of a trachyandesitic dome is preserved 8 km to the south-southeast (sample R8–1-16-LN; Table 1). Lava in this dome is dense with jointing (meters spacing) that is suggestive of having originated as shear lobe boundaries in some cases. The trachyandesite is fine grained, and contains plagioclase, sanidine (both <5 mm), and clinopyroxene (<3 mm) phenocrysts often with glomeroporphyritic texture. Isolated exposures on the hillside below the dome expose as much as 50 m of trachyandesitic tuff breccia that might record a talus apron around the dome. This small dome grew on a south-facing slope that was underlain by coarse basaltic talus derived from episode 1 volcanoes on an ∼200 m paleohigh just 800 m to the north. An additional trachyte lava flow is present beneath episode 2 lavas in the area, but it is only exposed in cross section as an ∼10-m-thick massive lava with minor basal flow breccia, and is not discussed further.

Eruptive Products and Processes

Lavas. All lavas that we observed were either emplaced primarily as single ‘a‘ā flows or as compound ‘a‘ā fields. Flow distances for most episode 3 and 4 lavas range from ∼1 to 5 km, and volumes are ∼3 × 107 to 2 × 108 m3 (estimated only for flows whose margins and thicknesses have not been obscured or eroded). Older episode 3 and most episode 1 and 2 lavas are either partly buried or eroded, and have sufficient surface modification to overprint many emplacement-related features.

Marcath lava. The Marcath lava field (Fig. 22) has a bulk volume of ∼108 m3 (based upon an average thickness of 15 m), similar to that of the cone and rampart. This corresponds to a DRE volume of ∼8 × 107 m3, assuming an average of 20% vesicularity of the lava. The flow field preserves many features related to its emplacement despite its ca. 35 ka age, reflecting the slow initial pace of surface modification of lavas in this high desert setting (e.g., Wells et al., 1985; Valentine and Harrington, 2006; Valentine et al., 2007). Lava emanated from the Marcath cone, which is situated on the edge of the aforementioned ∼100-m-high escarpment that steps down to the west toward Big Sand Springs Valley. The lavas flowed westward from the crater into a saddle between the escarpment and a cluster of episode 3 cones, forming the flat-topped lava upland (Figs. 22 and 23A). Lapilli and bombs that blanket the surface, and interstitial eolian sediments, obscure detailed surface textures of the upland. Depressed areas that occur near the edges of the lava upland are typically a few tens of meters across, and are several meters lower than the main upland surface. Flows first broke out from the northern side of the lava upland, forming an ∼120-m-wide channel that fed lava down the escarpment and around the older cones onto the flat valley floor, spreading out to form the northern and central lobes. Subsequent lavas flowed southwestward in similar channels, as narrow as 40 m, to feed the southern lobe, which abuts and conforms to the edge of the central lobe (Fig. 22). The lava surfaces are rubbly ‘a‘ā (Fig. 23B) and preserve arcuate ridges (also known as ogives, without genetic connotation), ∼1 m high, that widen downflow on the major lobes. These ridges are especially clear on the ∼700-m-wide central portion of the southern lobe, and an ∼500-m-wide central part of the northern lobe (Fig. 22), but vestiges of such large-scale wrinkles are found toward the outer edges of the lobes as well. Linear depressions in the flow surfaces parallel the lateral edges of the arcuate-ridged portions of the field. Such depressions also occur in the distal hundreds of meters of the lobes, where they roughly parallel flow directions. Smaller lobes around the edge of the lava field extend tens of meters outward onto the desert floor. The lobate flow margins rise 5–15 m from the desert floor (Fig. 23C).

We infer that lavas first ponded in the saddle to form the lava uplands, possibly inflating, as suggested by the flat surface, as a pāhoehoe sheet (Self et al., 1998). Initial northbound lava may have been sourced from the inflated interior of the lava upland, but eventually lava channels connected directly to the vent area in the crater of the Marcath cone, to which they can be traced (this is especially evident toward the head of the southern lobe). Lavas accelerated down the escarpment through these narrow, 40–120-m-wide channels, promoting their development as ‘a‘ā flows. The channels widened to several hundred meters when advancing lavas encountered the relatively flat valley floor. Arcuate ridges suggest broad fluid cores across the channels, and linear depressions along channel edges might result from shear zones along nearly solidified margins (similar to block lavas; Harris and Rowland, 2015). The major lobes grew outward by advancement of smaller lobes around their margins. Lava was fed to these small lobes from the broad fluid cores. In a way, lava emplacement was hybrid between pāhoehoe and ‘a‘ā, in that the flows advanced across flat surfaces by advancement of lobes that stacked and interfingered and eventually merged as the broad fluid cores reached a given location; however, the surface textures are those of ‘a‘ā.

All three major lobes of Marcath lava have rafts of agglomerate on their surfaces, but they are especially abundant on the southern lobe. The rafts are several meters to 20 m wide and as high as 5 m. They range from coherent bodies of bedded lapilli and bombs with different degrees of welding (Fig. 24A) to crudely bedded or massive, variably welded agglomerate (Fig. 24B), to loose mounds of lapilli and bombs. Some rafts have poorly welded or nonwelded interiors but are partially coated with a decimeters-thick layer of lava (Fig. 24C). We interpret the rafts to be material derived from the cone that was carried away on top of the lavas. The rafts shifted in orientation to varying degrees during transport, and some of them acquired lava coatings as they bobbed and rotated on the flow (accretionary lava balls of Lockwood and Hazlett, 2010). Most rafts are in early stages of degradation, surrounded by talus and coherent chunks that disaggregated from the raft during transport and subsequent weathering.

Rafting records a competition between growth of a cone by accumulation of pyroclasts, and destruction by failure and conveyance of failed material out onto a lava field (e.g., Holm, 1987; Luhr and Simkin, 1993; Sumner, 1998; Gutmann, 2002; Valentine and Gregg, 2008). The final configuration of the Marcath cone features a crater that is open to the west. Rafts might have originated from failure of the western crater wall if the wall existed during phases of eruptive activity. Alternatively, if the crater were open for the duration of activity, the rafts might have originated from periodic collapse of its oversteepened inner walls. Mounds of variably oriented, crudely bedded agglomerate in the crater near the breached side (Figs. 22 and 24D) appear to have been stranded just before exiting the crater, as lava flows waned. These mounds consist of abundant vesicular, fluidal bombs as much as 5 m long, and abundant subspherical dense bombs. The latter have diameters ranging from decimeters to ∼1.5 m, and broken examples exhibit concentric layers of densely welded spatter. The cores of these bombs are in many cases simply tightly wound welded spatter, and in other cases are loose scoria lapilli or small bombs; many of them have radial prismatic joints. The subspherical bombs are interpreted to result from rapid accumulation of spatter on steep, unstable slopes. The spatter rolled back toward the vent, accreting material as it traveled (accretionary bombs; Heiken, 1978; Valentine and Cortés, 2013). The presence of large fluidal bombs and accretionary bombs in these stranded rafts is consistent with an origin of the rafts as domains of failed, steep crater walls proximal to the vent.

Comparison with other LCVF lavas. Prior to the Marcath eruption there appears to have been an ∼100 k.y. gap in eruptive activity; the next youngest lavas are associated with the Giggle Springs event (ca.130 ka; Table 1). (Note that the Marcath lavas partly buried two other lavas with relatively young surface characteristics, but we did not obtain ages for these.) Thus older lavas have undergone substantially more surface modification that has obscured many details that are present on the Marcath lavas. However, the broad characteristics of all the lavas seem similar, with steep flow fronts of several meters to ∼25 m high, lobate margins, and agglomerate rafts. Lavas that were erupted from topographic highs (e.g., the escarpment along the northwestern edge of the volcanic field) formed narrow, 40–100-m-wide channels on slopes of more than a few degrees, which broadened to several hundred meters width as the lavas spread onto flat (∼1° or less) valley floors. Lavas erupted from vents in low-relief areas such as Lunar Lake basin spread immediately upon leaving their sources to form broad lava platforms around the cones, similar to those in the Crater Flat basin 200 km to the south (Valentine et al., 2006). Most lava fields have at least a few agglomerate rafts, in various stages of disaggregation, on their surfaces. The Giggle Springs lava field in particular has numerous rafts and most of them have lava coatings. This lava field was fed by two narrow channels, ∼30 and 50 m wide, that curved between older cones on their way down the escarpment; perhaps relatively high speeds and turning flow promoted rotation of rafted domains such that they were coated with lava. Older episode 1 and 2 lavas are typically eroded around their edges and exhibit cross-sectional characteristics typical of ‘a‘ā lavas, with rubbly tops and bottoms and dense cores. Where exposed, lavas that are proximal to their vents often have abundant relict, flattened clast shapes that are commonly sheared and deformed (e.g., Valentine and Cortés, 2013; Harp and Valentine, 2015). These textures do not generally occur in medial and distal exposures.

The general similarities of lava fields in the LCVF may suggest similar lava flow dynamics for most of the eruptions, with relatively high effusion rates to promote ‘a‘ā flows but relatively short flow durations, probably on the order of months in most cases, reflecting the relatively small volumes, 3 × 107 to 2 × 108 m3, of the lavas (compare with the lava volume of 7 × 108 m3 for the historical Parícutin lavas, Mexico; Ort et al., 2008). Development of ‘a‘ā was also promoted where lavas flowed down slopes such as the escarpment near the Marcath cone, and the rugged paleotopography of the Reveille Range (Hintz and Valentine, 2012; Harp and Valentine, 2015). Note that effusion rate here is used as the lava volume flux feeding a flow, and does not imply that they were all fed by purely effusive activity. Many, if not most, of the lavas in the LCVF likely had significant clastogenic components fed by lava fountaining.

Proximal (cone, rampart) pyroclastic deposits. Pyroclastic deposits on inner crater walls of cones and ramparts are poorly bedded agglomerates and scoria lapilli. The deposits are dominated by juvenile clasts, many of which are recycled and/or composite bombs as evidenced by internal clastic textures (White and Houghton, 2006), while lithic clasts are generally sparse. Agglomerates include abundant fluidal and ragged bombs (e.g., Fig. 25A) with individual ribbon, spindle, and cowpat bombs as much as 5 m long. Welding ranges from nonwelded to densely welded, lava-like horizons a few meters thick and extending tens of meters laterally (Fig. 25B; Valentine and Cortés, 2013; Harp and Valentine, 2015). Small bomb and scoria lapilli beds (e.g., Figs. 25C, 25D) tend to be less welded than agglomerates. Deposits forming outer slopes of cones and ramparts are rarely exposed, but are generally better bedded with lenticular reverse-graded beds indicative of grain avalanches, and are mostly nonwelded and finer grained than the inner crater deposits. The eastern foot of the Marcath cone is littered with rounded bombs as much as 2 m in diameter that rolled down the slopes during cone growth. Many of these appear to be accretionary bombs, and their final size may not be indicative of the maximum bomb size that was ejected out of the crater. In general the characteristics of proximal cone and rampart deposits are consistent with Strombolian to Hawaiian eruptive styles driven by magmatic volatiles (e.g., Valentine and Connor, 2015).

Tephra fall. As with lavas, tephra associated with the Marcath eruption is the best preserved in the LCVF because younger volcanic deposits have not yet buried it, although eolian processes have modified the finer grained distal and upper portions of its deposits. The Marcath tephra extends away from the cone in two lobes, a minor one to the south (Fig. 26A) and the more voluminous one to the east. The relationships between the two lobes, such as whether they were originally part of a broad single sheet, or if one preceded the other, are uncertain due to erosion and topography in the intervening area. Johnson et al. (2014) documented the eastern lobe, which consists of a single massive bed of well-sorted scoria lapilli in its proximal and medial reaches. Many clasts in relatively proximal areas are jigsaw pieces of bombs that broke apart upon landing. The proximal thickness of the eastern lobe (not including the cone) exceeds 4 m. The deposit is quite homogeneous at any given location, with no apparent grain-size grading, and appears to have formed during a single period of steady fountaining at the cone. Although the deposit is essentially free of lithics, megacrysts of plagioclase, olivine, clinopyroxene, and amphibole exist as free dense, unbroken clasts, which allowed estimation of eruption column height at ∼7 km using the approach of Carey and Sparks (1986). This height is similar to the maximum modeled height of the 1959–60 Kilauea Iki eruption (Parfitt and Wilson, 1999), although the highest fountaining of that eruption as visually recorded was only reported as >660 m (it is not clear if this value includes a buoyant plume above the fountain; Richter et al., 1970). The bulk volume of the eastern Marcath tephra lobe is ∼2 × 107 m3. Its DRE volume is ∼6 × 106 m3, based upon a representative bulk deposit density of 750 kg/m3 (Johnson et al., 2014). This compares with the DRE volumes of ∼2 × 107 m3 and ∼8 × 107 m3 for the cone-rampart and lava field, respectively. Thus most of the erupted mass fed lavas or was deposited proximally from fountaining, while the component of the eruption column that lofted clasts for distribution in the tephra blanket represented ∼10% of the total mass. If the lavas are assumed to be mainly clastogenic (fountain fed) as well, this means that 90% or more of the erupted mass is associated with lava fountains; this aspect is also similar to the Kilauea Iki event.

Exposures of older tephra fall deposits are limited. A tephra quarry along Highway 6, near ephemeral streams that drain the Big Sand Springs Valley, exposes a >2-m-thick, well-sorted, massive scoria lapilli deposit (Fig. 26B). A bedded horizon of scoria and silt near its top suggests minor fluvial reworking prior to deposition of an overlying thinner tephra bed. The origins of the two tephra beds are uncertain, and could be from any of the several nearby scoria cones. They are overlain by bedded and cross-bedded deposits of silt, sand, lapilli, and bombs with cut-and-fill structures indicative of fluvial processes associated with the Big Sand Springs Valley drainage, and finally by a massive deposit of silt, sand, lapilli, and bombs that grades upward into a desert soil. This latter material is probably colluvium derived from degradation of a small cone that is only 300 m to the south. Such interlayering of tephra fall, fluvial deposits, and colluvium, along with other material such as playa deposits and lavas, is probably common in the subsurface of the Pancake Range. In the Reveille Range more advanced erosion has produced inverted topography, and tephra fall deposits are only rarely exposed beneath mesa-capping lavas (e.g., Fig. 26C). These locations were topographic lows at the time the tephras were emplaced, so like the previous example they are variably reworked and interbedded with fluvial or colluvial deposits. The exposures, although limited, suggest that tephra-producing eruptions have been common throughout the history of the LCVF.

Maar-related deposits. Tephra ring deposits around the maars consist of lapilli tuffs that are dominated by shallow-derived lithics at all four maars, and isolated outsized blocks (Fig. 27). Lapilli tuffs preserve evidence of emplacement from pyroclastic density currents and ballistic curtains (Graettinger et al., 2015a), including (1) massive, poorly sorted beds (e.g., Figs. 27A, 27B); (2) thickening of deposits across pre-maar topographic lows; (3) lenticular block concentrations, likely recording ejecta rays (Graettinger et al., 2015b); and (4) local low-angle cross-stratification. Isolated outsized blocks were ejected and emplaced as individual ballistics. Block size data from Lunar Crater maar were used to estimate possible ejection velocities, which ranged between 50 and 195 m/s; the number of uncertain parameters in the calculations limited more refined estimates (Valentine et al., 2011). The general thinning of individual lapilli tuff beds upward in the section at the Easy Chair maar could be attributed to the tendency of the developing crater to focus eruption jets vertically (Taddeucci et al., 2013; Graettinger et al., 2014; Valentine et al., 2015), thus reducing the amount of material that was emplaced outside the crater for a given blast.

The dominant lithic clast types at all of the maars are previously erupted basalts either from nearby volcanoes or earlier phases of the same eruptive episode, which are also the shallowest seated (Figs. 27A, 27B). Ignimbrite clasts are most abundant at Lunar Crater, where they are sourced from as little as 30 m below the paleosurface, and at Bea’s Crater (Figs. 27C, 27D). Amin (2013) estimated the depth to the top of the ignimbrite section to be ∼250 m, but it is possible that ignimbrites are shallower, given the proximity of Bea’s Crater to the inferred horst along which most episode 4 volcanoes erupted (see discussion of Relationships of Preexisting Structures, Vents, and Feeder Dikes in the Pancake Range). Ignimbrite lithics are sparse at Easy Chair and Middle Maar, where their depth is inferred to be 200–400 m (Valentine and Cortés, 2013). Traditionally, the deepest-seated lithic clasts in tephra rings have been interpreted to record the deepest level of subsurface phreatomagmatic explosions and therefore the bases of underlying diatremes (e.g., Lorenz, 1986; Jordan et al., 2013). Recent field, experimental, and modeling studies, however, indicate that these clasts only record a minimum depth of a diatreme beneath a maar (Valentine and White, 2012; LeFebvre et al., 2013; Graettinger et al., 2014; Valentine et al., 2014; Sweeney and Valentine, 2015). Deep (>200–250 m) explosions probably do not vent at the surface, but mix brecciated wall rock upward (and downward) within a developing diatreme until shallower explosions can eject the clasts. Thus the low abundances of relatively deep-seated ignimbrite clasts in Easy Chair and Middle Maar could be an artifact of: (1) the maars having undergone a relatively smaller number of subsurface explosions and less vertical mixing of the clasts, (2) generally shallower explosions, or (3) some combination of the two.

Geomorphic Evolution

The ∼6 m.y. range of ages of volcanoes in the LCVF provides an excellent opportunity to improve our understanding of landform evolution in high-altitude, arid environments. The apparent similarity of dominant eruption sizes, processes, and products across that timespan further strengthens that potential, as we can assume some similarity in the initial conditions of various surfaces and landforms. Here we give a brief, qualitative overview of time evolution of topographic features (landforms) and surfaces as a starting point for more detailed and quantitative studies in the future.

Burial and Inverted Topography

Landforms evolve slowly in the high desert environment, and can be strongly influenced by structural subsidence or uplift. A number of volcanoes in the Pancake Range occur within the Lunar Lake basin (Fig. 5). Many cones and ramparts in that domain do not have associated mapped lava fields. It is likely that in most cases this is due to subsidence and alluvial filling of the basin that buries all but the youngest lavas, and eventually the cones. Similarly, as mentioned herein, we expect that there are a number of buried episode 1 and 2 volcanoes in the valleys adjacent to the Reveille Range. Lavas sourced from vents on the upthrown side of normal faults, such as may be inferred on the northwestern edge of the Pancake Range part of the LCVF, can also result in cones that have no mapped associated lava fields because the valley-filling lavas have been detached and buried in the adjacent basin.

Volcanic fields such as the LCVF produce both high-standing, relatively easily eroded landforms such as cones and ramparts, and valley-filling, erosionally resistant deposits in the form of lava fields. A common relative age indicator in volcanic fields is the development of inverted topography, and the LCVF provides excellent examples of this. The young Marcath cone rampart overlies an escarpment, and its lavas spread upon the floor of Big Sand Springs Valley to the west (Fig. 28A). Dark Peak volcano, which is several million years old (episode 1), produced a similar cone atop a paleotopographic high and a lava field that spread westward onto a valley floor that was underlain by easily eroded volcaniclastic material. Posteruption subsidence of the surrounding basins has lowered the erosional base level, such that the former valley-filling, erosionally resistant lava field now caps a high-standing mesa (Fig. 28B). While only a remnant of the Dark Peak cone is preserved, the lava field is largely intact with apparently only minor losses of material around its margins that locally expose its internal ‘a‘ā structure.

Another form of inverted topography occurs when a pyroclastic cone or rampart is partly or fully surrounded by lavas; the lavas might be associated with the cone or rampart or might be sourced at surrounding volcanoes. Cone degradation will first produce small debris fans on top of the lava. However, continued preferential erosion of the cone or rampart can eventually produce a small basin that corresponds approximately to the original footprint of the construct, surrounded by higher lava surfaces. The result is a moat-like basin with lava forming its outer wall, and a remnant of the cone or rampart in its center; this inverted topography is most pronounced where lavas from surrounding volcanoes have partly buried an earlier construct (e.g., Figs. 28C, 28D). Eventually the only remnant of the original volcano might be a plug that corresponds to its original conduit, resulting in the “neck-in-a-basin” landform of Aranda-Gómez et al. (2010).

Cone and Rampart Degradation

LCVF cones and ramparts undergo a sequence of posteruptive surficial processes in this arid environment that is similar to those documented in the broader region (Valentine et al., 2006, 2007; Valentine and Keating, 2007). If the last eruptions at a cone cover it with nonwelded lapilli and bombs, degradation at first advances very slowly as slope creep that is manifested by the development of clast garlands (Fig. 29A). During this phase most rainwater simply infiltrates into the porous cone slopes with little or no runoff. Over a time period on the order of tens of thousands of years, eolian silt and sand accumulate on the cone surface and infiltrate the interclast pore spaces of the uppermost decimeters of lapilli and bombs, while simultaneously building sand ramps at the base of some cones (Fig. 29B). Eventually this reduces the permeability of the surficial deposits and promotes nucleation of runoff, rilling, and gullying with attendant transport of cone material to debris fans that accumulate around the base of the cone or rampart. Cones with welded agglomerate beds on their outer slopes can immediately nucleate runoff and begin the rilling phase, without the time lag to accumulate eolian sediment. The fans merge through time to form a continuous debris apron composed of sand, silt, lapilli, and bombs that thickens and grows outward as the cone erodes (Fig. 29C). Continued erosion on a million-year time frame eventually removes all but the most resistant material in the cones, such as densely welded agglomerate horizons and dikes that compose a cone’s internal plumbing (Fig. 29D), as discussed herein. In detail the rates of these cone degradation processes depend upon, in addition to the original pyroclastic facies (e.g., degree of welding), proximity to eolian sediment sources, posteruptive climate variations, and local topographic setting.

Lava Surface Evolution

Lava surfaces in the LCVF, like those in other desert settings, undergo simultaneous breaking down of roughness and topographically high points (e.g., squeeze ups, agglomerate rafts), and accumulation of eolian sand and silt. In the absence of pyroclastic deposits on top of a lava, eolian sand and silt accumulate in topographic lows, while highs gradually break down and provide clasts that become desert pavements (e.g., Valentine et al., 2007). Desert pavement in turn stabilizes eolian sediment and promotes soil development (McFadden et al., 1987). Thus lava surfaces that are initially very rough eventually evolve toward smooth surfaces (Fig. 30). During the initial tens of thousands of years after emplacement, the rates of these processes depend upon the local setting, in particular proximity to eolian sediment sources such as ephemeral streambeds and playas (cf. Figs. 30A and 23B). An additional factor is whether a lava surface starts its posteruptive evolution having been covered by tephra that provides a ready source of lapilli to form a desert pavement, beneath which eolian sediment is armored and soil formation can quickly begin (e.g., Fig. 23A; Valentine and Harrington, 2006). Carbonates precipitate in maturing soils and clay content increases with time, progressively reducing permeability and promoting runoff and development of drainage networks (Wells et al., 1985); this can be seen on many episode 3 and older lavas (Dohrenwend et al., 1987).

We have described and interpreted many individual aspects of the LCVF, from geochronology to geochemistry, spatial patterns, and physical volcanology. Here we return to the questions (posed in the Introduction) that integrate across multiple aspects of the field and have general implications for other volcanic fields.

Were the Mantle Magma Sources Homogeneous or Were They Heterogeneous in Time and Space?

The areal extent of the volcanic field–scale melting anomaly that feeds the LCVF appears to have been shrinking with time (Yogodzinski et al., 1996), although the progressive shrinking might be partly an artifact of the fact that we are only ∼300 k.y. into episode 4. Previous episodes tended to have durations of ∼1–1.5 m.y., thus it is possible that if episode 4 continues it will encompass a larger area. The apparent shrinking footprint of the volcanic field was accompanied by the north-northeast shifting of volcanism (Fig. 6; Foland and Bergman, 1992). The areas of episodes 2 and 3 partly overlap episodes 1 and 2, respectively, while episode 4, in its current level of development, is fully within the area of episode 3. Either (1) the different mantle source regions of episodes 1–3 only partially overlapped in plan view, or (2) if they did fully overlap their magmas were not evenly tapped in an areal sense to feed eruptions, or the magmas stalled in some areas in the lithosphere without erupting. We prefer the first, simpler interpretation, that the sources of episodes have been laterally offset from each other. This, along with the different source compositions for the episodes, implies lateral upper mantle heterogeneity on the spatial scales of the episodes (∼20–40 km along the axis of the volcanic field). The asthenosphere might have advected as much as a few tens of kilometers upward and/or laterally relative to the lithosphere during the ∼6 m.y. history of basaltic volcanism, distances of similar scale to the lateral shifts of episodes that could have affected the migration of volcanism. For example, different 20–40-km-scale compositional domains may have risen into the melting depth window (very loosely constrained to ∼60–130 km; Cortés et al., 2015), possibly through a shear-driven upwelling mechanism (e.g., Ballmer et al., 2015). Volcanoes within a given episode, although scattered over tens of kilometers, tapped broadly similar mantle sources. However, subtle differences between the volcanoes indicate heterogeneity within episodes, on a range of scales down to hundreds of meters (Cortés et al., 2015; Rasoazanamparany et al., 2015). This makes sense in the context of the geochemistry of the source materials, which are expected to have small-scale variability due to preferential fluid and magma pathways during metasomatization. In addition, this interpretation is consistent with the likely compositional origin of the seismically slow domain in the upper mantle beneath the LCVF (Fig. 2). The presence of an asthenospheric domain that is enriched in hydrous phases compared to surrounding upper mantle can also explain the unusual location of the volcanic field within the interior of the Basin and Range Province, although we do not speculate here about how the domain arrived in the area starting ca. 6 Ma.

Are There Relationships between Magma Sources, Clusters, and Preexisting Crustal Structures?

Volcano clusters and alignments in volcanic fields can be interpreted in two end-member ways in relation to their underlying mantle sources. The first end-member interpretation is that the clusters and alignments formed along preferred pathways for magma ascent. Implicit in this view is that melt production at depth may have been relatively uniformly dispersed beneath the volcanic field, but the distribution of volcanoes was determined by leaky crustal structures. This is perhaps the most common interpretation of clusters and alignments in the literature, especially when they are colocated with mappable structures (typically faults or joints) at the surface (e.g., Aranda-Gómez et al., 2003; Cebriá et al., 2011; Cimarelli et al., 2013). The leaky interpretation is problematic, however, because mechanical analyses show that it is more efficient for rising dikes to create their own paths than to occupy preexisting structures unless those structures are oriented within a narrow range of values relative to the principal stresses (see discussion of Relationships of Preexisting Structures, Vents, and Feeder Dikes in the Pancake Range; Delaney et al., 1986; Gaffney et al., 2007). This range of orientations narrows with depth; thus it is most likely for dikes to be captured within favorably oriented structures such as faults only in the uppermost hundreds of meters of the crust (Gaffney et al., 2007). In addition, to be captured by a preexisting fault at all, a rising dike must intersect it, or come within ∼200 m of it (Le Corvec et al., 2013b). Dikes that ascend more than ∼200 m from a fault can only follow their own paths and will not erupt along a fault. Thus, even if a fault influences the formation of a cluster, there must be a mechanism for sending dikes upward in close proximity to the fault.

A second end-member view is that clusters and alignments form over sites of relatively high melt fraction in the upper mantle; their length scales are related to the lateral length scales of these mantle source domains (so-called magmatic footprints; Valentine and Perry, 2006, 2007; Germa et al., 2013; Tadini et al., 2014). In this view the melt fraction is heterogeneous in space, and crustal structures provide secondary controls that affect rising dikes only in the uppermost kilometer or so of the crust. Our geochemical data support compositionally heterogeneous sources that could be manifested as melt-rich domains with length scales of hundreds to thousands of meters and that could be magma sources for individual volcanoes or clusters. However, we do not know of an a priori mechanism for a correlation between the map-view location of such sources in the asthenosphere and the locations of upper crustal structures at the surface. Individual volcano clusters in the LCVF had prolonged histories of several hundred thousand years or more, which has additional implications for magmatism in clusters and alignments. These time scales incorporate multiple episodes of volcanism within individual clusters. This raises questions with respect to the interpretation about source heterogeneity: if each episode was sourced from mantle with a unique geochemical signature that is internally heterogeneous (within the footprint of an episode), why would the smaller, relatively melt-rich domains that feed clusters share some locations from one episode to the next? Why would clusters be consistently elongated perpendicular to σ3 and be closely related to crustal structure?

We suggest that LCVF feeder dikes were sourced from small melt-rich domains within the larger, episode-scale melting anomalies. Individual small-volume magma batches temporarily ponded within lithospheric mantle at or just below the crust-mantle boundary (see Cortés et al., 2015). The magma source domains might have released melt batches randomly or clustered in time and space, and the overlying small reservoirs or sills were similarly distributed (Fig. 31). Ductile deformation at the ponding interval beneath fault zones might have preferentially mobilized, or activated, magmas and triggered upward dike propagation and formation of vent clusters. These dikes ascended close enough to faults to allow for dike-fault interaction in the brittle upper crust. In map view the activation of these ponded magmas was elongate along the trends of fault zones, promoting alignments of vents at the surface and elongation of clusters. Magmas that ponded beneath locations that did not undergo deformation also sometimes ascended to the surface, but less frequently, and were not necessarily closely associated with near-surface structures. In both situations some magma batches ascended straight from their deep reservoirs near the crust-mantle boundary (e.g., episode 3), while others underwent additional temporary ponding at mid-crustal levels (episode 4).

This overall conceptual model for magma ascent, volcano clustering, and alignments (Fig. 31) is a sort of hybrid of the two views of clustering and alignments summarized in the preceding two paragraphs, but rather than faults providing leaky pathways, they are associated with preferential mobilization of ponded magmas. The fault zones are long-lived features compared to the time scales of individual eruptions or magmatic episodes, and thus could promote magma mobilization and resulting volcano clustering across episodes. It is consistent with analysis of clusters and alignments in 37 volcanic fields in a range of tectonic settings that showed that “the mechanisms of magma generation are similar to those controlling its propagation through the lithosphere” (Le Corvec et al., 2013a, p. 112), with the modification that we are suggesting mechanical activation of magmas ponded in the lithosphere rather than associating the locations of clusters with deeper mantle source domains (i.e., magma generation).

What Are the Characteristics of the Transition from Crustal Magma Transport to Eruptions at the Surface?

Exposures in the LCVF indicate that each volcano was fed by a dike or dike set; individual dike thicknesses are ∼0.5–3 m. Estimated magma overpressures range from ∼5 to 20 MPa; the lower end of this range is probably the most realistic. The caveats mentioned above about overpressure calculations are very important to keep in mind: we are basing them on preserved dike widths and lengths, and the calculations do not account for the effects of multiple magma pulses, inelastic country-rock accommodation, or relaxation of country rock after overpressure has waned. It is not clear how meaningful the estimates are.

For magmatic-volatile–driven eruptions (Strombolian, Hawaiian, effusive), most eruptions began when a dike or dike segment intersected the surface, after which localized widening occurred that promoted focused flow and formation of individual vents. Eroded examples in the LCVF, although few in number, all show that this widening occurred only in the uppermost 10 m beneath a preeruptive paleosurface, similar to (but on the shallow side of) depths reported from other studies around the world. Magma and bubble ascent was mostly through narrow, planar dikes and dike sets until the last meters, where it was focused toward small (meters wide and deep) conduits. Genareau et al. (2010) suggested that progressive focusing of flow from a fissure vent to a central vent was responsible for increased magma speed, and therefore decreased gas segregation, and relatively sustained discharge with attendant tephra deposits at the Lathrop Wells cone, ∼200 km to the south of the LCVF. Some agglomerate ramparts in the northern LCVF, vented from fissures, appear to have been dominated by relatively weak Strombolian activity consistent with efficient gas segregation and formation of gas slugs in the feeder system. At Marcath, sustained fountaining and tephra deposition were related to central vent activity after an initial phase of Strombolian activity along a fissure, which produced partly welded to nonwelded agglomerate. This suggests processes similar to the Lathrop Wells case; however, it is difficult to determine whether they were consistent across other LCVF centers due to the nature (or lack) of exposures.

These essentially three-dimensional feeding systems supported Strombolian to Hawaiian explosive and effusive activity for weeks to years, consistent with the lifespans of historical scoria cone eruptive episodes (Wood, 1980; Luhr and Simkin, 1993). Such geometry, let alone more complex examples (e.g., Valentine and Krogh, 2006; Re et al., 2015, 2016; Muirhead et al., 2016), is not what we normally think of when we describe conduit models of volcanic systems. It is likely that the very localized dike to conduit transition may, for example, have an important role in gas segregation processes that determine eruption dynamics.

How Do Individual Volcanoes Grow and Feed Various Eruptive Processes?

Individual basaltic volcanoes in the LCVF had a range of eruptive processes. For example, the Marcath volcano was constructed by relatively weak (Strombolian) fountaining along a fissure, followed by sustained high fountaining (Hawaiian) similar to that observed at Kilauea Iki (Hawaii) in 1959–1960 (Richter et al., 1970; Johnson et al., 2014). Lava flows were apparently in part concurrent with high fountaining and were probably fed by it. Growth of the main cone was a competition between accumulation of coarse pyroclastic deposits and localized collapses that were rafted away by the lavas. Easy Chair also transitioned from fissure-fed activity to two main vents around which high central cones developed and from which lava was emitted (also clastogenic; Valentine and Cortés, 2013). All of this took place within the context of the feeding systems discussed here. Part of the agglomerate ramparts associated with the early eruptive phase was destroyed by subsequent phreatomagmatic, crater-excavating explosions. We assume that a small diatreme extends beneath the maar crater and completely overprints the feeder dike there, at least to the base of the inferred diatreme. Internal structures in eroded cones such as Kimana and Dark Peak (Hintz and Valentine, 2012; Harp and Valentine, 2015) show that lateral intrusions within the edifices fed bocas and lava fields. This was likely promoted by the presence of dense, degassed magma in the feeding systems. The lateral intrusions both weakened their parent cones, promoting rafting of material on the lavas during eruptions, and provided resistant material that partially preserves the cones during posteruptive degradation.

Are There Relationships between Source Characteristics (e.g., Scales of Heterogeneity), the Scales and Eruptive Styles of Individual Volcanoes, and the Time-Volume Behavior of the Volcanic Field?

Geochemical and petrologic data have shown that the sources and, to some extent, ascent histories in terms of temporary crustal ponding, of LCVF magmas varied between the four episodes of magmatism. In some volcanic fields, partial melting is crudely related to volumes of individual eruptions in that larger degrees of melting (e.g., indicated by trace elements) produced larger eruption volumes (e.g., southwest Nevada volcanic field; Valentine and Perry, 2007). We do not have enough information at the LCVF to constrain whether the degree of partial melting has varied systematically, but, to the extent that we can determine given that older volcanoes are partly eroded and/or buried, the range of eruption volumes appears to have been similar across the episodes. Similarly, eruption styles have been similar with the exception of the three small evolved centers during episode 2 that formed lava domes and small coulees. Thus in the case of the LCVF we do not see relationships between eruptive scales and styles, and the magma source and ascent characteristics.

Total Eruptive Volume and Time-Volume Relationships

A key parameter in defining the long-term dynamics of volcanic systems is their time-volume behavior, normally presented as a plot of cumulative eruptive volume versus time for a volcanic field (e.g., Bacon, 1982; Kuntz et al., 1986) that in turn requires well-constrained eruption ages and reliable eruption volume estimates. Existing age data (Table 1) cover only 22%–25% of the known volcanoes in the LCVF. Volumes can only be reliably constrained for ∼10 individual volcanoes; others have been partly buried by later volcanic products, are partly buried by basin-filling sediments, or, especially in the Reveille Range, have been significantly eroded. This is compounded by the possible occurrence of an unknown number of completely buried volcanoes (e.g., in the Railroad, Kawich, and Big Sand Springs Valleys, and in the Lunar Lake basin). Buried volcanoes also reduce the percentage of centers for which eruption ages have been determined, by some unknown amount.

Previous publications that compiled long-term eruptive volume fluxes for volcanic systems around the world listed the LCVF as having a total volume of 100 km3 (Crisp, 1984; White et al., 2006). This value is based on the total area of the volcanic field (∼400 km2) and an assumed average of the volcanic products of 250 m (Bergman, 1982); no basis is provided for that assumed thickness. Given that the Marcath DRE eruption volume totals ∼1.1 × 108 m3 (∼0.1 km3), and that it appears to be representative of larger eruptions in the volcanic field, which contains ∼142–161 mappable vents, the published total volume might be too large by a factor of ∼5 (but this does not include possible buried volcanoes).

Thus while we know that the LCVF must have a long-term eruptive flux that is between the two type end members of volcanic fields (southwest Nevada volcanic field and eastern Snake River Plain; Kuntz et al., 1986; Valentine and Perry, 2007), the age and volume uncertainties are such that we cannot quantify time-volume behavior beyond estimating that the long-term flux was probably between ∼3 and 5 km3/m.y. averaged over the lifetime of the field. Tadini et al. (2014) also noted this problem, but suggested that isolated volcanoes around the periphery of the LCVF (Pancake Range) might be driven by processes akin to tectonically controlled fields that have very low long-term magma fluxes.

Comparison with Nearby Volcanic Fields

Most studies of volcanic fields in the interior southwestern United States have focused on petrogenesis; only a few have also been the subject of physical volcanological studies. Here we very briefly compare magma source, ascent, and volcanology (where it has been studied) between the LCVF and volcanic fields in the region, keeping in mind that the LCVF is relatively isolated even from those fields.

The Big Pine volcanic field (Fig. 1) is located ∼200 km west-southwest of the LCVF, at the boundary between the Basin and Range Province and the Sierra Nevada. The Big Pine field includes ∼24 basaltic volcanoes with abundant mantle xenoliths, and one silicic center, with ages ranging from 1.2 to 0.035 Ma (Blondes et al., 2008). Many of the volcanoes occupy young normal faults. The basaltic volcanoes comprise scoria cones with lava fields and are similar in size to those of the LCVF. Gazel et al. (2012) concluded that the asthenospheric melting interval has shifted upward (from 60–80 km to 30–55 km) during volcanism in concert with the shallowing asthenosphere-lithosphere boundary (defined seismically); younger (post–500 ka) magmas also bear evidence of having interacted with lithospheric mantle. Gazel et al. (2012) suggested that relationships between magmatism, vents, and structures are related to lithospheric deformation (faulting at the surface), perhaps similar to our interpretation of at least some aspects of the LCVF (Fig. 31). Blondes et al. (2008) concluded from isotopic studies that the source mantle was heterogeneous. It is interesting that one of the eruptions captured the full range of source characteristics that are represented in the field as a whole; this might reflect vertical heterogeneity that could be variably captured as magmas collected within ∼20-km-thick melting intervals.

The Buffalo Valley volcanic field (Fig. 1; Cousens et al., 2013), like the LCVF ∼230 km to the south-southeast, is located in the interior of the Basin and Range Province rather than near its edge. It consists of ∼15–20 basaltic scoria cones and lava fields with compositions from alkali-rich basalt to trachybasalt. The field had two episodes of magmatism, both of which were sourced in asthenospheric mantle. The Pliocene episode (4.02–2.75 Ma) erupted magmas that originated at relatively shallow depths in the asthenosphere, indicated by the absence of a garnet signature in trace element compositions. Quaternary episode (1.99–0.07 Ma) products, however, have a compositional signature consistent with garnet in the source mantle. Thus melting appears to have deepened with time (Cousens et al., 2013). Volcanoes of Quaternary age are clearly colocated with faults that bound the Fish Creek Mountains, potentially supporting our conceptual model linking crustal faulting with melt collection at depth.

The Hurricane volcanic field is located ∼250 km east-southeast of the LCVF, just in the Colorado Plateau Province (Fig. 1; Smith et al., 1999). Colorado Plateau lithosphere and crust is much thicker (>120 km and 40–50 km, respectively; e.g., West et al., 2004) than in most of the Basin and Range because it has not been extended and thinned. The small field contains ∼10 scoria cone volcanoes with attendant lava fields, with individual volumes of 0.05–0.3 km3. Pyroclastic deposits are consistent with Hawaiian and Strombolian explosive styles, accompanied by lava effusion. Compositions are basanites and alkali basalts. Age dates range between 353 and 258 ka (Smith et al., 1999), although eroded remnants suggest some older volcanism as well. Hurricane magmas are inferred to have originated by partial melting of heterogeneous lithospheric mantle (Smith et al., 1999), differently from the LCVF asthenospheric sources. The lithospheric mantle source was compositionally heterogeneous at scales of ∼1 km (or less) in plan view. The volcanoes are closely associated with faults and one of them is cut by the Hurricane fault, a major normal fault associated with the transition from Colorado Plateau to Basin and Range tectonics to the west.

The southwest Nevada volcanic field (Fig. 1) includes 17 Pliocene–Quaternary monogenetic volcanoes (Valentine and Perry, 2007). Compositions include alkali basalts to trachyandesites, and are consistent with most of the volcanoes having been fed by magmas that ascended relatively directly from their source depths to the surface (one exception has evidence for storage and fractionation in the crust prior to eruption). Trace elements define two major magmatic episodes. The first one produced three volcanoes of relatively large volumes (∼2.5, 1.2, and 0.6 km3), and several smaller volume centers that are now buried beneath basin fill. Ages for this first episode are between 4.63 and 3 Ma. Younger volcanoes are all smaller than 0.1 km3 with 2 exceptions at 0.84 and 0.12 km3. Trace elements indicate that the older episode was related to a larger degree of partial melting than the younger one. All of the magmas were sourced in enriched, old (Precambrian), and compositionally heterogeneous lithospheric mantle, which extends to depths of ∼125 km (Schulte-Pelkum et al., 2011). Note that this area, referred to as the central Basin and Range, began extending significantly later than most of the province, such that its crust and lithosphere has not thinned as much. Apparently all of the feeder dikes were captured by north-south–trending, Miocene-age faults as they approached the surface. This volcanic field was proposed as the type example of a tectonically controlled field, where melt collection in relatively melt richer domains of lithospheric mantle was dominantly in response to regional tectonic deformation (Valentine and Perry, 2007). The Quaternary long-term eruptive magma flux is extremely small, ∼0.5 km3/m.y. Eruption styles included Hawaiian, Strombolian, and violent Strombolian explosive activity, as well as lava effusion.

Of the volcanic fields summarized here, two (Big Pine and Hurricane) are at or near the margins of the Basin and Range Province. The Buffalo Valley and southwest Nevada fields are in the interior of the province. Like the LCVF, all of these fields were fed by magmas sourced in heterogeneous upper mantle, be it asthenospheric or lithospheric, with heterogeneities extending to small scales such that closely spaced volcanoes have slightly different source compositions. As shown in McGee and Smith (2016), small-scale compositional heterogeneity of source mantle appears to be a common characteristic of small-volume volcanoes fed by low degrees of partial melting. These heterogeneities produce small-volume domains of partial melt that in turn feed small-volume volcanoes. In many cases the magmas ascend relatively quickly; even if they temporarily pond in the crust, their small thermal budgets work against significant assimilation of host rock (see Cortés et al., 2015). In this way the volcanoes are able to preserve the subtle source variability. All of the fields exhibit close links between surface faulting and volcano location, consistent with our conceptual model (Fig. 31).

The LCVF is similar to many volcanic fields in the interior western United States. It is dominated by monogenetic basaltic volcanoes (in this case with a very minor component of lava domes and flows with evolved compositions). The volcanoes number 142–161 and span ∼6 m.y. in age, with the most recent only ca. 35 ka, indicating that the field is likely to produce future eruptions. The most common volcanic landforms are scoria cones, agglomerate ramparts, and ‘a‘ā lava fields. The volcanoes were fed at depth by dikes that occur in en echelon sets and that preserve evidence of multiple pulses of magma. The dikes locally flared in the upper ∼10 m of the crust to form shallow conduits that fed eruptions. Volumes of individual eruptions were typically ∼0.1 km3 and smaller. Eruptive styles were dominantly Strombolian to Hawaiian, the latter producing tephra fallout blankets, with effusive activity (although many lavas were likely clastogenic and associated with lava fountains). Phreatomagmatic maar volcanoes compose a small percentage of the landform types; we are unable to identify any clear hydrologic or climatic driver for the phreatomagmatic activity, suggesting that intrinsic factors such as magma flux played an important role.

Volcanism in the LCVF occurred in four magmatic episodes; each successive episode shifted northward but partly overlapped the area of its predecessor. The areas encompassed by the episodes appear to be progressively shrinking, although this may be partly an artifact of episode 4 having not yet fully unfolded. Asthenospheric mantle was the source throughout LCVF volcanism, but each episode was derived from a specific source composition. Thus the episodes record large-scale (length scales of tens of kilometers) heterogeneities in the upper mantle. The longevity of the volcanic field is such that advection of different compositional domains into the melting region might have been important. Geochemical evidence suggests subtle source heterogeneity at the scales of hundreds of meters to kilometers within each episode-scale area of activity; episode 1 magmas also underwent contamination by Paleozoic carbonates during temporary ponding in the crust. The LCVF, like most volcanic fields, exhibits clustering of vents that we relate to domains of elevated partial melt content in the asthenosphere, consistent with the heterogeneous source compositions. Magmas underwent temporary ponding near the crust-mantle boundary, and we infer that these were preferentially mobilized by ductile deformation of host rocks; this is in turn represented by faulting at the surface. Thus clusters could persist from one episode to the next, and their locations are closely related to crustal structures. Clusters result in overlapping and colocated monogenetic volcanoes that are separated by variable amounts of time, but that record different mantle sources (between episodes) without having shared, sustained crustal reservoirs. The total volume and time-volume behavior of the LCVF cannot be determined by surface observations due to erosion and to burial by basin-fill sediments and subsequent eruptive products. However, previous estimates of a total volume of 100 km3 are likely too high by a factor of ∼5.

Monogenetic volcanic fields are areas of active investigation around the world. Linking physical volcanological, structural, geochemical, and petrologic data are key in advancing our understanding of these systems. High-resolution geophysical studies that identify crustal structures and potential kilometer-scale domains of elevated melt, the inferred sources of individual volcanoes and clusters, could provide important tests of the conceptual model we have suggested for relating surface expressions of volcanism with deep magma sources, as will numerical modeling. Such work is needed to improve our understanding of the full spectrum of volcanism on Earth and, by analogy, the other terrestrial planets.

This work was supported by U.S. National Science Foundation collaborative grants EAR-1016100 (Valentine), 1016042 (Widom), and 1016104 (Smith). We thank Duane Champion for making us aware of certain age dates provided by Brent Turrin, and Brandon Schmandt for providing the tomographic image. Dave Kuentz provided valuable support in sampling and geochemical analyses. We also thank the following students who assisted in the field: Nicole Shufelt, Peter Johnson, Amanda Hintz, Jamal Amin, Alessandro Tadini, Matthew Sweeney, Dayana Schonwalder, Marissa Mnich, Christine Sumner, Patrick Whelley, Mai Sas, Amber Ciravolo, and Hugo Belmontes. We thank Jan Lindsay, Marco Brenna, and Nicolas Le Corvec for their helpful reviews of the manuscript.

1Supplemental Table. Analytical data and calculations related to the 10Be age determinations at Lunar Crater maar. Please visit http://dx.doi.org/10.1130/GES01428.S1 or the full-text article on www.gsapubs.org to view the Supplemental Table.