He–Ar and Nd–Sr isotopic compositions of late Pleistocene felsic plutonic back arc basin rocks from Ulleungdo volcanic island, South Korea: Implications for the genesis of young plutonic rocks in a back arc basin

We report analyses of noble gases and Nd–Sr isotopes in mineral separates and whole rocks of late Pleistocene (< 0.2 Ma) monzonites from Ulleungdo, South Korea, a volcanic island within the back arc basin of the Japan island arc. A Rb–Sr mineral isochron age for the monzonites is 0.12 ± 0.01 Ma. K–Ar biotite ages from the same samples gave relatively concordant ages of 0.19 ± 0.01and 0.22 ± 0.01 Ma. ⁴⁰Ar/³⁹Ar yields a similar age of 0.29 ± 0.09 Ma. Geochemical characteristics of the felsic plutonic rocks, which are silica oversaturated alkali felsic rocks (av., 12.5 wt% in K₂O + Na₂O), are similar to those of 30 alkali volcanics from Ulleungdo in terms of concentrations of major, trace and REE elements. The initial Nd–Sr isotopic ratios of the monzonites (⁸⁷Sr/⁸⁶Sr = 0.70454–0.71264, ¹⁴³Nd/¹⁴⁴Nd = 0.512528–0.512577) are comparable with those of the alkali volcanics (⁸⁷Sr/⁸⁶Sr = 0.70466–0.70892, ¹⁴³Nd/¹⁴⁴Nd = 0.512521–0.512615) erupted in Stage 3 of Ulleungdo volcanism (0.24–0.47 Ma). The high initial ⁸⁷Sr/⁸⁶Sr values of the monzonites imply that seawater and crustally contaminated pre-existing trachytes may have been melted or assimilated during differentiation of the alkali basaltic magma.A mantle helium component (³He/⁴He ratio of up to 6.5 RA) associated with excess argon was found in the monzonites. Feldspar and biotite have preferentially lost helium during slow cooling at depth and/or during their transportation to the surface in a hot host magma. The source magma noble gas isotopic features are well preserved in fluid inclusions in hornblende, and indicate that the magma may be directly derived from subcontinental lithospheric mantle metasomatized by an ancient subduction process, or may have formed as a mixture of MORB-like mantle and crustal components. The radiometric ages, geochemical and Nd–Sr isotopic signatures of the Ulleungdo monzonites as well as the presence of mantle-derived helium and argon, suggests that these felsic plutonic rocks evolved from alkali basaltic magma that formed by partial melting of subcontinental lithospheric mantle beneath the back arc basin located along the active continental margin of the southeastern part of the Eurasian plate.     

Experimental solidification of anhydrous latitic and trachytic melts at different cooling rates: The role of nucleation kinetics

Two sets of cooling experiments were run at atmospheric conditions for two anhydrous starting latitic and trachytic melts: 1) five cooling rates (25, 12.5, 3, 0.5, and 0.125 °C/min) between 1300° and 800 °C, and 2) a 0.5 °C/min cooling rate from 1300 °C with quench temperatures at 1200°, 1100°, 1000° and 900 °C. Trachytic run-products are invariably glassy. Nucleation is also suppressed in the latitic run-products at the three highest cooling rates. Conversely, in the 0.5 and 0.125 °C/min runs, latites have a crystal content of  90 vol.%. The phases are: plagioclase, clinopyroxene, glass and iron-bearing oxide (in order of abundance). The variable quench temperatures, investigated by coupling experiments with Pt wire and Pt capsule sample containers in set 2, again did not produce crystallization of trachyte, whereas latitic samples are characterized by  10 vol.% of oxides, pyroxenes and plagioclase (in order of appearance), at temperature < 1000 °C. Effects of (preferential) heterogeneous nucleation on sample holders, of superheating degree, and chemical species loss during cooling are absent for both melt compositions. The difference of solidification paths between these two silicate melts can be ascribed only to their small chemical differences. In comparison with calculated equilibrium conditions all the experimental latitic and trachytic run-products revealed strong kinetic effects, interpretable in the light of the nucleation theory. The glass-forming ability (GFA) of trachyte is higher, whereas their critical cooling rate (Rc) is lower (< 0.125 °C/min), in comparison to latitic melts (Rc > 0.5 °C/min). The experimental results carried out in this study can be applied to lava flows and domes; trachytic lavas are able to flow for longer period with respect to latitic ones in a metastable condition. Glass-rich terrestrial lavas, i.e. obsidians, can be the result of sluggish nucleation kinetics due to the relative high polymerisation of evolved silicate melts.     

