Thermal evolution of n- and iso-alkanes in oils. Part 1: Pyrolysis model for a mixture of 78 alkanes (C1-C32) including 13,206 free radical reactions

A mechanistic model consisting of 13,206 lumped free radical reactions has been developed to describe the thermal evolution of a mixture of 78 alkanes: all n-alkanes from C₁ to C₃₂ and 46 branched alkane model compounds from C₄ to C₃₂. The mixture was meant to represent the major part of the saturated fraction of petroleum. The rate constants used are available from the literature. The lumping together procedure is described and the model validated on the basis of several experimental results from the literature and relating to pure alkanes. The model is also compared to the saturated fraction obtained from pyrolysis of Elgin oil at 372 °C for up to 1000 h. The cracking global activation energy of n-C₁₅ as well as iso-C₁₅ is close to 69 kcal/mol in the range 200-350 °C. The implications of the model for geological reservoirs will be discussed in a following paper.     

The separate production of H2S from the thermal reaction of hydrocarbons with magnesium sulfate and sulfur: Implications for thermal sulfate reduction

The yields and stable C and H isotopic composition of gaseous products from the reactions of pure n-C₂₄ with (1) MgSO₄; and (2) elemental S in sealed Au-tubes at a series of temperatures over the range 220–600 °C were monitored to better resolve the reaction mechanisms. Hydrogen sulfide formation from thermochemical sulfate reduction (TSR) of n-C₂₄ with MgSO₄ was initiated at 431 °C, coincident with the evolution of C₂–C₅ hydrocarbons. Whereas the yields of H₂S increased progressively with pyrolysis temperature, the hydrocarbon yields decreased sharply above 490 °C due to subsequent S consumption. Ethane and propane were initially very ¹³C depleted, but became progressively heavier with pyrolysis temperature and were more ¹³C enriched than the values of a control treatment conducted on just n-C₂₄ above 475 °C. TSR of MgSO₄ also led to progressively higher concentrations of CO₂ showing relatively low δ¹³C values, possibly due to input of isotopically light CO₂ derived from gaseous hydrocarbon oxidation (e.g., more depleted CH₄).     

Influence of physiology and climate on δD of leaf wax n-alkanes from C3 and C4 grasses

We measured hydrogen isotope compositions (δD) of high-molecular-weight n-alkanes (C₂₇-C₃₃) from grasses grown in greenhouses and collected from the US Great Plains. In both cases, n-alkanes from C₄ grasses are enriched in D by more than 20‰ relative to those from C₃ grasses. The apparent enrichment factor (ε_C29-GW) between C₂₉n-alkane and greenhouse water is −165 ± 12‰ for C₃ grasses and −140 ± 15‰ for C₄ grasses. For samples from the Great Plains, δD values of C₂₉n-alkanes range from −280 to −136‰, with values for C₄ grasses ca. 21‰ more positive than those for C₃ grasses from the same site. Differences in C₃ and C₄ grass n-alkane δD values are consistent with the shorter interveinal distance in C₄ grass leaves, and greater back-diffusion of enriched water from stomata to veins, than in C₃ grass leaves. Great Plains’ grass n-alkane isotopic ratios largely reflect precipitation δD values. However, the offset or apparent fractionation between n-alkanes and precipitation is not uniform and varies with annual precipitation and relative humidity, suggesting climatic controls on lipid δD values. The dryer sites exhibit smaller absolute apparent fractionation indicative of D-enrichment of source waters through transpiration and/or soil evaporation. To explore the relationship between climate and n-alkane δD values, we develop three models. (1) The ‘direct analog’ model estimates δD_C29 values simply by applying the apparent enrichment factors, ε_C29-GW, observed in greenhouse grasses to precipitation δD values from the Great Plains. (2) The ‘leaf-water’ model uses a Craig-Gordon model to estimate transpirational D-enrichment for both greenhouse and field sites. The transpiration-corrected enrichment factors between C₂₉ and bulk leaf-water, ε_C29-GW, calculated from the greenhouse samples (−181‰ for C₃ and −157‰ for C₄) are applied to estimate δD_C29 values relative to modeled bulk leaf-water δD values. (3) The ‘soil- and leaf-water’ model estimates the combined effects of soil evaporation, modeled by analogy with a flow-through lake, and transpiration on δD_C29 values. Predictions improve with the addition of the explicit consideration of transpiration and soil evaporation, indicating that they are both important processes in determining plant lipid δD values. D-enrichment caused by these evaporative processes is controlled by relative humidity, suggesting that important climatic information is recorded in leaf wax n-alkane δD values. Calibration studies such as this one provide a baseline for future studies of plant-water-deuterium systematics and form the foundation for interpretation of plant wax hydrogen isotope ratios as a paleo-aridity proxy.     

