Sequestration of As and Mo in uranium mill precipitates (pH 1.5–9.2): An XAS study

As- and Mo- bearing secondary mineral phases formed during the neutralization of uranium mill wastes require characterization. Previous studies indicate that arsenate and molybdate adsorbed to ferrihydrite are the dominant controls in the tailings materials. A lab-scale plant was employed to characterize secondary precipitates from a variety of ore blends. Through total elemental analysis of precipitates and As and Mo K-edge X-ray absorption spectroscopy, different ratios of contributing phases were determined for each pH stage (4.2, 6.5, and 9.2) of the neutralization process. Overall, arsenate adsorbed to ferrihydrite was the dominant As mineral phase regardless of pH or sample blend (53–77%), with fractional contribution from ferric arsenates, and adsorption to aluminum phases. Molybdate adsorbed to ferrihydrite was the dominant Mo mineral phase, with fractional contribution decreasing with increasing pH (100–69%). The characterization of these phases in the secondary precipitates provides further understanding of the contributing mineral species in tailing facilities.     

Preservation of As III and As V in surface and ground waters

The literature on preservation procedures for the stabilization of As redox species in waters is controversial. A paper by McCleskey et al. (Appl. Geochem., 2004) states that any field collection procedure that filters out microorganisms, adds a reagent t…     

An experimental study of dissolution–reprecipitation in fluorapatite: fluid infiltration and the formation of monazite

In a series of timed experiments, monazite inclusions are induced to form in the Durango fluorapatite using 1 and 2 N HCl and H₂SO₄ solutions at temperatures of 300, 600, and 900°C and pressures of 500 and 1,000 MPa. The monazite inclusions form only in reacted areas, i.e. depleted in (Y+REE)+Si+Na+S+Cl. In the HCl experiments, the reaction front between the reacted and unreacted regions is sharp, whereas in the H₂SO₄ experiments it ranges from sharp to diffuse. In the 1 N HCl experiments, Ostwald ripening of the monazite inclusions took place both as a function of increased reaction time as well as increased temperature and pressure. Monazite growth was more sluggish in the H₂SO₄ experiments. Transmission electron microscopic (TEM) investigation of foils cut across the reaction boundary in a fluorapatite from the 1 N HCl experiment (600°C and 500 MPa) indicate that the reacted region along the reaction front is characterized by numerous, sub-parallel, 10–20 nm diameter nano-channels. TEM investigation of foils cut from a reacted region in a fluorapatite from the 1 N H₂SO₄ experiment at 900°C and 1,000 MPa indicates a pervasive nano-porosity, with the monazite inclusions being in direct contact with the surrounding fluorapatite. For either set of experiments, reacted areas in the fluorapatite are interpreted as replacement reactions, which proceed via a moving interface or reaction front associated with what is essentially a simultaneous dissolution–reprecipitation process. The formation of a micro- and nano-porosity in the metasomatised regions of the fluorapatite allows fluids to permeate the reacted areas. This permits rapid mass transfer in the form of fluid-aided diffusion of cations to and from the growing monazite inclusions. Nano-channels and nano-pores also serve as sites for nucleation and the subsequent growth of the monazite inclusions.     

A combined IR and XRD study of natural well crystalline goethites (α-FeOOH)

Acta Geochimica - Goethite (α-FeOOH) is one of the most abundant minerals on the Earth surface, occurring in temperate, tropical and equatorial climates. Fe in goethite can be substituted by...     

A Compensated-Redlich-Kwong (CORK) equation for volumes and fugacities of CO2 and H2O in the range 1 bar to 50 kbar and 100–1600°C

We present a simple virial-type extension to the modified Redlich-Kwong (MRK) equation for calculation of the volumes and fugacities of H₂O and CO₂ over the pressure range 0.001–50 kbar and 100 to 1400°C (H₂O) and 100 to 1600°C (CO₂). This extension has been designed to: (a) compensate for the tendency of the MRK equation to overestimate volumes at high pressures, and (b) accommodate the volume behaviour of coexisting gas and liquid phases along the saturation curve. The equation developed for CO₂ may be used to derive volumes and fugacities of CO, H₂, CH₄, N₂, O₂ and other gases which conform to the corresponding states principle. For H₂O the measured volumes of Burnham et al. are significantly higher in the range 4–10 kbar than those presented by other workers. For CO₂ the volume behaviour at high pressures derived from published MRK equations are very different (larger volumes, steeper (P/T)_V, and hence larger fugacities) from the virial-type equations of Saxena and Fei. Our CORK equation for CO₂ yields fugacities which are in closer agreement with the available high pressure experimental decarbonation reactions.     

