Benzene solubility in ionic liquids: working toward an understanding of liquid clathrate formation

The solubility of benzene in 15 imidazolium, pyrrolidinium, pyridinium, and piperidinium ionic liquids has been determined; the resulting, benzene-saturated ionic liquid solutions, also known as liquid clathrates, were examined with (1) H and (19) F nuclear magnetic resonance spectroscopy to try and understand the molecular interactions that control liquid clathrate formation. The results suggest that benzene interacts primarily with the cation of the ionic liquid, and that liquid clathrate formation (and benzene solubility) is controlled by the strength of the cation-anion interactions, that is, the stronger the cation-anion interaction, the lower the benzene solubility. Other factors that were determined to be important in the final amount of benzene in any given liquid clathrate phase included attractive interactions between the anion and benzene (when significant), and larger steric or free volume demands of the ions, both of which lead to greater benzene solubility.

Is water structure around hydrophobic groups clathrate-like?

The term "clathrate structure" is quantified for solvation of nonpolar groups by enumerating hydrogen-bonded ring sizes both in the solvation shell and through the shell-bulk interface and comparing it to a bulk control using the ST4 water model. For clathrate-like structure to be evident, the distributions along the hydrophobic surface are expected to be dominated by pentagons, with significant depletion of hexagons and larger polygons. While the distribution in this region is indeed distinguished by a large number of pentagons, there are significant contributions from hexagons and larger rings as well. Calculated polygon distributions through the shell-bulk interface indicate that when water structure is highly cooperative along the hydrophobic surface, hydrogen-bonded pathways leading back into bulk are then reduced. These results are qualitatively consistent with the observation that hydrophobicity is proportional to the nonpolar solute surface area.

Molecular computations for reactions and phase transitions: applications to protein stabilization, hydrates and catalysis

In this work we have made significant contributions in three different areas of interest: therapeutic protein stabilization, thermodynamics of natural gas clathrate-hydrates, and zeolite catalysis. In all three fields, using our various computational techniques, we have been able to elucidate phenomena that are difficult or impossible to explain experimentally. More specifically, in mixed solvent systems for proteins we developed a statistical-mechanical method to model the thermodynamic effects of additives in molecular-level detail. It was the first method demonstrated to have truly predictive (no adjustable parameters) capability for real protein systems. We also describe a novel mechanism that slows protein association reactions, called the “gap effect.” We developed a comprehensive picture of methioine oxidation by hydrogen peroxide that allows for accurate prediction of protein oxidation and provides a rationale for developing strategies to control oxidation. The method of solvent accessible area (SAA) was shown not to correlate well with oxidation rates. A new property, averaged two-shell water coordination number (2SWCN) was identified and shown to correlate well with oxidation rates. Reference parameters for the van der Waals Platteeuw model of clathrate-hydrates were found for structure I and structure II. These reference parameters are independent of the potential form (unlike the commonly used parameters) and have been validated by calculating phase behavior and structural transitions for mixed hydrate systems. These calculations are validated with experimental data for both structures and for systems that undergo transitions from one structure to another. This is the first method of calculating hydrate thermodynamics to demonstrate predictive capability for phase equilibria, structural changes, and occupancy in pure and mixed hydrate systems. We have computed a new mechanism for the methanol coupling reaction to form ethanol and water in the zeolite chabazite. The mechanism at 400°C proceeds via stable intermediates of water, methane, and protonated formaldehyde.

The long-term stability of a possible aqueous ammonium sulfate ocean inside Titan

We model the thermal evolution of a subsurface ocean of aqueous ammonium sulfate inside Titan using a parameterized convection scheme. The cooling and crystallization of such an ocean depends on its heat flux balance, and is governed by the pressure-dependent melting temperatures at the top and bottom of the ocean. Using recent observations and previous experimental data, we present a nominal model which predicts the thickness of the ocean throughout the evolution of Titan; after 4.5 Ga we expect an aqueous ammonium sulfate ocean 56 km thick, overlain by a thick (176 km) heterogeneous crust of methane clathrate, ice I and ammonium sulfate. Underplating of the crust by ice I will give rise to compositional diapirs that are capable of rising through the crust and providing a mechanism for cryovolcanism at the surface. We have conducted a parameter space survey to account for possible variations in the nominal model, and find that for a wide range of plausible conditions, an ocean of aqueous ammonium sulfate can survive to the present day, which is consistent with the recent observations of Titan's spin state from Cassini radar data [Lorenz, R.D., Stiles, B.W., Kirk, R.L., Allison, M.D., del Marmo, P.P., Iess, L., Lunine, J.I., Ostro, S.J., Hensley, S., 2008. Science 319, 1649–1651].

A sublimation technique for high-precision measurements of d13CO2 and mixing ratios of CO2 and N2O from air trapped in ice cores