Molecular H2O in armenite, BaCa2Al6Si9O30·2H2O, and epididymite, Na2Be2Si6O15·H2O: Heat capacity, entropy and local-bonding behavior of confined H2O in microporous silicates

Armenite, ideal formula BaCa₂Al₆Si₉O₃₀·2H₂O, and its dehydrated analog BaCa₂Al₆Si₉O₃₀ and epididymite, ideal formula Na₂Be₂Si₆O₁₅·H₂O, and its dehydrated analog Na₂Be₂Si₆O₁₅ were studied by low-temperature relaxation calorimetry between 5 and 300 K to determine the heat capacity, C_p, behavior of their confined H₂O. Differential thermal analysis and thermogravimetry measurements, FTIR spectroscopy, electron microprobe analysis and powder Rietveld refinements were undertaken to characterize the phases and the local environment around the H₂O molecule.The determined structural formula for armenite is Ba_0.88(0.01)Ca_1.99(0.02)Na_0.04(0.01)Al_5.89(0.03)Si_9.12(0.02)O₃₀·2H₂O and for epididymite Na_1.88(0.03)K_0.05(0.004)Na_0.01(0.004)Be_2.02(0.008)Si_6.00(0.01)O₁₅·H₂O. The infrared (IR) spectra give information on the nature of the H₂O molecules in the natural phases via their H₂O stretching and bending vibrations, which in the case of epididymite only could be assigned. The powder X-ray diffraction data show that armenite and its dehydrated analog have similar structures, whereas in the case of epididymite there are structural differences between the natural and dehydrated phases. This is also reflected in the lattice IR mode behavior, as observed for the natural phases and the H₂O-free phases. The standard entropy at 298 K for armenite is S° = 795.7 ± 6.2 J/mol K and its dehydrated analog is S° = 737.0 ± 6.2 J/mol K. For epididymite S° = 425.7 ± 4.1 J/mol K was obtained and its dehydrated analog has S° = 372.5 ± 5.0 J/mol K. The heat capacity and entropy of dehydration at 298 K are Δ = 3.4 J/mol K and ΔS^rxn = 319.1 J/mol K and Δ = −14.3 J/mol K and ΔS^rxn = 135.7 J/mol K for armenite and epididymite, respectively. The H₂O molecules in both phases appear to be ordered. They are held in place via an ion-dipole interaction between the H₂O molecule and a Ca cation in the case of armenite and a Na cation in epididymite and through hydrogen-bonding between the H₂O molecule and oxygen atoms of the respective silicate frameworks. Of the three different H₂O phases ice, liquid water and steam, the C_p behavior of confined H₂O in both armenite and epididymite is most similar to that of ice, but there are differences between the two silicates and from the C_p behavior of ice. Hydrogen-bonding behavior and its relation to the entropy of confined H₂O at 298 K is analyzed for various microporous silicates.The entropy of confined H₂O at 298 K in various silicates increases approximately linearly with increasing average wavenumber of the OH-stretching vibrations. The interpretation is that decreased hydrogen-bonding strength between a H₂O molecule and the silicate framework, as well as weak ion-dipole interactions, results in increased entropy of H₂O. This results in increased amplitudes of external H₂O vibrations, especially translations of the molecule, and they contribute strongly to the entropy of confined H₂O at T < 298 K.     

A new application of headspace single-drop microextraction technique for compound specific carbon isotopic determination of gasoline range hydrocarbons

