Weak coupling between sulfate reduction and the anaerobic oxidation of methane in methane-rich seafloor sediments during ex situ incubation

We investigated coupling between sulfate reduction (SR) and anaerobic oxidation of methane (AOM) by quantifying pore water geochemical profiles, determining rates of microbial processes, and examining microbial community structure at two sites within Mississippi Canyon lease block 118 (MC118) in the Northern Gulf of Mexico. Sediments from the northwest seep contained high concentrations of methane while sediments from the southwest seep contained methane, gaseous n-alkanes and liquid hydrocarbons and had abundant surficial accumulations of gas hydrate. Volumetric (21.5 μmol cm⁻³ day⁻¹) and integrated (1429 mmol m⁻² day⁻¹) rates of SR at MC118 in ex situ incubations are the highest reported thus far for seafloor environments. AOM rates were small in comparison, with volumetric rates ranging from 0.1 to 12.6 nmol cm⁻³ day⁻¹. Diffusion cannot adequately supply the sulfate required to support these high SR rates so additional mechanisms, possibly biological sulfide oxidation and/or downward advection, play important roles in supplying sulfate at these sites. The microbial communities at MC118 included sulfate-reducing bacteria phylogenetically associated with Desulfobacterium anilini, which is capable of complex hydrocarbon degradation. Despite low AOM rates, the majority of archaea identified were phylogenetically related to previously described methane oxidizing archaea. To evaluate whether weak coupling between SR and AOM occurs in habitats lacking the complex hydrocarbon milieu present at MC118, we compiled available SR and AOM rates and found that the global median ratio of SR to AOM was 10.7:1 rather than the expected 1:1. The global median integrated AOM rate was used to refine global estimates for AOM rates at cold seeps; these new estimates are only 5% of the previous estimate.     

Different types of methane monooxygenases produce similar carbon and hydrogen isotope fractionation patterns during methane oxidation

We determined the stable carbon and hydrogen isotope fractionation factors for methane oxidation under oxic conditions using strains with known degradation pathways. The aerobic oxidation of methane can be initiated by two different forms of enzymes known as methane monooxygenases (MMO). The expression of these enzymes is type-specific and dependent upon the adjusted copper concentration in the medium (or environment). In this study, the expression of either the soluble MMO or the particulate MMO was supported by adjusting the copper concentrations in the growth medium. Taxonomically different aerobic methanotrophic strains, mainly belonging to the alpha- and gamma- classes of Proteobacteria, produced methane isotope enrichment factors (ε_bulk) ranging from −14.8 to −27.9‰ for carbon, and from −110.0 to −231.5‰ for hydrogen. The ratio of hydrogen versus carbon discrimination (Λ = (α_H⁻¹ − 1)/(α_C⁻¹ − 1) ≈ Δ(δ²H)/Δ(δ¹³C)) were similar for all tested cultures, and are also identical to values calculated from previously published enrichment factors for aerobic and anaerobic methane degradation. In contrast, Λ-values for the abiotic oxidation of methane with OH radicals (this process is considered as the main removal process for methane from the atmosphere) were significantly higher than the values derived from biotic oxidation. Due to the low variability of microbial methane isotope fractionation patterns, we propose that combined carbon and hydrogen isotope fractionation analyses can be used to monitor and assess the occurrence of microbial methane oxidation in marine or terrestrial environments. However, it is not possible to distinguish distinct aerobic or anaerobic methane-oxidation pathways by this approach.     

Carbon and hydrogen isotope fractionation associated with the aerobic microbial oxidation of methane, ethane, propane and butane

Carbon isotope fractionation factors associated with the aerobic consumption of methane (C1), ethane (C2), propane (C3), and n-butane (C4) were determined from incubations of marine sediment collected from the Coal Oil Point hydrocarbon seep field, located offshore Santa Barbara, CA. Hydrogen isotope fractionation factors for C1, C2 and C3 were determined concurrently. Fresh sediment samples from two seep areas were each slurried with sea water and treated with C1, C2, C3 or C4, or with mixtures of all four gases. Triplicate samples were incubated aerobically at 15 °C, and the stable isotope composition and headspace levels of C1-C4 were monitored over the course of the experiment. Oxidation was observed for all C1-C4 gases, with an apparent preference for C3 and C4 over C1 and C2 in the mixed-gas treatments. Fractionation factors were calculated using a Rayleigh model by comparing the δ¹³C and δD of the residual C1-C4 gases to their headspace levels. Carbon isotope fractionation factors (reported in ε or (α-1) × 1000 notation) were consistent between seep areas and were −26.5‰ ± 3.9 for C1, −8.0‰ ± 1.7 for C2, −4.8‰ ± 0.9 for C3 and −2.9‰ ± 0.9 for C4. Fractionation factors determined from mixed gas incubations were similar to those determined from individual gas incubations, though greater variability was observed during C1 consumption. In the case of C1 and C3 consumption, carbon isotope fractionation appears to decrease as substrate becomes limiting. Hydrogen isotope fractionation factors determined from the two seep areas differed for C1 oxidation but were similar for C2 and C3. Hydrogen isotope fractionation factors ranged from −319.9‰ to −156.4‰ for C1 incubations, and averaged −61.9‰ ± 8.3 for C2 incubations and −15.1‰ ± 1.9 for C3 incubations. The fractionation factors presented here may be applied to estimate the extent of C1-C4 oxidation in natural gas samples, and should prove useful in further studying the microbial oxidation of these compounds in the natural environment.     

