The formation of natural cryogenic brines

The source of salts in the Ca-chloridic, hypersaline brines (up to 190 g Cl L⁻¹) occurring in crystalline basement rocks in the Canadian, Fennoscandian and Bohemian Shields and their evolution have been investigated and reported. The Cl-Br-Na relationship indicates that these waters have been concentrated from seawater, by freezing during glacial times. The Na/Cl ratio (0.25 to 0.35) in the more saline fluids is compatible with cooling down to −30°C, where the most saline waters have been concentrated by a factor of 25 to 30 relative to the parent seawater.The brines formed from seawater within cryogenic troughs, along the subarctic continental margins, around ice sheets. The depressions within which the brines formed are the cryogenic analogues of the classic, evaporitic lagoon. One million years suffice to saturate with brine a 2000km-radius by 1km-depth rock volume at an H₂O removal rate of only 2.8 mm/yr. Density-induced brine migration on a continental scale takes place via fissures below the ice.Our calculations, that were performed on a hypothetical ice sheet with dimensions compatible with the Laurentide ice sheet, demonstrate that during 1m.y., a 60m thick cryogenic sediment section could have formed. However, the precipitated minerals (mirabilite and hydrohalite) are repeatedly dispersed by the advance and retreat of the ice sheet, dissolved by melt water-seawater mixtures, and eroded during postglacial uplift, leaving almost no trace in the geological record.The cryogenic brines formed intermittently during and between glacial periods. The repeating advance and retreat of the ice sheets exerted a major control on the direction and intensity of brine flow. The cryogenic concentration of seawater and the migration of brine towards the center of the glaciostatic depression occurred mainly during the build up of the ice sheet, while reversal of the water flow from the center of the cryogenic basin outwards happened upon deglaciation. The flow of the waters in the subsurface was, inevitably, accompanied by significant dilution with melt water from the ice sheets.Using a “granitic” U concentration of 4 ppm and a (Ca-Mg mass balance based) rock/water ratio anywhere between 3.4 and 6.8 kg L⁻¹, a few hundred thousand years of brine-rock interaction are sufficient for the growth of ¹²⁹I in the most saline Canadian Shield brine to its present concentration (3.4×10⁸ atoms ¹²⁹I L⁻¹). Hence, both the formation of the saline fluids and their emplacement in their present sites occurred most likely within the Pleistocene.The young age calculated for cryogenic brines in crystalline shields and the dynamic water flow therein should raise concern about the planning and construction of high-grade nuclear waste repositories in such rocks, which are already under way.     

The major-ion composition of Permian seawater

The major-ion (Mg²⁺, Ca²⁺, Na⁺, K⁺, SO₄²⁻, and Cl⁻) composition of Permian seawater was determined from chemical analyses of fluid inclusions in marine halites. New data from the Upper Permian San Andres Formation of Texas (274-272 Ma) and Salado Formation of New Mexico (251 Ma), analyzed by the environmental scanning electron microscopy (ESEM) X-ray energy-dispersive spectrometry (EDS) method, along with published chemical compositions of fluid inclusions in Permian marine halites from North America (two formations of different ages) and the Central and Eastern European basins (eight formations of four different ages) show that Permian seawater shares chemical characteristics with modern seawater, including SO₄²⁻ > Ca²⁺ at the point of gypsum precipitation, evolution into Mg²⁺-Na⁺-K⁺-SO₄²⁻-Cl⁻ brines, and Mg²⁺/K⁺ ratios ∼5. Permian seawater, however, is slightly depleted in SO₄²⁻ and enriched in Ca²⁺, although modeling results do not rule out Ca²⁺ concentrations close to those in present-day seawater. Na⁺ and Mg²⁺ in Permian seawater are close to (slightly below) their concentrations in modern seawater. Permian and modern seawater are both classified as aragonite seas, with Mg²⁺/Ca²⁺ ratios >2, conditions favorable for precipitation of aragonite and magnesian calcite as ooids and cements.The chemistry of Permian seawater was modeled using the chemical composition of brine inclusions for three periods: Lower Permian Asselian-Sakmarian (296-283 Ma), Lower Permian Artinskian-Kungurian (283-274 Ma), and Upper Permian Tatarian (258-251 Ma). Parallel changes in the chemistry of brine inclusions from equivalent age evaporites in North America, Central Europe, and Eastern Europe show that seawater underwent secular variations in chemistry over the 50 million years of the Permian. Modeled SO₄²⁻ concentrations are 20 mmol per kg H₂O (mmolal) and 19 mmolal in the Asselian-Sakmarian and Artinskian-Kungurian, with higher concentrations in the Upper Permian Tatarian (23 mmolal). Modeled Ca²⁺ is at or above its concentration in modern seawater throughout the Permian. Mg²⁺ is close to (slightly below) its concentration in modern seawater (55 mmolal) in the Asselian-Sakmarian (52 mmolal), and Tatarian (52 mmolal), but slightly higher than modern seawater in the Artinskian-Kungurian (60 mmolal). Mg²⁺/Ca²⁺ ratios are 3.5 (total range = 2.7 to 5.5) in the Lower Permian and rose slightly to 3.7 (total range = 3.1 to 5.8) in the Upper Permian, primarily due to decreases in Ca²⁺. These results are consistent with models that predict oscillations in the major-ion composition of Phanerozoic seawater on the basis of changes in the midocean ridge/river water flux ratio driven by changes in the rate of midocean ridge crust production.The Permian was characterized by low sea levels, icehouse conditions, and southern hemisphere glaciation. Such conditions, analogous to the present ice age, and the similarities between Permian seawater and modern seawater, all suggest that general Phanerozoic supercycles, driven by mantle convection and global volcanicity, also control the major-ion chemistry of seawater.     

