Biogeochemical and historical drivers of microbial community composition and structure in sediments from Mercer Subglacial Lake, West Antarctica.

Ice streams that flow into Ross Ice Shelf are underlain by water-saturated sediments, a dynamic hydrological system, and subglacial lakes that intermittently discharge water downstream across grounding zones of West Antarctic Ice Sheet (WAIS). A 2.06 m composite sediment profile was recently recovered from Mercer Subglacial Lake, a 15 m deep water cavity beneath a 1087 m thick portion of the Mercer Ice Stream. We examined microbial abundances, used 16S rRNA gene amplicon sequencing to assess community structures, and characterized extracellular polymeric substances (EPS) associated with distinct lithologic units in the sediments. Bacterial and archaeal communities in the surficial sediments are more abundant and diverse, with significantly different compositions from those found deeper in the sediment column. The most abundant taxa are related to chemolithoautotrophs capable of oxidizing reduced nitrogen, sulfur, and iron compounds with oxygen, nitrate, or iron. Concentrations of dissolved methane and total organic carbon together with water content in the sediments are the strongest predictors of taxon and community composition. δ¹³C values for EPS (-25 to -30‰) are consistent with the primary source of carbon for biosynthesis originating from legacy marine organic matter. Comparison of communities to those in lake sediments under an adjacent ice stream (Whillans Subglacial Lake) and near its grounding zone provide seminal evidence for a subglacial metacommunity that is biogeochemically and evolutionarily linked through ice sheet dynamics and the transport of microbes, water, and sediments beneath WAIS.

Global fjords as transitory reservoirs of labile organic carbon modulated by organo-mineral interactions.

The global carbon cycle is strongly modulated by organic carbon (OC) sequestration and decomposition. Whereas OC sequestration is relatively well constrained, there are few quantitative estimates of its susceptibility to decomposition. Fjords are hot spots of sedimentation and OC sequestration in marine sediments. Here, we adopt fjords as model systems to investigate the reactivity of sedimentary OC by assessing the distribution of the activation energy required to break OC bonds. Our results reveal that OC in fjord sediments is more thermally labile than that in global sediments, which is governed by its unique provenance and organo-mineral interactions. We estimate that 61 ± 16% of the sedimentary OC in fjords is degradable. Once this OC is remobilized and remineralized during glacial maxima, the resulting metabolic CO2 could counterbalance up to 50 ppm of the atmospheric CO2 decrease during glacial times, making fjords critical actors in dampening glacial-interglacial climate fluctuations through negative carbon cycling loops.

Mapping spatial and temporal variation of seafloor organic matter Δ14C and δ13C in the Northern Gulf of Mexico following the Deepwater Horizon Oil Spill.

Following the Deepwater Horizon oil spill of 2010, large amounts of biodegraded oil (petrocarbon) sank to the seafloor. Our objectives were to 1) determine post-spill isotopic values as the sediments approached a new baseline and 2) track the recovery of affected sediments. Sediment organic carbon δ13C and Δ14C reached a post-spill baseline averaging -21.2 ± 0.9‰ (n = 129) and -220 ± 66‰ (n = 95). Spatial variations in seafloor organic carbon baseline isotopic values, 13C and 14C, were influenced by river discharge and hydrocarbon seepage, respectively. Inverse Distance Weighting of surface sediment Δ14C values away from seep sites showed a 50% decrease in the total mass of petrocarbon, from 2010 to 2014. We estimated a rate of loss of -2 × 109 g of petrocarbon-C/year, 2-11% of the degradation rates in surface slicks. Despite the observed recovery in sediments, lingering residual material in the surface sediments was evident seven years following the blowout.

Microbial Fe(III) reduction as a potential iron source from Holocene sediments beneath Larsen Ice Shelf.

Recent recession of the Larsen Ice Shelf C has revealed microbial alterations of illite in marine sediments, a process typically thought to occur during low-grade metamorphism. In situ breakdown of illite provides a previously-unobserved pathway for the release of dissolved Fe2+ to porewaters, thus enhancing clay-rich Antarctic sub-ice shelf sediments as an important source of Fe to Fe-limited surface Southern Ocean waters during ice shelf retreat after the Last Glacial Maximum. When sediments are underneath the ice shelf, Fe2+ from microbial reductive dissolution of illite/Fe-oxides may be exported to the water column. However, the initiation of an oxygenated, bioturbated sediment under receding ice shelves may oxidize Fe within surface porewaters, decreasing dissolved Fe2+ export to the ocean. Thus, we identify another ice-sheet feedback intimately tied to iron biogeochemistry during climate transitions. Further constraints on the geographical extent of this process will impact our understanding of iron-carbon feedbacks during major deglaciations.