Differential aerobic and anaerobic oxidation of hydrocarbon gases discharged at mud volcanoes in the Nile deep-sea fan

The present study investigates hydrocarbon oxidation processes at Isis and Amon mud volcanoes (MV’s), in the eastern Nile deep-sea fan. In the water column, molecular and carbon isotopic signatures of light hydrocarbons indicate that gases rapidly dissolve in seawater and are partially oxidized.In the upper sediments, anaerobic oxidation of the light hydrocarbons takes place, as clearly shown by their molecular and isotopic composition. These processes lead to the presence of a distinct Sulfate-Hydrocarbon Interface at 120-145 cm and 20-50 cm below the seafloor, for Isis and Amon MV’s, respectively. In contrast to processes occurring in the water column, a clear preferential oxidation of methane, propane and n-butane over ethane and i-butane is observed in the anoxic sediments. Furthermore, for the first time, fractionation factors have been determined for the anaerobic oxidation of propane and butane, being respectively −4.80‰ and −0.7‰ for δ¹³C, and −43.3‰ for δ²H of propane.     

Pathways and regulation of carbon, sulfur and energy transfer in marine sediments overlying methane gas hydrates on the Opouawe Bank (New Zealand)

This study combines sediment geochemical analysis, in situ benthic lander deployments and numerical modeling to quantify the biogeochemical cycles of carbon and sulfur and the associated rates of Gibbs energy production at a novel methane seep. The benthic ecosystem is dominated by a dense population of tube-building ampharetid polychaetes and conspicuous microbial mats were unusually absent. A 1D numerical reaction-transport model, which allows for the explicit growth of sulfide and methane oxidizing microorganisms, was tuned to the geochemical data using a fluid advection velocity of 14 cm yr⁻¹. The fluids provide a deep source of dissolved hydrogen sulfide and methane to the sediment with fluxes equal to 4.1 and 18.2 mmol m⁻² d⁻¹, respectively. Chemosynthetic biomass production in the subsurface sediment is estimated to be 2.8 mmol m⁻² d⁻¹ of C biomass. However, carbon and oxygen budgets indicate that chemosynthetic organisms living directly above or on the surface sediment have the potential to produce 12.3 mmol m⁻² d⁻¹ of C biomass. This autochthonous carbon source meets the ampharetid respiratory carbon demand of 23.2 mmol m⁻² d⁻¹ to within a factor of 2. By contrast, the contribution of photosynthetically-fixed carbon sources to ampharetid nutrition is minor (3.3 mmol m⁻² d⁻¹ of C). The data strongly suggest that mixing of labile autochthonous microbial detritus below the oxic layer sustains high measured rates of sulfate reduction in the uppermost 2 cm of the sulfidic sediment (100-200 nmol cm⁻³ d⁻¹). Similar rates have been reported in the literature for other seeps, from which we conclude that autochthonous organic matter is an important substrate for sulfate reducing bacteria in these sediment layers. A system-scale energy budget based on the chemosynthetic reaction pathways reveals that up to 8.3 kJ m⁻² d⁻¹ or 96 mW m⁻² of catabolic (Gibbs) energy is dissipated at the seep through oxidation reactions. The microorganisms mediating sulfide oxidation and anaerobic oxidation of methane (AOM) produce 95% and 2% of this energy flux, respectively. The low power output by AOM is due to strong bioenergetic constraints imposed on the reaction rate by the composition of the chemical environment. These constraints provide a high potential for dissolved methane efflux from the sediment (12.0 mmol m⁻² d⁻¹) and indicates a much lower efficiency of (dissolved) methane sequestration by AOM at seeps than considered previously. Nonetheless, AOM is able to consume a third of the ascending methane flux (5.9 mmol m⁻² d⁻¹ of CH₄) with a high efficiency of energy expenditure (35 mmol CH₄ kJ⁻¹). It is further proposed that bioenergetic limitation of AOM provides an explanation for the non-zero sulfate concentrations below the AOM zone observed here and in other active and passive margin sediments.     

The role of sulfur in chemical weathering and atmospheric CO2 fluxes: Evidence from major ions, δCDIC, and δSSO4 in rivers of the Canadian Cordillera

