Chemical and isotopic constraints on water/rock interactions at the Lost City hydrothermal field, 30°N Mid-Atlantic Ridge

Low temperature vent fluids (<91 °C) issuing from the ultramafic-hosted hydrothermal system at Lost City, 30°N Mid-Atlantic Ridge, are enriched in dissolved volatiles (H₂,CH₄) while attaining elevated pH values, indicative of the serpentization processes that govern water/rock interactions deep in the oceanic crust. Here, we present a series of theoretical models to evaluate the extent of hydrothermal alteration and assess the effect of cooling on the systematics of pH-controlled B aqueous species. Peridotite-seawater equilibria calculations indicate that the mineral assemblage composed of diopside, brucite and chrysotile likely dictates fluid pH at moderate temperature serpentinization processes (<300 °C), by imposing constraints on the aCa⁺⁺/a²H⁺ ratios and the activity of dissolved SiO₂. Based on Sr abundances and the ⁸⁷Sr/⁸⁶Sr isotope ratios of vent fluids reported from Lost City, estimated water/rock mass ratios (w/r = 2-4) are consistent with published models involving dissolved CO₂ and alkane concentrations. Combining the reported δ¹⁸O values of vent fluids (0.7‰) with such w/r mass ratios, allows us to bracket subseafloor reaction temperatures in the vicinity of 250 °C. These estimates are in agreement with previous theoretical studies supporting extensive conductive heat loss within the upflow zones. Experimental studies on peridotite-seawater alteration suggest that fluid pH increases during cooling which then rapidly enhances boron removal from solution and incorporation into secondary phases, providing an explanation for the highly depleted dissolved boron concentrations measured in the low temperature but alkaline Lost City vent fluids. Finally, to account for the depleted ¹¹B composition (δ¹¹B ∼25-30‰) of vent fluids relative to seawater, isotopic fractionation between tetrahedrally coordinated aqueous boron species with BO₃-bearing mineral sites (e.g. in calcite, brucite) is proposed.     

Mineral control of arsenic content in thermal waters from volcano-hosted hydrothermal systems: Insights from island of Ischia and Phlegrean Fields (Campanian Volcanic Province, Italy)

This paper documents arsenic concentrations in 157 groundwater samples from the island of Ischia and the Phlegrean Fields, two of the most active volcano-hosted hydrothermal systems from the Campanian Volcanic Province (Southern Italy), in an attempt to identify the environmental conditions and mineral-solution reactions governing arsenic aqueous cycling. On Ischia and in the Phlegrean Fields, groundwaters range in composition from NaCl brines, which we interpret as the surface discharge of deep reservoir fluids, to shallow-depth circulating fluids, the latter ranging from acid-sulphate steam-heated to hypothermal, cold, bicarbonate groundwaters. Arsenic concentrations range from 1.6 to 6900 μg·l^− 1 and from 2.6 to 3800 μg·l^− 1 in the Phlegrean Fields and on Ischia, respectively. They increase with increasing water temperature and chlorine contents, and in the sequence bicarbonate groundwaters < steam-heated groundwaters < NaCl brines. According to thermochemical modeling, we propose that high As concentrations in NaCl brines form after prolonged water-rock interactions at reservoir T, fO₂ and fH₂S conditions, and under the buffering action of an arsenopyrite + pyrite + pyrrhotite rock assemblage. On their ascent toward the surface, NaCl brines become diluted by As-depleted meteoric-derived bicarbonate groundwaters, giving rise to hybrid water types with intermediate to low As contents. Steam-heated groundwaters give their intermediate to high As concentrations to extensive rock leaching promoted by interaction with As-bearing hydrothermal steam.     

Catabolic and anabolic energy for chemolithoautotrophs in deep-sea hydrothermal systems hosted in different rock types

