Sensitive high-resolution ion microprobe U-Pb dating of prograde and retrograde ultrahigh-temperature metamorphism as exemplified by Sri Lankan granulites

Ultrahigh-temperature (UHT) granulites of the central Highland Complex, Sri Lanka, underwent some of the highest known peak temperatures of crustal metamorphism. Zircon and monazite U-Pb systems in granulites near Kandy, the highest grade region (similar to 1050 degrees C; 0.9 GPa), preserve both a record of the timing of prograde and retrograde phases of UHT metamorphism and evidence for the ages of older protolith components. Zircon grains from a quartz-saturated granulite containing relics of the peak UHT assemblage have remnant detrital cores with dates of ca. 2.5-0.83 Ga. Date clusters of ca. 1.7 and 1.04-0.83 Ga record episodes of zircon growth in the source region of the protolith sediment. Two generations of overgrowths with contrasting Th/U record metamorphic zircon growth at 569 +/- 5 and 551 +/- 7 Ma, probably in the absence and presence of monazite, respectively. The age of coexisting metamorphic monazite (547 +/- 7 Ma) is indistinguishable from that of the younger, low-Th/U zircon overgrowths. Zircon from a quartz-undersaturated monazite-absent UHT granulite with a mainly retrograde assemblage is mostly metamorphic (551 +/- 5 Ma). The ca. 570 Ma zircon overgrowths in the quartz-saturated granulite probably record partial melting just before or at the metamorphic peak. The ca. 550 Ma zircon in both rocks, and the ca. 550 Ma monazite in the quartz-saturated sample, record post-peak isothermal decompression. A possible model for this pressure-temperature-time evolution is ultrahot collisional orogeny during the assembly of Gondwana, locally superheated by basaltic underplating, followed by fast extensional exhumation.

Tab. 2: U-Pb ion-microprobe ages from Taymyr peninsula