Water samples from the Fraser, Skeena and Nass River basins of the Canadian Cordillera were analyzed for dissolved major element concentrations (HCO₃⁻, SO₄²⁻, Cl⁻, Ca²⁺, Mg²⁺, K⁺, Na⁺), δ¹³C of dissolved inorganic carbon (δ¹³C_DIC), and δ³⁴S of dissolved sulfate (δ³⁴S_SO4) to quantify chemical weathering rates and exchanges of CO₂ between the atmosphere, hydrosphere, and lithosphere. Weathering rates of silicates and carbonates were determined from major element mass balance. Combining the major element mass balance with δ³⁴S_SO4 (−8.9 to 14.1‰_CDT) indicates sulfide oxidation (sulfuric acid production) and subsequent weathering of carbonate and to a lesser degree silicate minerals are important processes in the study area. We determine that on average, 81% of the riverine sulfate can be attributed to sulfide oxidation in the Cordilleran rivers, and that 25% of the total weathering cation flux can be attributed to carbonate and silicate dissolution by sulfuric acid. This result is validated by δ¹³C_DIC values (−9.8 to −3.7‰ _VPDB) which represents a mixture of DIC produced by the following weathering pathways: (i) carbonate dissolution by carbonic acid (−8.25‰) > (ii) silicate dissolution by carbonic acid (−17‰) ≈ (iii) carbonate dissolution by sulfuric acid derived from the oxidation of sulfides (coupled sulfide-carbonate weathering) (+0.5‰).δ³⁴S_SO4 is negatively correlated with δ¹³C_DIC in the Cordilleran rivers, which further supports the hypothesis that sulfuric acid produced by sulfide oxidation is primarily neutralized by carbonates, and that sulfide-carbonate weathering impacts the δ¹³C_DIC of rivers. The negative correlation between δ³⁴S_SO4 and δ¹³C_DIC is not observed in the Ottawa and St. Lawrence River basins. This suggests other factors such as landscape age (governed by tectonic uplift) and bedrock geology are important controls on regional sulfide oxidation rates, and therefore also on the magnitude of sulfide-carbonate weathering—i.e., it is more significant in tectonically active areas.Calculated DIC fluxes due to Ca and Mg silicate weathering by carbonic acid (38.3 × 10³ mol C · km⁻² · yr⁻¹) are similar in magnitude to DIC fluxes due to sulfide-carbonate weathering (18.5 × 10³ mol C · km⁻² · yr⁻¹). While Ca and Mg silicate weathering facilitates a transfer of atmospheric CO₂ to carbonate rocks, sulfide-carbonate weathering can liberate CO₂ from carbonate rocks to the atmosphere when sulfide oxidation exceeds sulfide deposition. This implies that in the Canadian Cordillera, sulfide-carbonate weathering can offset up to 48% of the current CO₂ drawdown by silicate weathering in the region.     

Characterization of calcium isotopes in natural and synthetic barite

The mineral barite (BaSO₄) accommodates calcium in its crystal lattice, providing an archive of Ca-isotopes in the highly stable sulfate mineral. Holocene marine (pelagic) barite samples from the major ocean basins are isotopically indistinguishable from each other (δ^44/40Ca = −2.01 ± 0.15‰) but are different from hydrothermal and cold seep barite samples (δ^44/40Ca = −4.13 to −2.72‰). Laboratory precipitated (synthetic) barite samples are more depleted in the heavy Ca-isotopes than pelagic marine barite and span a range of Ca-isotope compositions, Δ^44/40Ca = −3.42 to −2.40‰. Temperature, saturation state, , and aCa²⁺/aBa²⁺ each influence the fractionation of Ca-isotopes in synthetic barite; however, the fractionation in marine barite samples is not strongly related to any measured environmental parameter. First-principles lattice dynamical modeling predicts that at equilibrium Ca-substituted barite will have much lower ⁴⁴Ca/⁴⁰Ca than calcite, by −9‰ at 0 °C and −8‰ at 25 °C. Based on this model, none of the measured barite samples appear to be in isotopic equilibrium with their parent solutions, although as predicted they do record lower δ^44/40Ca values than seawater and calcite. Kinetic fractionation processes therefore most likely control the extent of isotopic fractionation exhibited in barite. Potential fractionation mechanisms include factors influencing Ca²⁺ substitution for Ba²⁺ in barite (e.g. ionic strength and trace element concentration of the solution, competing complexation reactions, precipitation or growth rate, temperature, pressure, and saturation state) as well as nucleation and crystal growth rates. These factors should be considered when investigating controls on isotopic fractionation of Ca²⁺ and other elements in inorganic and biogenic minerals.     

Can pelagic net heterotrophy account for carbon fluxes from eastern Canadian lakes?