Chemistry of fluids from a natural analogue for a geological CO2 storage site (Montmiral, France): Lessons for CO2–water–rock interaction assessment and monitoring

Chemical and isotope studies of natural CO₂ accumulations aid in assessing the chemical effects of CO₂ on rock and thus provide a potential for understanding the long-term geochemical processes involved in CO₂ geological storage. Several natural CO₂ accumulations were discovered during gas and oil exploration in France’s carbogaseous peri-Alpine province (south-eastern France) in the 1960s. One of these, the Montmiral accumulation at a depth of more than 2400 m, is currently being exploited. The chemical composition of the water collected at the wellhead has changed in time and the final salinity exceeds 75 g/L. These changes in time can be explained by assuming that the fraction of the reservoir brine in the recovered brine–CO₂–H₂O mixture varies, resulting in variable proportions of H₂O and brine in the sampled water. The proportions can be estimated in selected samples due to the availability of gas and water flowrate data. These data enabled the reconstruction of the chemical and isotope composition of the brine. The proportions of H₂O and brine can also be estimated from isotope (δ²H, δ¹⁸O) composition of collected water and δ¹⁸O of the sulfates or CO₂. The reconstituted brine has a salinity of more than 85 g/L and, according to its Br⁻ content and isotope (δ²H, δ¹⁸O, δ³⁴S) composition, originates from an evaporated Triassic seawater that underwent dilution by meteoric water. The reconstitution of the brine’s chemical composition enabled an evaluation of the CO₂–water–rock interactions based on: (1) mineral saturation indices; and (2) comparison with initial evaporated Triassic seawater. Dissolution of K- and SO₄-containing minerals such as K-feldspar and anhydrite, and precipitation of Ca and Mg containing minerals that are able to trap CO₂ (carbonates) are highlighted. The changes in concentration of these elements in the brine, which are attributed to CO₂ interactions, illustrate the relevance of monitoring the water quality at future industrial CO₂ storage sites.     

History of carbonate ion concentration over the last 100 million years

Instead of having been more or less constant, as once assumed, it is now apparent that the major ion chemistry of the oceans has varied substantially over time. For instance, independent lines of evidence suggest that calcium concentration ([Ca²⁺]) has approximately halved and magnesium concentration ([Mg²⁺]) approximately doubled over the last 100 million years. On the other hand, the calcite compensation depth, and hence the CaCO₃ saturation, has varied little over the last 100 My as documented in deep sea sediments. We combine these pieces of evidence to develop a proxy for seawater carbonate ion concentration ([CO₃²⁻]) over this period of time. From the calcite saturation state (which is proportional to the product of [Ca²⁺] times [CO₃²⁻], but also affected by [Mg²⁺]), we can calculate seawater [CO₃²⁻]. Our results show that [CO₃²⁻] has nearly quadrupled since the Cretaceous. Furthermore, by combining our [CO₃²⁻] proxy with other carbonate system proxies, we provide calculations of the entire seawater carbonate system and atmospheric CO₂. Based on this, reconstructed atmospheric CO₂ is relatively low in the Miocene but high in the Eocene. Finally, we make a strong case that seawater pH has increased over the last 100 My.     

Diel Aquatic CO<Subscript>2</Subscript> System Dynamics of a Bermudian Mangrove Environment

Mangrove ecosystems play an important, but understudied, role in the cycling of carbon in tropical and subtropical coastal ocean environments. In the present study, we examined the diel dynamics of seawater carbon dioxide (CO₂) and dissolved oxygen (DO) for a mangrove-dominated marine ecosystem (Mangrove Bay) and an adjacent intracoastal waterway (Ferry Reach) on the island of Bermuda. Spatial and temporal trends in seawater carbonate chemistry and associated variables were assessed from direct measurements of dissolved inorganic carbon, total alkalinity, dissolved oxygen (DO), temperature, and salinity. Diel pCO₂ variability was interpolated across hourly wind speed measurements to determine variability in daily CO₂ fluxes for the month of October 2007 in Bermuda. From these observations, we estimated rates of net sea to air CO₂ exchange for these two coastal ecosystems at 59.8 ± 17.3 in Mangrove Bay and 5.5 ± 1.3 mmol m⁻² d⁻¹ in Ferry Reach. These results highlight the potential for large differences in carbonate system functioning and sea-air CO₂ flux in adjacent coastal environments. In addition, observation of large diel variability in CO₂ system parameters (e.g., mean pCO₂: 390–2,841 μatm; mean pH_T: 8.05–7.34) underscores the need for careful consideration of diel cycles in long-term sampling regimes and flux estimates.     

Dissolved inorganic carbon (DIC) and hydrologic mixing in a subtropical riverine estuary,southwest Florida,USA