Active deep-sea hydrothermal vents are hosted by a range of different rock types, including basalt, peridotite, and felsic rocks. The associated hydrothermal fluids exhibit substantial chemical variability, which is largely attributable to compositional differences among the underlying host rocks. Numerical models were used to evaluate the energetics of seven inorganic redox reactions (potential catabolisms of chemolithoautotrophs) and numerous biomolecule synthesis reactions (anabolism) in a representative sampling of these systems, where chemical gradients are established by mixing hydrothermal fluid with seawater. The wide ranging fluid compositions dictate demonstrable differences in Gibbs energies (ΔG_r) of these catabolic and anabolic reactions in three peridotite-hosted, six basalt-hosted, one troctolite-basalt hybrid, and two felsic rock-hosted systems. In peridotite-hosted systems at low to moderate temperatures (<∼45 °C) and high seawater:hydrothermal fluid (SW:HF) mixing ratios (>10), hydrogen oxidation yields the most catabolic energy, but the oxidation of methane, ferrous iron, and sulfide can also be moderately exergonic. At higher temperatures, and consequent SW:HF mixing ratios <10, anaerobic processes dominate the energy landscape; sulfate reduction and methanogenesis are more exergonic than any of the aerobic respiration reactions. By comparison, in the basalt-hosted and felsic rock-hosted systems, sulfide oxidation was the predominant catabolic energy source at all temperatures (and SW:HF ratios) considered. The energetics of catabolism at the troctolite-basalt hybrid system were intermediate to these extremes. Reaction energetics for anabolism in chemolithoautotrophs—represented here by the synthesis of amino acids, nucleotides, fatty acids, saccharides, and amines—were generally most favorable at moderate temperatures (22-32 °C) and corresponding SW:HF mixing ratios (∼15). In peridotite-hosted and the troctolite-basalt hybrid systems, ΔG_r for primary biomass synthesis yielded up to ∼900 J per g dry cell mass. The energetics of anabolism in basalt- and felsic rock-hosted systems were far less favorable. The results suggest that in peridotite-hosted (and troctolite-basalt hybrid) systems, compared with their basalt (and felsic rock) counterparts, microbial catabolic strategies—and consequently variations in microbial phylotypes—may be far more diverse and some biomass synthesis may yield energy rather than imposing a high energetic cost.     

Rare earth element abundances in hydrothermal fluids from the Manus Basin, Papua New Guinea: Indicators of sub-seafloor hydrothermal processes in back-arc basins

Rare earth element (REE) concentrations are reported for a large suite of seafloor vent fluids from four hydrothermal systems in the Manus back-arc basin (Vienna Woods, PACMANUS, DESMOS and SuSu Knolls vent areas). Sampled vent fluids show a wide range of absolute REE concentrations and chondrite-normalized (REE_N) distribution patterns (La_N/Sm_N ∼ 0.6-11; La_N/Yb_N ∼ 0.6 - 71; ). REE_N distribution patterns in different vent fluids range from light-REE enriched, to mid- and heavy-REE enriched, to flat, and have a range of positive Eu-anomalies. This heterogeneity contrasts markedly with relatively uniform REE_N distribution patterns of mid-ocean ridge hydrothermal fluids. In Manus Basin fluids, aqueous REE compositions do not inherit directly or show a clear relationship with the REE compositions of primary crustal rocks with which hydrothermal fluids interact. These results suggest that the REEs are less sensitive indicators of primary crustal rock composition despite crustal rocks being the dominant source of REEs in submarine hydrothermal fluids. In contrast, differences in aqueous REE compositions are consistently correlated with differences in fluid pH and ligand (chloride, fluoride and sulfate) concentrations. Our results suggest that the REEs can be used as an indicator of the type of magmatic acid volatile (i.e., presence of HF, SO₂) degassing in submarine hydrothermal systems. Additional fluid data suggest that near-seafloor mixing between high-temperature hydrothermal fluid and locally entrained seawater at many vent areas in the Manus Basin causes anhydrite precipitation. Anhydrite effectively incorporates REE and likely affects measured fluid REE concentrations, but does not affect their relative distributions.     

Determination of fluid/melt partition coefficients by LA-ICPMS analysis of co-existing fluid and silicate melt inclusions: Controls on element partitioning

Analyses of co-existing silicate melt and fluid inclusions, entrapped in quartz crystals in volatile saturated magmatic systems, allowed direct quantitative determination of fluid/melt partition coefficients. Investigations of various granitic systems (peralkaline to peraluminous in composition, log fO₂ = NNO−1.7 to NNO+4.5) exsolving fluids with various chlorinities (1-14 mol/kg) allowed us to assess the effect of these variables on the fluid/melt partition coefficients (D). Partition coefficients for Pb, Zn, Ag and Fe show a nearly linear increase with the chlorinity of these fluid (D_Pb ∼ 6 ∗ m_Cl, D_Zn ∼ 8 ∗ m_Cl, D_Ag ∼ 4 ∗ m_Cl, D_Fe ∼ 1.4 ∗ m_Cl, where m_Cl is the molinity of Cl). This suggests that these metals are dissolved primarily as Cl-complexes and neither oxygen fugacity nor the composition of the melt affects significantly their fluid/melt partitioning. By contrast, partition coefficients for Mo, B, As, Sb and Bi are highest in low salinity (1-2 mol/kg Cl) fluids with maximum values of D_Mo ∼ 20, D_B ∼ 15, D_As ∼ 13, D_Sb ∼ 8, D_Bi ∼ 15 indicating dissolution as non-chloride (e.g., hydroxy) complexes. Fluid/melt partition coefficients of copper are highly variable, but highest between vapor like fluids and silicate melt (D_Cu ? 2700), indicating an important role for ligands other than Cl. Partition coefficients for W generally increase with increasing chlorinity, but are exceptionally low in some of the studied brines which may indicate an effect of other parameters. Fluid/melt partition coefficients of Sn show a high variability but likely increase with the chlorinity of the fluid (D_Sn = 0.3-42, D_W = 0.8-60), and decrease with decreasing oxygen fugacity or melt peraluminosity.     

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.     

