Petroleum and aqueous inclusions from deeply buried reservoirs: Experimental simulations and consequences for overpressure estimates

Synthetic hydrocarbon and aqueous inclusions have been created in the laboratory batch reactors in order to mimic inclusion formation or re-equilibration in deeply buried reservoirs. Inclusions were synthesized in quartz and calcite using pure water and Mexican dead oil, or n-tetradecane (C₁₄H₃₀), at a temperature and pressure of 150 °C and 1 kbar. One-phase hydrocarbon inclusions are frequently observed at standard laboratory conditions leading to homogenization temperatures between 0 and 60 °C. UV epifluorescence of Mexican oil inclusions is not uniform; blue and green-yellow colored inclusions coexist; however, no clear evidence of variations in fluid chemistry were observed. Homogenization temperatures were recorded and the maxima of T_h plotted on histograms are in good agreement with expected T_h in a range of 6 °C. Broad histograms were reconstructed showing non-symmetrical T_h distributions over a 20 °C temperature range centered on the expected T_h. This histogram broadening is due to the fragility of the fluid inclusions that were created by re-filling of pre-existing microcavities. Such T_h histograms are similar to T_h histograms recorded on natural samples from deeply buried carbonate reservoirs. T_h values lower than those expected were measured for hydrocarbon inclusions in quartz and calcite, and for aqueous inclusions in calcite. However, the results confirm the ability of fluid inclusions containing two immiscible fluids to lead to PT reconstructions, even in overpressured environments.     

Experimental determination of the Raman CH4 symmetric stretching (ν1) band position from 1-650 bar and 0.3-22 °C: Application to fluid inclusion studies

The position of the Raman methane (CH₄) symmetric stretching band (ν₁) over the range 1-650 bar and 0.3-22 °C has been determined using a high-pressure optical cell mounted on a Raman microprobe. Two neon emission lines that closely bracket the CH₄ band were collected simultaneously with each CH₄ spectrum. The peak position was determined after least squares fitting using a summed Gaussian-Lorentzian method, resulting in a precision of ≈±0.02 cm⁻¹ in peak position determination. The CH₄ν₁ band position shifts to lower wave number with increasing pressure. At a given pressure, the band shifts to lower wave number with decreasing temperature, and the magnitude of the temperature shift increases with increasing pressure. The relationship between the Raman CH₄ν₁ band position and temperature and pressure determined here may be used to estimate the internal pressure in natural or synthetic CH₄-bearing fluid inclusions. This information, in turn, may be used to determine the density of pure CH₄ fluid inclusions and the salinity of CH₄-bearing aqueous inclusions.     

Metamorphic fluid flow in the northeastern part of the 3.8-3.7 Ga Isua Greenstone Belt (SW Greenland): A re-evalution of fluid inclusion evidence for early Archean seafloor-hydrothermal systems

