Geology and Metallogenesis of the Sawayaerdun Gold Deposit in the Southwestern Tianshan Mountains, Xinjiang, China

The Sawayaerdun gold deposit, located in Wuqia County, Southwest Tianshan, China, occurs in Upper Silurian and Lower Devonian low‐grade metamorphic carbonaceous turbidites. The orebodies are controlled by a series of NE‐NNE‐trending, brittle–ductile shear zones. Twenty‐four gold mineralized zones have been recognized in the Sawayaerdun ore deposit. Among these, the up to 4‐km‐long and 200‐m wide No. IV mineralized zone is economically the most important. The average gold grade is 1–6 g/t. Gold reserves of the Sawayaerdun deposit have been identified at approximately 37 tonnes and an inferred resource of 123 tonnes. Hydrothermal alteration is characterized by silicification, pyritization, arsenopyritization, sericitization, carbonatization and chloritization. On the basis of field evidence and petrographic analysis, five stages of vein emplacement and hydrothermal mineralization can be distinguished: stage 1, early quartz stage, characterized by the occurrence of quartz veins; stage 2, arsenopyrite–pyrite–quartz stage, characterized by the formation of auriferous quartz veinlets and stockworks; stage 3, polymetallic sulfide quartz stage, characterized by the presence of auriferous polymetallic sulfide quartz veinlets and stockworks; stage 4, antimony–quartz stage, characterized by the formation of stibnite–jamesonite quartz veins; and stage 5, quartz–carbonate vein stage. Stages 2 and 3 represent the main gold mineralization, with stage 4 representing a major antimony mineralization episode in the Sawayaerdun deposit. Two types of fluid inclusion, namely H₂O–NaCl and H₂O–CO₂–NaCl types, have been recognized in quartz and calcite. Aqueous inclusions show a wide range of homogenization temperatures from 125 to 340°C, and can be correlated with the mineralization stage during which the inclusions formed. Similarly, salinities and densities of these fluids range for each stage of mineralization from 2.57 to 22 equivalent wt% NaCl and 0.76 to 1.05 g/cm³, respectively. The ore‐forming fluids thus are representative of a medium‐ to low‐temperature, low‐ to medium‐salinity H₂O–NaCl–CO₂–CH₄–N₂ system. The δ³⁴S_CDT values of sulfides associated with mineralization fall into a narrow range of ?3.0 to +2.6‰ with a mean of +0.1‰. The δ¹³C_PDB values of dolomite and siderite from the Sawayaerdun gold deposit range from ?5.4 to ?0.6‰, possibly reflecting derivation of the carbonate carbon from a mixed magmatic/sedimentary source. Changes in physico‐chemical conditions and composition of the hydrothermal fluids, water–rock exchange and immiscibility of hydrothermal fluids are inferred to have played important roles in the ore‐forming process of the Sawayaerdun gold–antimony deposit.     

Carbon dioxide in igneous petrogenesis: II. Fluid dynamics of mantle metasomatism

Petrographic, fluid inclusion, geochemical and isotopic evidence from xenoliths in alkali basalts suggests that low-viscosity fluids rich in O-H-C, dissolved silicates and especially the incompatible elements may ascend, decompress and precipitate crystalline phases and/or induce partial fusion in the upper mantle. Such mantle metasomatic fluids (MMF) may be important in generating isotopic heterogeneity and in transporting and focusing mantle heat. In order to model the movement of MMF, the ordinary differential equations governing the variation ofP, T, ascent velocity and fluid density of a compressible, viscous, single-phase (H₂O or CO₂) non-reacting fluid ascending through a vertical crack of constant width have been solved. A large number of numerical simulations were carried out in which the significant factors affecting flow behavior (thermodynamic and transport fluid properties, roughness and width of cracks, geothermal gradient, initial conditions, etc.) were systematically varied. The calculations show that: (1) MMF tends to move at uniform rates following a short period of rapid initial acceleration, (2) MMF ascends nearly isothermally, (3) MMF acts as an efficient heat transfer agent; numerical experiments show that transport of heat into regions undergoing metasomatism can lead to partial fusion. The heat transported by movement of MMF averaged over the age of the Earth is sufficient to generate about 0.1 km³ of basaltic magma per year, which is approximately equal to the production rate of alkaline magma. If an intense period of mantle degassing occured early in the history of the Earth, the transport of heat and mass (K, U, Rb, LREE) by migrating fluids might have been important.     