Physical and chemical parameters were measured in a subtropical estuary with a blind river source in southwest Florida, United States, to assess seasonal discharge of overland flow and groundwater in hydrologic mixing. Water temperature, pH, salinity, alkalinity, dissolved inorganic carbon (DIC), δ¹⁸O, and δ¹³C_DIC varied significantly due to seasonal rainfall and climate. Axial distribution of the physical and chemical parameters constrained by tidal conditions during sampling showed that river water at low tide was a mixture of freshwater from overland flow and saline ground-water in the wet season and mostly saline groundwater in the dry season. Relationships between salinity and temperature, δ¹⁸O, and DIC for both the dry and wet seasons showed that DIC was most sensitive to seawater mixing in the estuary as DIC changed in concentration between values measured in river water at the tidal front to the most seaward station. A salinity-δ¹³C_DIC model was able to describe seawater mixing in the estuary for the wet season but not for the dry season because river water salinity was higher than that of seawater and the salinity gradient between seawater and river water was small. A DIC-δ¹³C_DIC mixing model was able to describe mixing of carbon from sheet flow and river water at low tide, and river water and seawater at high tide for both wet and dry seasons. The DIC-δ¹³C_DIC model was able to predict the seawater end member DIC for the wet season. The model was not able to predict the seawater end member DIC for the dry season data due to secondary physical and biogeochemical processes that altered estuarine DIC prior to mixing with seawater. The results of this study suggest that DIC and δ¹³C_DIC can provide additional insights into mixing of river water and seawater in estuaries during periods where small salinity gradients between river water and seawater and higher river water salinities preclude the use of salinity-carbon models.     

Experimental determination of the stability and stoichiometry of sulphide complexes of silver(I) in hydrothermal solutions to 400°C

The solubility of silver sulphide (acanthite/argentite) has been measured in aqueous sulphide solutions between 25 and 400°C at saturated water vapour pressure and 500 bar to determine the stability and stoichiometry of sulphide complexes of silver(I) in hydrothermal solutions. The experiments were carried out in a flow-through autoclave, connected to a high-performance liquid chromatographic pump, titanium sampling loop, and a back-pressure regulator on line. Samples for silver determination were collected via the titanium sampling loop at experimental temperatures and pressures. The solubilities, measured as total dissolved silver, were in the range 1.0 × 10⁻⁷ to 1.30 × 10⁻⁴ mol kg⁻¹ (0.01 to 14.0 ppm), in solutions of total reduced sulphur between 0.007 and 0.176 mol kg⁻¹ and pH_T,p of 3.7 to 12.7. A nonlinear least squares treatment of the data demonstrates that the solubility of silver sulphide in aqueous sulphide solutions of acidic to alkaline pH is accurately described by the reactions0.5Ag₂S(s) + 0.5H₂S(aq) = AgHS(aq) K_s,1110.5Ag₂S(s) + 0.5H₂S(aq) + HS⁻ = Ag(HS)2− K_s,122Ag₂S(s) + 2HS⁻ = Ag₂S(HS)22− K_s,232where AgHS(aq) is the dominant species in acidic solutions, Ag(HS)2− under neutral pH conditions and Ag₂S(HS)22− in alkaline solutions. With increasing temperature the stability field of Ag(HS)2− increases and shifts to more alkaline pH in accordance with the change in the first ionisation constant of H₂S(aq). Consequently, Ag₂S(HS)22− is not an important species above 200°C. The solubility constant for the first reaction is independent of temperature to 300°C, with values in the range logK_s,111 = −5.79 (±0.07) to −5.59 (±0.09), and decreases to −5.92 (±0.16) at 400°C. The solubility constant for the second reaction increases almost linearly with inverse temperature from logK_s,122 = −3.97 (±0.04) at 25°C to −1.89 (±0.03) at 400°C. The solubility constant for the third reaction increases with temperature from logK_s,232 = −4.78 (±0.04) at 25°C to −4.57 (±0.18) at 200°C. All solubility constants were found to be independent of pressure within experimental uncertainties. The interaction between Ag⁺ and HS⁻ at 25°C and 1 bar to form AgHS(aq) has appreciable covalent character, as reflected in the exothermic enthalpy and small entropy of formation. With increasing temperature, the stepwise formation reactions become progressively more endothermic and are accompanied by large positive entropies, indicating greater electrostatic interaction. The aqueous speciation of silver is very sensitive to fluid composition and temperature. Below 100°C silver(I) sulphide complexes predominate in reduced sulphide solutions, whereas Ag⁺ and AgClOH⁻ are the dominant species in oxidised waters. In high-temperature hydrothermal solutions of seawater salinity, chloride complexes of silver(I) are most important, whereas in dilute hydrothermal fluids of meteoric origin typically found in active geothermal systems, sulphide complexes predominate. Adiabatic boiling of dilute and saline geothermal waters leads to precipitation of silver sulphide and removal of silver from solution. Conductive cooling has insignificant effects on silver mobility in dilute fluids, whereas it leads to quantitative loss of silver for geothermal fluids of seawater salinity.     

Remnants of premetamorphic fluid and oxygen isotopic signatures in eclogites and garnet clinopyroxenite from the Dabie-Sulu terranes, eastern China