Development of a 3D GIS and its application to karst areas

There is a growing interest in modeling and analyzing karst phenomena in three dimensions. This paper integrates geology, groundwater hydrology, geographic information system (GIS), database management system (DBMS), visualization and data mining to study karst features in Huaibei, China. The 3D geo-objects retrieved from the karst area are analyzed and mapped into different abstract levels. The spatial relationships among the objects are constructed by a dual-linker. The shapes of the 3D objects and the topological models with attributes are stored and maintained in the DBMS. Spatial analysis was then used to integrate the data in the DBMS and the 3D model to form a virtual reality (VR) to provide analytical functions such as distribution analysis, correlation query, and probability assessment. The research successfully implements 3D modeling and analyses in the karst area, and meanwhile provides an efficient tool for government policy-makers to set out restrictions on water resource development in the area.     

Groundwater flow and contaminant transport modelling at an air weapons range

Numerical modelling was done at the Cold Lake Air Weapons Range, Canada, to test whether the dissolved RDX and nitrate detected in groundwater come from the same sources, and to predict whether contamination poses a threat to the surface water receptors near the site. Military live fire training activities may indeed pose a risk of contamination to groundwater resources, however field investigations on military bases are quite recent, and little information is available on the long-term behaviour of munition residues related contaminants. Very limited information was available about the contaminant source zones, which were assigned based on our knowledge of current training activities. The RDX plume was well represented with the model, but the heterogeneous distribution of nitrate concentrations was more difficult to reproduce. It was nonetheless determined that both contaminants originate from the same areas. According to the model, both contaminants should reach the nearby river, but concentrations in the river should remain very low if the source zone concentration does not change. Finally, the model allowed the recommendation of a new location for the main bombing target, which would offer added protection to the river and the lake into which it flows.     

Magma Ascent and Crustal Accretion at Ultraslow-Spreading Ridges: Constraints from Plagioclase Ultraphyric Basalts from the Arctic Mid-Ocean Ridge

Plagioclase ultraphyric basalts (PUBs) with up to 54% plagioclasephenocrysts were dredged in the rift valley and adjacent flanksof the ultraslow-spreading Mohns and Knipovich ridges. The PUBsshow large variations in crystal morphologies and zoning. Thelarge variations suggest that single basalt samples containa mixture of plagioclase crystals that aggregated at differentlevels in the magma conduits. Resorbed crystals and repeatedreverse zones suggest that the magma reservoirs were replenishedand heated several times. Thin concentric zones with melt inclusions,and sharp reductions in the anorthite content of 3–7%,are common between the reverse zones. These zones, and skeletalcrystals with distinctly lower anorthite contents than massivecrystals, are interpreted to be the result of rapid crystalliztionduring strong undercooling. The changes between short periodsof cooling and longer periods with reheating are explained bymultiple advances of crystal-rich magma into cool regions followedby longer periods of gradual magma inflow and temperature increase.The porphyritic basalts are characterizd by more depleted andmore fractionated compositions than the aphyric basalts, withlower (La/Sm)_N, K₂O and Mg-numbers. This relationship, and theobservation that PUBs are sampled only close to segment centresalong these ridges, suggests that the PUBs formed by higherdegrees of melting and evolved in more long-lived magma reservoirs.We propose that the zoning patterns of plagioclase crystalsand crystal morphologies of these PUBs reflect the developmentand flow of magma through a stacked sill complex-like conduitsystem, whereas the aphyric equivalents represent later flowof magma through the conduit. The formation of voluminous higher-degreemelts may trigger the development of the magma conduits andexplain the generally depleted compositions of PUB magmas. KEY WORDS: basalt; mineral chemistry; MORB; magma mixing; magma chamber; major element     

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     