The effect of ionic interactions on the oxidation of metals in natural waters

The effect of ionic interactions of the major components of natural waters on the oxidation of Cu(I) and Fe(II) has been examined. The various ion pairs of these metals have been shown to have different rates of oxidation. For Fe(II), the chloride and sulfate ion pairs are not easily oxidized. The measured decrease in the rate constant at a fixed pH in chloride and sulfate solutions agrees very well with the values predicted. The effect of pH (6 to 8) on the oxidation of Fe(II) in water and seawater have been shown to follow the rate equation

$\text{-d in [Fe(II)]/dt = k}{}_1\text{β}{}_1\text{α}{}_{\mathrm{Fe}}\text{/[H}{}^{\mathrm{+}}\text{] + k}{}_2\text{β}{}_2\text{α}{}_{\mathrm{Fe}}\text{/[H}{}^{\mathrm{+}}\text{]}{}^2$

where k₁ and k₂ are the pseudo first order rate constants, β₁ and β₂ are the hydrolysis constants for Fe(OH)⁺ and Fe(OH)⁰. The value of α_FE is the fraction of free Fe²⁺. The value of k₁ (2.0 ±0.5 min^?1) in water and seawater are similar within experimental error. The value of k₂ (1.2 × 10⁵ min^?1) in seawater is 28% of its value in water in reasonable agreement with predictions using an ion pairing model.For the oxidation of Cu(I) a rate equation of the form

$\text{?d ln [Cu(I)]/dt = k}{}_0\text{α}{}_{\mathrm{Cu}}\text{+ k}{}_1\text{β}{}_1\text{α}{}_{\mathrm{Cu}}\text{[Cl]}$

was found where k₀ (14.1 sec^?1) and k₁ (3.9 sec^?1) are the pseudo first order rate constants for the oxidation of Cu⁺ and CuCl⁰, β₁ is the formation constant for CuCl⁰ and α_Cu is the fraction of free Cu⁺. Thus, unlike the results for Fe(II), Cu(I) chloride complexes have measurable rates of oxidation.     

The Gibbs method and Duhem's theorem: The quantitative relationships among P,T, chemical potential,phase composition and reaction progress in igneous and metamorphic systems

The Gibbs method permits simultaneous evaluation of the relationships among all intensive thermodynamic variables of a heterogeneous system in equilibrium. Addition of mass balance constraints permits simultaneous evaluation of both intensive and extensive variables so that changes in phase chemistry and modes may be monitored. Assumption of closed system behavior results in a system of equations with two degrees of freedom, regardless of the thermodynamic variance, as specified by Duhem's theorem. Open system behavior increases the number of degrees of freedom by the number of components to which the system is open. The methodology presented is therefore a formal statement of the constraints among the differential of all of the intensive and extensive variables of a heterogeneous system.Examples of the application of this formalism include contouring pressure-temperature space for mineral composition, modal changes and reaction progress; contouring reaction space with pressure, temperature and mineral composition; and calculation of compositional and modal changes of phases for prescribed changes in pressure and temperature, as, for example, in the calculation of synthetic garnet zoning profiles or liquid lines of descent in crystallizing magmas.     

Elemental and isotopic fractionation of Type B calcium-, aluminum-rich inclusions: experiments, theoretical considerations, and constraints on their thermal evolution

