A. A. Ariskin

Calculation of water-bearing primary basalt and estimation of source mantle conditions beneath arcs: PRIMACALC2 model for WINDOWS

We present a new method for estimating the composition of water-bearing primary arc basalt and its source mantle conditions. The PRIMACALC2 model uses a thermodynamic fractional crystallization model COMAGMAT3.72 and runs with an Excel macro to examine the mantle equilibrium and trace element calculations of a primary basalt. COMAGMAT3.72 calculates magma fractionation in 0–10 kb at various compositions, pressure, oxygen fugacity, and water content, but is only applicable for forward calculations. PRIMACALC2 first calculates the provisional composition of a primary basalt from an observed magma. The basalt composition is then calculated by COMAGMAT3.72 for crystallization. Differences in elemental concentrations between observed and the closest-match calculated magmas are then adjusted in the primary basalt. Further iteration continues until the calculated magma composition converges with the observed magma, resulting in the primary basalt composition. Once the fitting is satisfied, back calculations of trace elements are made using stepwise addition of fractionated minerals. Mantle equilibrium of the primary basalt is tested using the Fo-NiO relationship of olivine in equilibrium with the primary basalt, and thus with the source mantle. Source mantle pressure, temperature, and degree of melting are estimated using petrogenetic grids based on experimental data obtained in anhydrous systems. Mantle melting temperature in a hydrous system is computed by adjusting T with a parameterization for a water-bearing system. PRIMACALC2 can be used either in dry or water-bearing arc magmas and is also applicable to mid-ocean ridge basalts and nonalkalic ocean island basalts.

The Dovyren Intrusive Complex: Problems of Petrology and Ni Sulfide Mineralization

This paper presents a review of petrological-geochemical studies at the Yoko-Dovyren Massif with an emphasis on relations between parameters of the parental magma, a model for the genesis of the lower contact zone, and the nature of Ni sulfide ore mineralization, including the evaluation of the possible ore potential. Arguments are presented in support of the conclusion that the Dovyren magma brought much intratelluric olivine of the composition Fo85–87 into the chamber, and the composition of the initial melt corresponded to gabbronorite or moderately magnesian basite with no more than 10 wt % MgO. The probable temperature of the parental magma was approximately 1200–1250°C, and the sulfur solubility did not exceed 0.10–0.12 wt % (P = 1 kbar, WM buffer). The comparison of this estimate with the average S contents in the bottom plagioperidotites (0.12±0.06 wt %) indicates that the initial magma was saturated with a sulfide phase. For the first time the problem of the composition of contaminated dunites was formulated (these rocks occur in the Layered Series and contain more magnesian olivine Fo87–92). The reason for the increase in the mg# of olivine is thought to be the partial melting and compaction of the original cumulates due to the infiltration of intercumulus melt enriched in volatile components. The volatiles were presumably provided by the thermal decomposition of carbonate xenoliths, a process that resulted in an increase in the CO2 pressure and the transfer of calcite-magnesite components of carbonates into the melt. This follows from (1) the occurrence of magnesian skarn developing after carbonates, (2) high CaO contents in olivine form the contaminated dunite, (3) the appearance of olivine-bearing pyroxenites and wehrlites in the upper part of the dunite zone, (4) correlation between the olivine and chromite composition in the contaminated and uncontaminated dunites, (5) broad variations in the oxygen isotopic composition of olivine and plagioclase from rocks of the Layered Series, (6) experimental data on the dissolution of carbonates in alkali basalt melts, and (7) analogies with isotopic-geochemical characteristics of rocks from the Jinchuan ultramafic complex. Petrological implications of the interpretation of the Dovyren chamber are discussed with reference to closed and flow-through (during an initial stage) magmatic systems. A petrological-geological model is proposed for the genesis of the Synnyr-Dovyren volcanic-plutonic complex and related Ni sulfide ore mineralization. The potential resources of Cu-Ni sulfide ores in the plagioperidotites are evaluated with regard to the still-unexposed part of the massif.

Geochronology of the Dovyren intrusive complex, northwestern Baikal area, Russia, in the Neoproterozoic