Petrocarbon evolution: Ramped pyrolysis/oxidation and isotopic studies of contaminated oil sediments from the Deepwater Horizon oil spill in the Gulf of Mexico.

Hydrocarbons released during the Deepwater Horizon (DWH) oil spill weathered due to exposure to oxygen, light, and microbes. During weathering, the hydrocarbons' reactivity and lability was altered, but it remained identifiable as "petrocarbon" due to its retention of the distinctive isotope signatures (14C and 13C) of petroleum. Relative to the initial estimates of the quantity of oil-residue deposited in Gulf sediments based on 2010-2011 data, the overall coverage and quantity of the fossil carbon on the seafloor has been attenuated. To analyze recovery of oil contaminated deep-sea sediments in the northern Gulf of Mexico we tracked the carbon isotopic composition (13C and 14C, radiocarbon) of bulk sedimentary organic carbon through time at 4 sites. Using ramped pyrolysis/oxidation, we determined the thermochemical stability of sediment organic matter at 5 sites, two of these in time series. There were clear differences between crude oil (which decomposed at a lower temperature during ramped oxidation), natural hydrocarbon seep sediment (decomposing at a higher temperature; Δ14C = -912‰) and our control site (decomposing at a moderate temperature; Δ14C = -189‰), in both the stability (ability to withstand ramped temperatures in oxic conditions) and carbon isotope signatures. We observed recovery toward our control site bulk Δ14C composition at sites further from the wellhead in ~4 years, whereas sites in closer proximity had longer recovery times. The thermographs also indicated temporal changes in the composition of contaminated sediment, with shifts towards higher temperature CO2 evolution over time at a site near the wellhead, and loss of higher temperature CO2 peaks at a more distant site.

Compositional and beam-size-dependent effects on pressure baseline in clumped isotope mass spectrometry.

RATIONALE: The analysis of carbonate samples for the application of clumped isotopes to paleoclimate reconstruction necessitates smaller beam intensities. However, there is a relationship between beam intensity and pressure-dependent baseline (PBL), and therefore between beam intensity and the correction for PBL. Here we explain the relationship between PBL and beam intensity to develop a better correction protocol and an improved understanding of clumped isotope mass spectrometry. METHODS: We describe a beam size experiment using our Isoprime isotope ratio mass spectrometer in which samples of the carbonate standard IAEA-C1 were analyzed at 30, 50, and 70 nA to establish an optimal protocol and a new method to correct for PBL using the theoretical constraint of invariable Δ47 over a range of δ47 (bulk isotope composition) values. We also explore the effects of both over- and under-correction of PBL on equilibrated and heated gas samples to understand the effect of mis-correction of PBL. RESULTS: The results of our beam size experiments showed that a direct measurement of the baseline consistently introduced variability to measurements of the Δ47 of heated gases, equilibrated gases, and carbonate standards. These results necessitated a new protocol to account for PBL in our system. Our new approach flattens the reference frame line slope to 0 and, importantly, reduces the variability of data points about the heated gas line. We also describe, for the first time, an empirically derived description of the compositional effect of PBL. CONCLUSIONS: A seemingly small change in our isotope ratio mass spectrometer resulted in a better understanding of PBL, for which we have developed an empirically based correction protocol to apply. Our new protocol has the potential to reduce analytical time for laboratories measuring PBL, and supports the need for carbonate mineral-based clumped isotope standards.

A centuries-long delay between a paleo-ice-shelf collapse and grounding-line retreat in the Whales Deep Basin, eastern Ross Sea, Antarctica.

