Resistivity image beneath an area of active methane seeps in the west Svalbard continental slope

The Arctic continental margin contains large amounts of methane in the form of methane hydrates. The west Svalbard continental slope is an area where active methane seeps have been reported near the landward limit of the hydrate stability zone. The presence of bottom simulating reflectors (BSR) on seismic reflection data in water depths greater than 600 m suggests the presence of free gas beneath gas hydrates in the area. Resistivity obtained from marine controlled source electromagnetic (CSEM) data provides a useful complement to seismic methods for detecting shallow hydrate and gas as they are more resistive than surrounding water saturated sediments. We acquired two CSEM lines in the west Svalbard continental slope, extending from the edge of the continental shelf (250 m water depth) to water depths of around 800 m. High resistivities (5-12 Ωm) observed above the BSR support the presence of gas hydrate in water depths greater than 600 m. High resistivities (3-4 Ωm) at 390-600 m water depth also suggest possible hydrate occurrence within the gas hydrate stability zone (GHSZ) of the continental slope. In addition, high resistivities (4-8 Ωm) landward of the GHSZ are coincident with high-amplitude reflectors and low velocities reported in seismic data that indicate the likely presence of free gas. Pore space saturation estimates using a connectivity equation suggest 20-50% hydrate within the lower slope sediments and less than 12% within the upper slope sediments. A free gas zone beneath the GHSZ (10-20% gas saturation) is connected to the high free gas saturated (10-45%) area at the edge of the continental shelf, where most of the seeps are observed. This evidence supports the presence of lateral free gas migration beneath the GHSZ towards the continental shelf.

New insights into the formation and modification of carbonate-bearing minerals and methane gas in geological systems using multiply substituted isotopologues

This thesis describes the use of multiply-substituted stable isotopologues of carbonate minerals and methane gas to better understand how these environmentally significant minerals and gases form and are modified throughout their geological histories. Stable isotopes have a long tradition in earth science as a tool for providing quantitative constraints on how molecules, in or on the earth, formed in both the present and past. Nearly all studies, until recently, have only measured the bulk concentrations of stable isotopes in a phase or species. However, the abundance of various isotopologues within a phase, for example the concentration of isotopologues with multiple rare isotopes (multiply substituted or 'clumped' isotopologues) also carries potentially useful information. Specifically, the abundances of clumped isotopologues in an equilibrated system are a function of temperature and thus knowledge of their abundances can be used to calculate a sample’s formation temperature. In this thesis, measurements of clumped isotopologues are made on both carbonate-bearing minerals and methane gas in order to better constrain the environmental and geological histories of various samples.

Clumped-isotope-based measurements of ancient carbonate-bearing minerals, including apatites, have opened up paleotemperature reconstructions to a variety of systems and time periods. However, a critical issue when using clumped-isotope based measurements to reconstruct ancient mineral formation temperatures is whether the samples being measured have faithfully recorded their original internal isotopic distributions. These original distributions can be altered, for example, by diffusion of atoms in the mineral lattice or through diagenetic reactions. Understanding these processes quantitatively is critical for the use of clumped isotopes to reconstruct past temperatures, quantify diagenesis, and calculate time-temperature burial histories of carbonate minerals. In order to help orient this part of the thesis, Chapter 2 provides a broad overview and history of clumped-isotope based measurements in carbonate minerals.

In Chapter 3, the effects of elevated temperatures on a sample’s clumped-isotope composition are probed in both natural and experimental apatites (which contain structural carbonate groups) and calcites. A quantitative model is created that is calibrated by the experiments and consistent with the natural samples. The model allows for calculations of the change in a sample’s clumped isotope abundances as a function of any time-temperature history.

In Chapter 4, the effects of diagenesis on the stable isotopic compositions of apatites are explored on samples from a variety of sedimentary phosphorite deposits. Clumped isotope temperatures and bulk isotopic measurements from carbonate and phosphate groups are compared for all samples. These results demonstrate that samples have experienced isotopic exchange of oxygen atoms in both the carbonate and phosphate groups. A kinetic model is developed that allows for the calculation of the amount of diagenesis each sample has experienced and yields insight into the physical and chemical processes of diagenesis.

The thesis then switches gear and turns its attention to clumped isotope measurements of methane. Methane is critical greenhouse gas, energy resource, and microbial metabolic product and substrate. Despite its importance both environmentally and economically, much about methane’s formational mechanisms and the relative sources of methane to various environments remains poorly constrained. In order to add new constraints to our understanding of the formation of methane in nature, I describe the development and application of methane clumped isotope measurements to environmental deposits of methane. To help orient the reader, a brief overview of the formation of methane in both high and low temperature settings is given in Chapter 5.

In Chapter 6, a method for the measurement of methane clumped isotopologues via mass spectrometry is described. This chapter demonstrates that the measurement is precise and accurate. Additionally, the measurement is calibrated experimentally such that measurements of methane clumped isotope abundances can be converted into equivalent formational temperatures. This study represents the first time that methane clumped isotope abundances have been measured at useful precisions.

