Scholarly article on topic 'Deep Sediment-Sourced Methane Contribution to Shallow Sediment Organic Carbon: Atwater Valley, Texas-Louisiana Shelf, Gulf of Mexico'

Deep Sediment-Sourced Methane Contribution to Shallow Sediment Organic Carbon: Atwater Valley, Texas-Louisiana Shelf, Gulf of Mexico Academic research paper on "Environmental engineering"

CC BY
0
0
Share paper
Academic journal
Energies
OECD Field of science
Keywords
{""}

Academic research paper on topic "Deep Sediment-Sourced Methane Contribution to Shallow Sediment Organic Carbon: Atwater Valley, Texas-Louisiana Shelf, Gulf of Mexico"

Energies 2015, 8, 1561-1583; doi:10.3390/en8031561

OPEN ACCESS

energies

ISSN 1996-1073

www.mdpi.com/journal/energies

Article

Deep Sediment-Sourced Methane Contribution to Shallow Sediment Organic Carbon: Atwater Valley, Texas-Louisiana Shelf, Gulf of Mexico

Richard B. Coffin 1,t,ft, Christopher L. Osburn 2, Rebecca E. Plummer 3, Joseph P. Smith 4, Paula S. Rose 1,t and Kenneth S. Grabowski 1

1 Naval Research Laboratory, Washington, DC 20735, USA; E-Mails: paula.rose@tamucc.edu (P.S.R.); kenneth.grabowski@nrl.navy.mil (K.S.G.)

2 Marine, Earth, and Atmospheric Sciences, North Carolina State University, Raleigh, NC 27695, USA; E-Mail: closburn@ncsu.edu

3 Department of Geology, University of Maryland, College Park, MD 20742, USA; E-Mail: rplummer@umd.edu

4 Oceanography Department, US Naval Academy, Annapolis, MD 21402, USA; E-Mail: jpsmith@usna.edu

f Current Address: Department of Physical and Environmental Sciences, Texas A&M University—Corpus Christi, Corpus Christi, TX 78412, USA.

* Author to whom correspondence should be addressed; E-Mail: richard.coffin@tamucc.edu; Tel.: +1-361-825-2456.

Academic Editor: Enrico Sciubba

Received: 23 August 2014 /Accepted: 11 February 2015 /Published: 18 February 2015

Abstract: Coastal methane hydrate deposits are globally abundant. There is a need to understand the deep sediment sourced methane energy contribution to shallow sediment carbon relative to terrestrial sources and phytoplankton. Shallow sediment and porewater samples were collected from Atwater Valley, Texas-Louisiana Shelf, Gulf of Mexico near a seafloor mound feature identified in geophysical surveys as an elevated bottom seismic reflection. Geochemical data revealed off-mound methane diffusion and active fluid advection on-mound. Gas composition (average methane/ethane ratio ~11,000) and isotope ratios of methane on the mound (average S13CcH4(g) = -71.2%o; A14CcH4(g) = -961%o) indicate a deep sediment, microbial source. Depleted sediment organic carbon values on mound (S13Csoc = -25.8%0; A14Csoc = -930%o) relative to off-mound (513Csoc = -22.5%0;

A14Csoc = -629%o) suggest deep sourced ancient carbon is incorporated into shallow sediment organic matter. Porewater and sediment data indicate inorganic carbon fixed during anaerobic oxidation of methane is a dominant contributor to on-mound shallow sediment organic carbon cycling. A simple stable carbon isotope mass balance suggests carbon fixation of dissolved inorganic carbon (DIC) associated with anaerobic oxidation of hydrate-sourced CH4 contributes up to 85% of shallow sediment organic carbon.

Keywords: methane; advection; geochemistry; carbon isotopes; sediment carbon

1. Introduction

Sediment organic carbon (SOC) composition and provenance have been extensively studied in the Gulf of Mexico (GoM). Generally, inputs of terrestrially-derived organic carbon (OC) dominate near shore sediments and decrease with distance offshore [1-5]. Terrestrial-derived OC in surficial sediments of the GoM shifts from C3 plant material near shore to highly degraded, soil-derived material offshore [1-3]. Surface sediments on the slope (365-2270 m water depth) have A14C values between -309.1% and -228.6% corresponding to 64% to 78% modern carbon [3]. Mayer et al. [4] calculated that surficial SOC on the slope is 41% to 46% of marine origin. Carbon isotope signatures typical for marine phytoplankton have been measured in surface sediments at water column depths from 74 to 2250 m (mean 513C value = -20.8%; [6]). In the same study, depleted (513C = -24.1%) sedimentary organic matter was observed at a cold seep site (688 m water depth) that suggested a contribution of seep hydrocarbons to sediment organic matter.

Methane hydrate deposits are abundant in deep sediments along the Texas-Louisiana Shelf [7,8]. Hydrate-sourced CH4 may contribute to the sediment OC pool and should be considered in investigations of carbon cycling and regional carbon mass balance calculations. During anaerobic oxidation of methane (AOM), bacteria utilize seawater-sourced sulfate (SO42-) as a terminal electron acceptor to oxidize methane as an energy source, and produce sulfide (H2S) and dissolved inorganic carbon (DIC) in the process [9,10]. During this oxidation porewater DIC sourced from downward seawater diffusion, upward deep system advection or diffusion, and AOM is incorporated into microbial biomass [11-14]. This is a contrast to aerobic CH4 assimilation, where methane-derived carbon is directly incorporated into microbial biomass [15]. Deep sediment hydrocarbons, including CH4, have been shown to contribute to carbon cycling in shallow sediments and the water column [16-21]. However, the indirect contribution of hydrate-sourced CH4 to shallow sediment carbon pools through AOM and subsequent carbon dioxide (CO2) fixation has not been thoroughly evaluated.

Shallow sediment CH4 in Atwater Valley, Gulf of Mexico has been shown to be dominated by biogenic gas from deep sediments [17]. In this study, two sites in Atwater Valley were contrasted: A sediment mound with active fluid advection and high vertical CH4 flux and an off-mound site exhibiting steady-state CH4 diffusion [22,23]. Gas speciation and 813C values at both sites indicate that the shallow sediment gas is primarily biogenic-sourced CH4, absent of higher molecular weight gases [22]. In this study a carbon budget is constructed for each site using hydrocarbon gas concentrations and stable carbon isotope (813C) and radiocarbon isotope (A14C) signatures of the

organic and inorganic carbon pools in solid phase sediment and pore water. The hypothesis is that dissolved inorganic carbon (DIC) assimilated during AOM is a dominant contributor to on-mound shallow sediment organic carbon cycling.

2. Materials and Methods

2.1. Study Location

Atwater Valley is a shallow trough in the Mississippi Canyon, Gulf of Mexico (Figure 1a). This underwater trough is part of the Mississippi Fan Fold-belt with sediment features including basinward-verging anticlines and underlain southern verging thrust faults across a 300 km long and 50 km wide region. Fold strata were formed during Late Jurassic to Miocene geologic periods resulting in formation of substantial salt tongues and sheets [24]. Canyon fill through the gas hydrate stability zone (HSZ) is comprised of fine-grained sediments, mostly interbedded debris flows and hemipelagic sediments overlain by a fine Holocene pelagic drape [25]. The area is characterized by seafloor mounds and basins formed through vertical advection of hydrocarbons, minerals, and CH4-rich fluids [26].

W94°30' W92°30' W90°30' W88°30'

N 31°30 N 30°30 N 29°30' N 28°30' N 27°30'

■ if ' ' ! . j * ^Bf — ■■'• -sa '■fWt

'jm^Mmf^mm

- \ N \ N

\ -. -«A off moundx corc^ A 1 km t mound ij M\C7

O I J oo

s 1400 18mmm gams«

& 1500 >

™«JB

^Miffip

USGS Seismic LineAV65

Figure 1. (a) Atwater Valley, Gulf of Mexico, Texas-Louisiana Shelf coring location (provided by Google Earth); (b) Sediment floor contours through the sample region, with sites on- (cores 3,7,8) and off-mound (core 2), the solid line represents the seismic data transect through the coring locations; (c) Core sites selected for this study were over a strong sediment bottom reflection (BSR) off-mound and on the mound observed with an elevated BSR and vertical seismic blanking above the BSR mound.

A 4 m high seafloor mound, at a water column depth between 1296 and 1300 m, was chosen as a study site based on previous seismic and geochemical studies conducted in May 2004 [22,23,27]

(Table 1). A key mound feature (Figure 1b) is a shallowing bell-shaped bottom simulating reflection (BSR) that raises approximately 200 m relative to off-mound BSR at 240 mbsf (Figure 1c, Table 1). At this location, heat flow elevates from a background signature off-mound of 40 to 160 mWm-2 on-mound (Table 1). These geophysical conditions on mound were interpreted to be a thermal perturbation to the HSZ base, creating vertical fluid fluxes [23].

