Petrogenesis of olivine-phyric shergottite Yamato 980459, revisited

Primitive magmas provide critical information on mantle sources, but most Martian meteorites crystallized from fractionated melts. An olivine-phyric shergottite, Yamato 980459 (Y-980459), has been interpreted to represent a primary melt, because its olivine megacrysts have magnesian cores (Fo84-86) that appear to be in equilibrium with the Y-980459 whole-rock composition based on Fe-Mg partitioning. However, crystal size distribution (CSD) plots for Y-980459 olivines show a size gap, suggesting a cumulus origin for some megacrysts. Because melting experiments using the Y-980459 whole-rock composition have been used to infer the thermal structure and volatile contents of the Martian mantle, the interpretation that this rock is primitive should be scrutinized.We report major, minor and trace element compositions of Y-980459 olivines and compare them with results from melting experiments (both hydrous and anhydrous) and thermodynamic calculations. Cores of the olivine megacrysts have major and minor element contents identical to those of the most magnesian olivines from the experiments, but they differ slightly from those of thermodynamic calculations. This is probably because the Y-980459 whole-rock composition lies near the limit of the range of liquids used to calibrate these models. The megacryst cores (Fo80-85) exhibit minor and trace element (Mn-Ni-Co-Cr-V) characteristics distinct from other olivines (megacryst rims and groundmass olivines, Fo < 80), implying that the megacryst cores crystallized under more reduced conditions (∼IW + 1).Y-980459 contains pyroxenes with orthopyroxene cores mantled by pigeonite and augite. We also found some reversely zoned pyroxenes that have augite cores (low-Mg#) mantled by orthopyroxenes (high-Mg#), although they are uncommon. These reversely zoned pyroxenes are interpreted to have grown initially as atoll-like crystals with later crystallization filling in the hollow centers, implying disequilibrium crystallization at a moderate cooling rate (3-7 °C/h). The calculated REE pattern of a melt in equilibrium with normally zoned pyroxene is parallel to those of glass and the Y-980459 whole-rock as well as other depleted olivine-phyric shergottites, suggesting that Y-980459 was derived from a depleted mantle reservoir.Considering the CSD patterns of Y-980459 olivines, we propose that the olivine megacrysts are cumulus crystals which probably formed in a feeder conduit by continuous melt replenishment, and the parent melt composition was indistinguishable from the Y-980459 whole-rock with 0-2 wt% of H₂O and 0-5 wt% of CO₂. The final magma pulse entrained these cumulus olivines and then crystallized groundmass olivines and pyroxenes. Although Y-980459 contains small amounts of cumulus olivine (<∼6 vol%), we conclude that the Y-980459 whole-rock composition closely approximates a Martian primary melt composition.     

Nitrogen solubility in basaltic melt. Part I. Effect of oxygen fugacity

The role of the oxygen fugacity on the incorporation of nitrogen in basaltic magmas has been investigated using one atmosphere high temperature equilibration of tholeiitic-like compositions under controlled nitrogen and oxygen partial pressures in the [C-N-O] system. Nitrogen was extracted with a CO₂ laser under high vacuum and analyzed by static mass spectrometry. Over a redox range of 18 oxygen fugacity log units, this study shows that the incorporation of nitrogen in silicate melts follows two different behaviors. For log fO₂ values between −0.7 and −10.7 (the latter corresponding to IW − 1.3), nitrogen dissolves as a N₂ molecule into cavities of the silicate network (physical solubility). Nitrogen presents a constant solubility (Henry’s) coefficient of 2.21 ± 0.53 × 10⁻⁹ mol g⁻¹ atm⁻¹ at 1425°C, identical within uncertainties to the solubility of argon. Further decrease in the oxygen fugacity (log fO₂ between −10.7 and −18 corresponding to the range from IW − 1.3 to IW − 8.3) results in a drastic increase of the solubility of nitrogen by up to 5 orders of magnitude as nitrogen becomes chemically bounded with atoms of the silicate melt network (chemical solubility). The present results strongly suggest that under reducing conditions nitrogen dissolves in silicate melts as N³⁻ species rather than as CN⁻ cyanide radicals. The nitrogen content of a tholeiitic magma equilibrated with N₂ is computed from thermochemical processing of our data set as

     

Experimental study of radon production and transport in an analogue for the Martian regolith

The suggestion that radon could be used as a radioactive tracer of regolith-atmosphere exchanges and as a proxy for subsurface water on Mars, as well as its indirect detection in the Martian atmosphere by the rover Opportunity, have raised the need for a better characterization of its production process and transport efficiency in the Martian soil. More specifically, a proper estimation of radon exhalation rate on Mars requires its emanation factor and diffusion length to be determined. The dependence of the emanation factor as a function of pore water content (at 267 and 293 K) and the dependence of the adsorption coefficient on temperature, specific surface area and nature of the carrier gas (He, He + CO₂) have been measured on a Martian soil analogue (Hawaiian palagonitized volcanic ash, JSC Mars-1), whose radiometric analysis has been performed. An estimation of radon diffusion lengths on Mars is provided and is used to derive a global average emanation factor (2-6.5%) that accounts for the exhalation rate inferred from the ²¹⁰Po surface concentration detected on Martian dust and from the ²¹⁴Bi signal measured by the Mars Odyssey Gamma Ray Spectrometer. It is found to be much larger than emanation factors characterizing lunar samples, but lower than the emanation factor of the palagonite samples obtained under dry conditions. This result probably reflects different degrees of aqueous alteration and could indicate that the emanation factor is also affected by the current presence of pore water in the Martian soil. The rationale of the “radon method” as a technique to probe subsurface water on Mars, and its sensitivity to soil parameters are discussed. These experimental data are useful to perform more detailed studies of radon transport in the Martian atmosphere using Global Climate Models and to interpret neutron and gamma data from Mars Odyssey Gamma Ray Spectrometer.     

Calcite dissolution kinetics in the system CaCO3-H2O-CO2 at high undersaturation

Dissolution rates of limestone covered by a water film open to a CO₂-containing atmosphere are controlled by the chemical composition of the CaCO₃-H₂O-CO₂ solution at the water-mineral interface. This composition is determined by the Ca²⁺-concentration at this boundary, conversion of CO₂ into H⁺ and in the solution, and by diffusional mass transport of the dissolved species from and towards the water-limestone interface. A system of coupled diffusion-reaction equations for Ca²⁺, , and CO₂ is derived. The Ca²⁺ flux rates at the surface of the mineral are defined by the PWP-empirical rate law. These flux rates by the rules of stoichiometry must be equal to the flux rates of CO₂ across the air-water interface. In the solution, CO₂ is converted into H⁺ and . At low water-film thickness this reaction becomes rate limiting. The time dependent diffusion-reaction equations are solved for free drift dissolution by a finite-difference scheme, to obtain the dissolution rate of calcite as a function of the average calcium concentration in the water film. Dissolution rates are obtained for high undersaturation. The results reveal two regimes of linear dissolution kinetics, which can be described by a rate law F = α_i(m_ic_eq − c), where c is the calcium concentration in the water film, c_eq the equilibrium concentration with respect to calcite. For index i = 0, a fast rate law, which here is reported for the first time, is found with α₀ = 3 × 10⁻⁶ m s⁻¹ and m₀ = 0.3. For c > m₀c_eq, a slow rate law is valid with α₁ = 3 × 10⁻⁷ m  s⁻¹ and m₁ = 1, which confirms earlier work. The numbers given above are valid for film thickness of several tenths of a millimetre and at 20 °C. These rates are proven experimentally, using a flat inclined limestone plate covered by a laminar flowing water film injected at an input point with known flow rate Q and calcium concentration. From the concentration measured after flow distance x the dissolution rates are determined. These experiments have been performed at a carbon-dioxide pressure of 0.00035 atm and also of 0.01 atm. The results are in good agreement to the theoretical predictions.     

Raman spectroscopic measurements of CO2 density: Experimental calibration with high-pressure optical cell (HPOC) and fused silica capillary capsule (FSCC) with application to fluid inclusion observations

Raman spectroscopy is a powerful method for the determination of CO₂ densities in fluid inclusions, especially for those with small size and/or low fluid density. The relationship between CO₂ Fermi diad split (Δ, cm⁻¹) and CO₂ density (ρ, g/cm³) has been documented by several previous studies. However, significant discrepancies exist among these studies mainly because of inconsistent calibration procedures and lack of measurements for CO₂ fluids having densities between 0.21 and 0.75 g/cm³, where liquid and vapor phases coexist near room temperature.In this study, a high-pressure optical cell and fused silica capillary capsules were used to prepare pure CO₂ samples with densities between 0.0472 and 1.0060 g/cm³. The measured CO₂ Fermi diad splits were calibrated with two well established Raman bands of benzonitrile at 1192.6 and 1598.9 cm⁻¹. The relationship between the CO₂ Fermi diad split and density can be represented by: ρ = 47513.64243 − 1374.824414 × Δ + 13.25586152 × Δ² − 0.04258891551 × Δ³ (r² = 0.99835, σ = 0.0253 g/cm³), and this relationship was tested by synthetic fluid inclusions and natural CO₂-rich fluid inclusions. The effects of temperature and the presence of H₂O and CH₄ on this relationship were also examined.     

Microbial dissolution of calcite at T = 28 °C and ambient pCO2

This study used batch reactors to quantify the mechanisms and rates of calcite dissolution in the presence and absence of a single heterotrophic bacterial species (Burkholderia fungorum). Experiments were conducted at T = 28°C and ambient pCO₂ over time periods spanning either 21 or 35 days. Bacteria were supplied with minimal growth media containing either glucose or lactate as a C source, NH₄⁺ as an N source, and H₂PO₄⁻ as a P source. Combining stoichiometric equations for microbial growth with an equilibrium mass-balance model of the H₂O-CO₂-CaCO₃ system demonstrates that B. fungorum affected calcite dissolution by modifying pH and alkalinity during utilization of ionic N and C species. Uptake of NH₄⁺ decreased pH and alkalinity, whereas utilization of lactate, a negatively charged organic anion, increased pH and alkalinity. Calcite in biotic glucose-bearing reactors dissolved by simultaneous reaction with H₂CO₃ generated by dissolution of atmospheric CO₂ (H₂CO₃ + CaCO₃ → Ca²⁺ + 2HCO₃⁻) and H⁺ released during NH₄⁺ uptake (H⁺ + CaCO₃ → Ca²⁺ + HCO₃⁻). Reaction with H₂CO₃ and H⁺ supplied ∼45% and 55% of the total Ca²⁺ and ∼60% and 40% of the total HCO₃⁻, respectively. The net rate of microbial calcite dissolution in the presence of glucose and NH₄⁺ was ∼2-fold higher than that observed for abiotic control experiments where calcite dissolved only by reaction with H₂CO₃. In lactate bearing reactors, most H⁺ generated by NH₄⁺ uptake reacted with HCO₃⁻ produced by lactate oxidation to yield CO₂ and H₂O. Hence, calcite in biotic lactate-bearing reactors dissolved by reaction with H₂CO₃ at a net rate equivalent to that calculated for abiotic control experiments. This study suggests that conventional carbonate equilibria models can satisfactorily predict the bulk fluid chemistry resulting from microbe-calcite interactions, provided that the ionic forms and extent of utilization of N and C sources can be constrained. Because the solubility and dissolution rate of calcite inversely correlate with pH, heterotrophic microbial growth in the presence of nonionic organic matter and NH₄⁺ appears to have the greatest potential for enhancing calcite weathering relative to abiotic conditions.     

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.     

Abundance of mass 47 CO2 in urban air, car exhaust, and human breath

Atmospheric carbon dioxide is widely studied using records of CO₂ mixing ratio, δ¹³C and δ¹⁸O. However, the number and variability of sources and sinks prevents these alone from uniquely defining the budget. Carbon dioxide having a mass of 47 u (principally ¹³C¹⁸O¹⁶O) provides an additional constraint. In particular, the mass 47 anomaly (Δ₄₇) can distinguish between CO₂ produced by high temperature combustion processes vs. low temperature respiratory processes. Δ₄₇ is defined as the abundance of mass 47 isotopologues in excess of that expected for a random distribution of isotopes, where random distribution means that the abundance of an isotopologue is the product of abundances of the isotopes it is composed of and is calculated based on the measured ¹³C and ¹⁸O values. In this study, we estimate the δ¹³C (vs. VPDB), δ¹⁸O (vs. VSMOW), δ47, and Δ₄₇ values of CO₂ from car exhaust and from human breath, by constructing ‘Keeling plots’ using samples that are mixtures of ambient air and CO₂ from these sources. δ47 is defined as , where is the R⁴⁷ value for a hypothetical CO₂ whose δ¹³C_VPDB = 0, δ¹⁸O_VSMOW = 0, and Δ₄₇ = 0. Ambient air in Pasadena, CA, where this study was conducted, varied in [CO₂] from 383 to 404 μmol mol⁻¹, in δ¹³C and δ¹⁸O from −9.2 to −10.2‰ and from 40.6 to 41.9‰, respectively, in δ47 from 32.5 to 33.9‰, and in Δ₄₇ from 0.73 to 0.96‰. Air sampled at varying distances from a car exhaust pipe was enriched in a combustion source having a composition, as determined by a ‘Keeling plot’ intercept, of −24.4 ± 0.2‰ for δ¹³C (similar to the δ¹³C of local gasoline), δ¹⁸O of 29.9 ± 0.4‰, δ47 of 6.6 ± 0.6‰, and Δ₄₇ of 0.41 ± 0.03‰. Both δ¹⁸O and Δ₄₇ values of the car exhaust end-member are consistent with that expected for thermodynamic equilibrium at∼200 °C between CO₂ and water generated by combustion of gasoline-air mixtures. Samples of CO₂ from human breath were found to have δ¹³C and δ¹⁸O values broadly similar to those of car exhaust-air mixtures, −22.3 ± 0.2 and 34.3 ± 0.3‰, respectively, and δ47 of 13.4 ± 0.4‰. Δ₄₇ in human breath was 0.76  ± 0.03‰, similar to that of ambient Pasadena air and higher than that of the car exhaust signature.     

Melting of the Indarch meteorite (EH4 chondrite) at 1 GPa and variable oxygen fugacity: Implications for early planetary differentiation processes

In order to derive constraints on planetary differentiation processes, and ultimately the formation of the Earth, it is required to study a variety of meteoritic materials and to investigate their melting relations and elemental partitioning at variable pressures, temperatures, and oxygen fugacities (f_O2). This study reports the first high pressure (HP) and high temperature (HT) investigation of an enstatite chondrite (Indarch). Four series of experiments exploring various f_O2 conditions have been carried out at 1 GPa in a piston-cylinder apparatus using the EH4 chondrite Indarch. We show that temperature and redox conditions have important effects on the phase equilibria of the meteorite: the solidus and liquidus temperatures of the silicate portion increase with decreasing f_O2, and the stability fields of various phases are modified. Olivine and pyroxene are stable around 1.5 log f_O2 unit below the iron-wüstite buffer (IW−1.5), whereas quartz and pyroxene is the stable assemblage under the most reducing conditions, between IW−5.0 and IW−4.0, due to reduction of the silicate. While these changes are occurring in the silicate, the metal gains Si from the silicate, (Fe, Mg, Mn, Ca, Cr)-bearing sulfides are observed at f_O2 less than IW−4, and the partitioning of Ni and Mo are both affected by the presence of Si in Fe-S-C liquids. The f_O2 has also a significant effect on the liquid metal-liquid silicate partitioning behavior of Si and S, two possible light elements in planetary cores, and of the slightly siderophile elements Cr and Mn. With decreasing f_O2, S becomes increasingly lithophile, Si becomes increasingly siderophile, and Cr and Mn both become strongly siderophile and chalcophile. The partitioning behavior of these elements places new constraints on models of core segregation for the Earth and other differentiated bodies.     

Between carbonatite and lamproite—Diamondiferous Torngat ultramafic lamprophyres formed by carbonate-fluxed melting of cratonic MARID-type metasomes

New U-Pb perovskite ages reveal that diamondiferous ultramafic lamprophyre magmas erupted through the Archean crust of northern Labrador and Quebec (eastern Canada) between ca. 610 and 565 Ma, a period of strong rifting activity throughout contiguous Laurentia and Baltica. The observed Torngat carbonate-rich aillikite/carbonatite and carbonate-poor mela-aillikite dyke varieties show a large spread in Sr-Nd-Hf-Pb isotope ratios with pronounced correlations between isotope systems. An isotopically depleted component is identified solely within aillikites (⁸⁷Sr/⁸⁶Sr_i = 0.70323-0.70377; εNd_i = +1.2-+1.8; εHf_i = +1.4-+3.5; ²⁰⁶Pb/²⁰⁴Pb_i = 18.2-18.5), whereas some aillikites and all mela-aillikites range to more enriched isotope signatures (⁸⁷Sr/⁸⁶Sr_i = 0.70388-0.70523; εNd_i = −0.5 to −3.9; εHf_i = −0.6 to −6.0; ²⁰⁶Pb/²⁰⁴Pb_i = 17.8-18.2). These contrasting isotopic characteristics of aillikites/carbonatites and mela-aillikites, along with subtle differences in their modal carbonate, SiO₂, Al₂O₃, Na₂O, Cs-Rb, and Zr-Hf contents, are consistent with two distinctive metasomatic assemblages of different age in the mantle magma source region.Integration of petrologic, geochemical, and isotopic information leads us to propose that the isotopically enriched component originated from a reduced phlogopite-richterite-Ti-oxide dominated source assemblage that is reminiscent of MARID suite xenoliths. In contrast, the isotopically depleted component was derived from a more oxidized phlogopite-carbonate dominated source assemblage. We argue that low-degree CO₂-rich potassic silicate melts from the convective upper mantle were preferentially channelled into an older, pre-existing MARID-type vein network at the base of the North Atlantic craton lithosphere, where they froze to form new phlogopite-carbonate dominated veins. Continued stretching and thinning of the cratonic lithosphere during the Late Neoproterozoic remobilized the carbonate-rich vein material and induced volatile-fluxed fusion of the MARID-type veins and the cold peridotite substrate. Isotopic modelling suggests that only 5-12% trace element contribution from such geochemically extreme MARID-type material is required to produce the observed compositional shift from the isotopically most depleted aillikites/carbonatites towards enriched mela-aillikites.We conclude that cold cratonic mantle lithosphere can host several generations of contrasting vein assemblages, and that each may have formed during past tectonic and magmatic events under distinctively different physicochemical conditions. Although cratonic MARID-type and carbonate-bearing veins in peridotite can be the respective sources for lamproite and carbonatite magmas when present as the sole metasome, their concomitant fusion in a complex source region may give rise to a whole new variety of deep volatile-rich magmas and we suggest that orangeites (formerly Group 2 kimberlites), kamafugites, and certain types of ultramafic lamprophyre are formed in this manner.     

Determination of carbon fractionation factor between aqueous carbonate and CO2(g) in two-direction isotope equilibration

The experiments were conducted in the open CO₂ system to find out the equilibrium fractionation between the carbonate ion and CO₂(g). The existence of isotopic equilibrium was checked using the two-direction approach by passing the CO₂−N₂ gases with different δ¹³C compositions (− 1.5‰ and − 23‰) through the carbonate solution with δ¹³C = − 4.2‰. The Δ_{CO3T²⁻−CO2(g)} equilibrium fractionation is given as 6.03 ± 0.17‰ at 25 °C. Discussion is provided about the significance of carbonate complexing in determination of Δ_{CO3T²⁻−CO2(g)} and Δ_{HCO3T⁻−CO2(g)} fractionations. Finally, an isotope numerical model of flow and kinetics of hydration and dehydroxylation is built to predict the isotopic behaviour of the system with time.     

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.     

Surface heat flow and CO2 emissions within the Ohaaki hydrothermal field,Taupo Volcanic Zone,New Zealand

Carbon dioxide emissions and heat flow have been determined from the Ohaaki hydrothermal field, Taupo Volcanic Zone (TVZ), New Zealand following 20 a of production (116 MW_e). Soil CO₂ degassing was quantified with 2663 CO₂ flux measurements using the accumulation chamber method, and 2563 soil temperatures were measured and converted to equivalent heat flow (W m⁻²) using published soil temperature heat flow functions. Both CO₂ flux and heat flow were analysed statistically and then modelled using 500 sequential Gaussian simulations. Forty subsoil CO₂ gas samples were also analysed for stable C isotopes. Following 20 a of production, current CO₂ emissions equated to 111 ± 6.7 T/d. Observed heat flow was 70 ± 6.4 MW, compared with a pre-production value of 122 MW. This 52 MW reduction in surface heat flow is due to production-induced drying up of all alkali–Cl outflows (61.5 MW) and steam-heated pools (8.6 MW) within the Ohaaki West thermal area (OHW). The drying up of all alkali–Cl outflows at Ohaaki means that the soil zone is now the major natural pathway of heat release from the high-temperature reservoir. On the other hand, a net gain in thermal ground heat flow of 18 MW (from 25 MW to 43.3 ± 5 MW) at OHW is associated with permeability increases resulting from surface unit fracturing by production-induced ground subsidence. The Ohaaki East (OHE) thermal area showed no change in distribution of shallow and deep soil temperature contours despite 20 a of production, with an observed heat flow of 26.7 ± 3 MW and a CO₂ emission rate of 39 ± 3 T/d. The negligible change in the thermal status of the OHE thermal area is attributed to the low permeability of the reservoir beneath this area, which has limited production (mass extraction) and sheltered the area from the pressure decline within the main reservoir. Chemistry suggests that although alkali–Cl outflows once contributed significantly to the natural surface heat flow (∼50%) they contributed little (<1%) to pre-production CO₂ emissions due to the loss of >99% of the original CO₂ content due to depressurisation and boiling as the fluids ascended to the surface. Consequently, the soil has persisted as the major (99%) pathway of CO₂ release to the atmosphere from the high temperature reservoir at Ohaaki. The CO₂ flux and heat flow surveys indicate that despite 20 a of production the variability in location, spatial extent and magnitude of CO₂ flux remains consistent with established geochemical and geophysical models of the Ohaaki Field. At both OHW and OHE carbon isotopic analyses of soil gas indicate a two-stage fractionation process for moderate-flux (>60 g m⁻² d⁻¹) sites; boiling during fluid ascent within the underlying reservoir and isotopic enrichment as CO₂ diffuses through porous media of the soil zone. For high-flux sites (>300 g m⁻² d⁻¹), the δ¹³CO₂ signature (−7.4 ± 0.3‰ OHW and −6.5 ± 0.6‰ OHE) is unaffected by near-surface (soil zone) fractionation processes and reflects the composition of the boiled magmatic CO₂ source for each respective upflow. Flux thresholds of <30 g m⁻² d⁻¹ for purely diffusive gas transport, between 30 and 300 g m⁻² d⁻¹ for combined diffusive–advective transport, and ?300 g m⁻² d⁻¹ for purely advective gas transport at Ohaaki were assigned. δ¹³CO₂ values and cumulative probability plots of CO₂ flux data both identified a threshold of ∼15 g m⁻² d⁻¹ by which background (atmospheric and soil respired) CO₂ may be differentiated from hydrothermal CO₂.     

Diffuse gas emissions at the Ukinrek Maars,Alaska: Implications for magmatic degassing and volcanic monitoring

Diffuse CO₂ efflux near the Ukinrek Maars, two small volcanic craters that formed in 1977 in a remote part of the Alaska Peninsula, was investigated using accumulation chamber measurements. High CO₂ efflux, in many places exceeding 1000 g m⁻² d⁻¹, was found in conspicuous zones of plant damage or kill that cover 30,000–50,000 m² in area. Total diffuse CO₂ emission was estimated at 21–44 t d⁻¹. Gas vents 3-km away at The Gas Rocks produce 0.5 t d⁻¹ of CO₂ that probably derives from the Ukinrek Maars basalt based on similar δ¹³C values (∼−6‰), ³He/⁴He ratios (5.9–7.2 R_A), and CO₂/³He ratios (1–2 × 10⁹) in the two areas. A lower ³He/⁴He ratio (2.7 R_A) and much higher CO₂/³He ratio (9 × 10¹⁰) in gas from the nearest arc-front volcanic center (Mount Peulik/Ugashik) provide a useful comparison. The large diffuse CO₂ emission at Ukinrek has important implications for magmatic degassing, subsurface gas transport, and local toxicity hazards. Gas–water–rock interactions play a major role in the location, magnitude and chemistry of the emissions.     

Negative δO values in Allan Hills 84001 carbonate: Possible evidence for water precipitation on mars

The Martian meteorite ALH84001 contains ∼1% by weight of carbonate formed by secondary processes on the Martian surface or in the shallow subsurface. The major form of this carbonate is chemically and isotopically zoned rosettes which have been well documented elsewhere. This study concentrates upon carbonate regions ∼200 μm across which possess previously unobserved magnesium rich inner cores, interpreted here as rosette fragments, surrounded by a later stage cement containing rare Ca-rich carbonates (up to Ca₈₁Mg₀₇Fe₀₄Mn₀₇) intimately associated with feldspar. High spatial resolution ion probe analyses of Ca-rich carbonate surrounding rosette fragments have δ¹⁸O_V-SMOW values as low as −10‰. These values are not compatible with deposition from a global Martian atmosphere invoked to explain ALH84001 rosettes. The range of δ¹⁸O values are also incompatible with a fluid that has equilibrated with the Martian crust at high temperature or from remobilisation of carbonate of rosette isotopic composition. At Martian atmospheric temperatures, the small CO₂(gas)-CO₂(ice) fractionation makes meteoric CO₂ an unlikely source for −10‰ carbonates. In contrast, closed system Rayleigh fractionation of H₂O can generate δ¹⁸O_H2O −30‰, as observed at high latitudes on Earth. We suggest that atmospheric transport and precipitation of H₂O in a similar fashion to that on Earth provides a source of suitably ¹⁸O depleted water for generation of carbonate with δ¹⁸O_V-SMOW = −10‰.     

Gas seepage from Tokamachi mud volcanoes, onshore Niigata Basin (Japan): Origin, post-genetic alterations and CH4-CO2 fluxes

Methane and CO₂ emissions from the two most active mud volcanoes in central Japan, Murono and Kamou (Tokamachi City, Niigata Basin), were measured in from both craters or vents (macro-seepage) and invisible exhalation from the soil (mini- and microseepage). Molecular and isotopic compositions of the released gases were also determined. Gas is thermogenic (δ¹³C_CH4 from −32.9‰ to −36.2‰), likely associated with oil, and enrichments of ¹³C in CO₂ (δ¹³C_CO2 up to +28.3‰) and propane (δ¹³C_C3H8 up to −8.6‰) suggest subsurface petroleum biodegradation. Gas source and post-genetic alteration processes did not change from 2004 to 2010. Methane flux ranged within the orders of magnitude of 10¹-10⁴ g m⁻² d⁻¹ in macro-seeps, and up to 446 g m⁻² d⁻¹ from diffuse seepage. Positive CH₄ fluxes from dry soil were widespread throughout the investigated areas. Total CH₄ emission from Murono and Kamou were estimated to be at least 20 and 3.7 ton a⁻¹, respectively, of which more than half was from invisible seepage surrounding the mud volcano vents. At the macro-seeps, CO₂ fluxes were directly proportional to CH₄ fluxes, and the volumetric ratios between CH₄ flux and CO₂ flux were similar to the compositional CH₄/CO₂ volume ratio. Macro-seep flux data, in addition to those of other 13 mud volcanoes, supported the hypothesis that molecular fractionation (increase of the “Bernard ratio” C₁/(C₂ + C₃)) is inversely proportional to gas migration fluxes. The CH₄ “emission factor” (total measured output divided by investigated seepage area) was similar to that derived in other mud volcanoes of the same size and activity. The updated global “emission-factor” data-set, now including 27 mud volcanoes from different countries, suggests that previous estimates of global CH₄ emission from mud volcanoes may be significantly underestimated.     

Carbon dioxide degassing and thermal energy release in the Monte Amiata volcanic-geothermal area (Italy)

The quaternary volcanic complex of Mount Amiata is located in southern Tuscany (Italy) and represents the most recent manifestation of the Tuscan Magmatic Province. The region is characterised by a large thermal anomaly and by the presence of numerous CO₂-rich gas emissions and geothermal features, mainly located at the periphery of the volcanic complex. Two geothermal systems are located, at increasing depths, in the carbonate and metamorphic formations beneath the volcanic complex. The shallow volcanic aquifer is separated from the deep geothermal systems by a low permeability unit (Ligurian Unit). A measured CO₂ discharge through soils of 1.8 × 10⁹ mol a⁻¹ shows that large amounts of CO₂ move from the deep reservoir to the surface. A large range in δ¹³C_TDIC (−21.07 to +3.65) characterises the waters circulating in the aquifers of the region and the mass and isotopic balance of TDIC allows distinguishing a discharge of 0.3 × 10⁹ mol a⁻¹ of deeply sourced CO₂ in spring waters. The total natural CO₂ discharge (2.1 × 10⁹ mol a⁻¹) is slightly less than minimum CO₂ output estimated by an indirect method (2.8 × 10⁹ mol a⁻¹), but present-day release of 5.8 × 10⁹ mol a⁻¹ CO₂ from deep geothermal wells may have reduced natural CO₂ discharge. The heat transported by groundwater, computed considering the increase in temperature from the infiltration area to the discharge from springs, is of the same order of magnitude, or higher, than the regional conductive heat flow (>200 mW m⁻²) and reaches extremely high values (up to 2700 mW m⁻²) in the north-eastern part of the study area. Heat transfer occurs mainly by conductive heating in the volcanic aquifer and by uprising gas and vapor along fault zones and in those areas where low permeability cover is lacking. The comparison of CO₂ flux, heat flow and geological setting shows that near surface geology and hydrogeological setting play a central role in determining CO₂ degassing and heat transfer patterns.     

The origin of Zn isotope fractionation in sulfides

Isotope fractionation of Zn between aqueous sulfide, chloride, and carbonate species (Zn²⁺, Zn(HS)₂, , , ZnS(HS)⁻, ZnCl⁺, ZnCl₂, , and ZnCO₃) was investigated using ab initio methods. Only little fractionation is found between the sulfide species, whereas carbonates are up to 1‰ heavier than the parent solution. At pH > 3 and under atmospheric-like CO₂ pressures, isotope fractionation of Zn sulfides precipitated from sulfidic solutions is affected by aqueous sulfide species and the δ⁶⁶Zn of sulfides reflect these in the parent solutions. Under high P_CO2 conditions, carbonate species become abundant. In high P_CO2 conditions of hydrothermal solutions, Zn precipitated as sulfides is isotopically nearly unfractionated with respect to a low-pH parent fluid. In contrast, negative δ⁶⁶Zn down to at least −0.6‰ can be expected in sulfides precipitated from solutions with pH > 9. Zinc isotopes in sulfides and rocks therefore represent a potential indicator of mid to high pH in ancient hydrothermal fluids.     

Modeling CO2 chemistry, δC, and oxidation of organic carbon and methane in sediment porewater: Implications for paleo-proxies in benthic foraminifera

I present a numerical diffusion-advection-reaction model to simulate CO₂ chemistry, δ¹³C, and oxidation of organic carbon and methane in sediment porewater. The model takes into account detailed reaction kinetics of dissolved CO₂ compounds, H₂O, H⁺, OH⁻, boron and sulfide compounds. These reactions are usually assumed to be in local equilibrium, which is shown to be a good approximation in most cases. The model also includes a diffusive boundary layer across which chemical species are transported between bottom water and the sediment-water interface. While chemical concentrations and δ¹³C_TCO2 at these locations are frequently assumed equal, I demonstrate that they can be quite different. In this case, shells of benthic foraminifera do not reflect the desired properties of bottom water, even for species living at the sediment-water interface (z = 0 cm). Environmental conditions recorded in their shells are strongly influenced by processes occurring within the sediment. The model is then applied to settings in the Santa Barbara Basin and at Hydrate Ridge (Cascadia Margin), locations of strong organic carbon and methane oxidation. In contrast to earlier studies, I show that a limited contribution of methane-derived carbon to porewater TCO₂ in the Santa Barbara Basin cannot be ruled out. Simulation of methane venting shows that at oxidation rates greater than , the δ¹³C of porewater TCO₂ at z > 1 cm is depleted by more than 15‰ relative to bottom water. Depletions of this magnitude have not been observed in living benthic foraminifera, even at methane vents with much higher oxidation rates. This suggests that foraminifera at these sites either calcify at very shallow sediment depth or during times when oxidation rates are much lower than ∼50 μmol cm⁻² y⁻¹.