Lakes worldwide are commonly oversaturated with CO₂, however the source of this CO₂ oversaturation is not well understood. To examine the magnitude of the C flux to the atmosphere and determine if an excess of respiration (R) over gross primary production (GPP) is sufficient to account for this C flux, metabolic parameters and stable isotopes of dissolved O₂ and C were measured in 23 Québec lakes. All of the lakes sampled were oversaturated with CO₂ over the sampling period, on average 221 ± 25%. However, little evidence was found to conclude that this CO₂ oversaturation was the result of an excess of pelagic R over GPP. In lakes Croche and à l’Ours, where CO₂ flux, R and GPP were measured weekly, the annual difference between pelagic GPP and R, or net primary production (NPP), was not sufficient to account for the size of the CO₂ flux to the atmosphere. In Lac Croche average annual NPP was 14.4 mg C m⁻² d⁻¹ while the average annual flux of CO₂ to the atmosphere was 34 mg C m⁻² d⁻¹. In Lac à l’Ours average annual NPP was −9.1 mg C m⁻² d⁻¹ while the average annual flux of CO₂ to the atmosphere was 55 mg C m⁻² d⁻¹. In all of the lakes sampled, O₂ saturation averaged 104.0 ± 1.7% during the ice-free season and the isotopic composition of dissolved O₂ (δ¹⁸O_DO) was 22.9 ± 0.3‰, lower than atmospheric values and indicative of net autotrophy. Carbon evasion was not a function of R, nor did the isotopic signature of dissolved CO₂ in the lakes present evidence of excess R over GPP. External inputs of C must therefore subsidize the lake to explain the continued CO₂ oversaturation. The isotopic composition of dissolved inorganic C (δ¹³C_DIC) indicates that the CO₂ oversaturation cannot be attributed to in situ aerobic respiration. δ¹³C_DIC reveals a source of excess C enriched in ¹³C, which may be accounted for by anaerobic sediment respiration or groundwater inputs followed by kinetic isotope fractionation during degassing under open system conditions.     

Coupling between Pbex and organic matter in sediments of a nutrient-enriched lake: An example from Lake Chenghai, China

Sediment cores were collected from deep-water areas of Lake Chenghai, China in June 1997. The vertical profile of ¹³⁷Cs activity gives reliable geochronological results. The results also indicate that sediment accumulation rates in deep-water areas of Lake Chenghai were relatively constant in recent decades, averaging 0.43 g cm^− 2 y^− 1, despite a variable organic carbon influx. ²¹⁰Pb_eq (= ²²⁶Ra) activity was relatively constant also, with an average value of 54.3 ± 3.2 Bq kg^− 1. Vertical profiles of ²¹⁰Pb_ex (= ²¹⁰Pb_total − ²²⁶Ra) decreased exponentially, resulting in somewhat lower sediment accumulation rates (0.3 g cm^− 2 y^− 1). These lower rates are likely less reliable, as the relatively large fluctuations in ²¹⁰Pb_ex activities correlate closely to the organic carbon (C_org) content of the sediments. For example, the vertical profile of ²¹⁰Pb_ex activity displays peaks at mass depths of 3.7-4.7 g cm^− 2 (10-12 cm) and 10-11 g cm^− 2(25-28 cm), similar to the maxima in the vertical profile of C_org. This phenomenon must be related to the delivery of particulate organic matter (POM) from the water to the sediments, or to watershed soil erosion. Since the mean atomic ratios of H_org / C_org and C_org / N_org in Lake Chenghai sediments are 5.5 and 7.0, respectively, indicating that POM was predominantly derived from the remains of authigenic algae, this eliminates watershed erosion rates as a primary control on lake sedimentation rates as resolved by ²¹⁰Pb_ex. Sedimentation fluxes (F(C_org)) of particulate organic carbon since 1970 varied between 60 to 160 g m^− 2 y^− 1, and appeared to closely influence variations in ²¹⁰Pb_ex concentrations. For example, sedimentation fluxes of ²¹⁰Pb_ex (F(²¹⁰Pb_ex)) showed maxima in the years 1972-1974 and 1986-1989, likely reflecting historical variations of lake biological productivity or carbon preservation.     

A cryptic sulfur cycle driven by iron in the methane zone of marine sediment (Aarhus Bay, Denmark)