In this study, headspace single-drop microextraction (HS-SDME) coupled with gas chromatography-isotope ratio mass spectrometry (GC-IRMS), was employed to determine compound specific carbon isotopic values (δ¹³C) of gasoline range hydrocarbons. The reproducibility of the method was found to be satisfactory. By comparison with the δ¹³C values of the twelve target compounds determined using direct injection of their n-C₁₆ solution, no obvious isotopic fractionation was observed during the HS-SDME procedures. Some parameters that could affect the carbon isotopic fractionation, such as ionic strength of working solutions and inlet split ratio, were examined. The results also suggest that these factors had no significant effect on the carbon isotopic determination of gasoline range hydrocarbons. The application of HS-SDME to a crude oil sample proved that this method could be a promising tool for the determination of carbon isotopic values of gasoline range hydrocarbons in oils or aqueous samples.     

Na2O solubility in CaO-MgO-SiO2 melts

The sodium solubility in silicate melts in the CaO-MgO-SiO₂ (CMS) system at 1400 °C has been measured by using a closed thermochemical reactor designed to control alkali metal activity. In this reactor, Na_(g) evaporation from a Na₂O-xSiO₂ melt imposes an alkali metal vapor pressure in equilibrium with the molten silicate samples. Because of equilibrium conditions in the reactor, the activity of sodium-metal oxide in the molten samples is the same as that of the source, i.e., aNa₂O_(sample) = aNa₂O_(source). This design also allows to determine the sodium oxide activity coefficient in the samples. Thirty-three different CMS compositions were studied. The results show that the amount of sodium entering from the gas phase (i.e., Na₂O solubility) is strongly sensitive to silica content of the melt and, to a lesser extent, the relative amounts of CaO and MgO. Despite the large range of tested melt compositions (0 < CaO and MgO < 40; 40 < SiO₂ < 100; in wt%), we found that Na₂O solubility is conveniently modeled as a linear function of the optical basicity (Λ) calculated on a Na-free basis melt composition. In our experiments, γNa₂O_(sample) ranges from 7 × 10⁻⁷ to 5 × 10⁻⁶, indicating a strongly non-ideal behavior of Na₂O solubility in the studied CMS melts (γNa₂O_(sample) ? 1). In addition to showing the effect of sodium on phase relationships in the CMS system, this Na₂O solubility study brings valuable new constraints on how melt structure controls the solubility of Na in the CMS silicate melts. Our results suggest that Na₂O addition causes depolymerization of the melt by preferential breaking of Si-O-Si bonds of the most polymerized tetrahedral sites, mainly Q⁴.     

Mathematical modeling of the aquatic macrophyte inputs of mid-chain n-alkyl lipids to lake sediments: Implications for interpreting compound specific hydrogen isotopic records

We present a systematic study of chain-length distributions and D/H ratios of n-alkyl lipids (both n-alkanes and n-alkanoic acids) in a wide range of terrestrial and aquatic plants around and in Blood Pond, Massachusetts, USA. The primary goal is to establish a model to quantitatively assess the aquatic plant inputs of the mid-chain length n-alkyl lipids to lake sediments and to determine the average hydrogen isotopic ratios of these lipids in different plants. Our results show that middle-chain n-alkyl lipids (C₂₁-C₂₃n-alkanes and C₂₀-C₂₄n-alkanoic acids) are exceptionally abundant in floating and submerged aquatic plants, in contrast to the dominance of long-chain n-alkyl lipids (C₂₇-C₃₁n-alkanes and C₂₆-C₃₂n-alkanoic acids) in other plant types, which are consistent with previously published data from Mountain Kenya and the Tibetan Plateau. Combining available data in different environmental settings allows us to establish statistically robust model distributions of n-alkyl lipids in floating/submerged macrophytes relative to other plant types. Based on the model distributions, we established a multi-source mixing model using a linear algebra approach, in order to quantify the aquatic inputs of mid-chain n-alkyl lipids in lake sediments. The results show that ∼97% of the mid-chain n-alkyl lipids (C₂₃n-alkane and C₂₂n-acid (behenic acid)) in Blood Pond sediments are derived from floating and submerged macrophytes. In addition, D/H ratios of C₂₂n-acid and C₂₃n-alkane in the floating and submerged plants from Blood Pond display relatively narrow ranges of variation (−161 ± 16‰ and −183 ± 18‰, respectively). Our study demonstrates that mid-chain n-alkyl lipids such as C₂₃n-alkane and C₂₂n-acid could be excellent recorders of past lake water isotopic ratios in lakes with abundant floating and submerged macrophyte inputs.     