In order to provide high precision stable carbon isotope ratios (δ13CO2 or δ13C of CO2) from small bubbly, partially and fully clathrated ice core samples we developed a new method based on sublimation coupled to gas chromatography-isotope ratio mass spectrometry (GC-IRMS). In a first step the trapped air is quantitatively released from ~30 g of ice and CO2 together with N2O are separated from the bulk air components and stored in a miniature glass tube. In an off-line step, the extracted sample is introduced into a helium carrier flow using a minimised tube cracker device. Prior to measurement, N2O and organic sample contaminants are gas chromatographically separated from CO2. Pulses of a CO2/N2O mixture are admitted to the tube cracker and follow the path of the sample through the system. This allows an identical treatment and comparison of sample and standard peaks. The ability of the method to reproduce δ13C from bubble and clathrate ice is verified on different ice cores. We achieve reproducibilities for bubble ice between 0.05 ‰ and 0.07 ‰ and for clathrate ice between 0.05 ‰ and 0.09 ‰ (dependent on the ice core used). A comparison of our data with measurements on bubble ice from the same ice core but using a mechanical extraction device shows no significant systematic offset. In addition to δ13C, the CO2 and N2O mixing ratios can be volumetrically derived with a precision of 2 ppmv and 8 ppbv, respectively.

A ~ 32-70 K formation temperature range for the ice grains agglomerated by Comet 67P/Churyumov-Gerasimenko

Grand Canonical Monte Carlo simulations are used to reproduce the N₂/CO ratio ranging between 1.7 x 10⁻³ and 1.6 x 10⁻² observed in situ in the Jupiter-family comet 67 P/Churyumov-Gerasimenko (67 P) by the ROSINA mass spectrometer on board the Rosetta spacecraft. By assuming that this body has been agglomerated from clathrates in the protosolar nebula (PSN), simulations are developed using elaborated interatomic potentials for investigating the temperature dependence of the trapping within a multiple-guest clathrate formed from a gas mixture of CO and N₂ in proportions corresponding to those expected for the PSN. By assuming that 67 P agglomerated from clathrates, our calculations suggest the cometary grains must have been formed at temperatures ranging between ~ 31.8 and 69.9 K in the PSN to match the N₂/CO ratio measured by the ROSINA mass spectrometer. The presence of clathrates in Jupiter-family comets could then explain the potential N₂ depletion (factor of up to ~ 87 compared to the protosolar value) measured in 67 P/Churyumov-Gerasimenko.

Pits formation from volatile outgassing on 67P/Churyumov-Gerasimenko

We investigate the thermal evolution of comet 67P/Churyumov-Gerasimenko's subsurface in the Seth_01 region, where active pits have been observed by the ESA/Rosetta mission. Our simulations show that clathrate destabilization and amorphous ice crystallization can occur at depths corresponding to those of the observed pits in a timescale shorter than 67P/Churyumov-Gerasimenko's lifetime in the comet's activity zone in the inner solar system. Sublimation of crystalline ice down to such depths is possible only in the absence of a dust mantle, which requires the presence of dust grains in the matrix small enough to be dragged out by gas from the pores. Our results are consistent with both pits formation via sinkholes or subsequent to outbursts, the dominant process depending on the status of the subsurface porosity. A sealed dust mantle would favor episodic and disruptive outgassing as a result of increasing gas pressure in the pores, while high porosity should allow the formation of large voids in the subsurface due to the continuous escape of volatiles. We finally conclude that the subsurface of 67P/Churyumov-Gerasimenko is not uniform at a spatial scale of similar to 100-200 m.

Methane concentration and delta deuterium of methane from the NGRIP ice core

The causes of past changes in the global methane cycle and especially the role of marine methane hydrate (clathrate) destabilization events are a matter of debate. Here we present evidence from the North Greenland Ice Core Project ice core based on the hydrogen isotopic composition of methane [dD(CH4)] that clathrates did not cause atmospheric methane concentration to rise at the onset of Dansgaard-Oeschger (DO) events 7 and 8. Box modeling supports boreal wetland emissions as the most likely explanation for the interstadial increase. Moreover, our data show that dD(CH4) dropped 500 years before the onset of DO 8, with CH4 concentration rising only slightly. This can be explained by an early climate response of boreal wetlands, which carry the strongly depleted isotopic signature of high-latitude precipitation at that time.

Stable isotope and geochemical record of sediments from the Bering Sea