Porewater chloride (Cl-) concentrations on-mound averaged 934 ± 74 mM, well above the seawater background, indicating active vertical advection of higher salinity porewaters (Table 1) and a rise in the BSR (Figure 1c). This interpretation is supported with the observation of Late Jurassic to Miocene salt tongues and salt diapirism in this region that would reduce deep sediment hydrate stability [24]. An advective CH4 flux of 3250 mMm-2year-1 was estimated on the mound (Table 1). In contrast, 750 m off-mound (Figure 1b) a moderate diffusive CH4 flux of 20.6 mMm-2year-1 was observed (Table 1). Off-mound, sediment porewater SO42- concentrations were observed to decrease linearly from near seawater values at the sediment-water interface to below detection limits at 410 centimeters below the sea floor (cmbsf). Below 410 cmbsf, sediment headspace methane concentrations increased linearly, indicating a sulfate methane-transition (SMT) supported by the anaerobic oxidation of methane (AOM). Vertical fluid advection limited downward SO42- diffusion on-mound (Table 1).

Table 1. Review of previously published data for core sites on- and off-mound. Where data is general for the mound and not located at specific stations it is listed under C7.

Parameter Cores

Off Mound Mound

C2 C7 C3 C8 Ref.

Latitude 27°56'26.644 27°56'26.644 27°56'26.644 27°56'26.644 [27]

Longitude 89°16'49.696 89°16'49.696 89°16'49.696 89°16'49.696 [27]

Water Column Depth (m) 1300 1296 1301 1296 [27]

Core Pentetration (cmbsf) 470 867 515 370 [27]

SMT (cmbsf) 410 no SMT no SMT 59 [22]

Methane Flux (mM CH4 m-2year-1) 20 3250 ND 167 [22]

Average Cl- (mM) 551 ± 6 934 ± 44 911 ± 79 770 ±102 [22]

Heat Flow (mWm-2) 40 160 - - [22]

BSR depth (mbsf) 240 40 - - [23]

Close proximity in spatial variation of downward SO42- diffusion and inferred vertical CH4 flux observed in previous studies on this mound makes it a unique study site to investigate CH4 contribution to shallow sediment carbon cycling [22]. Porewater and sediment inorganic and organic carbon concentration and carbon isotope ratio data from an off-mound core (C2) are compared to an on-mound core (C7). Results from two additional sediment mound cores (C3 and C8) located near C7 are used to provide additional data for interpretation of CH4 contribution to sediment carbon pools (Table 1).

2.2. Sediment Core Collection and Processing

Piston cores C7, C3, and C8 were collected on-mound (water depth = 1296 m) and C2 was collected off-mound (water depth = 1300 m, Figure 1b,c). For on-mound cores C7 was located near the

center, while C3 and C8 were near the mound edge, 150 and 70 m, respectively, from C7. Sediment cores were collected and processed shipboard as described in Coffin et al. [22]. Briefly, sediment cores were obtained using a 10 m piston coring system with 2.75 polycarbonate core liners. Cores ranged in length from ~300 to 800 cm and were processed immediately onboard the ship.

Core liners were inspected for gas pockets and gas expansion voids. At void spaces, the liner was drilled and gas sampled with a 60 mL polypropylene syringe fitted with a modified 3-way stopcock. Gas samples were then transferred to 30 mL pre-evacuated, glass serum vials fitted with a gastight stopper and aluminum seal. Subsequently, sediment plugs were collected from regular intervals along the core using a 3 mL polypropylene syringe with the tip cut off, transferred to pre-weighed 20 mL serum vials, and capped with gastight stoppers and aluminum seals to determine sediment headspace light hydrocarbon concentrations; (CH4 through C3H8) as well as 513CCH4(g) ratios.

For additional sampling, core liners were removed and cut in 10 cm sections within an interval of 25-45 cm. Wet sediment from each section was frozen in snap-tight Petri-dishes for laboratory measurements of sediment porosity and percent organic carbon. Porewater pressed from sediment using 70 mL Reeburgh-style PVC press containers pressurized to 400 KPa (60 psi) by a low-pressure air on a latex sheet between core sections and press gas inflow was collected into 60-mL polypropylene syringes [28]. Porewater was filtered from syringes through ashed Whatman GF-F filters into ashed 20 mL vials and subsequently distributed into appropriate vials for each analysis; 2 mL in a 5 mL glass serum vial for [DIC]; 1 mL in a 2 mL glass serum vial for 813Cdic, and; 2 mL in a 5 mL glass screw-top vial for dissolved organic carbon concentration [DOC] and 813Cdoc. Pressed sediment for inorganic and organic carbon concentration and isotope analyses was wrapped in ashed aluminum foil, sealed in Whirlpack bags, and stored frozen at -20 °C for analyses at the land-based laboratory.

2.3. Shipboard Analyses

To extract volatile hydrocarbons from sediment into vial headspace for gas analysis, 3 mL of nitrogen sparged, deionized water was injected through the septum of the serum vial and vial was shaken for 3 min. After this extraction, the headspace sample was removed from the vial and injected into a sampling loop on a Shimadzu GC-14A gas chromatograph-flame ionization detector (GC-FID) with a Hayesep-Q packed column (Alltech, Deerfield, IL, USA) to measure CH4, ethane (C2H6), and propane (C3H8) concentrations. Sediment CH4 concentrations were corrected for atmospheric background in the vials (95% extraction efficiency was assumed). Core gas pocket C1-C3 alkane concentrations were also measured using the GC-FID. Analytical precision was within 0.1 mM, based on replicate analyses. The limit of detection for methane was 0.009 mM and where concentrations were lower data are presented as 0.0 mM.

Porewater DIC concentrations were determined using a UIC CO2 coulometer (UIC, Inc., Joliet, IL, USA) standardized to a certified seawater reference material (University of California, San Diego, CA, USA). Replicate variability was less than 0.15 mM.

2.4. Post-Cruise Laboratory Analyses

Sediment total carbon and OC (%TC, %SOC) concentrations and 813C values were determined on a Fisons EA 1108 C/H/N analyzer in line with a Thermo Electron Delta Plus XP Isotope Ratio Mass

Spectrometer (IRMS) interface via a Conflo II (Thermo Scientific, Waltham, MA, USA). Pressed sediment was dried at 80 °C, ground with a mortar and pestle, then 15 to 20 mg of sediment was weighed in tin capsules for TC analysis. For SOC analysis, sub-samples were weighed in silver capsules, treated with an excess of 10% HCl and dried in an oven at 70 °C overnight to remove inorganic carbon. Sediment inorganic carbon (%CaCO3) concentrations were determined from the difference between TC and SOC. A concentration calibration curve for carbon concentration analysis was generated daily by analyzing an acetanilide standard. For sediment 813Csoc values, IAEA-C8 (oxalic acid), IAEA-CH-6 (sucrose) and USGS 40 (l-glutamic acid) were used as calibration standards. Acetanilide standards (USGS-40 and IAEA-C8) were also used as check standards during analysis. All 813C data presented in this work are in per mil units (%o) and referenced to the Vienna Pee Dee Belemnite (VPDB) scale. Errors were based on triplicate runs. Error for % SOC was within ±0.03%C, %CaCO3 varied by less 0.2%, and 813Csoc varied by less than 0.2%.

Sediment 813CcaCO3, pore water 813Cdic, and gas pocket and sediment S13Cch4 ratios were determined using a Thermo Electron Trace GC equipped with a Varian Porapak-Q column and GC-CIII combustion interface in-line with the Delta Plus XP IRMS (Thermo Scientific, Waltham, MA, USA) [22]. For 813CcaCO3 analysis, 250 mg of sediment in a serum vial was treated with 2 mL of 10% HCl. For 813Cdic analysis, 2 mL porewater samples were treated with 200 pL of 85% H3PO4. In both cases, CO2 was extracted from the vial headspace and injected into the GC via a split/splitless inlet in split mode. All 813CcaCO3 and 813Cdic values were normalized through analysis of CO2 and C1-C5 alkanes in NIST RM 8560 (natural gas, petroleum origin). Samples for S13Cch4 analysis were introduced via an in-line cryogenic focusing system according to the method of Plummer et al. [29]. A separate 813C normalization curve was generated for C1-C4 alkanes and used to normalize S13Cch4 data. Replicate 813CcaCO3 values varied by less than 0.2%, 813CDIC by less than 0.5%, and 513CCH4 by less than 1.0%.

Porewater DOC was measured on an OI Analytical 1010 total organic carbon analyzer (OI Analytical, College Station, TX, USA) using a heated persulfate oxidation method modified for seawater analyses [30]. The samples were kept frozen until ready for analysis, then acidified and nitrogen sparged in the lab prior to analysis to remove DIC. The total organic carbon analyzer was interfaced with the Delta Plus XP IRMS. A DOC calibration curve was generated using standards of potassium hydrogen phthalate. Measured 813Cdoc values were normalized to the VPDB scale by analyzing solutions of IAEA-CH-6 (= -10.449 ± 0.0339%) and USGS 40 (-26.389 ± 0.042) standards. Using triplicate analyses, DOC concentrations varied by less than 0.02 mM, and 813Cdoc varied by less than 0.2%.

2.5. Radiocarbon Isotope Analysis

Graphite sample preparation for A14C analysis of CH4 and SOC are described in detail by Pohlman et al. [31]. Targets were prepared in the U.S. Naval Research Laboratory (NRL) Graphite Lab. Targets were then analyzed at the NRL accelerator mass spectrometer (AMS) facility using an AMS equipped with a high intensity cesium sputter source for 14C analysis [32]. Data analysis was according to standard procedures described in Tumey et al. [33]. Samples were measured against OX II standards and blanks distributed throughout loadings of the sample wheel. On average, each wheel contained 2 AMS blanks for tuning the accelerator, three processing blanks appropriate for the samples on the wheel, and 7 OX II standards. 813C results were obtained from GC-IRMS analysis of the samples.