The role of metal sulfates in thermochemical sulfate reduction (TSR) of hydrocarbons: Insight from the yields and stable carbon isotopes of gas products

The mechanism of thermochemical sulfate reduction (TSR) was investigated by separately heating n-C₂₄ with three different sulfates (CaSO₄, Na₂SO₄, MgSO₄) in sealed gold tubes at 420 °C and measuring the stable carbon isotope values of hydrocarbon (C₁-C₅) and non-hydrocarbon (CO₂) products. Extensive TSR was observed with the MgSO₄ reactant as reflected by increasing concentrations of H₂S, ¹³C depleted CO₂ and relatively low concentrations of H₂ (compared to the control). H₂S yields were already very high at the first monitoring time (12 h) when the temperature had just reached 420 °C, suggesting that TSR had commenced well prior to this temperature. Only trace amounts of n-C₂₄ and secondary C₃-C₅ alkanes were detected at 12 h, reflecting the efficient TSR utilization of the reactant and lower molecular weight alkane products. Ethane levels were still relatively high at 12 h, but declined thereafter as it was subject to TSR in the absence of higher molecular weight alkanes which had already been utilized. Methane yields were consistently high throughout the 48 h MgSO₄ treatment. The temporal decrease in the concentrations of alkanes available for TSR may also contribute to the sharp enhancement of CO₂ after 36 h. Absence or dampening of the molecular and isotopic trends of MgSO₄ TSR was observed with Na₂SO₄ and CaSO₄ respectively, directly reflecting the levels of TSR reached using these sulfate treatments.For all treatments, the δ¹³C values of C_1-5n-alkanes showed an increase with both molecular weight and treatment time. MgSO₄ TSR led to a 5-10‰ increase in the δ¹³C values of the C₁-C₅ hydrocarbons and a 20‰ decrease in the δ¹³C value of CO₂. The significant ¹³C depletion of the CO₂ may be due to co-production of ¹³C enriched MgCO₃, although this remains unproven as the δ¹³C of MgCO₃ was not measured. The difference in the δ¹³C values of ethane and propane (Δδ¹³C_EP) increased in magnitude with the degree of TSR, and this trend could be used to help evaluate the occurrence and extent of TSR in subsurface gas reservoirs.     

Residence times for protons bound to three oxygen sites in the AlO4Al12(OH)24(H2O)12 polyoxocation

Although, the kinetic reactivity of a mineral surface is determined, in part, by the rates of exchange of surface-bound oxygens and protons with bulk solution, there are no elementary rate data for minerals. However, such kinetic measurements can be made on dissolved polynuclear clusters, and here we report lifetimes for protons bound to three oxygen sites on the AlO₄Al₁₂(OH)₂₄(H₂O)₁₂⁷⁺ (Al₁₃) molecule, which is a model for aluminum-hydroxide solids in water. Proton lifetimes were measured using ¹H NMR at pH ∼ 5 in both aqueous and mixed solvents. The ¹H NMR peak for protons on bound waters (η-H₂O) lies near 8 ppm in a 2.5:1 mixture of H₂O/acetone-d₆ and broadens over the temperature range −20 to −5 °C. Extrapolated to 298 K, the lifetime of a proton on a η-H₂O is τ²⁹⁸ ∼ 0.0002 s, which is surprisingly close to the lifetime of an oxygen in the η-H₂O (∼0.0009 s), but in the same general range as lifetimes for protons on fully protonated monomer ions of trivalent metals (e.g., Al(H₂O)₆³⁺). The lifetime is reduced somewhat by acid addition, indicating that there is a contribution from the partly deprotonated Al₁₃ molecule in addition to the fully protonated Al₁₃ at self-buffered pH conditions. Proton lifetimes on the two distinct sets of hydroxyls bridging two Al(III) (μ₂-OH) differ substantially and are much shorter than the lifetime of an oxygen at these sites. The average lifetimes for hydroxyl protons were measured in a 2:1 mixture of H₂O/dmso-d₆ over the temperature range 3.7-95.2 °C. The lifetime of a hydrogen on one of the μ₂-OH was also measured in D₂O. The τ²⁹⁸ values are ∼0.013 and ∼0.2 s in the H₂O/dmso-d₆ solution and the τ²⁹⁸ value for the μ₂-OH detectable in D₂O is τ²⁹⁸ ∼ 0.013 s. The ¹H NMR peak for the more reactive μ₂-OH broadens slightly with acid addition, indicating a contribution from an exchange pathway that involves a proton or hydronium ion. These data indicate that surface protons on minerals will equilibrate with near-surface waters on the diffusional time scale.     