A study of outburst activity of Z and in 2000–2010

We present the results of our study of the classical symbiotic star Z And during its period of activity in 2000–2010. In this period, the system experienced a series of six outbursts, for three of which high-resolution spectra were obtained and analyzed. These observations provided information about the system’s behavior during the entire activity period, rather than during an individual outburst. In particular, we found a fundamental difference between the first outburst, which initiated the activity period, and subsequent outbursts, namely, the presence of bipolar collimated optical outflows for some of the outbursts. We propose a model that can explain all the spectroscopic phenomena observed during this series of outbursts, as well as previous series of outbursts of Z And, and suggest that similar scenarios may be valid for other classical symbiotic stars.     

A Raman spectroscopic study on the structural disorder of monazite–(Ce)

This study addresses whether Raman spectra can be used to estimate the degree of accumulated radiation damage in monazite-(Ce) samples whose chemical composition was previously determined. Our results indicate that the degree of disorder in monazite–(Ce), as observed from increasing Raman band broadening, generally depends on both the structural state (i.e., radiation damage) and the chemical composition (i.e., incorporation of non-formula elements). The chemical effects were studied on synthetic orthophosphates grown using the Li-Mo flux method, and non radiation-damaged analogues of the naturally radiation-damaged monazite–(Ce) samples, produced by dry annealing. We found that the “chemical” Raman-band broadening of natural monazite–(Ce) can be predicted by the empirical formula, $$ {\hbox{FWHM}} {\hbox{[c}}{{\hbox{m}}^{ - {1}}}{]} = {3}{.95} + {26}{.66} \times {\hbox{(Th}} + {\hbox{U}} + {\hbox{Ca}} + {\hbox{Pb)}} {\hbox{[apfu]}} $$ where, FWHM = full width at half maximum of the main Raman band of monazite–(Ce) (i.e., the symmetric PO₄ stretching near 970?cm^?1), and (Th+U+Ca+Pb) = sum of the four elements in apfu (atoms per formula unit). Provided the chemical composition of a natural monazite–(Ce) is known, this “chemical band broadening” can be used to estimate the degree of structural radiation damage from the observed FWHM of the ν₁(PO₄) band of that particular sample using Raman spectroscopy. Our annealing studies on a wide range of monazite–(Ce) reference materials and other monazite–(Ce) samples confirmed that this mineral virtually never becomes highly radiation damaged. Potential advantages and the practical use of the proposed method in the Earth sciences are discussed.     

A photometric study of the polar MT draconis in 2005–2009

We have obtained data on various brightness states of the polar MT Dra over five years of observations, including the first multicolor photometry for this close binary. We confirm the known orbital period, which has remained constant over 17 years, which is unusual for cataclysmic binaries. Our observations in October 2006 demonstrated a transition of the polar to its low brightness state within one day.     

An experimental study of the partitioning of a rare earth element (Gd) in the system diopside—aqueous vapour

The partitioning of Gd in the experimental system diopside-aqueous vapor as a function of temperature, pressure, composition of the phases, time, grain size, solid-liquid ratio and Gd concentration has been investigated. A radioactive tracer measurement was used to determine Gd concentration in the separated phases. Diposides were reacted with aqueous vapor containing tracer Gd and reversibility was tested by reacting Gd-doped diopsides with pure aqueous vapor. Equilibration of Gd between the bulk of the diopside and the liquid was found to be limited by the slow rate of Gd diffusion in diopside, maximum value of D = 2 × 10^?15cm²sec^?1 at 800°C and 1 kb. Depending on whether the diopside was previously synthesized or synthesized from an oxide mix during the experiment, Gd concentrations were zoned in the crystal such that higher concentrations existed at the edges or center, respectively. Equilibrium is difficult to achieve in these experiments, but at the optimum experimental conditions for equilibration, the Gd diopside-aqueous vapor distribution coefficient is 20 ± 6 (800°C, 1 kb) in approximate agreement with previous results of 55 ± 23. Changing the composition of the aqueous vapor indicated that possible mechanisms for Gd substitution included coupling of Gd³⁺ with H⁺ or Na⁺ replacing 2Ca²⁺, or substitution of 2Gd³⁺ for 3Ca²⁺ with formation of a cation vacancy.     

