Geochemical and Isotopic Heterogeneities along an Island Arc-Spreading Ridge Intersection: Evidence from the Lewis Hills, Bay of Islands Ophiolite, Newfoundland

This study focuses on the origin of magma heterogeneity andthe genesis of refractory, boninite-type magmas along an arc–ridgeintersection, exposed in the Lewis Hills (Bay of Islands Ophiolite).The Lewis Hills contain the fossil fracture zone contact betweena split island arc and its related marginal oceanic basin. Threetypes of intrusions, which are closely related to this narrowtectonic boundary, have been investigated. Parental melts inequilibrium with the ultramafic cumulates of the PyroxeniteSuite are inferred to have high MgO contents and low Al₂O₃,Na₂O and TiO₂ contents. The trace element signatures of thesePyroxenite Suite parental melts indicate a re-enriched, highlydepleted source with 0·1 x mid-ocean ridge basalt (MORB)abundances of the heavy rare earth elements (HREE). Initial_Nd values of the Pyroxenite Suite range from -1·5 to+0·6, which overlap those observed for the island arc.Furthermore, the Pyroxenite Suite parental melts bear strongsimilarities to boninite-type equilibrium melts from islandarc-related pyroxenitic dykes and harzburgites. Basaltic dykessplit into two groups. Group I dykes have 0·6 x MORBabundances of the HREE, and initial _Nd values ranging from +5·4to +7·5. Thus, they have a strong geochemical affinitywith basalts derived from the marginal basin spreading ridge.Group II dykes have comparatively lower trace element abundances(0·3 x MORB abundances of HREE), and slightly lower initial_Nd values (+5·4 to +5·9). The geochemical characteristicsof the Group II dykes are transitional between those of GroupI dykes and the Pyroxenite Suite parental melts. Cumulates fromthe Late Intrusion Suite are similarly transitional, with _Ndvalues ranging from +2·9 to +4·6. We suggest thatthe magma heterogeneity observed in the Lewis Hills is due tothe involvement of two compositionally distinct mantle sources,which are the sub-island lithospheric mantle and the asthenosphericmarginal basin mantle. It is likely that the refractory, boninite-typeparental melts of the Pyroxenite Suite result from remeltingof the sub-arc lithospheric mantle at an arc–ridge intersection.Furthermore, it is suggested that the thermal-dynamic conditionsof the transtensional transform fault have provided the prerequisitefor generating magma heterogeneity, as a result of mixing relationshipsbetween arc-related and marginal basin-related magmas. KEY WORDS: Bay of Islands ophiolite; transform (arc)–ridge intersection; boninites; rare earth elements, Nd isotopes     

Petrological, Geochemical and Isotopic Constraints on the Origin of the Harzburg Intrusion, Germany

We present mineralogical, petrological and geochemical datato constrain the origin of the Harzburg mafic–ultramaficintrusion. The intrusion is composed mainly of mafic rocks rangingfrom gabbronorite to quartz diorite. Ultramafic rocks are veryrare in surface outcrops. Dunite is observed only in deepersections of the Flora I drill core. Microgranitic (fine-grainedquartz-feldspathic) veins found in the mafic and ultramaficrocks result from contamination of the ultramafic magmas bycrustal melts. In ultramafic and mafic compositions cumulatetextures are widespread and filter pressing phenomena are obvious.The order of crystallization is olivine pargasite, phlogopite,spinel plagioclase, orthopyroxene plagioclase, clinopyroxene.Hydrous minerals such as phlogopite and pargasite are essentialconstituents of the ultramafic cumulates. The most primitiveolivine composition is Fo_89·5 with 0·4 wt % NiO,which indicates that the olivine may have been in equilibriumwith primitive mantle melts. Coexisting melt compositions estimatedfrom this olivine have mg-number = 71. The chemical varietyof the rocks constituting the intrusion and the mg-number ofthe most primitive melt allow an estimation of the approximatecomposition of the mantle-derived primary magma. The geochemicalcharacteristics of the estimated magma are similar to thoseof an island-arc tholeiite, characterized by low TiO₂ and alkalisand high Al₂O₃. Geochemical and Pb, Sr and Nd isotope data demonstratethat even the most primitive rocks have assimilated crustalmaterial. The decoupling of Sr from Nd in some samples demonstratesthe influence of a fluid that transported radiogenic Sr. Leadof crustal origin from two isotopically distinct reservoirsdominates the Pb of all samples. The ultramafic rocks and thecumulates best reflect the initial isotopic and geochemicalsignature of the parent magma. Magma that crystallized in theupper part of the chamber was more strongly affected by assimilatedmaterial. Petrographic, geochemical and isotope evidence demonstratesthat during a late stage of crystallization, hybrid rocks formedthrough the mechanical mixing of early cumulates and melts withstrong crustal contamination from the upper levels of the magmachamber. KEY WORDS: Harzburg mafic–ultramafic intrusion; Sr–Nd–Pb isotopes; magma evolution; crustal contamination     

The Petrology and Geochemistry of the Aniakchak Caldera-forming Ignimbrite, Aleutian Arc, Alaska