Use of biomarkers indices in a sediment core to evaluate potential pollution sources in a subtropical reservoir in Brazil

The presence of PAHs, n-alkanes, pristane, and phytanes in core sediment from the Vossoroca reservoir (Parana, southern Brazil) was investigated. The total concentration of the 16 PAHs varied from 15.5 to 1646 μg kg⁻¹. Naphthalene was present in all layers (3.34–74.0 μg kg⁻¹). The most abundant and dominant n-alkanes were n-C₁₅ and n-C₃₆, with average concentrations of 198.1 ± 46.8 and 522.9 ± 167.7 μg kg⁻¹, respectively. Lighter n-alkanes were distributed more evenly through the layers and showed less variation, specially n-C₉, n-C₁₂, and n-C₁₈, with average concentrations of 14.6 ± 3.0, 31.6 ± 1.9, and 95.0 ± 5.2 μg kg⁻¹, respectively; heavier n-alkanes were more unevenly distributed.     

Application of the Pitzer ion interaction model to isopiestic data for the Fe2(SO4)3-H2SO4-H2O system at 298.15 and 323.15 K

Recent isopiestic studies of the Fe₂(SO₄)₃-H₂SO₄-H₂O system at 298.15 K are represented with an extended version of Pitzer’s ion interaction model. The model represents osmotic coefficients for aqueous {(1 − y)Fe₂(SO₄)₃ + yH₂SO₄} mixtures from 0.45 to 3.0 m at 298.15 K and 0.0435 ? y ? 0.9370. In addition, a slightly less accurate representation of a more extended molality range to 5.47 m extends over the same y values, translating to a maximum ionic strength of 45 m. Recent isopiestic data for the system at 323.15 K are represented with the extended Pitzer model over a limited range in molality and solute fraction. These datasets are also represented with the usual “3-parameter” version of Pitzer’s model so that it may be incorporated in geochemical modeling software, but is a slightly less accurate representation of thermodynamic properties for this system. Comparisons made between our ion interaction model and available solubility data display partial agreement for rhomboclase and significant discrepancy for ferricopiapite. The comparisons highlight uncertainty remaining for solubility predictions in this system as well as the need for additional solubility measurements for Fe³⁺-bearing sulfate minerals. The resulting Pitzer ion interaction models provide an important step toward an accurate and comprehensive representation of thermodynamic properties in this geochemically important system.     

Possible ecological and environmental significance of the predominance of even-carbon number C20-C30n-alkanes

Even-carbon number C₂₂-C₂₈n-alkanes, are present in higher concentrations than their C₂₁-C₂₉ odd-carbon number homologs in an organic-rich rock of Mississippian (ca. 310 × 10⁶-yr-old) age, which contains algae-like remnants. Environmental indicators associated with this organic-rich rock permitted environmental comparison with ancient and modern environments that contain n-alkanes which exhibit similar distributions. The results of these comparisons suggest that even-carbon number n-alkanes are preferentially produced in highly saline, carbonate environments where aerobic and anaerobic bacteria have subsisted on the remains of blue-green algae.     

Solubility of grossular, Ca3Al2Si3O12, in H2O-NaCl solutions at 800 °C and 10 kbar, and the stability of garnet in the system CaSiO3-Al2O3-H2O-NaCl

