In ancient or partially eroded arc sections, a protracted history of tectonism and deformation makes interpretation of local volcanic-plutonic relationships challenging. The fragmentary preservation of volcanic rocks relative to the extensive plutonic record in upper-crustal arc sections also suggests that a broader-scale approach that includes volcanic-hypabyssal-plutonic “fields” is useful. In this context, studies of hypabyssal intrusions emplaced at the intersection of volcanic and plutonic fields provide additional physical and chemical constraints on shallow-level magmatic processes. New mapping, U-Pb zircon geochronology, and geochemistry at Tioga Pass, in the central Sierra Nevada arc section, document the physical and chemical evolution of the Tioga Pass hypabyssal complex, a ca. 100 Ma system that includes an intrusive dacite-rhyolite porphyry unit and comagmatic Tioga Lake quartz monzodiorite. We interpret these units as a Cretaceous subvolcanic magma feeder system intruding a package of tectonically displaced Triassic and Jurassic volcanic and sedimentary rocks, rather than the previous interpretation of a Triassic caldera. The Tioga Pass magmatic system is a well-exposed example of a hypabyssal complex with meso- to micro-scale structures that are consistent with rapid cooling and emplacement between 0–6 km depth and compositions suggestive of extensive fractionation of largely mantle-derived magma. The Tioga Pass porphyry unit is one of many hypabyssal intrusions scattered along a ~50-kilometer-wide belt of the east-central Sierra Nevada that are spatially associated with coeval volcanic and plutonic rocks due to tectonic downward transfer of arc crust. They provide a valuable perspective of shallow magmatic processes that may be used to test upper-crustal plutonic-volcanic links in tectonically reorganized arc sections.

Studying the volcanic-plutonic connection in ancient and partially eroded arc sections presents some unique challenges, particularly from a temporal and geochemical perspective. Exposed sections, often biased towards deeper plutonic levels, may provide a somewhat time-averaged (yet longer-lived) record of arc activity, while volcanic units are susceptible to the effects of hydrothermal alteration and metamorphism (e.g., [1, 2]). Syn- to post-arc deformation and tectonism overprints the magmatic record and can reorganize and displace local and regional stratigraphy, masking true field relationships.

However, tectonic shuffling is also one of the primary reasons to study these crustal sections. Tectonism often results in the exposure of multiple crustal levels, providing information on the physical and spatial relationships of both plutons and contemporaneous volcanic rocks and can additionally reveal the physical structure (e.g., shape and size) of the subvolcanic magma plumbing system. Restoring local- to regional-scale structures allows for inclusion and evaluation of magmatic-tectonic relationships through mapped spatial and temporal associations of structures and magmatic rocks. Together, these factors suggest that ancient arc crustal sections provide a different perspective to understanding plutonic-volcanic relationships compared to modern systems.

In upper-crustal arc sections, hypabyssal intrusive complexes are important features as they are exposed in both volcanic- and plutonic-dominated arc domains (e.g., [3, 4]). They represent a physical connection between these environments, which is key to reconciling existing plutonic-volcanic observations and models (e.g., [58]). As hypabyssal complexes are emplaced at shallow levels, they have the potential to provide a snapshot of magmatic processes that relate to both plutons and volcanic rocks, such as crystal accumulation, melt extraction, and eruption (e.g., [9]).

In the Sierra Nevada arc section, several studies have drawn attention to the Mesozoic volcanic-plutonic connection, either by regional-scale comparisons between mid- to upper-crustal plutons and volcanic stratigraphic sections (e.g., [1016]) or by the study of local volcanic features, such as calderas (e.g., [1719]). More information on shallow levels of the plumbing system, accessed through the hypabyssal record, improves our understanding of the physical development and evolution of upper-crustal magma plumbing systems and has the potential to link observations made across wide spatial scales, from a single intrusion to a regional-scale magmatic “field.”

At Tioga Pass, in the eastern-central Sierra Nevada, we reinterpret a proposed Triassic caldera [19] as a Late Cretaceous hypabyssal complex. We present new field mapping, U-Pb zircon ages, and whole-rock element and isotopic geochemical data from Tioga Pass to characterize magmatic and tectonic histories from Triassic to Late Cretaceous time. We demonstrate that both the Tioga Lake quartz monzonite and intrusive dacite-rhyolite porphyry are ca. 100 Ma and represent comagmatic parts of a punched laccolith-shaped body (steep-sided and flat-topped, with vertical walls discordant to host rock structure; [20]) that intruded blocks of faulted and tilted Triassic and Jurassic volcanic and sedimentary strata. The Tioga Pass system shares compositional and structural affinity with other Late Cretaceous intrusive and volcanic rocks in the east-central Sierra Nevada region and is one of several hypabyssal intrusions found across a ~50-kilometer-wide belt in the central Sierra Nevada arc section. These intrusions are key structural markers and may represent feeders to volcanic eruptions or stalled late melts from plutons.

2.1. Regional Background

The Mesozoic Sierra Nevada Batholith is an upper- to lower-crustal composite arc section, dominated by magnesian, calc-alkaline, metaluminous tonalite, hornblende-biotite granodiorite, and locally peraluminous granite. In the central Sierra Nevada, emplacement depths of Triassic-Cretaceous intrusive suites range from 6–10 km based on Al-in-hornblende barometry [2123]. The 120–85 Ma Cretaceous flare-up resulted in the emplacement of voluminous upper-crustal intrusive complexes (e.g., [13, 24]). Cretaceous volcanic sections are locally exposed in several pendants of the Sierra Nevada Batholith, including the Ritter Range, Saddlebag Lake, and Piute pendants (Figure 1) [2529], although most of the Cretaceous volcanic cover is interpreted to have been eroded. Calculated magma addition rates, using Cretaceous plutonic and volcanic rock areas, resulted in plutonic-volcanic ratios between 20 : 1 to 30 : 1 during flare-up magmatism [24].

2.2. Geology of the Saddlebag Lake Pendant

The Saddlebag Lake pendant exposes Paleozoic deep marine framework rocks and a sequence of Mesozoic arc volcanic and sedimentary packages [11, 2831] (Figures 1 and 2). It is intruded by the Triassic Scheelite Intrusive Suite in the east and the Cretaceous (95-85 Ma) Tuolumne Intrusive Complex to the west [10, 13, 32]. Regional metamorphism pervasively altered rocks to greenschist facies [31, 33]; the prefix “meta-” is hereafter omitted from rock descriptions where the protolith is known. Unconformable contacts between westward younging major stratigraphic packages were subsequently reactivated by strike-slip faulting [28, 30, 34]. Tioga Pass is located at the southeastern end of the Saddlebag Lake pendant (Figures 1 and 2), marking the transition of the Saddlebag Lake pendant into the northern Ritter Range pendant.