The paper reports newly obtained data on the geochronology of the Dovyren intrusive complex and associated metarhyolites of the Inyaptuk Formation in the Synnyr Range. The data were obtained by local LA-ICPMS analysis of zircons in samples. The U-Pb age of olivine-free gabbronorite from near the roof of the Yoko-Dovyren Massif is 730 ± 6 Ma (MSWD = 1.7, n = 33, three samples) is close to the estimated age of 731 ± 4 Ma (MSWD = 1.3, n = 56, five samples) of a 200-m-thick sill beneath the pluton. These data overlap the age of recrystallized hornfels found within the massif (“charnockitoid”, 723 ± 7 Ma, MSWD = 0.12, n = 10) and a dike of sulfidated gabbronorite below the bottom of the massif (725 ± 8 Ma, MSWD = 2.0, n = 15). The estimates are also consistent with the age of albite hornfels (721 ± 6 Ma, MSWD = 0.78, n = 12), which was produced in a low-temperature contact metamorphic facies of the host rocks. The average age of the Dovyren Complex is 728.4 ± 3.4 Ma (MSWD = 1.8, n = 99) based on data on the sill, near-roof gabbronorite, and “charnockitoid”) and is roughly 55 Ma older than the estimate of 673 ± 22 Ma (Sm-Nd; [13]). The U-Pb system of zircon in two quartz metaporphyre samples from the bottom portion of the Inyaptuk volcanic formation in the northeastern part of the Yoko-Dovyren Massif turns out to be disturbed. The scatter of the data points can be explained by the effect of two discrete events. The age of the first zircon population is then 729 ± 14 Ma (MSWD = 0.74, n = 8), and that of the second population is 667 ± 14 Ma (MSWD = 1.9, n = 13). The older value pertains to intrusive rocks of Dovyren, and the age of the “rejuvenated” zircon grains corresponds to the hydrothermal-metasomatic processes, which affected the whole volcano-plutonic sequence and involved the serpentinization of the hyperbasites. This is validated by the results of Rb-Sr isotopic studies with the partial acid leaching of two serpentinized peridotite samples from the Verblyud Sill. These studies date the overprinted processes at 659 ± 5 Ma (MSWD = 1.3, n = 3).

Crystallization differentiation of intrusive magmatic melt: Development of a convection-accumulation model

The paper summarizes the principal results obtained over the past three decades at the Vernadsky Institute and the Department of Geochemistry of the Moscow State University by the computer simulation of basaltic magma differentiation in magma chambers. The processes of diffusion-controlled mass transfer in a chamber are demonstrated to be principally limited by the heat resources of the cooling magma and cannot play any significant role during the large-scale partitioning of melt components. The leading mass-transfer mechanism is the settling of crystals from convecting magma in the form of suspension flows that are enriched and depleted in the solid phase. The physical prerequisite for the onset of this concentration convection is the existence of boundary layers, which are characterized by volume crystallization and a gradient distribution of the suspended phases. Considered in detail are the principles used in the development of algorithms with regard for the occurrence of a boundary layer and the “instantaneous” stirring of the crystallizing magma that does not hamper the settling of mineral grains forming the cumulus. The plausibility of the convection-accumulation model is illustrated by the example of the reconstructed inner structure of differentiated Siberian traps. In application to these bodies, it is demonstrated that the solutions of the forward and inverse simulation problems with the use of geochemical thermometry techniques are identical. This is a convincing argument for the predominance of convection-accumulation processes during the formation of thin tabular magmatic bodies. The further development of the computer model for the differentiation dynamics should involve the processes of compositional convection related to the migration and reactivity of the intercumulus melt.

Test of the Ballhaus–Berry–Green Ol–Opx–Sp oxybarometer and calibration of a new equation for estimating the redox state of melts saturated with olivine and spinel

Testing the Ballhaus–Berry–Green Ol–Opx–Sp oxybarometer (BBG) on independent experimental data indicates that it overestimates the oxygen fugacity by 0.6–1.3 log units under mildly reduced conditions (near the C–CO buffer) and by as much as 2–3 log units under reduced conditions (at the IW buffer and below it). A newly developed oxibarometer is suggested to minimize this effect and enhance the capabilities of redoxometry of low-pressure mineral associations, including magmatic melts undersaturated with respect to orthopyroxene (Opx). The new empirical equation of the oxybarometer is applicable to a wide range of mafic–ultramafic magmas of normal alkalinity, including terrestrial, lunar, and meteoritic systems under pressures of 0.001–25 kbar and oxygen fugacity ranging from IW–3 to NNO + 1. The derived regression fits the ΔQFM values of the calibration dataset (154 experiments) accurate to ~0.5 log units. The new oxybarometer eliminates systematic errors when redox parameters are evaluated for the reduced region (from IW–3 to C–CO) and for crystallization of magmas without Opx on the liquidus. The efficiency of the suggested model is demonstrated by its application to natural rocks: (1) low-Ti lunar basalts, (2) tholeiites from the Shatsky Rise, (3) Siberian flood basalts, (4) rocks of the layered series of the Yoko-Dovyren intrusion, and (5) mantle xenoliths collected in southern Siberia, Mongolia, China, and the southern Russian Far East. The values yielded by such oxybarometers for intrusive rocks, which underwent long-lasting cooling and postcumulus reequilibration, should be regarded with reserve.