here is controversy over the role of marine methane hydrates in atmospheric methane concentrations and climate change during the last glacial period. In this study of two sediment cores from the southeast Bering Sea (700 m and 1467 m water depth), we identify multiple episodes during the last glacial period of intense methane flux reaching the seafloor. Within the uncertainty of the radiocarbon age model, the episodes are contemporaneous in the two cores and have similar timing and duration as Dansgaard-Oeschger events. The episodes are marked by horizons of sediment containing 13C-depleted authigenic carbonate minerals; 13C-depleted archaeal and bacterial lipids, which resemble those found in ANME-1 type anaerobic methane oxidizing microbial consortia; and changes in the abundance and species distribution of benthic foraminifera. The similar timing and isotopic composition of the authigenic carbonates in the two cores is consistent with a region-wide increase in the upward flux of methane bearing fluids. This study is the first observation outside Santa Barbara Basin of pervasive, repeated methane flux in glacial sediments. However, contrary to the "Clathrate Gun Hypothesis" (Kennett et al., 2003), these coring sites are too deep for methane hydrate destabilization to be the cause, implying that a much larger part of the ocean's sedimentary methane may participate in climate or carbon cycle feedback at millennial timescales. We speculate that pulses of methane in these opal-rich sediments could be caused by the sudden release of overpressure in pore fluids that builds up gradually with silica diagenesis. The release could be triggered by seismic shaking on the Aleutian subduction zone caused by hydrostatic pressure increase associated with sea level rise at the start of interstadials.

Bromine and iodine concentrations in sediments and pore fluids of ODP Leg 112

At the Peruvian convergent margin, two distinct pore fluid regimes are recognized from differences in their Cl- concentrations. The slope pore fluids are characterized by low Cl- concentrations, but elevated Br- and I- concentrations due to biogenic production. The shelf pore fluids exhibit elevated Cl- and Br- concentrations due to diffusive mixing with an evaporitic brine. In the slope pore fluids, the Br-, I-, and NH4+ concentrations are elevated following bacterial decomposition of organic matter, but the I- concentrations are in excess of those expected based on mass balance calculations using NH4+ and Br- concentrations. The slope sediment organic matter, which is enriched in iodine from oxidationreduction processes at the oxygenated sediment-water interface, is responsible for this enrichment. The increases in dissolved I- and the I- enrichments relative to NH4+ and Br- correlate well with sedimentation rates because of differential trapping following regeneration. The pore-fluid I-/Br- ratios suggest that membrane ion fiitration is not a major cause of the decreases in Cl- concentrations. Other possible sources for low Cl- water, including meteoric water, clathrate dissociation, and/or mineral dehydration reactions, imply that the diluting component of the slope low-Cl- fluids has flowed at least 1 km through the sediment. The low bottom-water oxygenation in the shelf is responsible for the low (if any) enrichment of iodine in the shelf sediments. Fluctuations in bottom-water oxygen concentrations in the past, however, may be responsible for the observed variations in the sediment I/Br ratios. Comparison of Na+/Cl- and Br-/Cl- molar ratios in the pore fluids shows that the shelf high-Cl- fluid formed from mixing with a brine that formed from seawater concentrated by twelve to nineteen times and probably was modified by halite dissolution. This dense brine, located below the sediment sections drilled, appears to have flowed a distance >500 km through the sediment.

Stable carbon and oxygen isotope record of Neogloboquadrina pachyderma and Cibicides lobatulus from the southwestern Greenland Sea

A core transect across the southwestern Greenland Sea reveals coeval events of extremely negative planktic and benthic delta13C excursions between 40 and 87 ka. The most pronounced event, event 1, began at peak Dansgaard-Oeschger stadial 22 (85 ka) with a duration of 18 k.y. During this episode, incursions of Atlantic Intermediate Water caused a bottom-water warming of up to 8 °C. The amplitude, timing, and geographic pattern of the delta13C events suggest that this bottom-water warming triggered clathrate instability along the East Greenland slope and a methane-induced depletion of delta13CDIC (DIC- dissolved inorganic carbon). Since delta13C event 1 matches a major peak in atmospheric CH4 concentration, this clathrate destabilization may have contributed to the rise in atmospheric CH4 and thus to climate warming over marine isotope stage 5.1.

CO2 concentration and stable isotope ratios of three Antarctic ice cores covering the period from 149.4 - 1.5 kyr before 1950

We present new d13C measurements of atmospheric CO2 covering the last glacial/interglacial cycle, complementing previous records covering Terminations I and II. Most prominent in the new record is a significant depletion in d13C(atm) of 0.5 permil occurring during marine isotope stage (MIS) 4, followed by an enrichment of the same magnitude at the beginning of MIS 3. Such a significant excursion in the record is otherwise only observed at glacial terminations, suggesting that similar processes were at play, such as changing sea surface temperatures, changes in marine biological export in the Southern Ocean (SO) due to variations in aeolian iron fluxes, changes in the Atlantic meridional overturning circulation, upwelling of deep water in the SO, and long-term trends in terrestrial carbon storage. Based on previous modeling studies, we propose constraints on some of these processes during specific time intervals. The decrease in d13C(atm) at the end of MIS 4 starting approximately 64 kyr B.P. was accompanied by increasing [CO2]. This period is also marked by a decrease in aeolian iron flux to the SO, followed by an increase in SO upwelling during Heinrich event 6, indicating that it is likely that a large amount of d13C-depleted carbon was transferred to the deep oceans previously, i.e., at the onset of MIS 4. Apart from the upwelling event at the end of MIS 4 (and potentially smaller events during Heinrich events in MIS 3), upwelling of deep water in the SO remained reduced until the last glacial termination, whereupon a second pulse of isotopically light carbon was released into the atmosphere.

