Email updates

Keep up to date with the latest news and content from Geochemical Transactions and Chemistry Central.

Open Access Highly Accessed Research article

Geochemically induced shifts in catabolic energy yields explain past ecological changes of diffuse vents in the East Pacific Rise 9°50'N area

Michael Hentscher* and Wolfgang Bach

Author Affiliations

Department of Geosciences, University of Bremen, Klagenfurter Straße, 28359 Bremen, Germany

For all author emails, please log on.

Geochemical Transactions 2012, 13:2  doi:10.1186/1467-4866-13-2

The electronic version of this article is the complete one and can be found online at: http://www.geochemicaltransactions.com/content/13/1/2


Received:28 June 2011
Accepted:27 January 2012
Published:27 January 2012

© 2012 Hentscher and Bach; licensee Chemistry Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The East Pacific Rise (EPR) at 9°50'N hosts a hydrothermal vent field (Bio9) where the change in fluid chemistry is believed to have caused the demise of a tubeworm colony. We test this hypothesis and expand on it by providing a thermodynamic perspective in calculating free energies for a range of catabolic reactions from published compositional data. The energy calculations show that there was excess H2S in the fluids and that oxygen was the limiting reactant from 1991 to 1997. Energy levels are generally high, although they declined in that time span. In 1997, sulfide availability decreased substantially and H2S was the limiting reactant. Energy availability dropped by a factor of 10 to 20 from what it had been between 1991 and 1995. The perishing of the tubeworm colonies began in 1995 and coincided with the timing of energy decrease for sulfide oxidizers. In the same time interval, energy availability for iron oxidizers increased by a factor of 6 to 8, and, in 1997, there was 25 times more energy per transferred electron in iron oxidation than in sulfide oxidation. This change coincides with a massive spread of red staining (putative colonization by Fe-oxidizing bacteria) between 1995 and 1997.

For a different cluster of vents from the EPR 9°50'N area (Tube Worm Pillar), thermodynamic modeling is used to examine changes in subseafloor catabolic metabolism between 1992 and 2000. These reactions are deduced from deviations in diffuse fluid compositions from conservative behavior of redox-sensitive species. We show that hydrogen is significantly reduced relative to values expected from conservative mixing. While H2 concentrations of the hydrothermal endmember fluids were constant between 1992 and 1995, the affinities for hydrogenotrophic reactions in the diffuse fluids decreased by a factor of 15 and then remained constant between 1995 and 2000. Previously, these fluids have been shown to support subseafloor methanogenesis. Our calculation results corroborate these findings and indicate that the 1992-1995 period was one of active growth of hydrogenotrophic communities, while the system was more or less at steady state between 1995 and 2000.

Introduction

Microorganisms have the ability to gain energy for their metabolism by promoting a large range of redox reactions. Well-known energy sources are for example aerobic oxidation of methane or hydrogen sulfide, methanogenesis, fermentation, and sulfate reduction under anaerobic conditions [1]. In habitats like hydrothermal systems or mines, lacking sunlight and organic carbon sources, the primary production depends on electron donors that are released by water-rock reactions. High-temperature (> 400°C) processes of water-rock interaction determine the composition of seawater-derived hydrothermal fluids that are equilibrated with rocks at depths as much as several kilometers (Figure 1). Upon upwelling, these fluids cool (conductively and/or adiabatically) and mix with cold seawater to varying extents. High temperature fluids, venting focused via black smoker chimneys, often show little evidence for subseafloor mixing and are typically used as "hydrothermal endmember" compositions. Commonly, sites of diffuse venting are developed around the black smokers, and the temperature-composition relations of the fluids issuing through the seafloor there indicate that the diffuse fluids formed by subseafloor cooling and mixing of hot hydrothermal fluids with cold seawater. The seafloor underneath these diffuse vent sites is a particularly favorable environment for a variety of chemosynthetic microorganisms in terms of suitable temperature and large energy availability (Figure 1). The composition of the upwelling hydrothermal fluids in these diffuse vent sites imposes a major control on the metabolic diversity in the colonizing ecosystem. Because of this tight relation between vent ecosystem and fluid compositions, chemical changes in the fluid may directly influence the ecosystem.

thumbnailFigure 1. Sketch of idealized fluid flow within a hydrothermal system and potential catabolic reactions in different environments (chimney wall, plume, recharge zone, and subseafloor mixing zone). Upwelling hot, reducing hydrothermal fluids mix with entrained cold, oxygenated seawater in subseafloor mixing zones. For these zones, the affinities of the catabolic reactions provided in the inset are examined in this paper.

Thermodynamic calculations based on geochemical compositions of waters in these habitats provide insights into the energy availability and can determine possible reactions that can support primary production in these systems [2-5]. Tight relations between the availability of geochemical energy and microbial processes have been demonstrated for a variety of submarine hydrothermal environments, including chimney walls, diffuse fluids, and vent mussels [6-10].

In this study we use geochemical data from two hydrothermally active vents in the East Pacific Rise 9°50'N area to show that thermodynamic modeling can help interpret the microbial metabolism in such systems. For the first area, our calculations provide clues to the biological evolution of a vent site influenced by dynamic changes in fluid chemistry and, consequently, catabolic energy. The other case shows that microbial processes in the subseafloor may be deciphered by determining and comparing free energies of reactions for catabolic reactions of hypothetical fluids derived from conservative mixing of seawater and hydrothermal fluid with the diffuse fluid actually sampled.