Global Correlations of Ocean Ridge Basalt Chemistry with Axial Depth: a New Perspective

The petrological parameters Na₈ and Fe₈, which are Na₂O andFeO contents in mid-ocean ridge basalt (MORB) melts correctedfor fractionation effects to MgO = 8 wt%, have been widely usedas indicators of the extent and pressure of mantle melting beneathocean ridges. We find that these parameters are unreliable.Fe₈ is used to compute the mantle solidus depth (P_o) and temperature(T_o), and it is the values and range of Fe₈ that have led tothe notion that mantle potential temperature variation of T_P= 250 K is required to explain the global ocean ridge systematics.This interpreted T_P = 250 K range applies to ocean ridges awayfrom ‘hotspots’. We find no convincing evidencethat calculated values for P_o, T_o, and T_P using Fe₈ have anysignificance. We correct for fractionation effect to Mg^# = 0·72,which reveals mostly signals of mantle processes because meltswith Mg^# = 0·72 are in equilibrium with mantle olivineof Fo_89·6 (vs evolved olivine of Fo_{88·1–79·6}in equilibrium with melts of Fe₈). To reveal first-order MORBchemical systematics as a function of ridge axial depth, weaverage out possible effects of spreading rate variation, local-scalemantle source heterogeneity, melting region geometry variation,and dynamic topography on regional and segment scales by usingactual sample depths, regardless of geographical location, withineach of 22 ridge depth intervals of 250 m on a global scale.These depth-interval averages give Fe₇₂ = 7·5–8·5,which would give T_P = 41 K (vs 250 K based on Fe₈) beneathglobal ocean ridges. The lack of Fe₇₂–Si₇₂ and Si₇₂–ridgedepth correlations provides no evidence that MORB melts preservepressure signatures as a function of ridge axial depth. We thusfind no convincing evidence for T_P > 50 K beneath globalocean ridges. The averages have also revealed significantcorrelations of MORB chemistry (e.g. Ti₇₂, Al₇₂, Fe₇₂,Mg₇₂, Ca₇₂, Na₇₂ and Ca₇₂/Al₇₂) with ridge axial depth. Thechemistry–depth correlation points to an intrinsic linkbetween the two. That is, the 5 km global ridge axial reliefand MORB chemistry both result from a common cause: subsolidusmantle compositional variation (vs T_P), which determines themineralogy, lithology and density variations that (1) isostaticallycompensate the 5 km ocean ridge relief and (2) determine thefirst-order MORB compositional variation on a global scale.A progressively more enriched (or less depleted) fertileperidotite source (i.e. high Al₂O₃ and Na₂O, and low CaO/Al₂O₃)beneath deep ridges ensures a greater amount of modal garnet(high Al₂O₃) and higher jadeite/diopside ratios in clinopyroxene(high Na₂O and Al₂O₃, and lower CaO), making a denser mantle,and thus deeper ridges. The dense fertile mantle beneath deepridges retards the rate and restricts the amplitude of the upwelling,reduces the rate and extent of decompression melting, givesway to conductive cooling to a deep level, forces melting tostop at such a deep level, leads to a short melting column,and thus produces less melt and probably a thin magmatic crustrelative to the less dense (more refractory) fertile mantlebeneath shallow ridges. Compositions of primitive MORB meltsresult from the combination of two different, but geneticallyrelated processes: (1) mantle source inheritance and (2) meltingprocess enhancement. The subsolidus mantle compositional variationneeded to explain MORB chemistry and ridge axial depth variationrequires a deep isostatic compensation depth, probably in thetransition zone. Therefore, although ocean ridges are of shalloworigin, their working is largely controlled by deep processesas well as the effect of plate spreading rate variation at shallowlevels. KEY WORDS: mid-ocean ridges; mantle melting; magma differentiation; petrogenesis; MORB chemistry variation; ridge depth variation; global correlations; mantle compositional variation; mantle source density variation; mantle potential temperature variation; isostatic compensation     

Crystal-Melt Separation and the Development of Isotopic Heterogeneities in Hybrid Magmas