These values were used to calculate the 8C fractionation correction for each sample. The A14C data were calculated as:

A14 C =

x 1000

where R is 14C/13C or 14C/12C, Rsn represents the sample with normalization to a standardized 513C and Rabs represents oxalic acid standards normalized to a 513C standard and the standardized atmospheric 14C level in 1950 [34-36]. To compare the influence of C cycling with aging of deposited sediment, A14C was converted to conventional radiocarbon age (CRA) where:

t = - 8033

f A14C ^ +1

v 1000

+ X (y - 1950)

In this equation, t is the CRA in years before present (1950), using the Libby half-life for 14C (t1/2 = 5568 years). To obtain actual calendar years, a standardized conversion was applied. The A14C term is the value obtained after the 13C correction is applied, X is the 14C constant (1/8267 years), and y is the year of measurement. Thus, for A14C of 0, t would be 0 for a measurement performed in 1950 and 1955 for 2007, when the measurement was performed [33].

2.6. d3C Data Interpretation

The relative contribution of CH4 to the shallow sediment carbon pools on-mound with a high advection and off-mound with a moderate diffusion assumes microbial assimilation of DIC during AOM [37]. Methane contribution at these different flux rates is summarized at specific core depths or averages through the sulfate methane transition zone. Accounting for isotope fractionation during assimilation, DIC contribution to organic carbon pools is estimated using a simple, two end-member isotope mass balance [38]:

Rx = ^PD^PD + ^dic^DIC (3)

where Rx represents 813C of SOC or DOC and the isotopic composition of marine phytodetritus (PD) and DIC are represented by Rpd and Rdic, respectively. The variables Cpd and Cdic represent the corresponding fractional contributions of marine PD and DIC such that:

CPD + CDIC = 1 (4)

The percent contribution of DIC to each carbon pool (%X) can then be derived from Equations (3) and (4):

Rr fipn

= X 100 (5)

where Rs represents the isotopic composition of DIC. Consideration of isotope fractionation during DIC assimilation is presented in the Discussion.

3. Results

Results focus on a comparison of off-mound (C2) and on-mound (C7) CH4 contribution to organic carbon in gas, porewater and sediment samples. Core C3 provides supplementary data for assessing the

on-mound CH4 flux and cycling. Core C8, located near the edge of the mound, shown to have active vertical CH4 diffusion [22], is also used in data interpretation.

3.1. Radiocarbon

The A14CCH4(g) values measured in gas pockets of near-mound cores C3 and C8 taken close to the primary mound core are presented in Figure 2 and Table 2. The gas pockets in core C8 were slightly enriched in 14C relative to gas pockets in core C3. On-mound SOC in core C7 was substantially more 14C-depleted than SOC in off-mound core C2 (Figure 2). In each core, the shallowest sediments had a more modern radiocarbon age, however, there was a large difference in the A14Csoc data between the cores. To estimate a change in apparent sediment age caused by the presence of CH4, a sedimentation rate of 0.037 ± 0.022 cm-year-1 [3] and natural radiocarbon decay starting at shallow sediment off-mound A14C value (-283%) was applied to data to construct a conservative age line (Figure 3). CRA in off-mound and on-mound SOC did not conform to the predicted conservative age line. In the shallow section of the off-mound core, carbon age was older than predicted aging line down to approximately 300 cm and was younger below 300 cm. On-mound sediment CRA was older than the age line through the core.

A14C (%o)

-1000 -800 -600 -400 -200

0 ■ cP □ o o o

a o 200 ■ 400 ■ « □ □ v □ □ V • V Core 3 CH4 Core 8 CH4

£ □ □ mound SOC

ft <u 600 - o off mound SOC

Q 800 - □ □ □ □

1000 J

Figure 2. A14Cch4 of free gas and SOC at on-mound and off-mound sites.

Age (Years)

0 5000 10000 15000 20000 25000 30000

Figure 3. Conventional radiocarbon age of SOC for cores taken off- and on-mound. The conservative aging line shown was calculated using estimates of annual sedimentation rate for this region.

3.2. Solid Phase Sediment Profiles

An overview of all carbon pool concentration and S13C data is presented in Table 2. A complete data set is available in Supplementary. Sediment OC concentrations were higher in off-mound C2 core than on-mound C7 core (Figure 4a, Table 2). On-mound S13Csoc values showed strong 13C-depletion near the seawater interface (SWI, Figure 4b). Off-mound SOC was more 13C-enriched with higher values near surface above the SMT of 410 cmbsf (Figure 4b, Table 1). Profiles of CaCO3 (Figure 4c) and S13CCaCO3 (Figure 4d) also were substantially different when compared on- and off-mound. CaCO3 as a percentage of total sediment mass was substantially higher on-mound (average = 11.4 ± 1.1%, n = 16) than off-mound (average = 4.1 ± 3.0%, n = 19). Sediment S13CCaCO3 values were generally 13C-enriched on-mound (average = -2.7 ± 2.8%, n = 16) relative to off-mound (average = -4.5 ± 2.4%, n = 19). However, at the SWI, on-mound S13CCaCO3 value was 13C-depleted while off-mound S13CCaCO3 was 13C-enriched.