Aniakchak caldera, Alaska, produced a compositionally heterogeneousignimbrite 3400 years ago, which changes from rhyodacitic atthe base to andesitic at the top of the eruptive sequence. Interpretationsof compositionally heterogeneous ignimbrites typically includeeither in situ fractional crystallization of mafic magma andgeneration of a stratified magma body or replenishment of asilicic magma chamber by mafic inputs. Another possibility,silicic replenishment of a more mafic chamber, exists. Geochemicalcharacteristics of the caldera-forming rhyodacite and severallate pre-caldera rhyodacites indicate independent origins foreach, within a maximum of 5000 years prior to caldera formation.Isotopic considerations preclude derivation of the caldera-formingrhyodacite from the caldera-forming andesite. However, the caldera-formingrhyodacite can be explained as the residual liquid of a mostlycrystallized basalt, with addition of crustal material. TheAniakchak andesite probably formed in a shallow chamber by successivemixing events involving small volumes of basalt and rhyodacite,together with contamination. The pre-caldera rhyodacites representerupted portions of intruding silicic magma, whereas anotherportion homogenized with the resident mafic magma. The caldera-formingevent reflects a large influx of rhyodacite, which erupted beforesignificant mixing occurred and also triggered draining of muchof the andesitic magma from the chamber. KEY WORDS: Aniakchak; caldera-forming eruption; geochemistry; ignimbrite; silicic replenishment     

Magmatic Evolution of the Skye Igneous Centre, Western Scotland: Modelling of Assimilation, Recharge and Fractional Crystallization

The Skye igneous centre, forming part of the British Tertiarymagmatic province, developed over a 7 Myr period (61–54Ma) and is characterized by a complex suite of lavas, hypabyssaland intrusive rocks of picritic to granitic composition. Theintrusion of magma from mantle to crust at 2x10^–3km³/yr(6 Mt/yr) advected magmatic heat of roughly 0·2 GW averagedover the period of magmatism supporting an ‘excess’heat flux of about 130 mW/m², or about twice the present-dayaverage continental heat flow. The volume of new crust generatedat Skye (15000 km³) spread over the present-day area of Skyecorresponds to 9 km of new crust. The geochemical evolutionof the Skye magmatic system is constrained using the Energy-ConstrainedRecharge, Assimilation, and Fractional Crystallization (EC-RAFC)model to understand variations in the Sr- and Pb-isotopic andSr trace-element composition of the exposed magmatic rocks withtime. The character (composition and specific enthalpy) of bothassimilant and recharge magma appears to change systematicallyup-section, suggesting that the magma reservoirs migrated toprogressively shallower levels as the system matured. The modelof the magma transport system that emerges is one in which magmabatches are stored initially at lower-crustal levels, wherethey undergo RAFC evolution. Residual magma from this stagethen migrates to shallower levels, where mid-crustal wall rockis assimilated; the recharge magma at this level is characterizedby an increasingly crustal signature. For some of the stratigraphicallyyoungest rocks, the data suggest that the magma reservoirs ascendedinto, and interacted with, upper-crustal Torridonian metasediments. KEY WORDS: assimilation; EC-RAFC model; geochemical modelling; magma recharge; Skye magmatism     

Isotopic Studies of Processes in Mafic Magma Chambers: I. The Kiglapait Intrusion, Labrador

The initial Nd and Sr isotopic ratios of cumulate rocks andminerals of the Kiglapait intrusion are pertinent as indicatorsof the processes that affected the Kiglapair magma while itwas resident in the crust. A Sm–Nd mineral isochron indicatesthat the crystallization age of the intrusion is 1305?22 millionyears. The initial _Nd values range from—1?6 in the LowerZone to 6 in the Upper Zone and correlate with the anorthitecontent of plagioclase. The initial ₈₇Sr/⁸⁶Sr ratios are 0.70407to 0.70433 in the Lower Zone, and increase monotonically up-section(decreasing plagioclase anorthite content) to 0.7068 in theUpper Zone. These variations are attributed to assimilationof roof rock concurrently with crystallization. It is evidentthat replenishment of the chamber with uncontaminated magmawas important during the accumulation of the lower zone rocks,but did not occurduring crystallization of the Upper Zone. Amathematical model relating isotope ratio shifts to the relativerates of crystallization, assimilation, and replenishment ispresented. It is estimated that the rate of assimilation wasbetween 0.01 and 0.04 of the crystallization rate, and thatduring the accumulation of the Lower Zone, the rate of replenishmentwas about half of the crystallization rate. A formulation for‘assimilation efficiency’is presented that relatedthe actual amount of assimilation (determined from the isotopedata) to the thermodynamic maximum allowable amount. The assimilationefficiency of the Kiglapait intrusion is not tightly constrained,but appears to be of the order of 0.01 near the end of crystallizationand perhaps as high as 0.07 near the beginning. Further quantitativeestimates of assimilation efficiency in different intrusionsemplaced under different conditions may aid in understandingassimilation mechanisms and in assessing the role of assimilationin magmatic evolution in general. Many layered intrusions studiedto date have initial Nd isotope ratios close to the chondriticvalue (_Nd= 0), suggestive of mantle magma sources that are chemicallyundepleted, and thus different from the sources of mid-oceanridge basalt and island-arec basalt.     

Pressures of Crystallization of Icelandic Magmas