Mantle and crustal sources of Archean anorthosite: a combined in situ isotopic study of Pb–Pb in plagioclase and Lu–Hf in zircon

Isotopic analyses of ancient mantle-derived magmatic rocks are used to trace the geochemical evolution of the Earth’s mantle, but it is often difficult to determine their primary, initial isotope ratios due to the detrimental effects of metamorphism and secondary alteration. We present in situ analyses by LA-MC-ICPMS for the Pb isotopic compositions of igneous plagioclase (An_75–89) megacrysts and the Hf isotopic compositions of BSE-imaged domains of zircon grains from two mantle-derived anorthosite complexes from south West Greenland, Fiskenæsset and Nunataarsuk, which represent two of the best-preserved Archean anorthosites in the world. In situ LA-ICPMS U–Pb geochronology of the zircon grains suggests that the minimum crystallization age of the Fiskenæsset complex is 2,936 ± 13 Ma (2σ, MSWD = 1.5) and the Nunataarsuk complex is 2,914 ± 6.9 Ma (2σ, MSWD = 2.0). Initial Hf isotopic compositions of zircon grains from both anorthosite complexes fall between depleted mantle and a less radiogenic crustal source with a total range up to 5 ε_Hf units. In terms of Pb isotopic compositions of plagioclase, both anorthosite complexes share a depleted mantle end member yet their Pb isotopic compositions diverge in opposite directions from this point: Fiskenæsset toward a high-μ, more radiogenic Pb, crustal composition and Nunataarsuk toward low-μ, less radiogenic Pb, crustal composition. By using Hf isotopes in zircon in conjunction with Pb isotopes in plagioclase, we are able to constrain both the timing of mantle extraction of the crustal end member and its composition. At Fiskenæsset, the depleted mantle melt interacted with an Eoarchean (~3,700 Ma) mafic crust with a maximum ¹⁷⁶Lu/¹⁷⁷Hf ~0.028. At Nunataarsuk, the depleted mantle melt interacted with a Hadean (~4,200 Ma) mafic crust with a maximum ¹⁷⁶Lu/¹⁷⁷Hf ~0.0315. Evidence from both anorthosite complexes provides support for the long-term survival of ancient mafic crusts that, although unidentified at the surface to date, could still be present within the Fiskenæsset and Nunataarsuk regions.     

A combined EMPA and LA-ICP-MS study of Li-bearing mica and Sn–Ti oxide minerals from the Qiguling topaz rhyolite (Qitianling District,China): The role of fluorine in origin of tin mineralization

The Sn-rich Qiguling topaz rhyolite dike intrudes the Qitianling biotite granite of the Nanling Range in southern China; the granite hosts the large Furong Sn deposit. The rhyolite dike is typically peraluminous, volatile-enriched, and highly evolved. Whole-rock F and Sn concentrations attain 1.9 wt.% and 2700 ppm, respectively. The rhyolite consists of a fine-grained matrix formed by quartz, feldspar, mica and topaz, enclosing phenocrysts of quartz, feldspar and mica; it is locally crosscut by quartz veinlets. Lithium-bearing micas in both phenocrysts and the groundmass can be classified as primary zinnwaldite, “Mus-Ann” (intermediate member between annite and muscovite), and secondary Fe-rich muscovite. Topaz is present in the groundmass only; common fluorite occurs in the groundmass and also in a specific cassiterite, rutile and fluorite (Sn–Ti–F) assemblage. Cassiterite and rutile are the only Sn and Ti minerals; both cassiterite and Nb-rich rutile are commonly included in the phenocrysts. The Sn–Ti–F assemblage is pervasive, and contains spongy cassiterite in some cases; cassiterite also occurs in quartz veinlets which cut the groundmass. Electron microprobe and LA-ICP-MS compositions were used to study the magmatic and hydrothermal processes and the role of F in Sn mineralization. The presence of zinnwaldite and “Mus-Ann”, which are respectively representative of early and late mica crystallization during magma differentiation, also suggests a significant decrease in f(HF)/f(H₂O) of the system. Cassiterite included in the zinnwaldite phenocrysts is suggested to have crystallized from the primary magma at high temperature. Within the Sn–Ti–F aggregates, rutile crystallized as the earliest mineral, followed by fluorite and cassiterite. Spongy cassiterite containing inclusions of the groundmass minerals indicate a low viscosity of the late fluid. The cassiterite in the quartz veinlets crystallized from low-temperature hydrothermal fluids, which possibly mixed with meteoric water. In general, cassiterite precipitated during both magmatic and hydrothermal stages, and over a range of temperatures. The original fluorine and tin enrichments, f(HF)/f(H₂O) change in the residual magma, formation of Ca,Sn,F-rich immiscible fluid, decrease of the f(HF) during groundmass crystallization, and mixing of magma-derived fluids with low-saline meteoric water during the late hydrothermal stage, are all factors independently or together responsible for the Sn mineralization in the Qiguling rhyolite.     