(Table 1) Stable isotope records of clathrite and water

Aragonitic clathrites are methane-derived precipitates that are found at sites of massive near-seafloor gas hydrate (clathrate) accumulations at the summit of southern Hydrate Ridge, Cascadia margin. These platy carbonate precipitates form inside or in proximity to gas hydrate, which in our study site currently coexists with a fluid that is highly enriched in dissolved ions as salts are excluded during gas hydrate formation. The clathrites record the preferential incorporation of 18O into the hydrate structure and hence the enrichment of 16O in the surrounding brine. We measured d18O values as high as 2.27 per mil relative to Peedee belemnite that correspond to a fluid composition of -1.18 per mil relative to standard mean ocean water. The same trend can be observed in Ca isotopes. Ongoing clathrite precipitation causes enrichment of the 44Ca in the fluid and hence in the carbonates. Carbon isotopes confirm a methane source for the carbonates. Our triple stable isotope approach that uses the three main components of carbonates (Ca, C, O) provides insight into multiple parameters influencing the isotopic composition of the pore water and hence the isotopic composition of the clathrites. This approach provides a tool to monitor the geochemical processes during clathrate and clathrite formation, thus recording the evolution of the geochemical environment of gas hydrate systems.

Stable carbon isotope ratios of carbon dioxide from EDC and Berkner Island ice cores for 40-50 ka BP

The stable carbon isotopic signature of carbon dioxide (d13CO2) measured in the air occlusions of polar ice provides important constraints on the carbon cycle in past climates. In order to exploit this information for previous glacial periods, one must use deep, clathrated ice, where the occluded air is preserved not in bubbles but in the form of air hydrates. Therefore, it must be established whether the original atmospheric d13CO2 signature can be reconstructed from clathrated ice. We present a comparative study using coeval bubbly ice from Berkner Island and ice from the bubble-clathrate transformation zone (BCTZ) of EPICA Dome C (EDC). In the EDC samples the gas is partitioned into clathrates and remaining bubbles as shown by erroneously low and scattered CO2 concentration values, presenting a worst-case test for d13CO2 reconstructions. Even so, the reconstructed atmospheric d13CO2 values show only slightly larger scatter. The difference to data from coeval bubbly ice is statistically significant. However, the 0.16 per mil magnitude of the offset is small for practical purposes, especially in light of uncertainty from non-uniform corrections for diffusion related fractionation that could contribute to the discrepancy. Our results are promising for palaeo-atmospheric studies of d13CO2 using a ball mill dry extraction technique below the BCTZ of ice cores, where gas is not subject to fractionation into microfractures and between clathrate and bubble reservoirs.

Methane Clathrates in the Solar System

We review the reservoirs of methane clathrates that may exist in the different bodies of the Solar System. Methane was formed in the interstellar medium prior to having been embedded in the protosolar nebula gas phase. This molecule was subsequently trapped in clathrates that formed from crystalline water ice during the cooling of the disk and incorporated in this form into the building blocks of comets, icy bodies, and giant planets. Methane clathrates may play an important role in the evolution of planetary atmospheres. On Earth, the production of methane in clathrates is essentially biological, and these compounds are mostly found in permafrost regions or in the sediments of continental shelves. On Mars, methane would more likely derive from hydrothermal reactions with olivine-rich material. If they do exist, martian methane clathrates would be stable only at depth in the cryosphere and sporadically release some methane into the atmosphere via mechanisms that remain to be determined. In the case of Titan, most of its methane probably originates from the protosolar nebula, where it would have been trapped in the clathrates agglomerated by the satellite's building blocks. Methane clathrates are still believed to play an important role in the present state of Titan. Their presence is invoked in the satellite's subsurface as a means of replenishing its atmosphere with methane via outgassing episodes. The internal oceans of Enceladus and Europa also provide appropriate thermodynamic conditions that allow formation of methane clathrates. In turn, these clathrates might influence the composition of these liquid reservoirs. Finally, comets and Kuiper Belt Objects might have formed from the agglomeration of clathrates and pure ices in the nebula. The methane observed in comets would then result from the destabilization of clathrate layers in the nuclei concurrent with their approach to perihelion. Thermodynamic equilibrium calculations show that methane-rich clathrate layers may exist on Pluto as well. Key Words: Methane clathrate-Protosolar nebula-Terrestrial planets-Outer Solar System.