The role of hydrologic regimes on dissolved organic carbon composition in an agricultural watershed

Willow Slough, a seasonally irrigated agricultural watershed in the Sacramento River valley, California, was sampled weekly in 2006 in order to investigate seasonal concentrations and compositions of dissolved organic carbon (DOC). Average DOC concentrations nearly doubled from winter baseflow (2.75 mg L⁻¹) to summer irrigation (5.14 mg L⁻¹), while a concomitant increase in carbon-normalized vanillyl phenols (0.11 mg 100 mg OC⁻¹ increasing to 0.31 mg 100 mg OC⁻¹, on average) indicates that this additional carbon is likely vascular plant-derived. A strong linear relationship between lignin concentration and total suspended sediments (r² = 0.79) demonstrates that agricultural management practices that mobilize sediments will likely have a direct and significant impact on DOC composition. The original source of vascular plant-derived DOC to Willow Slough appears to be the same throughout the year as evidenced by similar syringyl to vanillyl and cinnamyl to vanillyl ratios. However, differing diagenetic pathways during winter baseflow as compared to the rest of the year are evident in acid to aldehyde ratios of both vanillyl and syringyl phenols. The chromophoric dissolved organic matter (CDOM) absorption coefficient at 350 nm showed a strong correlation with lignin concentration (r² = 0.83). Other CDOM measurements related to aromaticity and molecular weight also showed correlations with carbon-normalized yields (e.g. specific UV absorbance at 254 nm (r² = 0.57) and spectral slope (r² = 0.54)). Our overall findings suggest that irrigated agricultural watersheds like Willow Slough can potentially have a significant impact on mainstem DOC concentration and composition when scaled to the entire watershed of the main tributary.     

The use of vein fluorite as probe for paleofluid REE and Sr-Nd isotope geochemistry: The Santa Catarina Fluorite District, Southern Brazil

Fluorite can be used as a probe for the source of Sr and REE, as well as for the Sr and Nd isotope systematics of mineralizing solutions, allowing characterization of the composition, oxidation state and sources of the fluids. The ⁸⁷Sr / ⁸⁶Sr ratios in vein fluorite from the Santa Catarina Fluorite District, southern Brazil, are low (0.720 to 0.745) relative to those of the majority of host granites at the time of mineralization (90 Ma), but are similar to those of less abundant and less evolved Sr- and Ca-rich granites and plagioclases of the heterogeneous Pedras Grandes granite association. Major contributions of Sr from the unradiogenic Parana Basin rocks (⁸⁷Sr / ⁸⁶Sr_{90 Ma} = 0.705 to 0.718) are unlikely, considering the radiogenic character of the lower ⁸⁷Sr / ⁸⁶Sr end-member in fluorite mixing lines. Estimated fluorite fluid partition coefficients (K_d^Sr-Ca = 0.019 and D^Sr ≈ 600) indicate a Sr / Ca ratio in the fluorite-forming solution of 0.012, and Sr contents of 0.05 to 0.25 ppm, which are similar to those of present-day granitic geothermal waters. Initial Nd isotopic compositions of the vein fluorites (0.5120 to 0.512) are similar to those of the Pedras Grandes granites. The ¹⁴³Nd / ¹⁴⁴Nd_{90 Ma} of the evolved granites of the Tabuleiro granite association, their accessory fluorites and the Parana Basin rocks are considerably more radiogenic (0.5120 to 0.5127) and these are thus considered to be unlikely sources of the fluids. The REE patterns of vein fluorites, normalized to upper continental crust, show a range of LREE-depleted patterns, with highly variable positive and negative Eu anomalies. The host Pedras Grandes granites show flat to slightly depleted UCC normalized LREE patterns with strong negative Eu anomalies. Depletion of the LREE in fluorites resulted from the mobility of HREE fluoride complexes during fluid migration. A REE fractionation model based on ionic potential ratios indicates that Eu³⁺ was stable during fluid migration and fluorite precipitation. The coexistence of pyrite and Eu³⁺ in the mineralizing fluids is consistent with low pH and oxygen fugacities near the hematite-magnetite buffer.     

Geochemistry of hydrothermal fluids from the PACMANUS, Northeast Pual and Vienna Woods hydrothermal fields, Manus Basin, Papua New Guinea