Organic Carbon (%

0.2 0.4 0.6

off mound mound

"I 400 £

& 600 H Q

SOC (%0)

off mound mound

CaCO3 (%)

g 400 £

& 600 Q

1000 J

8 10 12 14 16

off mound mound

CaCO3 1 -14 -12 -10 -8 -6 -4

■2 0

& 600 Q

800 1000

Figure 4. Sediment profiles at on-mound and off-mound sites: (a) Percent sediment organic carbon (SOC); (b) 813Csoc; (c) percent CaCO3; and (d) 813Ccaco3.

Table 2. Summary of C pool concentrations and S13C, with minimum, maximum and average values for each core. The full data set for each core is available in Supplementary.

Gas Porewater Sediment Core ID - - -

_813Cch4 A14Cch4 CH4 (mM) 813Cch4 DIC (mM) 813Cdic DOC (mM) 813Cdoc % SOC 513Csoc A14Csoc % CaCO3 813Ctic

Min ND

Max ND

AVG ± SD ND C3

Min -72.5 -963

Max -72 -957

AVG ± SD -72.3 ± 0.2 -960 ± 3 C7 Min Max AVG ± SD

Min -72.3 -905 0.0 -85.2 5.5 -44.3 ND ND ND ND ND ND ND

Max -65.8 -901 13.0 -71.5 13.8 6.7 ND ND ND ND ND ND ND

AVG ± SD -69.7 ± 3.4 -903 ± 2 5.64 ± 3.62 -76.6 ± 3.9 10.6 ± 2.0 -15.0 ± 21.7 ND ND ND ND ND ND ND

ND LOD -86.9 2.8 -47.1 0.82

ND 12.5 -48.8 20 -9.4 3.64

ND 5.3 ± 5.5 -74.1 ±11.5 10.2 ± 5.3 -27.6 ± 11.5 2.06 ± 1.18

-26.2 0.91 -23 -679 2.4 -6.1

-20.1 0.64 -20.5 -283 13.04 1.8

-23.9 ± 1.5 0.72 ± 0.07 -22.5 ± 0.6 -629 ± 88 4.4 ± 3.2 -4.6 ± 2.4

0.1 -75.4 6.7 -42.9 ND ND ND ND ND ND ND

7.9 -72.8 13.6 3.7 ND ND ND ND ND ND ND

4.2 ± 1.8 -73.7 ± 0.9 9.7 ± 2.3 -3.0 ± 12.6 ND ND ND ND ND ND ND

-71.4 ND 3.3 -82 5.6 -48.1 2.05 -27.6 0.57 -28.2 -951 9.89 -12.6

-71.2 ND 10.6 -69 12.1 4.6 3.92 -24.8 0.66 -25 -890 13.04 -1.1

-71.3 ± 0.3 ND 6.2 ± 1.8 -70.4 ± 3.3 7.3 ± 1.6 -3.4 ± 14.4 3.55 ± 0.42 -25.6 ± 0.8 0.60 ± 0.03 -25.8 ± 0.8 -929 ± 22 11.4 ± 1.1 -2.7 ± 2.8

3.3. Gas Sources and Sediment and Porewater Carbon Profiles

Through the discussion we assume that shallow CH4 originates from deep sediment CH4 fluxes; deep CH4 is based on on-mound core liner gas pocket data. Assuming void gas from core liner pockets is a deeper source is supported with depleted A14Cch4 data (see Section 3.1. Radiocarbon), coupled with observation of elevated Cl- profiles (Table 1) indicating deep vertical advection to the shallow system. Further support for this interpretation is presented below with on-mound porewater DIC data. This 813Cch4 value is compared to sediment headspace methane (CH4(g)) to assess shallow sediment cycling (Table 2). Measured 813CcH4(g) values from on-mound core C7 and near-mound core C3 showed little variation within and between cores. Variation of 813CcH4(g) was greater and moderately depleted in 13C in near-mound core C8. For all cores with gas pockets in which ethane was detected, the C1/C2 ratios were high ([27]; average 11,000, n = 13). As a note, all CH4 concentrations shown in Table 2 are relative (headspace) measurements and underestimate actual sediment methane concentrations as a result of pressure changes during core retrieval from the ocean floor resulting in degassing [21].

Off-mound (C2) CH4 concentrations were highest below the apparent SMT at 410 cmbsf (Table 2) and near the limit of detection above the SMT (Figure 5a). On-mound CH4 concentrations were higher in shallow sediments (C3, C7, and C8), relative to off-mound, and showed a general decrease toward the sediment-water interface (SWI). On-mound sediment 813Cch4 was relatively uniform through the profile except for 13C depletion at the SWI (Figure 5b). Mound core C3 had a similar profile with moderate 13C depletion through the profile. More variation was observed in on-mound core C8, with 13C-depletion observed at 50 cmbsf (Figure 5b). Off-mound sediment CH4 was 13C-enriched in shallow sediments up to 260 cmbsf and was depleted deeper than 300 cmbsf (Figure 5b). Off-mound 813Cch4 data are not presented above 260 cmbsf because CH4 concentrations were below the limits of detection for carbon isotope analyses.

On-mound (C7) porewater DIC concentrations increased from the SWI to 34 cmbsf, and subsequently decreased towards the core base (Figure 6a). Off-mound (C2) porewater DIC concentrations gradually increased from a SWI minimum to a maximum at 350 cmbsf and then declined rapidly toward the core base. On-mound, 813Cdic values decreased to most depleted 13C value at 34 cmbsf where concentration increased; below this point 813Cdic values increased and remained uniform (Figure 6b). Off-mound 813Cdic values showed a similar patter with the minimum value at 330 cmbsf.

While on-mound and off-mound porewater DOC concentration ranges were similar, notable differences in vertical profiles were observed (Figure 7a). Off-mound DOC concentrations in porewaters were lowest in near surface sediments and generally increased with depth. On-mound porewater DOC concentrations were relatively consistent throughout the profile, except for a low value near surface. On-mound porewater 813Cdoc values showed little variation, however, porewater DOC was substantially 13C-depleted throughout the core with a minimum value at 34 cmbsf (Figure 7b). In contrast, off-mound porewater 813Cdoc values varied substantially and were generally 13C-enriched; S13Cdoc values were elevated in samples taken less than 100 cmbsf and depleted below 100 cmbsf (Figure 7b).

CH4 (mM)

0.0001 0.001 0.01

-•— off mound

— mound core 3 -r— mound

— mound core 8

s13cch4

off mound mound core 3 mound mound core 8

Figure 5. (a) Sediment CH4 concentrations; and (b) 813Cch4 values measured for on-mound (C7) and off-mound (C2) core locations. On-mound cores include cores 3 and 8 that are used for estimating the C isotope values for source CH4 and interpretation of mound C cycling.

DIC (mM) a S13CDIc (o/oo) b

0 5 10 15 20 25 -60 -50 -40 -30 -20 -10 0 10

S 400 -

H 600 -

800 -1000 -

Figure 6. (a) Sediment pore water dissolved inorganic carbon (DIC) concentrations; and (b) S13Cdic values measured on-mound and off-mound.

1000 J

DOC (mM) 2 3

off mound mound

DOC (%o)

-24 -22

Figure 7. (a) Vertical profiles of sediment pore water dissolved organic carbon (DOC) concentrations compared on-mound and off-mound; and (b) Vertical profiles of 813Cdoc values measured on-mound and off-mound.

4. Discussion

4.1. Shallow Sediment Carbon

Shallow sediment carbon cycling in Atwater Valley was investigated assuming marine phytodetritus and AOM carbon fixation are the primary sources to SOC and DOC. The difference in A14Csoc values between cores on- and off-mound is likely a result of AOM incorporation of isotopically-depleted DIC into OC [11,12,37], with a substantially greater CH4 influence on-mound (Figures 2 and 3). There is a pronounced 14C-depletion in SOC on-mound (-955% to -890%) with a slightly more modern signature in surface sediments. Elevation in A14C observed in on-mound (C8) surface sediment results from more modern seawater DIC fixation during AOM, further discussed below. Note that on mound 813Cdic was relatively uniform near 0% (Figure 6b) and represents a deep sourced signature associated with the fluid advection observed in the mound 813Cch4 (Figure 5b). An alternate interpretation of on-mound aged carbon data could be sediment erosion resulting in uncovering relic SOC. While erosion could contribute to the aged carbon pattern on-mound, elevated CaCO3 concentrations relative to off-mound values indicate active AOM, resulting in over saturated

DIC and subsequent precipitation (Figure 4). In contrast to on-mound data, off-mound SOC was substantially more modern (-678% to -283%, Figure 3) implying a lower contribution of AOM DIC fixation to carbon cycling. Gordon and Goni [3] reported surface sediment A14Csoc ranging from -309.1% to -228.6% in the same general sampling area, with water column depths of 365 to 2270 m. In this study, similar A14Csoc on- and off-mound profiles with more 14C-depletion observed on-mound does indicate depleted radiocarbon DIC contribution to carbon cycling (Figures 2 and 3). The near-constant vertical CRA of the sediments sampled in this study, especially off-mound, suggests sediment mixing or rapid deposition, perhaps created by a shelf slump (Figure 3). Subsequent diagenesis of organic matter with different contributions of 14C-depleted CH4 between core locations is evident with differences in CRA between on-mound and off-mound cores. 813C data in this study provides an estimate of the DIC fixation driven by the AOM of deep sourced CH4 in shallow sediment carbon cycling. Other investigators have shown that Mississippi River particulate organic carbon (POC) changes from a terrigenous (513C = -28% to -26%) to phytoplankton source (513C = -19%) near the river mouth [1,3,39,40]. In a previous study, POC-S13C reported well offshore near Atwater Valley ranged from -22.5% to -18.7%, [41], in the range of phytodetritus [42,43]. For our study off-mound S13Csoc values (-20.5% at the surface to -23% near the core bottom, Figure 4) were characteristic of SOC dominated by phytodetritus. In contrast, on-mound S13Csoc ranged from -25.0% down core to -28.2% near surface, indicating an alternate carbon source.

A common interpretation in 13C-depleted isotope signatures in coastal waters is carbon sourced from terrestrial plants [39]. However, more recent studies show 13C-depleted S13Csoc values observed in these sediments likely results from carbon assimilated into bacterial biomass during CH4 cycling [22,37]. In anoxic sediments bacterial biomass is incorporated and preserved in solid phase sediments [44,45] and through time can constitute a significant portion of the SOC pool [46,47]. While our data indicate SOC is initially derived from marine phytoplankton, once deposited it is subject to diagenesis which includes incorporation of deep sediment CH4 into solid phase sediment and pore water carbon pools.

Several additional observations pertaining to sediment carbon pools support interpretation of CH4 contribution to carbon cycling. In the presence of sediment anaerobic CH4 oxidation DIC oversaturation results in formation of CaCO3 [48,49]. On-mound where 13C depleted SOC was observed, CaCO3 concentrations were 3 to 4 times higher than off-mound (Figure 4c). On-mound and off-mound data also suggest CH4 contribution to pore water DOC. On-mound core porewater S13Cdoc values ranged from -27.6% to -24.8% with 13C-depletion in shallower sediments, though porewater DOC concentrations showed little variation (Figure 7). The shift to lower S13Cdoc values in shallow on-mound sediments coincided with a decrease in pore water SO42- concentrations and apparent AOM [22]. Off-mound porewater DOC concentrations below the SMT were similar to those on-mound and had depleted S13Cdoc (-26.2%; Figure 7) indicating a CH4 contribution to porewater DOC. In off-mound surface sediments above the SMT, porewater DOC concentrations were lower with an increase in S13Cdoc (-20.1%) more typical of a marine phytoplankton source. Here it is assumed that DO13C depletion is not a result of diagenesis or selective metabolism during heterotrophic microbial carbon cycling that would result in isotopic fractionation of the total DOC pool [50,51].

4.2. Methane Source and Cycling

Composition and isotopic signatures of gas pocket samples on-mound indicate a microbial CH4 source with average C1/C2 = 11,000 (methane/ethane) and ô13CcH4(g) = -71.2% (Table 2) [8,52-54]. While data interpretation is contingent on piston core depths, there are indications that on-mound CH4 originates from deep sediment; (1) A14CCH4(g) depletion measured in on-mound (C3) (-961%) is a value characteristic of CH4 that originates in from a deep system [55]; and (2) SO42- not observed below the SWI in cores C3 and C7 suggests active advection of deep sediment CH4 to the surface [22]. In contrast, 14C-enrichment indicated by A14CCH4(g) values in on-mound (C8, -903%, Figure 2) and presence of a shallow SMT (59 cmbsf) suggests modern seawater DIC is being reduced to CH4 during methanogenesis. A biogenic CH4 source, absent of higher molecular weight gases and oil products that are present in a thermogenic source [53], presents an opportunity to study deep sediment CH4 contribution to shallow sediment carbon cycling. On mound core C3 gas pocket carbon isotope ratios is used as a CH4(g) end-member to estimate CH4 contribution to organic carbon pools, presented below.

Strong difference in shallow sediment CH4 cycling was observed within and between cores. Measured on-mound sediment headspace ô13CcH4 values (-70.4%) deeper than 100 cmbsf were similar to the ô13CCH4(g) value indicating little microbial CH4 consumption or production (Figure 5b). In contrast, in the top 100 cm on-mound, sediment headspace ô13CcH4 values were depleted 11% relative to ô13CCH4(g) values (Figure 5, Table 1). The 13C-enrichment indicated by ô13CCH4(g) values in C8 (-65.8%) at 240 cmbsf likely results from isotopic fractionation during AOM [22,56]. Where the upward CH4 flux due to active vertical fluid advection (such as C7 and C3) impedes downward SO42-diffusion into sediment, AOM or methanogenesis will not occur and there is not a change in the ô13CCH4(g) ratio. Depleted ô13CcH4 values on-mound near surface and corresponding enriched 813Cdic values do indicate active methanogenesis [57,58]. Methanogenesis typically occurs in marine sediments below AOM [59,60]. However, with advection on-mound, 513Cdic data indicate shallow methanogenesis near the SWI (Figure 5b). While enriched ô13CcH4 and an associated SMT was not observed there is potential for seawater SO42- to support AOM at the SWI (Figure 5b).

Active CH4 cycling was also observed in the off-mound core. Above the BSR (Figure 1) there was greater variation in the ô13CcH4 profile with lowest value (-86.9%), potentially resulting from methanogenesis at a mid-core depth (Figure 5b), just below the SMT (Table 1). Note in a recent study carbon isotope equilibration at the point of sulfate-limited anaerobic oxidation was observed and may account for this depleted ô13CcH4 [61]. Above the SMT, decreased CH4 concentrations and enriched sediment headspace ô13CcH4 values up to -48.8% result from AOM through the SO42- gradient (Figure 5). A corresponding depletion in 513Cdic values to -47.1% (330 cmbsf) and a subsequent increase to -10.5% towards the core base (Figure 6) is characteristic of AOM through the SMT located above the zone of methanogenesis [57,62].

4.3. Estimation of CH4 Contribution to the Shallow Organic Carbon

The relative contributions of phytodetritus and AOM to shallow sediment carbon cycling on-mound and off-mound can be estimate using a simple carbon budget calculation (Equation (5)). We assume sediment CH4 is primarily oxidized during AOM for energy and DIC is assimilated into cellular

biomass [11-13] with an isotopic fraction of 3.75% , a value recently applied in a study on the Hikurangi Margin [37] and an intermediate (2.0% to 5.5%) for isotopic fractionation during DIC assimilation in the reversed tricarboxylic acid cycle [63-65]. For estimating the DIC contribution to SOC during AOM it is assumed 813Cpd end-member (Equation (5)) was -20.5%, the off-mound signature in surface sediments with the most modern A14Csoc value (Figure 2). Use of this 813Cpd end-member is supported by other studies on the TX-LA Shelf where offshore to nearshore sediment 813C values ranged from -21.7% to -19.7% [1,3,40]. On the Gulf of Mexico abyssal plain, with no substantial terrestrial inputs, 813Csoc measured values were -20.6% [66]. The average porewater 813Cdic value measured on-mound, through the SMT (between 34 and 144 cm) was -29.5%. Adjusting this value for fractionation during DIC fixation (3.75%, [37]) provides a 813Cdic end member of -25.8%. On-mound, with a 813Csoc average of -25.0%, results in an AOM-DIC contribution to SOC of 85% (Equation (5), Figure 8). We assume, with on-mound advection, AOM occurs at the sediment water column interface where SO42- is available and this signature is buried through time. This observation of high CH4 contribution to shallow sediment carbon cycling is supported with a comparison of regional sedimentation rates and on-mound vertical CH4 flux. Phytodetritus contribution to SOC in this region is estimated on the basis of a sedimentation rate near this location of 0.005 cm-year-1 [1], sediment density of 2.6 gmcm-3 [67], and SOC mineralization of shallow off-mound sediment relative to deep on-mound sediment of 37% [22,68] to be 40.8 mMCm-2year-1. This suggests a small contribution of carbon relative to CH4 at 3250 mMCm-2year-1 (Table 1). Applying this approach to off-mound sediment 813Csoc values below 35 cmbsf averaged -22.7% and results in a 42% contribution of AOM-derived DIC incorporation to SOC (Figure 8).

Porewater DOC concentrations and 813Cdoc also indicate different carbon sources between on- and off-mound cores (Figure 7). In this analysis we do assume that porewater bacterial production, cell DOC excretion, and degradation of SOC contribute to porewater DOC. Assuming DIC is fixed into bacterial carbon and subsequently cycled to DOC (Figure 8), the contribution of AOM to DOC (Equation (5)) can be estimated using the phytodetritus end-member 813Cpd of -20.5% and 813Cdic of -25.8%. The 813Cdoc average off-mound (-24.3%0), suggests 72% of DOC below the SMT is AOM derived (Figure 8). The on-mound average 813Cdoc value of -25.6% through the entire C7 core indicates an indirect contribution of AOM to porewater DOC of 96%. The production of DOC via AOM in sediment may result in significant contribution relative to phytoplankton derived carbon. It is important to note that patterns in 813Cdic and 813Cdoc depletion were similar in the vertical core profiles and with high porewater DOC concentrations ogranoclastic degradation of DOC and subsequent production from light DIC from further AOM can result in an overestimate of AOM related DOC production. Furthermore acetogenesis can contribute to the observation of depleted 13C in the DOC [69]. Also, it is interesting that this observation is similar to a recent study off the coast of New Zealand where AOM contribution to porewater DOC was estimated to be up to 71% at a location with lower vertical CH4 flux rates [37].

Water Colum

Phytoplankton Sedimentation

Deep Sediment

S13CDIC = -25.8%o

CH4 Diffusion

CH4 Advection

off mound

Figure 8. A summary of dissolved organic carbon (DIC) and phytoplankton biomass contributions to shallow SOC and DOC at the on-mound and off-mound sites using 813C values.

5. Conclusions

Diffusive and advective fluxes of deep sediment CH4 contribute significantly to shallow sediment carbon cycling in Atwater Valley, Gulf of Mexico. Deep salt diapirs at this study location are predicted to result in destabilization of CH4 hydrates and result in an elevated vertical methane flux [22,23]. This prediction is supported in the observation of a substantial rise in the BSR (Figure 1c).The upward flux of deep sediment CH4 via fluid advection on-mound supports a substantial fraction of sediment carbon production (up to 85%), whereas marine phytoplankton sourced carbon contributes ~50% where there is a low diffusive CH4 flux (Figure 8). Consistent 13C-enriched 813Cdic values deeper than 100 cmbsf and the 14C-depleted A14Csoc and A14Cch4 values measured in gas pockets on-mound can be interpreted as long-term mound formation driven by active advection that supported shallow sediment carbon cycling driven by AOM near the SWI where SO42- is abundant. The cycling of CH4 via AOM is reflected in the organic and inorganic carbon pools. The variation in CH4 cycling in an advection dominated on-mound sediment and a diffusion dominated sediment off-mound is evident in the same sediment carbon pools. Observed 813Cdic profiles followed changes in 813Cch4 values, as did 813CcaCO3 values, indicating oxidation of CH4 during AOM was occurring in the sediment pore fluids on the mound. Mass balances using 813C of CH4 and SOC showed a large difference in the estimated contribution of deep sediment CH4 to the shallow sediment carbon pools off and on-mound relative to marine phytoplankton. The DOC pool in sediment porewaters also reflected a significant contribution (up to 96%) from deep sediment CH4. In another study on the TX-LA Shelf, petroleum seeps were found to contribute 40% to 60% to the total organic C in a shallower slope region [17].

This study suggests a need for global consideration of the distribution of deep sediment CH4, especially hydrate bound methane, and the flux of carbon from this globally-significant carbon pool to shallow sediments. More estimates of deep sediment CH4 and petroleum contributions to shallow sediments and the water column will help refine marine carbon cycling models and budget estimates and improve predictions of the impacts of climate change on these reservoirs. The current estimate for the world coastal CH4 hydrate distribution is 21 x 1015 m3 of CH4 at standard temperature and pressure [8]. A more thorough understanding of the potential for fixation of CH4 into organic carbon and carbonate phases in shallow sediments will assist in evaluating deep sediment sourced CH4 fluxes through the shallow sediments to the water column and the atmosphere.

Supplementary Materials

Supplementary materials can be accessed at: http://www.mdpi.com/1996-1073/8/3/1561/s1. Acknowledgments

Warren Wood, NRL-Stennis Space Center, MS provided seismic profiles and seafloor topography. Ross Downer, Milbar Hydo-Test, Inc. (Shreveport, LA, USA) was lead for all coring operations. We appreciate the technical discussions and reviews of this manuscript from Paula Rose, Jeff Chanton, Thomas Boyd, and Leila Hamdan. David Knies contributed to radiocarbon analyses. We also appreciate the support by the crew of the RV Gyre. Finally, reviewers of this manuscript have provided excellent input to research and the presentation. This research was supported by Department of Energy-National Energy Technology Laboratory, Office of Naval Research and the US Naval Research Laboratory.

Author Contributions

All authors have contributed to writing and revisions of this manuscript. Conflicts of Interest

The authors declare no conflict of interest. References

1. Goni, M.A.; Ruttenberg, K.C.; Eglington, T.I. A reassesment of the sources and importance of land-derived organic matter in surface sediment from the Gulf of Mexico. Geochem. Cosmochim. Acta 1998, 62, 3055-3075.

2. Bianchi, T.S.; Mitra, S.; McKee, B.A. Sources of terrestrially-derived organic carbon in lower Mississippi River and Louisiana shelf sediments: Implications for differential sedimentation and transport at the coastal margin. Mar. Chem. 2002, 77, 211-223.

3. Gordon, E.S.; Goni, M.A. Controls on the distribution and accumulation of terrigenous organic matter in sediments from the Mississippi and Atchafalaya river margin. Mar. Chem. 2004, 92, 331-352.

4. Mayer, L.M.; Schick, L.L.; Allison, M.A.; Rutthenberg, K.C.; Bentley, S.J. Marine vs. terrigenous organic matter in Louisiana coastal sediments: The uses of bromine organic carbonratios. Mar. Chem. 2007, 107, 244-254.

5. Rabalais, N.N.; Atilla, N.; Normandeau, C.; Turner, R.E. Ecosystem history of Mississippi River influenced continental shelf revealed through preserved phytoplankton pigments. Mar. Pollut. Bull. 2004, 49, 537-547.

6. Ruttenberg, K.C.; Goni, M.A. Phosphourus distribution, C:N:P ratios, and 813Coc in arctic, temperate, and tropical coastal sediments: Tools for characterizing bulk sedimentary organic matter. Mar. Geol. 1997, 139, 123-145.

7. Anderson, R.K.; Scalan, R.S.; Parker, P.L.; Behrens, E.W. Seep oil and gas in Gulf of Mexico slope sediment. Science 1998, 222, 619-621.

8. Milkov, A.V. Molecular and stable isotope compositions of natural gas hydraes: A revised global data set and basic interpretations in the context of geological settings. Org. Geochem. 2005, 36, 681-702.

9. Borowski, W.S.; Paull, C.K.; Ussler, W., III. Marine porewater sulfate profiles indicate in situ methane flux from underlying gas hydrate. Geology 1996, 24, 655-658.

10. Pancost, R.D.; Damste, J.S.S.; de Lint, S.; van der Maarel, M.J.E.C.; Gottschal, J.C.; the Medinaut Shipboard Scientific Party. Biomarker evidence for widespread anaerobic methane oxidation in Mediterranean sediments by a consortium of methanogenic archaea and bacteria. Appl. Environ. Microbiol. 2000, 66, 1126-1132.

11. Knittel, K.; Boetius, A. Anaerobic oxidation of methane: Progress with an unknown process. Annu. Rev. Microbiol. 2009, 63, 311-334.

12. Alperin, M.J.; Hoehler, T.M. Anaerobic methane oxidation by archaea/sulfate-reducing bacteria aggregates: 2. Isotopic constraints. Am. J. Sci. 2009, 309, 958-984.

13. Wegener, G.; Niemann, H.; Elvert, M.; Hinrichs, K.-U.; Boetius, A. Assimilation of methane and inorganic carbon by microbial communities mediating the anaerobic oxidation of methane. Environ. Microbiol. 2008, 10, 2287-2298.

14. Kellermann, M.; Wegener, G.; Elvert, M.; Yoshinga, M.Y.; Lin, Y.-S.; Holler, T.; Mollar, X.P.; Knittel, K.; Hinrichs, K.-U. Autotrophy as a predominant mode of carbon fixation in anaerobic methane-oxidizing mincrobial communities. PNAS 2012, 109, 19321-19326.

15. Leak, D.J.; Dalton, H. Growth yields of methanotrophs. Appl. Microbiol. Biotechnol. 1986, 23, 470-476.

16. Kelley, C.A.; Coffin, R.B.; Cifuentes, L.A. Stable isotope evidence for alternate carbon sources in the Gulf of Mexico. Limnol. Oceanogr. 1998, 43, 1962-1969.

17. Wang, X.-C.; Chen, R.F.; Whelan, J.; Eglinton, T. Contribution of "old" carbon from natural marine hydrocarbon seeps to sedimentary and dissolved organic carbon pools in the Gulf of Mexico. Geophys. Res. Lett. 2001, 28, 3313-3316.

18. Joye, S B.; Boetius, A.; Orcutt, B.N.; Montoya, R.P.; Schulz, H.N.; Ericson, M.J.; Lugo, S.K. The anaerobic oxidation of methane and sulfate reduction in sediments from Gulf of Mexico cold seeps. Chem. Geo. 2004, 205, 219-238.

19. Paull, C.K.; Ussler, W., III; Lorenson, T.; Winters, W.; Dougherty, J. Geochemical constraints on the distribution of gas hydrates in the Gulf of Mexico. Geo. Mar. Lett. 2005, 25, 273-280.

20. Ruppel, C.; Dickens, G.R.; Castellini, D.G.; Gilhooly, W.; Lisarralde, D. Heat and salt inhibition of gas hydrae formation in the northern Gulf of Mexico. Geophys. Res. Lett. 2005, 3, L04605, doi:1029/2004GL021909.

21. Lapham, L.L.; Chanton, J.P.; Martens, C.S.; Sleeper, S.; Woolsey, J.R. Microbial activity in surficial sediments overlying acoustic wipeout zones at a Gulf of Mexico cold seep. Geochem. Geophys. Geosyst. 2008, 9, doi:10.1029/2008GC001944.

22. Coffin, R.B.; Hamdan, L.; Plummer, R.; Smith, J.; Gardner, J.; Wood, W.T. Analysis of methane and sulfate flux in methane charged sediments from the Mississippi Canyon, Gulf of Mexico. Mar. Pet. Geol. 2008, 25, 977-987.

23. Wood, W.T.; Hart, P.E.; Hutchinson, D.R.; Dutta, N.; Snyder, F.; Coffin, R.B.; Gettrust, J.F. Gas and gas hydrate distribution around seafloor seeps in Mississippi Canyon, Northern Gulf of Mexico, using multi-resolution seismic imagery. Mar. Pet. Geol. 2008, 9, 952-959.

24. Weimer, P.; Buffler, R.T. Structural geology and evolution of the Mississippi fan fold belt, deep Gulf of Mexico. AAPG Bull. 1992, 76, 225-251.

25. Goodwin, R.H.; Prior, D.B. Geometry and depositional sequences of the Mississippi Canyon, Gulf of Mexico. J. Sediment. Res. 1989, 59, 318-329.

26. Ellwood, B.B.; Balsam, W.L.; Roberts, H.H. Gulf of Mexico sediment sources and sediment transport trends from magnetic susceptibility measurements of surface samples. Mar. Geol. 2006, 230, 237-248.

27. Coffin, R.B.; Gardner, J.; Pohlman, J.; Downer, R.; Wood, W. Methane Hydrate Exploration, Atwater Valley, Texas-Louisiana Shelf: Geophysical And Geochemical Profiles; Naval Research Laboratory: Washington, DC, USA, 2006.

28. Reeburgh, W.S. An improved interstitial water sampler. Limnol. Oceanogr. 1967, 12, 163-165.

29. Plummer, R.E.; Pohlman, J.; Coffin, R.B. Compound-Specific Stable Carbon Isotope Analysis of Low-Concentration Complex Hydrocarbon Mixtures from Natural Gas Hydrate Systems; American Geophysical Union: Washington, DC, USA, 2005.

30. Osburn, C.L.; St-Jean, G. The use of wet chemical oxidation with high-amplification isotope ratio mass spectrometry (WCO-IRMS) to measure stable isotope values of dissolved organic carbon in seawater. Limnol. Oceanogr. Methods 2007, 5, 296-308.

31. Pohlman, J.W.; Knies, D.L.; Grabowski, K.S.; DeTurck, T.M.; Treacy, D.J.; Coffin, R.B. Sample distillation/graphitization system for carbon pool analysis by accelerator mass spectrometry (AMS). Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms 2000, 172, 428-433.

32. Grabowski, K.S.; Knies, D.L.; DeTurck, T.M.; Treacy, D.J.; Pohlman, J.W.; Coffin, R.B.; Hubler, G.K. A report on the Naval Research Laboratory AMS facility. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms 2000, 172, 34-39.

33. Tumey, S.T.; Grabowski, K.S.; Knies, D.L.; Mignerey, A.C. Data collection, filtering and analysis at the Naval Research Laboratory trace element accelerator mass spectrometry facility. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. Atoms 2004, 176, 428-433.

34. Stuiver, M.; Polach, H.A. Discussion: Reporting of 14C data. Radiocarbon 1977, 19, 355-363.

35. Stuiver, M. International agreements and the use of the new oxalic acid standard. Radiocarbon 1983, 25, 793-795.

36. Donahue, D.J.; Linick, T.W.; Jull, A.J.T. Isotope-ratio and background corrections for accelerator mass spectrometry radiocarbon measurements. Radiocarbon 1990, 32, 135-142.

37. Coffin, R.B.; Hamdan, L.J.; Smith, J.P.; Rose, P.S.; Plummer, R.E.; Yoza, B.; Pecher, I.; Montgomery, M.T. Contribution of vertical methane flux to shallow sediment carbon pools across porangahau ridge, New Zealand. Energies 2014, 7, 5332-5356, doi:10.3390/en7085332.

38. Macko, S.A.; Ostrom, N.E. Pollution studies using stable isotopes. In Stable Isotopes in Ecology and Environmental Science; Lajtha, K., Michener, R., Eds.; Blackwell Scientific Publications: Hoboken, NJ, USA, 1994; pp. 45-62.

39. Fry, B. Using stable isotopes to monitor watershed influences on aquatic trophodynamics. Can. J. Fish. Aquat. Sci. 1999, 56, 2167-2171.

40. Gordon, E.S.; Goni, M.A. Source and distribution of terrigenous organic matter delivered by the Atchafalaya River to sediments in the northern Gulf of Mexico. Geochem. Cosmochim. Acta 2003, 67, 2359-2375.

41. Wang, X.; Chen, R.F.; Gardner, G.B. Sources and transport of dissolved and particulate organic carbon in the Mississippi River estuary and adjacent coastal waters of the northern Gulf of Mexico. Mar. Chem. 2004, 89, 241-256.

42. Peterson, B.J.; Howarth, R.W.; Garritt, R.H. Multiple stable isotopes used to trace the flow of organic matter in estuarine food webs. Science 1985, 277, 1361-1363.

43. Chen, X.; Lohrenz, S.E.; Wiesenburg, D.A. Distribution and controlling mechanisms of primary production on the Louisianan-Texas continental shelf. J. Mar. Syst. 2000, 25, 179-207.

44. Burdige, D.J.; Martens, C.S. Biogeochemical cycling in an organic-rich coastal marine basin: The sedimentary cycling of dissolved, free amino acids. Geochim. Cosmochim. Acta 1990, 54, 3033-3052.

45. Lee, C. Controls on organic carbon preservation: The use of stratified water bodies to compare intrinsic rates of decomposition in oxic and anoxic systems. Geochim. Cosmochim. Acta 1992, 56, 3323-3335.

46. Canuel, E.A.; Martens, C.S. Seasonal variations in the sources and alteration of organic matter associated with recently-deposited sediments. Org. Geochem. 1993, 20, 563-577.

47. Gong, C.; Hollander, D.J. Differential contribution of bacteria to sedimentary organic matter in oxic and anoxic environments, Santa Monica Basin, California. Org. Geochem. 1997, 26, 545-563.

48. Orphan, V.J.; Hinrichs, K.U.; Ussler, W., III; Paull, C.K.; Taylor, L.T.; Sylva, S.P.; Hayes, J.M.; Delong, E.F. Comparative analysis of methane-oxidizing archaea and sulfate-reducing bacteria in anoxic marine sediments. Appl. Environ. Microbiol. 2001, 67, 1922-1934.

49. Reeburgh, W.S. Oceanic methane biogeochemistry. Chem. Rev. 2007, 107, 486-513.

50. Goevert, D.; Conrad, R. Effect of substrate concentration on carbon isotope fractionation during acetoclastic methanogenesis by Methanosarcina barkeri and M. acetivorans and in rice field soil. Appl. Environ. Microbiol. 2009, 75, 2605-2612.

51. Szynkiewicz, A.M.; Jedrysek, O.; Kurasiewicz, M. Carbon isotope effects during precipitation of barium carbonate: Implications for environmental studies. Environ. Chem. Lett. 2006, 4, 29-35.

52. Martens, C.S.; Chanton, J.P.; Paull, C.K. Fossil biogenic methane at the Florida escarpment. Geology 1991, 19, 851-854.

53. Sassen, R.; MacDonald, I.R. Hydrocarbons of experimental and natural gas hydrates, Gulf of Mexico continental slope. Org. Geochem. 1997, 26, 289-293.

54. Sassen, R.; Sweet, S.T.; DeFreitas, D.A.; Morelos, J.A.; Milkov, A.V. Gas hydrate and crude oil from the Mississippi Fan Foldbelt, downdip Gulf of Mexico Salt Basin: Significance to petroleum system. Org. Geochem. 2001, 32, 999-1008.

55. Pohlman, J.W.; Kaneko, M.; Heuer, V.B.; Coffin, R.B.; Whiticar, M. Methane soucrces and production in the north Cascadia Margin gas hydrate system. Earth Planet. Sci. Lett. 2009, 287, 504-512.

56. Whiticar, M.J. Carbon and hydrogen isotope systematics of bacterial formation and oxidation of methane. Chem. Geol. 1999, 16, 291-314.

57. Boetius, A.; Ravenschlag, K.; Schubert, C.J.; Rickert, D.; Widdel, F.; Gleske, A.; Amann, R.; Jorgensen, B.B.; Witte, U.; Pfannkuche, O. A marine microbial consortium apparently mediating anaerobic oxidation of methane. Nature 2000, 407, 623-626.

58. Orphan, V.J.; House, C.H.; Hinrichs, K.-U.; McKeegan, K.D.; DeLong, E.F. Methane-consuming archea revealed by directly coupled isotopic and phylogenetic analysis. Science 2001, 293, 484-487.

59. Alperin, M.J.; Blair, N.E.; Albert, D.B.; Hoehler, T.M.; Martens, C.S. Factors that control the stable carbon isotopic composition of methane produced in an anoxic marine sediment. Glob. Biogeochem. Cycles 1992, 6, 271-291.

60. Hoehler, T.M.; Alperin, M.J.; Albert, D.B.; Martens, C.S. Field and laboratory studies of methane oxidation in an anoxic marine sediment: Evidence for a methanogen-sulfate reducer consortium. Glob. Biogeochem. Cycles 1994, 8, 451-463.

61. Yoshinaga, M.Y.; Holler, T.; Goldhammer, T.; Wegener, G.; Pohlman, J.W.; Brunner, B.; Kuypers, M.M.M.; Hinrichs, K-U.; Elvert, M. Carbon isotope equilibration during sulphate-limited anaerobic oxidation of methane. Nat. Geosci. 2014, 7, doi:10.1038/NGE02069.

62. Borowski, W.S.; Paull, C.K.; Ussler, W., III. Global and local variations of interstitial sulfate gradients in the deep-water, continental margin sediments: Sensitivity to underlying methane and gas hydrates. Mar. Geol. 1999, 159, 131-154.

63. House, C.H.; Schopf, J.W.; Stetter, K.O. Carbon isotopic fractionation by Archaens and other thermophilic prokaryotes. Org. Geochem. 2003, 34, 345-356.

64. Zhang, C.L.; Fouke, B.W.; Bonheyo, G.T.; Peackock, A.D.; White, D.C.; Huang, Y.; Romanek, C.S. Lipid biomarkers and carbon-isotopes of modern travertine deposits (Yellowstone Natioinal Park, USA): Implications for biogeochemical dynamics in hot-spring systems. Geochem. Cosmochim. Acta 2004, 68, 3157-3169.

65. Quandt, L.; Gottschalk, G.; Ziegler, H.; Stichler, W. Isotope discrimination by photosynthetic bacteria. REMS Microbiol. Lett. 1977, 1, 125-128.

66. Morse, J.W.; Beazley, M.J. Organic matter in deepwater sediments of the Northern Gulf of Mexico and its relationship to the distribution of benthic organisms. Deep Sea Res. II 2008, 55, 2563-2571.

67. Burdige, D.J. Geochemistry of Marine Sediments; Princeton University Press: Princeton, NJ, USA, 2006.

68. Martens, C.S.; Klump, J.V. Biogeochemical cycling of an organic-rich coastal marine basin 4. An organic carbon budget for sediments dominated by sulfate reduction and methanogenesis. Geochem. Cosmochim. Acta 1984, 48, 1987-2004.

69. Hoehler, T.M.; Albert, D.B.; Alperin, M.J.; Martens, C.S. Acetogenesis from CO2 in an anoxic marine sediment. Limnol. Oceanogr. 1999, 44, 662-667.

© 2015 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/4.0/).