Iceland lies astride the Mid-Atlantic Ridge and was createdby seafloor spreading that began about 55 Ma. The crust is anomalouslythick (20–40 km), indicating higher melt productivityin the underlying mantle compared with normal ridge segmentsas a result of the presence of a mantle plume or upwelling centeredbeneath the northwestern edge of the Vatnajökull ice sheet.Seismic and volcanic activity is concentrated in 50 km wideneovolcanic or rift zones, which mark the subaerial Mid-AtlanticRidge, and in three flank zones. Geodetic and geophysical studiesprovide evidence for magma chambers located over a range ofdepths (1·5–21 km) in the crust, with shallow magmachambers beneath some volcanic centers (Katla, Grimsvötn,Eyjafjallajökull), and both shallow and deep chambers beneathothers (e.g. Krafla and Askja). We have compiled analyses ofbasalt glass with geochemical characteristics indicating crystallizationof ol–plag–cpx from 28 volcanic centers in the Western,Northern and Eastern rift zones as well as from the SouthernFlank Zone. Pressures of crystallization were calculated forthese glasses, and confirm that Icelandic magmas crystallizeover a wide range of pressures (0·001 to 1 GPa), equivalentto depths of 0–35 km. This range partly reflects crystallizationof melts en route to the surface, probably in dikes and conduits,after they leave intracrustal chambers. We find no evidencefor a shallow chamber beneath Katla, which probably indicatesthat the shallow chamber identified in other studies containssilica-rich magma rather than basalt. There is reasonably goodcorrelation between the depths of deep chambers (> 17 km)and geophysical estimates of Moho depth, indicating that magmaponds at the crust–mantle boundary. Shallow chambers (<7·1 km) are located in the upper crust, and probablyform at a level of neutral buoyancy. There are also discretechambers at intermediate depths (11 km beneath the rift zones),and there is strong evidence for cooling and crystallizing magmabodies or pockets throughout the middle and lower crust thatmight resemble a crystal mush. The results suggest that themiddle and lower crust is relatively hot and porous. It is suggestedthat crustal accretion occurs over a range of depths similarto those in recent models for accretionary processes at mid-oceanridges. The presence of multiple stacked chambers and hot, porouscrust suggests that magma evolution is complex and involvespolybaric crystallization, magma mixing, and assimilation. KEY WORDS: Iceland rift zones; cotectic crystallization; pressure; depth; magma chamber; volcanic glass     

A Textural Record of Solidification and Cooling in the Skaergaard Intrusion, East Greenland

The clinopyroxene–plagioclase–plagioclase dihedralangle, _cpp, in gabbroic cumulates records the time-integratedthermal history in the sub-solidus and provides a measure oftextural maturity. Variations in _cpp through the Layered Seriesof the Skaergaard intrusion, East Greenland, demonstrate thatthe onset of crystallization of clinopyroxene (within LZa),Fe–Ti oxides (at the base of LZc) and apatite (at thebase of UZb) as liquidus phases in the bulk magma is recordedby a stepwise increase in textural maturity, related to an increasein the contribution of latent heat to the total heat loss tothe surroundings and a reduction in the specific cooling rateat the crystallization front of the intrusion. The onset ofboth liquidus Fe–Ti oxide and apatite crystallizationis marked by a transient increase in textural maturity, probablylinked to overstepping before nucleation. Textural maturationat pyroxene–plagioclase–plagioclase triple junctionseffectively ceases in the uppermost parts of the Layered Seriesas a result of the entire pluton cooling below the closure temperaturefor dihedral angle change, which is 1075°C. Solidificationof the Layered Series of the Skaergaard intrusion occurred viathe upwards propagation of a mush zone only a few metres thick. KEY WORDS: magma; partial melting; asthenosphere; olivine; mantle     

The Geochemical Evolution of Anorthosite Residual Magmas in the Laramie Anorthosite Complex, Wyoming

Olivine- and pyroxene-bearing Fe-enriched dioritic rocks inthe 1434 Ma Laramie anorthosite complex are interpreted to representvariably fractionated and contaminated magmas residual afterthe crystallization of anorthosite. Geochemical characteristicsof this suite include the following: high contents of TiO₂,, and P₂O₅; high incompatibletrace element contents; rare earth element patterns with a largerange of Eu anomalies; and isotopic compositions that reflectthe geographic location of individual samples, with I_Sr increasingand _Nd decreasing from south to north. After extraction fromanorthosite, fractionation of ferrodioritic residual magmasresulted in secondary residual monzodioritic melts and complementaryoxide-rich ferrodiorite cumulates. Geographic trends in isotopiccomposition reflect an increasing Archean crustal componentfrom south to north. Dioritic dikes and cumulates with isotopiccompositions similar to associated anorthosites were derivedlocally. Large isotopic discrepancies between some dioritesand their hosting anorthosites reflect preferential contaminationof residual magma during ascent and emplacement of mantle-derivedplagioclase-rich diapirs, followed by subsequent extractionand isolation of Fe-enriched interstitial melt. Strong isotopiccontrasts between anorthosite and associated Fe-enriched rocksin anorthosite complexes do not preclude a direct relationshipbetween them and reflect the diversity and complexity of processesduring their petrogenesis. KEY WORDS: anorthosite; ferrodiorite; geochemistry; Laramie anorthosite complex; residual magma ^*Corresponding author: Present address: Ceramic Science and Technology, Mail Stop K762, Los Alamos National Laboratory, Los Alamos, NM 87545, USA. Phone: (505) 665-3934. Fax: (505) 665-3935. E-mail: jeremy{at}lanl.gov.     