Methods

Calculation of affinity

Free energy for catabolic reactions is available only if the system is out of geochemical equilibrium. Disequilibrium prevails when the properties of the system change at rates faster than the rates at which the thermodynamically favored reactions proceed. The abiotic rates of many redox reactions are sluggish, in particular at temperatures conducive of life (< 120°C)[5]. Microbes use enzymes to catalyze these redox reactions and harness the free energy by controlling electron transfers and converting a sizable fraction of the catabolic energy in ATP production for their anabolic metabolism [11]. The maximum quantity of free energy that microorganisms can catabolize (ΔrG) is given by the Gibbs energy at a reference state (ΔrG° = -RTlnKr) representing the intensive parameters (P, T) and an extensive term (RTlnQ) that captures the compositions of the vent solutions (equation 1):

<a onClick="popup('http://www.geochemicaltransactions.com/content/13/1/2/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.geochemicaltransactions.com/content/13/1/2/mathml/M1">View MathML</a>

where R is the universal gas constant and T the temperature in Kelvin. Kr is the calculated equilibrium for the temperature and pressure of interest, and Qr expresses the activities of species participating a specific reaction. Qr is evaluated through equation (2):

<a onClick="popup('http://www.geochemicaltransactions.com/content/13/1/2/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.geochemicaltransactions.com/content/13/1/2/mathml/M2">View MathML</a>

where ai represents the activity of the chemical species in the reaction, vir denotes the stoichiometric coefficient for the ith chemical species in the reaction, which is positive for products and negative for reactants. If ΔrG for a reaction is negative, then the reaction should proceed from left to right; if it is greater than zero, the reaction will proceed in the opposite direction. By convention a negative sign indicates that the reaction should take place spontaneously and energy can be gained by microbes catalyzing this reaction.

Commonly, affinity is used instead of ΔrG for a reaction [12]. Affinities express the change of the Gibbs energy with reaction progress (ξ) (equation 3):

<a onClick="popup('http://www.geochemicaltransactions.com/content/13/1/2/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.geochemicaltransactions.com/content/13/1/2/mathml/M3">View MathML</a>

It follows that the reaction is favorable if the affinity is positive. Combining equations 1 and 3, the affinity can be evaluated through equation (4):

<a onClick="popup('http://www.geochemicaltransactions.com/content/13/1/2/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.geochemicaltransactions.com/content/13/1/2/mathml/M4">View MathML</a>

This relation demonstrates that, if Kr > Qr, then Ar > 0 and the reaction may proceed while free energy is released [12].

Two types of computations were employed in this study: (1) calculation of concentrations and activities of dissolved species in diffuse fluids, and (2) calculation of affinities of potential catabolic reactions in these fluids. These calculations were conducted for actual diffuse fluids sampled and analyzed by Von Damm and Lilley [13] and for hypothetical mixtures of endmember vent fluids [14] and ambient seawater. It is assumed here that the diffuse fluids form by subseafloor mixing of ascending hydrothermal fluids with seawater. The endmember hydrothermal fluid composition is taken from black smoker vent fluids issuing within a few meters of the diffuse vent site [13,14]. The percentage of hydrothermal fluid is estimated using a simple mass balance for silica:

<a onClick="popup('http://www.geochemicaltransactions.com/content/13/1/2/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.geochemicaltransactions.com/content/13/1/2/mathml/M5">View MathML</a>

Silica is known to precipitate slowly at low temperatures from mildly acidic fluids [15] and can be assumed to behave conservatively at the time scales of fluid mixing [16].

Geochemist Workbench® (GWB) was used to conduct the thermodynamic calculations [17]. A Log Kr database was created, covering temperatures from 0 to 350°C at a pressure of 25 MPa, using SUPCRT92 [18] and the thermodynamic database OBIGT [19], and including all speciation reactions in an aqueous system with Na, Ca, Mg, Fe, Sr, K, SiO2, Cl, sulfate, sulfide, oxygen, hydrogen and carbon dioxide. Likewise, equilibrium constants for the following catabolic reactions were calculated.

Aerobic sulfide oxidation

<a onClick="popup('http://www.geochemicaltransactions.com/content/13/1/2/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.geochemicaltransactions.com/content/13/1/2/mathml/M6">View MathML</a>

(1)

Aerobic methane oxidation

<a onClick="popup('http://www.geochemicaltransactions.com/content/13/1/2/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.geochemicaltransactions.com/content/13/1/2/mathml/M7">View MathML</a>

(2)

Aerobic iron oxidation

<a onClick="popup('http://www.geochemicaltransactions.com/content/13/1/2/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.geochemicaltransactions.com/content/13/1/2/mathml/M8">View MathML</a>

(3)

Aerobic hydrogen oxidation

<a onClick="popup('http://www.geochemicaltransactions.com/content/13/1/2/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.geochemicaltransactions.com/content/13/1/2/mathml/M9">View MathML</a>

(4)

Hydrogenotrophic sulfate reduction