Supplementary Information

Table S1. Sediment and pore water geochemistry data from C2, C3, C7 and C8. All isotope values are in units of per mil (%o) as described in the manuscript.

Core ID Pore Water Depth Below Sea Floor (cm) DIC (mM) 513Cdic %OC 513Csoc %CaCO3 513Ctic DOC (mM) 513Cdoc A14Csoc Gas Sample Depth Below Sea Floor (cm) CH4 (mM) Sediment S13CCH4 Gas S13CCH4 C1/C2 Gas A14Cch4

C2 5 2.8 -9.4 0.91 -20.5 13.04 1.8 0.95 -22.9 -283 0 6.00 X 10-4 ND ND ND ND

35 3.3 -14.6 0.67 -21.1 11.55 1.6 0.91 -20.8 -561 30 9.34 x 10-4 ND ND ND ND

70 4.2 -19.9 0.52 -22.7 8.82 -2.2 0.82 -20.1 -665 65 8.13 x 10-4 ND ND ND ND

105 5.0 -23.3 0.70 -22.8 2.71 -4.9 0.85 -23.2 -652 100 8.73 x 10-4 ND ND ND ND

140 5.2 -25.8 0.76 -22.8 2.4 -5.8 1.30 -25.3 -650 135 1.01 x 10-3 ND ND ND ND