If a magma is a hybrid of two (or more) isotopically distinctend-members, at least one of which is partially crystalline,separation of melt and crystals after hybridization will leadto the development of isotopic heterogeneities in the magmaas long as some of the pre-existing crystalline material (antecrysts)retains any of its original isotopic composition. This holdstrue whether the hybridization event is magma mixing as traditionallyconstrued, bulk assimilation, or melt assimilation. Once a magma-scaleisotopic heterogeneity is formed by crystal–melt separation,it is essentially permanent, persisting regardless of subsequentcrystallization, mixing, or equilibration events. The magnitudeof the isotopic variability resulting from crystal–meltseparation can be as large as that resulting from differentialcontamination, multiple isotopically distinct sources, or insitu isotopic evolution. In one model, a redistribution of one-thirdof the antecryst cargo yielded a crystal-enriched sample with⁸⁷Sr/⁸⁶Sr of 0·7058, whereas the complementary crystal-poorsample has ⁸⁷Sr/⁸⁶Sr of 0·7068. In other models, crystal-richsamples are enriched in radiogenic Sr. Isotopic heterogeneitiescan be either continuous (controlled by the modal distributionof crystals and melt) or discontinuous (when there is completeseparation of crystals and liquid). The first case may be exemplifiedby some isotopically zoned large-volume rhyolites, formed bythe eruptive inversion of a modally zoned magma chamber. Inthe latter case, the isotopic composition of any (for example)interstitial liquid will be distinct from the isotopic compositionof the bulk crystal fraction. The separation of such an interstitialliquid may explain the presence of isotopically distinct late-stageaplites in plutons. Crystal–melt separation provides anadditional option for the interpretation of isotopically zonedor heterogeneous magmas. This option is particularly attractivefor systems whose chemical variation is otherwise explicableby fractionation-dominated processes. Non-isotopic chemicalheterogeneities can also develop in this fashion. KEY WORDS: isotopic heterogeneity; zoning; hybrid magma; crystal separation; Sr isotopes; aplite; rhyolite     

Sedimentary petrology and geochemistry of siliciclastic rocks from the upper Jurassic Tordillo Formation (Neuquén Basin, western Argentina): Implications for provenance and tectonic setting

The Upper Jurassic Tordillo Formation is exposed along the western edge of the Neuquén Basin (west central Argentina) and consists of fluvial strata deposited under arid/semiarid conditions. The pebble composition of conglomerates, mineralogical composition of sandstones and pelitic rocks, and major- and trace-element geochemistry of sandstones, mudstones, and primary pyroclastic deposits are evaluated to determine the provenance and tectonic setting of the sedimentary basin. Conglomerates and sandstones derived almost exclusively from volcanic sources. The stratigraphic sections to the south show a clast population of conglomerates dominated by silicic volcanic fragments and a predominance of feldspathic litharenites. This framework composition records erosion of Triassic–Jurassic synrift volcaniclastic rocks and basement rocks from the Huincul arch, which was exhumed as a result of Late Jurassic inversion. In the northwestern part of the study area, conglomerates show a large proportion of mafic and acidic volcanic rock fragments, and sandstones are characterised by a high content of mafic volcanic rock fragments and plagioclase. These data suggest that the source of the sandstones and conglomerates was primarily the Andean magmatic arc, located west of the Neuquén Basin. The clay mineral assemblage is interpreted as the result of a complex set of factors, including source rock, climate, transport, and diagenesis. Postdepositional processes produced significant variations in the original compositions, especially the fine-grained deposits. The Tordillo sediments are characterised by moderate SiO₂ contents, variable abundances of K₂O and Na₂O, and a relatively high proportion of ferromagnesian elements. The degree of chemical weathering in the source area, expressed as the chemical index of alteration, is low to moderate. The major element geochemistry and Th/Sc, K/Rb, Co/Th, La/Sc, and Cr/Th values point to a significant input of detrital volcanic material of calcalkaline felsic and intermediate composition. However, major element geochemistry is not useful for interpreting the tectonic setting. Discrimination plots based on immobile trace elements, such as Ti, Zr, La, Sc, and Th, show that most data lie in the active continental margin field. Geochemical information is not sufficiently sensitive to differentiate the two different source areas recognized by petrographic and modal analyses of conglomerates and sandstones.