<a onClick="popup('http://www.geochemicaltransactions.com/content/13/1/2/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.geochemicaltransactions.com/content/13/1/2/mathml/M10">View MathML</a>

(5)

Hydrogenotrophic methanogenesis

<a onClick="popup('http://www.geochemicaltransactions.com/content/13/1/2/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.geochemicaltransactions.com/content/13/1/2/mathml/M11">View MathML</a>

(6)

Anaerobic oxidation of methane

<a onClick="popup('http://www.geochemicaltransactions.com/content/13/1/2/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.geochemicaltransactions.com/content/13/1/2/mathml/M12">View MathML</a>

(7)

Published compositions of endmember vent fluids [13,14] issuing from black smoker chimneys in proximity (few meters) to the diffuse vent site were used in the calculations of affinities for these reactions (Table 1). In determining Qr (equation2), the extended Debye-Hückel equation was used to calculate activity coefficients with extended parameters and hard core diameters for each species from Wolery and Jove-Colon [20]. Dissolved neutral species were assigned an activity coefficient of one, except non-polar species for which CO2 activity coefficients were used [21]. Reported pH values of hydrothermal vents (measured at 25°C) were used in determining the in situ pH (Table 1) by re-speciating the fluids at the temperatures of venting [22]. The percentage of hydrothermal endmember fluid in the diffuse fluids derived from the silica mass balance was used to calculate idealized mixed fluids, assuming conservative behavior of all elements. These hypothetical fluids were also speciated and compared with actual compositions of diffuse fluids in terms of concentrations and affinities (Tables 1 and 2). Deviations from conservative behavior in the diffuse fluids indicate that removal or release processes take place in the subseafloor mixing zones in which the diffuse fluids are formed. In the calculations, the activities of species in the hypothetical diffuse fluids (ideal conservative mixing) were determined in a batch mixing model simulating titration of hot hydrothermal endmember fluid into cold seawater and tracking the chemical speciation changes in the mixed fluid. In these calculations, redox reactions were suppressed, while kinetically fast reactions like protonation of bases, dissociation of acids, and complex formation are allowed to take place spontaneously. Redox reactions were suppressed, because these reactions are not expected to proceed at the low temperatures of the diffuse fluids and on the short time scales of the mixing process [23]. This procedure has the advantage that disequilibria formed during mixing can be determined and the affinities of selected redox reactions may be calculated. GWB also allows suppressing the Knallgas reaction, so elevated concentrations of both O2 and H2 in the mixed fluids could be accounted for [9].

Table 1. Compositions of discrete and diffuse vent fluids

Table 2. Affinities for different catabolic reactions in kJ and normalized affinities in J per e- and Kg Vent-fluid at the Southern Area (TWP)

The precipitation of minerals was also suppressed. The thermodynamically stable Fe-minerals in the diffuse fluids are hematite and pyrite. If these phases were allowed to precipitate, Fe concentrations would drop to extremely low values in the hypothetical mixed fluids. The measured diffuse fluids have Fe concentrations that are many orders of magnitude higher than values corresponding to pyrite and hematite solubility. They are hence strongly oversaturated with respect to pyrite and hematite and indicate that precipitation of these minerals was largely inhibited. The concentrations of Fe and H2S calculated for the hypothetical mixed fluids represent maximum values.

The affinities calculated (Table 2) represent the maximum energy content for the different catabolic reactions, disregarding the fact that limiting electron donors and acceptors, which appear in several reactions, can still only be used once within the ecosystem [24]. Moreover, comparisons of the raw affinities do not reflect differences in the numbers of electrons transferred in these reactions. This is problematic, because a given quantity of proton motive force driving chemiosmosis is generated by a set number of electrons transferred. The fact that the reactions considered have between one and eight electron transferred therefore skews a comparison of the affinities of different reactions. We hence report the affinities in values per electron transferred (Table 2). Furthermore, the energy flux into the system is controlled by the concentration of the limiting reactant in the upwelling fluid. To examine these combined effects on energy availability, we normalized affinity to kg vent fluid by multiplying the energy with the concentration of the limited reactant, and divided by the fraction of endmember vent fluid in the mix [24]. These normalized affinities provide us with a meaningful parameter for assessing the fluxes of energy for different catabolic reactions into a system. While affinities expressed in both notations are reported in Table 2, the following discussion will primarily use normalized affinity, i.e., energy flux.

Results and Discussion

Case studies

The sample locations are situated in the axial summit caldera of the fast spreading (11 cm/yr full rate) ridge EPR at 9°50'N at a water depth of 2500 meters. In both case studies we use data from time series studies conducted in the 1990s. Venting temperatures are ≤55°C, and fluids issue from cracks in the seafloor or from lava pillars that are fissured near the base. Compositions of vent fluids from the two sites are reported in Von Damm and Lilley [13] and Von Damm [11]. These publications also present a detailed description of the geological setting and vent field characteristics, so we here highlight only the key features of these localities.

The northern area is characterized by the high temperature vents Bio9 and Bio9' and the associated diffuse flow sites BM9Riftia (BM9R), BM91o and BM12. The data set for this system (hereinafter referred to as Bio9 area) is the most detailed, because the site was the target of long-term measurements of fluid composition and temperature [25] and seismic activity [26]. Sohn et al. [26] documented a seismic swarm in 1995 in this area, followed by a temperature increase in the Bio9 vent with a delay of a few days [25]. Temperatures of the diffuse fluid samples range between 22.0 and 33.3°C (Table 1).