165 5.9 -28.6 0.72 -22.7 2.36 -5.8 1.16 -26.2 -647 160 1.10 x 10-3 ND ND ND ND

190 7.3 -31.9 0.71 -22.6 2.42 -6.1 1.23 -25.4 -647 185 1.00 x 10-3 ND ND ND ND

215 8.2 -32.7 0.73 -22.4 2.52 -5.9 0.88 -24.4 -657 210 1.31 x 10-3 ND ND ND ND

240 8.4 -33.8 0.75 -22.4 2.54 -6 1.13 -24.6 -643 235 1.45 x 10-3 ND ND ND ND

265 9.5 -36.2 0.72 -22.4 2.33 -5.9 1.47 -24.2 -651 260 1.38 x 10-3 -48.8 ND ND ND

290 12.0 -40.4 0.73 -22.4 2.47 -6.2 1.88 -24.0 -645 285 2.86 x 10-3 -62.3 ND ND ND

310 15.6 -43.7 0.64 -22.9 2.53 -4.3 2.24 -23.9 -679 305 1.13 x 10-2 -66.8 ND ND ND

330 19.5 -47.1 0.66 -22.9 2.68 -4 2.78 -25.3 -673 325 3.46 x 10-2 -81.1 ND ND ND

350 20.0 -41.5 0.78 -22.7 4.13 -5.3 3.59 -23.8 -624 345 1.39 x 10° -86.9 ND ND ND