The distinctive lithologic units exposed west of Tioga Lake were previously mapped and studied by Kistler [35], Brook [31], and Schweickert and Lahren [19] (Figure 2). Interbedded sandstone and tuff units, discordant to regionally NW-striking volcanic units, could not be stratigraphically correlated to nearby Saddlebag Lake or northern Ritter Range pendants. Geologic interpretations from field relations varied widely, from hypabyssal pluton emplacement along a normal fault, a shallow volcanic vent structure, and a caldera structure [19, 31]. Schweickert and Lahren [19] interpreted a Triassic age for the proposed caldera collapse event, caldera fill sequence, and resurgent intrusions based on contact relationships between an intracaldera vent structure and the dated 222 Ma tuff of Saddlebag Lake (TRsl; Figure 2).

The Dana sequence is a ~1 km thick layered package of volcanic and sedimentary units located southeast of Tioga Pass, in the vicinity of Mt. Dana, and has been linked to units at Tioga Lake (Figure 2; [19, 36, 37]). Russell [38] found basal shallow marine interbedded volcanic and sedimentary deposits, grading upwards into terrestrial ash-flow tuff, cross-bedded sandstone, and conglomerate. Intrusive units mapped at Tioga Pass, such as the quartz monzodiorite and dacite porphyry unit, were also found at Mt. Dana [35, 36]. Several intrusions and faults mask original contacts of the Dana sequence with surrounding units. Kistler [35] and Russell [38] assigned a lower Jurassic age to the Dana sequence, based on lithologic correlations to the Dunlop Formation (Nevada) and nearby Ritter Range pendant. Schweickert and Lahren [19] reinterpreted this section as part of the Triassic caldera-fill package.

We remapped the Tioga Pass area at 1:10,000 scales, compiled map data from the literature [19, 35, 36, 38] and collected U-Pb zircon laser ablation-inductively coupled plasma mass spectrometry (LA-ICP-MS) data to temporally constrain regional stratigraphy. We dated tuff units at the base and the top of volcanic packages, defined by regional unconformities or faults, to bracket the range of possible ages for coherent volcanic sequences, and dated sedimentary units to estimate maximum depositional ages. We dated cross-cutting intrusions to provide crystallization age estimates. Detailed 1:1,000 scale mapping, petrography, and whole-rock geochemistry focused on the Tioga Pass intrusive units. Stereographic projections were created using the program Stereonet 9.8 [39].

3.1. U-Pb Zircon Geochronology

LA-ICP-MS U-Pb zircon analyses were performed at the University of Arizona Laserchron Center using conventional lab procedures (e.g., [4042]). Full analytical methods and data tables are included in Supplementary Materials (available here). Sri Lanka, FC-1, and R33 zircon grains were used as primary age standards [41]. For igneous samples, 20–30 zircons were analyzed where possible, and a weighted mean age calculated, reported at the 2σ level. We included in our weighted mean calculation only concordant zircons with U/Th < 10 as defined by Gehrels et al. [42] for igneous zircons. This excluded a small number of zircons (n=17 grains) from the total sample set of igneous samples (n=199 grains). For detrital zircon ages from sedimentary units approximately 100 zircons were analyzed, where possible. The location of the youngest zircon peak, made of at least 3 zircons with overlapping uncertainty, is interpreted as the maximum depositional age of the sample [43]. For published LA-ICP-MS ages where sample data are available (e.g., [28]), we reinterpreted age estimates using the above protocol in order to improve consistency across samples.

3.2. Whole-Rock Geochemistry

Sixteen samples were collected for whole-rock elemental geochemistry to characterize the porphyry (n=10) and the quartz monzodiorite unit (n=6). Fourteen samples were analyzed by X-ray fluorescence (XRF) at Pomona College using lab procedures outlined in Lackey et al. [44]. Five samples were analyzed by inductively coupled plasma-mass spectrometry (ICP-MS) for major oxides and/or trace elements at Activation Laboratories, Ancaster, under the “4LithoResearch” package. Three samples were analyzed for trace elements by laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) at Texas Tech University. Two samples were selected for whole-rock isotopic analysis performed at the University of Arizona, following methods in Otamendi et al. [45]. The standard error for isotope samples is between 0.0007%-0.0008% for 87Sr/86Sr and 0.0009-0.0018% for 143Nd/144Nd, and between 0.003-0.008% for Pb isotopes. Detailed analytical methods are included in Supplementary Materials.

4.1. Stratigraphy of the Study Area

The wide-ranging geologic interpretations at Tioga Pass highlight the importance of combining a detailed stratigraphic and field-based record with modern geochronology to resolve outstanding uncertainty. Stratigraphic, field, and geochronologic relationships are described below for the intrusive complex and the surrounding host rocks. We use these to establish which units are spatially and/or temporally associated, restore tectonic effects, and examine the contact relationships between units. We use previously defined regional terms, such as the Triassic “Koip sequence” and Triassic-Jurassic “Dana sequence,” to describe and group temporally and spatially related strata [13, 35].

The upper Paleozoic deep marine section forms the eastern boundary of the Tioga Pass complex, where fine-grained sandstones are interbedded with chert and phyllite (Pzs; Figures 2 and 3(a)–3(c)). A detrital zircon maximum depositional age of 435 Ma was attained from a fine-grained sandstone sample northwest of Tioga Lake (sample TP-8; Table 1, Figures 2, 4, and 5). The lithology and age of overlying Triassic Koip sequence metavolcanic rocks in the Saddlebag Lake pendant have been described by Brook [31], Barth et al. [10], Paterson and Memeti [28], and Cao et al. [34]. At Tioga Pass, the rhyolite tuff of Black Mountain (TRbm; Figures 2, 3(a), and 4) defines the base of the Koip sequence and yielded an age of 229.6±1.8Ma (sample TP-14, Figure 5). This age estimate is overlapping with a U-Pb zircon sensitive high-resolution ion microprobe (SHRIMP) age from Barth et al. [10] of 232±2Ma. The overlying Cooney Lake conglomerate (TRcl; Figures 2, 3(a) and 3(b), and 4) contains meter-wide pebbly layers interbedded with coarse sandstone. Detrital zircon ages for Cooney Lake conglomerate samples form a narrow unimodal distribution with minimum peaks at 219-220 Ma (samples KA-9, TP-7, TP-1; Figure 5). The overlying 224±1Ma rhyolite tuff of Saddlebag Lake (TRsl; Figures 2, 3(a) and 3(b), and 4) [10] separates the Cooney Lake conglomerate from a package of crystal-rich flows of andesite-dacite composition (TRan; Figures 2, 3(a), and 4). The uppermost section of the Koip sequence (TRuv; Figures 2, 3(a) and 3(b), and 4) includes layers of clastic rhyolite-dacite tuff, crystal-rich andesite flows, and volcanic breccia, thinly interbedded with volcaniclastic conglomerate and sandstone. One sample of a rhyolite breccia east of the Maul Lake fault yielded an age of 220.0±2.2Ma and constrains the top of the preserved section of Koip sequence metavolcanic rock units at Tioga Pass (Figures 4 and 5; sample KA-6).

