Skip Navigation



Journal of Petrology Advance Access published online on June 19, 2008

Journal of Petrology, doi:10.1093/petrology/egn030
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Data
Right arrow All Versions of this Article:
49/7/1397    most recent
egn030v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Lackey, J. S.
Right arrow Articles by Stockli, D. F.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2008. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oxfordjournals.org

Dynamic Magma Systems, Crustal Recycling, and Alteration in the Central Sierra Nevada Batholith: the Oxygen Isotope Record

Jade Star Lackey1,*, John W. Valley2, James H. Chen3 and Daniel F. Stockli4

1Geology Department, Pomona College, Claremont, CA 91711, USA
2Department of Geology and Geophysics, University of Wisconsin, Madison, WI 53706, USA
3Science Division, Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA
4Department of Geology, University of Kansas, Lawrence, KS 66045, USA

Received January 6, 2007; Revised typescript accepted May 27, 2008


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL BACKGROUND
 SAMPLE SELECTION, PREPARATION,...
 RESULTS
 DISCUSSION
 CONCLUDING REMARKS
 SUPPLEMENTARY DATA
 REFERENCES
 
Values of {delta}18O of zircon from the central Sierra Nevada batholith (SNB), California, yield fresh insight into the magmatic evolution and alteration history of this classic convergent margin batholith. Direct comparison of whole-rock and zircon (Zrc) {delta}18O provides evidence for modest (0.5{per thousand}), but widespread, alteration, which has complicated interpretation in previous whole-rock {delta}18O studies. Four discrete belts of {delta}18O values are recognized in the central Sierra. A small belt of plutons with relatively low {delta}18O(Zrc) values (5·2–6·0{per thousand}) intrudes the foothills, with a sharp increase of {delta}18O revealing the concealed Foothills Suture; high {delta}18O(Zrc) values (7·0–8·5{per thousand}) dominate the rest of the western SNB. East of the axis of the Sierra, {delta}18O is distinctly lower (6·75–5·75{per thousand}), and decreases monotonically to the Sierra Crest. A sharp 1{per thousand} increase of {delta}18O in the eastern Sierra reveals a second crustal boundary, with the fourth belt hosted in high-{delta}18O North American crust in the White Mountains and Owens and Long Valleys. Correlated O, Sr, and Pb isotope ratios reveal differences in magma generation between the western and eastern Sierra. The western Sierra experienced massive crustal recycling, with substantial melting and mobilization of accreted oceanic and volcanic arc rocks; crustal contamination affects many western SNB plutons. In contrast, the eastern Sierra was dominated by voluminous recycling of the lithospheric mantle and lower crust, with minimal crustal contamination. Batholith-wide shifts in {delta}18O occur between pulses of Cretaceous magmatism that may be linked to tectonic reorganizations of magma sources. Within intrusive suites, {delta}18O may be unchanged (Tuolumne); increase (Sonora and Whitney); or decrease (Sequoia and John Muir) with time. These trends show stable long-lived sources, or those where recycling and contamination may increase or decrease with time. Overall, {delta}18O reveals diverse magma system behavior at a range of scales in the Sierran arc.

KEY WORDS: zircon; crustal growth; granitoids; supracrustal; magma systems; Sierra Nevada


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL BACKGROUND
 SAMPLE SELECTION, PREPARATION,...
 RESULTS
 DISCUSSION
 CONCLUDING REMARKS
 SUPPLEMENTARY DATA
 REFERENCES
 
The Sierra Nevada batholith (SNB, Fig. 1) has been intensely studied for decades to deduce the petrogenesis of individual plutons and intrusive suites and to understand the integrated history of the arc that formed it (see Bateman, 1992Go; Moore, 2000Go). Isotope studies have been particularly helpful in identifying significant crustal and mantle ‘reservoirs’ that contributed to the diversity of magmas of the SNB (Doe & Delevaux, 1973Go; Kistler & Peterman, 1973Go, 1978Go; DePaolo, 1980Go, 1981Go; Masi et al., 1981Go; Saleeby et al., 1987Go; Kistler, 1990Go; Kistler & Ross, 1990Go; Chen & Tilton, 1991Go; Kistler & Fleck, 1994Go; Lee et al., 2000Go; Ducea, 2001Go; Ratajeski et al., 2001Go; Wenner & Coleman, 2004Go; Lackey et al., 2005Go, 2006Go). Reservoirs include sub-arc mantle, subducted ocean crust, Phanerozoic sedimentary rocks, and Precambrian continental crust. The relative contributions of these reservoirs to SNB magmas are largely unknown. For example, studies of mantle xenoliths show considerable geochemical heterogeneity in the sub-arc mantle (Mukhopadhyay & Manton, 1994Go; Ducea, 1998Go, 2001Go, 2002Go; Lee et al., 2000Go), which suggests that the lithospheric mantle was heavily contaminated during subduction in the Jurassic and Triassic. Also, the amount of recycling of young rocks in the arc is uncertain. Recent experimental studies have reproduced characteristic major and trace element patterns of SNB granitoids by partially melting gabbroic compositions at pressure, temperature, and fluid conditions inferred for the magma-generating regions of the arc (Ratajeski et al., 2005Go; Sisson et al., 2005Go). These findings suggest significant remelting of mantle-derived basaltic sources at the base of the arc.