The solubility and stability of synthetic grossular were determined at 800 °C and 10 kbar in NaCl-H₂O solutions over a large range of salinity. The measurements were made by evaluating the weight losses of grossular, corundum, and wollastonite crystals equilibrated with fluid for up to one week in Pt capsules and a piston-cylinder apparatus. Grossular dissolves congruently over the entire salinity range and displays a large solubility increase of 0.0053 to 0.132 molal Ca₃Al₂Si₃O₁₂ with increasing NaCl mole fraction (X_NaCl) from 0 to 0.4. There is thus a solubility enhancement 25 times the pure H₂O value over the investigated range, indicating strong solute interaction with NaCl. The Ca₃Al₂Si₃O₁₂ mole fraction versus NaCl mole fraction curve has a broad plateau between X_NaCl = 0.2 and 0.4, indicating that the solute products are hydrous; the enhancement effect of NaCl interaction is eventually overtaken by the destabilizing effect of lowering H₂O activity. In this respect, the solubility behavior of grossular in NaCl solutions is similar to that of corundum and wollastonite. There is a substantial field of stability of grossular at 800 °C and 10 kbar in the system CaSiO₃-Al₂O₃-H₂O-NaCl. At high Al₂O₃/CaSiO₃ bulk compositions the grossular + fluid field is limited by the appearance of corundum. Zoisite appears metastably with corundum in initially pure H₂O, but disappears once grossular is nucleated. At X_NaCl = 0.3, however, zoisite is stable with corundum and fluid; this is the only departure from the quaternary system encountered in this study. Corundum solubility is very high in solutions containing both NaCl and CaSiO₃: Al₂O₃ molality increases from 0.0013 in initially pure H₂O to near 0.15 at X_NaCl = 0.4 in CaSiO₃-saturated solutions, a >100-fold enhancement. In contrast, addition of Al₂O₃ to wollastonite-saturated NaCl solutions increases CaSiO₃ molality by only 12%. This suggests that at high pH (quench pH is 11-12), the stability of solute Ca chloride and Na-Al ± Si complexes account for high Al₂O₃ solubility, and that Ca-Al ± Si complexes are minor. The high solubility and basic dissolution reaction of grossular suggest that Al may be a very mobile component in calcareous rocks in the deep crust and upper mantle when migrating saline solutions are present.     

An experimental study of the solubility of baddeleyite (ZrO2) in fluoride-bearing solutions at elevated temperature

The solubility of baddeleyite (ZrO₂) and the speciation of zirconium have been investigated in HF-bearing aqueous solutions at temperatures up to 400 °C and pressures up to 700 bar. The data obtained suggest that in HF-bearing solutions zirconium is transported mainly in the form of the hydroxyfluoride species ZrF(OH)₃° and ZrF₂(OH)₂°. Formation constants determined for these species (Zr⁴⁺ + nF⁻ + mOH⁻ = ZrF_n(OH)_m°) range from 43.7 at 100 °C to 46.41 at 400 °C for ZrF(OH)₃°, and from 37.25 at 100 °C to 43.88 at 400 °C for ZrF₂(OH)₂°.Although the solubility of ZrO₂ is retrograde with respect to temperature, the measured concentrations of Zr are orders of magnitude higher than those predicted from theoretical extrapolations based on simple fluoride species (ZrF³⁺-ZrF₆²⁻). Model calculations performed for zircon show that zirconium can be transported by aqueous fluids in concentrations sufficient to account for the concentration of this metal at conditions commonly encountered in fluoride-rich natural hydrothermal systems.     

Solubility of Fe-ettringite (Ca6[Fe(OH)6]2(SO4)3 · 26H2O)

The solubility of Fe-ettringite (Ca₆[Fe(OH)₆]₂(SO₄)₃ · 26H₂O) was measured in a series of precipitation and dissolution experiments at 20 °C and at pH-values between 11.0 and 14.0 using synthesised material. A time-series study showed that equilibrium was reached within 180 days of ageing. After equilibrating, the solid phases were analysed by XRD and TGA while the aqueous solutions were analysed by ICP-OES (calcium, sulphur) and ICP-MS (iron). Fe-ettringite was found to be stable up to pH 13.0. At higher pH-values Fe-monosulphate (Ca₄[Fe(OH)₆]₂(SO₄) · 6H₂O) and Fe-monocarbonate (Ca₄[Fe(OH)₆]₂(CO₃) · 6H₂O) are formed. The solubilities of these hydrates at 25 °C are:      