The ten-meter-wide, brittle Maul Lake fault (NW corner of Figure 2) juxtaposes Koip sequence units (TRuv) with the ca. 95 Ma Kuna Crest granodiorite (Kkc; Figures 2, 3(a), and 4; [32]). Jurassic metasedimentary rocks and Cretaceous volcanic rocks form part of the stratigraphy further north in the Saddlebag Lake pendant (e.g., Sawmill Canyon, [28]; Virginia Canyon, [34]) (Figures 1 and 4), and south in the Ritter Range pendant [46]; however, they are not exposed at Tioga Pass. They were likely entirely removed during intrusion of the Tuolumne Intrusive Complex (Figure 4) [28, 47].

4.2. Tioga Lake and Dana Sequence Sections

4.2.1. Tioga Lake Section

A coherent, NW-striking, and SW- to NE-dipping section of volcanic and volcaniclastic strata forms an approximately 800 x 400 m block exposed west of Tioga Lake (TRvs, TRan, and TRrh; Figure 2). This volcanic block is intruded on all exposed sides by a dacite-rhyolite porphyry unit and quartz monzodiorite unit (Kdp, Ktm; Figures 2, 3(b), and 4). Units include volcaniclastic sandstone, pebbly monomict conglomerate, and andesite and rhyolite composition lava flows and tuff. The section has a total thickness of ~850 m. Meter-wide dacite porphyry dikes intrude these units at a high angle. A maximum depositional age of 217 Ma from a volcaniclastic sandstone in the intra-porphyry package (sample TCL 13-5; Figures 4 and 5) is similar to the youngest Koip sequence units in the Saddlebag Lake pendant.

4.2.2. Dana Sequence Section

The Dana sequence is a fault-bounded package of interbedded metavolcanic and metasedimentary units that is restricted to the area around Mt. Dana. It is surrounded by the Paleozoic chert-argillite unit. The package grades from mudstones and intermediate to felsic tuffs and flows (TRuv and TRtm; Figures 2 and 3(c)) into a ~200-meter-wide rhyolite tuff (TRrh; Figures 2 and 3(c)) and overlying volcaniclastic sandstone unit (TRvs; Figures 2 and 3(c)), consistent with the descriptions from previous studies [36, 38].

An andesite tuff from the exposed base of the Dana sequence yielded a weighted mean age of 212.8±6.4Ma; however, this sample yielded only 7 zircons, 6 of which were concordant (sample MD-8C; Figures 2, 4, and 5). A rhyolite flow ~200 m below the summit of Mt. Dana yielded an age of 195.1±2.1Ma (sample L101; Figures 2, 4, and 5). A fault-bounded rhyolite tuff yielded an age of 221.7±2.5Ma and could be a part of either the Koip or Dana sequence (sample MD-9; Figures 2, 3(c), 4, and 5).

4.3. Tioga Pass Intrusive Units

4.3.1. Quartz Monzodiorite

The Tioga Lake quartz monzodiorite (Ktm; Figure 2) crops out north, west, and east of Tioga Lake and has an elongate shape in map view (Figure 2). It was previously named the granodiorite of Tioga Lake ([13, 19], “unassigned granitic rocks” by [35]) and ranges in composition from quartz monzodiorite to granodiorite [48]. It is an equigranular biotite- (~5-10%) and hornblende-rich (~10%), plagioclase dominated (60–70%) unit with quartz (15-30%); grain size varies between 0.5 and 4 mm. The unit is finer-grained and porphyritic towards the western intrusive contact with the porphyry unit (Figures 2 and 6). The quartz monzodiorite typically grades into the porphyry intrusion close to the contact and contains abundant rounded to subangular porphyritic enclaves throughout (Figure 7(a)).

A ~50-meter-wide mutually intrusive contact zone separates the quartz monzodiorite unit from the adjacent porphyry; however, in some outcrops, the units have a sharp intrusive boundary (Figures 6 and 7(b)). The quartz monzodiorite intrudes Paleozoic metasediments along a near-vertical contact on northern and eastern exposures and contains stoped blocks of metasediments and metavolcanic tuff from 1 to 100-meter scale (Figures 3(a) and 3(b)). Close to the contact with the Paleozoic quartzite, plutonic samples contain clots of biotite.

Petrography. Plagioclase is the dominant mineral in the hypidiomorphic quartz monzodiorite unit. Distinct, resorbed cores, and well-defined oscillatory rim zoning are observed throughout the pluton (Figure 7(c)). Plagioclase crystals have euhedral to subhedral habit, show frequent plagioclase-plagioclase contacts, and are almost always sutured with other plagioclase crystals (Figures 7(c) and 7(d)). Euhedral to subhedral hornblende crystals contain inclusions of magnetite. Biotite is interstitial, <0.5 mm in size, with minor replacement by chlorite. Quartz occupies interstitial spaces between feldspar crystal clusters and shows undulose extinction and recrystallization indicating low-temperature crystal plastic deformation.

4.3.2. Porphyry Intrusion

The dacite-rhyolite porphyry (Kdp; Figure 2) crops out as several large NW-striking elongate intrusive bodies across the Tioga Pass area. The largest continuous body is the exposed west of Tioga Lake. The porphyry has an aphanitic quartz-biotite groundmass and euhedral plagioclase and alkali feldspar (~75–80%), biotite (~10–15%), and hornblende (~5%) phenocrysts. The unit grades from a low silica dacite to high silica rhyolite to the west.