Experiments exposing Type B calcium-, aluminum-rich inclusion (CAI)-like melts at high temperatures to high vacuum or reducing hydrogen-rich gas mixtures were used to determine the rates and consequences of elemental and isotopic fractionation by evaporation. Silicon and magnesium were found to evaporate much faster than calcium and aluminum, and the resulting residual liquid trajectories in composition space are reproduced via a thermodynamic model for the saturation vapor pressure of the evaporating species. Isotopic fractionations associated with evaporation were measured for magnesium. The resulting relationship between fraction of magnesium lost and enrichment of the residue in the heavy isotopes of magnesium follows a Rayleigh fractionation curve with a fractionation factor that is close to, but not exactly, the theoretically expected value. The rate of evaporation is found to be a strong function of temperature, oxygen fugacity, and melt composition, which can be understood and modeled in terms of the dependence of the saturation vapor pressures on these variables. The relationship between evaporation rate, which we measure, and calculated saturation vapor involves empirical evaporation coefficients that we find to be significantly less than one (∼0.1). Analytical and numerical models are used to characterize how diffusion in both the melt and in the surrounding gas affects evaporation rates and the degree of chemical and isotopic fractionation. The experimental data and theoretical considerations are combined to give a parameterization of the rates and consequences of evaporation of Type B CAI-like liquids, which is then used to translate the measured isotopic fractionation of Type B CAIs into constraints on their thermal history. Cooling rates of the order of 10°C per hour are indicated.     

Gold partitioning in melt-vapor-brine systems

We used laser-ablation inductively coupled plasma mass spectrometry to measure the solubility of gold in synthetic sulfur-free vapor and brine fluid inclusions in a vapor + brine + haplogranite + magnetite + gold metal assemblage. Experiments were conducted at 800°C, oxygen fugacity buffered at Ni-NiO (NNO), and pressures ranging from 110 to 145 MPa. The wt% NaCl eq. of vapor increases from 2.3 to 19 and that of brine decreases from 57 to 35 with increasing pressure. The composition of the vapors and brines are dominated by NaCl + KCl + FeCl₂ + H₂O. Gold concentrations in vapor and brine decrease from 36 to 5 and 50 to 28 μg/g, respectively, and the calculated vapor:brine partition coefficients for gold decrease from 0.72 to 0.17 as pressure decreases from 145 to 110 MPa. These data are consistent with the thermodynamic boundary condition that the concentration of gold in the vapor and brine must approach a common value as the critical pressure is approached along the 800°C isotherm in the NaCl-KCl-FeCl₂-HCl-H₂O system.We use the equilibrium constant for gold dissolution as AuOH⁰, extrapolated from lower temperature and overlapping pressure range, to calculate expected concentrations of AuOH⁰ in our experimental vapors. These calculations suggest that a significant quantity of gold in our experimental vapors is present as a non-hydroxide species. Possible chloridogold(I) species are hypothesized based on the positively correlated gold and chloride concentrations in our experimental vapors. The absolute concentration of gold in our synthetic vapor, brine, and melt and calculated mass partition coefficients for gold between these physicochemically distinct magmatic phases suggests that gold solubility in aqueous fluids is a function of aqueous phase salinity, specifically total chloride concentration, at magmatic conditions. However, though we highlight here the effect of salinity, the combination of our data with data sets from lower temperatures evinces a significant decrease in gold solubility as temperature drops from 800°C to 600°C. This decrease in solubility has implications for gold deposition from ascending magmatic fluids.     

Scattered Data Approximation on the Bisphere and Application to Texture Analysis

The paper deals with the approximation and optimal interpolation of functions defined on the bisphere \mathbb S²×\mathbb S²\mathbb {S}^{2}\times \mathbb {S}^{2} from scattered data. We demonstrate how the least square approximation to the function can be computed in a stable and efficient manner. The analysis of this problem is based on Marcinkiewicz–Zygmund inequalities for scattered data which we present here for the bisphere. The complementary problem of optimal interpolation is also solved by using well-localized kernels for our setting. Finally, we discuss the application of the developed methods to problems of texture analysis in material science.     

Holocene environments and climate in the Mongolian Altai reconstructed from the Hoton-Nur pollen and diatom records: a step towards better understanding climate dynamics in Central Asia