The southern area (Tube Worm Pillar, TWP) features high temperature venting through an 11-m high sulfide structure on top of a lava pillar. Discrete venting of 351°C fluid is restricted to the top of the chimney, while leakage of diffuse fluids is observed from around the base of the chimney. Eponymous for the site name, a large tubeworm colony inhabits the area of diffuse venting. The associated diffuse fluid samples were retrieved from Y vent, an adjacent broken-off lava pillar that issued fluids of temperatures between 20 and 25°C in 1992-1995, dropping to 18°C in 1997 and finally to 12°C in 2000.

Case study 1 - Bio9 area

Shank et al. [27] studied the change in the vent community during the time period from 1991 to 1995 at the vents in the Bio9 area. These authors report of a magmatic event in 1991, followed by venting of fluids high in hydrogen sulfide. These conditions boosted the establishment of a strong population of the tubeworm Riftia. During the following cruises in 1994, Shank et al. [27] observed the development of rusty spots that appeared within the Riftia colonies. In 1995, the rusty spots had spread and covered large areas of the Riftia population. In 1997, the Riftia population had broken down largely, while the rust had extended to cover much of the Riftia patch [13].

The temporal evolution of the fluid compositions in that time span reveals a decrease in hydrogen sulfide concentrations over the entire period after 1992 with a slight increase in November of 1995 (Figure 2). Before March of 1994, soluble iron follows the hydrogen sulfide concentration; afterwards the iron content increased and reached maximum concentrations during November of 1995. In November of 1997, the Fe concentration had dropped slightly, but was still much higher than during the beginning of the time series.

thumbnailFigure 2. Temporal changes in concentrations of dissolved iron and hydrogen sulfide in diffuse fluids from the Bio9 area. The increase in iron coincides with the appearance of rusty spots in the tubeworm colony (black arrow). The two horizontal lines represent the maximum concentrations of sulfide (dashed) and iron (continuous) that can be oxidized by seawater with an oxygen concentration of 100 μM. Iron is the limiting reactant over the whole time period, in contrast to sulfide, which is oxygen limited, except for the conditions in November 97.

It has been suggested that the biological development of this area depends on the bioavailability of iron and H2S [13]. This interpretation is plausible, because Riftia live in symbiosis with sulfide oxidizing bacteria [28] and depend on the energy associated with sulfide oxidation. Also, Fe-oxidizing bacteria oxidize Fe2+ in the fluids to ferric hydroxide [29]. So the "rust" in the study area is an indicator that these microorganisms are thriving.

We determined the affinities for both catabolic pathways for the time period of critical geochemical and ecological changes (1991-1997) to improve the understanding of the biological evolution of the vent ecosystems. The calculations make use of the measured concentrations of iron, H2S and oxygen. Unfortunately concentrations of oxygen and the pH for diffuse fluids are not available; therefore, these values are estimated from conservative mixing. Calculated pH values for the fluids show a narrow range of 5.3 to 5.7; likewise, small variations are predicted for oxygen concentrations (92 - 96 μM). Both pH and O2 concentrations reflect the large fraction of seawater calculated from the silica mass balance. Depending on the mixing ratio of vent fluid and seawater, either one of the electron donors (Fe2+, H2S) or oxygen is the limiting reactant determining the amount of energy available per unit vent fluid based (Figure 2). Figure 2 shows the upper limit for iron and H2S oxidation based on an oxygen concentration of the East Pacific bottom seawater of circa 100 μM oxygen [30]. For H2S oxidation, O2 is the limiting reactant, while Fe-oxidation is limited by the availability of iron. An exception is the fluid sampled last in the time series; it exhibits exceptionally low sulfide concentrations and H2S is the compound limiting energy availability.

Figure 3 illustrates the normalized affinities for both reactions. It shows the consequence of limitation; sulfide oxidation has the highest affinity when the fraction of vent fluid in the mixture is lowest (Table 1), because then oxygen contents are greatest. In contrast, the normalized affinities for iron oxidation more closely mirror the iron concentration in the fluid. But affinities are also dependent on the vent fluid fraction, as increased pH favors ferric hydroxide precipitation from the mixed fluids.

thumbnailFigure 3. Normalized affinities for the oxidation of Fe2+ and H2S in diffuse fluids from the Bio9 area during the period from April 1991 until November 1997. The generally high affinities for iron and hydrogen sulfide oxidation support life catabolizing these reactions. In October of 1994, affinities for both reactions are high, so that tubeworms (H2S-oxidizers) and iron oxidizing microorganism (rusty staining) can grow simultaneously. The demise of the Riftia population in November of 1997 coincides with a sudden drop in the affinity of H2S oxidation.

