skip to main content


Title: Open-system Evolution of a Crustal-scale Magma Column, Klamath Mountains, California
Abstract This study addresses the question of how and where arc magmas obtain their chemical and isotopic characteristics. The Wooley Creek batholith and Slinkard pluton are a tilted, mid- to upper-crustal part of a vertically extensive, late-Jurassic, arc-related magmatic system in the Klamath Mountains, northern California. The main stage of the system is divided into an older lower zone (c. 159 Ma) emplaced as multiple sheet-like bodies, a younger upper zone (c. 158–156 Ma), which is gradationally zoned upward from mafic tonalite to granite, and a complex central zone, which represents the transition between the lower and upper zones. Xenoliths are common and locally abundant in the lower and central zones and preserve a ghost stratigraphy of the three host terranes. Bulk-rock Nd isotope data along with ages and Hf and oxygen isotope data on zircons were used to assess the location and timing of differentiation and assimilation. Xenoliths display a wide range of εNd (whole-rock) and εHf (zircon), ranges that correlate with rocks in the host terranes. Among individual pluton samples, zircon Hf and oxygen isotope data display ranges too large to represent uniform magma compositions, and very few data are consistent with uncontaminated mantle-derived magma. In addition, zoning of Zr and Hf in augite and hornblende indicates that zircon crystallized at temperatures near or below 800 °C; these temperatures are lower than emplacement temperatures. Therefore, the diversity of zircon isotope compositions reflects in situ crystallization from heterogeneous magmas. On the basis of these and published data, the system is interpreted to reflect initial MASH-zone differentiation, which resulted in elevated δ18O and lowered εHf in the magmas prior to zircon crystallization. Further differentiation, and particularly assimilation–fractional crystallization, occurred at the level of emplacement on a piecemeal (local) basis as individual magma batches interacted with partial melts from host-rock xenoliths. This piecemeal assimilation was accompanied by zircon crystallization, resulting in the heterogeneous isotopic signatures. Magmatism ended with late-stage emplacement of isotopically evolved granitic magmas (c. 156 Ma) whose compositions primarily reflect reworking of the deep-crustal MASH environment.  more » « less
Award ID(s):
2004618
NSF-PAR ID:
10334916
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Journal of Petrology
Volume:
62
Issue:
12
ISSN:
0022-3530
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract The 119 Ma Dinkey Dome pluton in the central Sierra Nevada Batholith is a peraluminous granite and contains magmatic garnet and zircon that are complexly zoned with respect to oxygen isotope ratios. Intracrystalline SIMS analysis tests the relative importance of magmatic differentiation processes vs. partial melting of metasedimentary rocks. Whereas δ18O values of bulk zircon concentrates are uniform across the entire pluton (7.7‰ VSMOW), zircon crystals are zoned in δ18O by up to 1.8‰, and when compared to late garnet, show evidence of changing magma chemistry during multiple interactions of the magma with wall rock during crustal transit. The evolution from an early high-δ18O magma [δ18O(WR) = 9.8‰] toward lower values is shown by high-δ18O zircon cores (7.8‰) and lower δ 18O rims (6.8‰). Garnets from the northwest side of the pluton show a final increase in δ18O with rims reaching 8.1‰. In situ REE measurements show zircon is magmatic and grew before garnets. Additionally, δ18O in garnets from the western side of the pluton are consistently higher (avg = 7.3‰) relative to the west (avg = 5.9‰). These δ18O variations in zircon and garnet record different stages of assimilation and fractional crystallization whereby an initially high-δ18O magma partially melted low-δ18O wallrock and was subsequently contaminated near the current level of emplacement by higher δ18O melts. Collectively, the comparison of δ18O zoning in garnet and zircon shows how a peraluminous pluton can be constructed from multiple batches of variably contaminated melts, especially in early stages of arc magmatism where magmas encounter significant heterogeneity of wall-rock assemblages. Collectively, peraluminous magmas in the Sierran arc are limited to small <100 km2 plutons that are intimately associated with metasedimentary wall rocks and often surrounded by later and larger metaluminous tonalite and granodiorite plutons. The general associations suggest that early-stage arc magmas sample crustal heterogeneities in small melt batches, but that with progressive invigoration of the arc, such compositions are more effectively blended with mantle melts in source regions. Thus, peraluminous magmas provide important details of the nascent Sierran arc and pre-batholithic crustal structure. 
    more » « less
  2. PhD Dissertation Abstract: The imposing andesite stratovolcano is the characteristic expression of subduction zone magmatism, posing hazards to coastal populations and bearing insight into deep Earth processes. On a map of a typical volcanic arc, one can easily distinguish the approximately linear alignment and regular spacing of these major edifices that stand out from a diffuse distribution of mafic volcanoes (e.g. the Quaternary Cascades; Hildreth, 2007). The andesitic composite volcanoes have a reputation for being complex, open systems: crystal zoning “stratigraphies,” diverse crystal cargoes including antecrysts or xenocrysts, quenched magmatic inclusions, and variations in isotopic signatures are among the many lines of evidence that these systems involve a variety of igneous processes and melt sources. To investigate the development and evolution of such transcrustal magma factories, I have conducted a detailed temporal, spatial, and geochemical characterization of a long-lived arc volcanic center in the southern Washington Cascades, the Goat Rocks volcanic complex. Results from ⁴⁰Ar/³⁹Ar and U/Pb geochronology constrain the lifespan of the Goat Rocks volcanic complex from ~3.1 Ma to ~100 ka. During this time, four major composite volcanoes were built (as well as several smaller volcanoes). From oldest to youngest, these are Tieton Peak, Bear Creek Mountain, Lake Creek volcano, and Old Snowy Mountain. Four volcanic stages are defined based on the lifespans of these centers and distinct compositional changes that occur from one to the next: Tieton Peak stage (3.1-2.6 Ma), Bear Creek Mountain stage (1.6-1.1 Ma), Lake Creek stage (1.1 Ma to 456 ka), and Old Snowy Mountain stage (440 ka to 115 ka). Two lava flow remnants also have ages in the interim between Tieton Peak stage and Bear Creek Mountain stage (2.3 Ma and 2.1 Ma), and their sources are not yet identified. The ages of the Bear Creek Mountain and Lake Creek stages in fact overlap, and the gap between Lake Creek stage and Old Snowy Mountain stage is only on the order of 10⁴ years. Based on supporting compositional evidence, the Bear Creek Mountain, Lake Creek, and Old Snowy Mountain stage volcanoes are considered to be the migrating surface expressions of a continuous magmatic system that was active over at least ~1.5 million years. It remains uncertain whether the gaps between the Tieton Peak stage, scattered early Pleistocene andesites, and Bear Creek Mountain stage are due to incomplete exposure/sampling or real quiescent periods earlier in the development of the Goat Rocks volcanic complex. Throughout the construction of the andesitic complex, mafic volcanoes were active on its periphery. These include the Miriam Creek volcano (~3.6-3.1 Ma), Devils Washbasin volcano (3.0-2.7 Ma), Hogback Mountain (1.1 Ma – 891 ka), Lakeview Mountain (194 ka), and Walupt Lake volcano (65 ka). Two basalt and basaltic andesite units (Qob₁ and Qob₂, 1.4 and 1.3 Ma; Hammond, 2017) also erupted from the Goat Rocks area, likely an older incarnation of Hogback Mountain. The suite of mafic magmas erupted in this region are all calcalkaline basalt (or basaltic andesite; CAB), but two compositional groups emerge from the trace element and isotopic data. Group 1 is LILE and LREE-enriched, with higher ⁸⁷Sr/⁸⁶Sr isotopes, and includes compositions from Devils Washbasin, Lower Hogback Mountain, and Lakeview Mountain. Group 2 is less enriched in LILE and LREE and lower in ⁸⁷Sr/⁸⁶Sr, and includes the compositions of Miriam Creek, Qob1, Upper Hogback Mountain, Walupt Lake, and Coleman Weedpatch. The two groups are recurrent through time and with no geographic distinction; in fact, both types were tapped by the Hogback Mountain volcano. Together both of these groups, alongside CABs from Mount Adams and various basalts from Mount St. Helens, form a compositional array between the basalts of the High Cascades and the intraplate-type basalts (IPB) of Mount Adams and Simcoe volcanic field. These results lead to three conclusions. 1) Variably subduction-modified mantle is distributed across the region, perhaps either as stratified layers or a web-like network of fluid pathways amongst less metasomatized mantle. 2) Transitional compositions between the IPBs and typical “High Cascades” CAB/HAOT signature suggest a broader influence of the mantle domain that feeds IPBs—if asthenospheric mantle through a slab window, as suggested by Mullen et al. (2017), then perhaps it bleeds in smaller quantities over a broader area. This compositional trend solidifies the interpretation of the southern Washington Cascades as a unique and coherent “segment” of the arc (the Washington segment of Pitcher and Kent, 2019). 3) The recurrence of variable mafic magma types through time, and with no geographic boundaries, indicates that the compositional evolution of the Goat Rocks volcanic complex was not likely driven by a change in mafic input. Indeed, the Sr, Nd, Hf, and Pb isotope ratios of the intermediate to felsic suite are closely aligned with the local basalts and suggest a limited role of crustal assimilation. Importantly, several mineral thermometers (zircon, ilmenite-magnetite pairs, and amphibole) align in recording higher crystallization temperatures in Bear Creek Mountain to early Lake Creek time, a cooling trend through the Lake Creek stage, and a more diverse range of temperatures in the transition to Old Snowy Mountain stage. Thus, it is suggested that the compositional evolution at Goat Rocks represents a thermal cycle of waxing and waning magmatic flux: where the period of Bear Creek Mountain to early Lake Creek volcanism was the climactic phase of a vertically extensive magma homogenization factory, then the system waned and cooled, ultimately losing its ability to filter, homogenize, and enrich magmas. 
    more » « less
  3. Zircon U-Pb, and garnet Sm-Nd and Lu-Hf dates provide important constraints on local and orogenic scale processes in lower-crustal rocks. However, in high-temperature metamorphic rocks these isotopic systems typically yield significant ranges reflecting both igneous and metamorphic processes. Therefore, linking dates to specific aspects of rock history can be problematic. In Fiordland, New Zealand, granulite-facies orthogneiss is cut by leucosomes that are bordered by garnet clinopyroxene reaction zones (garnet reaction zones). In both host orthogneiss and garnet reaction zones, zircon are typically anhedral with U-Pb dates ranging from 118.30 ± 0.13 to 115.70 ± 0.18 Ma (CA-ID-TIMS) and 121.4 ± 2.0 to 109.8 ± 1.8 Ma (SHRIMP-RG). Zircon dates in host and garnet reaction zone do not define distinct populations. In addition, the dates cannot be readily grouped based on external morphology or internal CL zoning. Zircon trace-element concentrations indicate two distinct crystallization trends, clearly seen in Th and U. Garnet occurs in selvages to the leucosome veins and in the adjacent garnet reaction zones. In selvages and host orthogneiss, garnet is generally 0.5 to 1 cm diameter and euhedral and is 0.1 to 0.5 cm diameter and subhedral in garnet reaction zones. Garnet Sm-Nd and Lu-Hf dates range from ca. 115 to 101 Ma (including uncertainties) and correlate with grain size. We interpret the CA-ID-TIMS zircon dates to record the age of magma emplacement and the SHRIMP-RG dates to record a range from igneous crystallization to metamorphic dissolution and reprecipitation and/or local Pb loss. Zircon compositional trends within the garnet reaction zone and host are compatible with locally isolated melt and/or separate intrusive magma batches for the two samples described here. Dates for the largest, ~1 cm, garnet of ~113 Ma record growth during metamorphism, while the smaller grains with younger dates reflect high-temperature intracrystalline diffusion and isotopic closure during cooling. The comprehensive geochronological data set for a single location in the Malaspina Pluton illustrates a complex and protracted geologic history common in granulite facies rocks, estimates lower crustal cooling rates of ~20 °C/m.y., and underlines the importance of multiple chronometers and careful textural characterization for assigning meaningful ages to lower-crustal rocks. Numerous data sets from single locations, like the one described here, are needed to evaluate the spatial extent and variation of cooling rates for Fiordland and other lower crustal exposures. 
    more » « less
  4. Interpretation of geochronological and petrological data from partially-melted granulite is challenging. However, integration of multiple chronometers and mineral assemblage diagrams (MAD) can be used to estimate the nature and duration of processes. Excellent lower-crustal exposures of garnet granulite from the Malaspina Pluton, Fiordland New Zealand provide an ideal place to employ this kitchen sink approach. We use zircon U-Pb ages from LA-ICPMS, SHRIMP-RG, and CA-TIMS, garnet Lu-Hf and Sm-Nd ages, and MAD in order to evaluate local partial melting vs. melt injection, equilibrium volumes, P-T conditions, and the duration of lower crustal thermal events. Host diorite (H), garnet-clinopyroxene reaction zones (GRZ), coarse garnet selvages, and tonalite veins provide a record of intrusion and granulite facies partial melting. Zircon U-Pb ages range from 123 to 107 Ma (all); LA-ICPMS ages contain the entire range; CA-TIMS ages range from 118.30±0.13 to 115.7±0.18 Ma; and SHRIMP-RG ages range from 121.4±2 to 109.8±1.8 Ma. The latter two techniques are interpreted to indicate primary igneous crystallization from ~119 to ~116 Ma and the youngest ~110 Ma ages are interpreted as metamorphic zircon growth. Garnet ages for ~1 cm grains are ~113 Ma (Lu-Hf & Sm-Nd) and record metamorphic growth, and <0.3 mm grains with Sm-Nd ages from 113 to 104 Ma reflect high temperature intracrystalline diffusion and isotopic closure during cooling to amphibolite facies. Zircon trace-element compositions indicate 2 distinct crystallization trends reflecting evolution of primary magma batches. MAD indicate that garnet was not in equilibrium with sampled rock compositions. Instead, garnet shows apparent equilibrium with a modeled mixture of the GRZ and the H and grew in equilibrium with an effective bulk composition that shifted toward the leucosome. This would produce the observed increase in garnet grossular content. We conclude that: Malaspina rocks from Crooked Arm preserve evidence for 2 igneous layers which evolved as discrete magmas, igneous crystallization lasted 2 to 3 m.y., granulite metamorphism peaked ~ 3 m.y. after intrusion, metamorphism lasted ≥3 m.y., cooling occurred at ~20°C/m.y., and granulite minerals equilibrated with a mixture of solid phases and melt at ~14 kbar and 920°C (based on garnet compositions and MAD). 
    more » « less
  5. null (Ed.)
    Abstract Recovering the time-evolving relationship between arc magmatism and deformation, and the influence of anisotropies (inherited foliations, crustal-scale features, and thermal gradients), is critical for interpreting the location, timing, and geometry of transpressional structures in continental arcs. We investigated these themes of magma-deformation interactions and preexisting anisotropies within a middle- and lower-crustal section of Cretaceous arc crust coinciding with a Paleozoic boundary in central Fiordland, New Zealand. We present new structural mapping and results of Zr-in-titanite thermometry and U-Pb zircon and titanite geochronology from an Early Cretaceous batholith and its host rock. The data reveal how the expression of transpression in the middle and lower crust of a continental magmatic arc evolved during emplacement and crystallization of the ∼2300 km2 lower-crustal Western Fiordland Orthogneiss (WFO) batholith. Two structures within Fiordland’s architecture of transpressional shear zones are identified. The gently dipping Misty shear zone records syn-magmatic oblique-sinistral thrust motion between ca. 123 and ca. 118 Ma, along the lower-crustal WFO Misty Pluton margin. The subhorizontal South Adams Burn thrust records mid-crustal arc-normal shortening between ca. 114 and ca. 111 Ma. Both structures are localized within and reactivate a recently described >10 km-wide Paleozoic crustal boundary, and show that deformation migrated upwards between ca. 118 and ca. 114 Ma. WFO emplacement and crystallization (mainly 118–115 Ma) coincided with elevated (>750 °C) middle- and lower-crustal Zr-in-titanite temperatures and the onset of mid-crustal cooling at 5.9 ± 2.0 °C Ma−1 between ca. 118 and ca. 95 Ma. We suggest that reduced strength contrasts across lower-crustal pluton margins during crystallization caused deformation to migrate upwards into thermally weakened rocks of the mid-crust. The migration was accompanied by partitioning of deformation into domains of arc-normal shortening in Paleozoic metasedimentary rocks and domains that combined shortening and strike-slip deformation in crustal-scale subvertical, transpressional shear zones previously documented in Fiordland. U-Pb titanite dates indicate Carboniferous–Cretaceous (re)crystallization, consistent with reactivation of the inherited boundary. Our results show that spatio-temporal patterns of transpression are influenced by magma emplacement and crystallization and by the thermal structure of a reactivated boundary. 
    more » « less