A combined oxygen‐isotope and fluid‐inclusion study has been carried out on high‐ and ultrahigh‐pressure metamorphic (HP/UHPM) eclogites and garnet clinopyroxenite from the Dabie‐Sulu terranes in eastern China. Coesite‐bearing eclogites/garnet clinopyroxenite and quartz eclogites have a wide range in whole‐rock δ¹⁸O_VSMOW, from 0 to 11‰. The high‐T oxygen‐isotope fractionations preserved between quartz and garnet preclude significant retrograde isotope exchange during exhumation, and the wide range in whole‐rock oxygen‐isotope composition is thought to be a presubduction signature of the precursors. Aqueous fluids with variable salinities and gas species (N₂‐, CO₂‐, or CH₄‐rich), are trapped as primary inclusions in garnet, omphacite and epidote, and in quartz blebs enclosed within eclogitic minerals. In high‐δ¹⁸O HP/UHPM rocks from Hujialin and Shima, high‐salinity brine and/or N₂ inclusions occur in garnet porphyroblasts, which also contain inclusions of coesite, Cl‐rich blue amphibole and dolomite. In contrast, in low‐δ¹⁸O eclogites from Qinglongshan and Huangzhen, the Cl concentrations in amphibole are very low, < 0.2 wt.%, and low‐salinity aqueous inclusions occur in quartz inclusions in epidote porphyroblasts and in epidote cores. These low‐salinity fluid inclusions are believed to be remnants of meteoric water, although the fluid composition was modified during pre‐ and syn‐peak HP/UHPM. Eclogites at Houshuichegou and Hetang contain CH₄‐rich fluid inclusions, coexisting with high‐salinity brine inclusions. Methane was probably formed under the influence of CO₂‐rich aqueous fluids during serpentinisation of mantle‐derived peridotites prior to or during plate subduction. Remnants of premetamorphic low‐ to high‐salinity aqueous fluid with minor N₂ and/or other gas species preserved in the Dabie‐Sulu HP/UHPM eclogites and garnet clinopyroxenite indicate a great diversity of initial fluid composition in the precursors, implying very limited fluid–rock interaction during syn‐ and post‐peak HP/UHPM.     

Mercury in the Southern Ocean

We present here the first mercury speciation study in the water column of the Southern Ocean, using a high-resolution south-to-north section (27 stations from 65.50°S to 44.00°S) with up to 15 depths (0-4440 m) between Antarctica and Tasmania (Australia) along the 140°E meridian. In addition, in order to explore the role of sea ice in Hg cycling, a study of mercury speciation in the “snow-sea ice-seawater” continuum was conducted at a coastal site, near the Australian Casey station (66.40°S; 101.14°E). In the open ocean waters, total Hg (Hg_T) concentrations varied from 0.63 to 2.76 pmol L⁻¹ with “transient-type” vertical profiles and a latitudinal distribution suggesting an atmospheric mercury source south of the Southern Polar Front (SPF) and a surface removal north of the Subantartic Front (SAF). Slightly higher mean Hg_T concentrations (1.35 ± 0.39 pmol L⁻¹) were measured in Antarctic Bottom Water (AABW) compared to Antarctic Intermediate water (AAIW) (1.15 ± 0.22 pmol L⁻¹). Labile Hg (Hg_R) concentrations varied from 0.01 to 2.28 pmol L⁻¹, with a distribution showing that the Hg_T enrichment south of the SPF consisted mainly of Hg_R (67 ± 23%), whereas, in contrast, the percentage was half that in surface waters north of PFZ (33 ± 23%). Methylated mercury species (MeHg_T) concentrations ranged from 0.02 to 0.86 pmol L⁻¹. All vertical MeHg_T profiles exhibited roughly the same pattern, with low concentrations observed in the surface layer and increasing concentrations with depth up to an intermediate depth maximum. As for Hg_T, low mean MeHg_T concentrations were associated with AAIW, and higher ones with AABW. The maximum of MeHg_T concentration at each station was systematically observed within the oxygen minimum zone, with a statistically significant MeHg_Tvs Apparent Oxygen Utilization (AOU) relationship (p < 0.001). The proportion of Hg_T as methylated species was lower than 5% in the surface waters, around 50% in deep waters below 1000 m, reaching a maximum of 78% south of the SPF. At Casey coastal station Hg_T and Hg_R concentrations found in the “snow-sea ice-seawater” continuum were one order of magnitude higher than those measured in open ocean waters. The distribution of Hg_T there suggests an atmospheric Hg deposition with snow and a fractionation process during sea ice formation, which excludes Hg from the ice with a parallel Hg enrichment of brine, probably concurring with the Hg enrichment of AABW observed in the open ocean waters. Contrastingly, MeHg_T concentrations in the sea ice environment were in the same range as in the open ocean waters, remaining below 0.45 pmol L⁻¹. The MeHg_T vertical profile through the continuum suggests different sources, including atmosphere, seawater and methylation in basal ice. Whereas Hg_T concentrations in the water samples collected between the Antarctic continent and Tasmania are comparable to recent measurements made in the other parts of the World Ocean (e.g., Soerensen et al., 2010), the Hg species distribution suggests distinct features in the Southern Ocean Hg cycle: (i) a net atmospheric Hg deposition on surface water near the ice edge, (ii) the Hg enrichment in brine during sea ice formation, and (iii) a net methylation of Hg south of the SPF.     

An evaporated seawater origin for the ore-forming brines in unconformity-related uranium deposits (Athabasca Basin, Canada): Cl/Br and δCl analysis of fluid inclusions