The dynamic changes in the normalized affinities of sulfide and iron oxidation (Figure 3 and Additional File 1) can fully explain the ecological changes within the system. The incipient occurrence of rusty staining in November of 1994 correlates with an increased normalized affinity of iron oxidation, while the normalized affinity for sulfide oxidation remains at the same level. In November of 1995 a further increase in iron concentration in the fluid explains the continued spreading of the iron oxide staining. Tied to this change, the normalized affinity for Fe-oxidation almost quadrupled. The normalized affinity for hydrogen sulfide oxidation was only slightly decreased relative to 1994, which explains why the tubeworm colonies were still thriving, despite the increased development of rusty staining. Apparently, both metabolic pathways were favorable and were being exploited at that stage of system evolution. After 1995, the normalized affinity for hydrogen sulfide oxidation dropped as a consequence of the strongly decreased sulfide concentration in the fluid. Because of this drop in the affinity of sulfide oxidation, the tubeworm population, relying on favorable energetics for H2S oxidation, collapsed. Unlike sulfide oxidation, the normalized affinity for iron oxidation remains high, so organisms with the ability to gain energy from iron oxidation can still thrive. Since both reactions depend on oxygen, the reactions are in competition for that electron acceptor and the calculated affinities (Figure 3) are the predicted maxima.

Additional file 1. Table: Calculated normalized affinities for Aerobic sulfide oxidation and iron oxidation in J per e- and kg Vent-fluid at the Northern Area (Bio9). Additional Data to Figure 3

Format: XLS Size: 16KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

The thermodynamic calculations presented here validate the interpretation by Von Damm and Lilley [13] and confirm that the ecological changes are driven by changes in fluid composition.

Case study 2 - Tube Worm Pillar (TWP)

The fluid compositions of the diffuse fluids issuing in TWP area have been proposed to reveal insights in the redox reactions in the subseafloor [13,31]. Increased methane concentrations in the diffuse fluids led Von Damm and Lilley [13] to propose that hydrogenotrophic methanogenesis takes place in the subseafloor. Proskurowski et al. [31] could confirm this interpretation through carbon stable isotope measurements of methane and CO2 demonstrating that the carbon isotope ratios are consistent with active microbial carbon cycling in this area.

The compositional changes of diffuse fluid compositions relative to the concentrations predicted from conservative mixing are depicted in Figure 4. Throughout the time series, H2 concentrations are decreased by 1.5 to 2 orders of magnitude relative to the concentrations expected from conservative mixing (cf. Table 1). Methane, in contrast, is enriched by a factor of ten relative to the value predicted from conservative mixing in February-March of 1992. In October 1994 and November 1995 this enrichment is about 3-fold. By 2000, measured methane corresponds to those predicted from conservative mixing, and no methane excess can be observed (Figure 4). The methane excess in 1992-1995 is consistent with the decrease in hydrogen, and ratios of H2 depletion to CH4 excess between 3 and 6 are consistent with the stoichiometry of the hydrogenotrophic methanogenesis reaction, from which that ratio would be predicted to be 4. In 1997 and 2000, however, methane excess was minimal and H2 depletion was still significant, suggesting that other hydrogen-consuming reactions may have also played a role.

thumbnailFigure 4. Predicted and measured concentrations of iron, hydrogen sulfide, hydrogen, and methane for diffuse fluids in the Tube Worm Pillar area. Hydrogen is strongly depleted over the entire period. Methane is enriched in the diffuse fluids, which may show methanogenesis in the subseafloor [13]. Until 1997, iron and H2S concentrations are generally lower than predicted from conservative mixing. In November of 1997 predicted and measured concentration are similar to each other. In April of 2000 measured H2S concentrations also correspond to the predictions from conservative mixing, but Fe-concentrations are higher than predicted in the diffuse fluid. Loss of Fe2+ and H2S may be associated with precipitation of minerals. The surplus of measured Fe2+ in 2000 could indicate hydrogenotrophic iron reduction.

While methane enrichment and depletion of hydrogen are indicators for methanogenesis, some of the methane may be metabolized shallower in the system prior to venting by aerobic or anaerobic respiration (Figure 1). There are indications from the Guaymas Basin that anaerobic oxidation of methane (AOM) may take place in vent settings and at temperatures > 30°C [32]. Our calculations suggest that AOM is energetically feasible, so a loss of methane through AOM may be possible. If this reaction took place in the subseafloor, depletions of methane should be associated with increased hydrogen sulfide concentrations. This trend is not observed. AOM may still be taking place, but the rates are too small to affect the compositions of the diffuse fluids.

Affinity calculations for hydrogenotrophic methanogenesis and sulfate reduction in the modeled fluid (Figure 5) indicate strong driving forces for both reactions. The affinity for the reaction is strongly controlled by the hydrogen activity, which has a power of 4 in the relevant mass action equations. Hydrogen endmember concentrations increase from 6460 μM in February/March of 1992 to 8910 μM in March of 1994. Hydrogen concentrations then decrease in the following years to 2700 μM in April of 2000 (Table 1). Unfortunately, the time series contains three points with data lacking for either the diffuse or endmember fluid. In March of 1994, the highest hydrogen concentration in the endmember fluid was measured, but no diffuse fluid was sampled. Hence, in our calculations, the sample collected in October of 1994 has the highest predicted hydrogen content (Figure 4) and also the highest normalized affinity for methanogenesis with 34.0 Joule per kg vent fluid and electron transferred in reaction (J/kg e-) or of 41.5 J/kg e- for sulfate reduction in the predicted diffuse fluid (Figure 5). The fluids sampled in April of 2000 have the smallest fraction of vent fluid and the lowest hydrogen endmember concentration of 12.5 μM, yielding normalized affinities for methanogenesis of 8.9 J/kg e- and 11.0 J/kg e- for sulfate reduction. In these diffuse fluids, hydrogen concentrations are still lower than predicted from conservative mixing and they do not correlate with the endmember concentration or extent of mixing with seawater.