Thermodynamic modeling of binary CH4-H2O fluid inclusions

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

Strain response and re-equilibration of CH4-rich synthetic aqueous fluid inclusions in calcite during pressure drops

Aqueous fluids in sedimentary basins often contain dissolved methane, particularly in petroleum environments. PVTX (Pressure-Volume-Temperature-Composition) reconstructions performed using fluid inclusion data are largely based on the assumption that inclusions do not change from the time of trapping until the present. Many authors, however, consider that fluid inclusions can re-equilibrate, particularly in fragile minerals like calcite. In order to understand this re-equilibration phenomenon in the metamorphic domain, previous experiments have been performed under high PT conditions, but few have been performed at low to medium PT conditions such as those associated with sedimentary burial diagenesis, and no previous studies have examined CH₄-bearing aqueous inclusions in calcite.An experimental study of the preservation/modification of CH₄-rich synthetic fluid inclusions in calcite during isothermal decompression was conducted. An autoclave was used for accurate PTX control allowing equilibrium between liquid and vapour in the CH₄-H₂O system. PTX conditions were maintained at four stages of decreasing pressure, with each stage held for 7 days to simulate an isothermal pressure drop. In order of decreasing pressure, the pressure-temperature conditions monitored were 276 ± 10 bar at 180 ± 7 °C, 176 ± 10 bar at 180 ± 7 °C, 76 ± 10 bar at 180 ± 7 °C and 10 ± 3 bar at 180 ± 15 °C. At the end of the experiment, the calcite was recovered and analyzed by microthermometry and Raman microspectroscopy for PTX reconstruction. A careful procedure was adopted to limit re-equilibration of inclusions during analytical procedures. Four types of inclusion shapes and four types of strain patterns were differentiated. Classification of the petrographic strain patterns was carried out. These strain patterns were associated with inclusion stretching and/or leakage regarding CH₄, T_h and P_h compared to experimental conditions. Factors controlling the preservation or acquisition of strain patterns included the initial shape and size of the inclusion, and the pressure differential (ΔP) between the confining pressure (P_cf) and the internal pressure (P_i) within the inclusion. Most fluid inclusions seemed to be trapped during the first 7 days of the experiment, although few (4%) of these preserved the initial PT conditions of 276 ± 10 bar, whereas 8% preserved the second and third run of PT conditions. Overall, the majority of inclusions (88%) did not reflect accurately the PTX trapping conditions. A petrographic guide to the inclusions is presented here that allows strain identification for PVT reconstructions. Re-equilibration patterns and evidence for preferential methane leakage from aqueous inclusions in calcite are important findings revealed by this study, and may be useful for the reconstruction of post-trapping events in investigations of natural samples, and in other experiments using synthetic inclusions in calcite.     

Indoor and outdoor urban atmospheric CO2: Stable carbon isotope constraints on mixing and mass balance

From July to November 2009, concentrations of CO₂ in 78 samples of ambient air collected in 18 different interior spaces on a university campus in Dallas, Texas (USA) ranged from 386 to 1980 ppm. Corresponding δ¹³C values varied from −8.9‰ to −19.4‰. The CO₂ from 22 samples of outdoor air (also collected on campus) had a more limited range of concentrations from 385 to 447 ppm (avg. = 408 ppm), while δ¹³C values varied from −10.1‰ to −8.4‰ (avg.=-9.0‰). In contrast to ambient indoor and outdoor air, the concentrations of CO₂ exhaled by 38 different individuals ranged from 38,300 to 76,200 ppm (avg. = 55,100 ppm), while δ¹³C values ranged from −24.8‰ to −17.7‰ (avg. = −21.8‰). The residence times of the total air in the interior spaces of this study appear to have been on the order of 10 min with relatively rapid approaches (∼30 min) to steady-state concentrations of ambient CO₂ gas. Collectively, the δ¹³C values of the indoor CO₂ samples were linearly correlated with the reciprocal of CO₂ concentration, exhibiting an intercept of −21.8‰, with r² = 0.99 and p < 0.001 (n = 78). This high degree of linearity for CO₂ data representing 18 interior spaces (with varying numbers of occupants), and the coincidence of the intercept (−21.8‰) with the average δ¹³C value for human-exhaled CO₂ demonstrates simple mixing between two inputs: (1) outdoor CO₂ introduced to the interior spaces by ventilation systems, and (2) CO₂ exhaled by human occupants of those spaces. If such simple binary mixing is a common feature of interior spaces, it suggests that the intercept of a mixing line defined by two data points (CO₂ input from the local ventilation system and CO₂ in the ambient air of the room) could be a reasonable estimate of the average δ¹³C value of the CO₂ exhaled by the human occupants. Thus, such indoor spaces appear to constitute effective “sample vessels” for collection of CO₂ that can be used to determine the average proportions of C₃ and C₄-derived C in the diets of the occupants. For the various groups occupying the rooms sampled in this study, C₄-derived C appears to have constituted ∼40% of the average diet.     