Processes controlling the composition of seafloor hydrothermal fluids in silicic back-arc or near-arc crustal settings remain poorly constrained despite growing evidence for extensive magmatic-hydrothermal activity in such environments. We conducted a survey of vent fluid compositions from two contrasting sites in the Manus back-arc basin, Papua New Guinea, to examine the influence of variations in host rock composition and magmatic inputs (both a function of arc proximity) on hydrothermal fluid chemistry. Fluid samples were collected from felsic-hosted hydrothermal vent fields located on Pual Ridge (PACMANUS and Northeast (NE) Pual) near the active New Britain Arc and a basalt-hosted vent field (Vienna Woods) located farther from the arc on the Manus Spreading Center. Vienna Woods fluids were characterized by relatively uniform endmember temperatures (273-285 °C) and major element compositions, low dissolved CO₂ concentrations (4.4 mmol/kg) and high measured pH (4.2-4.9 at 25 °C). Temperatures and compositions were highly variable at PACMANUS/NE Pual and a large, newly discovered vent area (Fenway) was observed to be vigorously venting boiling (358 °C) fluid. All PACMANUS fluids are characterized by negative δD_H2O values, in contrast to positive values at Vienna Woods, suggesting substantial magmatic water input to circulating fluids at Pual Ridge. Low measured pH (25 °C) values (∼2.6-2.7), high endmember CO₂ (up to 274 mmol/kg) and negative δ³⁴S_H2S values (down to −2.7‰) in some vent fluids are also consistent with degassing of acid-volatile species from evolved magma. Dissolved CO₂ at PACMANUS is more enriched in ¹³C (−4.1‰ to −2.3‰) than Vienna Woods (−5.2‰ to −5.7‰), suggesting a contribution of slab-derived carbon. The mobile elements (e.g. Li, K, Rb, Cs and B) are also greatly enriched in PACMANUS fluids reflecting increased abundances in the crust there relative to the Manus Spreading Center. Variations in alkali and dissolved gas abundances with Cl at PACMANUS and NE Pual suggest that phase separation has affected fluid chemistry despite the low temperatures of many vents. In further contrast to Vienna Woods, substantial modification of PACMANUS/NE Pual fluids has taken place as a result of seawater ingress into the upflow zone. Consistently high measured Mg concentrations as well as trends of increasingly non-conservative SO₄ behavior, decreasing endmember Ca/Cl and Sr/Cl ratios with increased Mg indicate extensive subsurface anhydrite deposition is occurring as a result of subsurface seawater entrainment. Decreased pH and endmember Fe/Mn ratios in higher Mg fluids indicate that the associated mixing/cooling gives rise to sulfide deposition and secondary acidity production. Several low temperature (?80 °C) fluids at PACMANUS/NE Pual also show evidence for anhydrite dissolution and water-rock interaction (fixation of B) subsequent to seawater entrainment. Hence, the evolution of fluid compositions at Pual Ridge reflects the cumulative effects of water/rock interaction, admixing and reaction of fluids exsolved from silicic magma, phase separation/segregation and seawater ingress into upflow zones.     

Fluid-rock interactions in a geothermal Rotliegend/Permo-Carboniferous reservoir (North German Basin)

Comprehensive data on the chemical composition of reservoir rocks and geothermal brines from the geothermal well doublet Groβ Schönebeck (North German Basin) drilled into a Rotliegend sedimentary and Permo-Carboniferous volcanic rock reservoir were sampled over the past years. They were characterized with respect to their major and minor elemental composition including various isotope ratios. The study considered the impact of drilling and reservoir operations on fluid composition and aimed at determining fluid–rock interactions to gain information on fluid origin and hydraulic pathways.The highly saline fluids (up to 265 g/L TDS) show δ ¹⁸O and δD of water (2.7–5.6 and −3.1–15, respectively) as well as δ ³⁴S of sulfate (3.6–5), and ⁸⁷Sr/⁸⁶Sr ratios (0.715–0.716) that resemble Rotliegend brines from an area located around 200 km in the west (the Altmark). Halogen ratios indicated that brines developed predominantly by evaporation of meteoric water (primary brine) together with halite dissolution brine (secondary brine). Indication for mixing with Zechstein brine or with younger meteoric water was not found.No geochemical distinction was possible between fluids deriving from different rock formations (dacites or sedimentary rocks, respectively). This is due to the evolution of the sediments from the effusive rocks resulting in a similar mineralogical and chemical composition and due to a hydraulic connectivity between the two types of rock. This connection existed probably already before reservoir stimulation as indicated by a set of faults identified in the area that could connect the Rotliegend formation with both, the volcanic rocks and the lower units of the Zechstein. Additional geochemical indication for a hydraulic connectivity is given by (1) the very high heavy metal contents (mainly Cu and Pb) in fluids and scaling that derive from the volcanic rocks and were that were also found in increased amounts up at the Zechstein border (Kupferschiefer formation). (2) The ⁸⁷Sr/⁸⁶Sr isotope ratios of fluid samples correspond to the ratios determined for the sedimentary rocks indicating that initially the fluids developed in the sedimentary rocks and circulated later, when faults structures were created by tectonic events into the volcanic rocks.     