Fluid inclusions in quartz globules and quartz veins of a 3.8-3.7 Ga old, well-preserved pillow lava breccia in the northeastern Isua Greenstone Belt (IGB) were studied using microthermometry, Raman spectrometry and SEM Cathodoluminescence Imaging. Petrographic study of the different quartz segregations showed that they were affected by variable recrystallization which controlled their fluid inclusion content. The oldest unaltered fluid inclusions found are present in vein crystals that survived dynamic and static recrystallization. These crystals contain a cogenetic, immiscible assemblage of CO₂-rich (+H₂O, +graphite) and brine-rich (+CO₂, +halite, +carbonate) inclusions. The gas-rich inclusions have molar volumes between 44.8 and 47.5 cm³/mol, while the brine inclusions have a salinity of ∼33 eq. wt% NaCl. Modeling equilibrium immiscibility using volumetric and compositional properties of the endmember fluids indicates that fluid unmixing occurred at or near peak-metamorphic conditions of ∼460 °C and ∼4 kbar. Carbonate and graphite were precipitated cogenetically from the physically separated endmember fluids and were trapped in fluid inclusions.In most quartz crystals, however, recrystallization obliterated such early fluid inclusion assemblages and left graphite and carbonate as solid inclusions in recrystallized grains. Intragranular fluid inclusion trails in the recrystallized grains of breccia cementing and crosscutting quartz veins have CO₂-rich assemblages, with distinctly different molar volumes (either between 43.7 and 47.5 cm³/mol or between 53.5 and 74.1 cm³/mol), and immiscible, halite-saturated H₂O-CO₂-NaCl(-other salt) inclusions. Later intergranular trails have CH₄-H₂ (X_H2 up to ∼0.3) inclusions of variable density (ranging from 48.0 to >105.3 cm³/mol) and metastable H₂O-NaCl(-other salt?) brines (∼28 eq. wt% NaCl). Finally, the youngest fluid inclusion assemblages are found in non-luminescent secondary quartz and contain low-density CH₄ (molar volume > 105.33 cm³/mol) and low-salinity H₂O-NaCl (0.2-3.7 eq. wt% NaCl). These successive fluid inclusion assemblages record a retrograde P-T evolution close to a geothermal gradient of ∼30 °C/km, but also indicate fluid pressure variations and the introduction of highly reducing fluids at ∼200-300 °C and 0.5-2 kbar. The quartz globules in the pillow fragments only contain sporadic CH₄(+H₂) and brine inclusions, corresponding with the late generations present in the cementing and crosscutting veins. We argue that due to the large extent of static recrystallization in quartz globules in the pillow breccia fragments, only these relatively late fluid inclusions have been preserved, and that they do not represent remnants of an early, seafloor-hydrothermal system as was previously proposed.Modeling the oxidation state of the fluids indicates a rock buffered system at peak-metamorphic conditions, but suggests a change towards fluid-graphite disequilibrium and a logf_H2/f_H2O above the Quartz-Fayalite-Magnetite buffer during retrograde evolution. Most likely, this indicates a control on redox conditions and on fluid speciation by ultramafic rocks in the IGB.Finally, this study shows that microscopic solid graphite in recrystallized metamorphic rocks from Isua can be deposited inorganically from a fluid phase, adding to the complexity of processes that formed reduced carbon in the oldest, well-preserved supracrustal rocks on Earth.     

Thermodynamic modeling of binary CH4-H2O fluid inclusions

This work reports the application of thermodynamic models, including equations of state, to binary (salt-free) CH₄-H₂O fluid inclusions. A general method is presented to calculate the compositions of CH₄-H₂O inclusions using the phase volume fractions and dissolution temperatures of CH₄ hydrate. To calculate the homogenization pressures and isolines of the CH₄-H₂O inclusions, an improved activity-fugacity model is developed to predict the vapor-liquid phase equilibrium. The phase equilibrium model can predict methane solubility in the liquid phase and water content in the vapor phase from 273 to 623 K and from 1 to 1000 bar (up to 2000 bar for the liquid phase), within or close to experimental uncertainties. Compared to reliable experimental phase equilibrium data, the average deviation of the water content in the vapor phase and methane solubility in the liquid phase is 4.29% and 3.63%, respectively. In the near-critical region, the predicted composition deviations increase to over 10%. The vapor-liquid phase equilibrium model together with the updated volumetric model of homogenous (single-phase) CH₄-H₂O fluid mixtures (Mao S., Duan Z., Hu J. and Zhang D. (2010) A model for single-phase PVTx properties of CO₂-CH₄-C₂H₆-N₂-H₂O-NaCl fluid mixtures from 273 to 1273 K and from 1 to 5000 bar. Chem. Geol.275, 148-160), is applied to calculate the isolines, homogenization pressures, homogenization volumes, and isochores at specified homogenization temperatures and compositions. Online calculation is on the website: http://www.geochem-model.org/.     

Magnesium-isotope fractionation during low-Mg calcite precipitation in a limestone cave - Field study and experiments

The chemical and isotopic composition of speleothem calcite and particularly that of stalagmites and flowstones is increasingly exploited as an archive of past environmental change in continental settings. Despite intensive research, including modelling and novel approaches, speleothem data remain difficult to interpret. A possible way foreword is to apply a multi-proxy approach including non-conventional isotope systems. For the first time, we here present a complete analytical dataset of magnesium isotopes (δ²⁶Mg) from a monitored cave in NW Germany (Bunker Cave). The data set includes δ²⁶Mg values of loess-derived soil above the cave (−1.0 ± 0.5‰), soil water (−1.2 ± 0.5‰), the carbonate hostrock (−3.8 ± 0.5‰), dripwater in the cave (−1.8 ± 0.2‰), speleothem low-Mg calcite (stalactites, stalagmites; −4.3 ± 0.6‰), cave loam (−0.6 ± 0.1‰) and runoff water (−1.8 ± 0.1‰) in the cave, respectively. Magnesium-isotope fractionation processes during weathering and interaction between soil cover, hostrock and solute-bearing soil water are non-trivial and depend on a number of variables including solution residence times, dissolution rates, adsorption effects and potential neo-formation of solids in the regolith and the carbonate aquifer. Apparent Mg-isotope fractionation between dripwater and speleothem low-Mg calcite is about 1000lnα_Mg-cc-Mg(aq) = −2.4‰. A similar Mg-isotope fractionation (1000lnα_Mg-cc-Mg(aq) ≈ −2.1‰) is obtained by abiogenic precipitation experiments carried out at aqueous Mg/Ca ratios and temperatures close to cave conditions. Accordingly, ²⁶Mg discrimination during low-Mg calcite formation in caves is highly related to inorganic fractionation effects, which may comprise dehydration of Mg²⁺ prior to incorporation into calcite, surface entrapment of light isotopes and reaction kinetics. Relevance of kinetics is supported by a significant negative correlation of Mg-isotope fractionation with the precipitation rate for inorganic precipitation experiments.     