Sulfate reduction and sulfur-iron geochemistry were studied in 5-6 m deep gravity cores of Holocene mud from Aarhus Bay (Denmark). A goal was to understand whether sulfate is generated by re-oxidation of sulfide throughout the sulfate and methane zones, which might explain the abundance of active sulfate reducers deep below the main sulfate zone. Sulfate penetrated down to 130 cm where methane started to build up and where the concentration of free sulfide peaked at 5.5 mM. Below this sulfate-methane transition, sulfide diffused downwards to a sulfidization front at 520 cm depth, below which dissolved iron, Fe²⁺, accumulated in the pore water. Sulfate reduction rates measured by ³⁵S-tracer incubations in the sulfate zone were high due to high concentrations of reactive organic matter. Within the sulfate-methane transition, sulfate reduction was distinctly stimulated by the anaerobic oxidation of methane. In the methane zone below, sulfate remained at positive “background” concentrations of <0.5 mM down to the sulfidization front. Sulfate reduction decreased steeply to rates which at 300-500 cm depth were 0.2-1 pmol SO₄²⁻ cm⁻³ d⁻¹, i.e., 4-5 orders of magnitude lower than rates measured near the sediment surface. The turn-over time of sulfate increased from 3 years at 12 cm depth to 100-1000 years down in the methane zone. Sulfate reduction in the methane zone accounted for only 0.1% of sulfate reduction in the entire sediment column and was apparently limited by the low pore water concentration of sulfate and the low availability of organic substrates. Amendment of the sediment with both sulfate and organic substrates immediately caused a 10- to 40-fold higher, “potential sulfate reduction” which showed that a physiologically intact community of sulfate reducing bacteria was present. The “background” sulfate concentration appears to be generated from the reaction of downwards diffusing sulfide with deeply buried Fe(III) species, such as poorly-reactive iron oxides or iron bound in reactive silicates. The oxidation of sulfide to sulfate in the sulfidic sediment may involve the formation of elemental sulfur and thiosulfate and their further disproportionation to sulfide and sulfate. The net reaction of sulfide and Fe(III) to form pyrite requires an additional oxidant, irrespective of the formation of sulfate. This could be CO₂ which is reduced with H₂ to methane. The methane subsequently diffuses upwards to become re-oxidized at the sulfate-methane transition and thereby removes excess reducing power and enables the formation of excess sulfate. We show here how the combination of these well-established sulfur-iron-carbon reactions may lead to the deep formation of sulfate and drive a cryptic sulfur cycle. The iron-rich post-glacial sediments underlying Holocene marine mud stimulate the strong sub-surface sulfide reoxidation observed in Aarhus Bay and are a result of the glacial to interglacial history of the Baltic Sea area. Yet, processes similar to the ones described here probably occur widespread in marine sediments, in particular along the ocean margins.     

Kinetics of microbial sulfate reduction in estuarine sediments

Kinetic parameters of microbial sulfate reduction in intertidal sediments from a freshwater, brackish and marine site of the Scheldt estuary (Belgium, the Netherlands) were determined. Sulfate reduction rates (SRR) were measured at 10, 21, and 30 °C, using both flow-through reactors containing intact sediment slices and conventional sediment slurries. At the three sites, and for all depth intervals studied (0-2, 2-4, 4-6 and 6-8 cm), the dependence of potential SRR on the sulfate concentration followed the Michaelis-Menten rate equation. Apparent sulfate half-saturation concentrations, K_m, measured in the flow-through reactor experiments were comparable at the freshwater and marine sites (0.1-0.3 mM), but somewhat higher at the brackish site (0.4-0.9 mM). Maximum potential SRR, R_max, in the 0-4 cm depth interval of the freshwater sediments were similar to those in the 0-6 cm interval of the marine sediments (10-46 nmol cm⁻³ h⁻¹ at 21 °C), despite much lower in situ sulfate availability and order-of-magnitude lower densities of sulfate-reducing bacteria (SRB), at the freshwater site. Values of R_max in the brackish sediments were lower (3.7-7.6 nmol cm⁻³ h⁻¹ at 21 °C), probably due to less labile organic matter, as inferred from higher C_org/N ratios. Inflow solutions supplemented with lactate enhanced potential SRR at all three sites. Slurry incubations systematically yielded higher R_max values than flow-through reactor experiments for the freshwater and brackish sediments, but similar values for the marine sediments. Transport limitation of potential SRR at the freshwater and brackish sites may be related to the lower sediment porosities and SRB densities compared to the marine site. Multiple rate controls, including sulfate availability, organic matter quality, temperature, and SRB abundance, modulate in situ sulfate-reducing activity along the estuarine salinity gradient.     

Silicate weathering in anoxic marine sediments