The porphyry is internally structurally complex at the micro- and meso-scale (Figure 6). Quartzite, pelitic schist, rhyolite tuff, basaltic-andesite, and andesite lava flows make up displaced blocks (up to tens of meters scale) that are incorporated into, and intruded by, the porphyry. Meter-scale rounded inclusions of the quartz monzodiorite unit are also observed within the porphyry intrusion (Figure 6). The abundance of phenocrysts varies from <5% to 50%, while phenocryst size ranges from 1-4 mm. Crystal-rich (45-50%) dacite is observed near the contact with the Tioga Lake quartz monzodiorite, as well as in dispersed zones throughout the porphyry unit (Figure 7(e)). Crystal-poor rhyodacite to high-silica rhyolite (<10% phenocrysts) forms the western part of the intrusion, with extensive mingling zones of intermediate composition magma and flow-banded rhyolite (Figures 6 and 7(f)). Grain-size generally decreases to the west. Massive and layered flow-banded zones of the porphyry unit are interpreted to be devitrified obsidian, from the outcrop and microscale observations of phenocryst and matrix grain size, composition, and local extent of the layering. Volcanic clasts are prominent in some layers of the unit (Figures 6 and 7(g)). Clasts of crystal-rich dacite are observed within the rhyolite portions of the porphyry (Figure 7(h)).

The porphyry has a discordant intrusive contact with Paleozoic metasediments and Lower Koip sequence units at Tioga Lake. In the east, the porphyry intrudes the Dana sequence and Paleozoic metasediments. Contacts are sharp and near-vertical (Figures 2 and 3(a) and 3(b)). The western contact of the porphyry with Saddlebag Lake pendant is offset<100m by the Gaylor Peak fault. The contact between the porphyry and quartz monzodiorite is typically gradational, except where the porphyry cross-cuts the quartz monzodiorite unit as dikes or sheets.

Petrography. In thin section, the porphyry has a microcrystalline groundmass (Figures 7(i)7 to (l)). Most groundmass appears devitrified to quartz, feldspar, mica, and oxides, although glass remains intact as inclusions or as amorphous aggregates 0.2–0.5 mm (Figure 7(l)). In flow-banded layers, the groundmass wraps around aligned phenocrysts forming a magmatic foliation (Figures 7(i) and 7(l)). Hornblende phenocrysts, observed in samples of intermediate composition, are euhedral-subhedral in shape and 1–3 mm in size. Some hornblende phenocrysts are replaced by biotite at their rims. Biotite phenocrysts are often euhedral and tabular mm-size grains, in some cases deformed (Figure 7(j)), while groundmass biotite is microcrystalline (up to 0.1 mm). Hornblende and biotite decrease in abundance from the contact zone towards the western exposed edge of the porphyry unit.

Plagioclase phenocrysts are ubiquitous in the porphyry unit (Figures 7(i)–7(k)). Plagioclase is the most abundant feldspar in the section of the porphyry unit closest to the contact, at approximately a 3 : 1 ratio of plagioclase to alkali feldspar. Westwards, the proportion of alkali feldspar phenocrysts increases to a 1 : 4 ratio of plagioclase to alkali feldspar (sample KA13; Figure 6). Alkali feldspar phenocrysts in sample KA13, a flow-banded rhyolite portion of the porphyry (Figures 6 and 7(i)) are ~1 mm in size and are coarser in the dark flow bands relative to the light bands. Feldspars have oscillatory zoning preserved throughout the unit, and there is a population of grains identified by truncated zones and contact melting points between two grains. Some phenocrysts are fractured and infilled with melt (Figure 7(l)). Plagioclase phenocrysts are often arranged in clusters (Figure 7(k)). Signs of local alteration by fluids include quartz-epidote veins and chlorite lenses in the groundmass (Figure 7(l)).

4.3.3. U-Pb Zircon Ages of Intrusive Units

The Tioga Lake quartz monzodiorite yielded a U-Pb zircon weighted mean age of 101.0±1.5Ma (sample TIOGA; Figure 8(a)), older than the 96 Ma U-Pb zircon age (n=3 samples) reported by Schweickert and Lahren [19]. The porphyry unit was dated at the northern and southern extents of the unit and yielded ages of 100.0±1.3Ma and 99.5±1.6Ma (samples KA17 and TP17), indistinguishable in age within uncertainty from the Tioga Lake quartz monzodiorite (Figures 8(b) and 8(c)). Low MSWD values for these samples indicate the overestimation of uncertainty [49] (Table 1). The 99.1±1.4Ma age estimate of an outcropping granite unit (sample 3-3; Kg on Figure 2) at the Kuna Crest granodiorite margin also overlaps in age with the quartz monzodiorite and dacite porphyry unit. In this sample, the age of the youngest peak, at 98-99 Ma, is interpreted as the unit age, as the high MSWD value indicates that the current weighted-mean model does not suitably fit the data [49, 50] (Figure 8(d), Table 1). This age estimate for sample 3-3 is older than the 93-95 Ma isotope-dilution thermal ionization mass spectrometry (ID-TIMS) age estimates for the Kuna Crest margin ~10 km south of Tioga Pass [32]. In summary, the granite is contemporaneous with the dacite porphyry and quartz monzodiorite, but the petrogenetic relationship of the granite to these units remains uncertain.

4.4. Structural Features of the Study Area

Structural data from Mesozoic volcanic and sedimentary rocks and Cretaceous intrusive rocks are used to evaluate the timing and impacts of local tectonism in deforming and displacing units related to the intrusions and to provide information on the three-dimensional structure of the intrusive complex.

A bedding-parallel foliation and cleavage is well-exposed in Paleozoic and Mesozoic rocks, with an average foliation of strike/dip = 160/85 in the mapped area (Figure 2). Mineral stretching lineation (measured in biotite and quartz) plunges between 70–90°. The Triassic Tioga Lake block (Figures 2, 3(b), and 4) displays bedding rotated ~30° counterclockwise relative to the average orientation of Saddlebag Lake pendant units in the mapped area, with an average orientation of 133/83, also seen in foliation measurements (average orientation 323/81; Figure 2). Dana sequence units are eastward-dipping and gently to tightly folded, with an average orientation of 318/79 (Figure 2). Magmatic foliations (147/83) and well-defined flow bands (141/71) in the porphyry unit are oriented in a NW-SE direction (Figures 2 and 6). Along a SW-NE transect, layering of flow-bands, clastic layers, and mingling zones generally dip between 55–70°, shallower than the regional structural grain (Figures 2 and 6).

The Sierra Crest Shear Zone is a ~300 km long dextral-transpressive shear zone in the eastern Sierra Nevada, recording ductile deformation from at least 95 Ma to 84 Ma based on pluton ages and biotite cooling ages, switching to brittle behavior between 84-80 Ma [28, 34, 5153]. The porphyry intrusion locally has a solid-state foliation defined by deformed quartz grains and biotite phenocrysts that roughly parallels the earlier magmatic foliation. This solid-state foliation occurs at the western margin of the porphyry, where it intersects with the eastern margin of the Sierra Crest Shear Zone (“Cretaceous ductile shear zone” in Figure 2). The eastern margin of the Sierra Crest Shear Zone also deforms Triassic andesite flows and breccias in the Saddlebag Lake pendant. The western edge of the shear zone extends to the margin of the Kuna Crest granodiorite of the Tuolumne Intrusive Complex. Shear bands, boudinaged quartz veins, and σ-porphyroclasts are common structures within the sheared parts of the porphyry intrusion and generally indicate dextral shear.