Fluid and gas migration in the North German Basin: fluid inclusion and stable isotope constraints

Fluid inclusions have been studied in minerals infilling fissures (quartz, calcite, fluorite, anhydrite) hosted by Carboniferous and Permian strata from wells in the central and eastern part of the North German Basin in order to decipher the fluid and gas migration related to basin tectonics. The microthermometric data and the results of laser Raman spectroscopy reveal compelling evidence for multiple events of fluid migration. The fluid systems evolved from a H₂O–NaCl±KCl type during early stage of basin subsidence to a H₂O–NaCl–CaCl₂ type during further burial. Locally, fluid inclusions are enriched in K, Cs, Li, B, Rb and other cations indicating intensive fluid–rock interaction of the saline brines with Lower Permian volcanic rocks or sediments. Fluid migration through Carboniferous sediments was often accompanied by the migration of gases. Aqueous fluid inclusions in quartz from fissures in Carboniferous sedimentary rocks are commonly associated with co-genetically trapped CH₄–CO₂ inclusions. P–T conditions estimated, via isochore construction, yield pressure conditions between 620 and 1,650 bar and temperatures between 170 and 300°C during fluid entrapment. The migration of CH₄-rich gases within the Carboniferous rocks can be related to the main stage of basin subsidence and stages of basin uplift. A different situation is recorded in fluid inclusions in fissure minerals hosted by Permian sandstones and carbonates: aqueous fluid inclusions in calcite, quartz, fluorite and anhydrite are always H₂O–NaCl–CaCl₂-rich and show homogenization temperatures between 120 and 180°C. Co-genetically trapped gas inclusions are generally less frequent. When present, they show variable N₂–CH₄ compositions but contain no CO₂. P–T reconstructions indicate low-pressure conditions during fluid entrapment, always below 500 bar. The entrapment of N₂–CH₄ inclusions seems to be related to phases of tectonic uplift during the Upper Cretaceous. A potential source for nitrogen in the inclusions and reservoirs is C_org-rich Carboniferous shales with high nitrogen content. Intensive interaction of brines with Carboniferous or even older shales is proposed from fluid inclusion data (enrichment in Li, Ba, Pb, Zn, Mg) and sulfur isotopic compositions of abundant anhydrite from fissures. The mainly light δ³⁴S values of the fissure anhydrites suggest that sulfate is either derived through oxidation and re-deposition of biogenic sulfur or through mixing of SO₄²⁻-rich formation waters with variable amounts of dissolved biogenic sulfide. An igneous source for nitrogen seems to be unlikely since these rocks have low total nitrogen content and, furthermore, even extremely altered volcanic rocks from the study area do not show a decrease in total nitrogen content.     

The Fe-rich liquidus in the Fe-FeS system from 1 bar to 10 GPa

The composition and evolution of a metallic planetary core is determined by the behavior with pressure of the eutectic and the liquidus on the Fe-rich side of the Fe-FeS eutectic. New experiments at 6 GPa presented here, along with existing experimental data, inform a thermodynamic model for this liquidus from 1 bar to at least 10 GPa. Fe-FeS has a eutectic that becomes more Fe-rich but remains constant in T up to 6 GPa. The 1 bar, 3 GPa, and 6 GPa liquidi all cross at a pivot point at 1640 ± 5 K and FeS_37 ± 0.5. This liquid/crystalline metal equilibrium is T-x-fixed and pressure independent through 6 GPa. Models of the 1 bar through 10 GPa experimental liquidi show that with increasing P there is an increase in the T separation between the liquidus and the crest of the metastable two-liquid solvus. The solvus crest decreases in T with increasing P. The model accurately reproduces all the experimental liquidi from 1 bar to 10 GPa, as well as reproducing the 0-6 GPa pivot point. The 14 GPa experimental liquidus ( [Chen et al., 2008a] and Chen et al., 2008b) deviates sharply from the lower pressure trends indicating that the 0-10 GPa model no longer applies to this 14 GPa data.     