Two sediment cores retrieved at the northern slope of Sakhalin Island, Sea of Okhotsk, were analyzed for biogenic opal, organic carbon, carbonate, sulfur, major element concentrations, mineral contents, and dissolved substances including nutrients, sulfate, methane, major cations, humic substances, and total alkalinity. Down-core trends in mineral abundance suggest that plagioclase feldspars and other reactive silicate phases (olivine, pyroxene, volcanic ash) are transformed into smectite in the methanogenic sediment sections. The element ratios Na/Al, Mg/Al, and Ca/Al in the solid phase decrease with sediment depth indicating a loss of mobile cations with depth and producing a significant down-core increase in the chemical index of alteration. Pore waters separated from the sediment cores are highly enriched in dissolved magnesium, total alkalinity, humic substances, and boron. The high contents of dissolved organic carbon in the deeper methanogenic sediment sections (50-150 mg dm⁻³) may promote the dissolution of silicate phases through complexation of Al³⁺ and other structure-building cations. A non-steady state transport-reaction model was developed and applied to evaluate the down-core trends observed in the solid and dissolved phases. Dissolved Mg and total alkalinity were used to track the in-situ rates of marine silicate weathering since thermodynamic equilibrium calculations showed that these tracers are not affected by ion exchange processes with sediment surfaces. The modeling showed that silicate weathering is limited to the deeper methanogenic sediment section whereas reverse weathering was the dominant process in the overlying surface sediments. Depth-integrated rates of marine silicate weathering in methanogenic sediments derived from the model (81.4-99.2 mmol CO₂ m⁻² year⁻¹) are lower than the marine weathering rates calculated from the solid phase data (198-245 mmol CO₂ m⁻² year⁻¹) suggesting a decrease in marine weathering over time. The production of CO₂ through reverse weathering in surface sediments (4.22-15.0 mmol CO₂ m⁻² year⁻¹) is about one order of magnitude smaller than the weathering-induced CO₂ consumption in the underlying sediments. The evaluation of pore water data from other continental margin sites shows that silicate weathering is a common process in methanogenic sediments. The global rate of CO₂ consumption through marine silicate weathering estimated here as 5-20 Tmol CO₂ year⁻¹ is as high as the global rate of continental silicate weathering.     

Aeration to degas CO2, increase pH,and increase iron oxidation rates for efficient treatment of net alkaline mine drainage

Passive treatment systems for mine drainage use no energy other than gravity, but they require greater area than active treatment systems. Researchers are considering “hybrid” systems that have passive and active components for increased efficiency, especially where space limitations render passive-only technology ineffective. Flow-through reactor field experiments were conducted at two large net-alkaline anthracite mine discharges in central Pennsylvania. Assuming an Fe removal rate of 20 g m⁻² day⁻¹ and Fe loading from field data, 3.6 × 10³ and 3.0 × 10⁴ m² oxidation ponds would be required for the passive treatment of Site 21 and Packer 5 discharges, respectively. However, only a small area is available at each site. This paper demonstrates aeration to drive off CO₂, increase pH, and increase Fe(II) oxidation rates, enabling treatment within a small area compared to passive treatment methods, and introduces a geochemical model to accurately predict these rates as well as semi-passive treatment system sizing parameters. Both net-alkaline discharges were suboxic with a pH of ≈5.7, Fe(II) concentration of ≈16 mg L⁻¹, and low Mn and Al concentrations. Flow rates were ≈4000 L min⁻¹ at Site 21 and 15,000 L min⁻¹ at Packer 5. Three-h aeration experiments with flow rates scaled to a 14-L reactor resulted in pH increases from 5.7 to greater than 7, temperature increases from 12 to 22 °C, dissolved O₂ increases to saturation with respect to the atmosphere, and Fe(II) concentration decreases from 16 to <0.05 mg L⁻¹. A 17,000-L pilot-scale reactor at Site 21 produced similar results although aeration was not as complete as in the smaller reactor. Two non-aerated experiments at Site 21 with 13 and 25-h run times resulted in pH changes of ?0.2 and Fe(II) concentration decreases of less than 3 mg L⁻¹.     

Gas seepage from Tokamachi mud volcanoes, onshore Niigata Basin (Japan): Origin, post-genetic alterations and CH4-CO2 fluxes

Methane and CO₂ emissions from the two most active mud volcanoes in central Japan, Murono and Kamou (Tokamachi City, Niigata Basin), were measured in from both craters or vents (macro-seepage) and invisible exhalation from the soil (mini- and microseepage). Molecular and isotopic compositions of the released gases were also determined. Gas is thermogenic (δ¹³C_CH4 from −32.9‰ to −36.2‰), likely associated with oil, and enrichments of ¹³C in CO₂ (δ¹³C_CO2 up to +28.3‰) and propane (δ¹³C_C3H8 up to −8.6‰) suggest subsurface petroleum biodegradation. Gas source and post-genetic alteration processes did not change from 2004 to 2010. Methane flux ranged within the orders of magnitude of 10¹-10⁴ g m⁻² d⁻¹ in macro-seeps, and up to 446 g m⁻² d⁻¹ from diffuse seepage. Positive CH₄ fluxes from dry soil were widespread throughout the investigated areas. Total CH₄ emission from Murono and Kamou were estimated to be at least 20 and 3.7 ton a⁻¹, respectively, of which more than half was from invisible seepage surrounding the mud volcano vents. At the macro-seeps, CO₂ fluxes were directly proportional to CH₄ fluxes, and the volumetric ratios between CH₄ flux and CO₂ flux were similar to the compositional CH₄/CO₂ volume ratio. Macro-seep flux data, in addition to those of other 13 mud volcanoes, supported the hypothesis that molecular fractionation (increase of the “Bernard ratio” C₁/(C₂ + C₃)) is inversely proportional to gas migration fluxes. The CH₄ “emission factor” (total measured output divided by investigated seepage area) was similar to that derived in other mud volcanoes of the same size and activity. The updated global “emission-factor” data-set, now including 27 mud volcanoes from different countries, suggests that previous estimates of global CH₄ emission from mud volcanoes may be significantly underestimated.     