Geochemical constraints on petrogenic processes on Venus

Abundances of three incompatible elements (K, U, Th) have been determined by robotic spacecraft in five surface materials on Venus. We present these data normalized to terrestrial normal mid-ocean ridge basalt (NMORB) derived from the depleted mantle. Relative to NMORB, all of the Venus materials studied are enriched in all of these elements. The moderately enriched (Venera 9 and 10, Vega 1 and 2) basaltic rocks are similar to one another in the sense of their enrichment trends (U N > K N ) but differ from the highly enriched Venera 8 material, where the trend is K N > U N . This difference implies that the enrichment pattern for the basaltic materials is not controlled by crystallization of an NMORB-like magma or by contamination of such a magma by highly enriched Venera 8 material within the crust, and the Venera 8 material cannot have evolved from the magma of any of the basaltic rocks. Our calculations show the K-U-Th pattern for any of the Venus rocks analyzed was not controlled by batch partial melting of primitive mantle. The Venera 8 material could be produced as a partial melt from eclogitic tholeiite, but none of the basaltic materials could be because the sense of their U N > K N trends is opposite to the enrichment trends in calculated models (K N > U N ). The Venus basalts differ from fresh terrestrial rocks (NMORBs and oceanic island-arc volcanics) in having nearly constant K/U ratios, while terrestrial rocks have nearly constant Th/U ratios. This may suggest an unusual composition of mantle source(s) of the Venus basalts and/or unusual fractionation process(es) on that planet.

Cu–Ni–PGE fertility of the Yoko-Dovyren layered massif (northern Transbaikalia, Russia): thermodynamic modeling of sulfide compositions in low mineralized dunite based on quantitative sulfide mineralogy

The geology and major types of sulfide mineralization in the Yoko-Dovyren layered massif (northern Transbaikalia, Russia) are presented. This study focuses on the structure, mineralogy, and geochemistry of poorly mineralized plagiodunite and dunite in the lower part of the intrusion. Assuming these rocks contain key information on the timing of sulfide immiscibility in the original cumulate pile, we apply a novel approach which combines estimates of the average sulfide compositions in each particular rock with thermodynamic modeling of the geochemistry of the original sulfide liquid. To approach the goal, an updated sulfide version of the COMAGMAT-5 model was used. Results of simulations of sulfide immiscibility in initially S-undersaturated olivine cumulates demonstrate a strong effect of the decreasing fraction of the silicate melt, due to crystallization of silicate and oxide minerals, on the composition of the intercumulus sulfide liquid. Comparison of the observed and modeled sulfide compositions indicates that the proposed modeling reproduces well the average concentrations of Cu, Cd, Ag, and Pd in natural sulfides. This suggests the sulfide control on the distribution of these elements in the rocks. Conversely, data for Pt and Au suggest that a significant portion of these elements could present in a native form, thus depleting the intercumulus sulfide melt at an early stage of crystallization.

Simulation of Molecular Mass Distributions and Evaluation of O 2- Concentrations in Polymerized Silicate Melts