Calibration of the calcite-water oxygen-isotope geothermometer at Devils Hole, Nevada, a natural laboratory

The δ¹⁸O of ground water (−13.54 ± 0.05 ‰) and inorganically precipitated Holocene vein calcite (+14.56 ± 0.03 ‰) from Devils Hole cave #2 in southcentral Nevada yield an oxygen isotopic fractionation factor between calcite and water at 33.7 °C of 1.02849 ± 0.00013 (1000 ln α_{calcite-water} = 28.09 ± 0.13). Using the commonly accepted value of ∂(α_{calcite-water})/∂T of −0.00020 K⁻¹, this corresponds to a 1000 ln α_{calcite-water} value at 25 °C of 29.80, which differs substantially from the current accepted value of 28.3. Use of previously published oxygen isotopic fractionation factors would yield a calcite precipitation temperature in Devils Hole that is 8 °C lower than the measured ground water temperature. Alternatively, previously published fractionation factors would yield a δ¹⁸O of water, from which the calcite precipitated, that is too negative by 1.5 ‰ using a temperature of 33.7 °C. Several lines of evidence indicate that the geochemical environment of Devils Hole has been remarkably constant for at least 10 ka. Accordingly, a re-evaluation of calcite-water oxygen isotopic fractionation factor may be in order.Assuming the Devils Hole oxygen isotopic value of α_{calcite-water} represents thermodynamic equilibrium, many marine carbonates are precipitated with a δ¹⁸O value that is too low, apparently due to a kinetic isotopic fractionation that preferentially enriches ¹⁶O in the solid carbonate over ¹⁸O, feigning oxygen isotopic equilibrium.     

A thermodynamic model for calculating methane solubility, density and gas phase composition of methane-bearing aqueous fluids from 273 to 523 K and from 1 to 2000 bar

A thermodynamic model is presented to calculate methane solubility, liquid phase density and gas phase composition of the H₂O-CH₄ and H₂O-CH₄-NaCl systems from 273 to 523 K (possibly up to 573 K), from 1 to 2000 bar and from 0 to 6 mol kg⁻¹ of NaCl with experimental accuracy. By a more strict theoretical approach and using updated experimental data, this model made substantial improvements over previous models: (1) the accuracy of methane solubility in pure water in the temperature range between 273 and 283 K is increased from about 10% to about 5%, but confirms the accuracy of the Duan model [Duan Z., Moller N., Weare J.H., 1992a. Prediction of methane solubilities in natural waters to high ionic strength from 0 to 250 °C and from 0 to 1600 bar. Geochim. Cosmochim. Acta56, 1451-1460] above 283 K up to 2000 bar; (2) the accuracy of methane solubility in the NaCl aqueous solutions is increased from >12% to about 6% on average from 273 K and 1 bar to 523 K and 2000 bar; (3) this model is able to calculate water content in the gas phase and liquid phase density, which cannot be calculated by previous models; and (4) it covers a wider range of temperature and pressure space. With a simple approach, this model is extended to predict CH₄ solubility in other aqueous salt solutions containing Na⁺, K⁺, Mg²⁺, Ca²⁺, Cl⁻ and , such as seawater and geothermal brines, with excellent accuracy. This model is also able to calculate homogenization pressure of fluid inclusions (CH₄-H₂O-NaCl) and CH₄ solubility in water at gas-liquid-hydrate phase equilibrium. A computer code is developed for this model and can be downloaded from the website: www.geochem-model.org/programs.htm.     

Sorption and catalytic oxidation of Fe(II) at the surface of calcite