In Chapter 7, the methane clumped isotope method is applied to natural samples from a variety of settings. These settings include thermogenic gases formed and reservoired in shales, migrated thermogenic gases, biogenic gases, mixed biogenic and thermogenic gas deposits, and experimentally generated gases. In all cases, calculated clumped isotope temperatures make geological sense as formation temperatures or mixtures of high and low temperature gases. Based on these observations, we propose that the clumped isotope temperature of an unmixed gas represents its formation temperature — this was neither an obvious nor expected result and has important implications for how methane forms in nature. Additionally, these results demonstrate that methane-clumped isotope compositions provided valuable additional constraints to studying natural methane deposits.

Oxygen isotopes and hydrocarbon composition of gas hydrate and hydrate water from ODP Leg 165 sites

Ocean Drilling Program (ODP) Leg 164 recovered a number of large solid gas hydrate from Sites 994, 996, and 997 on the Blake Ridge. Sites 994 and 997 samples, either nodular or thick massive pieces, were subjected to laboratory analysis and measurements to determine the structure, molecular and isotopic composition, thermal conductivity, and equilibrium dissociation conditions. X-ray computed tomography (CT) imagery, X-ray diffraction, nuclear magnetic resonance (NMR), and Raman spectroscopy have revealed that the gas hydrates recovered from the Blake Ridge are nearly 100% methane gas hydrate of Structure I, cubic with a lattice constant of a = 11.95 ± 0.05 angström, and a molar ratio of water to gas (hydration number) of 6.2. The d18O of water is 2.67 per mil to 3.51 per mil SMOW, which is 3.5-4.0 heavier than the ambient interstitial waters. The d13C and dD of methane are -66 per mil to -70 per mil and -201 per mil to -206 per mil, respectively, suggesting that the methane was generated through bacterial CO2 reduction. Thermal conductivity values of the Blake Ridge hydrates range from 0.3 to 0.5 W/(m K). Equilibrium dissociation experiments indicate that the three-phase equilibrium for the specimen is 3.27 MPa at 274.7 K. This is almost identical to that of synthetic pure methane hydrate in freshwater.

Geochemistry of hydrate gas and water at DSDP Hole 84-570

Molecular and isotopic measurements of gas and water obtained from a gas hydrate at Site 570, DSDP Leg 84, are reported. The hydrate appeared to be Structure I and was composed of a solid framework of water molecules enclosing methane and small amounts of ethane and carbon dioxide. Carbon isotopic values for the hydrate-bound methane, ethane, and carbon dioxide were -41 to about -44, -27, and -2.9 per mil, respectively. The d13C-C1 values are consistent with void gas values that were determined to have a biogenic source. A significant thermogenic source was discounted because of high C1/C2 ratios and because the d13C-CO2 values in these sections were also anomalously heavy (or more positive) isotopically, suggesting that the methane was formed biogenically by reduction of heavy CO2 . The isotopically heavy hydrate d13C-C2 is also similar to void gas isotopic compositions and is either a result of low-temperature diagenesis producing heavy C2 in these immature sediment sections or upward migration of deeper thermogenic gas. The salinity of the hydrate water was 2.6 per mil with dDH2O and d18OH2O values of +1 and +2.2 per mil, respectively.

Modelling carbon membranes for gas and isotope separation

Molecular modelling has become a useful and widely applied tool to investigate separation and diffusion behavior of gas molecules through nano-porous low dimensional carbon materials, including quasi-1D carbon nanotubes and 2D graphene-like carbon allotropes. These simulations provide detailed, molecular level information about the carbon framework structure as well as dynamics and mechanistic insights, i.e. size sieving, quantum sieving, and chemical affinity sieving. In this perspective, we revisit recent advances in this field and summarize separation mechanisms for multicomponent systems from kinetic and equilibrium molecular simulations, elucidating also anomalous diffusion effects induced by the confining pore structure and outlining perspectives for future directions in this field.

Catalyst size effects on the growth of single-walled nanotubes in neutral and plasma systems

The results of large-scale (∼109 atoms) numerical simulations of the growth of different-diameter vertically-aligned single-walled carbon nanotubes in plasma systems with different sheath widths and in neutral gases with the same operating parameters are reported. It is shown that the nanotube lengths and growth rates can be effectively controlled by varying the process conditions. The SWCNT growth rates in the plasma can be up to two orders of magnitude higher than in the equivalent neutral gas systems. Under specific process conditions, thin SWCNTs can grow much faster than their thicker counterparts despite the higher energies required for catalyst activation and nanotube nucleation. This selective growth of thin SWCNTs opens new avenues for the solution of the currently intractable problem of simultaneous control of the nanotube chirality and length during the growth stage.

Robust design optimization of gas turbine compression systems

Gas turbine compression systems are required to perform adequately over a range of operating conditions. Complexity has encouraged the conventional design process for compressors to focus initially on one operating point, usually the most commonor arduous, to draw up an outline design. Generally, only as this initial design is refined is its offdesign performance assessed in detail. Not only does this necessarily introduce a potentially costly and timeconsuming extra loop in the design process, but it also may result in a design whose offdesign behavior is suboptimal. Aversion of nonintrusive polynomial chaos was previously developed in which a set of orthonormal polynomials was generated to facilitate a rapid analysis of robustness in the presence of generic uncertainties with good accuracy. In this paper, this analysis method is incorporated in real time into the design process for the compression system of a three-shaft gas turbine aeroengine. This approach to robust optimization is shown to lead to designs that exhibit consistently improved system performance with reduced sensitivity to offdesign operation.