A Strontium and Neodymium Isotopic Investigation of the Fongen--Hyllingen Layered Intrusion, Norway

The Fongen—Hyllingen Intrusion, situated 60 km SE of Trondheim,Norway, is a synorogenic layered mafic intrusion of Caledonianage . The intrusion is divided into four evolutionary stages based on cryptic variations: StageI—a basal reversal; Stage II—unchanged mineral chemistryor slight normal evolution; Stage III—a gradual regression;Stage IV— a strong normal fractionation trend Magma replenishmentdominated during most of the crystallization, i.e. during StagesI, II and III Replenishing magma was more dense than resident,evolved magma, and continuing influx eventually caused a compositionallystratified magma column to form. Cryptic lateral variation isan important feature in the southern part of the complex andformed by in situ crystallization from a stratified magma alongan inclined floor, where modal layering formed parallel to thecrystallization front. Initial Sr- and Nd-isotopic ratios inthe cumulates vary as a result of assimilation of country rockand subsequent mixing between uncontaminated, replenishing magmaand contaminated, resident magma. The parental magma had a moderatelydepleted isotope composition, relative to Bulk Earth, with _Nd=584and Sr_i=070308, whereas the main contaminant was a partialmelt of metapelitic country rock with _Nd=-874 and Sr_i=07195(Sr_i is the initial ⁸⁷Sr/⁸⁶Sr). Sr_i in the analysed cumulatewhole-rock samples ranges from 070308 to 070535 and initial_Nd ranges from. 158 to 584. There is a strong correlationbetween mineralogical composition and isotopic trends in mostof the cumulates: the most primitive samples are the least contaminated,as reflected by relatively high ed and low Sr,, and more evolvedsamples have progressively lower eNi and higher Sry A gradualregression of several hundred metres thickness characterizesStage III; stratigraphically upwards mineral compositions becomemore primitive and isotope compositions more depleted (higher_Nd and lower Sr_i), implying a process of. progressive mixing-inof replenishing, primitive and uncontaminated magma. Magma influxin Stage III took place by fountaining, whereas magma additionwas more tranquil in the earlier stages. The fountaining influxentrained resident, relatively evolved and contaminated magma,resulting in a hybrid magma which ponded at the floor. Duringprolonged magma addition with concomitant crystallization, thelowermost magma layer was replaced by progressively more primitivehybrid magma, creating a gradual regression in the crystallizingcumulate sequence. A detailed two-dimensional study revealslateral variations in mineral compositions both at the baseand top of Stage III, whereas lateral variations in Sr- andNd-isotopic compositions are present at the top, but not atthe base. This implies that the lowest crystallizing part ofthe magma column was essentially isotopically homogeneous, butcompositionally stratified, before influx in Stage III. Isotopicgradients in the magma were strong close to the roof, wheremost of the assimilation occurred, and decreased downwards,merging into isotopically homogeneous magma. This stratifiedsystem was destroyed by turbulent mixing between replenishingand resident magma during fountaining influx in Stage III, anda new stratification was established with both an isotopic anda compositional gradient. After the final influx, crystallizationcontinued in an essentially closed system, in which the remainingmagma column eventually became homogenized, as magma layersmixed when their densities converged owing to release of buoyant,residual liquid during fractional crystallization. Corresponding author     

Diverse Origins of Xenoliths from Seamounts at the Continental Margin, Offshore Central California

A diverse assemblage of small mafic and ultramafic xenolithsoccurs in alkalic lava from Davidson and Pioneer seamounts locatedat the continental margin of central California. Based on mineralcompositions and textures, they form three groups: (1) mantlexenoliths of lherzolite, pyroxenite, and dunite with olivineof >Fo₉₀; (2) ocean crust xenoliths of dunite with olivine<Fo₉₀, troctolite, pyroxene-gabbro, and anorthosite withlow-K₂O plagioclase; (3) cumulates of seamount magmas of alkalicgabbro with primary amphibole and biotite and anorthosites withhigh-K₂O plagioclase. The alkalic cumulates are geneticallyrelated to, but more evolved than, their host lavas and probablycrystallized at the margins of magma reservoirs. Modeling andcomparison with experimentally derived phases suggest an originat moderate pressures (0·5–0·9 GPa). Thehigh volatile contents of the alkalic host lavas may have pressurizedthe magma chambers and helped to propel the xenolith-bearinglavas directly from deep storage at the base of the lithosphereto the eruption site on the ocean floor, entraining fragmentsof the upper mantle and ocean crust cumulates from the underlyingabandoned spreading center. KEY WORDS: basaltic magmatism; continental margin seamounts; geothermobarometry; mineral chemistry; xenoliths     

Rates of Thermal and Chemical Evolution of Magmas in a Cooling Magma Chamber: a Chronological and Theoretical Study on Basaltic and Andesitic Lavas from Rishiri Volcano, Japan