An experimental study of the Fe oxidation states in garnet and clinopyroxene as a function of temperature in the system CaO–FeO–Fe2O3–MgO–Al2O3–SiO2: implications for garnet–clinopyroxene geothermometry

Samples with eclogitic composition in the system CaO–FeO–Fe₂O₃–MgO–Al₂O₃–SiO₂ were produced from various kinds of starting materials held in graphite-lined Pt capsules at a pressure of 2.5–3.0 GPa and temperatures of 800–1,300 °C using a piston-cylinder or Belt apparatus. Garnets and clinopyroxenes were characterized by analytical transmission electron microscopy and electron probe micro-analysis (EPMA). Fe³⁺/ΣFe ratios determined by electron energy-loss spectroscopy (EELS) decrease in clinopyroxene from 22.2 ± 3.4 % at 800 °C to 13.3 ± 5.4 % at 1,300 °C, while in garnet, they vary between 10.8 ± 1.5 and 15.4 ± 4.7 %, respectively. Temperature estimates according to Krogh (Contrib Mineral Petrol 99:44–48, 1988) reproduce the experimental temperature to ±60 °C without systematic deviations if total iron is used in the calculation. If only the Fe²⁺ content is used, which was obtained by combining EPMA and EELS results, the experimental temperature is underestimated by 33 °C on average at 800–1,200 °C and overestimated by 77 °C on average at 1,300 °C. These systematic deviations can be explained by the temperature-dependent ratio of Fe²⁺/ΣFe in garnet divided by that in clinopyroxene. Since the difference between the calculated and experimental temperature is relatively small, a Fe²⁺-based recalibration of the thermometer appears not to be necessary for the investigated system in the range of pressure, temperature and composition covered by the experiments of this study.     

Crustal structure of the Alpine–Pannonian transition zone: a combined seismic and gravity study

 The crustal structure of the transition zone between the Eastern Alps and the western part of the Pannonian depression (Danube basin) is traditionally interpreted in terms of subvertical Tertiary strike-slip and normal faults separating different Alpine tectonic units. Reevaluation of approximately 4000-km-long hydrocarbon exploration reflection seismic sections and a few deep seismic profiles, together with data from approximately 300 wells, suggests a different structural model. It implies that extensional collapse of the Alpine orogene in the Middle Miocene was controlled by listric normal faults, which usually crosscut Alpine nappes at shallow levels, but at depth merge with overthrust planes separating the different Alpine units. The alternative structural model was tested along a transect across the Danube basin by gravity model calculations, and the results show that the model of low-angle extensional faulting is indeed viable. Regarding the whole lithosphere of the western Pannonian basin, gravity modelling indicates a remarkable asymmetry in the thickness minima of the attenuated crust and upper mantle. The approximately 160 km lateral offset between the two minima suggests that during the Miocene extension of the Pannonian basin detachment of the upper crust from the mantle lithosphere took place along a rheologically weak lower crust. Received: 13 July 1998 / Accepted: 18 March 1999     

Optimization of As(V) adsorption on Fe-RH-MCM-41-immobilized GAC using Box–Behnken Design: Effects of pH,loadings, and initial concentrations

This work focused on the removal from water of arsenate in the 100 to 500 μg/L range using Fe-RH-MCM-41 immobilized on GAC by a colloidal impregnation method. A Response Surface Method (RSM) based on the Box–Behnken Design (BBD) was used to evaluate the effects of pH, adsorbent loading, and initial arsenate concentrations. The thermal stability of Fe-RH-MCM-41/GAC was greater than 460 °C because Fe-RH-MCM-41 and epoxy resin covered the GAC, preventing GAC weight loss. The rate and capacity for arsenate adsorption on Fe-RH-MCM-41/GAC were higher than for GAC. The effects of pH, adsorbent loading, and initial arsenate concentration showed significant effects on efficiency, in terms of percent arsenate removal. Significant interaction occurred as positive effects between arsenate and arsenate concentrations, and pH and arsenate concentrations, and as negative effects between pH and pH, adsorbent loading and adsorbent loading, and adsorbent loading and arsenate concentrations.     