Figure 1
View larger version (46K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 1. Generalized geological map of the Sierra Nevada batholith. Study area (Fig. 3) is denoted by solid box, with Mariposa and Fresno 1° x 2° map sheets for reference. The initial 87Sr/86Sr = 0·706 line of Kistler (1990Go) is shown as a dashed line. Map after Jennings et al. (1977Go), Moore & Sisson (1987Go), Ross (1987Go), Saleeby et al. (1987Go), Bateman (1992Go), and Moore (2000Go).

 
Oxygen isotopes are well suited to clarifying ambiguity of magma source and contamination processes in the SNB. Low-temperature fractionation of 18O/16O, especially in the presence of water, imparts unambiguous {delta}18O ‘supracrustal’ signatures in sedimentary rocks that form at the Earth's surface, or in igneous and metamorphic rocks that are hydrothermally altered there. When such rocks melt following subduction or burial, the {delta}18O of these melts contrasts markedly with the restricted {delta}18O range of the primitive mantle (5·7 ± 0·3{per thousand}; Taylor & Sheppard, 1986Go).

Despite the utility of oxygen isotopes, their application to study the SNB has long been regarded with skepticism. Masi et al. (1981Go) conducted the first large-scale oxygen isotope study of the SNB. They found that, unlike the coeval Peninsular Ranges batholith (Fig. 2a; Taylor & Silver, 1978Go; Silver et al., 1979Go), the variation of {delta}18O with Sri in the SNB did not define a simple binary mixing pattern. Subsequent studies have shown an increasingly complex data array (Fig. 2b). Furthermore, hydrogen isotope analyses have revealed low {delta}D values, –85 and lower, indicating hydrothermal alteration of many samples.


Figure 2
View larger version (43K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 2. Oxygen and strontium isotopes in the Sierra Nevada and Peninsular Ranges batholiths. (a) The Peninsular Ranges trend relative to recognized reservoirs. (b) Sierran isotopic values by geographical area. Sources include: xenoliths (Ducea, 1998Go); altered ocean crust (McCulloch et al., 1980Go); altered greywacke (Magaritz & Taylor, 1976Go); Sierran metamorphic wall-rocks (Kistler & Peterman, 1973Go; DePaolo, 1981Go; Ross, 1983aGo, 1983bGo; Zeng et al., 2005Go). Data in (b) from Taylor & Silver (1978Go), Masi et al. (1981Go), Saleeby et al. (1987Go), Clemens-Knott (1992Go), Kistler (1993Go), Kistler & Fleck (1994Go), and Truschel (1996Go).

 
Our study provides fresh insight into the origin and evolution of magmas in the SNB. We present a comprehensive laser fluorination study of {delta}18O of zircon (Zrc) in the central Sierra (Fig. 3). Zircon was selected for the study because it is highly retentive of magmatic {delta}18O (Valley, 2003Go; Page et al., 2007Go). Prior study showed that analysis of {delta}18O(Zrc) in the SNB resolves magmatic {delta}18O in rocks that have undergone subsolidus exchange and hydrothermal alteration (Lackey et al., 2001Go, 2003Go, 2005Go, 2006Go). In this study, over 100 rocks were analyzed for both {delta}18O(WR) (where WR is whole-rock) and {delta}18O(Zrc); direct comparison of {delta}18O is presented to illustrate the magnitude and frequency of alteration of {delta}18O(WR). Also, we detail how whole-rock SiO2 content and {delta}18O(Zrc) allow calculation of magmatic {delta}18O(WR). These findings elucidate the scale and nature of magma systems in the arc.


Figure 3
View larger version (91K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 3. Geology of the central Sierra Nevada. Sample locations, plutons, the 0·706 line, and the Panthalassan–North American (PA/NA) break of Kistler (1990Go) are shown, along with the location of the transects of Figs 5 and 6. Details of Sequoia region samples 1–76 have been given by Chen & Moore (1982Go). Map after Jennings et al. (1977Go), Robinson & Kistler (1986Go), Kistler (1990Go), Bateman (1992Go), Clemens-Knott (1992Go), Saleeby & Busby (1993Go), and Moore (2000Go).

 

    GEOLOGICAL BACKGROUND
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL BACKGROUND
 SAMPLE SELECTION, PREPARATION,...
 RESULTS
 DISCUSSION
 CONCLUDING REMARKS
 SUPPLEMENTARY DATA
 REFERENCES
 
The central Sierra Nevada batholith is a collage of plutons and intrusive suites (Fig. 3) with pronounced spatial variation in composition, age, and isotope chemistry. Cretaceous granitoids are exposed over c. 35 000 km2 (Fig. 1). The estimated original thickness of granitic crust in the batholith is 35 km (Fliedner et al., 2000Go; Saleeby et al., 2003Go); thus the Sierran arc must have generated at least 1 x 106 km3 of granitic magma between 120 and 85 Myr ago. The batholith has long been recognized to have lateral (west to east) gradients in bulk composition (Moore, 1959Go; Bateman, 1992Go), pluton age (Kistler & Peterman, 1973Go, 1978Go; Stern et al., 1981Go; Chen & Moore, 1982Go), major and trace element geochemistry (Ague & Brimhall, 1988bGo; Bateman, 1992Go), and isotope ratios (Doe & Delevaux, 1973Go; Kistler & Peterman, 1973Go, 1978Go; DePaolo, 1980Go, 1981Go; Masi et al., 1981Go; Saleeby et al., 1987Go; Kistler, 1990Go; Kistler & Ross, 1990Go; Chen & Tilton, 1991Go; Kistler & Fleck, 1994Go; Ducea, 2001Go; Wenner & Coleman, 2004Go). Besides lateral variations in magma chemistry, recent geochemical studies of xenoliths (Fig. 2a) and deep crustal exposures suggest pronounced vertical zonation in the geochemistry of the batholith as well (Mukhopadhyay & Manton, 1994Go; Ducea, 1998Go, 2001Go, 2002Go; Lee et al., 2000Go; Saleeby et al., 2003Go; Lackey et al., 2005Go).

Various studies have interpreted transverse changes in SNB geochemistry as resulting from heterogeneity in the pre-batholith lithosphere and sub-arc mantle. The classic geochemical boundary in the SNB, the Sri = 0·706 line (Figs 1 and 3), delineates the inferred boundary between accreted Phanerozoic rocks in the western Sierra and Proterozoic continental lithosphere in the eastern Sierra (Kistler, 1990Go). A second boundary recognized inboard of the 0·706 line is the ‘PA/NA break’ of Kistler (1990Go). This break runs through the central Sierra and delineates a cryptic pre-batholithic domain boundary between oceanic ‘Panthalassan’ (Kistler, 1990Go) lithosphere to the west and North American lithosphere to the east (Fig. 3). The break is based on the location of the 0·706 line, variations of {delta}18O(WR), and the location of offset pendants or plutons with similar lithology, age, and chemistry (Kistler, 1990Go, 1993Go).

Intrusive suites in the Sierra are the ‘microcosms’ of the batholith. Suites are often defined from field relations (Bateman, 1992Go), geochemical, and isotopic zonation (Bateman & Chappell, 1979Go; Kistler et al., 1986Go; Bateman, 1992Go), and similar ages or age progressions in member plutons, such as inward decreasing age in zoned plutonic complexes (Chen & Moore, 1982Go; Coleman et al., 2004Go). Because trends in geochemistry and age in some suites mimic trends of those same parameters across the SNB, we conducted detailed sampling in several suites to determine how their magmatic systems relate to the batholith as a whole.


    SAMPLE SELECTION, PREPARATION, AND ANALYSIS
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL BACKGROUND
 SAMPLE SELECTION, PREPARATION,...
 RESULTS
 DISCUSSION
 CONCLUDING REMARKS
 SUPPLEMENTARY DATA
 REFERENCES
 
Zircon was concentrated by standard crushing, density, and magnetic separation techniques. Least magnetic fractions were concentrated with a Frantz separator. Most samples were originally collected for geochronology (Chen & Moore, 1979Go, 1982Go; Ratajeski et al., 2001Go) and thermochronology (House et al., 1997Go, 1998Go, 2001Go; Stockli et al., 2003Go) studies in the Sierra Nevada and White Mountains. The first author collected samples with ‘1S’ and ‘3S’ prefixes in 2001 and 2003, respectively (see Table 1). Electronic Appendix 1 (available for downloading at http://www.petrology.oxfordjournals.org/) provides details of location, lithological description, and researcher for all samples. Detailed discussion of sampling and mineral separation techniques for the earlier studies is given in the original papers.


View this table:
[in this window]
[in a new window]

 
Table 1: Oxygen isotope ratios of zircon, quartz, and whole-rock in the central Sierra Nevada

 
Zircon was prepared for {delta}18O analysis by sequential acid leaching in HNO3, HF, and HCl, as described by Lackey et al. (2006Go). HF preferentially dissolves radiation-damaged domains from zircon, which are subject to alteration of {delta}18O, but does not affect crystalline zircon which preserves pristine {delta}18O (King & Valley, 2001Go; Valley, 2003Go). Handpicking under a binocular microscope further purified zircon separates. For each analysis, comprising tens to hundreds of crystals, zircon was powdered using a boron carbide mortar and pestle before laser fluorination. Powdering creates uniform size and increases fluorination efficiency.

Mineral separates and whole-rock powders were analyzed in the University of Wisconsin Stable Isotope Laboratory by laser fluorination. Isotope ratios were measured on a dual inlet gas source Finnigan MAT 251 mass spectrometer. Whole-rock powders (2–2·5 mg) were analyzed with an airlock sample chamber (Spicuzza et al., 1998Go). Oxygen was liberated from silicates with BrF5 using a 30W CO2 laser ({lambda} = 10·6 µm), then purified cryogenically, passed through hot Hg to remove any residual F2, and finally converted to CO2 with a hot carbon rod (Valley et al., 1995Go). All analyses were standardized daily based on four or more analyses of UWG-2, Gore Mountain garnet standard, and sample {delta}18O values were corrected to the long-term accepted value of 5·80{per thousand} for UWG-2. The average raw {delta}18O of UWG-2 for 40 days of analyses (n = 217) in this study is 5·71 ± 0·13{per thousand} [1SD; standard error (1{sigma} = 1SD/(n2) is ±0·009{per thousand}]. Average day-to-day precision of UWG-2 is ± 0·06{per thousand}; daily corrections average 0·08{per thousand}. Average {delta}18O of NBS-28 run on 7 days was 9·45 ± 0·13{per thousand} after correction, n = 20, and {Delta}(NBS-28–UWG-2) averaged 3·65 ± 0·12{per thousand}, in excellent agreement with the UW laboratory long-term {Delta}(NBS-28–UWG-2) of 3·70 ± 0·06{per thousand} (1SD) (Valley et al., 1995Go). During 9 days of whole-rock analyses using the airlock chamber, UWG-2 averaged 5·67 ± 0·11{per thousand} (n = 41); replicate analyses of whole-rock powders averaged ±0·11{per thousand}. Replicate analyses give uncertainties of less than ±0·04% for zircon and ±0·09% for quartz.

Major elements and selected trace elements (Rb, Sr, Y, Nb, Zr) were determined for 96 samples by X-ray fluorescence at XRAL Laboratories, Canada. These data and discussion of analytical techniques are provided in Electronic Appendix 2.

Comparison of the {delta}18O values of mineral pairs with equilibrium fractionation factors is used in several sections of the paper to test if minerals have exchanged oxygen, and the degree to which they approach isotopic equilibrium. Isotopic fractionation factors are expressed as


Formula

Equilibrium fractionation factors are calculated for a particular temperature using the expression


Formula

where T is temperature (K) and A is an experimentally or empirically determined coefficient. The following A factors are used for oxygen isotopes: quartz–zircon = 2·64 (Valley et al., 2003Go); quartz–plagioclase = 1· 20 (Clayton et al., 1989Go); quartz–K-feldspar = 1· 00 (Clayton et al., 1989Go); quartz–biotite = 2·16 (Chacko et al., 2001Go); quartz–hornblende = 3·01 by combining the fractionation factors quartz–almandine (Valley et al., 2003Go) and garnet–hornblende (Kohn & Valley, 1998Go).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL BACKGROUND
 SAMPLE SELECTION, PREPARATION,...
 RESULTS
 DISCUSSION
 CONCLUDING REMARKS
 SUPPLEMENTARY DATA
 REFERENCES
 
Oxygen isotope ratios
Zircon
Zircon {delta}18O in the study area ranges from 4·2 to 9·2{per thousand} (Table 1). Variation of {delta}18O is presented both as a color-coded map of the central Sierra (Fig. 4), and as two east–west transects (Figs 5 and 6) that divide the data approximately between the Fresno (Fig. 5; 36–37°N) and Mariposa (Fig. 6; 37–38°) 1/250 000 sheets (Fig. 1).


Figure 4
View larger version (71K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 4. {delta}18O(Zrc) in the central Sierra Nevada. Four prominent high- and low-{delta}18O belts (Zones 1–4) are recognized (see inset). Zoning of {delta}18O in the Mount Givens pluton shown as {delta}18O zones along its west side. The significance of the Eagle Peak pluton, Boyden Cave, and Mineral King pendants is discussed in the text. P, peraluminous pluton; IBB-2, Intrabatholithic Break; 2. John Muir Suite includes: Kin, Inconsolable Granodiorite; Klk, Lamark Granodiorite; Kle, Lake Edison Granodiorite; Krv, Round Valley Peak Granodiorite; Kmo, Mono Creek Granite.

 

Figure 5
View larger version (32K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 5. Fresno sheet transect showing (a) {delta}18O(Zrc). (b) Measured and calculated {delta}18O(WR). (c) Initial 87Sr/86Sr. (Note location of 0·706 line, Foothills Break, and PA/NA break.) Error bars on calculated {delta}18O are 1SD and reflect greater uncertainty of values of samples for which SiO2 was estimated. Published data are from sources listed in Fig. 2 and Chen & Tilton (1991Go), Sisson et al. (1996Go), Coleman & Glazner (1997Go), Wenner & Coleman (2004Go), and Lackey et al. (2005Go, 2006Go).

 

Figure 6
View larger version (33K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 6. Mariposa sheet transect showing (a) {delta}18O(Zrc), (b) measured and calculated {delta}18O(WR), and (c) initial 87Sr/86Sr. Major features and boundaries are as in Fig. 3. Published data are from sources listed in Fig. 2 and Hill et al. (1988Go), Bateman et al. (1991Go), Coleman et al. (1992Go), Coleman & Glazner (1997Go), Ratajeski et al. (2001Go), Ernst & Rumble (2003Go), Ernst et al. (2003Go), Wenner & Coleman (2004Go), and Lackey et al. (2005Go, 2006Go).

 
Values span most of the terrestrial range for {delta}18O of igneous zircon in plutonic rocks (Valley et al., 2005Go); 94% are above the mantle value for zircon (5·3 ± 0·3{per thousand}, 1SD; Valley et al., 1998Go), and the rest are in or slightly below the mantle range. Comparison of {delta}18O variation between different grain-size fractions of zircon (see Table 1) reveals no consistent pattern of size versus {delta}18O within samples (Lackey, 2005Go). Intra-pluton {delta}18O variations typically are small, <0·3{per thousand} (Table 1), and comparable with the range of {delta}18O(Zrc) from volcanic systems with relatively simple zircon crystallization histories. For example, the Bishop Tuff varies by only 0·21{per thousand} between early and late erupted units (Bindeman & Valley, 2002Go). The Mount Givens, Bass Lake, and Giant Forest Plutons are zoned in {delta}18O (Figs 3 and 4).

Whole-rocks
Measured {delta}18O(WR) across the central Sierra varies from 6 to 11·7{per thousand}, with one value as low as –0·5{per thousand} (Table 1). The range of {delta}18O(WR) values is larger than previously reported in the central Sierra (Masi et al., 1981Go; Kistler, 1990Go). Overall, this study adds considerably to the coverage of {delta}18O(WR) in the Fresno sheet (Fig. 5b), and doubles the data for the Mariposa sheet (Fig. 6b), which increases the known range of {delta}18O(WR).

Quartz
Values of {delta}18O(Qtz) were analyzed from rocks throughout the central Sierra. The range of {delta}18O(Qtz) is 8·8–12·0{per thousand}. Distribution of oxygen isotopes between coexisting quartz and zircon yields apparent temperatures of 575–650°C, which are significantly lower than are expected for magmatic crystallization (850–650°, Fig. 7). Similar {Delta}18O(Qtz–Zrc) disequilibrium is recognized in the Idaho batholith (King & Valley, 2001Go) and is explained by closed-system resetting of {delta}18O(Qtz) during sub-solidus cooling, owing to relatively fast oxygen diffusion in quartz (Farver & Yund, 1991Go; Sharp et al., 1991Go), which yields a lower closure temperature (Dodson, 1973Go) for oxygen exchange in quartz than in zircon.


Figure 7
View larger version (33K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 7. {delta}{delta} plot for quartz and zircon. Fractionations yield apparent temperatures that are lower than expected for high-temperature magmatic equilibrium. Isotherms calculated from Valley et al. (2003Go).

 
Oxygen isotope mapping of the Central Sierra Nevada
High sampling density, homogeneity of {delta}18O in most plutons, and thorough mapping allow us to present the results as a detailed isotopic map of the central SNB (Fig. 4). This representation both shows regional {delta}18O variations and resolves isotopic zoning of intrusive suites with detail beyond simple contouring of {delta}18O. Except where {delta}18O varies by >0·5{per thousand}, values were averaged for plutons having multiple samples. An overarching pattern emerges with four north–south running belts of rocks having {delta}18O(Zrc) that is greater or less than 6·5{per thousand} (inset, Fig. 4). Nonetheless, there is substantial heterogeneity within the belts, including on-strike changes of {delta}18O (see ‘Regional {delta}18O belts’ below).


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL BACKGROUND
 SAMPLE SELECTION, PREPARATION,...
 RESULTS
 DISCUSSION
 CONCLUDING REMARKS
 SUPPLEMENTARY DATA
 REFERENCES
 
Relating {delta}18O(Zrc) to {delta}18O(WR)
Correlation with SiO2
To use zircon as a benchmark for evaluating the degree of contamination or sub-solidus alteration, the fractionation of {delta}18O(Zrc) and {delta}18O(WR) at magmatic temperatures must be calibrated. Overall, measured {delta}18O(WR) correlates with {delta}18O(Zrc), but {delta}18O(WR) is more variable and 0·5–3·0{per thousand} higher (Fig. 8a; Table 1). Values of {delta}18O(WR) positively correlate with SiO2 content and {delta}18O(Zrc) has a slight positive correlation. The high variance in Fig. 8a reflects differences of magmatic {delta}18O throughout the central Sierra Nevada.


Figure 8
View larger version (22K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 8. Correlation of {delta}18O and silica content. (a) {delta}18O(Zrc) and {delta}18O(WR) vs SiO2 for all samples. (b) Pattern of {delta}18O(Zrc) and {delta}18O(WR) vs SiO2 in the Tuolumne Intrusive Suite. Trends of fixed {delta}18O(Zrc) but increasing {delta}18O(WR) with SiO2 should be noted. One anomalously high {delta}18O(Zrc) value was omitted from the linear best fit of the data. (c) Calculated {delta}18O(WR–Zrc) vs SiO2 using a compilation of Sierran rocks with published SiO2 and modal data (see Electronic Appendix 3). (d) Direct comparison of calculated and measured {delta}18O(WR). The grey field bounding the equilibrium line is 3SD of analytical uncertainty wide.

 
Values of {delta}18O(Zrc) are independent of wt % SiO2 in the Tuolumne Intrusive Suite, with one outlier, whereas {delta}18O(WR) is positively correlated and varies by over 2·0{per thousand} from 57 to 76 wt % SiO2 (Fig. 8a). Mafic magmas undergoing closed-system differentiation typically see {delta}18O(WR) increase by 0·5{per thousand} for each 10% increase in SiO2 (Taylor & Sheppard, 1986Go); however, results from the Tuolumne suite indicate that {delta}18O(WR) in calc-alkaline systems increases by ~1· 0{per thousand} for each 10% increase of SiO2. In general, increasing SiO2 in calc-alkaline igneous rocks correlates with increasing modal abundance of high-{delta}18O minerals such as quartz and feldspar, and with decreasing amounts of lower {delta}18O minerals such as amphibole, biotite, and magnetite. The Tuolumne suite exemplifies how differentiation creates variability of {delta}18O(WR) in a pluton or co-genetic suites, even if minerals have constant {delta}18O, and why zircon best tracks magmatic {delta}18O.

Various approaches can be used to predict equilibrium fractionations ({Delta}18O) between igneous rocks, melts, and their various constituent minerals. Valley et al. (1994Go) used correlation of {Delta}18O(WR–Zrc) in four model rock compositions of varying SiO2 content to derive an approximate relationship of {Delta}18O(WR–Zrc) as a function of wt % SiO2. By assuming all minerals were in equilibrium at 800°C, they calculated whole-rock {delta}18O values that were mass balanced for the mode of each mineral in a rock. Two recent studies have used normative mineralogy to calculate {Delta}18O(melt–mineral) values (Appora et al., 2003Go; Zhao & Zheng, 2003Go). This approach uses the principles of the Garlick Index (Garlick, 1966Go), which considers the degree of polymerization in a magma. Felsic magmas contain greater amounts of Al and Si, stabilizing minerals whose crystal chemistry favors stronger oxygen bonding, and therefore higher 18O/16O ratios. The disadvantage of using normative mineralogy includes omission of hydrous phases, poorly known fractionation factors for some normative phases, and mismatch between normative mineralogy and actual rock modes.

{Delta}18O(WR–Zrc) in the Sierra Nevada
Modeling of {Delta}18O(WR–Zrc) uses a compilation of 297 SNB rocks ranging from gabbro to granite that the US Geological Survey analyzed for mineral modes and SiO2 (see Electronic Appendix 3). The silica content of each rock is directly related to {Delta}18O(WR–Zrc), which is calculated from the modal abundance and the equilibrium fractionation factors for quartz, K-feldspar, plagioclase, biotite, hornblende, and zircon (see ‘Methods’). Solidus temperatures appropriate for magma compositions ranging from gabbro to high-silica granite were chosen (Fig. 8c), with the justification that both solidus and zircon saturation temperatures decrease with increasing SiO2, such that {Delta}18O(WR–Zrc) represents different temperatures across a range of compositions. Recent Ti-in-zircon thermometry of a representative suite of these rocks confirms higher average Ti-in-zircon temperatures in diorites and gabbros than in peraluminous granites (Fu et al., 2008Go). Solidus temperatures were assumed to have decreased linearly from 900°C (gabbro, PH2O = 3·0 kbar, 45 wt % SiO2) to 650°C (high-silica granite, PH2O = 3·0 kbar, 78 wt % SiO2). Calculated {Delta}18O(WR–Zrc) values show a strong (r2 = 0·96) positive correlation with SiO2 (Fig. 8c). A linear fit yields the following equation:


Formula 1

(1)
Equation (1) allows {delta}18O(WR) values to be calculated based on measured {delta}18O(Zrc) and whole-rock silica content. Possible sources of error on calculated {delta}18O come from uncertainty in SiO2 content and {delta}18O(Zrc). If SiO2 and {delta}18O(Zrc) uncertainties are 0·5% and 0·1{per thousand}, respectively, the error on the calculated {delta}18O(WR) is ±0·12{per thousand}.

Comparison of measured and calculated {delta}18O(WR) values shows that many Sierran rocks deviate from the predicted values (Fig. 8d). To estimate the number of rocks that are hydrothermally altered, the average analytical precision (0·12{per thousand}, 1SD, Fig. 8d) of {delta}18O(WR) analyses was used as a cutoff. In cases where calculated and measured {delta}18O(WR) values vary by greater than ±0·12{per thousand} of the calculated WR uncertainty and 3SD of analytical precision (>0·36{per thousand}), the {delta}18O(WR) is considered altered. By these criteria, approximately 28% of the {delta}18O(WR) values in the central Sierra are altered. Many values deviating by less than 0·36{per thousand} are probably altered as well; thus 28% is a minimum estimate. Textural indicators of alteration include widespread sericitization of feldspar and chlorite and epidote veins and overgrowths.

Inspection of Fig. 8d shows that most of the measured {delta}18O values showing alteration are higher than their expected magmatic value. Closed-system sub-solidus exchange in a rock will result in shifts in {delta}18O among the minerals in a rock relative to one another, but the overall {delta}18O(WR) will remain unchanged. Alteration causing a shift in {delta}18O(WR) typically results from sub-solidus exchange with externally derived fluids, possibly at a wide range of temperatures and fluid compositions. The positive shift of {delta}18O(WR) observed for most samples is inconsistent with infiltration of meteoric water; however, three samples with low-{delta}18O(WR) values, 2–5{per thousand}, confirm some sub-solidus infiltration of heated meteoric water in the eastern Sierra (Fig. 8d). The samples come from areas of the eastern Sierra where igneous barometry indicates shallow (<2 kbar) depths of crystallization (Ague & Brimhall, 1988aGo; Ague, 1997Go). Hydrogen isotope data from these same areas confirm infiltration of meteoric water as {delta}D(WR) values are commonly less than –85{per thousand} (Godfrey, 1962Go; Masi et al., 1981Go). In addition, samples from the peraluminous Grant Grove pluton have measured {delta}18O(WR) values that are 1–2{per thousand} greater than predicted (Fig. 8d). These samples are from the margin of the pluton and record contamination by high-{delta}18O wall-rocks after zircon crystallized (Lackey et al., 2006Go).

Although generally slight, alteration is enough to obscure isotopic zoning in plutons or intrusive suites and can hide detail such as {delta}18O discontinuities. Figures 5a, b and 6a, b show direct comparison of measured and calculated {delta}18O(WR) across the SNB. In these two transects, additional calculated {delta}18O values are shown, although they are less precise because SiO2 was not measured for those samples and was estimated from published SiO2 data (see Electronic Appendix 4). This direct comparison along transects shows local zones of alteration. For instance, the White Mountains have both low and variable measured {delta}18O, and calculated values that are higher and less variable, suggesting additional areas of meteoric water infiltration, one of which is confirmed below.

Two case studies distinguishing magmatic and alteration history
The Barcroft pluton
In the Jurassic Barcroft pluton, northern White Mountains (Fig. 3), pervasive secondary mineralization and a wide range of {delta}18O(WR) values (6–10{per thousand}) indicate hydrothermal alteration of the pluton (Ernst & Rumble, 2003Go). Three widely spaced zircon samples in the pluton show that magmatic {delta}18O was in fact very homogeneous: 7·60{per thousand} ± 0·08 (Fig. 9a). Using the zircon data from this study and measured SiO2 of Ernst & Rumble (2003Go), magmatic {delta}18O(WR) is calculated to vary less than 1{per thousand} in the pluton. All measured {delta}18O(WR) values deviate over 1{per thousand} from calculated values, and three of four measured values are lower than calculated (Fig. 9a). Although none of the {delta}18O(WR) values that Ernst & Rumble (2003Go) reported is exceptionally low, the tendency toward lower measured values confirms some meteoric water infiltration during hydrothermal alteration of the Barcroft pluton.


Figure 9
View larger version (18K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 9. Alteration in the Barcroft pluton and Saddlebag Lake meta-rhyolite. (a) SiO2 vs {delta}18O of zircon and whole-rock (measured and calculated) values in the Jurassic Barcroft pluton. The small range of {delta}18O(Zrc) and calculated {delta}18O(WR) relative to measured values should be noted. Barcroft {delta}18O(WR) and SiO2 values from Ernst & Rumble (2003Go). (b) {delta}18O of quartz and whole-rock in the Saddlebag Lake meta-rhyolite are altered and higher than expected according to {delta}18O(Zrc). Disequilibrium {delta}18O between metamorphic garnet and zircon indicates that garnet crystallized after alteration. Calculated whole-rock and quartz {delta}18O relative to both garnet and zircon emphasize systematic differences between the meta-rhyolite and its precursor.

 
Saddlebag Lake meta-rhyolite
A Triassic rhyolite in the Saddlebag Lake pendant, east of the Tuolumne Suite (Fig. 3), has a complex magmatic, alteration, and metamorphic history. The rhyolite has a typical {delta}18O(Zrc) value of 6·2{per thousand}. In contrast, the {delta}18O values of the whole-rock and quartz phenocrysts are 11· 6{per thousand} and 12·6{per thousand}, respectively; both are elevated well above typical igneous values and confirm alteration (Fig. 9b). If the original SiO2 content of the rhyolite was 70 wt %, then the calculated {delta}18O(WR) is c. 8·3{per thousand} (Fig. 9b), much lower than the measured value.

Metamorphic garnet in the meta-rhyolite yields another constraint on its alteration history. Because garnet has slow oxygen diffusion (Couglan, 1990Go; Vielzeuf et al., 2005Go), like zircon, it records the {delta}18O of a rock at the time it crystallizes. For Fe–Mn garnets, equilibrium fractionation with zircon is ~0{per thousand} at magmatic temperatures (Valley et al., 2003Go). Therefore, the 9{per thousand} garnet is not equilibrated with the 6·2{per thousand} zircon. As such, a chronology of events affecting the rhyolite is: (1) the rhyolite was erupted; (2) hydrothermal alteration elevated {delta}18O(WR); (3) Cretaceous metamorphism of the altered rhyolite led to garnet crystallization and re-equilibration of all phenocrysts except zircon.

Regional {delta}18O belts
Zone 1: Foothills belt
Mantle-like {delta}18O(Zrc) values occur in and around mafic ring dike complexes in the western Sierran foothills (Fig. 4, inset). This low-{delta}18O domain abruptly changes to high-{delta}18O rocks 20 km into the batholith (Fig. 5a). The extent of interaction between basement rock and magmas in this belt is unclear, although detailed O, Sr, and Nd isotope investigation of the ring dike complexes indicates a nonradiogenic, depleted mantle source with melting at pluton–wall-rock contacts causing some localized contamination (Clemens-Knott, 1992Go). Mantle-like {delta}18O(Zrc) values appear to rule out widespread contamination; however, the average {delta}18O(WR) of the wall-rocks, a mélange of the Kings–Kaweah ophiolite belt (Saleeby, 1992Go; Saleeby & Busby, 1993Go), is 7·1 ± 1·3{per thousand} (Lackey et al., 2006Go). Thus, there is minimal magma–wall-rock isotopic contrast, which would obscure source or contamination contributions from melting of ophiolitic rocks.

Zone 2: Fine Gold–Sequoia belt
In the Fine Gold Suite {delta}18O decreases from west to east, and is negatively correlated with Sri (Figs 4 and 6a); rocks with high {delta}18O(Zrc), up to 8·7{per thousand}, have low Sri (~0·704). The pattern of these isotope ratios requires substantial input of hydrothermally altered ocean crust or volcanic arc sediments in the western parts of the Fine Gold suite (e.g. Fig. 2a). In ocean crust, circulation of heated water increases the {delta}18O(WR) of upper pillow basalt and sheeted dike sections considerably from mantle values (altered ~10{per thousand}, unaltered 5·4–6·0{per thousand}; McCulloch et al., 1980Go; Muehlenbachs, 1998Go). Low Sri values result because 87Sr/86Sr is only modestly elevated from mantle values by alteration (McCulloch et al., 1980Go; Muehlenbachs, 1998Go). Furthermore, low Sri indicates that the rocks are relatively young because insufficient time has passed for significant ingrowth of 87Sr. Similar {delta}18O and Sri would be expected for arc volcanic rocks, such as greywacke, that have been hydrothermally altered (Magaritz & Taylor, 1976Go).

East of the Fine Gold rocks, {delta}18O and Sri show positive correlation in a series of plutons that lead up to the PA/NA break (Figs 3 and 6). In the Fresno Sheet, rocks of Zone 2 initially have high {delta}18O and sub-0·706 Sri values (Fig. 5; 38–55 km), indicating sources with juvenile supracrustal rock, but less than the Fine Gold Suite. Eastward increasing {delta}18O and Sri in the belt, into the Sequoia region, mimics the pattern in plutons east of the Fine Gold Suite, with the PA/NA break again marking the transition to lower {delta}18O.

Overall, the geochemistry of the plutons in Zone 2 reflects the influence of the exposed wall-rocks. For instance, rocks of the Fine Gold Suite with the highest {delta}18O are in direct contact with volcaniclastic rocks of the Jurassic–Triassic island arc belt and outboard accretionary sequence of the Sierran foothills (Fig. 4; Snow & Scherer, 2006Go). These rocks have the exact isotopic composition, with low Sri (0·704–0·707; Kistler & Peterman, 1973Go) and high {delta}18O(WR) (8–14{per thousand}; Böhlke & Kistler, 1986Go), to produce the local contamination patterns in the Fine Gold Suite. Negative correspondence of Sri and {delta}18O on the east side of the Fine Gold Suite suggests differences of magma source and contamination (Fig. 6). Permo-Triassic chert–argillite mélange of the Calaveras Complex is in fault contact with the Jurassic–Triassic metamorphic rocks belt near the 0·706 line (Fig. 4). Pendants and septa show southeastern continuation of both metamorphic belts to their termination near the boundary of the Mariposa and Fresno sheets (Fig. 4). Therefore the Bass Lake Tonalite intrudes across the fault contact into both metamorphic belts. Because Calaveras rocks contain greater amounts of continental sediments than rocks to the west, eastward increases of Sri in the Bass Lake Tonalite indicate greater proportions of Calaveras Complex rock in magma sources or as a contaminant (Fig. 6a and b). Decreasing {delta}18O may reflect greater depleted mantle contributions to the east or lower degrees of partial melting of refractory chert and argillite-rich Calaveras rocks. Minor partial melting, especially involving labile phases such as micas, can radically shift magmatic Sri without markedly changing {delta}18O. Thus, the 0·706 line may reflect crustal architecture as well as wall-rock fertility.

Another example of wall-rock influence on magma chemistry in the Zone 2 high-{delta}18O belt concerns the Kings Sequence, which contaminates magmas throughout the Sequoia region (Fig. 4; Ague & Brimhall, 1988bGo; Kistler, 1990Go; Chen & Tilton, 1991Go). Kings Sequence rocks and some from the Calaveras Complex host the Strongly Contaminated and Reduced plutons of Ague & Brimhall (1988bGo). These plutons, including several that are peraluminous (Fig. 3), intrude Kings Sequence rocks (Fig. 3). Ague & Brimhall (1988bGo) showed increased Fe/(Fe + Mg) in biotite in these rocks and proposed that graphite-rich Kings Sequence rocks led to magmas being relatively reduced, 0–2 log fO2 units below QFM (quartz–fayalite–magnetite). High {delta}18O(Zrc) > 7·5{per thousand}, but relatively low Sri (0·704–0·707), in these plutons shows that they contain high proportions of young supracrustal rocks in their sources; comparison of pluton margin and interior geochemistry and mineralogy reveals localized melting and contamination by the Kings Sequence (Lackey et al., 2006Go).

Zone 3: Sierra Crest and eastern Sierra
East of Zones 1 and 2 there is a monotonic decrease in {delta}18O toward the Sierra Crest (Fig. 5, 55–104 km; Fig. 6, 85–130 km). Many of the values in the eastern Sierra fall in the range of mantle {delta}18O(Zrc) values (5·3 ± 0·3{per thousand}, Valley et al., 1998Go). Calculated {delta}18O(WR) values for granodiorites and granites in this belt are 7–8{per thousand}. These relatively low values agree with studies proposing that significant portions of the late Cretaceous SNB were generated directly from melting of enriched (Sri = 0·706; {varepsilon}Nd = –4·5) lithospheric mantle beneath the eastern Sierra (Coleman et al., 1992Go; Coleman & Glazner, 1997Go). Although pendants of metavolcanic and metasedimentary rock are exposed in the eastern Sierra (Fig. 3), {delta}18O values do not indicate significant magmatic incorporation of these or other crustal rocks.

Zone 4: Owens Valley–White Mts. belt
Triassic, Jurassic, and Cretaceous plutons in the Owens Valley, Long Valley re-entrant, and White Mountains typically have {delta}18O(Zrc) > 6·5{per thousand} (Fig. 4). Values of Sri in these rocks often exceed > 0·709 (Fig. 6c), indicating magma sources containing aged sediments or crust. Wall-rock contamination effects are pronounced in many of the solitary plutons in the White Mountains that intrude Proterozoic to Paleozoic continental-derived sedimentary rocks. Values of Sri and {varepsilon}Nd at the margin and cupola regions of the plutons are more crustally influenced (Sri > 0·710, {varepsilon}Nd < –9·0) than their interiors (Sri < 0·706, {varepsilon}Nd > –1· 5; Ernst et al., 2003Go), indicating ‘veneers’ of contamination.

Crustal boundaries
Projection of lateral changes of {delta}18O in the SNB shows mostly gradual changes (Figs 5 and 6); however, localized, sharp, step-function changes of {delta}18O suggest steeply dipping trans-crustal breaks. Abrupt lateral shifts in {delta}18O suggest that magmas generated on either side of a break inherit the {delta}18O signature of the terrane in which they form, with a steep geometry of the terrane boundary preventing significant averaging of values across the break (see Peck et al., 2004Go). Below we discuss two significant lithospheric breaks that are resolved by {delta}18O(Zrc), and the position of the PA/NA break.

Foothills Break
The abrupt increase of 1· 5{per thousand} in {delta}18O(Zrc) that separates Zones 1 and 2 in the western Fresno sheet (20 km, Fig. 5a), is not discernible from {delta}18O(WR) or Sri (Fig. 5), but may be the geochemical expression of the Foothills suture (Saleeby, 1992Go), which is poorly exposed. Saleeby (1992Go) proposed northward projection of the Foothills suture into the Bass Lake Tonalite, but the limited extent of Zone 1 low-{delta}18O rocks to the north means that the break terminates before the tonalite. Greater wall-rock heterogeneity to the north and different magmatic style may conceal the break.

Eastern Sierra Break
The sharp 1{per thousand} west-to-east increase of {delta}18O(Zrc) defining the boundary between Zones 3 and 4 (Figs 4 and 6a, 128 km), defines a feature in the eastern Sierra and Owens and Long Valleys that has long been inferred to be a major crustal boundary (Saleeby et al., 1986Go; Kistler, 1993Go). This discontinuity corresponds to Intrabatholithic Break 3 (IBB3) of Saleeby et al. (1986Go), and the Tinemaha Fault of Stevens et al. (1997Go). This boundary was originally recognized by displacement of Jurassic and Triassic rocks along a lineament that intersects a major right-lateral offset in the 0·706 line in Long Valley (Fig. 4). The boundary is also inferred from structural and stratigraphic relationships where SW-trending Neoproterozoic to Triassic sedimentary shelf sequences in the White–Inyo Mountains are juxtaposed against the NW-trending SNB (Stevens et al., 1997Go). A steeply dipping geometry for this break favors tectonic truncation as the mode of origin, and the expression of the {delta}18O step in Triassic and younger plutons means that the boundary is at least pre-Jurassic. In addition, recent dating of dikes offset in the Owens Valley indicates that the 65 km of dextral offset along this break occurred since 83 Ma (Kylander-Clark et al., 2005Go).

PA/NA Break
Values of {delta}18O(Zrc) decrease gradually at the PA/NA break (Figs 5a and 6a) and generally correspond to its originally proposed location (Fig. 4). The lack of a sharp step in {delta}18O suggests a diffuse transition, rather than a trans-crustal break; yet changes in isotopic gradients, and the preponderance of data from the high- to low-{delta}18O zones that define it, indicate that this is a major lithospheric discontinuity. Potentially, the integrated effects of lithospheric mantle differences as well as an upper crustal overprint have blurred the break; however, such complex controls cannot be evaluated with certainty. Although the PA/NA break remains enigmatic, {delta}18O(Zrc) values show that some adjustments of the position of the PA/NA break are needed.

In the Fresno sheet, we propose two refinements. First, the trace of the break appears to be west of the Mineral King pendant, rather than through it (Fig. 4). Whereas Kistler (1993Go) proposed linking the PA/NA break to faults in the Mineral King pendant, the low {delta}18O(Zrc) of the pluton to the west of the pendant (granodiorite of Castle Creek) suggests that the break is ~10 km to the west (Fig. 4). Whereas the Mineral King pendant has lithologies similar to the Kings Sequence, volcanic rocks in the pendant have affinities with eastern terranes (Saleeby & Busby, 1993Go). The second adjustment is a 5 km eastward shift of the break at the Boyden Cave pendant (Fig. 4). The pendant clearly juxtaposes rocks from the Kings Sequence and the Goddard terrane to the east; however, refinement based on {delta}18O may better reflect the deep orientation.

North of the Mount Givens pluton, in the Mariposa sheet, the original PA/NA break strikes NW and joins the 0·706 line; however, high- and low-{delta}18O domains are continuously juxtaposed to the north (Fig. 4). If the bulbous westward protrusions of the Mount Givens pluton and Tuolumne Intrusive Suite are ignored as distortional effects (Fig. 4), the {delta}18O-defined PA/NA break projects directly through Yosemite to the offset 0·706 line in the Sonora Intrusive Suite (Fig. 4). Whereas Kistler (1993Go) favored the joining of the PA/NA break and 0·706 line to include the Melones fault as a northern extension of the PA/NA break, Seleeby et al. (1986Go) proposed that a straight intrabatholithic break, labeled IBB2 in Fig. 4, continues from the PA/NA trace into Yosemite (Kistler, 1993Go). The pattern of {delta}18O supports the northward trend into Yosemite.

Direct comparison of O, Sr, and Pb isotopes
Chen & Tilton (1991Go) reported Sr and Pb isotope ratios on a subset of samples analyzed in this study from the Fresno sheet (see Chen samples in Fig. 3). Therefore, correlation of O, Sr, and Pb isotope ratios can be evaluated for a representative suite of rocks across the SNB (Fig. 10). Plots depict Sr and Pb isotope variation relative to {delta}18O of zircon (Fig. 10a–d) and whole-rocks (Fig. 10e–h) to evaluate the effect of differing silica among samples on isotopic correspondence. The samples are divided into western and eastern samples; western samples have Sri < 0·706; eastern samples occur east of the PA/NA break. A third group of samples are from the Sequoia region (Fig. 3).

Oxygen and strontium isotopes
Co-variation of {delta}18O(Zrc) with Sri is geographically distinct for the western and eastern Sierra (Fig. 10a). Calculated {delta}18O(WR) comparisons yield more scatter but similar trends (Fig. 10e). The array of {delta}18O and Sri for western samples intersects a mantle value at the low-Sri end (Clemens-Knott, 1992Go) and trends toward an evolved (high-Sri) reservoir; overall there is resemblance to the Peninsular Ranges trend (e.g. Fig. 2b). One sample from the margin of the Grant Grove pluton lies well off the western trend, recording contamination typical of the Sequoia region (Fig. 10a).


Figure 10
View larger version (33K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 10. Correspondence of {delta}18O to Sr and Pb isotope ratios across the Fresno sheet. (a) {delta}18O(Zrc) vs 87Sr/86Sr initial. (b–d) {delta}18O(Zrc) vs initial Pb isotope ratios. (e) {delta}18O calculated whole-rock (WR) vs 87Sr/86Sr initial. (f–h) {delta}18O calculated WR vs initial Pb isotope ratios. Trend lines represent least-squares best fits to the data. The high-Sri sample from the contaminated margin of the Grant Grove pluton was omitted from the best fit of western samples; Sequoia samples were not fitted with a trendline because of overall poor correlation.

 
Eastern samples also exhibit a distinct correlation, with a less steep positive slope, forming a flat array of high Sri values throughout the eastern Sierra. Both low- and high-{delta}18O plutons indicate considerable heterogeneity in the crustal content of these magmas, with low {delta}18O values reflecting old lithospheric mantle (e.g. Coleman & Glazner, 1997Go); high {delta}18O values record melting of crust with variable proportions of supracrustal-derived rock.

Sequoia region rocks have distinctly higher Sri values relative to the eastern and western series, but the data are more scattered (Fig. 10a). The variable {delta}18O and Sri of this sample group exemplifies the widespread but heterogeneous contamination by the Kings Sequence.

Oxygen and lead isotopes
Lead isotope ratios positively correlate with {delta}18O(Zrc) and western and eastern trends are distinct (Figs 10b–d). Calculated whole-rock {delta}18O again shows more scatter, but, nevertheless, consistent patterns (Figs 10f–h). All initial Pb isotope ratios (206Pb/204Pb; 207Pb/204Pb; 208Pb/204Pb) become more radiogenic with increasing {delta}18O(Zrc). The lowest {delta}18O values of both eastern and western samples tend to have similar initial 206Pb/204Pb and 207Pb/204Pb, close to the primitive mantle value at 100 Ma (Chen & Tilton, 1991Go). Chen & Tilton (1991Go) suggested that eastern samples have higher 208Pb/204Pb because U/Pb was lower than Th/Pb in the source rocks, an indication of preferential extraction of U by fluid escape during granulite-facies metamorphism.

The increasing variability of initial Pb isotope composition with {delta}18O in western samples produces a ‘fan-like’ spread of the data, which suggests a common low-{delta}18O mantle reservoir input (Fig. 10b–d), but highlights source or contamination complexity in high-{delta}18O rocks that is not evident from Sri (Fig. 10a). Most western samples document input of aged crust or sediments thereof; however, a group of four samples from the westernmost rocks change little in initial Pb, whereas {delta}18O(Zrc) ranges from 6·0 to 7· 0{per thousand}; these same samples trend toward increasing Sri (Fig. 10a). Not surprisingly, the O–Pb–Sr systematics of these samples are consistent with incorporation of altered ocean crust or arc volcanic rocks. Elevated Sri values but low initial Pb values show decoupling of these radiogenic isotope systems during supracrustal alteration of the source rocks.

Sequoia samples are also scattered with respect to 207Pb/204Pb and 208Pb/204Pb (Fig. 10c and d); 207Pb/204Pb values overlap considerably with those for western and eastern rocks (Fig. 10c). The highly variable {delta}18O–Pb pattern, and the fact that many samples fall between the western and eastern trends, suggests a diverse source and contamination history for this group of plutons.

Large-scale magma systems
Positive correlation of O, Pb, and Sr isotopes indicates both mantle and crustal inputs in the western and eastern Sierra magmas. Each sample defining these trends is a complex integration of source characteristics and contamination, magma mixing, differentiation, and assimilation, obscuring the origin. Nevertheless, correlations for eastern SNB samples are reasonably strong in most cases and distinct in slope from western SNB rocks. This is a remarkable pattern considering that samples in each group vary considerably in age, up to several million years, and location within their respective domains (Fig. 3). Therefore, we hypothesize that two distinct magma systems operated.

It is widely acknowledged that different ‘styles’ of magmatism occur in the western and eastern Sierra (e.g. Chen & Tilton, 1991Go; Coleman & Glazner, 1997Go; Tikoff & de Saint Blanquat, 1997Go; Wenner & Coleman, 2004Go), and the Peninsular Ranges batholith bears some similarity (Gromet & Silver, 1987Go; Tulloch & Kimbrough, 2003Go; Lee et al., 2007Go). The isotopic patterns imply unique controls on isotope ratios on each side of the batholith that permitted largely continuous trans-batholith gradients in major and minor element chemistry (e.g. Bateman, 1992Go).

Differing lithosphere compositions in the SNB undoubtedly provide a major control on the isotopic composition of the magmas (Kistler, 1990Go); yet to average out heterogeneity of crust and mantle reservoirs and produce well-defined isotopic trends requires a long-lived and stable source region, such as is envisaged by the melting, assimilation, storage, and homogenization (MASH) model of Hildreth & Moorbath (1988Go). Ducea (2001Go, 2002Go) and Saleeby et al. (2003Go) adopted the MASH model to explain magmatic averaging that results in complementary geochemistry between batholithic crust and xenoliths thought to represent residues of melting at the base of the arc crust. Annen et al. (2006Go) proposed the existence of deep crustal ‘hot zones’ beneath large magma systems in which the geochemical character of the magmas is largely imparted in the source but textural diversity is acquired in the crust during crystallization. Clearly, other models may be applicable, but the MASH and ‘hot zone’ models are apt in that they permit large-scale magma systems that could support long-term geochemical connectivity, but allow for progressive emplacement of texturally diverse plutons.

In addition to basement composition effects, the two magma systems probably varied in depth of melting. In particular, eastern Sierra {delta}18O values indicate mostly mantle sources. This finding concurs with a growing consensus of studies invoking direct melting of the lithospheric mantle and remelting of underplated juvenile magmas as the origin of most magmas in the eastern Sierra (Coleman & Glazner, 1997Go; Wenner & Coleman, 2004Go; Ratajeski et al., 2005Go; Sisson et al., 2005Go). Broadly characterized, this magma system was driven by massive infracrustal recycling at a deep locus of melting. Any crustal melting that did happen would involve lower crust composed of orthogneiss or amphibolite metamorphosed to granulite facies (Chen & Tilton, 1991Go); increases of {delta}18O from melting of such rocks would be minimal (Kistler, 1990Go).

Except for the Foothills belt, which represents the nascent stages of magmatism in the Sierran arc, greater crustal melting, recycling, and contamination, including abundant melting of supracrustal rocks, defines the western magma system. The high-{delta}18O magmas of the Fine Gold–Sequoia belt, and evidence of wall-rock controls on {delta}18O are consistent with a shallow locus of magma generation that melted and reworked the crust considerably. The high content of greywacke and hydrous, disaggregated ophiolite lithologies present in the basement rocks of the western Sierra potentially proved more fertile for melting and facilitated recycling (Montel & Vielzeuf, 1997Go; Vielzeuf & Schmidt, 2001Go).

Small-scale magma systems: plutons and intrusive suites
Beyond evidence for two large-scale magmatic domains in the SNB, we observe intriguing spatial and temporal patterns of {delta}18O variation within intrusive suites and individual plutons. Examples of these smaller-scale systems are discussed below and we evaluate the causes of their isotopic diversity and implications for magma sources, contamination, and batholith construction.

Mount Givens Granodiorite
The Mount Givens Granodiorite exhibits a gradual 1{per thousand} increase in {delta}18O(Zrc) on its west side (Fig. 11a) where the pluton abuts the older Dinkey Creek Granodiorite (102 Ma, Tobisch et al., 1995Go) at the PA/NA boundary (Fig. 6a). Compared with zircon, the calculated {delta}18O(WR) blurs this pattern because inward changes in lithology from a tonalitic margin to a granitic interior result in highly variable calculated {delta}18O(WR) values (Fig. 11b). Because textures and mineral proportions vary considerably (Bateman, 1992Go; McNulty et al., 2000Go), as does the range of U–Pb zircon ages (87·9–92·8 Ma), the pluton appears to be constructed by multistage growth (McNulty et al., 2000Go). The trend of {delta}18O(Zrc) suggests that the western margin of the pluton has a geochemical affinity with the PA domain, whereas the {delta}18O of the younger eastern areas reflects the NA domain. Similarly, the Eagle Peak pluton, just to the west of the Mount Givens pluton (Fig. 4), shows zoning in O, Sr, and Nd isotopes that hints at variable contamination of the magma as it ascended along the PA/NA break (Hill et al., 1988Go).


Figure 11
View larger version (14K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 11. Variation of {delta}18O(Zrc) in the Mount Givens pluton. (a) Increasing {delta}18O toward the western margin of the pluton approaches values of the Dinkey Creek Granodiorite across the PA/NA break. (b) Calculated {delta}18O(WR) values poorly resolve {delta}18O zoning because SiO2 content varies considerably in the pluton; curve fits are power law.

 
Sequoia Intrusive Suite
The Sequoia Intrusive Suite (Fig. 3) exhibits distinct {delta}18O zonation. The outer and largest pluton, the Giant Forest Granodiorite, has the highest {delta}18O(Zrc) values of the suite and greater overall variation in {delta}18O; the interior granites of Big Meadows and Weaver Lake have lower average {delta}18O values that are less variable (Fig. 12a). Average Sri also decreases inward in the suite: Giant Forest = 0·70799 ± 0·00204, n = 6; Big Meadows = 0·70733 ± 0·00011, n = 2; Weaver Lake = 0·70657 ± 0·00032, n = 4 (Chen & Tilton, 1991Go; Wenner & Coleman, 2004Go). Thus, the margin of the suite exhibits greater supracrustal contamination. In detail, samples from the Giant Forest Granodiorite have higher {delta}18O in the southern half of the pluton where it abuts the Kings Sequence rocks of the Sequoia pendant (Fig. 3, Table 1). Such localized effects indicate upper crustal contamination and explain the {delta}18O heterogeneity in the pluton. The Giant Forest pluton appears to have shielded the ‘core’ of the suite from contamination. Therefore, the lower isotope ratios of the interior plutons better represent the suite's source composition.


Figure 12
View larger version (23K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 12. {delta}18O(Zrc) of intrusive suites. (a) Sequoia Suite rocks show less variability in younger, interior members. (b) The Tuolumne Suite shows variation in {delta}18O among members, but the average {delta}18O is relatively constant. (c) {delta}18O increases inward in Sonora Suite rocks. (d) John Muir Suite plutons show complex patterns of {delta}18O with decreasing age. (e) Variation of {delta}18O vs SiO2 in Muir Suite rocks shows distinct northern and southern trends. Kin, Inconsolable Granodiorite; Klk, Lamark Granodiorite; Kle, Lake Edison Granodiorite; Krv, RoundValley Peak Granodiorite; Kmo,Mono Creek Granite.

 
Sierra Crest Suites
Evolution of magmatic {delta}18O(Zrc) from oldest to youngest units is apparent in some plutonic suites exposed in the Sierra Crest (Fig. 12b–d). Relatively constant {delta}18O(Zrc) values occur in the Tuolumne Intrusive Suite (Fig. 12b), but others show progressive {delta}18O changes (e.g. Sonora, John Muir, Fig. 12c and d). In addition, {delta}18O is often homogeneous within each pluton of a suite, but varies between plutons.

The John Muir Suite
Plutons of the John Muir Suite (JMS, Fig. 3) display intriguing spatial and temporal trends of {delta}18O(Zrc). A subset of the JMS, the Mono Pass Suite, is approximately coeval with the Tuolumne Intrusive Suite, and includes, from oldest (95·3 Ma) to youngest (87· 0 Ma), the gabbro of Rock Creek, Lake Edison Granodiorite, Round Valley Peak Granodiorite, and Mono Creek Granite (Gaschnig et al., 2006Go). The Lamarck Granodiorite (94–92 Ma, J. Gracely, personal communication) and the Evolution Basin Alaskite (94–92 Ma, Wenner & Coleman, 2004Go) are exposed to the south of the Mono Pass Suite, and a new U–Pb zircon age for the Inconsolable Granodiorite (95 Ma, J. Gracely, personal communication) suggests that it may be a low-silica end-member of the southern JMS.

Values of {delta}18O(Zrc) are reported for all the plutons of the JMS except the Rock Creek gabbro and Evolution Basin Alaskite (Table 1). Relative to pluton age, {delta}18O(Zrc) decreases with age in the two oldest plutons. It displays considerable variation in the Lake Edison pluton; {delta}18O(Zrc) is highest in the Round Valley Peak Granodiorite, and decreases in the Mono Creek Granite (Fig. 12d). Except for the Lake Edison Granodiorite, most plutons are relatively homogeneous in {delta}18O. Variation of {delta}18O vs SiO2 for John Muir samples also reveals two parallel trends of decreasing {delta}18O with increasing SiO2 (Fig. 12e). One trend, the ‘southern trend’ is defined by the two southern plutons (Inconsolable and Lamarck), whereas the Round Valley Peak, Mono Creek, and some of the Lake Edison samples define the ‘northern trend’.

To explain the offset {delta}18O trends in the JMS, we propose that the southern and northern trend plutons were produced from two distinct sources. The Lake Edison Granodiorite, with its wide range of {delta}18O, suggests a hybrid pluton produced by mixing between the two sources. Such behavior supports models for magma mixing as the origin of some granodiorites in the Sierra Nevada (Sisson et al., 1996Go; Wenner & Coleman, 2004Go). In addition, higher {delta}18O in the older plutons indicates that the early formed magmas contained higher proportions of recycled supracrustal rocks than later magmas. Decreasing {delta}18O in the younger plutons is consistent with exhaustion of fusible crustal rocks in the magma sources, increased mantle contributions, or possibly conditioning of magma conduits (Scoates & Frost, 1996Go).

Temporal patterns of {delta}18O
Considering changes in {delta}18O(Zrc) over time for all rocks gives a more complete sense of the evolving arc magma system (Fig. 13a). In addition, magma production rates deduced from measurement of pluton area (Fig. 13b) reveal that the Cretaceous magma system experienced several episodic magma ‘pulses’. These pulses are resolved at finer time-scales than the magmatic ‘flare-ups’ of Ducea (2001Go), and they correspond to changes in {delta}18O.


Figure 13
View larger version (35K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 13. Temporal variation of {delta}18O(Zrc) (a) and pluton area (b) in the central Sierra from 125 to 80 Ma. Gradients of {delta}18O in (a) are discontinuous between different pulses of magmatism. The Bass Lake Tonalite can be subdivided according to the ages of lobes within the body (105–116 Ma; Stern et al., 1981Go), or a single age of 114 Ma may be assigned (Bateman, 1992Go). The Sonora Intrusive Suite is not included in area calculations because incomplete mapping and geochronology west of the suite prevents estimation of the full intrusive area across that part of the batholith. Image analysis was employed to calculate the area of individual plutons that have been dated, including some that were not analyzed for {delta}18O. (See Electronic Appendix 1 for age references.)

 
The first magmatic pulse, marked by increasing {delta}18O, resembles the Peninsular Ranges trend, and shows more supracrustal input with time. This period includes the transition from magmatism in the Foothills belt to the Fine Gold Suite and western Sequoia region. Subsequent pulses at 103–102 Ma and 99–97 Ma in the western and axial potions of the SNB saw widening of the range of {delta}18O to both higher and lower values. A temporary decrease in magma production at around 95 Ma (Fig. 13b) led to voluminous magmatism in the eastern SNB. The pronounced downward shift of {delta}18O(Zrc) emphasizes the difference of the eastern Sierra magma system.

In arc settings, magmatism, heat flow, and development of MASH domains or ‘hot zones’ should facilitate greater crustal recycling with time. A secular change to increased recycling of supracrustal rocks is demonstrated globally in a study of {delta}18O of igneous zircon spanning most of Earth's magmatic history (Valley et al., 2005Go), and increasing ‘supracrustal oxygen’ with time is seen during magmatic refinement of the Gondwana supercontinent (Kemp et al., 2006Go) and Superior craton (Moser et al., 2008Go). In the SNB, individual magmatic pulses generally record greater crustal recycling with time (Fig. 13), including the eastern Sierra where {delta}18O(Zrc) increases slightly. Some short-lived pulses (e.g. 103–102 Ma) are isotopically homogeneous, suggesting insufficient time for the onset of melting of crustal material. Overall, the SNB demonstrates that arcs may not record the same supracrustal oxygen increase as occurs during long-term crustal growth, but that such patterns are resolvable at the time-scale of individual magmatic episodes.

Evolution of the Sierran arc
We now integrate spatial and temporal patterns of SNB magma geochemistry and present a tectonic model in the form of a scheme portraying three critical times during construction of the batholith (Fig. 14).


Figure 14
View larger version (48K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 14. Model for the origin of the distinct {delta}18O belts of the central Sierra Nevada. (a) Early Cretaceous setting before voluminous magmatism, showing distinct basement and lithospheric mantle domains, overlying sediments, and older plutons. The Foothills and eastern Sierra crustal breaks are shown as trans-crustal; ‘?’ symbols indicate uncertainty of the geometry and depth of the PA/NA break. The transition from oceanic to North American lithospheric mantle is also uncertain, although an east-dipping contact is inferred. (b) By 110 Ma, increased magmatism mobilized accreted supracrustal rocks into the western SNB high-{delta}18O belt, inboard of the low-{delta}18O belt. A relatively shallow western MASH zone resulted in widespread crustal melting. (c) At 90 Ma, transpression and increased rates of subduction facilitated emplacement of voluminous eastern SNB magmas from lithospheric mantle sources. Intra-arc shortening by lithospheric under-thrusting temporarily stifled magmatism, possibly inserting Proterozoic lower crust into the sub-arc mantle. ‘Eclogitic root’ refers to eclogitic residues expelled from the MASH zone. Transect line corresponds to Fig. 5.

 
Following the Nevadan Orogeny, the future Sierra was underlain by basement oceanic terranes accreted to North America and overlying sedimentary cover sequences (Fig. 14a). Subduction and underplating, or underthrusting in the Triassic and Jurassic probably contaminated the sub-arc lithospheric mantle. The eastern Sierra inherited Proterozoic North American crust and complementary lithospheric mantle, as shown by the {delta}18O break recorded in Jurassic and Triassic plutons that stitch across the eastern Sierra Break (Fig. 14a), and the patterns of supracrustal contamination in these and Cretaceous plutons.

Subduction of the Farallon Plate initiated major magmatism by 120 Ma and produced significant batholith volume by 110 Ma (Figs 13a and 14b). Largely depleted mantle sources produced the plutons in the Foothills Zone, with minimal crustal interaction. The Fine Gold Suite, however, documents the onset of massive supracrustal recycling. Increased involvement of the wall-rocks of the Calaveras Complex, Kings Sequence, and cryptic contributions from altered ocean crust and Proterozoic rocks led to more complex isotopic signatures in magmas towards the PA/NA break.

The later stages of magmatism in the SNB highlight the role of tectonic reconfigurations in magma evolution. The deformation field of the Sierran arc is believed to have switched from contraction to transpression and dextral shearing because of changing convergence of the Farallon Plate in the late Cretaceous (Tobisch & Cruden, 1995Go; Tobisch et al., 1995Go; Tikoff & de Saint Blanquat, 1997Go; Fig. 14c). Recent work on the proto-Kern Canyon fault shows evidence for intra-arc shortening at 95 Ma in response to this convergence switch (Nadin & Saleeby, 2008Go). Shortening potentially suppressed voluminous magmatism, as is seen in the lull at 95 Ma (Fig. 13b). However, post 95 Ma transpression established a positive feedback between magmatism and deformation, facilitating emplacement of the Sierra Crest magmas (Glazner, 1991Go; Tikoff & de Saint Blanquat, 1997Go). Shear-enhanced magma emplacement of eastern Sierra magmas possibly facilitated intrusion from MASH zone sources into the upper crust, thereby reducing crustal contamination. Likewise, transpressional conduits may have affected the constancy of {delta}18O seen within some units. Homogeneous {delta}18O in some plutons in the Tuolumne Suite emplaced over 2–3 Myr (Coleman et al., 2004Go) favors enduring magma sources and conduits. Spatially variable {delta}18O in plutons such as the Mount Givens indicates heterogeneous magmas, possibly the resulting of switching of sources (and conduits) during deformation, varying amounts of contamination, or magma mixing in evolving magma plumbing systems. In the southern Sierra, greater recycling of crustal rocks with time (Lackey et al., 2005Go) corresponds to increased ductile deformation, including both vertical and lateral transport of wall-rocks in response to elevated geothermal gradients (Saleeby et al., 2003Go). Thus, tectonic ‘re-shuffling’ in the arc, such as downward flow of wall-rocks, or intra-arc thrusting (Fig. 14c), could possibly have rejuvenated magma sources with crustal rocks.


    CONCLUDING REMARKS
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL BACKGROUND
 SAMPLE SELECTION, PREPARATION,...
 RESULTS
 DISCUSSION
 CONCLUDING REMARKS
 SUPPLEMENTARY DATA
 REFERENCES
 
The data presented in this study and by Lackey et al. (2005Go, 2006Go) put in place a modern and comprehensive oxygen isotope evaluation of the Sierra Nevada batholith. Values of {delta}18O(Zrc) provide clear records of changing magmatic {delta}18O in the batholith, and calculation of {delta}18O(WR) values using zircon as a magmatic benchmark can be used to circumvent widespread resetting of the {delta}18O of whole-rocks and single minerals such as quartz. Moreover, examples from the Barcroft pluton and Saddlebag Lake meta-rhyolite illustrate how {delta}18O analysis of a refractory mineral such as zircon is a powerful means to deconvolve magmatic and metamorphic histories previously unobtainable from pervasively altered rocks.

Recognition of distinct belts of {delta}18O values reveals both transverse and longitudinal domains within the batholith and a markedly more complex system than in the coeval Peninsular Ranges batholith. Large-scale pre-batholith boundaries resolved by {delta}18O mapping better define the locations of previously inferred terrane boundaries, including cryptic boundaries such as the Foothills Break. Furthermore, the pre-batholith arrangement of basement and supracrustal rocks, structural reconfiguration during the evolution of the batholith, and different magma source depths, led to the overall patterns of {delta}18O. Lastly, the discrete magma systems that we resolve at different scales in time and space, from batholith to pluton scale, emphasizes the dynamic nature of the Sierran arc.


    SUPPLEMENTARY DATA
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL BACKGROUND
 SAMPLE SELECTION, PREPARATION,...
 RESULTS
 DISCUSSION
 CONCLUDING REMARKS
 SUPPLEMENTARY DATA
 REFERENCES
 
Supplementary data for this paper are available at Journal of Petrology online.


    ACKNOWLEDGEMENTS
 
All isotope analyses were completed by J.S.L. in the University of Wisconsin Stable Isotope Laboratory with support from DOE 93ER14389 and NSF EAR99-02973 and EAR02-07340. A UW Dean Morgridge Distinguished Graduate Fellowship, GSA and Sigma Xi grants, and the UW Department of Geology and Geophysics Weeks Fund supported portions of this work. Pomona College Geology Department faculty research funds provided valuable support for publication of this study. We thank Mike Spicuzza for assistance with stable isotope analysis, Brian Hess for making thin sections and grain mounts, and Clark Johnson and Brian Beard for the use of the rock crushing and mineral separation laboratories. Clark Johnson, Brad Singer, Elizabeth King, William Peck, Ilya Bindeman, Cory Clechenko, Aaron Cavosie, and Tom Lapen have added valuable discussion of portions of this research. We thank Martha House for many of the samples. Jason Saleeby, Kent Ratajeski, and Basil Tikoff are thanked for their helpful discussions and samples. Drew Coleman, Rich Gaschnig, and John Gracely, at UNC, kindly provided unpublished U–Pb ages and helpful discussions of the geochronology of the eastern Sierra. Ron Kistler graciously shared unpublished isotope data. Hans Hinke assisted in the field. David Graber and Jan van Wagtendonk helped us obtain permits to sample in Sequoia, Kings Canyon, and Yosemite National Parks. Detailed and encouraging reviews by Allen Glazner, Tom Sisson, Ron Frost, and Charlotte Allen helped improve this paper.


*Corresponding author. Telephone: (909) 621-8677. Fax: (909) 621-8552. E-mail: JadeStar.Lackey{at}pomona.edu


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL BACKGROUND
 SAMPLE SELECTION, PREPARATION,...
 RESULTS
 DISCUSSION
 CONCLUDING REMARKS
 SUPPLEMENTARY DATA
 REFERENCES
 
Ague J. J. Thermodynamic calculation of emplacement pressures for batholithic rocks, California: Implications for the aluminum-in-hornblende barometer. Geology (1997) 25:563–566.[Abstract/Free Full Text]

Ague J. J., Brimhall G. H. Magmatic arc asymmetry and distribution of anomalous plutonic belts in the batholiths of California: effects of assimilation, crustal thickness, and depth of crystallization. Geological Society of America Bulletin (1988a) 100:912–927.[Abstract/Free Full Text]

Ague J. J., Brimhall G. H. Regional variations in bulk chemistry, mineralogy, and the compositions of mafic and accessory minerals in the batholiths of California: with Suppl. Data 88–13. Geological Society of America Bulletin (1988b) 100:891–911.[Abstract/Free Full Text]

Annen C, Blundy JD, Sparks RSJ. The genesis of intermediate and silicic magmas in deep crustal hot zones. Journal of Petrology (2006) 47:505–539.[Abstract/Free Full Text]

Appora I., Eiler J. M., Matthews A., Stolper E. M. Experimental determination of oxygen isotope fractionations between CO2 vapor and soda-melilite melt. Geochimica et Cosmochimica Acta (2003) 67:459–471.[CrossRef][Web of Science]

Bateman P. C. Plutonism in the central part of the Sierra Nevada Batholith, California. US Geological Survey, Professional Papers (1992) 1483:1–186.

Bateman P. C., Chappell B. W. Crystallization, fractionation, and solidification of the Tuolumne Intrusive Series, Yosemite National Park, California. Geological Society of America Bulletin (1979) 90:465–482.[Abstract/Free Full Text]

Bateman P. C., Dodge F. C. W., Kistler R. W. Magnetic susceptibility and relation to initial 87Sr/86Sr for granitoids of the central Sierra Nevada, California. Journal of Geophysical Research, B, Solid Earth and Planets (1991) 96:19555–19568.[CrossRef]

Bindeman I. N., Valley J. W. Oxygen isotope study of the Long Valley magma system, California: isotope thermometry and convection in large silicic magma bodies. Contributions to Mineralogy and Petrology (2002) 144:185–205.[Web of Science]

Böhlke J. K., Kistler R. W. Rb–Sr, K–Ar, and stable isotope evidence for the ages and sources of fluid components of gold-bearing quartz veins in the northern Sierra Nevada foothills metamorphic belt, California. Economic Geology (1986) 81:296–322.[Abstract/Free Full Text]

Chacko T., Cole D. R., Horita J. Equilibrium oxygen, hydrogen, and carbon isotope fractionation factors applicable to geological systems. In: Stable Isotope Geochemistry. Mineralogical Society of America/Geochemical Society, Reviews in Mineralogy and Geochemistry —Valley J. W., Cole D. R., eds. (2001) 43:1–81.

Chen J. H., Moore J. G. Late Jurassic Independence dike swarm in eastern California. Geology (1979) 7:129–133.[Abstract/Free Full Text]

Chen J. H., Moore J. G. Uranium–lead isotopic ages from the Sierra Nevada batholith, California. Journal of Geophysical Research B (1982) 87:4761–4784.[CrossRef]

Chen J. H., Tilton G. R. Applications of lead and strontium isotopic relationships to the petrogenesis of granitoid rocks, central Sierra Nevada batholith, California. Geological Society of America Bulletin (1991) 103:437–447.

Clayton R. N., Goldsmith J. R., Mayeda T. K. Oxygen isotope fractionation in quartz, albite, anorthite and calcite. Geochimica et Cosmochimica Acta (1989) 53:725–733.[CrossRef][Web of Science]

Clemens-Knott D. Geologic and isotopic investigations of the Early Cretaceous Sierra Nevada batholith, Tulare County, California, and the Ivrea Zone, Northwest Italian Alps: examples of interaction between mantle-derived magma and continental crust. In: Ph.D. thesis (1992) Pasadena: California Institute of Technology.

Coleman D. S., Glazner A. F. The Sierra Crest magmatic event: rapid formation of juvenile crust during the Late Cretaceous in California. International Geology Review (1997) 39:768–787.

Coleman D. S., Frost T. P., Glazner A. F. Evidence from the Lamarck Granodiorite for rapid Late Cretaceous crust formation in California. Science (1992) 258:1924–1926.[Abstract/Free Full Text]

Coleman D. S., Gray W., Glazner A. F. Rethinking the emplacement and evolution of zoned plutons: geochronologic evidence for incremental assembly of the Tuolumne Intrusive Suite, California. Geology (2004) 32:433–436.[Abstract/Free Full Text]

Couglan RAN. Studies in diffusional transport: Grain boundary transport of O in feldspars, diffusion of O, strontium, and the REEs in garnet and thermal histories of granitic intrusions in south–central Maine using O isotopes. Ph.D. thesis (1990) Brown University, Providence, RI.

DePaolo D. J. Sources of continental crust: neodymium isotope evidence from the Sierra Nevada and Peninsular Ranges. Science (1980) 209:684–687.[Abstract/Free Full Text]

DePaolo D. J. A neodymium and strontium isotopic study of the Mesozoic calc-alkaline granitic batholiths of the Sierra Nevada and Peninsular Ranges, California. Journal of Geophysical Research B (1981) 86:10470–10488.[CrossRef]

Dodson M. H. Closure temperature in cooling geochronological and petrologic systems. Contributions to Mineralogy and Petrology (1973) 40:259–274.[CrossRef][Web of Science]

Doe B. R., Delevaux M. H. Variations in lead-isotopic compositions in Mesozoic granitic rocks of California: a preliminary investigation. Geological Society of America Bulletin (1973) 84:3513–3526.[Abstract/Free Full Text]

Ducea MN. A petrologic investigation of deep-crustal and upper-mantle xenoliths from the Sierra Nevada, California: constraints on lithospheric composition beneath continental arcs and the origin of Cordilleran batholiths. In: Ph.D. thesis (1998) Pasadena: California Institute of Technology.

Ducea M. N. The California Arc: thick granitic batholiths, eclogitic residues, lithospheric-scale thrusting, and magmatic flare-ups. GSA Today (2001) 11:4–10.

Ducea M. N. Constraints on the bulk composition and root foundering rates of continental arcs: a California arc perspective. Journal of Geophysical Research, B, Solid Earth and Planets (2002) 107:ECV 15.1–15.13.

Ernst W. G., Rumble D. III. Oxygen isotopic study of late Mesozoic cooling of the Mount Barcroft area, central White Mountains, eastern California. Contributions to Mineralogy and Petrology (2003) 144:639–651.[Web of Science]

Ernst W. G., Coleman D. S., Van de Ven C. M. Petrochemistry of granitic rocks in the Mount Barcroft area: implications for arc evolution, central White Mountains, easternmost California. Geological Society of America Bulletin (2003) 115:499–512.[Abstract/Free Full Text]

Farver J. R., Yund R. A. Oxygen diffusion in quartz: dependence on temperature and water fugacity. Chemical Geology (1991) 90:55–70.[CrossRef][Web of Science]

Fliedner M. M., Klemperer S. L., Christensen N. I. Three-dimensional seismic model of the Sierra Nevada Arc, California, and its implications for crustal and upper mantle composition. Journal of Geophysical Research, B, Solid Earth and Planets (2000) 105:10899–10921.[CrossRef]

Fu B., Page F. Z., Cavosie A. J., Fournelle J. H., Kita N. T., Lackey J. S., Wilde S. A., Valley J. W. Ti-in-zircon thermometry: applications and limitation. Contributions to Mineralogy and Petrology (2008) doi:10.1007/s00410-008-0281-5.

Garlick G. D. Oxygen isotope fractionation in igneous rocks. Earth and Planetary Science Letters (1966) 1:361–368.[CrossRef]

Gaschnig R. M., Coleman D. S., Glazner A. F. Twin of the Tuolumne: new geochronology from the Mono Pass intrusive suite. Geological Society of America, Abstracts with Programs (2006) 38:559.

Glazner A. F. Plutonism, oblique subduction, and continental growth: an example from the Mesozoic of California. Geology (1991) 19:784–786.[Abstract/Free Full Text]

Godfrey J. D. The deuterium content of hydrous minerals from the east–central Sierra Nevada and Yosemite National Park. Geochimica et Cosmochimica Acta (1962) 26:1215–1245.[CrossRef][Web of Science]

Gromet P. L., Silver L. T. REE variations across the Peninsular Ranges batholith: implications for batholithic petrogenesis and crustal growth in magmatic arcs. Journal of Petrology (1987) 28:75–125.[Abstract/Free Full Text]

Hildreth W., Moorbath S. Crustal contributions to arc magmatism in the Andes of central Chile. Contributions to Mineralogy and Petrology (1988) 98:455–489.[CrossRef][Web of Science]

Hill M., O’Neil J. R., Noyes H., Frey F. A., Wones D. R. Sr, Nd, and O isotope variations in compsitionally zoned and unzoned plutons in the central Sierra Nevada batholith. American Journal of Science (1988) 288-A:213–214.

House M. A., Wernicke B. P., Farley K. A., Dumitru T. A. Cenozoic thermal evolution of the central Sierra Nevada, California, from (U–Th)/He thermochronometry. Earth and Planetary Science Letters (1997) 151:167–179.[CrossRef][Web of Science]

House M. A., Wernicke B. P., Farley K. A. Dating topography of the Sierra Nevada, California, using apatite (U–Th)/He ages. Nature (1998) 396:66–69.[CrossRef]

House M. A., Wernicke B. P., Farley K. A. Paleo-geomorphology of the Sierra Nevada, California, from (U–Th)/He ages in apatite. American Journal of Science (2001) 301:77–102.[Abstract/Free Full Text]

Jennings CW, Strand RG, Rogers TH. Geologic map of California (1977) Washington, DC: California Division of Mines and Geology. Williams and Heintz Map Corporation, 1 sheet, 1:750 000.

Kemp A. I. S., Hawkesworth C. J., Paterson B. A., Kinny P. D. Episodic growth of the Gondwana supercontinent from hafnium and oxygen isotopes in zircon. Nature (2006) 439:580–583.[CrossRef][Medline]

King E. M., Valley J. W. The source, magmatic contamination, and alteration of the Idaho batholith. Contributions to Mineralogy and Petrology (2001) 142:72–88.[Web of Science]

Kistler R. W. Two different lithosphere types in the Sierra Nevada, California. In: The Nature and Origin of Cordilleran Magmatism. Geological Society of America, Memoirs —Anderson L., ed. (1990) 174:271–281.

Kistler RW. Mesozoic Intrabatholithic faulting, Sierra Nevada, California. In. In: Mesozoic Paleography of the Western United States II—Dunne G, McDougall K, eds. (1993) Los Angeles: Pacific Section, SEPM. 247–262.

Kistler R. W., Fleck R. J. Field guide for a transect of the central Sierra Nevada, California: geochronology and isotope geology. US Geological Survey, Open-File Report (1994) 94-0267:50.

Kistler R. W., Peterman Z. E. Variations in Sr, Rb, K, Na, and initial Sr87/Sr86 in Mesozoic granitic rocks and intruded wall rocks in central California. Geological Society of America Bulletin (1973) 84:3489–3512.[Abstract/Free Full Text]

Kistler R. W., Peterman Z. E. Reconstruction of crustal blocks of California on the basis of initial strontium isotopic compositions of Mesozoic granitic rocks. US Geological Survey, Professional Papers (1978) 1071:17.

Kistler R. W., Ross D. C. A strontium isotopic study of plutons and associated rocks of the southern Sierra Nevada and vicinity, California. US Geological Survey Bulletin (1990) 1920:20.

Kistler R. W., Chappell B. W., Peck D. L., Bateman P. C. Isotopic variation in the Tuolumne intrusive suite, central Sierra Nevada, California. Contributions to Mineralogy and Petrology (1986) 94:205–220.[CrossRef][Web of Science]

Kohn M. J., Valley J. W. Effects of cation substitutions in garnet and pyroxene on equilibrium oxygen isotope fractionations. Journal of Metamorphic Geology (1998) 16:625–639.[CrossRef][Web of Science]

Kylander-Clark A. R. C., Coleman D. S., Glazner A. F., Bartley J. M. Evidence for 65 km of dextral slip across Owens Valley, California, since 83 Ma. Geological Society of America Bulletin (2005) 117:962–968.[Abstract/Free Full Text]

Lackey JS. The magmatic and alteration history of the Sierra Nevada batholith as recorded by oxygen isotope ratios in zircon, titanite, garnet, and quartz. In: PhD thesis (2005) Madison: University of Wisconsin.

Lackey J. S., Valley J. W., Chen J. H. Correlated O–Sr–Pb isotope ratios in the West–Central Sierra Nevada batholith, California. Geological Society of America, Abstracts with Programs (2001) 33:295.

Lackey J. S., Valley J. W., Stockli D. F., House M. A. Magmatic processes in the central Sierra Nevada batholith: a cryptic pre-Jurassic boundary in Long Valley, California. Geological Society of America, Abstracts with Programs (2003) 35:92.

Lackey J. S., Valley J. W., Saleeby J. B. Supracrustal input to magmas in the deep crust of Sierra Nevada batholith: evidence from high-{delta}18O zircon. Earth and Planetary Science Letters (2005) 235:315–330.[CrossRef][Web of Science]

Lackey J. S., Valley J. W., Hinke H. J. Deciphering the source and contamination history of peraluminous magmas using {delta}18O of accessory minerals: examples from garnet-bearing granitoids of the Sierra Nevada batholith. Contributions to Mineralogy and Petrology (2006) 151:20–44.[CrossRef][Web of Science]

Lee C. T., Yin Q., Rudnick R. L., Chesley J. T., Jacobsen S. B. Osmium isotopic evidence for Mesozoic removal of lithospheric mantle beneath the Sierra Nevada, California. Science (2000) 289:1912–1916.[Abstract/Free Full Text]

Lee C.-T., Morton D. M., Kistler R. W., Baird A. K. Petrology and tectonics of Phanerozoic continent formation: from island arcs to accretion and continental arc magmatism. Earth and Planetary Science Letters (2007) 263:370–387.[CrossRef][Web of Science]

Magaritz M., Taylor H. P. Jrs. Oxygen, hydrogen and carbon isotope studies of the Franciscan Formation, Coast Ranges, California. Geochimica et Cosmochimica Acta (1976) 40:215–234.[CrossRef][Web of Science]

Masi U., O’Neil J. R., Kistler R. W. Stable isotope systematics in Mesozoic granites of central and northern California and southwestern Oregon. Contributions to Mineralogy and Petrology (1981) 76:116–126.[CrossRef][Web of Science]

McCulloch M. T., Gregory R. T., Wasserburg G. J., Taylor H. P. A neodymium, strontium, and oxygen isotope study of the Cretaceous Samail Ophiolite and implications for the petrogenesis and seawater–hydrothermal alteration of oceanic crust. Earth and Planetary Science Letters (1980) 46:201–211.[CrossRef][Web of Science]

McNulty B. A., Tobisch O. T., Cruden A. R., Gilder S. Multi-stage emplacement of the Mt. Givens pluton, central Sierra Nevada, California, Geological Society of America Bulletin (2000) 112:103–119.

Montel J. M., Vielzeuf D. Partial melting of metagreywackes. 2. Compositions of minerals and melts. Contributions to Mineralogy and Petrology (1997) 128:176–196.[CrossRef][Web of Science]

Moore J. G. The quartz diorite boundary line in the western United States. Journal of Geology (1959) 67:198–210.[Web of Science]

Moore J. G. Exploring the Highest Sierra (2000) Palo Alto, CA: Stanford University Press. 427.

Moore, J. G. & Sisson TW. Preliminary geologic map of Sequoia and Kings Canyon national parks, California. US Geological Survey, Open-File Report (1987) 87-0651:1. Sheet, 1:250 000.

Moser D. E., Bowman J. R., Wooden J. L., Valley J. W., Mazdab F., Kita N. Creation of a continent recorded in zircon zoning. Geology (2008) 36:239–242.[Abstract/Free Full Text]

Muehlenbachs K. The oxygen isotopic composition of the oceans, sediments and the seafloor. Chemical Geology (1998) 145:263–273.[CrossRef][Web of Science]

Mukhopadhyay B., Manton W. I. Upper-mantle fragments from beneath the Sierra Nevada Batholith: partial fusion, fractional crystallization, and metasomatism in a subduction-related ancient lithosphere. Journal of Petrology (1994) 35:1417–1450.[Abstract/Free Full Text]

Nadin E. S., Saleeby J. B. Disruption of regional primary structure of the Sierra Nevada batholith by the Kern Canyon fault system, California. Geological Society of America, Special Papers (2008) 438:429–454.

Page F. Z., Ushikubo T., Kita N. T., Riciputi L. R., Valley J. W. High-precision oxygen isotope analysis of picogram samples reveals 2 µm gradients and slow diffusion in zircon. American Mineralogist (2007) 92:1772–1775.[Abstract/Free Full Text]

Peck W. H., Valley J. W., Corriveau L., Davidson A., McLelland J., Farber D. A. Oxygen-isotope constraints on terrane boundaries and origin of 1· 18–1· 13 Ga granitoids in the southern Grenville Province. Geological Society of America, Memoirs (2004) 197:163–183.

Ratajeski K., Glazner A. F., Miller B. V. Geology and geochemistry of mafic to felsic plutonic rocks in the Cretaceous intrusive suite of Yosemite Valley, California. Geological Society of America Bulletin (2001) 113:1486–1502.[Abstract/Free Full Text]

Ratajeski K., Glazner A. F., Sisson T. W. Experimental and geochemical evidence for derivation of the El Capitan Granite, California, by partial melting of hydrous gabbroic lower crust. Contributions to Mineralogy and Petrology (2005) 149:713–734.[CrossRef][Web of Science]

Robinson, A. C. & Kistler RW. Maps showing isotopic dating in the Walker Lake 1 degrees by 2 degrees Quadrangle, California and Nevada. US Geological Survey, Miscellaneous Field Studies Map MF-1382-N (1986) 49. 3 sheets, 1:250 000.

Ross DC. Generalized geologic map of the southern Sierra Nevada, California, showing the location of basement samples for which whole rock 18O has been determined. US Geological Survey, Open-File Report (1983a) 83-0904, 1 sheet, 1:250 000.

Ross DC. Generalized geologic map of the southern Sierra Nevada, California, showing the location of samples for which K–Ar radiometric age data and Rb/Sr data have been determined. US Geological Survey, Open-File Report (1983b) 83-0231, 1 sheet, 1:250 000.

Ross D. C. Generalized geologic map of the basement rocks of the southern Sierra Nevada, California. US Geological Survey, Open-File Report (1987) 87-0276:30.

Saleeby J. B. Petrotectonic and paleogeographic settings of U.S. Cordilleran ophiolites. In: The Cordilleran Orogen: Conterminous US: The Geology of North America—Burchfiel B. C., Lipman P. W., Zoback M. L., eds. (1992) Boulder, CA: Geological Society of America. 653–682.

Saleeby JB, Busby C. Paleogeographic and tectonic setting of axial and western metamorphic framework rocks of the southern Sierra Nevada, California. In:. In: Mesozoic Paleogeography of the Western United States: II—Dunne GC, McDougall K, eds. (1993) Los Angeles, CA: Pacific Section 71, SEPM. 197–225.

Saleeby JB, Speed RC, Blake MC, Allmendinger RW, Gans PB, Kistler RW, Ross DC, Stauber DA, Zoback ML, Griscom A, McCulloch DS. Centennial Continent/Ocean Transect 10, Corridor C-2: Central California offshore to the Colorado Plateau (1986) Boulder, CA: Geological Society of America, 63 pp. 2 sheets, 1:500 000.

Saleeby J. B., Sams D. B., Kistler R. W. U/Pb zircon, strontium, and oxygen isotopic and geochronological study of the southernmost Sierra Nevada Batholith, California. Journal of Geophysical Research, B, Solid Earth and Planets (1987) 92:10443–10466.[CrossRef]

Saleeby J., Ducea M., Clemens Knott D. Production and loss of high-density batholithic root, southern Sierra Nevada, California. Tectonics (2003) 22. doi:10.1029/2002TC001374.

Scoates J. S., Frost C. D. A strontium and neodymium isotopic investigation of the Laramie anorthosites, Wyoming, USA: Implications for magma chamber processes and the evolution of magma conduits in Proterozoic anorthosites. Geochimica et Cosmochimica Acta (1996) 60:95–107.[CrossRef][Web of Science]

Sharp Z. D., Giletti B. J., Yoder H. S. Jr. Oxygen diffusion rates in quartz exchanged with CO2. Earth and Planetary Science Letters (1991) 107:339–348.[CrossRef][Web of Science]

Silver L. T., Taylor H. P. Jr, Chappell B. Some petrological, geochemical and geochronological observations of the Peninsular Ranges Batholith near the international border of the U.S.A. and Mexico. In: Mesozoic Crystalline Rocks: Peninsular Ranges Batholith and Pegmatites, Point Sal Ophiolite—Abbott P. L., Todd V., eds. (1979) San Diego, CA: San Diego State University. 83–110.

Sisson T. W., Grove T. L., Coleman D. S. Hornblende gabbro sill complex at Onion Valley, California, and a mixing origin for the Sierra Nevada Batholith. Contributions to Mineralogy and Petrology (1996) 126:81–108.[CrossRef][Web of Science]

Sisson T. W., Ratajeski K., Hankins W. B., Glazner A. F. Voluminous granitic magmas from common basaltic sources. Contributions to Mineralogy and Petrology (2005) 148:635–661.[CrossRef][Web of Science]

Snow, C. A. & Scherer HH. Terranes of the western Sierra Nevada Foothills metamorphic belt, California: A critical review. International Geology Review (2006) 48:46–62.[Web of Science]

Spicuzza M. J., Valley J. W., McConnell V. S. Oxygen isotope analysis of whole rock via laser fluorination: An air-lock approach. Geological Society of America, Abstracts with Programs (1998) 30:80.

Stern T. W., Bateman P. C., Morgan B. A., Newell M. F., Peck D. L. Isotopic U–Pb ages of zircon from the granitoids of the central Sierra Nevada, California. US Geological Survey, Professional Papers (1981) 1185:17.

Stevens C. H., Stone P., Dunne G. C., Greene D. C., Walker J. D., Swanson B. J. Paleozoic and Mesozoic evolution of east–central California. International Geology Review (1997) 39:788–829.

Stockli D. F., Dumitru T. A., McWilliams M. O., Farley K. A. Cenozoic Tectonic Evolution of the White Mountains, California and Nevada. Geological Society of America Bulletin (2003) 115:788–816.[Abstract/Free Full Text]

Taylor H. P., Sheppard S. M. F. Igneous rocks: I. Processes of isotopic fractionation and isotope systematics. In: Stable Isotopes in High Temperature Geological Processes. Mineralogical Society of America, Reviews in Mineralogy —Valley J. W., Taylor H. P. Jr, O’Neil J. R., eds. (1986) 16:227–271.

Taylor H. P., Silver L. T. Oxygen isotope relationships in plutonic igneous rocks of the Peninsular Ranges Batholith, southern and Baja California. US Geological Survey, Open-File Report (1978) 78-0701:423–426.

Tikoff B., de Saint Blanquat M. Transpressional shearing and strike-slip partitioning in the Late Cretaceous Sierra Nevada magmatic arc, California. Tectonics (1997) 16:442–459.[CrossRef][Web of Science]

Tobisch O. T., Cruden A. R. Fracture-controlled magma conduits in an obliquely convergent continental magmatic arc. Geology (1995) 23:941–944.[Abstract/Free Full Text]

Tobisch O. T., Saleeby J. B., Renne P. R., McNulty B. A., Tong W. Variations in deformation fields during development of a large-volume magmatic arc, central Sierra Nevada, California. Geological Society of America Bulletin (1995) 107:148–166.[Abstract/Free Full Text]

Truschel JP. Petrogenesis of the Fine Gold intrusive suite, Sierra Nevada Batholith, California. In: M.S. thesis (1996) California: California State University, Northridge.

Tulloch A. J., Kimbrough D. L. Paired plutonic belts in convergent margins and the development of high Sr/Y magmatism: Peninsular Ranges Batholith of Baja California and Median Batholith of New Zealand. Geological Society of America, Special Papers (2003) 374:275–295.

Valley J. W. Oxygen isotopes in zircon. In: Zircon. Mineralogical Society of America/Geochemical Society, Reviews in Mineralogy and Geochemistry —Hanchar J. M., Hoskin P. W. O., eds. (2003) 53:343–385.

Valley J. W., Chiarenzelli J. R., McLelland J. M. Oxygen isotope geochemistry of zircon. Earth and Planetary Science Letters (1994) 126:187–206.[CrossRef][Web of Science]

Valley J. W., Kitchen N., Kohn M. J., Niendorf C. R., Spicuzza M. J. UWG-2, a garnet standard for oxygen isotope ratios: strategies for high precision and accuracy with laser heating. Geochimica et Cosmochimica Acta (1995) 59:5223–5231.[CrossRef][Web of Science]

Valley J. W., Kinny P. D., Schulze D. J., Spicuzza M. J. Zircon megacrysts from kimberlite: oxygen isotope variability among mantle melts. Contributions to Mineralogy and Petrology (1998) 133:1–11.[CrossRef][Web of Science]

Valley J. W., Bindeman I. N., Peck W. H. Empirical calibration of oxygen isotope fractionation in zircon. Geochimica et Cosmochimica Acta (2003) 67:3257–3266.[CrossRef][Web of Science]

Valley J. W., Lackey J. S., Cavosie A. J., Clechenko C., Spicuzza M. J., Basei M. A. S., Bindeman I. N., Ferreira V. P., Sial A. N., King E. M., Peck W. H., Sinha A. K., Wei C. S. 4·4 billion years of crustal maturation from oxygen isotope ratios of zircon. Contributions to Mineralogy and Petrology (2005) 150:561–580.[CrossRef][Web of Science]

Vielzeuf D., Schmidt M. W. Melting relations in hydrous systems revisited: application to metapelites, metagreywackes and metabasalts. Contributions to Mineralogy and Petrology (2001) 141:251–267.[Web of Science]

Vielzeuf D., Veschambre M., Brunet F. Oxygen isotope heterogeneities and diffusional profiles in composite metamorphic/magmatic garnets from the Pyrenees. American Mineralogist (2005) 90:463–472.[Abstract/Free Full Text]

Wenner J. M., Coleman D. S. Magma mixing and Cretaceous crustal growth: geology and geochemistry of granites in the central Sierra Nevada batholith, California. International Geology Review (2004) 46:880–903.[CrossRef][Web of Science]

Zeng L., Ducea M., Saleeby J. B. Geochemical characteristics of crustal anatexis during the formation of migmatite at the southern Sierra Nevada, California. Contributions to Mineralogy and Petrology (2005) 150:386–402.[CrossRef][Web of Science]

Zhao Z., Zheng Y. Calculation of oxygen isotope fractionation in magmatic rocks. Chemical Geology (2003) 193:59–80.[CrossRef][Web of Science]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
GeosphereHome page
C. G. Barnes, T. Prestvik, Y. Li, L. McCulloch, A. S. Yoshinobu, and C. D. Frost
Growth and zoning of the Hortavaer intrusive complex, a layered alkaline pluton in the Norwegian Caledonides
Geosphere, June 1, 2009; 5(3): 286 - 301.
[Abstract] [Full Text] [PDF]


Home page
GeologyHome page
J. W. Valley
The Origin of Habitats
Geology, November 1, 2008; 36(11): 911 - 912.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Data
Right arrow All Versions of this Article:
49/7/1397    most recent
egn030v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Lackey, J. S.
Right arrow Articles by Stockli, D. F.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?