Rates of magmatic processes in a cooling magma chamber wereinvestigated for alkali basalt and trachytic andesite lavaserupted sequentially from Rishiri Volcano, northern Japan, bydating of these lavas using ²³⁸U–²³⁰Th radioactive disequilibriumand ¹⁴C dating methods, in combination with theoretical analyses.We obtained the eruption age of the basaltic lavas to be 29·3± 0·6 ka by ¹⁴C dating of charcoals. The eruptionage of the andesitic lavas was estimated to be 20·2 ±3·1 ka, utilizing a whole-rock isochron formed by U–Thfractionation as a result of degassing after lava emplacement.Because these two lavas represent a series of magmas producedby assimilation and fractional crystallization in the same magmachamber, the difference of the ages (i.e. 9 kyr) is a timescaleof magmatic evolution. The thermal and chemical evolution ofthe Rishiri magma chamber was modeled using mass and energybalance constraints, as well as quantitative information obtainedfrom petrological and geochemical observations on the lavas.Using the timescale of 9 kyr, the thickness of the magma chamberis estimated to have been about 1·7 km. The model calculationsshow that, in the early stage of the evolution, the magma cooledat a relatively high rate (>0·1°C/year), and thecooling rate decreased with time. Convective heat flux fromthe main magma body exceeded 2 W/m² when the magma was basaltic,and the intensity diminished exponentially with magmatic evolution.Volume flux of crustal materials to the magma chamber and rateof convective melt exchange (compositional convection) betweenthe main magma and mush melt also decreased with time, from 0·1 m/year to 10^–3 m/year, and from 1 m/yearto 10^–2 m/year, respectively, as the magmas evolved frombasaltic to andesitic compositions. Although the mechanism ofthe cooling (i.e. thermal convection and/or compositional convection)of the main magma could not be constrained uniquely by the model,it is suggested that compositional convection was not effectivein cooling the main magma, and the magma chamber is consideredto have been cooled by thermal convection, in addition to heatconduction. KEY WORDS: convection; magma chamber; heat and mass transport; timescale; U-series disequilibria     

Geochemistry of Picritic and Associated Basalt Flows of the Western Emeishan Flood Basalt Province, China

Picritic lava flows near Lijiang in the late Permian Emeishanflood basalt province are associated with augite-phyric basalt,aphyric basalt, and basaltic pyroclastic units. The dominantphenocryst in the picritic flows is Mg-rich olivine (up to 91·6%forsterite component) with high CaO contents (to 0·42wt %) and glass inclusions, indicating that the olivine crystallizedfrom a melt. Associated chromite has a high Cr-number (73–75).The estimated MgO content of the primitive picritic liquidsis about 22 wt %, and initial melt temperature may have beenas high as 1630–1690°C. The basaltic lavas appearto be related to the picritic ones principally by olivine andclinopyroxene fractionation. Age-corrected Nd–Sr–Pbisotope ratios of the picritic and basaltic lavas are indistinguishableand cover a relatively small range [e.g. _Nd(t) = –1·3to +4·0]. The higher _Nd(t) lavas are isotopically similarto those of several modern oceanic hotspots, and have ocean-island-likepatterns of alteration-resistant incompatible elements. Heavyrare earth element characteristics indicate an important rolefor garnet during melting and that the lavas were formed byfairly small degrees of partial melting. Rough correlationsof isotope ratios with ratios of alteration-resistant highlyincompatible elements (e.g. Nb/La) suggest modest amounts ofcontamination involving continental material or a relativelylow-_Nd component in the source. Overall, our results are consistentwith other evidence suggesting some type of plume-head originfor the Emeishan province. KEY WORDS: Emeishan; flood basalts; picrites; mantle plumes; late Permian     

The Geochemistry of Ultramafic to Mafic Volcanics from the Belingwe Greenstone Belt, Zimbabwe: Magmatism in an Archean Continental Large Igneous Province

The evolution of the late Archean Belingwe greenstone belt,Zimbabwe, is discussed in relation to the geochemistry of theultramafic to mafic volcanic rocks. Four volcanic types (komatiite,komatiitic basalt, D-basalt and E-basalt) are distinguishedin the 2·7 Ga Ngezi volcanic sequence using a combinationof petrography and geochemistry. The komatiites and D-basaltsare rocks in which isotopic systems and trace elements are depleted.Chemical variations in komatiites and D-basalts can be explainedby fractional crystallization from the parental komatiite. Incontrast, komatiitic basalts and E-basalts are siliceous anddisplay enriched isotopic and trace element compositions. Theirchemical trends are best explained by assimilation with fractionalcrystallization (AFC) from the primary komatiite. AFC calculationsindicate that the komatiitic basalts and E-basalts are derivedfrom komatiites contaminated with 20% and 30% crustal material,respectively. The volcanic stratigraphy of the Ngezi sequence,which is based on field relationships and the trace elementcompositions of relict clinopyroxenes, shows that the leastcontaminated komatiite lies between highly contaminated komatiiticbasalt flows, and has limited exposure near the base of thesuccession. Above these flows, D- and E-basalts alternate. Thekomatiite appears to have erupted on the surface only in theearly stages, when plume activity was high. As activity decreasedwith time, komatiite magmas may have stagnated to form magmachambers within the continental crust. Subsequent komatiiticmagmas underwent fractional crystallization and were contaminatedwith crust to form D-basalts or E-basalts. KEY WORDS: komatiite; crustal assimilation; Belingwe greenstone belt; continental flood basalt; plume magmatism     

Lower Crustal Magma Genesis and Preservation: a Stochastic Framework for the Evaluation of Basalt-Crust Interaction