thumbnailFigure 5. Normalized affinities for sulfate reduction and methanogenesis in the Tube Worm Pillar area. Affinities are high in the hypothetical fluids calculated from conservative mixing of seawater and hydrothermal fluid. In the measured fluids the affinity is strongly decreased. Notably, affinities drop markedly in the first three years, which reflects the decrease in H2 concentration in that time span (Figure 4). The removal of H2 and the lowering of affinity reflect the exploitation of H2 in fueling catabolic activity. During the last five years of the time series, normalized affinities had plateaued, possibly indicating a steady-state between hydrothermal energy supply and microbial utilization of energy.

Overall, hydrogen concentrations decrease linearly within the first four years of the time series from 14.9 μM to 0.7 μM. Hydrogen concentrations then remain fairly constant in the range of 0.7 to 0.2 μM until 2000. Figure 5 shows the affinities of these reactions per mol electrons in the reaction and normalized to kg of vent fluid. The normalized affinities of sulfate reduction decrease from 1.4 J/kg e- to 0.05 J/kg e- in 1995 and rebounds to 0.08-0.09 J/kg e- in the following years. Hydrogenotrophic methanogenesis has an affinity of 1.0 J/kg e- in 1992 and decreases to 0.05 to 0.06 J/kg e- in 1995 - 2000 (Table 2).

The differences in the affinities (Table 2) reflect the variability in hydrogen concentrations, but the magnitude of these differences is quite small, because the intensive term (ΔrG°) in the Gibbs energy calculation is very large for both reactions. The calculated affinities for methanogenesis and sulfate reduction (Table 2) lie above the estimated energy limit of microbial metabolism 10 kJ/8 e- [33]. Minimum H2 concentrations required for microbial harnessing of hydrogen at the temperatures of diffuse venting (11.9 to 24.7°C) is on the order of 10-8 M [33], which is considerably lower than the measured concentrations (> 2 × 10-7 M). This result indicates that the microbial communities consume hydrogen but do not control its abundance. Apparently, hydrothermally driven influx of H2 into the system is overall greater than the rate at which H2 is metabolized. The constant hydrogen concentrations between 1995 and 2000 probably indicate some sort of steady-state between influx of hydrogen from below and hydrogen consumption in the subseafloor ecosystem. The early phase of decreasing hydrogen concentrations in the diffuse fluids is not related to changes in endmember compositions (steady between 250 and 350 μM), but may instead reflect the growth a hydrogenotrophic microbial community and increasing rates of consumption of H2 advected into the system by hydrothermal flow. In that early stage, there was a relation between H2 depletion and methane production, indicating that methanogenesis was responsible for both. In 1997 and 2000, H2 was still consumed, but the methane excess had disappeared. Instead, there was excess Fe in the fluids, suggesting that Fe-reduction was taking place, perhaps because it requires lower H2 activities than methanogenesis [33].

Conclusions

Thermodynamic calculations of energy yields of catabolic reactions from geochemical data of diffuse fluids facilitate an assessment of microbial metabolism in vent settings. This has been demonstrated in two case studies, both from the EPR 9°50'N region, where published geochemical data [13,14] where used in systematic calculations of affinities of different catabolic reactions.

In the Bio9 area, affinities for sulfide oxidation strongly decrease, which is in accordance with the dying Riftia population. At the same time, an increase in the affinity for iron oxidation corresponds to a massive spread of red staining in the area, which is likely evidence for Fe-oxidizing bacteria. The results of the energy calculations verify the idea that the sudden change in vent fauna is a result of changes in fluid chemistry.

The example from the Bio9 area is more relevant to subseafloor processes. Enrichment of methane in diffuse fluids points to methanogenesis in the mixing and cooling zone. Our calculations confirm that hydrogenotrophic catabolic reactions have large energy yields throughout the duration of the time series (1992-2000). A large discrepancy in the amount of H2 predicted from conservative mixing and the measured H2 concentrations indicate effective scrubbing of H2 by subseafloor hydrogenotrophic microorganisms. During the first three years of the time series, affinities for hydrogenotrophic reactions decreased despite continued high H2 concentrations in the endmember fluids. This is interpreted to indicate the development of a hydrogenotrophic-based microbial ecosystem in the subseafloor. Between 1995 and 2000, the affinities remained constant and low (about an order of magnitude above the biological energy quantum). Apparently, influx of hydrogen from below and consumption of hydrogen within the subseafloor had reached a steady state. In 1997 and 2000, methane excesses were minimal, but the fluids showed pronounced enrichment of Fe relative to the concentrations predicted from conservative mixing. This finding may indicate a switch within the system from methanogenesis to Fe-reduction.

