Geochronology of the lpanema layered mafic-ultramafic complex, Minas gerais, Brazil: Evidence of extension at the meso-neoproterozoic time boundary

The mafic/ultramafic Ipanema Layered Complex (ILC), Minas Gerais Brazil, consists of seven individual bodies. These units crosscut polyphase orthogneisses and interlayered paragneisses of the Paleoproterozoic Juiz de Fora Complex. Intrusive granitoids tectonically related to [lie Neoproterozoie Aracuai orogen are also present in the study area.A Sm-Nd whole-rock linear array for seven samples metapyoxenites, metaperidotiles, metagabbro. and meta-anorthosite) from the Santa Cruz massif, the largest body of the ILC. suggest that it was emplaced at 1104 +/- 78 Ma the original magma was derived from a depleted mantle source (epsilon(Ndt)= +3.8). U-Pb single-grain zircon stud of a meta-anorthosite yields all upper intercept age of 1719 +/- 4 Ma, which is interpreted to represent inheritance. The lower intercept at 630+/-3 Ma indicates (hat a Neoproterozoic tectonothermal episode overprinted the ILC, this event occurred under upper-amphiolite-, to granulite-facies conditions. The 630 Ma episode is consistent with the timing of regional metamorphism and deformation of the adjacent Aracuai orogen (Brasiliano collage). Emplacement of the ILC and other coeval metamafies and meta-ultramafics (of alkaline affinity) in the re, oil is attributed to early extension tectonics, accompanying accretion of the Rodinia super- continent during the Mesoproterozoic-Neoproterozoic time boundary.

Evolução geológica dos batólitos granitóides neoproterozóicos do sudeste do estado de são paulo

The Neoproterozoic granitogenesis related to the Central Mantiqueira Province comprise the calc alkaline to alkaline granitoid complexes of Sorocaba, San Francisco, São Roque, Ibiúna and Piedade. These complexes occur in a ruptil tectonic to tardi (Sn+3) event. The emplacement of the different facies in transtractives structures of the pull-apart type are characterized in the area by the main transcurrent shear zones of Taxaquara-Pirapora, Itu-Jundiuvira, Moreiras, Cangüera and Caucáia of ENE-WSW general direction. The massifs present complex internal architecture characterized by intrusions in restrict initial phase of intermediate equigranular nature. Also present a main phase of porfiroid monzo and sienogranite that fragments the previous phase, followed by lateral accretion of equi to inequigranular material, and in some cases by the accretion of late phases of circular bodies of porfiroid rapakivi granites, and a late to final phase of aplitic to pegmatitic composition. This magmatism grew with the intrusions of successive magmatic pulses, partially controlled by many reactivations of the shear zones. The REE also suggest that the magmatic phases are similar, synchronous and repetitive in four of the complexes in both domains, present in the São Francisco Complex. The crystallization starts from accretion processes, but compositionally quite different from the others. The variation in compositions and ages (TDM) for these granites reflect the derivation from different sources developed under different magmatic conditions, followed by processes of contamination that frequently occur in the crust.

Crustal growth in the 3.4-2.7Ga São José de Campestre Massif, Borborema Province, NE Brazil

The Archean (3.45-2.70Ga) rocks of the São José do Campestre Massif (SJCM) in the Borborema Province (NE Brazil) make up a small area (~6000km2) and are composed of granitoids and metasupracrustal rocks that define a complex magmatic and deformational history. The massif provides the opportunity to study mantle- and crustal-derived magmas generated since the Palaeoarchean. The orthogneisses of the SJCM are composed of: (1) tonalite to granodiorite with diorite enclaves (Bom Jesus gneiss, 3412±8Ma; TDM Nd model ages from 4.1 to 3.5Ga and negative epsilon Nd values); (2) biotite and ferroan-diopside monzogranite (Presidente Juscelino complex, 3356±21Ma and 3251±44Ma; TDM model ages range from 4.1 to 3.4Ga and epsilon Nd values that are slightly positive to negative); (3) hornblende tonalite to granodiorite (Brejinho complex, 3333±77Ma and 3187±8Ma; dominantly positive epsilon Nd values and TDM ages from 3.6 to 3.2Ga); (4) biotite monzogranite (São Pedro do Potengi gneiss, 3120±22Ma; TDM =3.5Ga; negative epsilon Nd value); (5) ferroan-diopside-grossular anorthosite and metagabbro (Senador Elói de Souza complex, 3033±3Ma); and (6) quartz diorite to syenogranite (São José do Campestre complex; 2685±9Ma and 2655±4Ma; negative epsilon Nd values and TDM ages from 3.9 to 3.3Ga). The orthogneisses are subalkaline to faintly alkaline, magnesian to ferroan, M- and I-type granitoids that follow either the K-enrichment or the trondhjemite trends. Each group has a subset with REE characteristics similar to Archean TTG and another that is analogous to Phanerozoic granitoids. They have negative Ta-Nb and Ti anomalies and have trace element contents of granitoids from subduction zones. Geochemical and Nd isotope data suggest that subducted oceanic crust and a depleted and metasomatised mantle wedge both acted as the magma sources. We propose a convergent tectonic model in which hybridisation of the upper mantle occurs through interactions with adakitic or trondhjemitic melts and recycling of earlier crust. The results imply that both the subducted oceanic crust and the mantle wedge played major roles in continent formation throughout successive episodes of arc accretion in Palaeo- and Mesoarchean times. The Archean rocks of the SJCM shares some similarities with the Pilbara, Kaapvaal, West African, and São Francisco cratons. However, the most reliable comparisons with the SJCM are with the neighbouring basement of the Nigeria and Cameroon shields. © 2012 Elsevier B.V.

Petrogenesis of syntectonic granites emplaced at the transition from thrusting to transcurrent tectonics in post-collisional setting: Whole-rock and Sr-Nd-Pb isotope geochemistry in the Neoproterozoic Quatro Ilhas and Mariscal Granites, Southern Brazil

The Neoproterozoic post-collisional period in southern Brazil (650-580 Ma) is characterized by substantial volumes of magma emplaced along the active shear zones that compose the Southern Brazilian Shear Belt. The early-phase syntectonic magmatism (630-610 Ma) is represented by the porphyritic, high-K, metaluminous to peraluminous Quatro Ilhas Granitoids and the younger heterogranular, slightly peraluminous Mariscal Granite. Quatro II has Granitoids include three main petrographic varieties (muscovite-biotite granodiorite mbg; biotite monzogranite - bmz: and leucogranite - lcg) that, although sharing some significant geochemical characteristics, are not strictly comagmatic, as shown by chemical and Sr-Nd-Pb isotope data. The most primitive muscovite-biotite granodiorite was produced by contamination of more mafic melts (possibly with some mantle component) with peraluminous crustal melts; the biotite monzogranite, although more felsic, has higher Ca, MgO,TiO2 and Ba, and lower K2O, FeOt, Sr and Rb contents, possibly reflecting some mixing with coeval mafic magmas of tholeiitic affinity; the leucogranite may be derived from pure crustal melts. The Mariscal Granite is formed by two main granite types which occur intimately associated in the same pluton, one with higher K (5-6.5 wt.% K2O) high Rb and lower CaO, Na2O, Ba and Zr as compared to the other (3-5 wt.% of K2O). The two Mariscal Granite varieties have compositional correspondence with fine-grained granites (fgg) that occur as tabular bodies which intruded the Quatro Ilhas Granoitoids before they were fully crystallized, and are inferred to correspond to the Mariscal Granite feeders, an interpretation that is reinforced by similar U-Pb zircon crystallization ages. The initial evolution of the post-collisional magmatism, marked by the emplacement of the Quatro Ilhas Granitoids varieties, activated sources that produced mantle and crustal magmas whose emplacement was controlled both by flat-lying and transcurrent structures. The transition from thrust to transcurrent-related tectonics coincides with the increase in the proportion of crustal-derived melts. The transcurrent tectonics seems to have played an essential role in the generation of mantle-derived magmas and may have facilitated their interaction with crustal melts which seem to be to a large extent the products of reworking of orthogneiss protoliths. (C) 2012 Elsevier B.V. All rights reserved.

Ages (U-Pb SHRIMP and LA ICPMS) and stratigraphic evolution of the Neoproterozoic volcano-sedimentary successions from the extensional Camaqua Basin, Southern Brazil

During the Ediacaran, southern Brazil was the site of multiple episodes of volcanism and sedimentation, which are best preserved in the 3000 km(2) Camaqua Basin. The interlayered sedimentary and volcanic rocks record tectonic events and paleoenvironmental changes in a more than 10 km-thick succession. In this contribution, we report new U-Pb and Sm-Nd geochronological constraints for the 605 to 580 Ma Born Jardim Group, the 570 Ma Acampamento Velho Formation, and a newly-recognized 544 Ma volcanism. Depositional patterns of these units reveal the transition from a restricted, fault-bounded basin into a wide, shallow basin. The expansion of the basin and diminished subsidence rates are demonstrated by increasing areal distribution and compressed isopachs and increasing onlap of sediments onto the basement to the west. The Sm-Nd isotopic composition of the volcanic rocks indicates mixed sources, including crustal rocks from the adjacent basement. Both Neoproterozoic and Paleoproterozoic sources are indicated for the western part of the basin, whereas only the older Paleoproterozoic signature can be discerned in the eastern part of the basin. (C) 2011 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

The pre-alpine evolution of the basement of the Pelagonian Zone and the Vardar Zone, Greece

The Pelagonian Zone and the Vardar Zone in Greece represent the western part of the Hellenide hinterland (Internal Hellenides). While the Pelagonian Zone comprises predominantly crystalline basement and sedimentary cover rocks, the Vardar Zone has long been regarded as an ophiolite-decorated suture zone separating the Pelagonian Zone from the Serbo-Macedonian Massif to the east. Felsic basement rocks from both areas, with the main focus put on the Pelagonian Zone, were dated in order to identify the major crust-forming episodes and to improve the understanding of the evolutionary history of the region. The interpretation of the single-zircon geochronology results was aided by geochemical investigations. The majority of the basement rocks from the Pelagonian Zone yielded Permo-Carboniferous intrusion ages around 300 Ma, underlining the importance of this crust-forming event for the Internal Hellenides of Greece. Geochemically these basement rocks are classified as subduction-related granitoids, which formed in an active continental margin setting. An important result was the identification of a Precambrian crustal unit within the crystalline basement of the Pelagonian Zone. Orthogneisses from the NW Pelagonian Zone yielded Neoproterozoic ages of c. 700 Ma and are so far the oldest known rocks in Greece. These basement rocks, which are also similar to active margin granitoids, were interpreted as remnants of a terrane, the Florina Terrane, which can be correlated to a Pan-African or Cadomian arc. Since the gneisses contain inherited zircons of Middle to Late Proterozoic ages, the original location of the Florina Terrane was probably at the northwestern margin of Gondwana. In the Vardar Zone an important phase of Upper Jurassic felsic magmatism is documented by igneous formation ages ranging from 155 to 164 Ma. The chemical and isotopic composition of these rocks is also in accord with their formation in a volcanic-arc setting at an active continental margin. Older continental material incorporated in the Vardar Zone is documented by 319-Ma-old gneisses and by inherited zircons of mainly Middle Palaeozoic ages. The prevalence of subduction-related igneous rocks indicates that arc formation and accretion orogeny were the most important processes during the evolution of this part of the Internal Hellenides. The geochronological results demonstrate that most of the Pelagonian Zone and the Vardar Zone crystalline basement formed during distinct pre-Alpine episodes at c. 700, 300 and 160 Ma with a predominance of the Permo-Carboniferous magmatic phase.

Estrutura crustal e mantélica da província Borborema através de funções do receptor e dispersão de ondas superficiais

The Borborema Province, located in northeastern Brazil, has a basement of Precambrian age and a tectonic framework structured at the Neoproterozoic (740-560 Ma). After separation between South America and Africa during the Mesozoic, a rift system was formed, giving rise to a number of marginal and inland basins in the Province. After continental breakup, episodes of volcanism and uplift characterized the evolution of the Province. Plateau uplift was initially related to magmatic underplating of mafic material at the base of the crust, perhaps related to the generation of young continental plugs (45-7 Ma) along the Macau-Queimadas Alignment (MQA), due to a small-scale convection at the continental edge. The goal of this study is to investigate the causes of intra-plate uplift and its relationship to MQA volcanism, by using broadband seismology and integrating our results with independent geophysical and geological studies in the Borborema Province. The investigation of the deep structure of the Province with broadband seismic data includes receiver functions and surface-wave dispersion tomography. Both the receiver functions and surface-wave dispersion tomography are methods that use teleseismic events and allow to develop estimates of crustal parameters such as crustal thickness, Vp/Vs ratio, and S-velocity structure. The seismograms used for the receiver function work were obtained from 52 stations in Northeast Brazil: 16 broadband stations from the RSISNE network (Rede Sismográfica do Nordeste do Brasil), and 21 short-period and 6 broadband stations from the INCT-ET network (Instituto Nacional de Ciência e Tecnologia – Estudos Tectônicos). These results add signifi- cantly to previous datasets collected at individual stations in the Province, which include station RCBR (GSN - Global Seismic Network), stations CAUB and AGBL (Brazilian Lithosphere Seismic Project IAG/USP), and 6 other broadband stations that were part of the Projeto Milênio - Estudos geofísicos e tectônicos na Província Borborema/CNPq. For the surface-wave vii tomography, seismograms recorde at 22 broadband stations were utilized: 16 broadband stations from the RSISNE network and 6 broadband stations from the Milênio project. The new constraints developed in this work include: (i) estimates of crustal thickness and bulk Vp/Vs ratio for each station using receiver functions; (ii) new measurements of surfassewave group velocity, which were integrated to existing measurementes from a continental-scale tomography for South America, and (iii) S-wave velocity models (1D) at various locations in the Borborema Province, developed through the simultaneous inversion of receiver functions and surface-wave dispersion velocities. The results display S-wave velocity structure down to the base of the crust that are consistent with the presence of a 5-7.5 km thick mafic layer. The mafic layer was observed only in the southern portion of the Plateau and absent in its northern portion. Another important observation is that our models divide the plateau into a region of thin crust (northern Plateau) and a region of thick crust (southern Plateau), confirming results from independent refraction surveys and receiver function analyses. Existing models of plateau uplift, nonetheless, cannot explain all the new observations. It is proposed that during the Brazilian orogeny a layer of preexisting mafic material was delaminated, as a whole or in part, from the original Brasiliano crust. Partial delamination would have happened in the southern portion of the plateau, where independent studies found evidence of a more resistant rheology. During Mesozoic rifting, thinning of the crust around the southern Plateau would have formed the marginal basins and the Sertaneja depression, which would have included the northern part of the Plateau. In the Cenozoic, uplift of the northern Plateau would have occurred, resulting in a northern Plateau without mafic material at the base of the crust and a southern Plateau with partially delaminated mafic layer.

New insights into Crustal Contributions to Large-volume Rhyolite Generation in the Mid-Tertiary Sierra Madre Occidental Province, Mexico, revealed by U Pb Geochronology

Voluminous (≥3·9 × 105 km3), prolonged (∼18 Myr) explosive silicic volcanism makes the mid-Tertiary Sierra Madre Occidental province of Mexico one of the largest intact silicic volcanic provinces known. Previous models have proposed an assimilation–fractional crystallization origin for the rhyolites involving closed-system fractional crystallization from crustally contaminated andesitic parental magmas, with <20% crustal contributions. The lack of isotopic variation among the lower crustal xenoliths inferred to represent the crustal contaminants and coeval Sierra Madre Occidental rhyolite and basaltic andesite to andesite volcanic rocks has constrained interpretations for larger crustal contributions. Here, we use zircon age populations as probes to assess crustal involvement in Sierra Madre Occidental silicic magmatism. Laser ablation-inductively coupled plasma-mass spectrometry analyses of zircons from rhyolitic ignimbrites from the northeastern and southwestern sectors of the province yield U–Pb ages that show significant age discrepancies of 1–4 Myr compared with previously determined K/Ar and 40Ar/39Ar ages from the same ignimbrites; the age differences are greater than the errors attributable to analytical uncertainty. Zircon xenocrysts with new overgrowths in the Late Eocene to earliest Oligocene rhyolite ignimbrites from the northeastern sector provide direct evidence for some involvement of Proterozoic crustal materials, and, potentially more importantly, the derivation of zircon from Mesozoic and Eocene age, isotopically primitive, subduction-related igneous basement. The youngest rhyolitic ignimbrites from the southwestern sector show even stronger evidence for inheritance in the age spectra, but lack old inherited zircon (i.e. Eocene or older). Instead, these Early Miocene ignimbrites are dominated by antecrystic zircons, representing >33 to ∼100% of the dated population; most antecrysts range in age between ∼20 and 32 Ma. A sub-population of the antecrystic zircons is chemically distinct in terms of their high U (>1000 ppm to 1·3 wt %) and heavy REE contents; these are not present in the Oligocene ignimbrites in the northeastern sector of the Sierra Madre Occidental. The combination of antecryst zircon U–Pb ages and chemistry suggests that much of the zircon in the youngest rhyolites was derived by remelting of partially molten to solidified igneous rocks formed during preceding phases of Sierra Madre Occidental volcanism. Strong Zr undersaturation, and estimations for very rapid dissolution rates of entrained zircons, preclude coeval mafic magmas being parental to the rhyolite magmas by a process of lower crustal assimilation followed by closed-system crystal fractionation as interpreted in previous studies of the Sierra Madre Occidental rhyolites. Mafic magmas were more probably important in providing a long-lived heat and material flux into the crust, resulting in the remelting and recycling of older crust and newly formed igneous materials related to Sierra Madre Occidental magmatism.

Neoproterozoic glacial deposits of the Kimberly Region and northwestern Northern Territory, Australia

Neoproterozoic glacigenic formations are preserved in the Kimberley region and northwestern Northern Territory of northern Australia. They are distributed in the west Kimberley adjacent to the northern margins of the King Leopold Orogen, the Mt Ramsay area at the junction of the King Leopold and Halls Creek Orogens, and the east Kimberley, adjacent to the eastern margin of the Halls Creek Orogen. Small outlier glacigenic deposits are preserved in the Litchfield Province, Northern Territory (Uniya Formation) and Georgina Basin, western Queensland (Little Burke Formation). Glacigenic strata comprise diamictite, conglomerate, sandstone and pebbly mudstone and characterize the Walsh, Landrigan and Fargoo/Moonlight Valley formations. Thin units of laminated dolomite sit conformably at the top of the Walsh, Landrigan and Moonlight Valley formations. Glacigenic units are also interbedded with the carbonate platform deposits of the Egan Formation and Boonall Dolomite. δ13C data are available for all carbonate units. There is no direct chronological constraint on these successions. Dispute over regional correlation of the Neoproterozoic succession has been largely resolved through biostratigraphic, chemostratigraphic and lithostratigraphic analysis. However, palaeomagnetic results from the Walsh Formation are inconsistent with sedimentologically based correlations. Two stratigraphically defined glaciations are preserved in northwestern Australia: the ‘Landrigan Glaciation’, characterized by southwest-directed continental ice-sheet movement and correlated with late Cryogenian glaciation elsewhere in Australia and the world; and, the ‘Egan Glaciation’, a more localized glaciation of the Ediacaran Period. Future research focus should include chronology, palaeomagnetic constraint and tectonostratigraphic controls on deposition.

Discriminating stromatolite formation modes using rare earth element geochemistry: Trapping and binding versus in situ precipitation of stromatolites from the Neoproterozoic Bitter Springs Formation, Northern Territory, Australia

Stromatolites consist primarily of trapped and bound ambient sediment and/or authigenic mineral precipitates, but discrimination of the two constituents is difficult where stromatolites have a fine texture. We used laser ablation-inductively coupled plasma-mass spectrometry to measure trace element (rare earth element – REE, Y and Th) concentrations in both stromatolites (domical and branched) and closely associated particulate carbonate sediment in interspaces (spaces between columns or branches) from bioherms within the Neoproterozoic Bitter Springs Formation, central Australia. Our high resolution sampling allows discrimination of shale-normalised REE patterns between carbonate in stromatolites and immediately adjacent, fine-grained ambient particulate carbonate sediment from interspaces. Whereas all samples show similar negative La and Ce anomalies, positive Gd anomalies and chondritic Y/Ho ratios, the stromatolites and non-stromatolite sediment are distinguishable on the basis of consistently elevated light REEs (LREEs) in the stromatolitic laminae and relatively depleted LREEs in the particulate sediment samples. Additionally, concentrations of the lithophile element Th are higher in ambient sediment samples than in stromatolites, consistent with accumulation of some fine siliciclastic detrital material in the ambient sediment but a near absence in the stromatolites. These findings are consistent with the stromatolites consisting dominantly of in situ carbonate precipitates rather than trapped and bound ambient sediment. Hence, high resolution trace element (REE + Y, Th) geochemistry can discriminate fine-grained carbonates in these stromatolites from coeval non-stromatolitic carbonate sediment and demonstrates that the sampled stromatolites formed primarily from in situ precipitation, presumably within microbial mats/biofilms, rather than by trapping and binding of ambient sediment. Identification of the source of fine carbonate in stromatolites is significant, because if it is not too heavily contaminated by trapped ambient sediment, it may contain geochemical biosignatures and/or direct evidence of the local water chemistry in which the precipitates formed.

Role of kilometer-scale weak circular heterogeneities on upper crustal deformation patterns: Evidence from scaled analogue modeling and the Sudbury Basin, Canada

The Sudbury Basin is a non-cylindrical fold basin occupying the central portion of the Sudbury Impact Structure. The impact structure lends itself excellently to explore the structural evolution of continental crust containing a circular region of long-term weakness. In a series of scaled analogue experiments various model crustal configurations were shortened horizontally at a constant rate. In mechanically weakened crust, model basins formed that mimic several first-order structural characteristics of the Sudbury Basin: (1) asymmetric, non-cylindrical folding of the Basin, (2) structures indicating concentric shortening around lateral basin termini and (3) the presence of a zone of strain concentration near the hinge zones of model basins. Geometrically and kinematically this zone corresponds to the South Range Shear Zone of the Sudbury Basin. According to our experiments, this shear zone is a direct mechanical consequence of basin formation, rather than the result of thrusting following folding. Overall, the models highlight the structurally anomalous character of the Sudbury Basin within the Paleoproterozoic Eastern Penokean Orogen. In particular, our models suggest that the Basin formed by pure shear thickening of crust, whereas transpressive deformation prevailed elsewhere in the orogen. The model basin is deformed by thickening and non-cylindrical synformal buckling, while conjugate transpressive shear zones propagated away from its lateral tips. This is consistent with pure shear deformation of a weak circular inclusion in a strong matrix. The models suggest that the Sudbury Basin formed as a consequence of long-term weakening of the upper crust by meteorite impact.

The causes of sinuous crustal-scale deformation patterns in hot orogens: Evidence from scaled analogue experiments and the southern Central Andes

The cause of upper-crustal segmentation into rhomb-shaped, shear zone-bound domains associated with contractional sedimentary basins in hot, wide orogens is not well understood. Here we use scaled multilayered analogue experiments to investigate the role of an orogen-parallel crustal-strength gradient on the formation of such structures. We show that the aspect ratio and size of domains, the sinuous character and abundance of transpressional shear zones vary with the integrated mechanical strength of crust. Upper-crustal deformation patterns and the degree of strain localization in the experiments are controlled by the ratio between the brittle and ductile strength in the model crust as well as gradients in tectonic and buoyancy forces. The experimental results match the first-order kinematic and structural characteristics of the southern Central Andes and provide insight on the dynamics of underlying deformation patterns in hot, wide orogens.

Heat-producing element enrichment in granitic rocks, the role of crustal composition and evolution

This study greatly enhanced our knowledge of the potential for geothermal energy development in Queensland as a viable clean energy source in the coming decades. Key outcomes of the project were understanding the first-order controls on the concentration of the heat-producing elements: uranium, thorium and potassium in granitic rocks, and constraining where rocks with the greatest heat-producing potential lie at depth in Queensland. Importantly, new temperature and heat flow maps for southwest Queensland were developed that will greatly assist future exploration efforts.

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.