Faults in a NW striking brittle fault system are locally quartz-filled or associated with the mineralization of tourmaline and epidote. Dextral strike-slip motion is the dominant fault style based on offset markers between stratigraphic units and measured slickenlines. The Gaylor Peak fault, a Jurassic thrust fault [19, 29], is one example of a reactivated through-going brittle fault in the area. The mapped fault trace cross-cuts the ca. 100 Ma porphyry unit (Figure 2), suggesting that the fault was reactivated in the Late Cretaceous. Subhorizontal slickenlines are found along the fault trace, indicating strike-slip movement on the fault. Several smaller, near-vertical strike-slip brittle faults cross-cut the porphyry unit with minor offsets (<100 m). These faults represent the transition from ductile to brittle behavior of the shear zone [51].

4.5. Whole-Rock Geochemistry of Tioga Pass Intrusive Rocks

Major and trace element whole-rock geochemistry from Tioga Pass intrusive rocks (Figure 9) illustrates wide inter- and intra-unit compositional variation. Fields of Cretaceous plutonic, hypabyssal, and volcanic rock samples, between 105 and 95 Ma, from the central Sierra Nevada (data sources listed in Figure 9 and Supplementary Materials) are shown for comparison with the Tioga Pass units. All four groups show similar trends with overlap in composition. Plutonic rocks span the widest range in composition, while hypabyssal samples show a narrow range in composition.

4.5.1. Major Elements

Quartz monzodiorite samples contain between 55 and 64 wt.% SiO2, whereas the silica content in the porphyry unit ranges between 63-73 wt.%. Al2O3, CaO, Fe2O3 (total Fe expressed as Fe2O3), TiO2, MnO, P2O5, and MgO decrease in abundance with increasing SiO2 contents in both units (Figure 9(a); Table 2). Samples of the porphyry unit are dominantly peraluminous with ASI between 1.02 and 1.15 (ASI = molar (Al2O3)/[(CaO-3.33P2O5) + Na2O + K2O)]). One metaluminous porphyry sample (sample KA20) with ASI 0.86 crops out at the contact between the pluton and the porphyry. The quartz monzodiorite has an ASI between 0.90 and 1.06. The single peraluminous sample in this unit (KA26) is also from the contact zone. These samples from the contact have similar major element compositions and also share a compositional affinity with porphyry and quartz monzodiorite samples collected from the southwest and northwest slopes of Mt. Dana, respectively.

One sample of the quartz monzodiorite (sample KA23; 56 wt.% SiO2) has high amounts of the major elements, e.g., 6.64 wt.% CaO, 9.42 wt.% Fe2O3, and 1.19 wt.% TiO2, compared to the other samples. A porphyritic mafic magmatic enclave (sample TP16-1a) with mm-size feldspar phenocrysts has a similar SiO2 content to KA23, but has higher MgO, Al2O3, and lower Fe2O3. Within the porphyry unit, there is a compositional trend in major elements linked to spatial location, as the porphyry becomes gradually more felsic towards the west (Figure 6; Figure 9 arrows); a pattern that is not observed in the quartz monzodiorite unit.

4.5.2. Alkali Enrichment

Samples in both units show a trend of increasing K2O and decreasing Na2O with increasing SiO2 (Figure 9(b)). Barium and Rb also increase with increasing K2O (Figure 9(c)). Sample KA13L contains anomalously high concentrations of K2O and Ba (7.41 wt.% and 3208 ppm, respectively). In addition, elevated K/Na, K/Al, and Ba/Ti in sample KA13L suggests that this sample experienced pervasive alkali metasomatism [2].

4.5.3. Minor and Trace Elements

Strontium decreases with increasing SiO2, and with decreasing CaO in the porphyry unit (Figure 9(d)), both decreasing towards the top of the mapped section (Figure 9 arrow). In the quartz monzodiorite, Sr varies from 500-600 ppm but does not vary with CaO. In both units V and Sc decrease with increasing SiO2 and increase with increasing TiO2 and CaO (Figures 9(e) and 9(f)). With increasing SiO2, Y, Zr, and Rb in the porphyry unit increase, but Y decreases slightly in the quartz monzodiorite unit; Zr and Rb in the quartz monzodiorite increase (Figure 9(g)). Compositional overlap between units is evident in all the above elements between 62-65 wt.% SiO2. Rare earth element (REE) patterns of the porphyry unit and the quartz monzodiorite display similar normalized La/Lu ratios of 7.3-11.1 and 8.5-15, respectively, although samples from the porphyry have a slight negative Eu anomaly that is absent in quartz monzodiorite samples (Figure 9(h)).

4.5.4. Isotopes

Initial strontium and neodymium isotope ratios (87Sr/86Sri and εNd) for the quartz monzodiorite sample (KA23) are 0.704997 and -0.29, respectively. In the porphyry sample (KA17), these ratios are 0.705494 and -1.33. εNd and 87Sr/86Sri values are within the mantle array [54] and are within the isotopic range of Cretaceous Sierra Nevada peridotite and pyroxenite xenoliths [55] (Figure 10(a)). Initial Pb isotopic ratios (207Pb/204Pb, 208Pb/204Pb, and 206Pb/204Pb) in the quartz monzodiorite are 15.67433, 38.97684 and 19.15353, respectively. The porphyry unit has initial Pb ratios (207Pb/204Pb, 208Pb/204Pb, and 206Pb/204Pb) of 15.65512, 38.67369, and 18.9539, respectively. Lead isotopes from the porphyry and quartz monzodiorite unit are also consistent with xenolith data [56] and are within the field defined by western US passive margin deep marine sediments [57] (Figure 10(b)).

The Tioga Pass area of the central Sierra Nevada contains a Cretaceous hypabyssal magmatic complex intruding into regionally extensive host-rock strata of the Saddlebag Lake and northern Ritter Range pendants, as well as the Tioga Lake and Dana sequence sections (Figure 11). Evidence for this new interpretation includes (1) ca. 100 Ma U-Pb zircon ages of the porphyry and quartz monzodiorite; (2) the gradational contact between these two units and overlap in compositions indicating comagmatic emplacement; and (3) field and petrographic features characteristic of shallow (subvolcanic) emplacement. Field evidence combined with element and isotope compositions are permissive of a “magma feeder system” model (Figure 11) where fractionation was significant in producing the compositional zoning and peraluminous rhyolite from mantle-derived magmas. Specifically, amphibole (ASI ~0.5; Zen [58]) and clinopyroxene (ASI ~0; Zen [58]) fractionation can drive the resulting melt to peraluminous compositions (e.g., Zen [58], Nandedkar et al. [59], Clemens et al. [60]).