A comparison study of the C–A and S–A models with singularity analysis to identify geochemical anomalies in covered areas

Fractal/multifractal modeling of geochemical data is an interesting topic in the field of applied geochemistry. Identification of weak anomalies for mineral exploration in covered areas is one of the most challenging tasks for utilization of geochemical data. In this study, three fractal models, consisting of the concentration–area (C–A), spectrum–area (S–A) and singularity index models were applied to identify geochemical anomalies in the covered area located in the Chaobuleng Fe polymetallic district, Inner Mongolia (China). The results show that (1) the grassland cover weakens the concentrations of geochemical elements; (2) the C–A model has a limitation to identify weak anomalies in covered areas; (3) the S–A model is a powerful tool to decompose mixed geochemical patterns into a geochemical anomaly map and a varied geochemical background map but suffers edge effects in an irregular shaped study area; and (4) the singularity index is a useful tool to identify weak geochemical anomalies.     

An experimental study of phase equilibria in the system H2O-CO2-NaCl at 800 °C and 9 kbar

Phase equilibria in the ternary system H₂O-CO₂-NaCl were studied at 800 °C and 9 kbar in internally heated gas pressure vessels using a modified synthetic fluid inclusion technique. The low rate of quartz overgrowth along the `b' and `a' axes of quartz crystals was used to avoid fluid inclusion formation during heating, prior to attainment of equilibrium run conditions. The density of CO₂ in the synthetic fluid inclusions was calibrated using inclusions in the binary H₂O-CO₂ system synthesised by the same method and measured on the same heating-freezing stage. In the two-phase field, two types of fluid inclusions with different densities of CO₂ were observed. Using mass balance calculations, these inclusions are used to constrain the miscibility gap and the orientation of two-phase tie-lines in the H₂O-CO₂-NaCl system at 800 °C and 9 kbar. The equation of state of Duan et al. (1995) approximately describes the P-T section of the ternary system up to about 40 wt% of NaCl. At higher NaCl concentrations the measured solubility of CO₂ in the brine is much smaller than predicted by the EOS. A “salting out” effect must be added to the equation of state to include coulomb interaction in the model of Anderko and Pitzer (1993) and Pitzer and Jiang (1996). The new experimental data together with published data up to 5 kbar (Shmulovich et al. 1995) encompass practically all subsolidus crustal P-T conditions. A feature of the new experimental results is the large compositional range in the H₂O-CO₂-NaCl system occupied by the stability fields of halite + CO₂-rich fluid ± H₂O-NaCl brine. The prediction of halite stability in equilibrium with CO₂-rich fluid in deep-crustal rocks is supported by recent petrological and fluid inclusion studies of granulites. Received: 29 June 1998 / Accepted: 17 March 1999     

Application of the content–area (C–A) fractal model and prediction–area (P–A) plot for mineral prospectivity modeling in the Luchun area of Yunnan Province,China

This method of assigning weights based on expert opinion introduces bias when we are evaluating the relative importance of evidence values. In this paper, we used a prediction–area (P–A) plot method and content–area (C–A) fractal model to estimate the weight of each evidence map. In this paper, we used the content region (C–A) fractal model to divide the evidence maps to the threshold of the corresponding dimensions. The P–A plot approach is an objective data-driven approach for evaluating map weights. Using geochemical layer and remote sensing, hydroxyl layers as weight evidence maps are the highlights of this study. We use the P–A method from which we can evaluate the predictive ability of each evidence map with respect to the known ore occurrences. We used the P–A plot for weighting each evidence map and choosing the appropriate threshold for predictor maps in the Luchun area of Yunnan Province, China. The method adopted in this paper can improve the prediction efficiency of ore prospecting.     

Age relationships in supracrustal sequences of the northern part of the Murchison Terrane,Archaean Yilgarn Craton,Western Australia: A combined field and zircon U–Pb study

