SEARCH · Engineering Papers
Results for “Geochemistry & Geophysics”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Experimental Observation of a New Attenuation Mechanism in hcp ‐Metals That May Operate in the Earth's Inner Core
Abstract Seismic observations show the Earth's inner core has significant and unexplained variation in seismic attenuation with position, depth and direction. Interpreting these observations is difficult without knowledge of the visco‐ or anelastic dissipation processes active in iron under inner core conditions. Here, a previously unconsidered attenuation mechanism is observed in zinc, a low pressure analog ofhcp‐iron, during small strain sinusoidal deformation experiments. The experiments were performed in a deformation‐DIA combined with X‐radiography, at seismic frequencies (∼0.003–0.1 Hz), high pressure and temperatures up to ∼80% of melting temperature. Significant dissipation (0.077 ≤ Q −1 (ω) ≤ 0.488) is observed along with frequency dependent softening of zinc's Young's modulus and an extremely small activation energy for creep (⩽7 kJ mol −1 ). In addition, during sinusoidal deformation the original microstructure is replaced by one with a reduced dislocation density and small, uniform, grain size. This combination of behavior collectively reflects a mode of deformation called “internal stress superplasticity”; this deformation mechanism is unique to anisotropic materials and activated by cyclic loading generating large internal stresses. Here we observe a new form of internal stress superplasticity, which we name as “elastic strain mismatch superplasticity.” In it the large stresses are caused by the compressional anisotropy. If this mechanism is also active inhcp‐iron and the Earth's inner‐core it will be a contributor to inner‐core observed seismic attenuation and constrain the maximum inner‐core grain‐size to ≲10 km.
High Sodium Solubility in Magnesiowüstite in Iron‐Rich Deep Lower Mantle
Abstract (Mg,Fe)O ferropericlase‐magnesiowüstite has been proposed to host the majority of Earth's sodium, but the mechanism and capacity for incorporating the alkali cation remain unclear. In this work, experiments in the laser‐heated diamond anvil cell and first‐principles calculations determine the solubility of sodium and favorability of sodium incorporation in iron‐rich magnesiowüstite relative to (Mg,Fe)SiO 3 bridgmanite. Reaction of Mg/(Mg + Fe) (Mg#) 55 and 28 olivine with NaCl at 33–128 GPa and 1600–3000 K produces iron‐rich magnesiowüstite containing several percent sodium, while iron‐rich bridgmanite contains little to no detectable sodium. In sodium‐saturated magnesiowüstite, sodium number [Na/(Na + Mg + Fe)] is 2–5 atomic percent at pressures below 60 GPa and drastically increases to 10–20 atomic percent at deep lower mantle pressures. For these two compositions, there is no significant dependence of the results on Mg#. Our calculations not only show consistent results with experiments but further indicate that such an increase in solubility and partitioning of Na into magnesiowüstite is driven by the spin transition in iron. These results provide fundamental constraints on the crystal chemistry of sodium at lower‐mantle conditions. If the sodium capacity of (Mg,Fe)O is not strongly dependent on Mg#, (Mg,Fe)O in the lower mantle may have the capacity to store the entire sodium budget of the Earth.
Multi‐Scale Seismic Imaging of the Ridgecrest, CA, Region With Waveform Inversion of Regional and Dense Array Data
Abstract We develop an inversion procedure for deriving multi‐scale velocity models with waveform inversions of earthquake and ambient noise data at multi‐frequency bands recorded by regional and dense sensor configurations. The method is applied for the area around the 2019 Ridgecrest earthquake rupture zones, utilizing data recorded by regional stations and dense 2D and 1D arrays with station spacings of ∼5 km and ∼100 m, respectively. Starting with regional Vp, Vs models and locations of Ridgecrest aftershocks, the velocity models and event locations are improved iteratively by inversions of waveforms recorded by regional stations and the 2D array, using a minimum spectral element size of ∼600 m. Waveforms from local events recorded by dense 1D arrays across the M7.1 rupture zone with frequencies of up to 10 Hz are used to resolve small‐scale features of the rupture zone and shallow crust with a local spectral element size of 80 m. The refined models provide self‐consistent descriptions of the rupture zone and the shallow crust embedded in the regional structures. The results reveal pronounced low Vs and high Vp/Vs in the M6.4 and M7.1 rupture zones coinciding with concentrations of seismicity, and also around the Garlock fault and in several local basins. We also observe clear velocity contrasts across the Garlock fault with polarity reversals along strike and with depth. The obtained multi‐scale velocity models can be used to improve derivations of earthquake source properties, simulations of dynamic ruptures and ground motions, and the understanding of fault and tectonic processes in the region.
Formation of a Large Cold Groundwater Mantle Helium Anomaly and High Temperature Geothermal Resources in Response to Bimodal Magmatism Near Roosevelt Hot Springs and Utah FORGE, Milford Valley, Southwest Utah
Abstract A large mantle helium anomaly and separate domains of high heat flow are the predominant manifestations of bimodal magmatic activity in the Milford valley. The mantle helium anomaly (1.9–2.6 R/Ra) covers 270 km 2 and is subdivided into two separated domains: a cold shallow groundwater regime and high temperature hydrothermal activity. The zone of anomalous heat flow covers >100 km 2 and is also subdivided into two adjacent domains, comprising hydrothermal activity at Roosevelt Hot Springs (RHS) (3–7 W/m 2 ) and conductive heat flow (100–180 mW/m 2 ). While the transfer of heat and mantle helium at RHS are coupled, heat and helium transfer are decoupled in the adjacent cold groundwater regime to the west. Both the mantle helium and geothermal anomalies are attributed to recent mafic‐felsic magmatic intrusions of >400 km 3 , however, the absence of volcanic eruptions <500,000 years indicates magmas stall before rising to shallow crustal level <10 km depth. Deep level magmatism produces a felsic composition melt, which is inferred to be responsible for the widespread and near uniform range of diluted mantle helium values. A thick and impermeable mass of crystalline granitic basement rock at the mid‐crustal level divides the ascent of mantle helium into separate flow paths. It may also impede the rise of buoyant magma trapping thermal energy that facilitates partial melting, slow cooling, and development of a thick thermal aureole. Partitioning of convective and conductive thermal regimes and independent flow paths supplying deeply derived helium characterize the development of a large long‐lived magma‐related geothermal system.
Geodetic Evidence for Distributed Shear Below the Brittle Crust of the Walker Lane, Western United States
Abstract Models of active deformation of the Earth's crust are predominantly represented with dislocations having a downdip continuation into the lower crust, where the fault slips continuously. This model predicts surface strain accumulation concentrated near the fault during the interseismic period. In an alternative model, faults do not extend beneath the elastic portion of the crust and are accompanied by a wide zone of distributed shear underneath, predicting a more constant strain rate lacking concentrations at the faults. We use high‐precision GPS data collected across the northern and central Walker Lane, USA— a region of complex faulting near the western edge of the Basin and Range Province to evaluate which model is appropriate. Despite the existence of dense continuous and semi‐continuous geodetic networks that have been surveyed for ∼20 years, the horizontal velocities reveal no evidence of localized strain accumulation across the fault surface expressions. Instead, deformation within the Walker Lane is uniformly linear, suggesting that the surface deformation reflects distributed shear within the ductile crust rather than focused deformation at faults. This suggests no downdip extension of the faults below the seismogenic layer. The shear zone is 172 ± 6 km wide in the northernmost Walker Lane narrowing to 116 ± 4 km in the central Walker Lane. The total velocity budget across the shear zone is 7.2 ± 0.1 mm/yr in the north, increasing to 10.1 ± 0.1 mm/yr in the central Walker Lane. We conclude that assuming the presence of lower crustal dislocations when estimating geodetic faults slip rates may be inappropriate.
Stress Field Dynamics and Fault Slip Potential in the Paradox Basin
Abstract The Paradox Basin, straddling Utah, Colorado, Arizona, and New Mexico is characterized by an intricate amalgamation of evaporites and clastic layers and is dominated by prominent salt walls and related subsurface structures. Our research offers a new examination of the stress distribution across the basin, deriving from continuous and discrete stress measurements conducted in boreholes in the region and focal mechanism analysis, emphasizing variations over salt structures. Integrating Coulomb failure criteria with probabilistic methods, we assess potential fault movements resulting from fluid pressure alterations. Our approach provides a comprehensive understanding of the Paradox Basin's state of stress, showing a continuous change of the maximum horizontal stress orientation from N‐S at the Wasatch Fault Zone to WNW‐ESE in the northern part of the Paradox Basin and to WSW‐ENE in the southern part of the basin. Further East, into the Colorado Plateau and the Uncompahgre Uplift, the S H max orientation becomes E‐W. Decoding stress orientation dynamics has enabled critical insights into fault slip potential, especially in the basin's northern region. The salt wall faults are less likely to slip, and the Paradox Formation's evaporite and clastic rock sequence can serve as a potential low seismic risk target for carbon storage and hydrocarbon extraction.
Velocity Dependence of Rate‐And‐State Friction in Granular Fault Gouge and Implications for Slow‐Slip Events
Abstract The standard rate‐and‐state friction (RSF) has extensively captured frictional behaviors, but it fails to explain the velocity dependence of frictional stability transition and widespread slow‐slip events (SSEs) in experiments and nature adequately. An alternative microphysical Chen‐Niemeijer‐Spiers (CNS) model can well describe the velocity dependence of frictional behaviors of granular gouges. Using the original CNS model, standard RSF parameters can be quantified microphysically. However, some micro‐parameters are not easy to estimate quantitatively, making it difficult to extrapolate to natural and experimental conditions. Here, we simplify the microphysically‐derived RSF parameters including direct effect a , evolution effect b , and critical slip distance D c , as well as equivalent values ( a eq , b eq , and D eq ). The simplified friction parameters directly illustrate their velocity dependence, namely the essentially constant a , a eq , and D c , negatively velocity‐dependent b and b eq , as well as varying D eq for different laws. They are roughly consistent with experimental results in various fault gouges. A modified CNS model is further derived from the original CNS model, establishing a direct link between the standard RSF and CNS models. The modified CNS model exhibits virtually identical frictional behaviors to the original CNS, but differs from the standard RSF at large velocity perturbations. Moreover, the linearized stability analysis indicates that the critical stiffness for the modified CNS model is velocity‐dependent. Compared with the standard RSF, the modified CNS model not only explains the velocity dependence of frictional stability transition, but also exhibits a more gradual transition for SSEs with a broader range of stiffness ratios.
Seismic and acoustic signals from the 2014 ‘interstellar meteor’
SUMMARY We conduct a thorough analysis of seismic and acoustic data purported to be from the so-called ‘interstellar meteor’ which entered the Earth’s atmosphere off the coast of Papua New Guinea on 2014 January 08. Previous work had suggested that this meteor may have been caused by an alien spacecraft burning up in the atmosphere. We conclude that both previously reported seismic signals are spurious—one has characteristics suggesting a local vehicular-traffic-based origin; whilst the other is statistically indistinguishable from the background noise. As such, previously reported localizations based on this data are unreliable. Analysis of acoustic data provides a best-fitting location estimate which is very far ($\sim$170 km) from the reported fireball location. Accordingly, we conclude that material recovered from the seafloor and purported to be from this event is almost certainly unrelated to it, and is likely of more mundane (non-interstellar) origin.
Graph-learning approach to combine multiresolution seismic velocity models
SUMMARY The resolution of velocity models obtained by tomography varies due to multiple factors and variables, such as the inversion approach, ray coverage, data quality, etc. Combining velocity models with different resolutions can enable more accurate ground motion simulations. Toward this goal, we present a novel methodology to fuse multiresolution seismic velocity maps with probabilistic graphical models (PGMs). The PGMs provide segmentation results, corresponding to various velocity intervals, in seismic velocity models with different resolutions. Further, by considering physical information (such as ray path density), we introduce physics-informed probabilistic graphical models (PIPGMs). These models provide data-driven relations between subdomains with low (LR) and high (HR) resolutions. Transferring (segmented) distribution information from the HR regions enhances the details in the LR regions by solving a maximum likelihood problem with prior knowledge from HR models. When updating areas bordering HR and LR regions, a patch-scanning policy is adopted to consider local patterns and avoid sharp boundaries. To evaluate the efficacy of the proposed PGM fusion method, we tested the fusion approach on both a synthetic checkerboard model and a fault zone structure imaged from the 2019 Ridgecrest, CA, earthquake sequence. The Ridgecrest fault zone image consists of a shallow (top 1 km) high-resolution shear-wave velocity model obtained from ambient noise tomography, which is embedded into the coarser Statewide California Earthquake Center Community Velocity Model version S4.26-M01. The model efficacy is underscored by the deviation between observed and calculated traveltimes along the boundaries between HR and LR regions, 38 per cent less than obtained by conventional Gaussian interpolation. The proposed PGM fusion method can merge any gridded multiresolution velocity model, a valuable tool for computational seismology and ground motion estimation.
SIMS Data Correction Procedure for Quasi‐Simultaneous Arrival ( QSA ) Under‐counting and Ramifications of Misapplication
Isotope geochemistry requires isotope ratios measured using secondary ion mass spectrometry (SIMS) to be made with optimal precision and accuracy. Under some analytical conditions when using electron multiplier detectors, secondary ions may be under‐counted because of quasi‐simultaneous arrival (QSA) at the first dynode. The relative magnitude of the associated QSA correction to raw measured isotopic ratios can be up to seventy permil or more. Therefore, not applying the correction, or misapplication of it could lead to significant inaccuracies in published isotope ratio data. Examples and ramifications of the latter are described in addition to a straightforward procedure for QSA under‐counting correction.
Compositional evidence for chondrule origins of low‐Ca pyroxenes in comet Wild 2 and a giant cluster IDP
Abstract A literature compilation of 1136 low‐Ca pyroxene compositions from chondrules from 12 primitive type 2–3 carbonaceous, ordinary and enstatite chondrite groups define unique regions on an Al 2 O 3 and Cr 2 O 3 diagram when compared to low‐Ca pyroxenes from equilibrated type 4‐6 chondrites. Measured compositions of 100 low‐Ca pyroxenes from comet Wild 2 and a giant cluster IDP of probable cometary origin are similar to each other and fall in the type 2–3 chondrite chondrule region suggesting that most of the pyroxenes likely formed in the solar nebula like conventional chondrules. The data imply that most low Ca‐pyroxenes from comet Wild 2 and the giant cluster IDP formed from igneous crystallization processes and did not experience significant thermal metamorphism, indicating that the low‐Ca pyroxenes were unlikely incorporated into large parent bodies prior to accretion in their respective comet bodies. An intriguing group of nine low‐Ca pyroxenes from comet Wild 2 with low Cr and Al that fall where type 4–6 chondrites are located are interpreted as products of condensation. The compositional data combined with previously measured oxygen isotopes on 17 low‐Ca pyroxenes support earlier conclusions that comet samples have links with carbonaceous, ordinary, and possibly enstatite chondrite groups. Our results provide additional evidence that comets accreted materials from multiple chondrule reservoirs throughout the solar nebula.
Advancing attenuation estimation through integration of the Hessian in multiparameter viscoacoustic full-waveform inversion
Accurate seismic attenuation models of subsurface structures not only enhance subsequent migration processes by improving fidelity, resolution, and facilitating amplitude-compliant angle gather generation but also provide valuable constraints on subsurface physical properties. Leveraging full-wavefield information, multiparameter viscoacoustic full-waveform inversion ( Q-FWI) simultaneously estimates seismic velocity and attenuation ( Q) models. However, a major challenge in Q-FWI is the contamination of crosstalk artifacts, where inaccuracies in the velocity model are mistakenly mapped to the inverted attenuation model. While incorporating the Hessian is expected to mitigate these artifacts, the explicit implementation is prohibitively expensive due to its formidable computational cost. In this study, we formulate and develop a Q-FWI algorithm via the Newton-conjugate gradient (CG) framework, where the search direction at each iteration is determined through an internal CG loop. In particular, the Hessian is integrated into each CG step in a matrix-free fashion using the second-order adjoint-state method. We find through synthetic experiments that our Newton-CG Q-FWI significantly mitigates crosstalk artifacts compared with the limited-memory Broyden-Fletcher-Goldfarb-Shanno method and the CG method, albeit with a notable computational cost. In the discussion of several key implementation details, we also determine the significance of the approximate Gauss-Newton Hessian, the second-order adjoint-state method, and the two-stage inversion strategy.
Temporal and spatial characterization of a thermogenic, fault-controlled gas hydrate system, Woolsey Mound, Gulf of Mexico
Woolsey Mound, located at Mississippi Canyon Lease Block 118 (MC118), is the site of the Gulf of Mexico hydrate research consortium’s seafloor observatory, where gas hydrates outcrop at the seafloor. The presence of gas hydrates in the mound is confirmed directly by coring and indirectly by 3D seismic reflection data. Craters, pockmarks, chemosynthetic communities, and authigenic carbonates populate the seafloor at Woolsey Mound. Each crater is characterized by a network of shallow crestal faults that connect the hydrate mound to the underlying allochthonous salt body. We characterize the temporal and spatial evolution of gas hydrates at Woolsey Mound under natural perturbations using four collocated 3D seismic reflection data sets that span over 14 years. Data acquisition differences embedded in the data sets arising from variation in geometry, sample rate, and phase are minimized using the “cross-equalization” method. Our results indicate that hydrate formation and dissociation vary temporally and spatially in close connection to the shallow crestal faults. Evidence of gas hydrate dissociation is observed over a period of three years (2000–2003), where major dissociation occurred along the southern portion of the crestal fault in the southeast crater. The dissociation is less prominent in the southwest crater. Evidence of methane venting is observed between 2000 and 2010, which is mostly concentrated in the southeast crater. The residual amplitude anomalies observed between 2000 and 2014 in the mound are mostly positive, implying that the methane venting had increased significantly. The positive anomalies are correlated with the methane seepage recorded in 2011. Our results indicate the evolution of a fault-controlled gas hydrate system in the northern Gulf of Mexico, which would aid in assessing its impact on the seafloor.
Seismic Monitoring near Ithaca, New York, Reveals Nonuniform Distribution of Microseismicity in an Intraplate Region
Abstract Cornell University intends to use a deep direct-use geothermal system to heat its Ithaca, New York, campus. In preparation for this project, the Cornell Seismic Network has been monitoring the background seismicity in this intraplate region since 2019. From January 2020 to June 2023, 95 events were detected within 20 km of the proposed geothermal well site, with local magnitudes ranging from −1.02 to 0.56. None of these events appear in regional or national catalogs. Events locate in a narrow geographic band, with one-fourth exhibiting multimodal hypocentral probability peaks both near the surface and at 1–4 km depth. We relocate events with a joint hypocenter and 1D velocity model inversion, in addition to a fully nonlinear method, and then compare observations with synthetic waveforms. Together, these approaches provide strong evidence for >95% of events locating at the surface or within the 3-km-thick sedimentary sequence. We explore how anthropogenic activity and regional topographic stress may contribute to frequent surficial events. This information is critical for characterizing the background microseismicity for comparison during future geothermal operations. Ithaca’s geology of Paleozoic sediments overlying Precambrian crystalline basement is typical of many continental interiors, so these results also provide insight into intraplate microseismicity patterns.
Sulfur speciation in dacitic melts using X-ray absorption near-edge structure spectroscopy of the S K -edge (S-XANES): Consideration of radiation-induced changes and the implications for sulfur in natural arc systems
Abstract The synchrotron technique of micro X-ray absorption near-edge structure spectroscopy at the sulfur K-edge (S-XANES) provides a unique opportunity to measure the proportion of different oxidation states of sulfur (S) in silicate glasses. Although applied extensively in the analysis of basaltic silicate glasses, few S-XANES studies have investigated variations in S oxidation states with fO2 in felsic silicate glasses. In addition, no study has systematically compared the S-XANES results obtained from the same samples at different photon flux densities to quantify the relationship between exposure time and changes in S speciation in silicate glass, as has been done for Fe and V. This study evaluates observed differences in S speciation measured in experimentally produced H2O-saturated dacitic glasses over a range of reducing to oxidizing conditions (from log fO2 = ΔFMQ-0.7 to ΔFMQ+3.3; FMQ is the fayalite-magnetite-quartz mineral redox buffer) and equilibrated at 1000 °C and 300 MPa. S-XANES spectra were collected at three different photon flux densities using three microspectroscopy beamlines. As is observed in S-XANES analyses of basaltic silicate glasses, beam-induced changes to the S6+/ΣS are observed as a function of photon flux density and beam exposure time. Our results demonstrate that silicate glasses of dacitic composition undergo beam-induced photo-reduction in samples equilibrated at ΔFMQ > +1.75 and photo-oxidation if equilibrated at ΔFMQ < +1. The time required to observe beam-induced changes in the spectra varies as a function of flux density, and our study establishes an upper photon density limit at ~1.0 × 1012 photons/µm2. The S6+/ΣS calculated from spectra collected below this absorbed photon limit at intermediate flux densities (~1–4 × 109 photons/s per µm2) are affected by beam damage, as no conditions were found to be completely free of beam-induced changes. However, the S6+/ΣS ratios calculated below the limit at intermediate flux densities are consistent with thermodynamic constraints, demonstrating that S6+/ΣS ratios calculated from S-XANES spectra can be considered reliable for estimating the oxygen fugacity. Our results carry important implications for the S budget of felsic magmas and dissolution mechanisms in evolved melts. While our results from all three flux densities show the presence of S4+ dissolved in relatively oxidized (ΔFMQ > +1.75) dacitic glass, even in the spectra exposed to the lowest photon densities, we are unable to rule out the possibility that the S4+ signal is the result of instantaneous X-ray irradiation induced beam damage using S-XANES alone. When our spectra are compared to S-XANES spectra from basaltic silicate glasses, important differences exist in the solubility of S2– and S6+ between dacitic silicate melts, pointing to differences in solubility mechanisms as melt composition changes. This study highlights the need for further investigation into beam damage systematics, presence of S4+, and the solubility mechanisms of different oxidation states of S as silicate melt composition changes.
Bobfinchite, Na[(UO2)8O3(OH)11]·10H2O, a new Na-bearing member of the schoepite family
Abstract The new mineral bobfinchite (IMA2020-082), Na[(UO2)8O3(OH)11]·10H2O, was found in the Burro mine, Slick Rock district, San Miguel County, Colorado, U.S.A., where it occurs as an oxidation product of uraninite on asphaltite matrix in intimate association with gypsum, natrozippeite, metaschoepite, and uranopilite. Bobfinchite crystals are transparent to translucent, yellow, lozenge-shaped disks up to 0.3 mm wide. Crystals are flattened on [100] and exhibit the forms {100}, {011}, {021}, {021}, and {011}. Bobfinchite has a pale-yellow streak and emits very dim yellow fluorescence under 365 nm ultraviolet illumination. The crystals are brittle with very good {100} cleavage and irregular, stepped fracture. The Mohs hardness is ca. 2 based on scratch tests. The calculated density is 5.044 g/cm3 based on the empirical formula and 5.036 g/cm3 for the ideal formula. Bobfinchite is optically biaxial (–), with α = 1.690(5), β = 1.7205(5), and γ = 1.730(5) (white light). The measured 2V, estimated from the interference figure, is 55(5)° and the calculated value is 59.1°. Dispersion is moderate, r > v; orientation: X = a, Y = b, Z = c; pleochroism: X nearly colorless, Y yellow, Z yellow; X < Y ≈Z. Electron microprobe analysis provided the empirical formula (Na0.99Pb0.02)[(UO2)7.99O3(OH)11]·10H2O. The five strongest X-ray powder diffraction lines are [dobs in Å(I)(hkl)]: 7.34(100)(200), 3.59(50)(024), 3.23(60)(224), 3.18(36)(240), and 2.01(23)(624,551,208,640,346). Bobfinchite is orthorhombic, Pbcn, a = 14.6249(9), b = 14.0389(10), c = 16.6923(10) Å, V = 3427.2(4) Å3, and Z = 4. The structure of bobfinchite (R1 = 0.0330 for 3770 I > 4σI) is built from uranyl oxide-hydroxide sheets that adopt the fourmarierite topology, with interlayer Na+ and H2O groups. Both the sheet and interlayer topology mimic those observed in natural and synthetic Na-metaschoepites studied previously, and as seen in other uranyl oxide hydrate minerals, charge balance is achieved at specific sites in the sheet through the substitution O2– ↔ (OH)–.
UT-GOM2-2 Preliminary Report: Terrebonne Basin Northern Gulf of Mexico, 30 July-28 September 2023
In the summer and fall of 2023, the Gulf of Mexico Deepwater Hydrate Coring Expedition (UT-GOM2-2) drilled, cored, made downhole measurements, and analyzed samples from the seafloor to the base of the gas hydrate stability zone in one location (Site H, WR313) in the Terrebonne basin, deepwater Gulf of Mexico. Analyses of data and samples from the expedition will inform biological, geochemical, and geomechanical models to constrain the role of gas hydrates in the carbon cycle and the potential for gas hydrates as an energy resource. Pressure and conventional cores were collected continuously to a depth of 155.1 meters below the seafloor (mbsf). At deeper depths, cores were taken periodically from hydrate-bearing sands and their bounding muds to a total depth of 861.3 mbsf. 162.6 m of conventional core and 54.8 m of pressure core were obtained. Twelve temperature measurements were made between 27.1 and 144.5 mbsf to determine the geothermal gradient. At the seafloor, more than 4 m of sandy silt of unknown origin was encountered. Beneath this sand, to a depth of ~200 mbsf, the section was composed of interbedded mud and biogenic carbonate ooze. The biogenic ooze correlated to low density and high porosity intervals observed in the previously acquired logging while drilling (LWD) data and as measured. Calcareous nannofossil biostratigraphy constrains the entire record to the Pleistocene (< 0.91 million years) with a pronounced increase in sedimentation rate with depth. Beneath 200 mbsf, the section was predominantly composed of mud with two thicker, hydrate-bearing coarse-grained intervals, which are commonly known as the Blue and Orange sands. The dissolved gas concentration was quantified from pressure cores. In the shallow section, dissolved methane concentration increased below the sulfate-methane transition zone (SMTZ) and reaches saturation (the limit of solubility for methane) at 147 mbsf. Gas expansion was very common in conventional and depressurized pressure (conventionalized) cores below the SMTZ. At deeper depths, the methane concentration within muds bounding the Blue and Orange reservoirs was generally found to be less than saturation. The dissolved and hydrate gas composition is consistent with a microbial source, containing greater than 99.99% methane and only trace concentrations of ethane, propane, and butane. The methane to ethane ratio (C 1 /C 2 ) and the methane to ethane plus propane (C 1 /(C 2 +C 3 )) decrease with depth down to at least 678 mbsf, mainly driven by the increase in ethane with depth. It is unclear if this trend continues through the Orange sand interval. The δ 13 C isotopic signature of methane ranges between -69.9 and -78.5 ‰ Vienna Pee Dee Belemnite (VPDB). Pressure core recovery of all sandy intervals was poor. However, pressure core logs of the Orange sand show intervals of low density and high velocity, which are indicative of high hydrate saturation. One core from within the Orange sand was composed of interbedded graded sandy silt and mud. The sandy silts from this core are composed of mainly quartz and feldspar with some lithics. Most of the recovered pressure core samples are maintained at near in-situ pressure and temperature (within the hydrate stability field) at the University of Texas Pressure Core Center awaiting analysis. In the shallow section, samples will be used to determine the flux of organic carbon through the basin system, find the rate at which that carbon was consumed, and understand the microbial population responsible for these processes. In the deeper section, samples from in and around the hydrate reservoirs will be used to determine the petrophysical properties of the reservoir and bounding seals in these systems.