Analyses of halogen concentration and stable chlorine isotope composition of fluid inclusions from hydrothermal quartz and carbonate veins spatially and temporally associated with giant unconformity-related uranium deposits from the Paleoproterozoic Athabasca Basin (Canada) were performed in order to determine the origin of chloride in the ore-forming brines. Microthermometric analyses show that samples contain variable amounts of a NaCl-rich brine (Cl concentration between 120,000 and 180,000 ppm) and a CaCl₂-rich brine (Cl concentration between 160,000 and 220,000 ppm). Molar Cl/Br ratios of fluid inclusion leachates range from ∼100 to ∼900, with most values between 150 and 350. Cl/Br ratios below 650 (seawater value) indicate that the high salinities were acquired by evaporation of seawater. Most δ³⁷Cl values are between −0.6‰ and 0‰ (seawater value) which is also compatible with a common evaporated seawater origin for both NaCl- and CaCl₂-rich brines.Slight discrepancies between the Cl concentration, Cl/Br, δ³⁷Cl data and seawater evaporation trends, indicate that the evaporated seawater underwent secondary minor modification of its composition by: (i) mixing with a minor amount of halite-dissolution brine or re-equilibration with halite during burial; (ii) dilution in a maximum of 30% of connate and/or formation waters during its migration towards the base of the Athabasca sandstones; (iii) leaching of chloride from biotites within basement rocks and (iv) water loss by hydration reactions in alteration haloes linked to uranium deposition.The chloride in uranium ore-forming brines of the Athabasca Basin has an unambiguous dominantly marine origin and has required large-scale seawater evaporation and evaporite deposition. Although the direct evidence for evaporative environments in the Athabasca Basin are lacking due to the erosion of ∼80% of the sedimentary pile, Cl/Br ratios and δ³⁷Cl values of brines have behaved conservatively at the basin scale and throughout basin history.     

Experimental study of aluminum speciation in fluoride-rich supercritical fluids

The solubility of the albite-paragonite-quartz mineral assemblage was measured as a function of NaCl and fluorine concentration at 400°C, 500 bars and at 450°C, 500 and 1000 bars. Decreasing Al concentrations with increasing NaCl molality in F-free fluids of low salinity (mNaCl < 0.01) demonstrates that Al(OH)₄⁻ dominates Al speciation and is formed according to the reaction 0.5 NaAl₃Si₃O₁₂H_2(cr)+2 H₂O = 0.5 NaAlSi₃O_8(cr)+Al(OH)₄⁻+H⁺. Log K results for this reaction are −11.28 ± 0.10 and −10.59 ± 0.10 at 400°C, 500 bars and 450°C, 1000 bars, respectively. Upon further salinity increase, Al concentration becomes constant (at 400°C, 500 bars) or even rises (at 450°C, 1000 bars). The observed Al behavior can be explained by the formation of NaAl(OH)₄⁰_(aq) or NaAl(OH)₃Cl_(aq)⁰. The calculated constant for the reaction Al(OH)₄⁻+Na⁺=NaAl(OH)₄⁰_(aq) expressed in log units is equal to 2.46 and 2.04 at 400°C, 500 bars and 450°C, 1000 bars, respectively. These values are in good agreement with the predictions given in Diakonov et al. (1996). Addition of fluoride at m(NaCl) = const = 0.5 caused a sharp increase in Al concentration in equilibrium with the albite-paragonite-quartz mineral assemblage. As fluid pH was also constant, this solubility increase indicates strong aluminum-fluoride complexation with the formation of NaAl(OH)₃F_(aq)⁰ and NaAl(OH)₂F₂⁰_(aq), according to 0.5 NaAl₃Si₃O₁₂H_2(cr)+Na⁺+HF_(aq)⁰+H₂O = 0.5 NaAlSi₃O_8(cr)+ NaAl(OH)₃F_(aq)⁰+H⁺, log K = −5.17 and −5.23 at 400°C and 450°C, 500 bars, respectively, and 0.5 NaAl₃Si₃O₁₂H_2(cr)+Na⁺+2 HF_(aq)⁰ = 0.5 NaAlSi₃O_8(cr)+NaAl(OH)₂F₂⁰_(aq)+H⁺, log K = −2.19 and −1.64 at the same P-T conditions. It was found that temperature increase and pressure decrease promote the formation of Na-Al-OH-F species. Stability of NaAl(OH)₂F₂⁰_(aq) in low-density fluids also increases relative to NaAl(OH)₃F_(aq)⁰. These complexes, together with Al(OH)₂F_(aq)⁰ and AlOHF₂⁰_(aq), whose stability constants were calculated from the corundum solubility measured by Soboleva and Zaraisky (1990) and Zaraisky (1994), are likely to dominate Al speciation in metamorphic fluids containing several ppm of fluorine.     

Study of shallow groundwater quality evolution under saline intrusion with environmental isotopes and geochemistry

Evolution of the shallow groundwater quality under saline intrusion in porous aquifer system has been studied with environmental isotopes and geochemistry in the Laizhou Bay area, China. Two campaigns of water sampling from various sources were carried out in spring and winter for environmental isotopic and chemical analyses. The origin of groundwater salinity from intrusion of both modern seawater and deep brine water was identified by analysing the correlations between ¹⁸O, D, T, Cl⁻, SO₄²⁻ and electrical conductivity. The results indicate that the brine is originated from evaporating and concentrating of intruded seawater and its δD and δ¹⁸O are different from modern seawater but similar to those of mixture of seawater with fresh groundwater. It is hard to distinguish the salinity origin in this area by the δD–δ¹⁸O relationship alone. The relations between δ¹⁸O and conductivity, Cl⁻ and SO₄²⁻ have been used to identify the salinity origin due to the distinct difference in salinity between the brine and seawater, conjunctively with use of T. A threshold of T = 12 TU was adopted to identify the origin of saline groundwater.     