Mapping carried out in the northern Murchison Terrane of the Archaean Yilgarn Craton, Western Australia, shows that correlation of units between isolated greenstone belts is very difficult and an informal stratigraphic subdivision is proposed where the greenstone sequences have been divided into a number of assemblages. The assemblages may not necessarily be time equivalent throughout the region. The lower units (Assemblages 1–3) consist of ultramafic, mafic and intermediate volcanic rocks deposited without significant breaks in volcanism. Felsic volcanic packages (Assemblage 4) are conformable with underlying units, but are spatially restricted. Discordant units of graphitic sedimentary rocks are developed along major crustal structures (Assemblage 5). SHRIMP and conventional U–Pb study of zircons reveal that felsic volcanic rocks of Assemblage 4 in the Dalgaranga Greenstone Belt were emplaced at 2747 ± 5 Ma, whereas those in the adjacent Meekatharra — Mt Magnet Greenstone Belt range in age from 2762 ± 6 to 2716 ± 4 Ma. The age of emplacement of a differentiated gabbro sill in the Dalgaranga Greenstone Belt at 2719 ± 6 Ma places a maximum age on major folding in the belt. The presence of 2.9–3.0 Ga inherited zircons in some of the felsic volcanic rocks indicates contamination with, or reworking of, underlying 3 Ga sialic crust. This distinguishes the Murchison Terrane from the central parts of the Eastern Goldfields terranes to the south, where there is no evidence for a 3 Ga imprint in zircons from volcanic or granitic rocks, and also from the Narryer Gneiss Terrane to the north and west, which is composed of older gneisses and granitoids. The ca 2.76–2.71 Ga felsic volcanism in the Murchison Terrane is significantly older than 2.71–2.67 Ga felsic volcanism in the Eastern Goldfields lending support to models advocating assemblage of the craton by terrane accretion.     

A combined fluid inclusion and S–Pb isotope study of the Neoproterozoic Pingshui volcanogenic massive sulfide Cu–Zn deposit,Southeast China

The Pingshui Cu–Zn deposit is located in the Jiangshan–Shaoxing fault zone, which marks the Neoproterozoic suture zone between the Yangtze block and Cathaysia block in South China. It contains 0.45 million tons of proven ore reserves with grades of 1.03 wt.% Cu and 1.83 wt.% Zn. This deposit is composed of stratiform, massive sulfide ore bodies, which contain more than 60 vol.% sulfide minerals. These ore bodies are hosted in altered mafic and felsic rocks (spilites and keratophyres) of the bimodal volcanic suite that makes up the Neoproterozoic Pingshui Formation. Metallic minerals include pyrite, chalcopyrite, sphalerite, tennantite, tetrahedrite and magnetite, with minor galena. Gangue minerals are quartz, sericite, chlorite, calcite, gypsum, barite and jasper. Three distinct mineralogical zones are recognized in these massive sulfide ore bodies: a distal zone composed of sphalerite + pyrite + barite (zone I); an intermediate zone characterized by a pyrite + sphalerite + chalcopyrite assemblages (zone II); and a proximal zone containing chalcopyrite + pyrite + magnetite (zone III). A thin, layer of exhalative jaspilite overlies the sulfide ore bodies except in the proximal zone. The volcanic rocks of the Pingshui Formation are all highly altered spilites and keratophyres, but their trace element geochemistry suggests that they were generated by partial melting of the depleted mantle in an island arc setting. Homogenization temperatures of the primary fluid inclusions in quartz from massive sulfide ores are between 217 and 328 °C, and their salinities range from 3.2 to 5.7 wt.% NaCl equivalent. Raman spectroscopy of the fluid inclusions showed that water is the dominant component, with no other volatile components. Fluid inclusion data suggest that the ore-forming fluids were derived from circulating seawater. The δ³⁴S values of pyrite from the massive sulfide ores range from − 3.6‰ to + 3.4‰, indicating that the sulfur was primarily leached from the arc volcanic rocks of the Pingshui Formation. Both pyrite from the massive sulfide ores and plagioclase from the spilites have similar lead isotope compositions, implying that the lead was also derived from the Pingshui Formation. The low lead contents of the massive sulfide ores and the geochemistry of their host rocks are similar to many VMS Cu–Zn deposits in Canada (e.g., Noranda) and thus can be classified as belonging to the bimodal-mafic subtype. The presence of magnetite and the absence of jaspilite and barite at the − 505 m level in the Pingshui deposit suggest that this level is most likely the central zone of the original lateral massive sulfide ore bodies. If this interpretation is correct, the deep part of the Pingshui Cu–Zn deposit may have significant exploration potential.