Natural seabed gas seeps as sources of atmospheric methane

Microbial and thermogenic methane migrates towards the seabed where some is utilised during microbially-mediated anaerobic oxidation. Excess methane escapes as gas seeps, which occur in a variety of geological contexts in every sea and ocean, from inter-tidal zones to deep ocean trenches. Some seeps are localised, gentle emanations; others are vigorous covering areas of >1 km²; the most prolific seeps reported (offshore Georgia) produce ~40 t CH₄ per year. Gas bubbles lose methane to the water as they rise, so deep water seeps are unlikely to contribute to the atmosphere. However, bubbles break the surface above some shallow water seeps. Estimates of the total methane contribution to the atmosphere are poorly constrained, largely because the data set is so small. 20 Tg yr^–1 is considered a realistic first approximation. This is a significant contribution to the global budget, particularly as methane from seeps is ¹⁴C-depleted. A seep measurement programme is urgently required.GEM     

Methane microseepage from different sectors of the Yakela condensed gas field in Tarim Basin,Xinjiang, China

Methane microseepage is the result of natural gas migration from subsurface hydrocarbon accumulations to the Earth’s surface, and it is quite common in commercial petroleum fields. While the role of microseepage as a pathfinder in petroleum exploration has been known for about 80 a, its significance as an atmospheric CH₄ source has only recently been studied, and flux data are currently available only in the USA and Europe. With the aim of increasing the global data-set and better understanding flux magnitudes and variabilities, microseepage is now being extensively studied in China. A static flux chamber method was recently applied to study microseepage emissions into the atmosphere in four different sectors of the Yakela condensed gas field in Tarim Basin, Xinjiang, China, and specifically in: (a) a faulted sector, across the Luntai fault systems; (b) an oil–water interface sector, at the northern margin of the field; (c) an oil–gas interface sector, in the middle of the field; (d) an external area, outside the northern gas field boundaries. The results show that positive CH₄ fluxes are pervasive in all sectors and therefore, only part of the CH₄ migrating from the deep oil–gas reservoirs is consumed in the soil by methanotrophic oxidation. The intensity of gas seepage seems to be controlled by subsurface geologic settings and lateral variabilities of natural gas pressure in the condensed gas field. The highest CH₄ fluxes, up to ∼14 mg m⁻² d⁻¹ (mean of 7.55 mg m⁻² d⁻¹) with higher spatial variability (standard deviation, σ: 2.58 mg m⁻² d⁻¹), occur in the Luntai fault sector. Merhane flux was lower in the oil–water area (mean of 0.53 mg m⁻² d⁻¹) and the external area (mean of 1.55 mg m⁻² d⁻¹), and at the intermediate level in the gas–oil sector (mean of 2.89 mg m⁻² d⁻¹). These values are consistent with microseepage data reported for petroleum basins in the USA and Europe. The build-up of methane concentration in the flux chambers is always coupled with an enrichment of ¹³C, from δ¹³C₁ of −46‰ to −42.5‰ (VPDB), which demonstrates that seeping methane is thermogenic, as that occurring in the deep Yakela reservoir. Daily variations of microseepage are very low, with minima in the afternoon, corresponding to higher soil temperature (and higher methanotrophic consumption), and maxima in the early morning (when soil temperatures are lowest). A preliminary and rough estimate of the total amount of CH₄ exhaled from the Yakela field is in the order of 10² tonnes a⁻¹.     

Methylmercury cycling in sediments on the continental shelf of southern New England