Formation of tonalite in island arcs by seawater-induced anatexis of mafic rocks; evidence from the Lyngen Magmatic Complex, North Norwegian Caledonides

The Lyngen Magmatic Complex (LMC) of North Norway, consists of a western suite of layered gabbros of normal-mid oceanic ridge basalt (N-MORB) affinity and an eastern suite of layered gabbronorites, quartz-bearing gabbros and diorites/quartz-diorites of IAT (island-arc tholeiitte) to boninitic affinity. The boundary between the suites is defined by a large-scale ductile shear zone, the Rypdalen shear zone (RSZ). In this shear zone anatectic tonalites were generated by partial melting of the gabbro in the presence of an H₂O bearing fluid phase.Quartz from the tonalites contains early secondary and secondary liquid-dominated inclusions (88-99 wt.% H₂O), with an average salinity of 18 wt.% (calculated as NaCl_eq). Combined gas and ion chromatography shows that the major ions in the fluid are Cl⁻, Ca²⁺, Na⁺ with smaller amounts of K⁺, Mg²⁺, Sr²⁺, Br⁻ and NO₃⁻. The dominant non-H₂O volatile species is N₂ (0.5-10%), and small amounts of CO₂, CH₄ and other hydrocarbons are also present.The cation concentrations in the fluid are variable, due to element exchange during interaction of the fluids with the tonalites, amphibolites and metagabbros of the RSZ. The fluid contributed Na⁺ and K⁺ to the melt and gained Ca²⁺ in exchange, explaining the variable Na⁺/Ca²⁺ ratio of the fluid. The Br⁻ and Cl⁻ contents of the fluid inclusions plot on the same line as evaporating sea water, which strongly suggests a seawater origin for the fluid phase, and a seawater source fits well with other geochemical signatures and the tectonic setting of the LMC.It is suggested that seawater escaped from a subducting slab and was channelled along the Rypdalen shear zone. This caused anatexis of the gabbro, generating tonalitic melts at 0.5-0.9 GPa and 680-800°C.     

Can pelagic net heterotrophy account for carbon fluxes from eastern Canadian lakes?

Lakes worldwide are commonly oversaturated with CO₂, however the source of this CO₂ oversaturation is not well understood. To examine the magnitude of the C flux to the atmosphere and determine if an excess of respiration (R) over gross primary production (GPP) is sufficient to account for this C flux, metabolic parameters and stable isotopes of dissolved O₂ and C were measured in 23 Québec lakes. All of the lakes sampled were oversaturated with CO₂ over the sampling period, on average 221 ± 25%. However, little evidence was found to conclude that this CO₂ oversaturation was the result of an excess of pelagic R over GPP. In lakes Croche and à l’Ours, where CO₂ flux, R and GPP were measured weekly, the annual difference between pelagic GPP and R, or net primary production (NPP), was not sufficient to account for the size of the CO₂ flux to the atmosphere. In Lac Croche average annual NPP was 14.4 mg C m⁻² d⁻¹ while the average annual flux of CO₂ to the atmosphere was 34 mg C m⁻² d⁻¹. In Lac à l’Ours average annual NPP was −9.1 mg C m⁻² d⁻¹ while the average annual flux of CO₂ to the atmosphere was 55 mg C m⁻² d⁻¹. In all of the lakes sampled, O₂ saturation averaged 104.0 ± 1.7% during the ice-free season and the isotopic composition of dissolved O₂ (δ¹⁸O_DO) was 22.9 ± 0.3‰, lower than atmospheric values and indicative of net autotrophy. Carbon evasion was not a function of R, nor did the isotopic signature of dissolved CO₂ in the lakes present evidence of excess R over GPP. External inputs of C must therefore subsidize the lake to explain the continued CO₂ oversaturation. The isotopic composition of dissolved inorganic C (δ¹³C_DIC) indicates that the CO₂ oversaturation cannot be attributed to in situ aerobic respiration. δ¹³C_DIC reveals a source of excess C enriched in ¹³C, which may be accounted for by anaerobic sediment respiration or groundwater inputs followed by kinetic isotope fractionation during degassing under open system conditions.     

The major-ion composition of silurian seawater

One-hundred fluid inclusions in Silurian marine halite were analyzed in order to determine the major-ion composition of Silurian seawater. The samples analyzed were from three formations in the Late Silurian Michigan Basin, the A-1, A-2, and B Evaporites of the Salina Group, and one formation in the Early Silurian Canning Basin (Australia), the Mallowa Salt of the Carribuddy Group. The results indicate that the major-ion composition of Silurian seawater was not the same as present-day seawater. The Silurian ocean had lower concentrations of Mg²⁺, Na⁺, and SO₄²⁻, and much higher concentrations of Ca²⁺ relative to the ocean’s present-day composition. Furthermore, Silurian seawater had Ca²⁺ in excess of SO₄²⁻. Evaporation of Silurian seawater of the composition determined in this study produces KCl-type potash minerals that lack the MgSO₄-type late stage salts formed during the evaporation of present-day seawater. The relatively low Na⁺ concentrations in Silurian seawater support the hypothesis that oscillations in the major-ion composition of the oceans are primarily controlled by changes in the flux of mid-ocean ridge brine and riverine inputs and not global or basin-scale, seawater-driven dolomitization. The Mg²⁺/Ca²⁺ ratio of Silurian seawater was ∼1.4, and the K⁺/Ca²⁺ ratio was ∼0.3, both of which differ from the present-day counterparts of 5 and 1, respectively. Seawaters with Mg²⁺/Ca²⁺ <2 facilitate the precipitation of low-magnesian calcite (mol % Mg < 4) marine ooids and submarine carbonate cements whereas seawaters with Mg²⁺/Ca²⁺ >2 (e.g., modern seawater) facilitate the precipitation of aragonite and high-magnesian calcite. Therefore, the early Paleozoic calcite seas were likely due to the low Mg²⁺/Ca²⁺ ratio of seawater, not the pCO₂ of the Silurian atmosphere.     