5.1. Tioga Pass Cretaceous Magmatic Complex

Weighted mean ages from the Tioga Lake quartz monzodiorite, granite, and porphyry unit are indistinguishable within uncertainty. As the contact between the granite and the other Cretaceous intrusions is concealed, interpretation of the relationship to the granite unit is hindered. However, field evidence from the quartz monzodiorite and the porphyry unit demonstrate that they are comagmatic across a contact that varies from gradational to sharp along strike. In addition, inclusions of the quartz monzodiorite are found within the porphyry and vice versa. The overlap in major and trace elements between units (at 62-65 wt.% SiO2) near the contact zone, and the overlap in REE patterns, is consistent with a comagmatic relationship.

Observations consistent with a hypabyssal interpretation include the discordant intrusive contact of the porphyry unit and quartz monzodiorite with the surrounding host rocks, the porphyritic texture, the glassy inclusions and clasts, as well as the glomerocrysts and fractured crystals. The glass inclusions in particular are indicative of rapidly cooled magma. Devitrification is interpreted where the glassy inclusions are found together with the recrystallized matrix [61, 62]. Layers including clastic material suggest interaction with the surface, or near-surface, environment, either during magma ascent or possible eruption. Plagioclase and alkali feldspar phenocrysts are clustered, forming the distinctive glomeroporphyritic structure of the porphyry unit. In this case, the preferential contacts of like minerals are observed (e.g., plagioclase-plagioclase clusters are abundant), which is a common feature of hypabyssal rocks (e.g., [63]). Melt extraction from a crystal mush and/or the tectonic disruption of the mush could produce feldspar glomerocrysts in porphyries and volcanic rocks [6, 6466]. Many of the phenocrysts are fractured, which could occur by a volcanic eruption, a decrease in pressure due to magma ascent, or the rapid “ungluing” of crystal groups [65, 6769].

The quartz monzodiorite transitions from an equigranular to porphyritic texture close to the contact with the porphyry unit, but the largest textural and structural variations are within the porphyry unit. The quartz monzodiorite and porphyry share similar mineralogy but differ in modal proportions of quartz, alkali feldspar, and biotite (more abundant in the western extent of the porphyry), and there is a decrease in hornblende and biotite phenocryst abundances westwards in the porphyry unit. The porphyry unit is stratified by rock type and structure, with layers uniformly dipping (Figure 6). Topographically, the porphyry grades from a crystal-rich dacite to a flow banded rhyolite at higher levels, and many of the major and trace elements are zoned from the base of the section to the top (Figures 9 and 11). This suggests that layers in the porphyry unit preserve their original orientation and have not been steeply tilted. Our tests of sample alteration (e.g., [2], alkali enrichment) suggest that this is a primary chemical gradation. Further, the compositional zoning implies that either the porphyry is made up of multiple pulses or gradually fractionated from a dacite to a high-silica rhyolite (Figure 12).

Plagioclase appears to be an important mineral defining the intermediate-felsic composition of the intrusive units. It is likely controlling the decrease in Sr, and the high concentration of Sr and CaO in some quartz monzodiorite and porphyry samples suggests that these samples accumulated plagioclase. The quartz monzodiorite does not show the negative Eu anomaly of the porphyry unit nor a positive anomaly, but retains a subtle signature of feldspar accumulation at the upper-crustal levels. This interpretation is consistent with observations of the elevated plagioclase content, plagioclase-plagioclase contacts in the quartz monzodiorite, and plagioclase glomerocrysts of the porphyry unit. Zoning of K2O and Rb across the mapped section suggests an important role for K-feldspar fractionation, except in sample KA13L.

The porphyry is distinct from the quartz monzodiorite in that it has a peraluminous composition (with the exception of one sample at the contact zone reported above). In addition to the role of fluids and volatiles (e.g., [70]), one explanation is that the porphyry may contain a higher proportion of highly fractionated melt (preserved in the rock as glass or recrystallized glass matrix). The metaluminous quartz monzodiorite on the other hand is likely a cumulate that lost melt and accumulated feldspar. Another possibility is that the porphyry unit assimilated some of the surrounding host rocks (Triassic volcanic rocks or Paleozoic chert/argillite unit), supported by field observations of stoped blocks (Figures 6 and 11). However, these blocks are found in both units, so assimilation of the blocks does not entirely explain the peraluminosity of the porphyry unit. Generating peraluminous rhyolite compositions in the porphyry unit from a mantle-derived magma (e.g., [71]; discussion below) suggest that fractionation played a significant role in generating the compositional variety observed in the porphyry unit.

A depleted mantle or oceanic crust (MORB) magma source (εNd>+7.5; [54]) does not fully explain the magmatic source of the Tioga Pass intrusions, which have negative εNd values, unless the depleted mantle source was progressively enriched and modified during earlier arc activity. In addition, the samples do not share isotopic affinity with enriched mantle reservoirs (EMI and EMII; Figure 10(a)). However, initial Sr, Nd, and Pb isotope ratios of the porphyry unit and quartz monzodiorite match Cretaceous lower crust and mantle lithosphere xenolith compositions [55], suggesting that the porphyry and quartz monzodiorite magmas are largely mantle-derived. Further, Lackey et al. [71] measured δ18O in zircon from the Tioga Lake quartz monzodiorite (5.98‰) and found values slightly elevated above the mantle range (5.3±0.3; [72]). Initial Pb isotopes indicate that the magmas incorporated deep-marine sediments of the western Cordilleran passive margin ([57]; Figure 10(b)), which is compatible with field observations of the Paleozoic chert-argillite unit in contact with the Tioga Pass intrusions and detrital zircon studies that show that the age range from this type of unit is regionally extensive in the eastern-central Sierra Nevada (e.g., [30] and references therein). Although crustal contamination was limited, based on Rb-Sr, Sm-Nd, and O isotope systems, the marine sediments (and the underlying basement/lithospheric mantle) represent one likely assimilant.