Exposure of humans to monomethylmercury (MMHg) occurs primarily through consumption of marine fish, yet there is limited understanding concerning the bioaccumulation and biogeochemistry of MMHg in the biologically productive coastal ocean. We examined the cycling of MMHg in sediments at three locations on the continental shelf of southern New England in September 2003. MMHg in surface sediments is related positively to inorganic Hg (Hg(II) = total Hg − MMHg), the geographical distribution of which is influenced by organic material. Organic matter also largely controls the sediment-water partitioning of Hg species and governs the availability of dissolved Hg(II) for methylation. Potential gross rates of MMHg production, assayed by experimental addition of ²⁰⁰Hg to intact sediment cores, are correlated inversely with the distribution coefficient (K_D) of Hg(II) and positively with the concentration of Hg(II), most probably as HgS⁰, in 0.2-μm filtered pore water of these low-sulfide deposits. Moreover, the efflux of dissolved MMHg to overlying water (i.e., net production at steady state) is correlated with the gross potential rate of MMHg production in surface sediments. These results suggest that the production and efflux of MMHg from coastal marine sediments is limited by Hg(II), loadings of which presumably are principally from atmospheric deposition to this region of the continental shelf. The estimated diffusive flux of MMHg from the shelf sediments averages 9 pmol m⁻² d⁻¹. This flux is comparable to that required to sustain the current rate of MMHg accumulation by marine fish, and may be enhanced by the efflux of MMHg from near-shore deposits contaminated more substantially with anthropogenic Hg. Hence, production and subsequent mobilization of MMHg from sediments in the coastal zone may be a major source of MMHg to the ocean and marine biota, including fishes consumed by humans.     

Transport and thermodynamics constrain belowground carbon turnover in a northern peatland

Rates of anaerobic respiration are of central importance for the long-term burial of carbon (C) in peatlands, which are a relevant sink in the global C cycle. To identify constraints on anaerobic peat decomposition, we determined detailed concentration depth profiles of decomposition end-products, i.e. methane (CH₄) and dissolved inorganic carbon (DIC), along with concentrations of relevant decomposition intermediates at an ombrotrophic Canadian peat bog. The magnitude of in situ net production rates of DIC and CH₄ was estimated by inverse pore-water modeling. Vertical transport in the peat was slow and dominated by diffusion leading to the buildup of DIC and CH₄ with depth (5500 μmol L⁻¹ DIC, 500 μmol L⁻¹ CH₄). Highest DIC and CH₄ production rates occurred close to the water table (decomposition constant k_d ∼ 10⁻³-10⁻⁴ a⁻¹) or in some distinct zones at depth (k_d ∼ 10⁻⁴ a⁻¹). Deeper into the peat, decomposition proceeded very slowly at about k_d = 10⁻⁷ a⁻¹. This pattern could be related to thermodynamic and transport constraints. The accumulation of metabolic end-products diminished in situ energy yields of acetoclastic methanogenesis to the threshold for microbially mediated processes (−20 to −25 kJ mol⁻¹ CH₄). The methanogenic precursor acetate also accumulated (150 μmol L⁻¹). In line with these findings, CH₄ was formed by hydrogenotrophic methanogenesis at Gibbs free energies of −35 to −40 kJ mol⁻¹ CH₄. This was indicated by an isotopic fractionation α_CO2-CH4 of 1.069-1.079. Fermentative degradation of acetate, propionate and butyrate attained Gibbs free energies close to 0 kJ mol⁻¹ substrate. Although methanogenesis was apparently limited by some other factor in some peat layers, transport and thermodynamic constraints likely impeded respiratory processes in the deeper peat. Constraints on the removal of DIC and CH₄ may thus slow decomposition and contribute to the sustained burial of C in northern peatlands.     

Contribution to characterisation of biochar to estimate the labile fraction of carbon

Different analytical techniques were used to find the most reliable and economic method for determining the labile fraction of C in biochar. Biochar was produced from pine, poplar and willow (PI, PO and WI, respectively) at two temperatures (400 and 550 °C) and characterised using spectroscopic techniques [solid state ¹³C nuclear magnetic resonance spectroscopy (NMR)], molecular markers [pyrolysis-gas chromatography-mass spectrometry (Py-GC-MS)], thermogravimetry (TG), elemental composition and wet oxidation (potassium permanganate and potassium dichromate). Short term incubation (110 h) of an A horizon from an Umbrisol amended with the biochar samples at two doses (7.5 and 15 t ha⁻¹) was also carried out to provide supplementary information on the influence of biochar-soil interaction on CO₂ evolution. Spectroscopic analysis demonstrated that the degree of biochar carbonisation was influenced by the type of feedstock and heating conditions and followed the order WI-400 < PI-400 ∼ WI-550 ∼ PO-400 < PO-550 < PI-550. The thermo-labile fraction of the biochar samples, estimated from TG, ranged between 21% and 49%. The fraction of total C oxidised with potassium permanganate (C_per/C_total) was <50 g kg⁻¹ in all cases, whereas potassium dichromate (C_dichro/C_total) oxidation efficiency ranged between 180 and 545 g kg⁻¹. For each type of feedstock, the highest values of either chemically or thermally degradable C corresponded to the biochar produced at low temperature. Results indicate that low cost methodologies, such as dichromate oxidation and TG, reflected the degree of biochar carbonisation, and could therefore be used to estimate the labile fraction of C in biochar.