Kump et al EPSL 2005.pdf - Bryn Mawr College
Kump et al EPSL 2005.pdf - Bryn Mawr College
Kump et al EPSL 2005.pdf - Bryn Mawr College
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Earth and Plan<strong>et</strong>ary Science L<strong>et</strong>ters 235 (2005) 654–662<br />
www.elsevier.com/locate/epsl<br />
Hydrotherm<strong>al</strong> Fe fluxes during the Precambrian: Effect of low<br />
oceanic sulfate concentrations and low hydrostatic pressure<br />
on the composition of black smokers<br />
Lee R. <strong>Kump</strong> a, *, William E. Seyfried Jr. b<br />
a Department of Geosciences and Astrobiology Research Center, The Pennsylvania State University, 535 Deike Bldg.,<br />
University Park, PA 16827, USA<br />
b Department of Geology and Geophysics, University of Minnesota, 310 Pillsbury Dr. SE, Minneapolis, MN 55455-0219, USA<br />
Received 18 June 2004; received in revised form 15 March 2005; accepted 25 April 2005<br />
Available online 15 June 2005<br />
Editor: E. Boyle<br />
Abstract<br />
Modern mid-ocean ridge hydrotherm<strong>al</strong> systems typic<strong>al</strong>ly release vent fluids to the ocean with dissolved H 2 S in excess of Fe.<br />
These fluids are the consequence of high-temperature interactions b<strong>et</strong>ween sulfate-rich seawater and mid-ocean ridge bas<strong>al</strong>t at<br />
conditions near the critic<strong>al</strong> point for seawater. The precipitation of FeS and FeS 2 during venting titrates most of the Fe from the<br />
fluid, significantly reducing the n<strong>et</strong> flux of Fe to the open ocean. Here we suggest that hydrotherm<strong>al</strong> fluids emanating from<br />
Precambrian seafloor systems older than ~1.8 Ga had Fe/H 2 S ratiosNN1, and with fH 2 higher than today, because seawater<br />
lacked its primary oxidant, dissolved sulfate ion. This predominance of Fe over H 2 S would have promoted the establishment of<br />
an iron-rich deep ocean and the deposition of banded iron formations (BIF). Accordingly, the end of BIF deposition at ~1.8 Ga<br />
was the result of the buildup of sulfate in seawater from oxidative weathering, and its r<strong>et</strong>urn at 750 Ma the result of reductions in<br />
seawater sulfate concentrations during Snowb<strong>al</strong>l Earth episodes, enhanced by elevated Fe concentrations during depressurization<br />
of hydrotherm<strong>al</strong> systems by large eustatic sea-level f<strong>al</strong>ls. Moreover, Precambrian chemosynth<strong>et</strong>ic vent communities may<br />
have been based on H 2 synthesis rather than on H 2 S oxidation, as they largely are today.<br />
D 2005 Elsevier B.V. All rights reserved.<br />
Keywords: hydrotherm<strong>al</strong>; banded iron formation; Precambrian; iron; sulfate<br />
1. Introduction<br />
* Corresponding author.<br />
E-mail address: lkump@psu.edu (L.R. <strong>Kump</strong>).<br />
Mid-ocean ridge vent fluids have been considered<br />
a likely source of Fe for Precambrian iron formations<br />
[1]. However, modern vent compositions are unlikely<br />
to generate BIF because the concentration of H 2 S<br />
0012-821X/$ - see front matter D 2005 Elsevier B.V. All rights reserved.<br />
doi:10.1016/j.epsl.2005.04.040
L.R. <strong>Kump</strong>, W.E. Seyfried Jr. / Earth and Plan<strong>et</strong>ary Science L<strong>et</strong>ters 235 (2005) 654–662 655<br />
gener<strong>al</strong>ly exceeds the Fe concentration such that Fe<br />
is titrated out of the vent fluid by the precipitation of<br />
FeS and FeS 2 [2,3]. W<strong>al</strong>ker and Brimblecombe [4]<br />
pointed out, however, that in the absence of seawater<br />
sulfate, vent fluids would be Fe dominated simply<br />
because bas<strong>al</strong>t has an abundance of Fe-miner<strong>al</strong>s<br />
relative to S-miner<strong>al</strong>s.<br />
In a similar vein, Canfield [5] suggested that the<br />
end of BIF deposition at ~1.8 Ga resulted from the<br />
accumulation of sulfate in the ocean following the rise<br />
in atmospheric oxygen abundance at ~2.3 Ga. In<br />
Canfield’s model, microbi<strong>al</strong> reduction of this sulfate<br />
led to the development of sulfidic conditions in the<br />
deep sea that severely reduced the concentration of<br />
ferrous iron because of the insolubility of FeS. An<br />
<strong>al</strong>ternative explanation, consistent with W<strong>al</strong>ker and<br />
Brimblecombe’s model, is that abiogenic sulfide associated<br />
with hydrotherm<strong>al</strong> venting of P<strong>al</strong>eoproterozoic<br />
sulfate-rich seawater-derived fluids led to the<br />
precipitation of iron-sulfide miner<strong>al</strong>s, creating conditions<br />
unfavorable for BIF production (with Fe/<br />
H 2 Sb1). According to either model, the r<strong>et</strong>urn of<br />
BIF deposition during the Neoproterozoic bSnowb<strong>al</strong>l<br />
EarthQ episodes would have been more related to the<br />
depl<strong>et</strong>ion of seawater sulfate during the ice-covered<br />
interv<strong>al</strong> than to the establishment of oceanic anoxia<br />
(the origin<strong>al</strong> argument put forth by Kirschvink [6] and<br />
Beukes and Klein [7]); anoxia <strong>al</strong>one does not affect<br />
the Fe/H 2 S ratio of vent fluids nor does it <strong>al</strong>low for<br />
the buildup of Fe because of the high concentrations<br />
of H 2 S that develop.<br />
In this paper we explore the W<strong>al</strong>ker and Brimblecombe’s<br />
model [4] using available thermodynamic<br />
data for seawater equilibria at elevated temperatures<br />
and pressures. Assuming the absence of seawater<br />
sulfate, we find that Precambrian vent fluids were<br />
considerably more reducing, with ratios of Fe/H 2 S<br />
that exceeded unity. Low hydrostatic pressures during<br />
Snowb<strong>al</strong>l Earth episodes of the Neoproterozoic and<br />
for much of the Archean, if mid-ocean ridge depths<br />
were sh<strong>al</strong>lower then than now [1], may have further<br />
enhanced Fe concentrations.<br />
2. Constraints on vent fluid composition<br />
After nearly 20 years of experimentation and theor<strong>et</strong>ic<strong>al</strong><br />
development, we now have a clearer understanding<br />
of the effects of temperature and pressure on<br />
the chemistry of mid-ocean ridge hydrotherm<strong>al</strong> systems<br />
[8,9]. The composition of fluids collected from<br />
black smoker chimneys can be explained for the most<br />
part in terms of fluid–miner<strong>al</strong> equilibria, established at<br />
elevated temperatures and pressures in subseafloor<br />
reaction zones near magma chambers below the<br />
chimneys. These solutions then rise convectively to<br />
the sea floor, mixing with cold seawater, precipitating<br />
insoluble miner<strong>al</strong>s, and in some cases, separating into<br />
vapor and brine [10,11].<br />
2.1. Effect of seawater sulfate on vent fluid<br />
composition<br />
In sharp contrast with modern hydrotherm<strong>al</strong> systems,<br />
low dissolved sulfate in the ocean [5,12,13]<br />
would likely have rendered Precambrian hydrotherm<strong>al</strong><br />
systems more reducing, enhancing dissolved Fe<br />
concentrations in coexisting vent fluids. We can illustrate<br />
this by means of a phase diagram for a portion of<br />
the FeO–Fe 2 O 3 –H 2 S–SiO 2 –CaO–H 2 O–HCl system at<br />
400 8C, 400 bars (Fig. 1). Reaction of bas<strong>al</strong>t/gabbro or<br />
even more reducing protolith with a fluid lacking<br />
sulfate would permit the inherent redox capacity of<br />
the rock to buffer the fluid resulting in relatively high<br />
H 2 /H 2 S ratios, consistent with constraints imposed by<br />
ferrous iron-bearing phases, such as might be approximated<br />
by the fay<strong>al</strong>ite–magn<strong>et</strong>ite–pyrrhotite-bearing<br />
system. Vent fluids impacted by magmatic degassing<br />
effects provide an indication of this. For example,<br />
hydrotherm<strong>al</strong> vent fluids from 98 to 108 N, East<br />
Pacific Rise and the Endeavour segment of the Juan<br />
de Fuca Ridge, indicate high H 2 /H 2 S ratios and high<br />
H 2 concentrations in the immediate aftermath of subseafloor<br />
magmatic intrusions [14]. Although uncertainties<br />
exist in terms of temperature, pressure and<br />
phase separation effects, the high H 2 in particular is<br />
an indicator of distinctly reducing conditions associated<br />
with vapor release at the magmatic–hydrotherm<strong>al</strong><br />
interface [14,15]. In <strong>al</strong>l cases, however, time series<br />
observations reve<strong>al</strong> sharp decreases in dissolved H 2<br />
and H 2 S, with changes in H 2 greater than H 2 S [14],<br />
suggesting more oxidizing conditions. Indeed, even<br />
modest rock–seawater interaction can be expected to<br />
cause the fluid to achieve saturation with respect to<br />
anhydrite, assuming coexistence of plagioclase feldspar<br />
and quartz (Fig. 1), which is in good agreement
656<br />
L.R. <strong>Kump</strong>, W.E. Seyfried Jr. / Earth and Plan<strong>et</strong>ary Science L<strong>et</strong>ters 235 (2005) 654–662<br />
-1.0<br />
T = 400 °C<br />
P = 400 bars<br />
Cl = 0.55 mol<strong>al</strong><br />
-1.5<br />
Pyrite<br />
Pyrrhotite<br />
log H 2 S(aq) (mol<strong>al</strong>)<br />
-2.0<br />
-2.5<br />
Hematite<br />
Anhydrite<br />
-2<br />
SO 4 - bearing<br />
seawater<br />
Magn<strong>et</strong>ite<br />
-2<br />
SO 4 - free seawater<br />
Fay<strong>al</strong>ite<br />
-3.0<br />
-4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0<br />
log H 2(aq) (mol<strong>al</strong>)<br />
Fig. 1. Phase equilibria in the FeO–Fe 2 O 3 –H 2 S–SiO 2 –CaO–H 2 O–HCl system at 400 8C, 400 bars. Reaction of fresh bas<strong>al</strong>t and/or ultramafic<br />
lithologies with sulfate-free aqueous fluid can be expected to result in relatively high H 2 and moderate H 2 S concentrations, respectively,<br />
an<strong>al</strong>ogous to constraints imposed by fay<strong>al</strong>ite–magn<strong>et</strong>ite–pyrrhotite–fluid equilibria. In contrast, the relative abundance of sulfate in modern<br />
seawater provides a powerful oxidizing agent that is capable of buffering H 2 at lower v<strong>al</strong>ues, <strong>al</strong>though H 2 S concentrations are gener<strong>al</strong>ly similar<br />
to the more reducing system. For plagioclase-bearing systems, these redox constraints ultimately result in anhydrite formation, a condition that<br />
is consistent with the chemistry of many modern vent fluids (diagon<strong>al</strong> boundary). The two distinctly different redox conditions—reducing/<br />
initi<strong>al</strong>ly sulfate-free fluid (fay<strong>al</strong>ite/magn<strong>et</strong>ite/pyrrhotite equilibria), and more oxidizing/initi<strong>al</strong>ly sulfate-bearing fluid (anhydrite/magn<strong>et</strong>ite<br />
equilibria) can exert a fundament<strong>al</strong> control on dissolved Fe concentrations (see below). Modern vent fluids unaffected of recent subseafloor<br />
magmatic activity reve<strong>al</strong> H 2 /H 2 S concentrations that are in good agreement with the plagioclase buffered, anhydrite-bearing system [8,14].<br />
Thermodynamic data used for the construction of the diagram are from Johnson <strong>et</strong> <strong>al</strong>. [41]. Activity–concentration relations for pyrrhotite and<br />
H 2 and H 2 S are from Barton and Skinner [42], Kishima [43] and Ding and Seyfried [44], respectively.<br />
with the vent fluid dissolved gas data [8,14,16]. Accordingly,<br />
modern vent fluids with dissolved chloride<br />
concentrations close to seawater v<strong>al</strong>ues have dissolved<br />
H 2 and H 2 S concentrations that typic<strong>al</strong>ly do<br />
not exceed 1 and 10 mmol/kg, respectively [3,14].<br />
The role of dissolved sulfate in the redox evolution of<br />
modern vent fluids, however, is <strong>al</strong>so manifest by d 34 S<br />
data, which reve<strong>al</strong> seawater and rock-derived sources<br />
of sulfur [17], while pen<strong>et</strong>ration of sufficient sulfate to<br />
render anhydrite stable in high-temperature hydrotherm<strong>al</strong><br />
reaction zones is consistent with d 34 S of sulfate in<br />
hydrotherm<strong>al</strong>ly <strong>al</strong>tered rocks [18].<br />
Redox constraints imposed by dissolved sulfate not<br />
only affect dissolved H 2 and H 2 S concentrations in<br />
vent fluids, but <strong>al</strong>so dissolved Fe. For example, c<strong>al</strong>culations<br />
for typic<strong>al</strong> seafloor hydrotherm<strong>al</strong> vent fluids<br />
(400 8C, 400 bars, and an in situ pH of 5; [19])<br />
indicate relatively high Fe concentrations for the fay<strong>al</strong>ite–pyrrhotite-bearing<br />
system, whereas more modest<br />
Fe concentrations are predicted for sulfate-bearing<br />
systems (Fig. 2a and b). Although dissolved Fe concentrations<br />
change greatly from one condition to the<br />
other, dissolved H 2 S remains relatively constant (see<br />
Fig. 1). Thus, the absence of sulfate in seawater would<br />
likely yield vent fluids with high Fe/H 2 S ratios,<br />
approaching v<strong>al</strong>ues imposed by the reducing nature<br />
of the un<strong>al</strong>tered rock. The high ratios depicted in (a)<br />
are likely even with modest increases in pH.<br />
A consequence of the high Fe /H 2 S ratio vent<br />
fluids would be more efficient delivery of Fe to the<br />
ocean. For example, c<strong>al</strong>culations show that mixing a<br />
fluid initi<strong>al</strong>ly buffered at 400 8C by fay<strong>al</strong>ite–magn<strong>et</strong>ite–pyrrhotite<br />
equilibria with sulfate-free seawater, as<br />
would occur on venting, fails to prevent the flux of Fe
L.R. <strong>Kump</strong>, W.E. Seyfried Jr. / Earth and Plan<strong>et</strong>ary Science L<strong>et</strong>ters 235 (2005) 654–662 657<br />
80<br />
70<br />
Fe<br />
a.<br />
Concentration (mmol/kg)<br />
60<br />
50<br />
40<br />
30<br />
20<br />
10<br />
0<br />
12<br />
11<br />
10<br />
H 2 S (aq)<br />
5.0 5.2 5.4 5.6 5.8 6.0<br />
pH<br />
H 2 S<br />
b.<br />
Concentration (mmol/kg)<br />
70<br />
60<br />
50<br />
40<br />
30<br />
20<br />
10<br />
0<br />
a.<br />
Fe<br />
H 2 S<br />
0 50 100 150 200 250 300 350 400<br />
Temperature °C<br />
Concentration (mmol/kg)<br />
9<br />
8<br />
7<br />
6<br />
5<br />
4<br />
3<br />
2<br />
1<br />
Fe<br />
Concentration (mmol/kg)<br />
10<br />
8<br />
6<br />
4<br />
2<br />
b.<br />
H 2 S<br />
Fe<br />
5.0 5.2 5.4 5.6 5.8 6.0<br />
pH<br />
0<br />
Fig. 2. The effect of pH on dissolved Fe for fay<strong>al</strong>ite–magn<strong>et</strong>ite–<br />
pyrrhotite–(quartz)–fluid equilibria (a) and anhydrite–magn<strong>et</strong>ite–<br />
(quartz)–plagioclase (An80) [8]. (b). C<strong>al</strong>culations were performed<br />
at 400 8C, 400 bars. At a pH of 5, which is typic<strong>al</strong> of modern hot<br />
spring vent fluids at mid-ocean ridges [15], high Fe and relatively<br />
low H 2 S are predicted (high Fe/H 2 S ratio) assuming equilibria with<br />
the more reducing assemblage (a), while the opposite is true for the<br />
anhydrite-bearing bmodernQ system (b). The relative absence of<br />
sulfate in ancient oceans (Archean, Neoproterozoic) would permit<br />
more reducing assemblages in host rocks to persist, enhancing Fe<br />
solubility. Model c<strong>al</strong>culations were performed using EQ3/6 [45]<br />
with thermodynamic data generated using SUPCRT92 [41,46],<br />
assuming dissolved chloride of 0.55 mol/kg. See Fig. 1 for addition<strong>al</strong><br />
sources of thermodynamic data.<br />
into the ancient ocean. Most of the loss in Fe is by<br />
dilution, <strong>al</strong>though minor Fe-miner<strong>al</strong>ization (pyrite,<br />
pyrrhotite and at sufficiently low temperature, hematite)<br />
<strong>al</strong>so occurs (Fig. 3a). Because of the high Fe/H 2 S<br />
ratio, however, compl<strong>et</strong>e remov<strong>al</strong> of H 2 S from the<br />
0 50 100 150 200 250 300 350 400<br />
Temperature °C<br />
Fig. 3. Reaction path model depicting the effect of mixing (cooling)<br />
on dissolved Fe and H 2 S initi<strong>al</strong>ly s<strong>et</strong> assuming fay<strong>al</strong>ite–magn<strong>et</strong>ite–<br />
pyrrhotite–(quartz)–fluid equilibria (a) and anhydrite–pyrite–magn<strong>et</strong>ite–(quartz)–plagioclase<br />
(An80) (b) (see Fig. 2). Concentrations<br />
of Fe and H 2 S at 400 8C, 400 bars were c<strong>al</strong>culated assuming pH=5,<br />
and 0.55 mol/kg dissolved chloride. Temperature change was c<strong>al</strong>culated<br />
assuming mixing with a NaCl fluid (0.55 mol/kg) at 25 8C<br />
(a), while modern seawater was the low temperature mix fluid for<br />
second simulation (b). Dilution effects and temperature dependent<br />
changes in sulfide miner<strong>al</strong> solubility (pyrite, pyrrhotite, hematite)<br />
and homogeneous equilibria (pH change) cause the predicted<br />
changes in Fe and H 2 S (a). These effects can result in the delivery<br />
of relatively high Fe and high Fe/H 2 S ratio fluids to the ancient<br />
ocean affecting BIF deposition. This is not the case for modern<br />
sulfate-bearing systems due to the initi<strong>al</strong>ly low Fe/H 2 S ratio of the<br />
predicted source fluid (b). Model c<strong>al</strong>culations were performed using<br />
EQ3/6 [45] with thermodynamic data generated using SUPCRT92<br />
[41,46], assuming dissolved chloride of 0.55 mol/kg. See Fig. 1 for<br />
addition<strong>al</strong> sources of thermodynamic data.
658<br />
L.R. <strong>Kump</strong>, W.E. Seyfried Jr. / Earth and Plan<strong>et</strong>ary Science L<strong>et</strong>ters 235 (2005) 654–662<br />
hydrotherm<strong>al</strong> fluid is predicted (Fig. 3a). In contrast,<br />
the low Fe/H 2 S ratio of the more oxidizing modern<br />
vent fluids results in compl<strong>et</strong>e remov<strong>al</strong> of Fe when<br />
mixed with modern sulfate-bearing seawater, a result<br />
largely due to pyrite precipitation (Fig. 3b). In actu<strong>al</strong>ity,<br />
vent miner<strong>al</strong>ogy is undoubtedly dominated by a<br />
complex series of kin<strong>et</strong>ic<strong>al</strong>ly mediated reactions [2],<br />
<strong>al</strong>though fundament<strong>al</strong> redox constraints will still dominate<br />
the temperature dependent abundance and sequence<br />
of miner<strong>al</strong> precipitates and corresponding<br />
changes in vent fluid chemistry.<br />
2.2. Effect of pressure on vent fluid composition<br />
concentration (mmol/kg)<br />
80 500 bars<br />
400 bars<br />
70<br />
60<br />
50<br />
40<br />
30<br />
20<br />
10<br />
It has long been known that in addition to redox,<br />
temperature and pressure play critic<strong>al</strong> roles in controlling<br />
Fe solubility in hydrotherm<strong>al</strong> systems. At<br />
elevated temperatures, relatively sm<strong>al</strong>l changes in<br />
pressure can result in large changes in dissolved<br />
Fe. Thermodynamic c<strong>al</strong>culations for the fay<strong>al</strong>ite–pyrrhotite–magn<strong>et</strong>ite–(quartz)-bearing<br />
system at pH=5<br />
show that a change in pressure from 500 to 400<br />
bars results in a doubling of Fe, a modest change<br />
in H 2 S, and correspondingly a substanti<strong>al</strong> increase in<br />
Fe/H 2 S (Fig. 4). Owing to the lack of thermodynamic<br />
data, it is not possible to perform similar<br />
c<strong>al</strong>culations at lower pressures, <strong>al</strong>though it can be<br />
inferred from experiment<strong>al</strong> data [20] that at temperatures<br />
in excess of 350 8C, a further decrease in<br />
pressure would result in addition<strong>al</strong> increases in dissolved<br />
Fe, particularly for hydrotherm<strong>al</strong> systems<br />
lacking sulfate. This sensitivity largely results from<br />
the reduction in pH that occurs as the two-phase<br />
boundary of seawater is approached, but is <strong>al</strong>so<br />
due to the enhanced stability of aqueous Fe-chlorocomplexes<br />
[19,20].<br />
Recent an<strong>al</strong>yses of fluids from modern vent systems<br />
further indicate that phase separation in response<br />
to pressure and/or temperature change can result in<br />
significant corresponding changes in dissolved Fe<br />
concentrations. For example, in spite of the low s<strong>al</strong>inity<br />
of the vapor phase fluids released from subseafloor<br />
reaction zones undergoing phase separation in<br />
connection with subseafloor magmatic activity at EPR<br />
9–108 N, Fe concentrations are unusu<strong>al</strong>ly high on an<br />
absolute and chloride norm<strong>al</strong>ized basis [11]. The high<br />
Fe concentrations of the vapors likely result from<br />
aqueous complexing effects, as noted earlier.<br />
Fe H 2 S Fe/H 2 S<br />
Fig. 4. The effect of pressure on dissolved Fe, H 2 S and Fe/H 2 S for a<br />
chemic<strong>al</strong> system buffered by fay<strong>al</strong>ite–quartz–magn<strong>et</strong>ite–pyrrhotite<br />
coexisting with 0.55 mol/kg chloride at a pH=5 and 400 8C (see<br />
text). The increase in dissolved Fe and Fe/H 2 S ratio predicted for<br />
the 100 bar pressure drop likely continues at still lower pressures<br />
owing to the effect of pressure on the stability of aqueous ferrous<br />
chloride complexes and miner<strong>al</strong> phase relations, as suggested by<br />
earlier experiment<strong>al</strong> data [20]. These data suggest that the decrease<br />
in hydrostatic pressure associated with glob<strong>al</strong> glaciations in the<br />
Neoproterozoic may enhance the flux of hydrotherm<strong>al</strong> Fe to the<br />
ancient ocean. Relatively sh<strong>al</strong>low ridges in the Archean tog<strong>et</strong>her<br />
with a sulfate-free ocean would similarly enhance Fe flux for BIF<br />
deposition at that time. Pressure sufficiently low that the two-phase<br />
boundary of the fluid (ancient or modern) is intersected would result<br />
in formation of Fe-bearing vapors and brines. Assuming both phases<br />
ultimately reach the seafloor, the Fe flux would be significant.<br />
Thermodynamic data sources are as discussed in earlier figures.<br />
The coexisting brine phase at depth is likely<br />
exceedingly rich in Fe, owing to its high chloride<br />
concentration and relatively low pH [8,10,21]. Although<br />
it is commonly believed that vapor loss<br />
causes an increase in the pH of the residu<strong>al</strong> brine<br />
phase [22], this is not the case, however, in the<br />
presence of silicate miner<strong>al</strong>s that can buffer pH at<br />
relatively low v<strong>al</strong>ues in response to the increase in<br />
chloride and charge b<strong>al</strong>ance constraints [23], especi<strong>al</strong>ly<br />
at high temperatures where such reactions<br />
occur rapidly. The Fe-bearing brines likely reside<br />
in the crust for some time, but ultimately mix with<br />
hydrotherm<strong>al</strong> seawater and vent as well, as indicated,<br />
for example, by the composition of vent fluids at the<br />
southern Cleft Segment of the Juan de Fuca Ridge<br />
[24–26]. If phase separation/segregation takes place<br />
very near the seafloor, as has been suggested for vent
L.R. <strong>Kump</strong>, W.E. Seyfried Jr. / Earth and Plan<strong>et</strong>ary Science L<strong>et</strong>ters 235 (2005) 654–662 659<br />
fluids from the Southern East Pacific Rise and East<br />
Pacific Rise at 9–108 N, venting of Fe-bearing brine<br />
is even more likely, especi<strong>al</strong>ly for conditions at or<br />
near the critic<strong>al</strong> point of the NaCl source fluid [11].<br />
Similar arguments apply to Precambrian vent systems,<br />
<strong>al</strong>though even higher Fe concentrations are<br />
likely owing to the lack of sulfate (more reducing<br />
conditions) in the hydrotherm<strong>al</strong> root zones from<br />
which the fluids are derived.<br />
3. Precambrian mid-ocean ridge systems<br />
3.1. Ridge depressurization and sulfate depl<strong>et</strong>ion during<br />
Snowb<strong>al</strong>l Earth<br />
Might the pressure decrease associated with the<br />
growth of continent<strong>al</strong> ice she<strong>et</strong>s cause an enhancement<br />
of the mid-ocean ridge Fe flux sufficient to<br />
stimulate BIF deposition? As emphasized above,<br />
redox constrains imposed by the lack of sulfate in<br />
Neoproterozoic seawater during glob<strong>al</strong> glaci<strong>al</strong> events<br />
likely play a key role in accounting for BIF deposition.<br />
Nevertheless, pressure changes associated with<br />
glaciation could have significant secondary effects on<br />
Fe flux. All of the Neoproterozoic iron formations are<br />
in glaciomarine sequences with associated glaci<strong>al</strong> diamictites<br />
[27]. Thus the coincidence of glaciation and<br />
iron deposition is on secure footing. Mid-ocean ridge<br />
fluxes represent a significant fraction of the tot<strong>al</strong><br />
delivery of soluble Fe to the ocean today [1,28,29],<br />
and on a Snowb<strong>al</strong>l Earth, would be the only major<br />
source of Fe to the ocean.<br />
Estimates of Neoproterozoic sea-level fluctuations<br />
range from z160 m [30] to more than 1000 m [31].<br />
At the low end of these estimates (equiv<strong>al</strong>ent to a drop<br />
in hydrostatic pressure of 16 bars), the enhancement<br />
of Fe fluxes would be minor. At the upper end (a 100<br />
bar depressurization), and with fluids in the seafloor<br />
reaction zone equilibrated at temperatures greater than<br />
400 8C, the enhancement of mid-ocean ridge Fe fluxes<br />
could be appreciable.<br />
Interestingly, <strong>al</strong>though the low-end estimated increase<br />
in Fe flux is sm<strong>al</strong>l for one-phase systems,<br />
even modest changes in pressure could enhance Fe<br />
delivery for systems sufficiently close to the critic<strong>al</strong><br />
point of the NaCl fluid, and this argument<br />
applies equ<strong>al</strong>ly well to the Quaternary. Even a<br />
relatively sm<strong>al</strong>l increase in Fe delivery, coupled<br />
with sh<strong>al</strong>lower ridge crests and perhaps more vigorous<br />
hydrotherm<strong>al</strong> circulation at near critic<strong>al</strong> conditions<br />
could presumably inject Fe-rich plumes<br />
more efficiently into the thermocline, stimulating<br />
glaci<strong>al</strong> biologic<strong>al</strong> productivity and a draw-down<br />
of atmospheric pCO 2 comparable to that estimated<br />
to result from intensified aerosol fluxes during the<br />
Quaternary [32].<br />
Recognizing that there were at least two and<br />
perhaps sever<strong>al</strong> Neoproterozoic glaciations [33,34],<br />
one might expect to find a correlation b<strong>et</strong>ween the<br />
extent of continent<strong>al</strong> glaciation (inferred magnitude<br />
of sea-level f<strong>al</strong>l) and abundance of iron formation. Of<br />
the 8 glaci<strong>al</strong>ly associated iron formations listed by<br />
Klein and Beukes [35], 6 are apparently bSturtianQ in<br />
age (~750 Ma; Kingston Peak, Death V<strong>al</strong>ley, C<strong>al</strong>ifornia;<br />
Rapitan, Mackenzie Mountains, Canada;<br />
Upper Tindir Group, Alaska; Adelaide Geosyncline,<br />
South Austr<strong>al</strong>ia; Numees Formation, South Africa<br />
and Namibia; and Chuos Formation, Namibia) and<br />
2 are bVarangian or MarinoanQ (~600 Ma; Bissokpabe<br />
Group, West Africa; and Jacadigo Group, Brazil)<br />
[36]. Thus, iron deposition was concentrated in<br />
the older glaci<strong>al</strong> interv<strong>al</strong>. But did the Sturtian event<br />
produce thicker and/or more extensive continent<strong>al</strong> ice<br />
she<strong>et</strong>s and thus a larger sea-level f<strong>al</strong>l? Unfortunately,<br />
the magnitude of the Neoproterozoic sea-level fluctuations<br />
remains controversi<strong>al</strong>. Thus, we are left with<br />
an untested hypothesis that the Sturtian glaciation<br />
involved larger continent<strong>al</strong> ice she<strong>et</strong>s than the Marinoan<br />
event.<br />
An <strong>al</strong>ternative explanation for the preponderance<br />
of BIF in the Sturtian glaci<strong>al</strong> deposits is that sulfate<br />
concentrations were lower during the Sturtian than<br />
during the Marinoan event. This explanation was<br />
forwarded by Hurtgen <strong>et</strong> <strong>al</strong>. [13], based on their<br />
an<strong>al</strong>ysis of the trace sulfate content and sulfur isotopic<br />
composition of cap carbonates deposited following<br />
these events. Of course, a more intense glaciation<br />
(presumably the Sturtian) would more effectively<br />
eliminate the weathering supply of sulfate. Addition<strong>al</strong>ly,<br />
the isolation of the ocean covered in sea ice<br />
would favor the consumption of what little sulfate<br />
initi<strong>al</strong>ly existed in the ocean by bacteri<strong>al</strong> sulfate reduction<br />
or hydrotherm<strong>al</strong> uptake, leading to a large<br />
increase in the Fe/S ratio and the absolute flux of<br />
Fe from mid-ocean ridges.
660<br />
L.R. <strong>Kump</strong>, W.E. Seyfried Jr. / Earth and Plan<strong>et</strong>ary Science L<strong>et</strong>ters 235 (2005) 654–662<br />
3.2. Archean and P<strong>al</strong>eoproterozoic BIF production<br />
The abundance of BIF in older (Archean and<br />
P<strong>al</strong>eoproterozoic) sedimentary rocks could <strong>al</strong>so be<br />
the result of a predominance of low oceanic sulfate<br />
concentrations and lower pressure hydrotherm<strong>al</strong> systems.<br />
The observation that Archean sh<strong>al</strong>es are<br />
enriched in Fe has been interpr<strong>et</strong>ed to indicate a<br />
larger hydrotherm<strong>al</strong> flux of Fe to the early ocean<br />
compared to today [37]. Isley [1] implicated higher<br />
heat flow (implying increased hydrotherm<strong>al</strong> circulation)<br />
and sh<strong>al</strong>lower ridges to account for this enhanced<br />
Fe delivery; sh<strong>al</strong>lower ridges would have<br />
facilitated the delivery of mid-ocean ridge Fe in<br />
hydrotherm<strong>al</strong> plumes to surface waters leading to<br />
BIF deposition. Upwelling of Fe-rich water to the<br />
surface may have supported anoxygenic photosynthesis<br />
based on Fe oxidation [38,39], or the Fe may<br />
have been oxidized anoxic<strong>al</strong>ly and inorganic<strong>al</strong>ly. The<br />
predictably high heat flow in the Archean and accelerated<br />
rates of hydrotherm<strong>al</strong> circulation at ridges,<br />
notwithstanding, Lowell and Keller [40] have<br />
shown that simple extrapolation of the modern hydrotherm<strong>al</strong><br />
Fe flux to the Archean using crust<strong>al</strong><br />
cooling models underestimates Fe in the large superior-type<br />
BIF by approximately an order of magnitude.<br />
This inconsistency, however, can be resolved in<br />
large part by taking explicit account of the absence<br />
of sulfate in the Archean ocean on hydrotherm<strong>al</strong><br />
<strong>al</strong>teration processes (Figs. 2 and 3).<br />
Archean vent fluids would have been relatively<br />
sulfide poor, but H 2 -rich, with H 2 concentrations predictably<br />
larger than modern vent fluids. Thus, the<br />
biota they supported likely was based on H 2 chemosynthesis<br />
(m<strong>et</strong>hanogenesis) than on sulfide oxidation<br />
as a chemosynth<strong>et</strong>ic mechanism, more prev<strong>al</strong>ent<br />
today.<br />
The end of Archean–P<strong>al</strong>eoproterozoic BIF deposition<br />
may have resulted from decreasing heat flow,<br />
deeper ridges [1], and perhaps most importantly, the<br />
transition to more oxidizing sulfate-bearing systems.<br />
Such a mechanism is not inconsistent with the suggestion<br />
that the increase in sulfate concentrations in<br />
the P<strong>al</strong>eoproterozoic ocean stimulated sulfate reduction<br />
and created a sulfidic deep ocean [5]. Such an<br />
ocean could conceivably have elevated concentrations<br />
of both sulfate and sulfide compared to those of the<br />
Archean ocean, but would be ineffective as a transfer<br />
agent for BIF deposition at points dist<strong>al</strong> from source<br />
regions at ridges.<br />
4. Conclusion<br />
Because they operate close to the critic<strong>al</strong> point and<br />
their redox state is s<strong>et</strong> largely by the availability of<br />
dissolved sulfate, the chemic<strong>al</strong> composition of midocean<br />
ridge vent fluids is particularly sensitive to<br />
changes in temperature, pressure, and seawater sulfate<br />
concentration. In particular, the Fe content and Fe/<br />
H 2 S ratio of mid-ocean ridge hydrotherm<strong>al</strong> fluids<br />
should have been significantly elevated if Precambrian<br />
seawater was depl<strong>et</strong>ed in sulfate and if ridges were<br />
sh<strong>al</strong>low during the Archean and during Snowb<strong>al</strong>l<br />
Earth episodes of eustatic sea-level draw-down. The<br />
abundance of BIF under these conditions would then<br />
have been largely the consequence of changes in<br />
hydrotherm<strong>al</strong> fluid composition rather than changes<br />
in the solubility of Fe in seawater d<strong>et</strong>ermined by the<br />
O 2 [28] or H 2 S content [5] of the deep sea.<br />
Acknowledgments<br />
This work was supported by the NASA Astrobiology<br />
Institute (Cooperative Agreements NCC2-0157<br />
and NNA04CC06A to LRK) and the Nation<strong>al</strong> Science<br />
Foundation. P. Hoffman offered v<strong>al</strong>uable<br />
insights and specific information on the geologic<strong>al</strong><br />
evidence for glacioeustatic sea-level f<strong>al</strong>l during the<br />
Neoproterozoic. H.D. Holland and N. Sleep provided<br />
careful and useful reviews of earlier versions of the<br />
manuscript.<br />
References<br />
[1] A.E. Isley, Hydrotherm<strong>al</strong> plumes and the delivery of iron to<br />
banded iron formation, J. Geol. 103 (1994) 169–185.<br />
[2] D.R. Janecky, W.E. Seyfried Jr., Formation of massive sulfide<br />
deposits on oceanic ridge crests: increment<strong>al</strong> reaction models<br />
for mixing b<strong>et</strong>ween hydrotherm<strong>al</strong> solutions and seawater,<br />
Geochim. Cosmochim. Acta 48 (1984) 2723–2738.<br />
[3] K.L. Von Damm, Controls on the chemistry and tempor<strong>al</strong><br />
variability of seafloor hydrotherm<strong>al</strong> fluids, in: S.E. Humphris,<br />
R.A. Zierenberg, L.S. Mullineaux, R.E. Thomson (Eds.),<br />
Seafloor Hydrotherm<strong>al</strong> Systems: Physic<strong>al</strong>, Chemic<strong>al</strong>, Biologic<strong>al</strong>,<br />
and Geologic<strong>al</strong> Interactions, Geophysic<strong>al</strong> Monograph,
L.R. <strong>Kump</strong>, W.E. Seyfried Jr. / Earth and Plan<strong>et</strong>ary Science L<strong>et</strong>ters 235 (2005) 654–662 661<br />
vol. 91. American Geophysic<strong>al</strong> Union, Washington, D.C.,<br />
1995, pp. 222–247.<br />
[4] J.C.G. W<strong>al</strong>ker, P. Brimblecombe, Iron and sulfur in the prebiologic<br />
ocean, Precambrian Res. 28 (1985) 205–222.<br />
[5] D.E. Canfield, A new model for Proterozoic ocean chemistry,<br />
Nature 396 (1998) 450–453.<br />
[6] J.L. Kirschvink, Late Proterozoic glob<strong>al</strong> glaciation: the Snowb<strong>al</strong>l<br />
Earth, in: J.W. Schopf, C. Klein (Eds.), The Proterozoic<br />
Biosphere: A Multidisciplinary Study, Cambridge University<br />
Press, Cambridge, UK, 1992, pp. 51–52.<br />
[7] N.J. Beukes, C. Klein, Models for iron-formation deposition,<br />
in: J.W. Schopf, C. Klein (Eds.), The Proterozoic Biosphere: A<br />
Multidisciplinary Study, Cambridge University Press, Cambridge,<br />
UK, 1992, pp. 147–151.<br />
[8] W.E. Seyfried Jr., K. Ding, Phase equilibria in subseafloor<br />
hydrotherm<strong>al</strong> systems; a review of the role of redox, temperature,<br />
pH and dissolved Cl on the chemistry of hot spring<br />
fluids at mid-ocean ridges, in: S.E. Humphris , R.A. Zierenberg,<br />
L.S. Mullineaux, R.E. Thomson (Eds.), Seafloor Hydrotherm<strong>al</strong><br />
Systems: Physic<strong>al</strong>, Chemic<strong>al</strong>, Biologic<strong>al</strong>, and<br />
Geologic<strong>al</strong> Interactions. Geophysysic<strong>al</strong> Monograph, vol. 91,<br />
Amer. Geop. Union, Washington, D.C., 1995, pp. 248–272.<br />
[9] K.L. Von Damm, Chemistry of hydrotherm<strong>al</strong> vent fluids from<br />
98–108 N, East Pacific Rise: btime zero,Q the immediate posteruptive<br />
period, J. Geophys. Res. 105 (2000) 11,203–11,222.<br />
[10] K.L. Von Damm, L.G. Buttermore, S.E. Oosting, A.M. Bray,<br />
D.J. Fornari, M.D. Lilley, W.C. Shanks III, Direct observation<br />
of the evolution of a seafloor bblack smokerQ from vapor to<br />
brine, Earth Plan<strong>et</strong>. Sci. L<strong>et</strong>t. 149 (1997) 101–111.<br />
[11] K.L. Von Damm, M.D. Lilley, W.C. Shanks, M. Brockington,<br />
A. Bray, K.M. O’Grady, E. Olson, A. Graham, G. Proskurowski,<br />
Extraordinary phase separation and segregation in vent<br />
fluids from the southern East Pacific Rise, Earth Plan<strong>et</strong>. Sci.<br />
L<strong>et</strong>t. 206 (2003) 365–378.<br />
[12] A.D. Anbar, A.H. Knoll, Proterozoic ocean chemistry and<br />
evolution; a bioinorganic bridge? Science 297 (2002)<br />
1137–1142.<br />
[13] M.T. Hurtgen, M.A. Arthur, N.S. Suits, A.J. Kaufman, The<br />
sulfur isotopic composition of Neoproterozoic seawater sulfate:<br />
implications for a Snowb<strong>al</strong>l Earth? Earth Plan<strong>et</strong>. Sci. L<strong>et</strong>t.<br />
203 (2002) 413–429.<br />
[14] M.D. Lilley, J.E. Lupton, D.A. Butterfield, E. Olson, Magmatic<br />
events produce rapid changes in hydrotherm<strong>al</strong> vent<br />
chemistry, Nature 422 (2003) 878–881.<br />
[15] K.L. Von Damm, Evolution of the hydrotherm<strong>al</strong> system at East<br />
Pacific Rise 9850V N: geochemic<strong>al</strong> evidence for changes in the<br />
upper oceanic crust, in: C.R. German, J. Linn, L.M. Parson<br />
(Eds.), Mid-Ocean Ridges: Hydrotherm<strong>al</strong> Interactions B<strong>et</strong>ween<br />
the Lithosphere and Oceans, American Geophysic<strong>al</strong><br />
Union, Washington, D.C., 2004, pp. 285–305.<br />
[16] M.K. Tivey, S.E. Humphris, G. Thompson, M.D. Hannington,<br />
P.A. Rona, Deducing patterns of fluid flow and mixing within<br />
the TAG active hydrotherm<strong>al</strong> mound using miner<strong>al</strong>ogic<strong>al</strong> and<br />
geochemic<strong>al</strong> data, J. Geophys. Res., Solid Earth 100 (B7)<br />
(1995) 12527–12555.<br />
[17] W.C. Shanks III, Stable isotopes in seafloor hydrotherm<strong>al</strong><br />
systems: vent fluids, hydrotherm<strong>al</strong> deposits, hydrotherm<strong>al</strong> <strong>al</strong>teration,<br />
and microbi<strong>al</strong> processes, in: J.W. V<strong>al</strong>ley , D.R. Cole<br />
(Eds.), Stable Isotope Geochemistry, Miner<strong>al</strong>ogic<strong>al</strong> Soci<strong>et</strong>y of<br />
America, Washington, D.C., 2001, pp. 469–526.<br />
[18] J.C. Alt, Sulfur isotopic profiles through the oceanic crust:<br />
sulfur mobility and seawater–crust<strong>al</strong> sulfur exchange during<br />
hydrotherm<strong>al</strong> <strong>al</strong>teration, Geology 23 (1995) 585–588.<br />
[19] K. Ding, W.E. Seyfried Jr., D<strong>et</strong>ermination of Fe–Cl complexing<br />
in the low-pressure supercritic<strong>al</strong> region (NaCl fluid)—iron<br />
solubility constraints on pH of subseafloor hydrotherm<strong>al</strong><br />
fluids, Geochim. Cosmochim. Acta 56 (1992) 3681–3692.<br />
[20] W.E. Seyfried Jr., D.R. Janecky, Heavy m<strong>et</strong><strong>al</strong> and sulfur<br />
transport during subcritic<strong>al</strong> and supercritic<strong>al</strong> hydrotherm<strong>al</strong><br />
<strong>al</strong>teration of bas<strong>al</strong>t: influence of fluid pressure and bas<strong>al</strong>t<br />
composition and cryst<strong>al</strong>linity, Geochim. Cosmochim. Acta<br />
49 (1985) 2245–2260.<br />
[21] D. Butterfield, I.R. Jonasson, G.J. Massoth, R.A. Feely, K.K.<br />
Roe, R.W. Embley, J.F. Holden, J.R. McDuff, M.D. Lilley,<br />
J.R. Delaney, Seafloor eruptions and evolution of hydrotherm<strong>al</strong><br />
fluid chemistry, Philos. Trans. R. Soc. Lond., A 355<br />
(1997) 369–386.<br />
[22] S.E. Drummond, H. Ohmoto, Chemic<strong>al</strong> evolution and miner<strong>al</strong><br />
deposition in boiling hydrotherm<strong>al</strong> systems, Econ. Geol. 80<br />
(1986) 126–147.<br />
[23] D. Butterfield, W.E. Seyfried Jr., M. Lilley, Composition and<br />
evolution of hydrotherm<strong>al</strong> fluids, in: P. H<strong>al</strong>bach , V. Tunnicliffe,<br />
J.R. Hein (Eds.), Energy and Mass Transfer in Marine<br />
Hydrotherm<strong>al</strong> Systems, Dahlem University Press, Berlin, Germany,<br />
2003, pp. 123–163.<br />
[24] K.L. Von Damm, J.L. Bischoff, Chemistry of hydrotherm<strong>al</strong><br />
solutions from the Southern Juan de Fuca Ridge, J. Geophys.<br />
Res. 92 (1987) 11334–11346.<br />
[25] D.A. Butterfield, G.J. Massoth, Geochemistry of the North<br />
Cleft segment vent fluids: tempor<strong>al</strong> changes in chlorinity and<br />
their possible relation to recent volcanism, J. Geophys. Res. 99<br />
(1994) 4951–4968.<br />
[26] S. Schoofs, U. Hansen, Depl<strong>et</strong>ion of a brine layer at the base of<br />
ridge–crest hydrotherm<strong>al</strong> systems, Earth Plan<strong>et</strong>. Sci. L<strong>et</strong>t. 180<br />
(2002) 341–353.<br />
[27] M.J. Hambrey, W.B. Harland, Earth’s Pre-Pleistocene Glaci<strong>al</strong><br />
Record, Cambridge University Press, Cambridge, UK, 1981,<br />
1004 pp.<br />
[28] H.D. Holland, The Chemic<strong>al</strong> Evolution of the Atmosphere<br />
and Oceans, Princ<strong>et</strong>on University Press, Princ<strong>et</strong>on, NJ, 1984,<br />
582 pp.<br />
[29] A.J. Watson, D.C.E. Bakker, A.J. Ridgwell, P.W. Boyd, C.S.<br />
Law, Effect of iron supply on southern ocean CO 2 uptake and<br />
implications for glaci<strong>al</strong> atmospheric CO 2 , Nature 407 (2000)<br />
730–733.<br />
[30] N. Christie-Blick, L.E. Sohl, M.J. Kennedy, Considering a<br />
Neoproterozoic Snowb<strong>al</strong>l Earth, Science 284 (1999) 1087.<br />
[31] W.T. Hyde, T.J. Crowley, S.K. Baum, W.R. Peltier, Neoproterozoic<br />
bSnowb<strong>al</strong>l EarthQ simulations with a coupled climate/<br />
ice-she<strong>et</strong> model, Nature 205 (2000) 425–429.<br />
[32] N. Lefevre, A.J. Watson, Modelling the geochemic<strong>al</strong> cycle of<br />
iron in the oceans and its impact on atmospheric carbon<br />
dioxide concentrations, Glob. Biogeochem. Cycles 13 (1999)<br />
727–736.
662<br />
L.R. <strong>Kump</strong>, W.E. Seyfried Jr. / Earth and Plan<strong>et</strong>ary Science L<strong>et</strong>ters 235 (2005) 654–662<br />
[33] P.F. Hoffman, A.J. Kaufman, G.P. H<strong>al</strong>verson, D.P. Schrag,<br />
A Neoproterozoic Snowb<strong>al</strong>l Earth, Science 281 (1998)<br />
1342–1346.<br />
[34] M.J. Kennedy, B. Runnegar, A.R. Prave, K.-H. Hoffmann,<br />
M.A. Arthur, Two or four Neoproterozoic glaciations? Geology<br />
26 (1998) 1059–1063.<br />
[35] C. Klein, N.J. Beukes, Sedimentology and geochemistry of the<br />
glaciogenic Late Proterozoic Rapitan iron-formation in Canada,<br />
Econ. Geol. 88 (1993) 542–565.<br />
[36] M.J. Hambrey, W.B. Harland, The Late Proterozoic glaci<strong>al</strong><br />
era, P<strong>al</strong>aeogeogr. P<strong>al</strong>aeoclim. P<strong>al</strong>aeoecol. 51 (1985)<br />
255–272.<br />
[37] L.R. <strong>Kump</strong>, H.D. Holland, Iron in Precambrian rocks: implications<br />
for the glob<strong>al</strong> oxygen budg<strong>et</strong> of the ancient Earth,<br />
Geochim. Cosmochim. Acta 56 (1992) 3217–3223.<br />
[38] F. Widdel, S. Schnell, S. Heising, A. Ehrenreich, B. Assmus,<br />
B. Schink, Ferrous iron oxidation by anoxygenic phototrophic<br />
bacteria, Nature 362 (1993) 834–836.<br />
[39] L.R. <strong>Kump</strong>, Bacteria forge a new link, Nature 362 (1993)<br />
790–791.<br />
[40] R.P. Lowell, S.M. Keller, High-temperature seafloor hydrotherm<strong>al</strong><br />
circulation over geologic time and Archean banded<br />
iron formations, Geophys. Res. L<strong>et</strong>t. 30 (2003) 1391,<br />
doi:10.1029/2002GL016536.<br />
[41] J.W. Johnson, E.H. Oelkers, H.C. Helgeson, SUPCRT92—a<br />
software package for c<strong>al</strong>culating the standard mol<strong>al</strong> thermodynamic<br />
properties of miner<strong>al</strong>s, gases, aqueous species, and<br />
reactions from 1-bar to 5000-bar and 0 8C to 1000 8C, Comput.<br />
Geosci. 18 (1992) 899–947.<br />
[42] P.B. Barton, B.J. Skinner, Sulfide miner<strong>al</strong> stabilities, in: H.L.<br />
Barnes (Ed.), Geochemistry of Hydrotherm<strong>al</strong> Ore Deposits,<br />
John Wiley & Sons, New York, 1979, pp. 278–390.<br />
[43] N. Kishima, A thermodynamic study on the pyrite–pyrrhotite–<br />
magn<strong>et</strong>ite–water system at 300–500 8C with relevance to the<br />
fugacity/concentration quotient of aqueous H 2 S, Geochim.<br />
Cosmochim. Acta 53 (1989) 2143–2155.<br />
[44] K. Ding, W.E.J. Seyfried Jr., Activity coefficients of H 2 and<br />
H 2 S in NaCl solutions at 300–425 8C, 300–500 bars with<br />
application to ridge crest hydrotherm<strong>al</strong> systems, EOS, Trans.<br />
Amer. Geophys. Union 71 (1990) 1680.<br />
[45] T.J. Wolery, S.A. Daveler, EQ6, a Computer Program for<br />
Reaction Path Modeling of Aqueous Geochemic<strong>al</strong> Systems:<br />
Theor<strong>et</strong>ic<strong>al</strong> Manu<strong>al</strong>, Users Guide, and Related Documentation<br />
(version 7.0), 1992, pp. 1 – 337. UCRL-MA-110662 PT IV.<br />
[46] D.A. Sverjensky, E.L. Shock, H.C. Helgeson, Prediction of the<br />
thermodynamic properties of aqueous m<strong>et</strong><strong>al</strong> complexes to<br />
1000 8C and 5 kb, Geochim. Cosmochim. Acta 61 (1997)<br />
1359–1412.