We present the initial results of a U-Th-Pb zircon ion-microprobe investigation on samples from the Central Belt of Taimyr, in order to constrain its tectono-magmatic evolution. The zircon samples are from a deformed twomica granite (Faddey Massif), deformed metamorphosed gabbroic dike entrained as pods and lenses within metamorphosed tholeiitic basalts of the Kunar-Mod volcanic suite (Klyaz'ma River region), a metamorphosed rhyolite of the same volcanic suite overlying the basic metavolcanic rocks, as well as an undeformed dolerite dike which intrudes the metamorphosed Kunar-Mod basic volcanic rocks. Preliminary results on zircons from the two-mica granite suggest a crystallization age of ~630 Ma for this rock, with inheritance from assimilated crust 840 Ma to 1.1 Ga in age. In the Klyaz'ma River region, zircons from the meta-rhyolite yield a concordant age of -630 Ma. Zircons from the entrained metagabbroic dikes have so far yielded an age of -615 Ma (1 grain), as well as Archean ages (5 grains, concordant at 2.6-2.8 Ga). It seems likely that the Archean grains represent assimilation of older crustal material. Zircons from the post-tectonic dolerite dike have a bimodal age distribution. A well-defined younger age of 281 ±9 Ma is interpreted to represent the crystallization age of the dike, while older, concordant ages of 2.6-2.9 Ga likely represent assimilation of Archean crust (Siberian craton at depth). Several important conclusions can be drawn from the data. (1) The mafic and felsic lithologies of the Kunar-Mod volcanic suite are genetically related and should be the same age. Ages of-630 Ma (meta-rhyolite) and -615 Ma (metagabbroic dikes representing the latest stage of mafic magmatism associated the Kunar-Mod suite) suggest that these lithologies may be the same age, but more data are required to confirm this hypothesis. (2) The 630 Ma two-mica granite is similar in age to the time of high-grade metamorphism, suggesting that syntectonic granite emplacement accompanied obduction of the accretionary Central Belt to the Siberian craton. (3) An Early Permian age is well defined for the undeformed dolerite dike. Dolerite dikes occur across the whole of Taimyr, but are deformed to the south. If related, this single magmatic event pre-dates Permo-Triassic Siberian trap magmatism. Furthermore, it suggests that deformation was localized to southeastern Taimyr.

A rare earth element study of complex zircons from early Archaean Amitsoq gneisses, Godthåbsfjord, south-west Greenland

We present high spatial resolution ion-microprobe rare earth element (REE) data for discrete growth phases of complex polyphase zircons from early Archaean Amitsoq gneisses, outer Godthabsfjord, SW Greenland. In Matsuda diagrams, the two major growth phases, >3.8 Ga cores and ca. 3.65 Ga rims, have steep positive slopes from La to Lu, prominent positive Ce anomalies and negative Eu anomalies that are consistent with growth in a melt. Exceptions to this are non-cathodolurnmescent zircon developed between the cores and rims, sometimes truncating zoning in the cores, and late Archaean prismatic tip overgrowths, both of which exhibit flatter light REE (LREE) patterns and have small or no Eu anomaly, which we interpret as the result of metamorphism and/or small-degree, isolated partial melting. Our data support previous interpretations that the ca. 3.65 Ga zircon phase was generated in a melt, with the >3.8 Ga phase representing either original protolith zircons in a large degree partial melt or inherited zircons in an introduced magma. Regardless which of these two interpretations is correct for these, and similar, rocks in the outer GodthAbsfjord, the 3.65 Ga event will have profoundly affected isotopic systems and obscured beyond recognition any earlier igneous features such as cross-cutting relationships, which may only be assigned a minimum 3.65 Ga age. (C) 2003 Elsevier Science B.V. All rights reserved.

Inheritance of early Archaean Pb-isotope variability from long-lived Hadean protocrust

Comparison of initial Pb-isotope signatures of several early Archaean (3.65-3.82 Ga) lithologies (orthogneisses and metasediments) and minerals (feldspar and galena) documents the existence of substantial isotopic heterogeneity in the early Archaean, particularly in the Pb-207/Pb-204 ratio. The magnitude of isotopic variability at 3.82-3.65 Ga requires source separation between 4.3 and 4.1 Ga, depending on the extent of U/Pb fractionation possible in the early Earth. The isotopic heterogeneity could reflect the coexistence of enriched and depleted mantle domains or the separation of a terrestrial protocrust with a U-238/Pb-204 (mu) that was ca. 20-30% higher than coeval mantle. We prefer this latter explanation because the high-p signature is most evident in metasediments (that formed at the Earth's surface). This interpretation is strengthened by the fact that no straightforward mantle model can be constructed for these high-mu lithologies without violating bulk silicate Earth constraints. The Pb-isotope evidence for a long-lived protocrust complements similar Hf-isotope data from the Earth's oldest zircons, which also require an origin from an enriched (low Lu/Hf) environment. A model is developed in which greater than or equal to3.8-Ga tonalite and monzodiorite gneiss precursors (for one of which we provide zircon U-Pb data) are not mantle-derived but formed by remelting or differentiation of ancient (ca. 4.3 Ga) basaltic crust which had evolved with a higher U/Pb ratio than coeval mantle in the absence of the subduction process. With the initiation of terrestrial subduction at, we propose, ca. 3.75 Ga, most of the greater than or equal to3.8-Ga basaltic shell (and its differentiation products) was recycled into the mantle, because of the lack of a stabilising mantle lithosphere. We argue that the key event for preservation of all greater than or equal to3.8-Ga terrestrial crust was the intrusion of voluminous granitoids immediately after establishment of global subduction because of complementary creation of a lithospheric keel. Furthermore, we argue that preservation of !3.8-Ga material (in situ rocks and zircons) globally is restricted to cratons with a high U/Pb source character (North Atlantic, Slave, Zimbabwe, Yilgarn, and Wyoming), and that the Pb-isotope systematics of these provinces are ultimately explained by reworking of material that was derived from ca. 4.3 Ga (i.e. Hadean) basaltic crust.

Assigning dates to thin gneissic veins in high-grade metamorphic terranes: A cautionary tale from Akilia, southwest Greenland

A granodiorite from Akilia, southwest Greenland, previously suggested to date putative life-bearing rocks to greater than or equal to3.84 Ga, is re-investigated using whole-rock major and trace-element geochemistry, and detailed cathodoluminescence image-guided secondary ion mass spectrometer analyses of zircon U-Th-Pb and rare earth elements. Complex zircon internal structure reveals three episodes of zircon growth and/or recrystallization dated to c. 3.84 Ga, 3.62 Ga and 2.71 Ga. Rare earth element abundances imply a significant role for garnet in zircon generation at 3.62 Ga and 2.71 Ga. The 3.62 Ga event is interpreted as partial melting of a c. 3.84 Ga grey gneiss precursor at granulite facies with residual garnet. Migration of this 3.62 Ga magma (or melt-crystal mush) away from the melt source places a maximum age limit on any intrusive relationship. These early Archaean relationships have been complicated further by isotopic reworking in the 2.71 Ga event, which could have included a further episode of partial melting. This study highlights a general problem associated with dating thin gneissic veins in polyphase metamorphic terranes, where field relationships may be ambiguous and zircon inheritance can be expected.

湘西地区上震旦统上部－下寒武统底部黑色岩系沉积学及地球化学

Widespread black chert-shales occur in the Ediacaran-Cambrian(E-C) boundary successions along the flank of Yangtze Platform, South China, remarkable changes in sedimentology, geochemistry and biology were recorded. Although extensive studies were carried out upon this boundary succession, the origin of black chert-shales still remain controversial. This paper focuses on the E-C black chert-shales in western Hunan, South China, upon which detailed depositional and geochemical changes are documented, accordingly a depositional model for black chert-shales is proposed. Stratigraphic anatomy across the depositional strike demonstrates that the shallow-water Dengying dolostone along the platform margin sharply pass basinward into the Liuchapo chert successions, which indicate syndepositional extensional faulting at depth could have occurred along the platform margin. The deep-water Niutitang phosphorite-rich black shales are either underlain by the Dengying dolostones on the platform margin toward platform interior or directly by the Liuchaopo chert successions farther basinwards. By detailed investigation, silica chimneys are firsly identified approximately in the chert along platform margin; two types of silica chimneys, including mounded and splayed/funnelized chert(generally brecciated) bodies are further sorted out. The mounded chert are exitbited by domed or hummocky surfaces on the top and irregular spongy to digitiform internal fabrics; within the silica mounds, abundant original vesicles/voids and/or channels were mostly plugged by initial chalcedony, quartze crystals with minor dolomite and bladed barite crystals. Splayed/funnelized brecciated chert “intrusion” cross-cut the uppermost dolostones capping to the horizon underneath, and are directly overlain by the Niutitang phosphorite-rich black shales. Their similarities to the silica chimneys reported from the oceanic spreading centres suggest a similar origin responsible for these unique silica bodies which is also supported by the microthermonmetric data and element geochemistry. High P, Ba, Fe contents and positive correlation between Fe and TOC concentrations in the Niutitang black shales indicate a high palaeo-productivity in the Early Cambrian ocean. The low Th/U and the high V/Cr, V/Sc, V/(V+Ni) ratios in the black shales suggest an anoxic water condition during this interval. Furthermore, Positive Eu anomalies and high Ba contents in the sediments also imply a hydrothermal influence on the formation of Niutitang black shales. To better constrain the placement of deep-water successions straddling the E-C boundary and the timing of hydrothermal silica chimneys, sensitive high-resoluton ion microprobe(SHRIMP) U-Pb dating of zircon grains from tuffs within the chert succession of Liuchapo Formation at Ganziping was conducted and yields a weighted-mean 206Pb/238Pb age at 536.6±5.5Ma, younger than E-C boundary age(542.0±0.3Ma). This age combined with carbon isotopic data is then proposed to correspond to the U-Pb age of zircons(538.2±1.5Ma) from the Zhongyicun member of Meishucun Formation at Meishucun in eastern Yunna, thus, the E-C boundary in Gazngziping was placed between the Dengying formations and Liuchapo formatioms. therefore, the silica chimneys took place at the beginning of the Cambrian period. The temporal coincidence of silica chimneys and negative excursions of δ13C and δ34Spy pairs suggest hydrothermal activities were likely responsible for the isotopic changes. Under such a circumstance, vast amounts of greenhouse gases(CO2, CH4, H2S), with highly 13C-depleted carbon and 34S-depleted sulfur would be released into the ocean and atmosphere. A positive shift in δ34Scas and Δ34S values from the late Ediacaran to the Early Cambrian could be a reflection of enhanced bacterial sulfate reduction(BSR), strengthened by the intensified oceanic anoxia stimulated by hydrothermal activities. Based on the analyses of sedimentology and geochemistry, a model- “oceanic anoxia induced by hydrothermal–volcanic activies” was proposed to responsible for the formation of black chert-shales during this E-C transition. Under this case, hydrothermal-volcanic activies could release large large amount of greenhouse into atmosphere and metal micronutrients into the ocean, which may lead to global warming, stratified ocean, thereby a high palaeoproductivity; on the other hand, the massive releasing of reduced hydrothermal fluids with abundant H2S, could have in turn enhanced the ocean anoxia. All of these were favourable the for preservation of organic matter, and subsequent extensive deposition of black silica-shales.

Dating prograde fluid pulses during subduction by in situ U–Pb and oxygen isotope analysis

Keywords High-pressure fluids · Whiteschists · U–Pb dating · Oxygen isotopes · Ion microprobe · Metasomatism Introduction The subduction of crustal material to mantle depths and its chemical modification during burial and exhumation contribute to element recycling in the mantle and the formation of new crust through arc magmatism. Crustal rocks that Abstract The Dora-Maira whiteschists derive from metasomatically altered granites that experienced ultrahighpressure metamorphism at ~750 °C and 40 kbar during the Alpine orogeny. In order to investigate the P–T–time– fluid evolution of the whiteschists, we obtained U–Pb ages from zircon and monazite and combined those with trace element composition and oxygen isotopes of the accessory minerals and coexisting garnet. Zircon cores are the only remnants of the granitic protolith and still preserve a Permian age, magmatic trace element compositions and δ18O of ~10 ‰. Thermodynamic modelling of Si-rich and Si-poor whiteschist compositions shows that there are two main fluid pulses during prograde subduction between 20 and 40 kbar. In Si-poor samples, the breakdown of chlorite to garnet + fluid occurs at ~22 kbar. A first zircon rim directly overgrowing the cores has inclusions of prograde phlogopite and HREE-enriched patterns indicating zircon growth at the onset of garnet formation. A second main fluid pulse is documented close to peak metamorphic conditions in both Si-rich and Si-poor whiteschist when talc + kyanite react to garnet + coesite + fluid. A second metamorphic overgrowth on zircon with HREE depletion was observed in the Si-poor whiteschists, whereas a single metamorphic overgrowth capturing phengite and talc inclusions was observed in the Si-rich whiteschists. Garnet rims, zircon rims and monazite are in chemical and isotopic equilibrium for oxygen, demonstrating that they all formed at peak metamorphism at 35 Ma as constrained by the age of monazite (34.7 ± 0.4 Ma) and zircon rims (35.1 ± 0.8 Ma). The prograde zircon rim in Si-poor whiteschists has an age that is within error indistinguishable from the age of peak metamorphic conditions, consistent with a minimum rate of subduction of 2 cm/year for the Dora-Maira unit. Oxygen isotope values for zircon rims, monazite and garnet are equal within error at 6.4 ± 0.4 ‰, which is in line with closed-system equilibrium fractionation during prograde to peak temperatures. The resulting equilibrium Δ18Ozircon-monazite at 700 ± 20 °C is 0.1 ± 0.7 ‰. The in situ oxygen isotope data argue against an externally derived input of fluids into the whiteschists. Instead, fluidassisted zircon and monazite recrystallisation can be linked to internal dehydration reactions during prograde subduction. We propose that the major metasomatic event affecting the granite protolith was related to hydrothermal seafloor alteration post-dating Jurassic rifting, well before the onset of Alpine subduction.

Geochemistry of fallout tephra from ODP Hole 125-782A

Very rare, halogen-rich andesite melt inclusions (HRA) in bytownitic plagioclase phenocrysts (An89-90) from tephra fallout of the Izu arc volcanic front (Izu VF) provide new insights into the processes of fluid release from slab trenchward to the volcanic front in a cool subduction zone. These HRA are markedly enriched in Cl, F and Li - by factors of up to 8 (Cl, F) and 1.5 (Li) - but indistinguishable with respect to the fluid-mobile large-ion lithophile elements (LILE; K, Sr, Rb, Cs, Ba, Pb, U), rare earths (REE) or high field strength elements (HFSE) from the low-K tholeiitic magmas of the Izu VF. We suggest that the chemical signature of the HRA reflects the presence of a fluid in the mantle source that originated from the serpentinized mantle peridotite above the metacrust. This "wedge serpentinite" presumably formed by fluid infiltration beneath the forearc and was subsequently down-dragged with the slab to arc front depths. The combined evidence from the Izu VF (?110 km above slab) and the outer forearc serpentinite seamounts (~25 to 30 km above slab) suggests that the slab flux of B and Cl is highest beneath the forearc, and decreases with increasing slab depths. In contrast, the slab flux of Li is minor beneath the forearc, but increases with depth. Fluorine may behave similarly to Li, whereas the fluid-mobile LILE appear to be largely retained in the slab trenchward from the Izu VF. Consequently, the chemical signatures of both Izu trench sediments and basaltic rocks appear preserved until arc front depths.

Integrated Pb- and S-isotope investigation of sulphide minerals from the early Archaean of southwest Greenland

We present high-spatial resolution secondary ion mass spectrometry (SIMS) measurements of Pb and S isotopes in sulphides from early Archaean samples at two localities in southwest Greenland. Secondary pyrite from a 3.71 Ga sample of magnetite-quartz banded iron formation in the Isua Greenstone Belt, which has previously yielded unradiogenic Pb consistent with its ancient origin, contains sulphur with a mass independently fractionated (MIF) isotope signature (Delta(33)S =+3.3 parts per thousand). This reflects the secondary mineralization of remobilized sedimentary S carrying a component modified by photochemical reactions in the early Archaean atmosphere. It further represents one of the most extreme positive excursions so far known from the early Archaean rock record. Sulphides from a quartz-pyroxene rock and an ultramafic boudin from the island of Akilia, in the Godth (a) over circle bsfjord, have heterogeneous and generally radiogenic Pb isotopic compositions that we interpret to represent partial re-equilibration of Pb between the sulphides and whole rocks during tectonothermal events at 3.6, 2.7 and 1.6 Ga. Both these samples have Delta(33)S=0 (within analytical error) and therefore show no evidence for MIF sulphur. These data are consistent with previous interpretations that the rock cannot be proven to have a sedimentary origin. Our study illustrates that SIMS S-isotope measurements in ancient rocks can be used to elucidate early atmospheric parameters because of the ability to obtain combined S and Pb-isotope data, but caution must be applied when using such data to infer protolith. When information from geological context, petrography and chronology (i.e. by Pb isotopes) is combined and fully evaluated, Delta(33)S signatures from sulphides and their geological significance can be interpreted with a higher degree of confidence. (c) 2005 Elsevier B.V All rights reserved.

Two cycles of voluminous pyroclastic volcanism and sedimentation related to episodic granite emplacement during the late Archean : Eastern Yilgarn Craton, Western Australia

The thick piles of late-Archean volcaniclastic sedimentary successions that overlie the voluminous greenstone units of the eastern Yilgarn Craton, Western Australia, record the important transition from the cessation in mafic-ultramafic volcanism to cratonisation between about 2690 and 2655 Ma. Unfortunately, an inability to clearly subdivide the superficially similar sedimentary successions and correlate them between the various geological terranes and domains of the eastern Yilgarn Craton has led to uncertainty about the timing and nature of the region's palaeogeographic and palaeotectonic evolution. Here, we present the results of some 2025 U–Pb laser-ablation-ICP-MS analyses and 323 Sensitive High-Resolution Ion Microprobe (SHRIMP) analyses of detrital zircons from 14 late-Archean felsic clastic successions of the eastern Yilgarn Craton, which have enabled correlation of clastic successions. The results of our data, together with those compiled from previous studies, show that the post-greenstone sedimentary successions include two major cycles that both commenced with voluminous pyroclastic volcanism and ended with widespread exhumation and erosion associated with granite emplacement. Cycle One commences with an influx of rapidly reworked feldspar-rich pyroclastic debris. These units, here-named the Early Black Flag Group, are dominated by a single population of detrital zircons with an average age of 2690–2680 Ma. Thick (up to 2 km) dolerite bodies, such as the Golden Mile Dolerite, intrude the upper parts of the Early Black Flag Group at about 2680 Ma. Incipient development of large granite domes during Cycle One created extensional basins predominantly near their southeastern and northwestern margins (e.g., St Ives, Wallaby, Kanowna Belle and Agnew), into which the Early Black Flag Group and overlying coarse mafic conglomerate facies of the Late Black Flag Group were deposited. The clast compositions and detrital-zircon ages of the late Black Flag Group detritus match closely the nearby and/or stratigraphically underlying successions, thus suggesting relatively local provenance. Cycle Two involved a similar progression to that observed in Cycle One, but the age and composition of the detritus were notably different. Deposition of rapidly reworked quartz-rich pyroclastic deposits dominated by a single detrital-zircon age population of 2670–2660 Ma heralded the beginning of Cycle Two. These coarse-grained quartz-rich units, are name here the Early Merougil Group. The mean ages of the detrital zircons from the Early Merougil Group match closely the age of the peak in high-Ca (quartz-rich) granite magmatism in the Yilgarn Craton and thus probably represent the surface expression of the same event. Successions of the Late Merougil Group are dominated by coarse felsic conglomerate with abundant volcanic quartz. Although the detrital zircons in these successions have a broad spread of age, the principal sub-populations have ages of about 2665 Ma and thus match closely those of the Early Merougil Group. These successions occur most commonly at the northwestern and southeastern margins of the granite batholiths and thus are interpreted to represent resedimented units dominted by the stratigraphically underlying packages of the Early Merougil Group. The Kurrawang Group is the youngest sedimentary units identified in this study and is dominated by polymictic conglomerate with clasts of banded iron formation (BIF), granite and quartzite near the base and quartz-rich sandstone units containing detrital zircons aged up to 3500 Ma near the top. These units record provenance from deeper and/or more-distal sources. We suggest here that the principal driver for the major episodes of volcanism, sedimentation and deformation associated with basin development was the progressive emplacement of large granite batholiths. This interpretation has important implication for palaeogeographic and palaeotectonic evolution of all late-Archean terranes around the world.

Low δ18O zircons from the Bruneau-­Jarbidge eruptive center: A key to crustal anatexis along the track of the Yellowstone hotspot

The Bruneau-Jarbidge eruptive center (BJEC) in the central Snake River Plain, Idaho, USA consists of the Cougar Point Tuff (CPT), a series of ten, high-temperature (900-1000°C) voluminous ignimbrites produced over the explosive phase of volcanism (12.8-10.5 Ma) and more than a dozen equally high-temperature rhyolite lava flows produced during the effusive phase (10.5-8 Ma). Spot analyses by ion microprobe of oxygen isotope ratios in 210 zircons demonstrate that all of the eruptive units of the BJEC are characterized by zircon δ¹⁸O values ≤ 2.5‰, thus documenting the largest low δ¹⁸O silicic volcanic province known on Earth (>10⁴ km³). There is no evidence for voluminous normal δ¹⁸O magmatism at the BJEC that precedes generation of low δ¹⁸O magmas as there is at other volcanic centers that generate low δ¹⁸O magmas such as Heise and Yellowstone. At these younger volcanic centers of the hotspot track, such low δ¹⁸O magmas represent ~45 % and ~20% respectively of total eruptive volumes. Zircons in all BJEC tuffs and lavas studied (23 units) document strong δ¹⁸O depletion (median CPT δ¹⁸OZrc = 1.0‰, post-CPT lavas = 1.5‰) with the third member of the CPT recording an excursion to minimum δ¹⁸O values (δ¹⁸OZrc= -1.8‰) in a supereruption > 2‰ lower than other voluminous low δ¹⁸O rhyolites known worldwide (δ¹⁸OWR ≤0.9 vs. 3.4‰). Subsequent units of the CPT and lavas record a progressive recovery in δ¹⁸OZrc to ~2.5‰ over a ~ 4 m.y. interval (12 to 8 Ma). We present detailed evidence of unit-to-unit systematic patterns in O isotopic zoning in zircons (i.e. direction and magnitude of Δcore-rim), spectrum of δ¹⁸O in individual units, and zircon inheritance patterns established by re-analysis of spots for U-Th-Pb isotopes by LA-ICPMS and SHRIMP. In conjunction with mineral thermometry and magma compositions, these patterns are difficult to reconcile with the well-established model for "cannibalistic" low δ¹⁸O magma genesis at Heise and Yellowstone. We present an alternative model for the central Snake River Plain using the modeling results of Leeman et al. (2008) for ¹⁸O depletion as a function of depth in a mid-upper crustal protolith that was hydrothermally altered by infiltrating meteoric waters prior to the onset of silicic magmatism. The model proposes that BJEC silicic magmas were generated in response to the propagation of a melting front, driven by the incremental growth of a vast underlying mafic sill complex, over a ~5 m.y. interval through a crustal volume in which a vertically asymmetric δ¹⁸OWR gradient had previously developed that was sharply inflected from ~ -1 to 10‰ at mid-upper crustal depths. Within the context of the model, data from BJEC zircons are consistent with incremental melting and mixing events in roof zones of magma reservoirs that accompany surfaceward advance of the coupled mafic-silicic magmatic system.