Isotopes from Tioga Pass samples are relatively primitive compared to the isotopic compositions of intrusive suites from the central Sierra Nevada, which are on average 87Sr/86Sri>0.706 and εNd<1.7 (e.g., [7375]). This signal is not unique to Tioga Pass, and other Late Cretaceous intrusions with a similar isotopic composition, emplaced south of Tioga Pass, include the ca. 97 Ma quartz monzodiorite of Rush Creek, which is comparable to the Tioga Lake quartz monzodiorite in composition and size ([4, 76]; Figure 10(a)), and the hypabyssal ca. 100 Ma Shellenbarger Lake granite porphyry ([18]; Figures 10(a) and 10(b)). To the north of Tioga Pass, the 97 Ma Solider Lake granodiorite [34] and a Cretaceous volcanic sample are isotopically similar to the Tioga dacite porphyry unit. The eastern margin of the Kuna Crest granodiorite and Kuna Crest lobe, emplaced at 95-94 Ma ([32], Memeti et al., in review), contain the most primitive compositions of the Tuolumne Intrusive Complex ([77]; Figures 412 in [78]; Memeti et al., in review) and represent much larger volumes of magma than the other examples. These intrusions and volcanic rocks are all located within an arc-parallel belt in the eastern-central Sierra Nevada that is also defined by δ18Ozrc approximating, or slightly higher than, mantle values between 5.5-6.5‰ [71]. Lackey et al. [71] interpreted this as a zone of minimal crustal contamination and extensive recycling of mantle lithosphere and upper-mantle (see also [79]). The similarity in isotope compositions suggests each of these intrusions had a similar mantle-derived magma source to the Tioga Pass system, including a minor component of deep marine sediments (Figure 10). It further indicates that the whole-rock Sr, Nd, and Pb isotopes across this eastern belt capture the basement composition of the arc framework rocks (e.g., [80]).

5.2. Alternate Interpretation of the Tioga Pass Magmatic Complex

An alternate interpretation of the field and geochronologic data is that the porphyry represents the subvolcanic (intrusive) roots of a lava dome (e.g., [81]). The lack of vesicles within this unit suggests that the magma was already degassed or the microstructure was overprinted by Cretaceous regional tectonic shortening [34, 8284]. Flow-banding structures are characteristic features of lava flows and domes, both in the upper parts of magma feeder conduits and at the surface [85, 86], a feature which is ubiquitous in the porphyry unit. Glomerocrysts have been recorded in extrusive lava flows and domes (e.g., [66]). The sharp, discordant intrusive contact between the intrusive units and older host rocks requires this to be a subvolcanic system but could have fed an extrusive structure such as a lava dome or stratovolcano (e.g., [4]).

5.3. Tectonic Implications

New geochronology from volcanic and sedimentary units at Tioga Pass expands our understanding of the regional early Mesozoic stratigraphy of the Saddlebag Lake and northern Ritter Range pendants. In combination with field relationships, this resolves differences between Triassic and Cretaceous tectonic histories and aids in refining the local structural history for this area. Triassic and Jurassic deposition is recorded in the Tioga Lake and Dana sequence sections, preserving the transition from Triassic emergent volcanism to shallow marine sedimentation. We suggest that correlative in situ units from the Saddlebag Lake and northern Ritter Range pendants are not observed because they were removed prior to, or during, emplacement of the 95-85 Ma Tuolumne Intrusive Complex, which truncates the ca. 220 Ma volcaniclastic unit at Tioga Pass. Some displacement of the Tioga Lake and Dana sequence sections may have been accommodated by the Gaylor Peak fault, as proposed by Greene [36] and Schweickert and Lahren [19]. If so, this must postdate Dana sequence formation (195 Ma) and precede the emplacement of the magmatic complex at 100 Ma, which intrudes these units. Subsequent strike-slip reactivation of this fault offsets the porphyry unit (Figure 11). Rotation of bedding and foliation in the Tioga Lake block could alternatively have occurred during the emplacement of the porphyry unit; rotated structures are also observed in <50 m stoped host-rock blocks within the porphyry unit.

Repeated episodes of deformation and tilting of strata occurred in the Triassic and Jurassic (e.g., [28, 34]); however, Cretaceous units of the Tioga Pass intrusive complex are less tilted and less strained than the older rocks in the same area. This may indicate that by 100 Ma (at the peak of the Cretaceous flare-up), tectonically driven steepening and shortening were waning. In intrusive rocks, the porphyry preserves evidence for ductile shear during Sierra Crest Shear Zone activity that occurred after porphyry emplacement at ca. 100 Ma, but prior to brittle faulting between 84-80 Ma (e.g., [51]).

The Kuna Crest granodiorite, Tioga Lake quartz monzodiorite, hypabyssal porphyry, and Cretaceous volcanic rocks are structurally juxtaposed at the present-day surface (Figures 1, 3(a)–3(c)). This tectonic association between plutons and volcanic rocks of similar age has been identified at other localities in the central Sierra Nevada (e.g., [14, 22]) and is attributed to downward transfer of host rock (e.g., [27, 47, 87, 88]). Using the minimum and maximum U-Pb zircon ages for each unit and Al-in hornblende barometry of the Kuna Crest granodiorite ([21]; Memeti et al., in review), the rate of downward transfer is estimated between 1.6–4 km/m.y., consistent with pendant-wide estimates from Cao et al. [47]. Whether downward transfer was episodic (over geologic timescales) or continuous remains to be tested.

5.4. Implications for Hypabyssal Intrusions in the Central Sierra Nevada and the Volcanic-Plutonic Connection

The Tioga Pass magmatic complex is one example of a shallowly emplaced, porphyritic magma feeder system in the central Sierra Nevada (Figure 11). Other Cretaceous hypabyssal intrusions scattered along a ~50-kilometer-wide belt of the central Sierra Nevada are recognized from a synthesis of published literature and our ongoing field studies using the field and textural criteria outlined in this study. Known hypabyssal intrusions are labeled in Figure 1 and are reported as far south as Kings Canyon National Park in the southern Sierra Nevada (not shown; [16, 89, 90]). Below we outline some of the general features of these intrusions in the context of our findings at Tioga Pass.

Hypabyssal intrusions across the central Sierra Nevada are typically small in areal extent (10 m–10 km in diameter) and have a range of 3D shapes, including punched laccoliths (flat-topped, steep-sided intrusions: [20], Figure 11), pipes, and porphyritic dikes [13, 17, 88]. They sometimes intrude through slightly older metavolcanic rocks; the ca. 100 Ma Shellenbarger Lake granite porphyry cuts across the Minarets caldera fill deposit and is interpreted as a hypabyssal resurgent dome to the Minarets caldera (Figure 1; [17, 86]).