Recent thinning and loss of Antarctic ice shelves has been followed by near synchronous acceleration of ice flow that may eventually lead to sustained deflation and significant contraction in the extent of grounded and floating ice. Here, we present radiocarbon dates from foraminifera that constrain the time elapsed between a previously described paleo-ice-shelf collapse and the subsequent major grounding-line retreat in the Whales Deep Basin (WDB) of eastern Ross Sea. The dates indicate that West Antarctic Ice Sheet (WAIS) grounding-line retreat from the continental shelf edge was underway prior to 14.7 ± 0.4 cal kyr BP. A paleo-ice-shelf collapse occurred at 12.3 ± 0.2 cal kyr BP. The grounding position was maintained on the outer-continental shelf until at least 11.5 ± 0.3 cal kyr BP before experiencing a 200-km retreat. Given the age uncertainties, the major grounding-line retreat lagged ice-shelf collapse by at least two centuries and by as much as fourteen centuries. In the WDB, the centuries-long delay in the retreat of grounded ice was partly due to rapid aggradational stacking of an unusually large volume of grounding-zone-wedge sediment as ice-stream discharge accelerated following ice-shelf collapse. This new deglacial reconstruction shows that ongoing changes to ice shelves may trigger complex dynamics whose consequences are realized only after a significant lag.

Hydrocarbon degradation and response of seafloor sediment bacterial community in the northern Gulf of Mexico to light Louisiana sweet crude oil.

The Deepwater Horizon (DWH) blowout resulted in the deposition to the seafloor of up to 4.9% of 200 million gallons of oil released into the Gulf of Mexico. The petroleum hydrocarbon concentrations near the wellhead were high immediately after the spill, but returned to background levels a few years after the spill. Microbial communities in the seafloor are thought to be responsible for the degradation of hydrocarbons, however, our knowledge is primarily based upon gene diversity surveys and hydrocarbon concentration in field sediment samples. Here, we investigated the oil degradation potential and changes in bacterial community by amending seafloor sediment collected near the DWH site with crude oil and both oil and Corexit dispersant. Polycyclic aromatic hydrocarbons were rapidly degraded during the first 30 days of incubation, while alkanes were degraded more slowly. With the degradation of hydrocarbons, the relative abundances of Colwelliaceae, Alteromonadaceae, Methylococales, Alcanivorax, Bacteriovorax, and Phaeobacter increased remarkably. However, the abundances of oil-degrading bacteria changed with oil chemistry. Colwelliaceae decreased with increasing oil degradation, whereas Alcanivorax and Methylococcales increased considerably. We assembled seven genomes from the metagenome, including ones belonging to Colwellia, Alteromonadaceae, Rhodobacteraceae, the newly reported genus Woeseia, and candidate phylum NC10, all of which possess a repertoire of genes for hydrocarbon degradation. Moreover, genes related to hydrocarbon degradation were highly enriched in the oiled treatment, suggesting that the hydrocarbons were biodegraded, and that the indigenous microflora have a remarkable potential for the natural attenuation of spilled oil in the deep-sea surface sediment.

Petroleum hydrocarbon persistence following the Deepwater Horizon oil spill as a function of shoreline energy.

An important aspect of oil spill science is understanding how the compounds within spilled oil, especially toxic components, change with weathering. In this study we follow the evolution of petroleum hydrocarbons, including n-alkanes, polycyclic aromatic hydrocarbons (PAHs) and alkylated PAHs, on a Louisiana beach and salt marsh for three years following the Deepwater Horizon spill. Relative to source oil, we report overall depletion of low molecular weight n-alkanes and PAHs in all locations with time. The magnitude of depletion, however, depends on the sampling location, whereby sites with highest wave energy have highest compound depletion. Oiled sediment from an enclosed bay shows high enrichment of high molecular weight PAHs relative to 17α(H),21ß(H)-hopane, suggesting the contribution from sources other than the Deepwater Horizon spill, such as fossil fuel burning. This insight into hydrocarbon persistence as a function of hydrography and hydrocarbon source can inform policy and response for future spills.

Sediment Core Extrusion Method at Millimeter Resolution Using a Calibrated, Threaded-rod.