This study presents the results of the palynological and diatom analyses of the sediment core recovered in Hoton-Nur Lake (48°37′18″N, 88°20′45″E, 2083 m) in 2004. Quantitative reconstruction of the Holocene vegetation and climate dynamics in the semiarid Mongolian Altai suggests that boreal woodland replaced the primarily open landscape of northwestern Mongolia at about 10 kyr BP (1 kyr = 1000 cal yr) in response to a noticeable increase in precipitation from 200–250 mm/yr to 450–550 mm/yr. A decline of the forest vegetation and a return to a predominance of open vegetation types occurred after 5 kyr BP when precipitation sums decreased to 250–300 mm/yr. Prior to 11.5 kyr BP diatom concentrations are relatively low and the lake is dominated by planktonic Cyclotella and small Fragilariaceae, suggesting the existence of a relatively deep and oligotrophic/mesotrophic lake. The great abundance of Staurosirella pinnata from the beginning of the record until 10.7 kyr BP might imply intensified erosion processes in the catchment and this is fully consistent with the presence of scarce and dry vegetation and the generally arid climate during this period. From about 10.7 kyr BP, more planktonic diatom taxa appeared and increased in abundance, indicating that the lake became more productive as diatom concentration increased. This change correlates well with the development of boreal woodland in the catchment. Decrease in precipitation and changes in the vegetation towards steppe are reflected by the rapid increase in Aulacoseira distans from about 5 kyr BP. The Holocene pollen and diatom records do not indicate soil and vegetation cover disturbances by the anthropogenic activities, implying that the main transformations of the regional vegetation occurred as a result of the natural climate change. Our reconstruction is in agreement with the paleomonsoon records from China, demonstrating an abrupt strengthening of the summer monsoon at 12 kyr BP and an associated increase in precipitation and in lake levels between 11 and 8 kyr BP, followed by the stepwise attenuation of the monsoon circulation and climate aridization towards the modern level. The records from the neighboring areas of Kazakhstan and Russia, situated west and north of Hoton-Nur, demonstrate spatially and temporally different Holocene vegetation and climate histories, indicating that the Altai Mountains as a climate boundary are of pivotal importance for the Holocene environmental and, possibly, habitation history of Central Asia.     

Trace element compositions of submicroscopic inclusions in coated diamond: A tool for understanding diamond petrogenesis

Trace element compositions of submicroscopic inclusions in both the core and the coat of five coated diamonds from the Democratic Republic of Congo (DRC, formerly Zaire) have been analyzed by Laser Ablation Inductively Coupled Mass Plasma Spectrometry (LA-ICP-MS). Both the diamond core and coat inclusions show a general 2-4-fold enrichment in incompatible elements relative to major elements. This level of enrichment is unlikely to be explained by the entrapment of silicate mantle minerals (olivine, garnet, clinopyroxene, phlogopite) alone and thus submicroscopic fluid or glass inclusions are inferred in both the diamond coat and in the gem quality diamond core. The diamond core fluids have elevated High Field Strength Element (Ti, Ta, Zr, Nb) concentrations and are enriched in U relative to inclusions in the diamond coats and relative to chondrite. The core fluids are also moderately enriched in LILE (Ba, Sr, K). Therefore, we suggest that the diamond cores contain inclusions of silicate melt. However, the Ni content and Ni/Fe ratio of the trapped fluid are very high for a silicate melt in equilibrium with mantle minerals; high Ni and Co concentrations in the diamond cores are attributed to the presence of a sulfide phase coexisting with silicate melt in the diamond core inclusions. Inclusions in the diamond coat are enriched in LILE (U, Ba, Sr, K) and La over the diamond core fluids and to chondrite. The coats have incompatible element ratios similar to natural carbonatite (coat fluid: Na/Ba ≈0.66, La/Ta≈130). The coat fluid is also moderately enriched in HFSE (Ta, Nb, Zr) when normalized to chondritic Al. LILE and La enrichment is related to the presence of a carbonatitic fluid in the diamond coat inclusions, which is mixed with a HFSE-rich hydrous silicate fluid similar to that in the core. The composition of the coat fluid is consistent with a genetic link to group 1 kimberlite.