375 16.9 -30.8 0.76 -22.6 4.47 -5.1 3.53 -23.7 -631 370 8.42 x 100 -82.4 ND 26555 ND

400 14.3 -22.4 0.77 -22.8 4.31 -5.3 3.60 -23.6 -640 395 8.91 x 100 -79.7 ND 26926 ND

420 13.0 -17.7 0.72 -22.9 4.61 -4.9 3.65 -24.0 -678 415 1.20 x 101 -78.3 ND 29430 ND

440 11.9 -13.4 0.71 -23.0 4.02 -5.1 3.59 -24.4 -659 435 1.25 x 101 -77.9 ND ND ND

465 10.3 -10.3 0.70 -22.9 4.07 -5.3 3.64 -25.1 -662 460 9.96 x 100 -76.5 ND 30494 ND

C3 5 25 13.6 15.6 -42.9 -16.6 ND ND ND ND ND ND ND ND ND ND ND ND ND ND 0 20 6.87 x 10-2 6.96 x 100 -75.4 -75.2 ND ND ND 6033 ND ND

Table S1. Cont.

Pore Water Core Depth Below ID Sea Floor (cm) DIC (mM) 813Cdic %OC 513Csoc %CaCO3 513Ctic DOC (mM) 513Cdoc A14Csoc Sample Depth Below Sea Floor (cm) CH4 (mM) Sediment 813Cch4 Gas 813Cch4 C1/C2 Gas A14Cch4