Our results show how thermodynamic calculations can be used to examine the relations between changes in fluid chemistry and seafloor biology. They are also a helpful tool in examining processes in the subseafloor and help highlight the tight relations and interdependencies between geochemistry and microbiology in vent systems.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

The study was developed jointly by both authors. MH conducted the thermodynamic calculations. Both authors were involved equally in the interpretations of the results and in writing and editing of the manuscript. All authors have read and approved the final manuscript.

Acknowledgements

We are grateful to two anonymous reviewers for numerous helpful comments and suggestions and to Greg Druschel for efficient editorial handling. This work was funded through the German Research Foundation (Grant BA1605/6-1). WB acknowledges additional support through the University of Bremen and the MARUM Research Center and Cluster of Excellence.

References

  1. Jannasch HW: Microbial interactions with hydrothermal fluids. In Seafloor Hydrothermal Systems: Physical, Chemical, Biological, and Geological Interactions. Volume 91. Edited by Humphris SE, et al. Washington, DC: American Geophysical Union, American Geophysical Union; 1995::273-273.

    Geophysical Monograph

    OpenURL

  2. McCollom TM, Shock EL: Geochemical constraints on chemolithoautotrophic metabolism by microorganisms in seafloor hydrothermal systems.

    Geochimica et Cosmochimica Acta 1997, 61:4375-4391. PubMed Abstract | Publisher Full Text OpenURL

  3. Amend JP, McCollom TM, Hentscher M, Bach W: Catabolic and anabolic energy for chemolithoautotrophs in deep-sea hydrothermal systems hosted in different rock types.

    Geochimica et Cosmochimica Acta 2011, 75:5736-5748. Publisher Full Text OpenURL

  4. Shock E, Holland ME: Geochemical energy sources that support the subsurface biosphere. In The Subseafloor Biosphere at Mid-Ocean Ridges. Edited by Wilcock WSD, DeLong EF, Kelley DS, Baross JA, Cary SC. Washington, DC: American Geophysical Union; 2004:153-165.

    Geophysical Monograph

    OpenURL

  5. Houghton JL, Seyfried WE Jr: An experimental and theoretical approach to determining linkages between geochemical variability and microbial biodiversity in seafloor hydrothermal chimneys.

    Geobiology 2010, 8:457-470. PubMed Abstract | Publisher Full Text OpenURL

  6. Huber JA, Butterfield DA, Baross JA: Bacterial diversity in a subseafloor habitat following a deep-sea volcanic eruption.

    FEMS Microbiology Ecology 2003, 43:393-409. PubMed Abstract | Publisher Full Text OpenURL

  7. Reysenbach A-L, Liu Y, Banta AB, Beveridge TJ, Kirshtein JD, Schouten S, Tivey MK, Von Damm KL, Voytek MA: A ubiquitous thermoacidophilic archaeon from deep-sea hydrothermal vents.

    Nature 2006, 442:444-447. PubMed Abstract | Publisher Full Text OpenURL

  8. Kormas KA, Tivey MK, Von Damm K, Teske A: Bacterial and archaeal phylotypes associated with distinct mineralogical layers of a white smoker spire from a deep-sea hydrothermal vent site (9°N, East Pacific Rise).

    Environmental Microbiology 2006, 8:909-920. PubMed Abstract | Publisher Full Text OpenURL

  9. Perner M, Bach W, Hentscher M, Koschinsky A, Garbe-Schönberg D, Streit WR, Strauss H: Short-term microbial and physico-chemical variability in low-temperature hydrothermal fluids near 5°S on the Mid-Atlantic Ridge.

    Environmental Microbiology 2009, 11:2526-2541. PubMed Abstract | Publisher Full Text OpenURL

  10. Petersen JM, Zielinski FU, Pape T, Seifert R, Moraru C, Amann R, Hourdez S, Girguis PR, Wankel SD, Barbe V, et al.: Hydrogen is an energy source for hydrothermal vent symbioses.

    Nature 2011, 476:176-180. PubMed Abstract | Publisher Full Text OpenURL

  11. Amend JP, Shock EL: Energetics of overall metabolic reactions of thermophilic and hyperthermophilic Archaea and Bacteria.

    FEMS Microbiology Reviews 2001, 25:175-243. PubMed Abstract | Publisher Full Text OpenURL

  12. Helgeson HC: Mass transfer among minerals and hydrothermal solutions.

    Geochemistry of hydrothermal ore deposits 1979, 2:568-610. OpenURL

  13. Von Damm KL, Lilley MD: Diffuse flow hydrothermal fluids from 9 50'N East Pacific Rise: evolution and biogeochemical controls. In The Subseafloor Biosphere at Mid-Ocean Ridges. Edited by Wilcock WSD, DeLong EF, Kelley DS, Baross JA, Cary SC. Washington, DC: American Geophysical Union; 2004:245-268.

    Geophysical Monograph

    OpenURL

  14. Von Damm KL: Evolution of the hydrothermal system at East Pacific Rise 9 50'N: Geochemical evidence for changes in the upper oceanic crust. In Mid-Ocean Ridges: Hydrothermal Interactions between the Lithosphere and Oceans. Edited by German CR, Lin J, Parson LM. Washington, DC: American Geophysical Union; 2004:285-304.

    Geophysical Monograph

    OpenURL

  15. Carroll S, Mroczek E, Alai M, Ebert M: Amorphous silica precipitation (60 to 120°C): comparison of laboratory and field rates.

    Geochimica et Cosmochimica Acta 1998, 62:1379-1396. Publisher Full Text OpenURL

  16. Edmond JM, Measures C, McDuff RE, Chan LH, Collier R, Grant B, Gordon LI, Corliss JB: Ridge crest hydrothermal activity and the balances of the major and minor elements in the ocean: The Galapagos data.

    Earth and Planetary Science Letters 1979, 46:1-18. Publisher Full Text OpenURL

  17. Bethke C: Geochemical reaction modeling: Concepts and applications. Oxford University Press, USA; 1996.

  18. Johnson JW, Oelkers EH, Helgeson HC: SUPCRT92: A software package for calculating the standard molal thermodynamic properties of minerals, gases, aqueous species, and reactions from 1 to 5000 bar and 0 to 1000°C.

    Computers & Geosciences 1992, 18:899-947. PubMed Abstract | Publisher Full Text OpenURL

  19. Dick J: Calculation of the relative metastabilities of proteins using the CHNOSZ software package.

    Geochemical Transactions 2008, 9. OpenURL

  20. Wolery T, Jove-Colon C: Qualification of thermodynamic data for geochemical modeling of mineral-water interactions in dilute systems. In Energy. Edited by US Department of Energy. Bechtel SAIC Company, Las Vegas, Nevadam, USA; 2004. OpenURL

  21. Drummond SE: Boiling and mixing of hydrothermal fluids. In Ph D thesis. University Microfilms Int.; 1981. OpenURL

  22. Tivey MK, Humphris SE, Thompson G, Hannington MD, Rona PA: Deducing patterns of fluid flow and mixing within the TAG active hydrothermal mound using mineralogical and geochemical data.

    J Geophys Res 1995, 100:12527-12555. Publisher Full Text OpenURL

  23. Foustoukos DI, Houghton JL, Seyfried WE Jr, Sievert SM, Cody GD: Kinetics of H2-O2-H2O redox equilibria and formation of metastable H2O2 under low temperature hydrothermal conditions.

    Geochimica et Cosmochimica Acta 2011, 75:1594-1607. Publisher Full Text OpenURL

  24. McCollom TM: Geochemical constraints on sources of metabolic energy for chemolithoautotrophy in ultramafic-hosted deep-sea hydrothermal systems.

    Astrobiology 2007, 7:933-950. PubMed Abstract | Publisher Full Text OpenURL

  25. Fornari DJ, Shank T, Von Damm KL, Gregg TKP, Lilley M, Levai G, Bray A, Haymon RM, Perfit MR, Lutz R: Time-series temperature measurements at high-temperature hydrothermal vents, East Pacific Rise 9°49'-51'N: evidence for monitoring a crustal cracking event.

    Earth and Planetary Science Letters 1998, 160:419-431. Publisher Full Text OpenURL

  26. Sohn RA, Fornari DJ, Von Damm KL, Hildebrand JA, Webb SC: Seismic and hydrothermal evidence for a cracking event on the East Pacific Rise crest at 9° 50' N.

    Nature 1998, 396:159-161. Publisher Full Text OpenURL

  27. Shank TM, Fornari DJ, Von Damm KL, Lilley MD, Haymon RM, Lutz RA: Temporal and spatial patterns of biological community development at nascent deep-sea hydrothermal vents (9°50'N, East Pacific Rise).

    Deep Sea Research Part II: Topical Studies in Oceanography 1998, 45:465-515. Publisher Full Text OpenURL

  28. Cavanaugh CM: Symbiotic chemoautotrophic bacteria in marine invertebrates from sulphide-rich habitats.

    Nature 1983, 302:58-61. Publisher Full Text OpenURL

  29. Toner BM, Santelli CM, Marcus MA, Wirth R, Chan CS, McCollom T, Bach W, Edwards KJ: Biogenic iron oxyhydroxide formation at mid-ocean ridge hydrothermal vents: Juan de Fuca Ridge.

    Geochimica et Cosmochimica Acta 2009, 73:388-403. Publisher Full Text OpenURL

  30. Betts JN, Holland HD: The oxygen content of ocean bottom waters, the burial efficiency of organic carbon, and the regulation of atmospheric oxygen.

    Palaeogeography, Palaeoclimatology, Palaeoecology 1991, 97:5-18. Publisher Full Text OpenURL

  31. Proskurowski G, Lilley MD, Olson EJ: Stable isotopic evidence in support of active microbial methane cycling in low-temperature diffuse flow vents at 9°50'N East Pacific Rise.

    Geochimica et Cosmochimica Acta 2008, 72:2005-2023. Publisher Full Text OpenURL

  32. Schouten S, Wakeham SG, Hopmans EC, Sinninghe Damste JS: Biogeochemical evidence that thermophilic archaea mediate the anaerobic oxidation of methane.

    Appl Environ Microbiol 2003, 69:1680-1686. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Hoehler TM: Biological energy requirements as quantitative boundary conditions for life in the subsurface.

    Geobiology 2004, 2:205-215. Publisher Full Text OpenURL