Genesis of fumarolic emissions as inferred by isotope mass balances: CO2 and water at Vulcano Island, Italy

We have developed a quantitative model of CO₂ and H₂O isotopic mixing between magmatic and hydrothermal gases for the fumarolic emissions of the La Fossa crater (Vulcano Island, Italy). On the basis of isotope balance equations, the model takes into account the isotope equilibrium between H₂O and CO₂ and extends the recent model of chemical and energy two-end-member mixing by Nuccio et al. (1999). As a result, the H₂O and CO₂ content and the δD, δ¹⁸O, and δ¹³C isotope compositions for both magmatic and hydrothermal end-members have been assessed. Low contributions of meteoric steam, added at a shallow depth, have been also recognized and quantified in the fumaroles throughout the period from 1988 to 1998. Nonequilibrium oxygen isotope exchange also seems to be occurring between ascending gases and wall rocks along some fumarolic conduits.The δ¹³C_CO2 of the magmatic gases varies around −3 to 1‰ vs. Peedee belemnite (PDB), following a perfect synchronism with the variations of the CO₂ concentration in the magmatic gases. This suggests a process of isotope fractionation because of vapor exsolution caused by magma depressurization. The hydrogen isotopes in the magmatic gases (−1 to −‰ vs. standard mean ocean water [SMOW]), as well as the above δ¹³C_CO2 value, are coherent with a convergent tectonic setting of magma generation, where the local mantle is widely contaminated by fluids released from the subducted slab. Magma contamination in the crust probably amplifies this effect.The computed isotope composition of carbon and hydrogen in the hydrothermal vapors has been used to calculate the δD and δ¹³C of the entire hydrothermal system, including mixed H₂O-CO₂ vapor, liquid water, and dissolved carbon. We have computed values of about 10‰ vs. SMOW for water and −2 to −6.5‰ vs. PDB for CO₂. On these grounds, we think that Mediterranean marine water (δD_H2O ≈ 10‰) feeds the hydrothermal system. It infiltrates at depth throughout the local rocks, reaching oxygen isotope equilibrium at high temperatures. Interaction processes between magmatic gases and the evolving seawater also seem to occur, causing the dissolution of isotopically fractionated aqueous CO₂ and providing the source for hydrothermal carbon. These results have important implications concerning fluid circulation beneath Vulcano and address the more convenient routine of geochemical surveillance.     

Multi-proxy approach (H/H, O/O, C/C and Sr/Sr) for the evolution of carbonate-rich groundwater in basalt dominated aquifer of Axum area, northern Ethiopia

Isotopic and chemical composition of groundwater from wells and springs, and surface water from the basalt-dominated Axum area (northern Ethiopia) provides evidence for the origin of water and dissolved species. Shallow (depth < 40 m) and deep groundwater are distinguished by both chemical and isotopic composition. Deep groundwater is significantly enriched in dissolved inorganic carbon up to 40 mmol l⁻¹ and in concentrations of Ca²⁺, Mg²⁺, Na⁺ and Si(OH)₄ compared to the shallow type.The δ²H and δ¹⁸O values of all solutions clearly indicate meteoric origin. Shifts from the local meteoric water line are attributed to evaporation of surface and spring water, and to strong water–rock interaction. The δ¹³C_DIC values of shallow groundwater between −12 and −7‰ (VPDB) display the uptake of CO₂ from local soil horizons, whereas δ¹³C_DIC of deep groundwater ranges from −5 to +1‰. Considering open system conditions with respect to gaseous CO₂, δ¹³C_DIC = +1‰ of the deep groundwater with highest PCO₂ = 10^−0.9 atm yields δ¹³C_CO2(gas) ≈ −5‰, which is close to the stable carbon isotopic composition of magmatic CO₂. Accordingly, stable carbon isotope ratios within the above range are referred to individual proportions of CO₂ from soil and magmatic origin. The uptake of magmatic CO₂ results in elevated cations and Si(OH)₄ concentrations. Weathering of local basalts is documented by ⁸⁷Sr/⁸⁶Sr ratios of the groundwater from 0.7038 to 0.7059. Highest values indicate Sr release from the basement rocks. Besides weathering of silicates, neoformation of solids has to be considered, which results in the formation of, e.g., kaolinite and montmorillonite. In several solutions supersaturation with respect to calcite is reached by outgassing of CO₂ from the solution leading to secondary calcite formation.     

四川盆地三叠纪古盐湖已达钾石盐析出阶段——来自石盐流体包裹体化学组成的约束