Isotopic evidence (Sr/Sr, δLi) for alteration of the oceanic crust at deep-rooted mud volcanoes in the Gulf of Cadiz, NE Atlantic Ocean

The chemical and isotopic composition of pore fluids is presented for five deep-rooted mud volcanoes aligned on a transect across the Gulf of Cadiz continental margin at water depths between 350 and 3860 m. Generally decreasing interstitial Li concentrations and ⁸⁷Sr/⁸⁶Sr ratios with increasing distance from shore are attributed to systematically changing fluid sources across the continental margin. Although highest Li concentrations at the near-shore mud volcanoes coincide with high salinities derived from dissolution of halite and late-stage evaporites, clayey, terrigenous sediments are identified as the ultimate Li source to all pore fluids investigated. Light δ⁷Li values, partly close to those of hydrothermal vent fluids (δ⁷Li: +11.9‰), indicate that Li has been mobilized during high-temperature fluid/sediment or fluid/rock interactions in the deep sub-surface. Intense leaching of terrigenous clay has led to radiogenic ⁸⁷Sr/⁸⁶Sr ratios (∼0.7106) in pore fluids of the near-shore mud volcanoes. In contrast, non-radiogenic ⁸⁷Sr/⁸⁶Sr ratios (∼0.7075) at the distal locations are attributed to admixing of a basement-derived fluid component, carrying an isotopic signature from interaction with the basaltic crust. This inference is substantiated by temperature constraints from Li isotope equilibrium calculations suggesting exchange processes at particularly high temperatures (>200 °C) for the least radiogenic pore fluids of the most distal location.Advective pore fluids in the off-shore reaches of the Gulf of Cadiz are influenced by successive exchange processes with both oceanic crust and terrigenous, fine-grained sediments, resulting in a chemical and isotopic signature similar to that of fluids in near-shore ridge flank hydrothermal systems. This suggests that deep-rooted mud volcanoes in the Gulf of Cadiz represent a fluid pathway intermediate between mid-ocean ridge hydrothermal vent and shallow, marginal cold seep. Due to the thicker sediment coverage and slower fluid advection rates, the overall geochemical signature is shifted towards the sediment-diagenetic signal compared to ridge flank hydrothermal environments.     

High integrated fluid/rock ratios during metamorphism at Naxos: evidence from carbon isotopes of calcite in schists and fluid inclusions

Calcite in schists of the metamorphic complex at Naxos is depleted both in ¹³C and in ¹⁸O with respect to massive marbles. This effect is attributed to isotope exchange with circulating CO₂-rich fluids, which had an >0.5 according to fluid inclusions. The carbon isotopic composition of the calcites is close to equilibrium with fluid inclusion CO₂ at metamorphic temperatures. Mass balance calculations assuming initial ¹³C values of 0 for calcite and –5 for the fluid, give integrated fluid/rock volume ratios between 0.1 and 2.0. Such high fluid/rock ratios are supported by observations on the distribution of CO₂/H₂O ratios of fluid inclusions, carbon isotopic compositions of fluid inclusion CO₂ and oxygen isotope systematics of silicates.     

Calcite dissolution kinetics in Na-Ca-Mg-Cl brines

Sedimentary basins can contain close to 20% by volume of pore fluids commonly classified as brines. These fluids can become undersaturated with respect to calcite as a result of migration, dispersive mixing, or anthropogenic injection of CO₂. This study measured calcite dissolution rates in geologically relevant Na-Ca-Mg-Cl synthetic brines (50-200 g L⁻¹ TDS). The dissolution rate dependency on brine composition, pCO₂ (0.1-1 bar), and temperature (25.0-82.5 °C) was modeled using the empirical rate equation

R=k(1-Ω)ⁿ     

Isotopic constraints on ice age fluids in active geothermal systems: Reykjanes, Iceland