The effect of sorption and coprecipitation of Fe(II) with calcite on the kinetics of Fe(II) oxidation was investigated. The interaction of Fe(II) with calcite was studied experimentally in the absence and presence of oxygen. The sorption of Fe(II) on calcite occurred in two distinguishable steps: (a) a rapid adsorption step (seconds-minutes) was followed by (b) a slower incorporation (hours-weeks). The incorporated Fe(II) could not be remobilized by a strong complexing agent (phenanthroline or ferrozine) but the dissolution of the outmost calcite layers with carbonic acid allowed its recovery. Based on results of the latter dissolution experiments, a stoichiometry of 0.4 mol% Fe:Ca and a mixed carbonate layer thickness of 25 nm (after 168 h equilibration) were estimated. Fe(II) sorption on calcite could be successfully described by a surface adsorption and precipitation model (Comans & Middelburg, GCA51 (1987), 2587) and surface complexation modeling (Van Cappellen et al., GCA57 (1993), 3505; Pokrovsky et al., Langmuir16 (2000), 2677). The surface complex model required the consideration of two adsorbed Fe(II) surface species, >CO₃Fe⁺ and >CO₃FeCO₃H⁰. For the formation of the latter species, a stability constant is being suggested. The oxidation kinetics of Fe(II) in the presence of calcite depended on the equilibration time of aqueous Fe(II) with the mineral prior to the introduction of oxygen. If pre-equilibrated for >15 h, the oxidation kinetics was comparable to a calcite-free system (t_1/2 = 145 ± 15 min). Conversely, if Fe(II) was added to an aerated calcite suspension, the rate of oxidation was higher than in the absence of calcite (t_1/2 = 41 ± 1 min and t_1/2 = 100 ± 15 min, respectively). This catalysis was due to the greater reactivity of the adsorbed Fe(II) species, >CO₃FeCO₃H⁰, for which the species specific rate constant was estimated.     

Hydrogen isotope fractionation by Methanothermobacter thermoautotrophicus in coculture and pure culture conditions

We grew a hydrogen-utilizing methanogen, Methanothermobacter thermoautotrophicus strain ΔH, in coculture and pure culture conditions to evaluate the hydrogen isotope fractionation associated with carbonate reduction under low (< several tens of μM; coculture) and high (>6 mM; pure culture) concentrations of H₂ in the headspace. In the cocultures, which were grown at 55 °C with a thermophilic butyrate-oxidizing syntroph, the hydrogen isotopic relationship between methane and water was well represented by the following equation:

δD_CH4=0.725(±0.003)·δD_H2O-275(±3),     

The partitioning behavior of silver in a vapor-brine-rhyolite melt assemblage

The partitioning of silver in a sulfur-free rhyolite melt-vapor-brine assemblage has been quantified at 800 °C, pressures of 100 and 140 MPa and f_O2≈NNO (nickel-nickel oxide). Silver solubility (±2σ) in rhyolite increases 5-fold from 105 ± 21 to 675 ± 98 μg/g as pressure increases from 100 to 140 MPa. Nernst-type partition coefficients describing the mass transfer of silver at 100 MPa between vapor and melt, brine and melt and vapor and brine are 32 ± 30, 1151 ± 238 and 0.026 ± 0.004, respectively. At 140 MPa, values for for vapor and melt, brine and melt, and vapor and brine are 32 ± 10, 413 ± 172 and 0.06 ± 0.03, respectively. Apparent equilibrium constant values (±2σ) describing the exchange of silver and sodium between vapor and melt, , at 100 and 140 MPa are 105 ± 68 and 14 ± 6. The average values (±2σ) for silver and sodium exchange between brine and melt, , at 100 and 140 MPa are 313 ± 288 and 65 ± 12. These data indicate that the mass transfer of silver from rhyolite melt to an exsolved volatile phase(s) is enhanced at 100 MPa relative to 140 MPa, suggesting that decompression increases the silver ore-generative potential of an evolving silicate magma. Model calculations using the new data suggest that the evolution of low-density, aqueous fluid (i.e., vapor) may be responsible for the the silver tonnage of many porphyry-type and perhaps epithermal-type ore deposits. For example, Halter et al. (Halter W. E., Pettke T. and Heinrich C. A. (2002) The origin of Cu/Au ratios in porphyry-type ore deposits. Science296, 1842-1844) used detailed silicate and sulfide melt inclusion and vapor and brine fluid inclusions analyses to estimate a melt volume on the order of 15 km³ to satisfy the copper budget at the Bajo de la Alumbrera copper-, gold-, silver-ore deposit. Using their melt volume estimate with the data presented here, model calculations for a 15-km³ felsic melt, saturated with pyrrhotite and magnetite, suggest that a low-salinity magmatic vapor may scavenge on the order of 7 × 10¹² g of silver from the melt. This quantity of silver exceeds the discovered 2 × 10⁹ g of Ag at Alumbrera. Calculated tonnages for numerous other deposits yield similar results. The excess silver in the vapor, remaining after porphyry formation, is then available to precipitate at lower PTconditions in the stratigraphically higher epithermal environment. These data suggest that silver, and perhaps other ore metals, in the porphyry-epithermal continuum may be derived solely from the time-integrated flux of dominantly low-salinity vapor exsolved from a series of sequential magma batches.     