Field observations of hypabyssal intrusions highlight the diversity of structures found at shallow levels and suggests that a range of magmatic conditions may be captured in these systems. In the central Sierra Nevada, hypabyssal intrusions are dominantly intermediate to felsic in composition, such as andesite, dacite, monzodiorite, granodiorite, and granite. Basaltic andesite compositions are locally found ([17]; S. Attia, pers. comm). Textures are largely porphyritic, with phenocrysts of alkali feldspar, plagioclase, and biotite common [18, 91, 92]. Phenocrysts are set in a groundmass that ranges from medium-grained (akin to typical plutonic equigranular textures) to fine-grained or microcrystalline (volcanic aphanitic textures). In the latter case, the only differentiating feature of a hypabyssal stock from a volcanic deposit is a discordant, intrusive contact with surrounding strata or older plutonic bodies [13, 90]. As exemplified at Tioga Pass, hypabyssal bodies often have complex macro- to microscale internal structures. In many instances they contain miarolitic cavities and granophyric microstructures, which are compatible with volatile saturation and shallow emplacement at approximately 1 kbar [93, 94], although miarolitic cavities are also found in the Tuolumne Intrusive Complex (6–10 km; [78]).

Compositions of the hypabyssal intrusions can be used to reconstruct source characteristics, test crystal accumulation processes, and explore the possible genetic associations between plutons, hypabyssal intrusions, and volcanic rocks in the upper crust. Cretaceous plutons in the central Sierra Nevada are magnesian, calc-alkaline, metaluminous to peraluminous granodiorite, and granite in composition [13]. Porphyry and volcanic fields overlap considerably with plutonic compositions (Figure 9). However, the hypabyssal samples trend towards peraluminous and alkali-calcic compositions and are generally more compositionally restricted than plutons, between 62–76 wt.% SiO2 (e.g., [18, 76, 95, 96]). Our findings at Tioga Pass extend the compositional range of known hypabyssal intrusions in the central Sierra Nevada (Figure 9), and in some cases, results are quite distinct from regional patterns (e.g., Yttrium; Figure 9). In addition to the Tioga Pass system, other Late Cretaceous plutons and hypabyssal intrusions emplaced along the eastern edge of the central Sierra Nevada have relatively primitive bulk isotopic compositions (Figure 10(a)) compared to average values for the arc section, illustrating a regional control on mantle and crustal components in the magma source that extends into the upper crust [71, 80].

Field relationships combined with chemical compositions of plutonic, hypabyssal, and volcanic rocks have been studied in other settings to characterize volcanic-plutonic systems and refine existing volcanic-plutonic models (Figure 12). For example, caldera structures that formed during the Cenozoic ignimbrite flare-up in western North America document a wide variety of structural and chemical volcanic-plutonic relationships, made possible due to the exposures of the magma plumbing from the surface to 9 km depths (e.g., [3, 6, 8, 97, 98]). Volcanic and shallow plutonic rocks within each system typically share a “compositional affinity” (e.g., [98]) and, in some cases, shallow plutons are interpreted as the residual magmas to erupted volcanic deposits (e.g., [9, 99]; Figure 12-magma feeder system model). Often, the exposed shallow plutons are resurgent into the erupted volcanic rocks, and thus did not directly feed volcanism (e.g., [3, 100, 101]; Figure 12-failed eruption model). In each example, shallow intrusions (0-6 km depth) provide windows into different levels of the magmatic system through time. The study of these intrusions is particularly important in areas where surface deposits and the uppermost crust are largely eroded, such as the Sierra Nevada.

(1) The Tioga Lake quartz monzodiorite and dacite-rhyolite porphyry are comagmatic intrusive units overlapping in age at ca. 100 Ma and are not related to a Triassic caldera system. Our favored interpretation is that the intrusions were emplaced at shallow crustal levels and represent part of the subvolcanic magmatic roots of once-extensive Cretaceous volcanic deposits, now partially preserved in host-rock pendants across the central Sierra Nevada. Thus, the only known Sierran caldera structures are the Cretaceous Minarets and Merced peak calderas.

(2) Whole-rock isotopes from Tioga Pass, and nearby intrusions and volcanic rocks along the eastern belt of the central Sierra Nevada, are sensitive to the composition of regional arc framework rocks, including the underlying basement and mantle sources. Mantle-derived magmas along this eastern belt contain a minor crustal component consistent with the composition of the host Paleozoic deep marine sediments. The Tioga Lake and Dana sequence sections represent ~2 km of Triassic and Jurassic strata that is otherwise not preserved at Tioga Pass, due to surface or magmatic erosion. Cretaceous ductile and brittle activity reactivated earlier structures, and rapid downward flow of rocks juxtaposed volcanic, hypabyssal, and plutonic rocks.

(3) Subvolcanic porphyry intrusions in the Cretaceous central Sierra Nevada are widespread and are also documented in other parts of the arc section (e.g., [16, 90]). Combining this dataset with temporally and spatially associated volcanic and plutonic rocks provides a multilevel view of the uppermost parts of an arc crustal column, owing to either differential erosion and/or downward transfer processes that transported volcanic rocks erupted at the surface to 6–10 km depth before intrusion of broadly coeval plutons. In this regard, porphyry intrusions represent key structural markers in tectonic studies, where they are constrained by a crystallization age and emplacement depth (between approximately 0–6 km). We posit that these hypabyssal intrusions are also significant in that they may physically and chemically relate to volcanic deposits and/or the late melts drained from plutons

Geochronology and geochemical datasets collected in this study are provided in Tables 1 and 2 in the manuscript, and in spreadsheets in the Supplementary File.

The authors declare that they have no conflicts of interest.

We acknowledge support from the National Science Foundation grants EAR-1624847 awarded to Scott Paterson, and EAR-1624854 awarded to Vali Memeti. We thank Mark Pecha and research scientists from the Arizona LaserChron Center (NSF EAR-1649254) for help with obtaining U-Pb zircon ages. We thank James McColl, Luke Ardill, and Maria Cuevas for field assistance, Alejandra Angulo for help preparing XRF beads, and Melanie Barnes and Kevin Werts for helping us obtain trace element data. Cal Barnes, Wenrong Cao, and Snir Attia are thanked for useful discussions and comments. Snir Attia provided samples of mafic hypabyssal rocks to include in our compilation. Helpful reviews by Jade Star Lackey, Kathryn Watts, Andrew Barth, Lily Claiborne, Joseph Colgan, Nancy Riggs, and Gerardo J Aguirre-Diaz improved the earlier versions of this manuscript. We thank Sarah Roeske for editorial handling.

Exclusive Licensee GeoScienceWorld. Distributed under a Creative Commons Attribution License (CC BY 4.0).