Aquatic sediment core subsampling is commonly performed at cm or half-cm resolution. Depending on the sedimentation rate and depositional environment, this resolution provides records at the annual to decadal scale, at best. An extrusion method, using a calibrated, threaded-rod is presented here, which allows for millimeter-scale subsampling of aquatic sediment cores of varying diameters. Millimeter scale subsampling allows for sub-annual to monthly analysis of the sedimentary record, an order of magnitude higher than typical sampling schemes. The extruder consists of a 2 m aluminum frame and base, two core tube clamps, a threaded-rod, and a 1 m piston. The sediment core is placed above the piston and clamped to the frame. An acrylic sampling collar is affixed to the upper 5 cm of the core tube and provides a platform from which to extract sub-samples. The piston is rotated around the threaded-rod at calibrated intervals and gently pushes the sediment out the top of the core tube. The sediment is then isolated into the sampling collar and placed into an appropriate sampling vessel (e.g., jar or bag). This method also preserves the unconsolidated samples (i.e., high pore water content) at the surface, providing a consistent sampling volume. This mm scale extrusion method was applied to cores collected in the northern Gulf of Mexico following the Deepwater Horizon submarine oil release. Evidence suggests that it is necessary to sample at the mm scale to fully characterize events that occur on the monthly time-scale for continental slope sediments.

Distributions and accumulation rates of polycyclic aromatic hydrocarbons in the northern Gulf of Mexico sediments.

Sediment samples collected from shelf, slope and interior basin of the northern Gulf of Mexico during 2011-2013, 1-3 years after the Deepwater Horizon (DWH) oil spill, were utilized to characterize PAH pollution history, in this region. Results indicate that the concentrations of surface ΣPAH43 and their accumulation rates vary between 44 and 160 ng g(-1) and 6-55 ng cm(-2) y(-1), respectively. ΣPAH43 concentration profiles, accumulation rates and Δ(14)C values are significantly altered only for the sediments in the immediate vicinity of the DWH wellhead. This shows that the impact of DWH oil input on deep-sea sediments was generally limited to the area close to the spill site. Further, the PAHs source diagnostic analyses suggest a noticeable change in PAHs composition from higher to lower molecular weight dominance which reflects a change in source of PAHs in the past three years, back to the background composition. Results indicate low to moderate levels of PAH pollution in this region at present, which are unlikely to cause adverse effects on benthic communities.

Using natural abundance radiocarbon to trace the flux of petrocarbon to the seafloor following the Deepwater Horizon oil spill.

In 2010, the Deepwater Horizon accident released 4.6­6.0 × 10(11) grams or 4.1 to 4.6 million barrels of fossil petroleum derived carbon (petrocarbon) as oil into the Gulf of Mexico. Natural abundance radiocarbon measurements on surface sediment organic matter in a 2.4 × 10(10) m(2) deep-water region surrounding the spill site indicate the deposition of a fossil-carbon containing layer that included 1.6 to 2.6 × 10(10) grams of oil-derived carbon. This quantity represents between 0.5 to 9.1% of the released petrocarbon, with a best estimate of 3.0­4.9%. These values may be lower limit estimates of the fraction of the oil that was deposited on the seafloor because they focus on a limited mostly deep-water area of the Gulf, include a conservative estimate of thickness of the depositional layer, and use an average background or prespill radiocarbon value for sedimentary organic carbon that produces a conservative value. A similar approach using hopane tracer estimated that 4­31% of 2 million barrels of oil that stayed in the deep sea settled on the bottom. Converting that to a percentage of the total oil that entered into the environment (to which we normalized our estimate) converts this range to 1.8 to 14.4%. Although extrapolated over a larger area, our independent estimate produced similar values.

Blank corrections for ramped pyrolysis radiocarbon dating of sedimentary and soil organic carbon.

Ramped pyrolysis (RP) targets distinct components of soil and sedimentary organic carbon based on their thermochemical stabilities and allows the determination of the full spectrum of radiocarbon ((14)C) ages present in a soil or sediment sample. Extending the method into realms where more precise ages are needed or where smaller samples need to be measured involves better understanding of the blank contamination associated with the method. Here, we use a compiled data set of RP measurements of samples of known age to evaluate the mass of the carbon blank and its associated (14)C signature, and to assess the performance of the RP system. We estimate blank contamination during RP using two methods, the modern-dead and the isotope dilution method. Our results indicate that during one complete RP run samples are contaminated by 8.8 ± 4.4 µg (time-dependent) of modern carbon (MC, fM ∼ 1) and 4.1 ± 5.5 µg (time-independent) of dead carbon (DC, fM ∼ 0). We find that the modern-dead method provides more accurate estimates of uncertainties in blank contamination; therefore, the isotope dilution method should be used with caution when the variability of the blank is high. Additionally, we show that RP can routinely produce accurate (14)C dates with precisions ∼100 (14)C years for materials deposited in the last 10,000 years and ∼300 (14)C years for carbon with (14)C ages of up to 20,000 years.