石盐的流体包裹体成分可提供古流体组成的物理化学信息,用以探查卤水组成变化及环境演化规律等。四川盆地位于上扬子地台,其中的早-中三叠纪沉积建造是中国海相找钾的有利层位之一。获取石盐沉积时期的卤水成分信息,是深刻认识四川盆地古海水蒸发浓缩程度的重要途径。文章利用激光剥蚀电感耦合等离子体质谱法,对采自川东地区长平3井嘉陵江组的石盐流体包裹体开展了化学组成分析,结果显示古卤水化学类型为Mg_SO4型;流体包裹体中的ρ(K~+)与现代海水浓缩到钾石盐析出阶段的ρ(K~+)基本一致,可能揭示了盆地三叠纪时期古卤水已达到钾石盐析出阶段,对四川盆地沉积环境演化及钾盐成矿规律研究等具有重要的理论意义。     

Ore metal redistribution by hydrocarbon–brine and hydrocarbon–halide melt phases, North Range footwall of the Sudbury Igneous Complex, Ontario, Canada

We report methane-dominant hydrocarbon (fluid) inclusions (CH₄±C₂H₆–C₂H₂, C₃H₈) coexisting with primary brine inclusions and secondary halide melt (solid NaCl) inclusions in Au–Pt-rich quartz-sulfide-epidote alteration veins associated with the footwall-style Cu–PGE (platinum-group element)–Au deposits at the Fraser Mine (North Range of the Sudbury Igneous Complex). Evidence for coentrapment of immiscible hydrocarbon–brine, and hydrocarbon–halide melt mixtures is demonstrated. A primary CH₄–brine assemblage was trapped during quartz growth at relatively low T (min. T _trapping∼145–315°C) and P (max. P _trapping∼500 bar), prior to the crystallization of sulfide minerals in the veins. Secondary inclusions contain solid halite and a mixture of CH₄, C₂H₆–C₂H₂ and C₃H₈ and were trapped at a minimum T of ∼710°C. The halite inclusions may represent halide melt that exsolved from crystallizing sulfide ores that texturally postdate (by replacement) early alteration quartz hosting the primary, lower T brine–CH₄ assemblage. Laser ablation ICP-MS analyses show that the brine, hydrocarbon and halide melt inclusions contain significant concentrations of Cu (0.1–1 wt% range), Au, Bi, Ag and Pt (all 0.1–10 ppm range). Cu:Pt and Cu:Au ratios in the inclusions are significantly (up to 4 log units) lower than in the host alteration veins and adjacent massive sulfide ore veins, suggesting either (1) early Cu loss from the volatiles by chalcopyrite precipitation or (2) enhanced Au and Pt solubilities relative to Cu at the temperatures of entrapment. Concentration ratios between coexisting brine and CH₄ inclusions are lower for Cu, Au, Bi and Ag than for other elements (Na, Ca, Fe, Mn, Zn, Pb) indicating that during interaction with the brine, the hydrocarbon phase was enriched in ore metals. The high concentrations of ore metals in hydrocarbon, brine and halide melt phases confirm that both aqueous and non-aqueous volatiles were carriers of precious metals in the Sudbury environment over a wide range of temperatures. Volatile evolution and magmatic sulfide differentiation were clearly part of a single, continuous process in the Sudbury footwall. The exsolution of H₂O-poor volatiles from fractionated sulfide liquid may have been a principal mechanism controlling the final distribution of PGE and Au in the footwall ore systems. The study reports the first measurements of precious metal concentrations in fluid inclusions from a magmatic Ni–Cu–PGE environment (the Sudbury district). Electronic Supplementary Material Supplementary material is available for this article at     

Chemical evolution of seawater during the Phanerozoic: Implications from the record of marine evaporites

The chemical evolution of seawater during the Phanerozoic is still a matter of debate. We have assembled and critically analyzed the available data for the composition of fluid inclusions in marine halite and for the mineralogy of marine evaporites. The composition of fluid inclusions in primary marine halite reveals two major long-term cycles in the chemistry of seawater during the past 600 myr. The concentration of Mg²⁺, Ca²⁺, and SO₄²⁻ has varied quite dramatically. The Mg²⁺ concentration in seawater during most of the early Paleozoic and Jurassic to Cretaceous was as low as 30 to 40 mmol/kg H₂O; it reached maximum values ≥50 mmol/kg H₂O during the Late Neoproterozoic and Permian. The Ca²⁺ concentration in seawater during the Phanerozoic has reached maximum values two to three times greater than the concentration in seawater today (10.6 mmol/kg H₂O), whereas SO₄²⁻ concentrations may have been as low as 5 to 10 mmol/kg H₂O (a third to a fifth of the modern value) during the Jurassic and Early Paleozoic. The Mg²⁺/Ca²⁺ ratio in seawater ranged from 1 to 1.5 during the early to middle Paleozoic and Jurassic-Cretaceous to a near-modern value of 5.2 during the Late Neoproterozoic and Permian. This change in seawater Mg²⁺/Ca²⁺ ratio is consistent with the notion of alternating “calcite-aragonite seas” recorded in oölites and marine carbonate cements.Several models have been proposed to explain the chemical evolution of seawater. These have invoked significant changes in one or more of the major geochemical processes that control the composition of seawater. The pattern and magnitude of the variations in the composition of seawater proposed in this study are similar to those proposed elsewhere that suggest that seawater fluxes through midocean ridges have played a major role in the evolution of seawater during the past 600 myr. Two Phanerozoic supercycles of the Earth’s exogenic processes were recognized in the literature that are caused by mantle convection and plate activity. The composition of seawater has apparently undergone dramatic secular changes in phase with these supercycles and as a consequence of biological evolution. Analyses of fluid inclusions containing unevaporated seawater and a better understanding of the processes that affect the composition of seawater are needed to refine our understanding of the history of Phanerozoic seawater.