Synthetic fluid inclusions XVIII: Experimental determination of the PVTX properties of H2O-CH4 to 500 °C, 3 kbar and XCH4 ?4 mol.%

The pressure-volume-temperature-composition (PVTX) properties of H₂O-CH₄ were determined from the bubble point curve to 500 °C and 3 kbar for compositions ?4 mol.% CH₄ using the synthetic fluid inclusion technique. H₂O-CH₄ inclusions were produced by loading known amounts of Al₃C₄ and H₂O into platinum capsules along with pre-fractured and inclusion-free quartz cores. During heating the Al₃C₄ and H₂O react to produce CH₄, and the H₂O-CH₄ homogeneous mixture was trapped as inclusions during fracture healing at elevated temperature and pressure. The composition of the fluid in the inclusion was confirmed using the weight loss technique after the experiment and by Raman spectroscopic analysis of the inclusions.Homogenization temperatures of the inclusions were determined and the results were used to construct iso-T_h lines, defined as a line connecting the formation temperature and pressure with the homogenization temperature and pressure. The pressure in the inclusion at the homogenization temperature was calculated from the Duan equation of state (EOS). The slope (ΔP/ΔT) of each iso-T_h line was calculated and the results fitted to a polynomial equation using step-wise multiple regression analysis to estimate the slope of the iso-T_h line as a function of the homogenization temperature and composition according to:

(ΔP/ΔT)=a+b·m+c·m⁴+d·(T_h)²+e·m·T_h+f·m·(T_h)⁴,     

Determination of methane concentrations in water in equilibrium with sI methane hydrate in the absence of a vapor phase by in situ Raman spectroscopy

Most submarine gas hydrates are located within the two-phase equilibrium region of hydrate and interstitial water with pressures (P) ranging from 8 to 60 MPa and temperatures (T) from 275 to 293 K. However, current measurements of solubilities of methane in equilibrium with hydrate in the absence of a vapor phase are limited below 20 MPa and 283.15 K, and the differences among these data are up to 30%. When these data were extrapolated to other P-T conditions, it leads to large and poorly known uncertainties. In this study, in situ Raman spectroscopy was used to measure methane concentrations in pure water in equilibrium with sI (structure one) methane hydrate, in the absence of a vapor phase, at temperatures from 276.6 to 294.6 (±0.3) K and pressures at 10, 20, 30 and 40 (±0.4%) MPa. The relationship among concentration of methane in water in equilibrium with hydrate, in mole fraction [X(CH₄)], the temperature in K, and pressure in MPa was derived as: X(CH₄) = exp [11.0464 + 0.023267 P − (4886.0 + 8.0158 P)/T]. Both the standard enthalpy and entropy of hydrate dissolution at the studied T-P conditions increase slightly with increasing pressure, ranging from 41.29 to 43.29 kJ/mol and from 0.1272 to 0.1330 kJ/K · mol, respectively. When compared with traditional sampling and analytical methods, the advantages of our method include: (1) the use of in situ Raman signals for methane concentration measurements eliminates possible uncertainty caused by sampling and ex situ analysis, (2) it is simple and efficient, and (3) high-pressure data can be obtained safely.     

Ore-forming fluids associated with granite-hosted gold mineralization at the Sanshandao deposit, Jiaodong gold province, China

The Sanshandao gold deposit, with total resources of more than 60 t of gold, is located in the Jiaodong gold province, the most important gold province of China. The deposit is a typical highly fractured and altered, disseminated gold system, with high-grade, quartz-sulphide vein/veinlet stockworks that cut Mesozoic granodiorite. There are four stages of veins that developed in the following sequence: (1) quartz-K-feldspar-sericite; (2) quartz-pyrite±arsenopyrite; (3) quartz-base metal sulfide; and (4) quartz-carbonate. Fluid inclusions in quartz and calcite in vein/veinlet stockworks contain C-O-H fluids of three main types. The first type consists of dilute CO₂–H₂O fluids coeval with the early vein stage. Molar volumes of these CO₂–H₂O fluid inclusions, ranging from 50–60 cm³/mol, yield estimated minimum trapping pressures of 3 kbar. Homogenization temperatures, obtained mainly from CO₂–H₂O inclusions with lower CO₂ concentration, range from 267–375 °C. The second inclusion type, with a CO₂–H₂O±CH₄ composition, was trapped during the main mineralizing stages. These fluids may reflect the CO₂–H₂O fluids that were modified by fluid/rock reactions with altered wallrocks. Isochores for CO₂-H₂O±CH₄ inclusions, with homogenization temperatures ranging from 204–325 °C and molar volumes from 55 to 70 cm³/mol, provide an estimated minimum trapping pressure of 1.2 kbar. The third inclusion type, aqueous inclusions, trapped in cross-cutting microfractures in quartz and randomly in calcite, are post-mineralization, and have homogenization temperatures between 143–228 °C and salinities from 0.71–7.86 wt% NaCl equiv. Stable isotope data show that the metamorphic fluid contribution is minimal and that ore fluids are of magmatic origin, most likely sourced from 120–126 Ma mafic to intermediate dikes. This is consistent with the carbonic nature of the fluid, and the cross-cutting nature of those deposits relative to the host Mesozoic granitoid.Editorial handling: R.J. Goldfarb     