The Reykjanes geothermal system is located on the landward extension of the Mid-Atlantic Ridge in southwest Iceland, and provides an on-land proxy to high-temperature hydrothermal systems of oceanic spreading centers. Previous studies of elemental composition and salinity have shown that Reykjanes geothermal fluids are likely hydrothermally modified seawater. However, δD values of these fluids are as low as −23‰, which is indicative of a meteoric water component. Here we constrain the origin of Reykjanes hydrothermal solutions by analysis of hydrogen and oxygen isotope compositions of hydrothermal epidote from geothermal drillholes at depths between 1 and 3 km. δD_EPIDOTE values from wells RN-8, -9, -10 and -17 collectively range from −60 to −78‰, and δ¹⁸O_EPIDOTE in these wells are between −3.0 and 2.3‰. The δD values of epidote generally increase along a NE trend through the geothermal field, whereas δ¹⁸O values generally decrease, suggesting a southwest to northeast migration of the geothermal upflow zone with time that is consistent with present-day temperatures and observed hydrothermal mineral zones. For comparative analysis, the meteoric-water dominated Nesjavellir and Krafla geothermal systems, which have a δD_FLUID of ∼ −79‰ and −89‰, respectively, show δD_EPIDOTE values of −115‰ and −125‰. In contrast, δD_EPIDOTE from the mixed meteoric-seawater Svartsengi geothermal system is −68‰; comparable to δD_EPIDOTE from well RN-10 at Reykjanes.Stable isotope compositions of geothermal fluids in isotopic equilibrium with the epidotes at Reykjanes are computed using published temperature dependent hydrogen and oxygen isotope fractionation curves for epidote-water, measured isotope composition of the epidotes and temperatures approximated from the boiling point curve with depth. Calculated δD and δ¹⁸O of geothermal fluids are less than 0‰, suggesting that fluids of meteoric or glacial origin are a significant component of the geothermal solutions. Additionally, δD_FLUID values in equilibrium with geothermal epidote are lower than those of modern-day fluids, whereas calculated δ¹⁸O_FLUID values are within range of the observed fluid isotope composition. We propose that modern δD_EPIDOTE and δD_FLUID values are the result of diffusional exchange between hydrous alteration minerals that precipitated from glacially-derived fluids early in the evolution of the Reykjanes system and modern seawater-derived geothermal fluids. A simplified model of isotope exchange in the Reykjanes geothermal system, in which the average starting δD_ROCK value is −125‰ and the water to rock mass ratio is 0.25, predicts a δD_FLUID composition within 1‰ of average measured values. This model resolves the discrepancy between fluid salinity and isotope composition of Reykjanes geothermal fluids, explains the observed disequilibrium between modern fluids and hydrothermal epidote, and suggests that rock-fluid interaction is the dominant control over the evolution of fluid isotope composition in the hydrothermal system.     

A synthetic fluid inclusion study of copper solubility in hydrothermal brines from 525 to 725 °C and 0.3 to 1.7 GPa

Fluid inclusions were synthesized in a piston-cylinder apparatus under mineral-buffered conditions over a range of Cl concentration (0.29 to 11.3 mol kg⁻¹), temperature (525 to 725 °C), and pressure (0.3 to 1.7 GPa). All fluids were buffered by the mineral assemblage native copper + cuprite + talc + quartz. In situ fluid composition was determined by analysing individual fluid inclusions by LA-ICPMS and independently analysing the quench solution. The solubility data provide basic information necessary to model the high temperature behaviour of Cu in magmatic-hydrothermal systems. Copper concentrations up to ∼15 wt% were measured at 630 °C and 0.34 GPa. These results give an upper limit for Cu in natural fluids and support field-based observations of similar high Cu concentrations in fluids at near-magmatic conditions. Experimental evidence indicates that Cu⁺ may form neutral chloride complexes with the general stoichiometry with n up to 4, though n ? 2 is typical for the majority of the experimental conditions. At high pressure (>∼0.5 GPa) there is evidence that hydroxide species, e.g., CuOH⁰, become increasingly important and may predominate over copper(I)-chloride complexes. The roles of fluid mixing, cooling and decompression in ore-forming environments are also discussed.     

Variations in Os isotope ratios of pyrrhotite as a result of water-rock and magma-rock interaction: Constraints from Virginia Formation-Duluth Complex contact zones