50 9.4 -1.8 ND ND ND ND ND ND ND 45 3.52 x 10° -75.3 ND 5244 ND

75 9.7 2.1 ND ND ND ND ND ND ND 70 7.87 x 100 -73.0 ND 7313 ND

100 9.1 3.0 ND ND ND ND ND ND ND 95 3.74 x 100 -73.1 ND 6122 ND

125 9.2 3.4 ND ND ND ND ND ND ND 120 3.36 x 100 -73.5 ND 6873 ND

155 10.5 2.7 ND ND ND ND ND ND ND 150 3.61 x 100 -73.1 -72.5 ± 0.1, n = 3 6021 -958

200 8.0 0.7 ND ND ND ND ND ND ND 195 3.80 x 100 -72.8 -72.0 6970 -963 ± 3, n = 5

225 8.9 -2.0 ND ND ND ND ND ND ND 220 4.28 x 100 -73.2 -72.4 6397 -960

250 9.2 3.3 ND ND ND ND ND ND ND 245 4.98 x 100 -73.6 -72.1 6686 -957

285 8.2 1.7 ND ND ND ND ND ND ND 280 3.69 x 100 -73.2 ND 7461 ND

310 8.3 0.8 ND ND ND ND ND ND ND 305 4.14 x 100 -73.4 ND 7033 ND

340 8.9 3.7 ND ND ND ND ND ND ND 335 4.37 x 100 -73.3 ND 7299 ND

365 6.7 0.1 ND ND ND ND ND ND ND 360 4.28 x 100 -73.7 ND 7380 ND

C7 5 5.8 3.2 0.63 -28.2 13.04 -12.6 3.83 -25.7 -890 0 8.34 x 100 -82.0 ND 9162 ND

34 12.1 -48.1 0.66 -26.4 10.36 -3.99 2.05 -27.6 -907 29 4.71 x 100 -73.0 ND 10602 ND

104 9.7 -25.4 0.64 -25.1 9.89 -1.93 3.01 -25.6 -939 99 3.29 x 100 -71.6 ND 7936 ND

144 8.3 -15.0 0.62 -25.0 9.93 -1.5 3.55 -25.4 -939 139 5.61 x 100 -69.3 ND 11753 ND

184 7.3 -0.7 0.58 -25.6 11.36 -2.4 3.67 -25.4 -950 179 5.85 x 100 -69.4 ND 8978 ND

224 7.1 3.0 0.57 -26.6 12.68 -1.7 3.77 -25.5 -955 219 4.82 x 100 -69.1 ND 7017 ND

322 7.3 3.1 0.57 -25.6 13.81 -2.4 3.76 -25.7 -950 317 7.75 x 100 -69.7 ND 8056 ND

372 7.4 2.5 0.58 -25.3 12.08 -1.1 3.69 -26.1 -947 367 6.06 x 100 -69.3 ND 5716 ND

422 7.1 4.5 0.63 -26.5 10.71 -3.9 3.79 -25.9 -899 417 4.89 x 100 -69.1 ND 7335 ND

Table S1. Cont.

Pore Water Core Depth Below ID Sea Floor (cm) DIC (mM) 813Cdic %OC 513Csoc %CaCO3 513Ctic DOC (mM) 513Cdoc A14Csoc Sample Depth Below Sea Floor (cm) CH4 (mM) Sediment 813Cch4 Gas 813Cch4 C1/C2 Gas A14Cch4

472 7.0 3.8 0.64 -25.2 10.32 -1.6 3.83 -26.0 -922 467 1.06 x 101 -69.4 ND 5210 ND

522 6.5 3.4 0.58 -25.3 11.71 -1.3 3.92 -26.7 -950 517 7.77 x 10° -68.5 ND 5836 ND

572 6.0 4.6 0.6 -26.3 10.69 -2.3 3.49 -25.2 -896 567 5.13 x 100 -68.5 -71.3 ± 0.1, n = 3 5967 ND

682 6.2 3.1 0.57 -25.5 11.01 -2.1 3.37 -24.8 -927 677 5.69 x 100 -69.0 ND 3106 ND

732 6.8 -1.5 0.58 -25.3 11.46 -1.6 3.73 -24.7 -929 727 5.29 x 100 -69.0 ND 6213 ND

782 6.1 3.0 0.58 -25.6 12.19 -1.5 3.67 -24.7 -951 777 6.79 x 100 -69.6 ND 5320 ND

842 5.6 2.5 0.63 -25.3 10.9 -1.2 3.41 -24.8 -918 837 6.57 x 100 -69.2 ND 3572 ND

C8 5 5.5 -38.3 ND ND ND ND ND ND ND 0 2.71 x 10-2 -73.3 ND ND ND

30 11.3 -54.7 ND ND ND ND ND ND ND 25 8.76 x 10-2 -78.9 ND ND ND

60 13.8 -44.3 ND ND ND ND ND ND ND 55 7.10 x 100 -85.2 ND 6777 ND

85 12.1 -37.9 ND ND ND ND ND ND ND 80 9.92 x 100 -81.8 ND 9108 ND

115 12.0 -29.3 ND ND ND ND ND ND ND 110 1.30 x 101 -80.1 ND 10562 ND

150 11.5 -18.1 ND ND ND ND ND ND ND 145 9.77 x 100 -78.6 ND 10021 ND

185 11.2 -7.4 ND ND ND ND ND ND ND 180 7.23 x 100 -76.9 ND 9691 ND

220 11.5 -4.7 ND ND ND ND ND ND ND 215 5.21 x 100 -76.0 -65.8 7994 -903

260 9.8 0.2 ND ND ND ND ND ND ND 255 5.32 x 100 -74.5 -70.9 6802 -905

310 10.8 2.0 ND ND ND ND ND ND ND 305 5.46 x 100 -75.3 ND 9054 ND

370 3.8 ND ND ND ND ND ND ND 365 4.96 x 100 -73.7 -72.3 8059 -901

420 10.1 5.8 ND ND ND ND ND ND ND 415 3.95 x 100 -73.6 ND 8036 ND

470 9.7 6.7 ND ND ND ND ND ND ND 465 2.85 x 100 -71.5 ND 9362 ND

510 8.6 6.1 ND ND ND ND ND ND ND 505 4.11 x 100 -72.5 ND 5450 ND

Copyright of Energies (19961073) is the property of MDPI Publishing and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.