Majoritic garnet: A new approach to pressure estimation of shock events in meteorites and the encapsulation of sub-lithospheric inclusions in diamond

A means for estimating pressures in natural samples based on both the coupled substitution (Na+)^[1+] (Ti + ^[VI]Si)^[4+] = (M)^[2+] (Al + Cr)^[3+], and the classic pyroxene-stoichiometry majorite-substitution into garnet at high-pressure, is derived for garnets with majoritic chemistry. The technique is based on a compilation of experimental data for different bulk compositions. It is compositionally and thermally robust and can be used to estimate pressures experienced by natural materials during formation of majoritic garnet. In addition, it can be used either retrospectively, or in new experimental studies to establish the pressures of crystallization of reaction products, and determine if disequilibrium is recorded by the chemistries of majoritic garnets. Pressures are calculated based on majoritic chemistries in chondritic meteorites and diamond inclusions. Majoritic garnets associated with Mg perovskite in shocked L chondrites (n = 4) yield uniform pressures of 23.8 ± 0.2 GPa that are slightly higher than pressures recorded by majoritic garnet in shock-derived melt veins in L chondrites (22.4 ± 0.6 GPa; n = 5). Similar pressures are also exhibited by shock-derived majoritic garnets in H chondrites (22.2 ± 1.1 GPa; n = 3). Diamond inclusions with eclogitic and peridotitic majoritic garnet chemistries exhibit mean pressures of 10.7 ± 2.7 GPa (n = 30) and 8.3 ± 1.6 GPa (n = 15) respectively, consistent with a sub-lithospheric origin. However, pressures defined by majoritic diamond inclusions from Jagersfontein (22.3 ± 0.8 GPa and 16.9 ± 1 GPa), Monastery (15.7 ± 7 GPa) and Kankan (15.5 ± 0.2 GPa) show that these inclusions originated from the mantle transition zone. Thus, this new single-phase method for pressure estimation has unmatched potential to map the depth of formation of garnets with majoritic chemistries that occur as diamond inclusions in all parageneses except those that include Ca silicate perovskite. The derived pressures confirm the sub-lithospheric origin of eclogitic majoritic diamond inclusions, and thus provide a more comprehensive picture of the important role of storage of oceanic lithosphere in the transition zone.     

Coupled phase and aqueous species equilibrium of the H2O-CO2-NaCl-CaCO3 system from 0 to 250 °C, 1 to 1000 bar with NaCl concentrations up to saturation of halite

A model is developed for the calculation of coupled phase and aqueous species equilibrium in the H₂O-CO₂-NaCl-CaCO₃ system from 0 to 250 °C, 1 to 1000 bar with NaCl concentrations up to saturation of halite. The vapor-liquid-solid (calcite, halite) equilibrium together with the chemical equilibrium of H⁺, Na⁺, Ca²⁺, , Ca(OH)⁺, OH⁻, Cl⁻, , , CO_2(aq) and CaCO_3(aq) in the aqueous liquid phase as a function of temperature, pressure, NaCl concentrations, CO_2(aq) concentrations can be calculated, with accuracy close to those of experiments in the stated T-P-m range, hence calcite solubility, CO₂ gas solubility, alkalinity and pH values can be accurately calculated. The merit and advantage of this model is its predictability, the model was generally not constructed by fitting experimental data.One of the focuses of this study is to predict calcite solubility, with accuracy consistent with the works in previous experimental studies. The resulted model reproduces the following: (1) as temperature increases, the calcite solubility decreases. For example, when temperature increases from 273 to 373 K, calcite solubility decreases by about 50%; (2) with the increase of pressure, calcite solubility increases. For example, at 373 K changing pressure from 10 to 500 bar may increase calcite solubility by as much as 30%; (3) dissolved CO₂ can increase calcite solubility substantially; (4) increasing concentration of NaCl up to 2 m will increase calcite solubility, but further increasing NaCl solubility beyond 2 m will decrease its solubility.The functionality of pH value, alkalinity, CO₂ gas solubility, and the concentrations of many aqueous species with temperature, pressure and NaCl_(aq) concentrations can be found from the application of this model. Online calculation is made available on www.geochem-model.org/models/h2o_co2_nacl_caco3/calc.php.     