We present a quantitative assessment of the thermal and dynamicresponse of an amphibolitic lower crust to the intrusion ofbasaltic dike swarms in an arc setting. We consider the effectof variable intrusion geometry, depth of intrusion, and basaltflux on the production, persistence, and interaction of basalticand crustal melt in a stochastic computational framework. Distinctmelting and mixing environments are predicted as a result ofthe crustal thickness and age of the arc system. Shallow crustal(30 km) environments and arc settings with low fluxes of mantle-derivedbasalt are likely repositories of isolated pods of mantle andcrustal melts in the lower crust, both converging on daciticto rhyodacitic composition. These may be preferentially rejuvenatedin subsequent intrusive episodes. Mature arc systems with thickercrust (50 km) produce higher crustal and residual basaltic meltfractions, reaching 0·4 for geologically reasonable basaltfluxes. The basaltic to basaltic andesite composition of bothcrustal and mantle melts will facilitate mixing as the networkof dikes collapses, and Reynolds numbers reach 10^–4–1·0in the interiors of dikes that have been breached by ascendingcrustal melts. This may provide one mechanism for melting, assimilation,storage and homogenization (MASH)-like processes. Residual mineralassemblages of crust thickened by repeated intrusion are predictedto be garnet pyroxenitic, which are denser than mantle peridotiteand also generate convective instabilities where some of thecrustal material is lost to the mantle. This reconciles thethinner than predicted crust in regions that have undergonea large flux of mantle basalt for a prolonged period of time,and helps explain the enrichment of incompatible elements suchas K₂O, typical of mature arc settings, without the associatedmass balance problem. KEY WORDS: crustal anatexis; delamination; lower crust; magma mixing; thermal model     

Contact Metamorphosed Ultramafic Rocks in the Western Sierra Nevada Foothills, California

Within the western Sierra Nevada metamorphic belt, linear bodiesof alpine-type ultramafic rock, now composed largely of serpentineminerals, parallel the regional strike and commonly coincidewith major fault zones. Within this metamorphic belt, east ofSacramento, California, ultramafic rocks near a large maficintrusion, the Pine Hill Intrusive Complex, have been emplacedduring at least two separate episodes. Those ultramafic rocks,evidently unaffected by the Pine Hill Intrusive Complex andcomposed largely of serpentine minerals, were emplaced alonga major fault zone after emplacement of the Pine Hill IntrusiveComplex. Those ultramafic rocks, contact metamorphosed by thePine Hill Intrusive Complex, show a zonation of mineral assemblagesas the igneous contact is approached: olivine+antigorite+chlorite+tremolite+Fe-Cr spinel olivine+talc+chlorite+tremolite+Fe-Crspinel olivine+anthophyllite+chlorite+tremolite+Fe-Cr spinel olivine+orthopyroxene+aluminous spinel+hornblende+Fe-Cr spinel.Superimposed on these mineral assemblages are abundant secondaryminerals (serpentine minerals, talc, chlorite, magnetite) whichformed after contact metamorphism. Correlation of observed mineralassemblages with the experimental systems, MgO-SiO₂-H₂O andMgO-Al₂O₃-SiO₂-H₂O suggests an initial contact temperature of775±25 °C for the Pine Hill Intrusive Complex assumingPtotal P_fluid PH₂O. The pressure acting on the metamorphic rockduring emplacement of the intrusion is estimated to be a minimumof 1.5 kb.     