Varying relative degradation rates of oil in different forms and environments revealed by ramped pyrolysis.

Degradation of oil contamination yields stabilized products by removing and transforming reactive and volatile compounds. In contaminated coastal environments, the processes of degradation are influenced by shoreline energy, which increases the surface area of the oil as well as exchange between oil, water, sediments, microbes, oxygen, and nutrients. Here, a ramped pyrolysis carbon isotope technique is employed to investigate thermochemical and isotopic changes in organic material from coastal environments contaminated with oil from the 2010 BP Deepwater Horizon oil spill. Oiled beach sediment, tar ball, and marsh samples were collected from a barrier island and a brackish marsh in southeast Louisiana over a period of 881 days. Stable carbon ((13)C) and radiocarbon ((14)C) isotopic data demonstrate a predominance of oil-derived carbon in the organic material. Ramped pyrolysis profiles indicate that the organic material was transformed into more stable forms. Our data indicate relative rates of stabilization in the following order, from fastest to slowest: high energy beach sediments > low energy beach sediments > marsh > tar balls. Oil was transformed most rapidly where shoreline energy and the rates of oil dispersion and exchange with water, sediments, microbes, oxygen, and nutrients were greatest. Still, isotope data reveal persistence of oil.

Measurement of multiply substituted isotopologues ('clumped isotopes') of CO2 using a 5 kV compact isotope ratio mass spectrometer: performance, reference frame, and carbonate paleothermometry.

RATIONALE: The burgeoning field of 'clumped isotope' paleothermometry, which has broad applications in geosciences, depends almost entirely on measurements made on two types of mass spectrometer with the same ion source design. Demonstration that these measurements can be carried out on a retrofitted mass spectrometer with a different type of ion source is important to the growing community of geoscientists considering employing this technique. METHODS: We summarize the sample preparation techniques, mass spectrometry, and data processing involved in establishing an absolute reference frame for the measurement of Δ47 in CO2 from carbonate minerals using a prototype 5 kV mass spectrometer. The prototype differs from previously marketed mass spectrometers only in the presence of six extra Faraday cups designed specifically to measure ion beams of masses 44-49. We employ pressure baseline corrections on the small mass 47 signal. RESULTS: The system was fully capable of establishing absolute reference frames for clumped isotope carbonate paleothermometry, despite alternative use of the system for other types of measurements (carbonate δ(18)O and δ(13)C values, water δ(18)O and δ(2)H values via equilibration, and δ(13) C and δ(15)N values of sedimentary organic material) in different reference frames. Repeated measurements of CO2 from carbonate standards yielded external precisions comparable with those from other laboratories and adequate for paleothermometry to a precision of near or below ±5 °C. The accuracy was demonstrated by reasonable temperatures being obtained for different types of repeatedly measured carbonates. Measurements of the pressure baseline with every sample increased analysis times by 20 min, but decreased the stabilization time of the reference frame by days to weeks, while improving precision. CONCLUSIONS: Simple modifications to existing instrumentation will allow more laboratories to perform clumped isotope measurements. The prototype mass spectrometer used herein provides precision comparable with all published data. We show that, where negatively charged species are associated with large ion beams, making pressure baseline corrections with every measurement is increasingly important no matter what instrument is used.

Accelerator mass spectrometry 14C determination in CO2 produced from laser decomposition of aragonite.

The determination of (14)C in aragonite (CaCO(3)) decomposed thermally to CO(2) using an yttrium-aluminum-garnet doped neodymium laser is reported. Laser decomposition accelerator mass spectrometry (LD-AMS) measurements reproduce AMS determinations of (14)C from the conventional reaction of aragonite with concentrated phosphoric acid. The lack of significant differences between these sets of measurements indicates that LD-AMS radiocarbon dating can overcome the significant fractionation that has been observed during stable isotope (C and O) laser decomposition analysis of different carbonate minerals. The laser regularly converted nearly 30% of material removed into CO(2) despite it being optimized for ablation, where laser energy breaks material apart rather than chemically altering it. These results illustrate promise for using laser decomposition on the front-end of AMS systems that directly measure CO(2) gas. The feasibility of such measurements depends on (1) the improvement of material removal and/or CO(2) generation efficiency of the laser decomposition system and (2) the ionization efficiency of AMS systems measuring continuously flowing CO(2).