长岭凹陷深层天然气藏中CO_2及CH_4充注次序的流体包裹体证据

对长岭凹陷深层天然气藏储层——营城组火山岩中发育的流体包裹体进行了详细研究,结果表明在火山岩发育的石英、方解石细网脉中均存在较多的碳质流体包裹体,单个包裹体激光拉曼光谱分析结果表明其主要为CO2及CH4两种类型的碳质包裹体。其中方解石细网脉体中发育有原生及次生CH4包裹体,而含CO2包裹体多以原生包裹体产于石英细网脉中。很多含CO2包裹体的石英细脉中发现了含CH4包裹体的方解石脉体的角砾,这就表明石英细脉形成晚于方解石细脉。营城组火山岩储层中CO2及CH4包裹体的产状特征研究表明,松辽盆地深层天然气藏的形成系火山岩成岩后CO2及CH4等气体不同期次充注的结果,CH4气的充注时间早于CO2气,火山岩中发育的原生孔隙及次生裂隙为上述气体的充注和聚集提供了重要通道。     

Thermal history of the unsaturated zone at Yucca Mountain,Nevada, USA

Secondary calcite, silica and minor amounts of fluorite deposited in fractures and cavities record the chemistry, temperatures, and timing of past fluid movement in the unsaturated zone at Yucca Mountain, Nevada, the proposed site of a high-level radioactive waste repository. The distribution and geochemistry of these deposits are consistent with low-temperature precipitation from meteoric waters that infiltrated at the surface and percolated down through the unsaturated zone. However, the discovery of fluid inclusions in calcite with homogenization temperatures (T_h) up to ∼80 °C was construed by some scientists as strong evidence for hydrothermal deposition. This paper reports the results of investigations to test the hypothesis of hydrothermal deposition and to determine the temperature and timing of secondary mineral deposition. Mineral precipitation temperatures in the unsaturated zone are estimated from calcite- and fluorite-hosted fluid inclusions and calcite δ¹⁸O values, and depositional timing is constrained by the ²⁰⁷Pb/²³⁵U ages of chalcedony or opal in the deposits. Fluid inclusion T_h from 50 samples of calcite and four samples of fluorite range from ∼35 to ∼90 °C. Calcite δ¹⁸O values range from ∼0 to ∼22‰ (SMOW) but most fall between 12 and 20‰. The highest T_h and the lowest δ¹⁸O values are found in the older calcite. Calcite T_h and δ¹⁸O values indicate that most calcite precipitated from water with δ¹⁸O values between −13 and −7‰, similar to modern meteoric waters.     

Retrograde fluid evolution at the Rampura Agucha Pb-Zn-(Ag) deposit, Rajasthan, India

Fluid inclusions in mineralized graphite-sillimanite-mica schist from the Rampura-Agucha Pb-Zn-(Ag) deposit, Rajasthan, northwest India, have been investigated by microthermometry and Raman microspectrometry. Three different main types of fluid inclusions in quartz can be distinguished: (1) gaseous (CO₂, partially mixed with CH₄-N₂), (2) low salinity aqueous inclusions (0–8 eq. wt% NaCl) and (3) high salinity aqueous inclusions (NaCl ± MgCl₂-CaCl₂). Low density CO₂-rich and low salinity H₂O inclusions are contemporaneous and occur, together with CH₄-N₂ inclusions, in close association with sulfide mineral inclusions. This indicates immiscibility between the gaseous and aqueous phase and participation of these fluids during the deposition or remobilization of the ore, which occurred over a wide P (1220 to 200 bar) and T (450 to 250 °C). Raman spectra of graphite indicate upper greenschist-facies metamorphic conditions, although host rocks have been metamorphosed at upper amphibolite-facies metamorphic conditions. This indicates that graphite re-equilibrated with the CO₂-rich phase during retrograde metamorphism.