Os isotope ratios in pyrrhotite-bearing pelitic rocks of the ∼1.85 Ga Virginia Formation are variable, with perturbations linked to the emplacement of the ∼1.1 Ga Duluth Complex. Pyrrhotite in footwall rocks of the contact aureole show evidence for a mixing event at 1.1 Ga involving a low ¹⁸⁷Os/¹⁸⁸Os fluid. However, because rocks with perturbed pyrrhotite Os isotope ratios occur 1½ km or more from the Duluth Complex, the fluid is unlikely to have been of magmatic origin. Fluid inclusions in layer-parallel quartz veins provide evidence of the involvement of a boiling fluid at temperatures between ∼300 and 400 °C. Analyses of fluid inclusions via LA-ICP-MS show that the fluids contain up to 1.7 wt% Na, 1.1 wt% K, 4330 ppm Fe, 2275 ppm Zn, and 415 ppm Mg. The veins also contain pyrite or pyrrhotite, plus minor amounts of chalcopyrite, bornite, pentlandite, and sphalerite. The Re-Os isotopic ratios of pyrite from the veins indicate that they crystallized from low ¹⁸⁷Os/¹⁸⁸Os fluids (<0.2). δ¹⁸O values of vein quartz range from 7.7‰ to 9.5‰, consistent with an origin involving fluid with a relatively low δ¹⁸O value between 2‰ and 5‰. Meteoric water with such a low δ¹⁸O value could have interacted with the igneous rocks of the Complex and would have acquired Os with a low ¹⁸⁷Os/¹⁸⁸Os ratio. Strongly serpentinized olivine-rich rocks of the Complex are commonly characterized by such low δ¹⁸O values and we propose that the fluid involved in serpentinization was also responsible for the perturbation of the Os isotopic system recorded by pyrrhotite in the Virginia Formation. Two important observations are that only pyrrhotite-bearing assemblages in the contact aureole show isotopic perturbation and that intervals showing Os exchange are spatially restricted, and not uniformly distributed. Os exchange and mixing has occurred only where temperatures were sufficient to convert pyrite to pyrrhotite, and where time-integrated water-rock ratios in the aureole were high enough to provide a supply of Os.Troctolitic and gabbroic rocks of the Partridge River Intrusion, Duluth Complex, are characterized by Os isotope ratios that are indicative of variable degrees of crustal contamination (γ_Os values of ∼0-543). Xenoliths of carbonaceous and sulfidic pelitic rocks of the Virginia Formation found in the igneous rocks provide evidence that Os was released by organic matter and pyrite in the sedimentary rocks and assimilated by mantle-derived magma. However, residual pyrrhotite produced as a result of pyrite breakdown in the xenoliths is characterized by ¹⁸⁷Os/¹⁸⁸Os ratios that are much lower than anticipated and similar to those of pyrrhotite in the contact aureole. The Os exchange and addition shown by pyrrhotite in the xenoliths highlight an unusual cycle of Re-Os liberation during devolatilization, kerogen maturation, and pyrite to pyrrhotite conversion (processes that contribute to magma contamination), followed by Os uptake by pyrrhotite during back reaction involving magma and/or fluid characterized by a relatively low ¹⁸⁷Os/¹⁸⁸Os ratio. The extreme Os uptake recorded by pyrrhotite in the xenoliths, as well as the lesser degree of uptake recorded by pyrrhotite in the contact aureole, is in line with the high Os diffusivity in pyrrhotite experimentally determined by Brenan et al. (2000). Our data confirm that Os isotope ratios in pyrrhotite-bearing rocks may be readily perturbed. For this reason caution should be exercised in the interpretation of Os isotope ratios in rocks where pyrrhotite may be the primary host of Os.     

Isotopic effects in fracture-dominated reactive fluid-rock systems

A mathematical model is presented that describes the effects of pore fluid aqueous diffusion and reaction rate on the isotopic exchange between fluids and rocks in reactive geo-hydrological systems where flow is primarily through fractures. The model describes a simple system with parallel equidistant fractures, and chemical transport in the matrix slabs between fractures by aqueous diffusion through a stagnant pore fluid. The solid matrix exchanges isotopes with pore fluid by solution-precipitation at a rate characterized by a time constant, R (yr⁻¹), which is an adjustable parameter. The effects of reaction on the isotopes of a particular element in the fracture fluid are shown to depend on the ratio of the diffusive reaction length for that element (L) to the fracture spacing (b). The reaction length depends on the solid-fluid exchange rate within the matrix, the partitioning of the element between the matrix pore fluid and the matrix solid phase, the porosity and density of the matrix, and the aqueous diffusivity. For L/b < 0.3, fluid-rock isotopic exchange is effectively reduced by a factor of 2L/b relative to a standard porous flow (single porosity) model. For L/b > 1, the parallel fracture model is no different from a porous flow model. If isotopic data are available for two or more elements with different L values, it may be possible to use the model with appropriate isotopic measurements to estimate the spacing of the primary fluid-carrying fractures in natural fluid-rock systems. Examples are given using Sr and O isotopic data from mid-ocean ridge (MOR) hydrothermal vent fluids and Sr isotopes in groundwater aquifers hosted by fractured basalt. The available data for MOR systems are consistent with average fracture spacing of 1-4 m. The groundwater data suggest larger effective fracture spacing, in the range 50-500 m. In general, for fractured rock systems, the effects of fracture-matrix diffusive exchange must be considered when comparing isotopic exchange effects for different elements, as well as for estimating water age using radioactive and cosmogenic isotopes.     

The system H2O-NaCl. Part I: Correlation formulae for phase relations in temperature-pressure-composition space from 0 to 1000 °C, 0 to 5000 bar, and 0 to 1 XNaCl