Open-System Magma Chamber Evolution: an Energy-constrained Geochemical Model Incorporating the Effects of Concurrent Eruption, Recharge, Variable Assimilation and Fractional Crystallization (EC-E'RA{chi}FC)

Significant petrogenetic processes governing the geochemicalevolution of magma bodies include magma Recharge (includingformation of ‘quenched inclusions’ or enclaves),heating and concomitant partial melting of country rock withpossible ‘contamination’ of the evolving magma body(Assimilation), and formation and separation of cumulates byFractional Crystallization (RAFC). Although the importance ofmodeling such open-system magma chambers subject to energy conservationhas been demonstrated, the effects of concurrent removal ofmagma by eruption and/or variable assimilation (involving imperfectextraction of anatectic melt from wall rock) have not been considered.In this study, we extend the EC-RAFC model to include the effectsof Eruption and variable amounts of assimilation, A. This model,called EC-E'RAFC, tracks the compositions (trace elements andisotopes), temperatures, and masses of magma body liquid (melt),eruptive magma, cumulates and enclaves within a composite magmaticsystem undergoing simultaneous eruption, recharge, assimilationand fractional crystallization. The model is formulated as aset of 4 + t + i + s coupled nonlinear differential equations,where the number of trace elements, radiogenic and stable isotoperatios modeled are t, i and s, respectively. Solution of theEC-E'RAFC equations provides values for the average temperatureof wall rock (T_a), mass of melt within the magma body (M_m),masses of cumulates (M_ct), enclaves (M_en) and wall rock () and the masses of anatectic melt generated () and assimilated (). In addition, t trace element concentrations and i + s isotopic ratios inmelt and eruptive magma (C_m, _m, _m), cumulates (C_ct, _m, _m), enclaves(C_en, , ) and anatectic melt (C_a, , ) as a function of magma temperature (T_m) are also computed. Input parametersinclude the (user-defined) equilibration temperature (T_eq),a factor describing the efficiency of addition of anatecticmelt () from country rock to host magma, the initial temperatureand composition of pristine host melt (, , , ), recharge melt (, , , ) and wall rock (, , , ), distribution coefficients (D_m, D_r, D_a) and their temperaturedependences (H_m, H_r, H_a), latent heats of transition (meltingor crystallization) for wall rock (h_a), pristine magma (h_m)and recharge magma (h_r) as well as the isobaric specific heatcapacity of assimilant (C_p,a), pristine (C_p,m) and recharge(C_p,r) melts. The magma recharge mass and eruptive magma massfunctions, M_r(T_m) and M_e(T_m), respectively, are specified apriori. M_r(T_m) and M_e(T_m) are modeled as either continuous orepisodic (step-like) processes. Melt productivity functions,which prescribe the relationship between melt mass fractionand temperature, are defined for end-member bulk compositionscharacterizing the local geologic site. EC-E'RAFC has potentialfor addressing fundamental questions in igneous petrology suchas: What are intrusive to extrusive ratios (I/E) for particularmagmatic systems, and how does this factor relate to rates ofcrustal growth? How does I/E vary temporally at single, long-livedmagmatic centers? What system characteristics are most profoundlyinfluenced by eruption? What is the quantitative relationshipbetween recharge and assimilation? In cases where the extractionefficiency can be shown to be less than unity, what geologiccriteria are important and can these criteria be linked to fieldobservations? A critical aspect of the energy-constrained approachis that it requires integration of field, geochronological,petrologic, and geochemical data, and, thus, the EC-ERAFC ‘systems’approach provides a means for answering broad questions whileunifying observations from a number of disciplines relevantto the study of igneous rocks. KEY WORDS: assimilation; energy conservation; eruption; open system; recharge     

Garnet-bearing Xenoliths from Salt Lake Crater, Oahu, Hawaii: High-Pressure Fractional Crystallization in the Oceanic Mantle

The focus of this study is a suite of garnet-bearing mantlexenoliths from Oahu, Hawaii. Clinopyroxene, olivine, and garnetconstitute the bulk of the xenoliths, and orthopyroxene is presentin small amounts. Clinopyroxene has exsolved orthopyroxene,spinel, and garnet. Many xenoliths also contain spinel-coredgarnets. Olivine, clinopyroxene, and garnet are in major elementchemical equilibrium with each other; large, discrete orthopyroxenedoes not appear to be in major-element chemical equilibriumwith the other minerals. Multiple compositions of orthopyroxeneoccur in individual xenoliths. The new data do not support theexisting hypothesis that all the xenoliths formed at 1 6–22GPa, and that the spinel-cored garnets formed as a consequenceof almost isobaric subsolidus cooling of a spinel-bearing assemblage.The lack of olivine or pyroxenes in the spinel–garnetreaction zones and the embayed outline of spinel grains insidegarnet suggest that the spinel-cored garnets grew in the presenceof a melt. The origin of these xenoliths is interpreted on thebasis of liquidus phase relations in the tholeiitic and slightlysilica-poor portion of the CaO–MgO–Al₂O₃–SiO₂(CMAS) system at pressures from 30 to 50 GPa. The phase relationssuggest crystallization from slightly silica-poor melts (ortransitional basaltic melts) in the depth range 110–150km beneath Oahu. This depth estimate puts the formation of thesexenoliths in the asthenosphere. On the basis of this study itis proposed that the pyroxenite xenoliths are high-pressurecumulates related to polybaric magma fractionation in the asthenosphere,thus making Oahu the only locality among the oceanic regionswhere such deep magmatic fractional crystallization processeshave been recognized. KEY WORDS: xenolith; asthenosphere; basalt; CMAS; cumulate; oceanic lithosphere; experimental petrology; mantle; geothermobarometry; magma chamber     

Petrogenesis of the Fiambal Gabbroic Intrusion, Northwestern Argentina, a Deep Crustal Syntectonic Pluton in a Continental Magmatic Arc

The Fiambala gabbroic intrusion is an Ordovician pluton, 30km² in outcrop area, situated in the Sierras Pampeanas of NWArgentina. It was emplaced at a depth of 21–24 km in anactively deforming continental magmatic are, and it representssome of the earliest are magmatism in the southern Cordillera.Its formation was accompanied by the emplacement of voluminousgabbroic sills and dikes, and was synchronous with regionalamphibolite facies metamorphism that is locally upgraded togranulite facies along the edges of the intrusion. The intrusionhas been tilted to the east and now stands almost verticallyon edge, striking NNW-SSE. It has a 200–500-m basal cumulatezone of dunite and websterite along its western margin, overlainstratigraphically to the east by 2.5 km of poorly layered two-pyroxenegabbro and homblende gabbro. Although olivine is present inthe basal ultramafic cumulates, nowhere does it coexist withplagioclase. The common differentiation indices, such as Ancontent of plagioclase, 100 Mg/(Mg+Fe) in the pyroxenes andthe whole rocks, and the concentrations of various incompatibletrace elements, all vary systematically through the intrusion.The parental magma was evidently a hydrous, high-Mg tholeiitewith >11 wt.% MgO, 100 Mg/(Mg+Fe)=73–74, and REE abundancesapproximately 10 chondrites. Its estimated trace element characteristicswere strikingly similar to those of primary island are magmasfrom occanic regions. Most of the gabbroic rocks represent orapproach liquid compositions; they are hypersthene-normative,almost silica-saturated tholeiites. Evolution from the mostmafic compositions comparable to primitive island are basaltsto derivative compositions with the typical characteristicsof continental are basalts occurred through the combined effectsof fractional crystallization, multiple magma batch injection,and crustal assimilation. There appears to be no need to invokea special subcontinental mantle source to produce the traceelement signatures that are so characteristic of continentalare magmas. These signatures can be developed by limited crystallizationand assimilation in the lower crust.     

Early-Middle Jurassic Dolerite Dykes from Western Dronning Maud Land (Antarctica): Identifying Mantle Sources in the Karoo Large Igneous Province

A suite of dolerite dykes from the Ahlmannryggen region of westernDronning Maud Land (Antarctica) forms part of the much moreextensive Karoo igneous province of southern Africa. The dykecompositions include both low- and high-Ti magma types, includingpicrites and ferropicrites. New ⁴⁰Ar/³⁹Ar age determinationsfor the Ahlmannryggen intrusions indicate two ages of emplacementat 178 and 190 Ma. Four geochemical groups of dykes have beenidentified in the Ahlmannryggen region based on analyses of60 dykes. The groups are defined on the basis of whole-rockTiO₂ and Zr contents, and reinforced by rare earth element (REE),⁸⁷Sr/⁸⁶Sr and ¹⁴³Nd/¹⁴⁴Nd isotope data. Group 1 were intrudedat 190 Ma and have low TiO₂ and Zr contents and a significantArchaean crustal component, but also evidence of hydrothermalalteration. Group 2 dykes were intruded at 178 Ma; they havelow to moderate TiO₂ and Zr contents and are interpreted tobe the result of mixing of melts derived from an isotopicallydepleted source with small melt fractions of an enriched lithosphericmantle source. Group 3 dyke were intruded at 190 Ma and formthe most distinct magma group; these are largely picritic withsuperficially mid-ocean ridge basalt (MORB)-like chemistry (flatREE patterns, ⁸⁷Sr/⁸⁶Sr_i 0·7035, Nd_i 9). However, theyhave very high TiO₂ (4 wt %) and Zr (500 ppm) contents, whichis not consistent with melting of MORB-source mantle. The Group3 magmas are inferred to be derived by partial melting of astrongly depleted mantle source in the garnet stability field.This group includes several high Mg–Fe dykes (ferropicrites),which are interpreted as high-temperature melts. Some Group3 dykes also show evidence of contamination by continental crust.Group 4 dykes are low-K picrites intruded at 178 Ma; they havevery high TiO₂–Zr contents and are the most enriched magmagroup of the Karoo–Antarctic province, with ocean-islandbasalt (OIB)-like chemistry. Dykes of Group 1 and Group 3 aresub-parallel (ENE–WSW) and both groups were emplaced at190 Ma in response to the same regional stress field, whichhad changed by 178 Ma, when Group 2 and Group 4 dykes were intrudedalong a dominantly NNE–SSW strike. KEY WORDS: flood basalt; depleted mantle; enriched mantle; Ahlmannryggen; Karoo dyke     

Phase Equilibria Constraints on the Chemical and Physical Evolution of the Campanian Ignimbrite

The Campanian Ignimbrite is a > 200 km³ trachyte–phonolitepyroclastic deposit that erupted at 39·3 ± 0·1ka within the Campi Flegrei west of Naples, Italy. Here we testthe hypothesis that Campanian Ignimbrite magma was derived byisobaric crystal fractionation of a parental basaltic trachyandesiticmelt that reacted and came into local equilibrium with smallamounts (5–10 wt%) of crustal rock (skarns and foid-syenites)during crystallization. Comparison of observed crystal and magmacompositions with results of phase equilibria assimilation–fractionationsimulations (MELTS) is generally very good. Oxygen fugacitywas approximately buffered along QFM + 1 (where QFM is the quartz–fayalite–magnetitebuffer) during isobaric fractionation at 0·15 GPa ( 6km depth). The parental melt, reconstructed from melt inclusionand host clinopyroxene compositions, is found to be basaltictrachyandesite liquid (51·1 wt% SiO₂, 9·3 wt%MgO, 3 wt% H₂O). A significant feature of phase equilibria simulationsis the existence of a pseudo-invariant temperature, 883 °C,at which the fraction of melt remaining in the system decreasesabruptly from 0·5 to < 0·1. Crystallizationat the pseudo-invariant point leads to abrupt changes in thecomposition, properties (density, dissolved water content),and physical state (viscosity, volume fraction fluid) of meltand magma. A dramatic decrease in melt viscosity (from 1700Pa s to 200 Pa s), coupled with a change in the volume fractionof water in magma (from 0·1 to 0·8) and a dramaticdecrease in melt and magma density acted as a destabilizingeruption trigger. Thermal models suggest a timescale of 200kyr from the beginning of fractionation until eruption, leadingto an apparent rate of evolved magma generation of about 10^–3km³/year. In situ crystallization and crystal settling in density-stratifiedregions, as well as in convectively mixed, less evolved subjacentmagma, operate rapidly enough to match this apparent volumetricrate of evolved magma production. KEY WORDS: assimilation; Campanian Ignimbrite; fractional crystallization; magma dynamics; phase equilibria