Carboxylic acid distribution patterns of temperate C3 and C4 crops

Agricultural grasses cover a major part of the land surface in temperate agro-ecosystems and contribute significantly to the formation of soil organic matter. Crop-derived lipids are assumed to be responsible for fast carbon turnover in soils. Differences in lipid distribution patterns between crops following C₃ and C₄ photosynthesis pathways have rarely been described, but could be useful for source apportionment of crop-derived input into soils or sediments. The distribution of long chain n-carboxylic acids (C₂₂, C₂₄, C₂₆) reveals significant differences between crop plants following either the C₃ or the C₄ photosynthetic carbon fixation pathway. The plant compartments leaves, stems and roots of C₄ plants contain relatively large proportions (> 40%) of n-C₂₄ carboxylic acid when compared to C₃ plants. These reveal larger relative proportions of n-C₂₂ and n-C₂₆ acids, whose relative abundance is subject to change between different plant compartments and during the growing season. The carboxylic acid ratio [CAR = n-C₂₄/(n-C₂₂ + n-C₂₆) carboxylic acids] provides distinct ratios for C₄ (> 0.67) and C₃ crops (< 0.67) and can thus be used as a molecular marker for the differentiation of crop plant biomass. In combination with the bulk stable carbon isotopic composition (δ¹³C) the CAR can be used as a tool for the estimation of the C₄ derived carbon proportion in soils or sediments.     

Hydrogen isotopic fractionation in lipid biosynthesis by H2-consuming Desulfobacterium autotrophicum

We report hydrogen isotopic fractionations between water and fatty acids of the sulfate-reducing bacterium Desulfobacterium autotrophicum. Pure cultures were grown in waters with deuterium (D) contents that were systematically varied near the level of natural abundance (−37‰ ? δD ? 993‰). H₂ of constant hydrogen isotope (D/H) ratio was supplied to the cultures. The D/H ratios of water, H₂, and specific fatty acids were measured by isotope-ratio mass spectrometry. The results demonstrate that D. autotrophicum catalyzes hydrogen isotopic exchange between water and H₂, and this reaction is conclusively shown to approach isotopic equilibrium. In addition, variation in the D/H ratio of growth water accounts for all variation in the hydrogen isotopic composition of fatty acids. The D/H ratios of fatty acids from cultures grown on H₂/CO₂ are compared with those from a separate set of cultures grown on D-enriched formate, an alternative electron donor. This comparison rules out H₂ as a significant source of fatty acid hydrogen. Grown on either H₂/CO₂ or formate, D. autotrophicum produces fatty acids in which all hydrogen originates from water. For specific fatty acids, biosynthetic fractionation factors are mostly in the range 0.60 ? α_FA-water ? 0.70; the 18:0 fatty acid exhibits a lower fractionation factor of 0.52. The data show that α_FA-water generally increases with length of the carbon chain from C₁₄ to C₁₇ among both saturated and unsaturated fatty acids. These results indicate a net fractionation associated with fatty acid biosynthesis in D. autotrophicum that is slightly smaller than in another H₂-consuming bacterium (Sporomusa sp.), but much greater than in most photoautotrophs.