Realistic simulations of fluid flow in geologic systems have severely been hampered by the lack of a consistent formulation for fluid properties for binary salt-water fluids over the temperature-pressure-composition ranges encountered in the Earth’s crust. As the first of two companion studies, a set of correlations describing the phase stability relations in the system H₂O-NaCl is developed. Pure water is described by the IAPS-84 equation of state. New correlations comprise the vapor pressure of halite and molten NaCl, the NaCl melting curve, the composition of halite-saturated liquid and vapor, the pressure of vapor + liquid + halite coexistence, the temperature-pressure and temperature-composition relations for the critical curve, and the compositions of liquid and vapor on the vapor + liquid coexistence surface. The correlations yield accurate values for temperatures from 0 to 1000 °C, pressures from 0 to 5000 bar, and compositions from 0 to 1 X_NaCl (mole fraction of NaCl). To facilitate their use in fluid flow simulations, the correlations are entirely formulated as functions of temperature, pressure and composition.     

The noble gas systematics of late-orogenic H2O-CO2 fluids, Mt Isa, Australia

The noble gases (He, Ne, Ar, Kr and Xe) are powerful geochemical tracers because they have distinctive isotopic compositions in the atmosphere, crust and mantle. This study illustrates how noble gases can be used to trace fluid origins in high-temperature metamorphic and mineralising environments; and at the same time provides new information on the composition of noble gases in deeper parts of the crust than have been sampled previously.We report data for H₂O and CO₂ fluid inclusions trapped at greenschist to amphibolite facies metamorphic conditions associated with three different styles of mineralisation and alteration in the Proterozoic Mt Isa Inlier of Australia. Sulphide fluid inclusions are dominated by crustal ⁴He. However, co-variations in fluid inclusion ²⁰Ne/²²Ne, ²¹Ne/²²Ne, ⁴⁰Ar/³⁶Ar and ¹³⁶Xe/¹³⁰Xe indicate noble gases were derived from three or more reservoirs. In most cases, the fluid inclusions elemental noble gas ratios (e.g. Ne/Xe) are close to the ranges expected in sedimentary and crystalline rocks. However, the elemental ratios have been modified in some of the samples providing evidence for independent pulses of CO_2, and interaction of CO₂ with high-salinity aqueous fluids.Compositional variation is attributed to mixing of: (i) magmatic fluids (or deeply sourced metamorphic fluids) characterised by basement-derived noble gases with ²⁰Ne/²²Ne ∼ 8.4, ²¹Ne/²²Ne ∼ 0.4, ⁴⁰Ar/³⁶Ar ∼ 40,000 and ¹³⁶Xe/¹³⁰Xe ∼ 8; (ii) basinal-metamorphic fluids with a narrow range of compositions including near-atmospheric values and (iii) noble gases derived from the meta-sedimentary host-rocks with ²⁰Ne/²²Ne ∼ 8-9.8, ²¹Ne/²²Ne < 0.1, ⁴⁰Ar/³⁶Ar < 2500 and ¹³⁶Xe/¹³⁰Xe ∼ 2.2.These data provide the strongest geochemical evidence available for the involvement of fluids from two distinct geochemical reservoirs in Mt Isa’s largest ore deposits. In addition the data show how noble gases in fluid inclusions can provide information on fluid origins, the composition of the crust’s major lithologies, fluid-rock interactions and fluid-fluid mixing or immiscibility processes.     

Bromine cycle in subduction zones through in situ Br monitoring in diamond anvil cells

The geochemical partitioning of bromine between hydrous haplogranitic melts, initially enriched with respect to Br and aqueous fluids, has been continuously monitored in situ during decompression. Experiments were carried out in diamond anvil cells from 890 °C to room temperature and from 1.7 GPa to room pressure, typically from high P, T conditions corresponding to total miscibility (presence of a supercritical fluid). Br contents were measured in aqueous fluids, hydrous melts and supercritical fluids. Partition coefficients of bromine were characterized at pressure and temperature between fluids, hydrous melts and/or glasses, as appropriate: D^Br_fluid/melt = (Br)_fluid/(Br)_melt, ranges from 2.18 to 9.2 ± 0.5 for conditions within the ranges 0.66-1.7 GPa, 590-890 °C; and D^Br_fluid/glass = (Br)_fluid/(Br)_glass ranges from 60 to 375 at room conditions. The results suggest that because high pressure melts and fluids are capable of accepting high concentrations of bromine, this element may be efficiently removed from the slab to the mantle source of arc magmas. We show that Br may be highly concentrated in subduction zone magmas and strongly enriched in subduction-related volcanic gases, because its mobility is strongly correlated with that of water during magma degassing. Furthermore, our experimental results suggest that a non negligible part of Br present in the subducted slab may remain in the down-going slab, being transported toward the transition zone. This indicates that the Br cycle in subduction zones is in fact divided in two related but independent parts: (1) a shallower one where recycled Br may leave the slab with a water and silica-bearing “fluid” leading to enriched arc magmas that return Br to the atmosphere. (2) A deeper cycle where Br may be recycled back to the mantle maybe to the transition zone, where it may be present in high pressure water-rich metasomatic fluids.