A new statistical model is proposed for the molecular mass distributions (MMD) of polymerized anions in silicate melts. The model is based on the known distribution of Qn species in the MeO-Me2O-SiO2 system. In this model, chain and ring complexes are regarded as a random series of Qn structons with various concentrations of bridging bonds (1 ≤ n ≤ 4, Q0 corresponds to SiO44−). This approach makes it possible to estimate the probability of formation of various ensembles of polymer species corresponding to the general formula (SiiO3i+1−j)2(i+1−j)−, where i is the size of the ion, and j is the cyclization number of intrachain bonds. The statistical model is utilized in the STRUCTON computer model, which makes use of the Monte Carlo method and is intended for the calculation of the composition and proportions of polyanions at a specified degree of polymerization of silicate melts (STRUCTON, version 1.2; 2007). Using this program, we simulated 1200 MMD for polyanions in the range of 0.52 ≤ p ≤98, where p is the fraction of nonbridging bonds in the silicon-oxygen matrix. The average number of types of anions in this range was determined to increase from three (SiO44−, Si2O76−, and Si3O108−) to 153, and their average size increases from 1 to 7.2. A special option of the STRUCTON program combines MMD reconstructions in silicate melts with the formalism of the Toop-Samis model, which enables the calculation of the mole fraction of the O2− ion relative to all anions in melts of specified composition. It is demonstrated that, with regard for the distribution and average size of anion complexes, the concentration of the O2− ion in the MeO-SiO2 system is characterized by two extrema: a minimum at 40–45 mol % SiO2, which corresponds to the initial stages of the gelenization of the polycondensated silicate matrix, and a maximum, which is predicted for the range of 60–80 mol % SiO2.

Inner Structure, Composition, and Genesis of the Chineiskii Anorthosite-Gabbronorite Massif, Northern Transbaikalia

The Chineiskii anorthosite-gabbronorite massif is the most typical layered intrusion in Russia, which is accompanied by large V and Cu deposits. This massif is first considered to be a component of the Proterozoic volcanic-plutonic system of the Kodar-Udokan district, whose largest massifs are Chineiskii and Lukturskii. This system also comprises numerous dikes (including the Main gabbronorite dike at the Udokan deposit, whose thickness reaches 200 m), which are likely the magmatic feeders of ancient volcanism. An intermediate position in the vertical section of the magmatic system is occupied by gabbroids, whose exposures occur in the peripheral part of the Lurbunskii granite massif. The intrusive rocks were proved to be genetically interrelated and show certain similar geochemical features: they bear elevated TiO2 concentrations and have similar trace element patterns and (La/Sm)N and (Gd/Yb)N ratios (1.5–2.3 and 1.87–2.06, respectively). The Chineiskii Massif is thought to have been formed by the successive emplacement of genetically similar basic magmas, which produced four rock groups with fine and coarse layering and cyclicity of variable rank (microrhythms, rhythms, units, and series). The results of cluster analysis indicate that the rocks can be classified into 13 petrochemical types. The phase and chemical characteristics of the parental melts of these compositions were simulated with the use of the COMAGMAT-3.5 computer model, which was also applied to evaluate the composition of the most primitive initial magma of the whole Chineiskii Massif. Our results indicate that the primitive magma was heterogeneous (olivine + plagioclase ± titanomagnetite + melt) at a temperature of approximately 1130°C. The initial melt had a ferrobasaltic composition and was close to saturation with magnetite at ∼NNO ± 0.5

Geochemical nature and age of the plagiogranite-gabbronorite association of the oceanic core complex of the Mid-Atlantic ridge at 5°10′S

Our newly obtained data on the geochemistry and age of plagiogranite-gabbronorite association in the oceanic core complex of the Mid-Atlantic Ridge (MAR) at 5°10′S suggest close genetic relations between these rocks in this segment of the ridge. The U/Pb zircon age of an oceanic plagiogranite (OPG) sample is 1.059 ± 0.055 Ma and is in good agreement with the zircon age of plutonic rocks in the oceanic core complex of northern MAR. A distinctive geochemical feature of the rocks is their unusually depleted 87Sr/86Sr and 143Nd/144Nd ratios, which suggest that the plutonic rocks of the gabbronorite-plagiogranite association in MAR at 5°10′S could be derived from the most strongly depleted mantle reservoir of all known to occur beneath the axial MAR zone. The COMAGMAT-5.2 numerical thermodynamic simulation of the possible crystallization links between the plagiogranite and gabbronorite from the MAR segment at 5°10′S led us to conclude that the leading role in the origin of the plagiogranite was played by a two-stage process: the partial melting of the gabbronorite and the subsequent fractionation of the newly generated melt. The regional differences between the isotopic-geochemical parameters of MAR plagiogranites can, perhaps, reflect local specifics of so-called hydrothermal anatexis, such as the geochemical features of the rocks involved in this process and the parameters of the hydrothermal process, for example, variations in the W/R ratio.