Molybdenum burial mechanism in sulfidic sediments: Iron-Sulfide

Mar 20, 2018 - Trent Vorlicek , George Rudolph Helz , Anthony Chappaz , Pakou Vue , Austin Vezina , and Wayland Hunter. ACS Earth Space Chem. , Just ...
0 downloads 0 Views 5MB Size
Subscriber access provided by Queen Mary, University of London

Molybdenum burial mechanism in sulfidic sediments: Iron-Sulfide pathway Trent Vorlicek, George Rudolph Helz, Anthony Chappaz, Pakou Vue, Austin Vezina, and Wayland Hunter ACS Earth Space Chem., Just Accepted Manuscript • DOI: 10.1021/ acsearthspacechem.8b00016 • Publication Date (Web): 20 Mar 2018 Downloaded from http://pubs.acs.org on March 25, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 52 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

1 2 3 4 5

Molybdenum burial mechanism in sulfidic sediments: Iron-Sulfide pathway

6 7

Trent P. Vorlicek1*, George R. Helz2, Anthony Chappaz3, Pakou Vue1, Austin Vezina1, Wayland

8

Hunter1

9 10 11 12 13 14 15 16 17 18 19 20

1

Department of Chemistry and Geology, Minnesota State University, Mankato, MN 56001, USA

2

Department of Chemistry and Biochemistry, University of Maryland, College Park, MD 20742, USA 3

Department of Earth and Atmospheric Sciences, Central Michigan University, Mount Pleasant, MI 48859, USA

21 22 23 24 25 26 27

* Corresponding author. Tel.: 507-389-1598

28

E-mail address: [email protected]

29 30 31 32 33 34 35 36 1 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 52

37

ABSTRACT: Relative to continental crust, sediments underlying sulfidic marine waters are

38

molybdenum-rich, a property preserved in the rock record and useful for characterizing

39

paleoenvironments. The enrichment mechanism is not agreed upon, but is attributed at least partly

40

to deposition of Fe-Mo-S compounds, which are as yet uncharacterized. Here, we determine the

41

composition and stability of colloidal Fe-Mo-S precipitates formed at mildly basic pH and H2S(aq)

42

> 10-5 M. The first product consists simply of FeMoS4, with Ksp = 10-14.95. Within hours, FeMoS4

43

irreversibly transforms by internal self-reduction to a Mo(IV) product of similar composition. The

44

reduced product is insoluble in 1 M HCl but soluble in concentrated HNO3, implying that it would

45

be recovered with pyrite in a common assay of sediments. XAFS data show that Mo(IV) in the

46

colloids is coordinated by a split first shell of about five sulfur atoms at average distances of 2.31

47

and 2.46 Å and in its second shell by an iron atom at about 2.80 Å. These properties resemble those

48

determined for Mo in modern anoxic lake sediments and in Phanerozoic black shales. The atomic

49

environment around Mo suggests that the colloidal products may be inorganic polymers containing

50

cuboid, Fe2Mo2S44+ cores. Such materials are so far unreported by mineralogists, although a rare

51

mineral, jordisite, may be a related, but more Mo-rich material. The low solubility of FeMoS4

52

makes it a feasible precipitate in euxinic waters like those in the modern Black Sea. We propose

53

that colloids similar to those studied here could account for Mo-enrichment in euxinic basin

54

sediments and black shales.

55 56

Keywords: Molybdenum, euxinic basins, iron sulfide, jordisite, colloids

57 58

2 ACS Paragon Plus Environment

Page 3 of 52 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

59 60

ACS Earth and Space Chemistry

1. INTRODUCTION Molybdenum (Mo) concentrations in ancient black shales and modern sediments shed light

61

on the redox state of Earth’s ocean-atmosphere system and on recent anoxic events in marine

62

basins.1,2 Molybdenum serves as a paleo-redox proxy for past sulfidic conditions because its

63

concentrations in black shales increase from detritic background (1-2 ppm) to exceptionally high

64

values (up to 200-300 ppm) when redox conditions shift from oxic to euxinic.3,4 Two pathways are

65

being considered to explain Mo sequestration under euxinic conditions:5,6 (1) the iron-sulfide

66

pathway7-10 and (2) organic matter pathway.6,11,12 Quantitative tests of these pathways are not

67

possible at present. The goal of this paper is to put the first on a quantitative footing.

68

The idea that Mo in euxinic basins follows an iron-sulfide pathway has a long history.13-17

69

Modern support for this pathway comes from X-ray Absorption Fine Structure (XAFS) studies: S

70

atoms are found in the first coordination shell and Fe atoms in the second coordination shell around

71

Mo atoms in some black shales and anoxic lake sediments.7,18,19 Based on field data from

72

Rogoznica Lake,8 a quantitative thermodynamic model of Mo deposition in sulfidic waters via an

73

iron sulfide pathway was proposed. The model postulated existence of an unknown mineral with an

74

approximate, empirical formula, FeMo0.6S2.8. This model predicted water-column Mo

75

concentrations in a half-dozen euxinic water bodies to within about a factor of three. In post facto

76

tests, it predicted Mo concentrations in the sulfidic water column of Green Lake and in the sulfidic

77

pore waters of Terrebonne Bay sediments to a similar level of accuracy.20,21 Nevertheless, a

78

weakness of this model is that compounds approaching the proposed mineral composition are

79

unknown to both mineralogists and inorganic chemists. Possibly this material has evaded discovery

80

in nature because its structure is unsuited to ordered crystallization, and therefore it is necessarily

81

colloidal and x-ray amorphous.8

82

Evidence has accumulated over several decades that suspended FeS particles, often of

83

colloidal size, are commonly found in sulfidic waters.22-30 In this paper, we determine properties of

3 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 84 3 4 85 5 6 86 7 8 9 87 10 11 88 12 13 89 14 15 16 90 17 18 91 19 20 92 21 22 23 93 24 25 94 26 27 95 28 29 96 30 31 32 97 33 34 98 35 36 99 37 38 39 100 40 41 101 42 43 102 44 45 46 103 47 48 104 49 50 105 51 52 106 53 54 55 107 56 57 108 58 59 60

Page 4 of 52

Fe-Mo-S precipitates formed by the reaction of colloidal FeS with thiomolybdate, and we propose that products of this reaction are significant Mo hosts in sediments of euxinic basins. 2. METHODS Materials. The following reagents were purchased and used as received: Na2B4O7∙10H20 (Sigma-Aldrich), orthorhombic sulfur (Sigma-Aldrich), Fe(NH4)2(SO4)2∙6H2O (Mallinckrodt), MgCl2∙6H2O (BDH), HgCl2 (J.T. Baker), Na2O2 (Merck), BaCl2∙2H2O (Fisher), Na2SO4 (Mallinckrodt), Na2CO3 (Fisher), 0.1000(±0.0001) M EDTA (Alfa-Aesar), LC-MS grade acetonitrile (EMS), and tetrabutylammonium hydroxide thirty hydrate (Sigma). Na2S∙9H2O (Sigma-Aldrich) crystals were individually rinsed with H2O until any yellowish oxidative impurities were no longer visible. Several authors have noted that commercial (NH4)2MoS4 reagents often contain substantial amounts of an unidentified impurity, giving them a mauve-brown color rather than their normal deep red color, which is like that of a full-bodied red wine.31,32 The specific impurity has not been identified, but it may be a dimeric Mo(V)-polysulfido compound (e.g., (NH4)2Mo2O2(μ-S)2(S2)2), arising from oxidation over long storage times.33 Consequently, we opted to avoid commercial thiomolybdate reagents, instead synthesizing fresh (NH4)2MoS4 by the reaction of MoO42- with sulfide. Test Solution Procedures. Solutions were prepared at room temperature using deoxygenated ultrapure H2O in an anaerobic chamber under a 95% N2:5% H2 atmosphere. Na2B4O7∙10H2O and Na2S∙9H2O were placed into a 1 L beaker. Following addition of ~800 mL of H2O, the pH was adjusted, and an aliquot of MoS42- stock solution previously standardized by UVvis absorbance34 was added. After diluting to 1L, the required volume of Fe2+ stock solution was added, resulting in formation of a stable, inky-black, colloidal suspension. Fourteen mL of this suspension was apportioned into each of a series of 15 mL plastic centrifuge tubes that were then capped and stored in the glove box. Upon aging ~100 hours, centrifuge tubes were opened within the glovebox and 350 μL of 4 M MgCl2 was added to promote flocculation of the colloids; tubes

4 ACS Paragon Plus Environment

Page 5 of 52 1 2 109 3 4 110 5 6 111 7 8 9 112 10 11 113 12 13 114 14 15 16 115 17 18 116 19 20 117 21 22 23 118 24 25 119 26 27 120 28 29 121 30 31 32 122 33 34 123 35 36 124 37 38 39 125 40 41 126 42 43 127 44 45 46 128 47 48 129 49 50 130 51 52 131 53 54 55 132 56 57 133 58 59 60

ACS Earth and Space Chemistry

were again capped and stored in the glovebox. Flocculation was unsatisfactory unless MgCl2 was added after some aging. Analyses. At various times, individual centrifuge tubes were opened and their contents filtered through Whatman 0.02 μm syringe filters. Filtered aliquots were immediately passed out of the anaerobic chamber for UV-vis analyses (Schimadzu UV-2450 double-beam spectrophotometer). We calculated MoS42- concentrations using extinction coefficients,35,36 attributing all absorption at 468 nm to MoS42-. Samples were treated with BrCl to oxidize sulfide and with aluminum (1000 ppm) as a matrix modifier prior to determination of total dissolved Mo by atomic absorption spectroscopy (Perkin-Elmer PinAAcle 900T spectrophotometer). At the end of a run, three samples were analyzed separately to obtain average total sulfide and pH. Total sulfide was titrated with standardized Hg2+ to a potentiometric endpoint (Ag/Ag2S vs. double junction reference electrodes) after adjusting the solution to pH ~13. This titration quantifies the sum of H2S and HS- plus any sulfide ligands associated with thiomolybdates and dissolved Fe-Mo-S complexes. Free sulfide (H2S+HS-) was calculated from the titration result after subtracting four times the measured MoS42- concentration and five times the excess of ΣMo over MoS42-. Under the experimental conditions, this excess is due mostly to dissolved [(FeS)2(MoS4)2]4-. Precipitate characterization. Precipitates used in the solid-phase analyses were obtained by filtering test solutions with 47 mm diameter, 0.45 μm pore size Tuffryn® or Metricel® membrane filters (Pall Life Sciences). The filtered solids were washed with several ~50 mL portions of ultrapure H2O and dried in the anaerobic chamber. The filter with filtrand was stored in the anaerobic chamber until analysis. Precipitates were freeze-dried prior to x-ray diffractometry (Rigaku Ultima-4 X-ray diffractometer, Cu Kα radiation, excitation energy 40 kV and 44 mA, scan rate 5°/min, step size 0.01°). Energy dispersive spectroscopy was performed on a ~0.5 cm2 piece of filter fixed to a nylon

5 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 134 3 4 135 5 6 136 7 8 9 137 10 11 138 12 13 139 14 15 16 140 17 18 141 19 20 142 21 22 23 143 24 25 144 26 27 145 28 29 146 30 31 32 147 33 34 148 35 36 149 37 38 39 150 40 41 151 42 43 152 44 45 46 153 47 48 154 49 50 155 51 52 156 53 54 55 157 56 57 158 58 59 60

Page 6 of 52

stub with carbon tape (Jeol 6510 LV scanning electron microscope equipped with a Thermo-Noran EDS System 7 and a Thermo-Noran Silicon Drift Detector; 20 kV accelerating voltage; 16 mm working distance; 0.066 μm spot size). More than 20,000 counts were obtained at each spot, implying counting errors of 1% or less; about 95 spots were assayed per sample. Peroxide digestion was used for elemental analyses. Precipitates were resuspended from filters by immersing the filters in water and shaking. The filters were removed and the water evaporated to dryness. Between 10 and 15 mg of dried precipitate was placed in a 50 mL Ni crucible with about 80 to 120 mg Na2O2. After capping the crucible, its contents were fused gingerly for 20 minutes over an air-CH4 flame. Upon cooling, the products in the crucible were leached with water, which was then filtered (0.45 µm Metricel® membrane filters). The filtrands were dissolved into 0.1 M HCl and analyzed by atomic absorption spectroscopy (Fe, Mo, Mg). Filtrate solutions were analyzed using reverse-phase ion pair chromatography (IPC: SO42-, MoO42-; Dionex IonPac NS1, conductivity detection). As a check, sulfate was also quantified turbidimetrically. Additional details concerning the turbidimetric and IPC methods are found in Supporting Information. X-ray Absorption Fine Structure spectroscopy was performed at beamline 13-BM-D at the Advanced Photon Source, Argonne National Laboratory. Powdered precipitates were packed into Teflon sample holders that were then sealed with Kapton tape within a glove box containing N2 atmosphere. During data collection, a constant flux of He was injected in a plastic bag containing our sample holders to prevent any oxidation. A Si (111) double crystal monochromator was used in conjunction with harmonic rejection mirrors. Fluorescence X-rays were measured using a Canberra 16 element Ge detector. The incident beam intensity was detuned by ∼20–30% to reject higherorder harmonic frequencies. Spectra were collected with energy and wave number resolution prior to the edge (19,850–19,980 eV); across the Mo K-edge (19,980–20,050 eV); and throughout the EXAFS region (20,050–20,700 eV) at 4 eV, 0.3 eV, and 0.05Å−1Å, respectively. At least three

6 ACS Paragon Plus Environment

Page 7 of 52 1 2 159 3 4 160 5 6 161 7 8 9 162 10 11 163 12 13 164 14 15 16 165 17 18 166 19 20 167 21 22 23 168 24 25 169 26 27 170 28 29 171 30 31 32 172 33 34 173 35 36 174 37 38 39 175 40 41 176 42 43 177 44 45 46 178 47 48 179 49 50 180 51 52 181 53 54 55 182 56 57 183 58 59 60

ACS Earth and Space Chemistry

spectra per sample were merged to improve the signal-to-noise ratio and resolution. Energy calibration was maintained by simultaneous measurements of a Mo(0) foil in transmission mode as an internal standard. The first and fourth peaks in the first derivative of the Mo(0) foil were assigned 20,000.0 and 20,039.1 eV. For data interpretation, five Mo standards were also characterized: Mo(IV)O2(s), Mo(IV)S2(s), Mo(VI)O3(s), Mo(VI)O42−(aq), and Mo(VI)S42−(aq). All data were processed and analyzed using the Demeter software package.37 Two features from the XANES region are particularly useful for delineating differences between Mo compounds: the pre-edge features (i.e. the presence, amplitude, and position of a peak before the edge) and the maximum peak (Em). A pronounced pre-edge feature is characteristic of terminal Mo=O double bonds (MoO42−) and an octahedral configuration. This pre-edge feature is less pronounced for Mo=S double bonds (MoS42−) and absent when only single bonds are present (i.e. MoO2(s) and MoS2(s)). Em increases with the Mo oxidation state.18,19 Thermodynamic Model. We calculated speciation in the solutions by solving the following mass balance equation for dissolved Fe: ΣFe = Fe2+ + FeOH+ + FeCl+ + Fe(HS)20 + [FeO(OH)MoS4]3- + 2[(Fe2S2)(MoS4)2]4-

(1)

Iron sulfide cluster complexes were excluded from this equation in the light of recent findings.38,39 Assuming chemical equilibrium, this equation can be solved for Fe2+ by replacing the term for each species by an expression that includes an analytically determined ligand concentration and a stability constant obtained from published sources.9,40,41 Activity coefficients were evaluated by the Davies equation. When Fe2+ had been calculated, the concentrations of all complexes then could be obtained by back substitution into the equilibrium constant expressions. In solving equation 1, the measured total Mo concentrations were not used. Therefore, the validity of this thermodynamic model could be tested by comparing the sum of calculated Mo species concentrations against the measured total Mo. Results of this test are shown in Figure 1, which reveals excellent agreement between measured and calculated values except for a small

7 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 184 3 4 185 5 6 186 7 8 9 187 10 11 188 12 13 189 14 15 16 190 17 18 191 19 20 192 21 22 23 193 24 25 194 26 27 195 28 29 196 30 31 32 197 33 34 198 35 36 199 37 38 39 200 40 41 201 42 43 202 44 45 46 203 47 48 204 49 50 205 51 52 206 53 54 55 207 56 57 58 59 60

Page 8 of 52

departure in the sample having the highest Mo concentration. This agreement verifies the model and the choices of stability constants, as well as the implicit assumption that all dissolved but not precipitated Mo is in the +VI oxidation state. 3. RESULTS Precipitate Formation. Table 1 summarizes both initial and final Mo, Fe and ΣS(-II) concentrations for all runs. Figure 2 shows a representative time course for an experiment. At mildly alkaline pH, no precipitation occurs from a solution containing only thiomolybdates and sulfide, but a stable colloidal suspension forms as soon as Fe(II) is added. The inset graph in Figure 2 compares optical absorption by the solution before and after colloids were flocculated with 0.1 M MgCl2 and removed by filtration. Prior to MgCl2 introduction, MoS42- can be quantified by intense absorption peaks at 318 and 468 nm, but colloids must be flocculated and removed in order to quantify dissolved Mo, Fe, and sulfide. The initially-formed colloids are likely to consist mainly of iron sulfide, which precipitates rapidly.42 The high ΣS-II/ΣFeII ratios and neutral pH of our solutions inhibit mackinawite crystal growth,43 preserving colloidal FeS in an unstable state as it reacts slowly with MoS42- (Figure 2). About half of the total MoS42- losses occur in the first day, but losses continue for many days. To allow precipitation to reach completion, we aged precipitates for 300 to 4000 h (4th column, Table 1), but we found no analytically significant trends with aging beyond 300 h. Curiously, the visual appearance of flocculated particles differed in solutions with high vs. low initial Fe/Mo mole ratios (NFe/Mo). At initial NFe/Mo>1, flocculated particles are large, wispy black sheets that settle appreciably within ~20 minutes after MgCl2 addition. In contrast, flocculated particles in solutions with initial NFe/Mo104, causing a resurgence of dissolved MoS42- absorbance from fresh precipitates. On the other hand, after only a few hours of aging, the figure shows that precipitates no longer release most of their MoS42-. Apparently, aging transforms the precipitates to much less soluble secondary products in this time period. Nevertheless, a few percent of precipitated MoS42- remain in an EDTA-reactive state, possibly as surface layers on the colloidal particles. Solution composition. At the pH and sulfide concentrations of these runs, optical absorption showed that MoS42- was overwhelmingly the predominant ion in the MoOnS4-n2- (n = 0 to 4) series, as anticipated from thermodynamics.34 Table 2 presents concentrations of selected other species obtained from the thermodynamic model. The fourth and fifth columns give the concentrations of the two most abundant Fe and Fe-Mo complexes. Not shown are the minor species that were considered (FeOH+, FeCl+ and [FeO(OH)MoS4]3-), but the sixth column shows the sum of Mo concentrations in all dissolved species in each run. This sum was plotted against the measured total Mo concentration in Figure 1 to verify the thermodynamic model. The last two columns in Table 2 list logarithms of activity products, QFeS and QMoS3, which we define as: QFeS = {Fe2+}{H2S} / 10-2pH

(2)

QMoS3 = {MoS42-} 10-2pH / {H2S}

(3)

(In this paper, braces denote activities of dissolved species; Q designates any activity product whereas K designates an activity product at equilibrium with a solid phase.) The Q values are proportional to the activities of FeS and MoS3 components (a component being one of the minimum number of independently variable constituents required to specify the composition of a phase). For

9 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 233 3 4 234 5 6 235 7 8 9 236 10 11 237 12 13 238 14 15 16 239 17 18 240 19 20 241 21 22 23 242 24 25 243 26 27 244 28 29 245 30 31 32 246 33 34 247 35 36 248 37 38 39 249 40 41 250 42 43 251 44 45 46 252 47 48 253 49 50 254 51 52 255 53 54 55 256 56 57 257 58 59 60

Page 10 of 52

example, we can write aFeS = QFeS/Cs where the constant of proportionality, Cs, defines the standard state and thus the aFeS scale. If we preferred to designate pure mackinawite (tetragonal FeS) as having aFeS = 1, then Cs would be equated to mackinawite’s conventional solubility product constant (Kmack = 103.21 at 25oC).44 If, on the other hand, we wished to designate nanocrystalline FeS as having aFeS = 1, then Cs would be KFeSam = 104.87.45 No choice of standard state is more correct than another, and for simplicity we will let Cs be unity, thus employing the Q values themselves as measures of the activities of their respective components. (In this case, the standard states for aFeS and aMoS3 are ideal aqueous solutions in which the component ions are present at unit concentrations.) In Figure 4, log QMoS3 has been plotted against log QFeS. The excellent inverse correlation reveals that the activities of FeS and MoS3 are linked by a solubility equilibrium. The GibbsDuhem equation specifies that the slope in this kind of plot is equal to minus the molar Fe/Mo ratio in the precipitate. Because log QFeS and log QMoS3 are linearly correlated, with an approximate slope of -1 (actual error-weighted regression value -0.94 ± 0.11, R2 0.95), we infer that over the > 2 log unit range of log QFeS in these experiments, the phase with which the solutions have equilibrated possesses a fixed composition with Fe/Mo = 1. Thus, the phase must be FeMoS4. If this phase were a solid solution having a significant range in Fe/Mo ratio, then the data in Figure 4 would display a non-linear trend. Employing the regression equation in Figure 4 and the definitions in equations (2) and (3), the solubility product constant of this FeMoS4 precipitate is given by: KFeMoS4 = (QMoS3)(QFeS) = {Fe2+}{MoS42-} = 10-14.95 ± 0.06

(4)

The standard free energy of formation from its elements of this colloidal precipitate is ΔGf0 = −1227.2 kJ/mol. For comparison, published data indicate that the corresponding solubility product for FeMoO4 is 10-10.48 and ΔGf0 = −975 kJ/mol.46 The large difference in solubility shows why converting MoO42- to MoS42- in sulfidic waters initiates Mo precipitation. The stability constant

10 ACS Paragon Plus Environment

Page 11 of 52 1 2 258 3 4 259 5 6 260 7 8 9 261 10 11 262 12 13 263 14 15 16 264 17 18 265 19 20 266 21 22 23 267 24 25 268 26 27 269 28 29 270 30 31 32 271 33 34 272 35 36 273 37 38 39 274 40 41 275 42 43 276 44 45 46 277 47 48 278 49 50 279 51 52 280 53 54 55 281 56 57 282 58 59 60

ACS Earth and Space Chemistry

and ΔGf0 values derived here refer to FeMoS4, not the secondary product formed after a few hours of aging. Because the secondary product has formed irreversibly, we can learn nothing of its stability from measurements of solution composition. Precipitate Properties. Figure 5 presents x-ray diffraction patterns for precipitates obtained at three initial Fe/Mo mole ratios. With no Mo in the initial solution (top trace), the pattern is similar to that of nanoparticulate FeS.47,48 The most prominent peaks of crystalline mackinawite can be recognized, but they are broadened and shifted to smaller 2θ (larger interplanar distances). In the absence of added Mo, these effects must be attributed to small particle size and the resulting lattice relaxation. Adding a small amount of MoS42- to the initial solution (middle trace) dilutes the mackinawite peaks and adding MoS42- in excess to Fe2+ (bottom trace) largely erases them. Significantly, no new peaks are created. Evidence for weak mackinawite peaks in the bottom trace may be due to survival of a small amount of mackinawite in cores of what are otherwise amorphous particles. As shown in Figure 6, we also characterized compositions of several precipitates by energy dispersive spectroscopy. The panels depict experiments that are arranged according to the ratio of initial Fe to initial Mo with the most Fe-rich experiment at the top. The horizontal axes in this figure are the ratios of counts obtained for the Fe peak at 6.40 keV to the sum of counts for Fe plus the counts for overlapping Mo and S peaks at ~2.30 keV. Overlap prevented obtaining independent Mo and S counts and, thus, quantifying elemental ratios. The precipitate in the top panel contains solely iron sulfide. The observed dispersion in count ratios probably is indicative of the instrument’s energy resolution with colloidal particles. In panels b and c, where Mo is present but the amount is not sufficient to convert all the FeS to FeMoS4, the observed dispersion in count ratios is much larger than in panel a. Apparently, the precipitate is chemically heterogeneous, and the electron beam is encountering the FeS and FeMoS4 phases in varying proportions as different spots are sampled. This evidence implies that the

11 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 283 3 4 284 5 6 285 7 8 9 286 10 11 287 12 13 288 14 15 16 289 17 18 290 19 20 291 21 22 23 292 24 25 293 26 27 294 28 29 295 30 31 32 296 33 34 297 35 36 298 37 38 39 299 40 41 300 42 43 301 44 45 46 302 47 48 303 49 50 304 51 52 305 53 54 55 306 56 57 307 58 59 60

Page 12 of 52

precipitates consist of discrete phases, not a homogeneous solid solution phase, but does not exclude limited Fe/Mo variation within either phase. In panels d and e, where enough Mo relative to Fe is initially available to convert all FeS to FeMoS4, the dispersion in precipitate compositions again resembles that in panel a, suggesting that a single FeMoS4 phase is present. Nevertheless, further increase in initial Mo produces a new phase still richer in Mo (panel f). Figure 7 presents XANES spectra for several standards as well as for precipitates generated by runs 6, 7, 9 and 12. A pre-edge peak near 20005 eV is seen in the standards containing Mo(VI) but is absent from the standards containing Mo(IV). This peak is also absent from the spectra of all precipitates. Additionally, the precipitates’ maximum absorption energies (Em) are similar to those for Mo(IV) standards (Figure 7 and Table 3). These features provide evidence that Mo has been reduced. Additionally, fits to EXAFS spectra (Table 4) indicate that the first two coordination shells around Mo are occupied by approximately 2 and 3 sulfur atoms at distances of 2.30 – 2.32 Å and 2.44 – 2.47 Å, respectively. These interatomic distances are clearly longer than interatomic distances for MoS42- (2.18 Å). If we consider the average of the shorter Mo–S interatomic distance for our precipitate (2.31 Å), it is similar to data for the type S sediment samples from Lake Cadagno (2.26–2.38 Å).18 The longer Mo-S distances were not observed in Cadagno, possibly because of greater disorder. When the Lake Cadagno EXAFS modeling was run, only one type of Mo–S bond was considered and tested for EXAFS modeling. The EXAFS evidence rules out MoS2(S), MoS3(S) and MoS42- as major constituents in the aged precipitates. Modeling of the third shell reveals that Mo is bound to an Fe atom at a distance of 2.79–2.81 Å. These Mo–Fe interatomic distances, although longer, remain comparable to those determined in Lake Cadagno type S samples (2.71– 2.73 Å). These comparisons suggest that the precipitates are structurally similar to products found in natural euxinic environments. Huerta-Diaz and Morse found that Mo was not appreciably dissolved from modern sulfidic sediments by 1 M HCl but was dissolved to a large extent by concentrated HNO317 an oxidizing

12 ACS Paragon Plus Environment

Page 13 of 52 1 2 308 3 4 309 5 6 310 7 8 9 311 10 11 312 12 13 313 14 15 16 314 17 18 315 19 20 316 21 22 23 317 24 25 318 26 27 319 28 29 320 30 31 32 321 33 34 322 35 36 323 37 38 39 324 40 41 325 42 43 326 44 45 46 327 47 48 328 49 50 329 51 52 330 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

acid. They interpreted this finding as evidence that Mo was mainly hosted by pyrite, which requires concentrated HNO3 for dissolution. We applied this test to a precipitate that formed under conditions of Run 1. Wet precipitate was placed in each of two centrifuge tubes. Air-saturated, 1 M HCl was added to one and concentrated HNO3 to the other. Visible solid disappeared immediately in HNO3, but persisted in HCl for >250 days; AAS analyses confirmed that only the HNO3 solution had appreciably dissolved Mo. This result implies that oxidation is required to dissolve the Mo(IV) contained in aged Fe-Mo-S precipitates. It also weakens support for the inference of Huerta-Diaz and Morse that Mo is hosted mainly by pyrite in modern sulfidic sediments. Because sulfidic sediments exposed to NaOH have been shown to leach a significant portion of particle-bound Mo, 50 workers have concluded that Mo must be hosted by humic materials. As a test, we added 1 M NaOH to a precipitate formed under the conditions of Run 6. After filtering, AAS analyses of the filtrate showed that >20% of Mo initially present in the precipitate had dissolved. These results imply that leaching of Mo by NaOH can not be used on its own as evidence of organic-bound Mo in sulfidic sediments. Sulfide as reductant. Reduction of Mo(VI) in solution to Mo(IV) in precipitates requires a reducing agent, which in these experiments can only be sulfide. The electron transfer could be internal, between Mo(VI) and coordinating S(-II) atoms, or external, between HS-(aq) and FeMoS4(s). These options can be represented as: FeMoS4(s)  FeMoS2(S2)(s)

(5)

FeMoS4(s) + 1/(n-1) HS-  FeMoS3(s) + 1/(n-1) Sn2- + 1/(n-1) H+

(6)

where (S2) represents a disulfide ligand in the precipitate structure and Sn2- represents a dissolved polysulfide ion.

13 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 331 3 4 332 5 6 333 7 8 9 334 10 11 335 12 13 336 14 15 16 337 17 18 338 19 20 339 21 22 23 340 24 25 341 26 27 342 28 29 343 30 31 32 344 33 34 345 35 36 346 37 38 39 347 40 41 348 42 43 349 44 45 46 350 47 48 351 49 50 352 51 52 353 53 54 55 354 56 57 355 58 59 60

Page 14 of 52

Molybdenum sulfides that include disulfide ligands are well known. Simply acidifying a MoS42- solution rapidly produces dissolved species such as MoV2S4(S2)22-.51 Even amorphous MoS3(s), which is nominally a Mo(VI) compound, is known actually to be MoIVS(S2).52-54 We approached the question of whether reaction 5 or 6 best describes the evolution of our precipitates in two ways. First, we looked for optical absorption evidence for dissolved polysulfide ions.55 This approach was obstructed in most runs by swamping optical absorption from MoS42and [(FeS)2(MoS4)2]4-. Only in one case, run 16, was interference from these sources low enough to obtain the concentration of dissolved S(0) in polysulfides. The result was about ⅓ mole dissolved S(0) per mole Mo precipitated, implying that at least ⅔ mole S(0) per mole Mo was retained by the precipitate. Retention could be higher if some O2 contamination were responsible for the dissolved S(0) that was observed. Second, we analyzed several aged precipitates. From three solutions initially containing NFe/Mo = 0.20, 0.13 and 0.10, we obtained precipitates with empirical formulas of Fe0.78MoS3.83, Fe0.83MoS3.75 and Fe0.63MoS3.60 respectively. These contained very small amounts of Mg2+ adsorbed from the electrolyte. Additionally, 45% to 50% of their masses was unexplained and presumed to be water, a common component of flocculated colloids. If all Fe in these formulas has a +2 charge and all molybdenum has a +4 charge (as required by x-ray spectroscopic data), then charge balance requires that anions in these formulas possess a total charge of -5.56, -5.66 and -5.26 respectively. If all of this negative charge is supplied by the sulfur (i.e. none from bound OH− or Cl−), then the implied final formulas would be Fe0.78MoS1.73(S2)1.05, Fe0.83MoS1.91(S2)0.92 and Fe0.63MoS1.66(S2)0.97. Apparently, these precipitates have retained about one mole of zero-valent sulfur (contained in S22-) per mole of Mo, suggesting that reaction 5 is better than reaction 6 as a representation of the secondary material’s formation. These analyses suggest that precipitates kept in Mo-rich solutions during weeks of aging contain Fe/Mo somewhat smaller than the value of unity inferred for the initial precipitates from Figure 4. Possibly, FeMoS4 takes on additional Mo or loses some Fe as it converts to its aged form.

14 ACS Paragon Plus Environment

Page 15 of 52 1 2 356 3 4 357 5 6 358 7 8 9 359 10 11 360 12 13 361 14 15 16 362 17 18 363 19 20 364 21 22 23 365 24 25 366 26 27 367 28 29 368 30 31 32 369 33 34 370 35 36 37 371 38 39 40 372 41 42 373 43 44 374 45 46 375 47 48 49 376 50 51 377 52 53 378 54 55 56 379 57 58 380 59 60

ACS Earth and Space Chemistry

Alternatively, the recovered products might contain a small amount of the additional, Mo-rich phase revealed by Figure 5f. 4. DISCUSSION Mo removal from sulfidic natural waters. To optimize experimental control, the precipitates reported here were produced in solutions containing higher Mo and S(-II) concentrations and higher pH than normally found in natural sulfidic waters. Can these precipitates form in nature? In order to investigate this question, we apply a thermodynamic model to the extensive data available for the Black Sea, the largest modern euxinic basin. This evaluation relies on an assumption that dissolved species reach aqueous-phase equilibrium during the centuries that deep water is resident in the Black Sea. From published analyses and stability constants, we first calculate the following activities in the Black Sea’s sulfidic waters: {Fe2+}BS, {H2S}BS and {MoS42-}BS (the BS subscript denotes activities in the Black Sea’s water). Sources of data and results are presented in Table S1.

Then,

we compare the {MoS42-}BS values to {MoS42-}eq values that would be required for saturation with respect to FeMoS4; {MoS42-}eq is calculated from a rearrangement of equation 4: {Mo𝑆42− }𝑒𝑞 =

10−14.95 {𝐹𝑒 2+ }

=

10−14.95 {𝐻2 𝑆}

(7

(𝑄𝐹𝑒𝑆 )10−2𝑝𝐻

If precipitation of FeMoS4 controls total Mo in the Black Sea, {MoS42-}BS should be approximately equal to {MoS42-}eq. It will not matter that FeMoS4 in the Black Sea would soon transform irreversibly to FeMoS2(S2). In an analogous situation, labile FeS phases control dissolved Fe in sulfidic waters even though these phases are transforming irreversibly to FeS2. Figure 8 presents the results of this assessment. At the time when the Black Sea was sampled, sulfide concentrations exceeding 1 μM appeared below a suboxic-to-anoxic transition at 90 m depth. Panel A shows that with increasing depth, the hydrogen sulfide activity, {H2S}BS, rose continuously all the way to the bottom. Although {Fe2+}BS rose with increasing depth in the uppermost sulfidic waters, owing to reductive dissolution of FeIII in sinking particles, this trend

15 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 381 3 4 382 5 6 383 7 8 9 384 10 11 385 12 13 386 14 15 16 387 17 18 388 19 20 389 21 22 23 390 24 25 391 26 27 392 28 29 393 30 31 32 394 33 34 395 35 36 396 37 38 39 397 40 41 398 42 43 399 44 45 46 400 47 48 401 49 50 402 51 52 53 403 54 55 404 56 57 405 58 59 60

Page 16 of 52

reversed below 150 m. The likely reason is that the water became saturated with FeS, forcing {Fe2+} to decline as {H2S} continued to increase. Panel B shows that total dissolved Mo in the uppermost sulfidic water column changed negligibly with respect to overlying oxic and suboxic waters, but it began to decline below about 150 m. In Panel C, the solid red triangles show that despite negligible variation in total dissolved Mo above 150 m, {MoS42─}BS rose sharply in this region; this was due to conversion of MoO42─ to MoS42─. In contrast, below 150 m, the figure shows that {MoS42─}BS ceased to change appreciably even though total dissolved Mo then began to decline very significantly. In this region, any additional MoS42- produced by rising sulfide was removed from the water, causing loss of total Mo at nearly constant {MoS42─}BS. Most of the loss of total dissolved Mo from the entire water column occurred in the zone where some removal process was fixing {MoS42-}BS at nearly constant values. The disconnect between behaviors of total dissolved Mo and {MoS42-}BS is striking. Is it plausible that saturation with respect to FeMoS4 is the process that fixed {MoS42─}BS? The answer to this question depends on the values of {Fe2+}BS, or alternatively, QFeS. The open diamond symbols in Panel C depict {MoS42-}eq values obtained from equation 7 by assuming that ambient {Fe2+}BS at each depth controlled {MoS42-}eq. The diamond symbols show that {MoS42─}eq is greater than {MoS42─}BS by more than two orders of magnitude. Seemingly, the Black Sea is immensely under saturated with respect to FeMoS4. If correct, then FeMoS4 precipitates cannot control dissolved Mo in the Black Sea. However, this difficulty can be circumvented if the operative {Fe2+}BS values are much higher than the ambient values in bulk water. (Bulk water is defined here simply as water in quantities usually taken for analysis.) Perhaps higher {Fe2+}BS exists in microenvironments where reductive dissolution of iron is taking place, for example at surfaces of Fe(III)-rich particles, or in pore spaces of composite particles and flocs. Support for this idea can be found by considering FeS. Below 150 m, log QFeS in bulk water averages 2.60 ± 0.17 (Table S1), which is much smaller

16 ACS Paragon Plus Environment

Page 17 of 52 1 2 406 3 4 407 5 6 408 7 8 9 409 10 11 410 12 13 411 14 15 16 412 17 18 413 19 20 414 21 22 23 415 24 25 416 26 27 417 28 29 418 30 31 32 419 33 34 420 35 36 421 37 38 39 422 40 41 423 42 43 424 44 45 46 425 47 48 426 49 50 427 51 52 428 53 54 55 429 56 57 58 59 60

ACS Earth and Space Chemistry

than 3.38 at mackinawite saturation or 4.03 at amorphous FeS saturation.56 Thus, the deep waters of the Black Sea are also immensely under saturated with respect to FeS, even though FeS and phases derived from it (greigite, pyrite) are known to exist there.25,28 Experiments show that nucleation of FeS nanoparticles occurs at log QFeS ≥ 4.87, 45 which is 2.3 log units larger than average log QFeS in the deep Black Sea’s bulk water. If we suppose that iron sulfide precursors form in microenvironments where log QFeS ≥ 4.87, then FeMoS4 could also form there. This is demonstrated in Figure 8C by open square symbols, which show {MoS42-}eq values calculated from equation 7 by assuming that log QFeS = 4.87. The resulting {MoS42─}eq values are in reasonable agreement with {MoS42-}BS values, fulfilling the feasibility criterion for FeMoS4 control of total dissolved Mo. If FeMoS4 colloids indeed form on or within Fe(III)-bearing particles, then they would be carried to the sediments with the particles. As FeMoS4 settles aboard its natal particles, it can experience elevated Fe(II) activity only until the Fe(III)-pool has been exhausted. The Fe(II) activity will then decline to ambient levels, dropping log QFeS to values permissive of FeMoS4 dissolution. How could FeMoS4 persist in these undersaturated waters and ultimately survive as a Mo burial product? In addition to the demonstrated irreversibility of its formation, FeMoS4 dissolution requires strongly oxidizing conditions. Given the persistently euxinic water column of the Black Sea, it is reasonable to assume FeMoS4 could survive to burial within Black Sea sediments. In euxinic basins like the Black Sea, Mo probably is trapped in sediments permanently as FeMoS2(S2) colloids. EXAFS patterns resembling those of FeMoS2(S2) colloids have been reported from Phanerozoic black shales.7 However, in environments subject to seasonal fluctuations between oxic and anoxic conditions, trapped Mo is apt to be partially or fully dissolved by exposure to O2; rapid oxidation has been demonstrated in resuspension experiments.57

17 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 430 3 4 431 5 6 432 7 8 9 433 10 11 434 12 13 435 14 15 16 436 17 18 437 19 20 438 21 22 23 439 24 25 440 26 27 441 28 29 442 30 31 32 443 33 34 444 35 36 445 37 38 39 446 40 41 447 42 43 448 44 45 46 449 47 48 450 49 50 451 51 52 452 53 54 55 453 56 57 454 58 59 60

Page 18 of 52

Some continental shelf sediments are considerably enriched in Mo even though negligible sulfide is found in the seawater above them. In these environments, Mo is believed to diffuse into the seafloor, where it encounters sulfide.58,59 Because reduction of Fe(III) is ongoing in these environments,60 the Mo-fixing mechanism described above is expected to operate in the same manner. Comparison to synthetic compounds. Clues to the atomic structure of the precipitates can be gleaned from known structures of synthetic Fe-Mo-S compounds by comparing Mo-S and MoFe interatomic distances. The Mo-S distances in Table 4 are conspicuously long and suggestive of cuboid compounds.61-64 Most Fe-Mo-S cuboids have been prepared in non-aqueous solvents and stabilized by organic ligands, but the cuboidal [FeMo3S4(H2O)10]4+ cation illustrated in Figure 9 has been prepared in 2 M aqueous HCl.65 The interatomic distances in this cation (see caption) are similar to those in our precipitates (Table 4). In contrast, most non-cuboid Fe-Mo-S compounds have substantially shorter Mo-S interatomic distances and slightly longer Mo-Fe distances. For example, the structure of the non-cuboid complex, [(FeS)2(MoS4)2]4-, which is the principal dissolved Fe-Mo-S complex in our solutions, has Mo-Sμ2 bridge interatomic distances of 2.24Å and Mo-Sterminal distances of 2.15Å; Mo-Fe distances are 2.78Å.66 Comparable distances have been reported for a number of similar, non-cuboidal complexes.67 Additionally, the formal oxidation state in [(FeS)2(MoS4)2]4- is MoVI, but is MoIV in our aged precipitates. These comparisons indicate that [(FeS)2(MoS4)2]4- and similar non-cuboid complexes are poor models for the building blocks of our aged precipitates. They also emphasize the point that Mo in the precipitates is chemically quite distinct from Mo in solution. Core compositions in known cuboid compounds form a complete series from Fe4S4 to Mo4S4.64 If we double the inferred empirical formula for our secondary product (reaction 5), the result is Fe2Mo2S4(S2)2. This suggests that our precipitates might have [Fe2Mo2S4]4+ cores coordinated at the Mo positions by disulfide ligands (illustrated in Figure 9). On the other hand,

18 ACS Paragon Plus Environment

Page 19 of 52 1 2 455 3 4 456 5 6 457 7 8 9 458 10 11 459 12 13 460 14 15 16 461 17 18 462 19 20 463 21 22 23 464 24 25 465 26 27 466 28 29 467 30 31 32 468 33 34 469 35 36 470 37 38 39 471 40 41 472 42 43 473 44 45 46 474 47 48 475 49 50 476 51 52 477 53 54 55 478 56 57 479 58 59 60

ACS Earth and Space Chemistry

panel f of Figure 5 hints that initial solutions having very low NFe/Mo ratios produce precipitates that are much less Fe-rich, possibly containing FeMo3S4 or even Mo4S4 cores. From solutions having extremely low initial NFe/Mo, precipitates with less than 1 mol % Fe have been obtained,11 but only in the presence of microbial cells or their non-living degradation products. The chemical mechanism by which these microbial products promote precipitation remains to be determined. The inferred disulfide ligands in our precipitates probably would be reactive. They might participate in cross-linking, polymerization reactions that could explain how a monomer, like the hypothetical one in Figure 9, aggregates to make colloids. Similarly, the disulfide ligands, might form covalent bridges to the S22- ions on pyrite surfaces or to organic polysulfides. Attachment of Mo to such materials has been proposed as scavenging mechanisms for Mo in nature.10,68 Jordisite—an example? Jordisite is a rare, dark gray-to-black sulfide that historically was believed to consist of x-ray amorphous MoS2. Its composition is uncertain. Its structure is known to differ from that of both molybdenite and synthetic MoS3, but is otherwise unknown.69 Its particle sizes are invariably in the colloidal range. Miners recognize it by its ready conversion to molybdenum blue compounds when exposed to urine, a property inconsistent with molybdenite.70 It forms in deep-seated groundwaters at temperatures up to 200oC and is often associated with uranium ores.70-75 It is also found in hypogene zones of weathered molybdenum ore bodies.76,77 Owing to its colloidal nature, no pure jordisite is available for analysis. The few reported analyses include significant amounts of constituents other than Mo and S (for example: iron, silica and water). Previous workers have attributed the iron to contamination by colloidal iron sulfide particles. In Table 5, we recalculate the available published jordisite analyses so that the Fe is included in the empirical formulas. The resulting formulas seem too regular to be attributed to arbitrary mixtures of Fe and Mo sulfides. On the other hand, Fe/Mo mole ratios in jordisites cluster below 0.6, making them distinctly different from our aged precipitates. The solutions from which jordisite precipitates must differ significantly in composition and temperature from those in this

19 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 480 3 4 481 5 6 482 7 8 9 483 10 11 484 12 13 485 14 15 16 486 17 18 487 19 20 488 21 22 23 489 24 25 490 26 27 491 28 29 492 30 31 32 493 33 34 494 35 36 495 37 38 39 496 40 41 497 42 43 498 44 45 46 499 47 48 500 49 50 501 51 52 502 53 54 55 503 56 57 504 58 59 60

Page 20 of 52

work. If jordisite colloids contain cuboid structures, then FeMo3S4 predominates over Fe2Mo2S4. The findings here and in Ref.11 hint that nature may harbor a whole class of Fe-Mo-S colloids, of which jordisite is merely a member. 5.

CONCLUSIONS When mixtures of dissolved Fe and Mo are precipitated by sulfide, two colloidal products

form. When Fe/Mo in solution exceeds unity, which is the usual case in nature, the early-formed colloidal products are composed of FeMoS4 and poorly crystalline mackinawite. When dissolved Fe/Mo is less than unity, very little iron monosulfide persists, because it is nearly all converted to FeMoS4. In either case, after only a few hours of aging, most FeMoS4 spontaneously reduces irreversibly to a non-crystalline, Mo(IV) material. The reduced product is soluble in an oxidizing acid, concentrated HNO3, but not in 1 M HCl. At the very lowest Fe/Mo investigated in this study, a further, very Fe-depleted colloidal phase forms. Precipitates formed in these experiments appear to share a number of characteristics with Mo contained within sulfidic sediments (although information about the latter is still very limited). Both contain Mo in similar atomic environments, as indicated by similar Mo-S and Mo-Fe interatomic distances. Both appear to contain Mo(IV) rather than Mo(VI), and consequently both are dissolved only under oxidizing conditions. On the other hand, deposition of Mo via a FeMoS4 precursor in the Black Sea requires {Fe2+}BS or QFeS values that are several orders of magnitude greater than exist in bulk water and therefore requires positing existence of such values in microenvironments. SUPPORTING INFORMATION Section S1, description of turbidimetric and IPC analyses of peroxide digestion solutions; Figure S1, EXAFS spectra of run 7 with k-representation in the left panel and the Fourier transform shown in the right panel; Table S1, Black Sea data for Figure 8 and associated citations.

20 ACS Paragon Plus Environment

Page 21 of 52 1 2 505 3 4 506 5 6 507 7 8 9 508 10 11 509 12 13 510 14 15 16 511 17 18 512 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

ACKNOWLEDGEMENTS The authors thank Tais Dahl, Florian Sholz, two anonymous reviewers, and AE Sumit Chakraborty for constructive comments that significantly improved an earlier version of the manuscript. The authors express gratitude to Chad Wittkop for invaluable assistance with XRD analyses. TPV and AC thank the US National Science Foundation (Awards # EAR-1503567 and EAR-1503596) for funding this research. TPV and AC acknowledge the Donors of the American Chemical Society Petroleum Research Fund (PRF # 52201-UR2 and ACS-PRF 54583-DNI2) for supporting this research.

21 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 513 3 4 514 5 515 6 7 516 8 517 9 518 10 519 11 12 520 13 521 14 522 15 523 16 17 524 18 525 19 526 20 527 21 528 22 23 529 24 530 25 531 26 532 27 28 533 29 534 30 535 31 536 32 33 537 34 538 35 539 36 540 37 38 541 39 542 40 543 41 544 42 43 545 44 546 45 547 46 548 47 48 549 49 550 50 551 51 552 52 53 553 54 554 55 555 56 556 57 557 58 59 558 60

Page 22 of 52

References 1) Lyons, T.W.; Reinhard, C.T.; Planavsky, N.J. Nature. The rise of oxygen in Earth’s early ocean and atmosphere. 2014, 506, 307-315. 2) Smedley, P.L.; Kinniburgh, D.G. Molybdenum in natural waters: A review of the occurrence, distributions and controls. Appl. Geochem. 2017, 84, 387-432. 3) Scott, C.; Lyons, T.W.; Bekker, A.; Shen, Y.; Poulton, S.W.; Chu, X.; Anbar, A.D. Tracing the stepwise oxygenation of the Proterozoic ocean. Nature. 2008, 452, 456-460. 4) Scott, C.; Lyons, T.W. Contrasting molybdenum cycling and isotopic properties in euxinic versus non-euxinic sediments and sedimentary rocks: Refining the paleoproxies. Chem. Geol. 2012, 324-325, 19-27. 5) Chappaz, A.; Lyons, T.W.; Gregory, D.D.; Reinhard, B.C.; Gill, C.; Large, R.R. Does pyrite act as an important host for molybdenum in modern and ancient euxinic sediments? Geochim. Cosmochim. Acta. 2014, 74, 203-214. 6) Wagner, M.; Chappaz, A.; Lyons, T.W. Molybdenum speciation and burial pathway in weakly sulfidic environments: Insights from XAFS. Geochim. Cosmochim. Acta. 2017, 206, 18-29. 7) Helz, G.R.; Miller, C.V.; Charnock, J.M.; Mosselmans, J.F.W.; Pattrick, R.A.D.; Garner, C.D.; Vaughan, D.J. Mechanism of molybdenum removal from the sea and its concentration in black shales: EXAFS evidence. Geochim. Cosmochim. Acta. 1996, 60, 3631-3642. 8) Helz, G.R.; Bura-Nakić, E.; Mikac, N.; Ciglenečki, I. New model for molybdenum behavior in euxinic waters. Chem. Geol. 2011, 284, 323-332. 9) Helz, G.R.; Erickson, B.E.; Vorlicek, T.P. Stabilities of thiomolybdate complexes of iron; implications for retention of essential trace elements (Fe, Cu, Mo) in sulfidic waters. Metalomics. 2014, 6, 1131-1140. 10) Freund, C.; Wishard, A.; Brenner, R.; Sobel, M.; Mizelle, J.; Kim, A.; Meyer, D.A.; Morford, J.A. The effect of a thiol-containing organic molecule on molybdenum adsorption onto pyrite. Geochim. Cosmochim. Acta. 2016, 174, 222-235. 11) Dahl, T.W.; Chappaz, A.; Hoek, J.; McKenzie, C.J.; Svane, S.; Canfield, D.E. Evidence of molybdenum association with particulate organic matter under sulfidic conditions. Geobiol. 2017, 15, 311-323. 12) Ardakani, O.H., Chappaz, A., Sanei, H. and Mayer, B. (2016) Effect of thermal maturity on remobilization of molybdenum in black shales. Earth Planet. Sci. Lett. 449, 311-320. 13) Korolev, D.F. The role of iron sulfides in the accumulation of molybdenum in sedimentary rocks of the reduced zone. Geochemistry. 1958, 4, 450-463. 14) Sugawara, K.; Okabe, S.; Tanaka, M. Geochemistry of molybdenum in natural waters (II). J. Earth Sci. Nagoya Univ. 1961, 9, 114-128. 15) Bertine, K.K. The deposition of molybdenum in anoxic waters. Mar. Chem. 1972, 1, 43-53. 16) Volkov, I.I.; Fomina, L.S. Influence of organic material and processes of sulfide formation on distribution of some trace elements in deep-water sediments of Black Sea. In: E.T. Degens and D.A. Ross (eds) The Black Sea: Geology, chemistry and Biology. Amer. Assoc. Petroleum Geol. Memoir 20, p. 456-476, 1974. 17) Huerta-Diaz, M.A.; Morse, J.W. Pyritization of trace metals in anoxic marine sediments. Geochim. Cosmochim. Acta. 1992, 56, 2681-2702. 18) Dahl, T.W.; Chappaz, A.; Fitts, J.P.; Lyons, T.W. Molybdenum reduction in a sulfidic lake: Evidence from x-ray absorption fine-structure spectroscopy and implications for the Mo paleoproxy. Geochim. Cosmochim. Acta. 2013, 103, 213-231.

22 ACS Paragon Plus Environment

Page 23 of 52 1 2 559 3 560 4 561 5 6 562 7 563 8 564 9 565 10 11 566 12 567 13 568 14 569 15 570 16 17 571 18 572 19 573 20 574 21 22 575 23 576 24 577 25 578 26 27 579 28 580 29 581 30 582 31 32 583 33 584 34 585 35 586 36 37 587 38 588 39 589 40 590 41 42 591 43 592 44 593 45 594 46 47 595 48 596 49 597 50 598 51 599 52 53 600 54 601 55 602 56 603 57 58 59 60

ACS Earth and Space Chemistry

19) Ardakani, O.H.; Chappaz, A.; Sanei, H.; Mayer B. Effect of thermal maturity on remobilization of molybdenum in black shales. Earth Planet. Sci. Lett. 2016, 449, 311-320. 20) Havig, J.R.; McCormic,k M.L.; Hamilton, T.L.; Kump, L.R. The behavior of biologically important trace elements across the oxic/euxinic transition of meromictic Fayetteville Green Lake, New York, USA. Geochim. Cosmochim. Acta. 2015, 165, 389-406. 21) Mohajerin, T.J.; Helz, G.R.; Johannesson, K.H. Tungsten-molybdenum fractionation in estuarine environments. Geochim. Cosmochim. Acta. 2016, 177, 105-119. 22) Buffle, J.; Zali, O.; Zumstein, J.; De Vitre, R. Analytical methods for the direct determination of inorganic and organic species: Seasonal changes of iron, sulfur, and pedogenic and aquagenic organic constituents in the eutrophic lake Bret, Switzerland. Sci. Total Environ. 1987, 664, 41-59. 23) Skei, J.M. Formation of framboidal iron sulfide in the water of a permanently anoxic fjord—Framvaren Fjord, south Norway. Mar. Chem. 1988., 3, 345-352. 24) Lewis, B.L; Landing, W.M. The biogeochemistry of manganese and iron in the Black Sea. Deep-Sea Res. 1991, 38, S773-S803. 25) Muramoto, J.A.; Honjo, S.; Fry, B.; Hay, B.J.; Howarth, R.W.; Cisne, J.L. Sulfur, iron and organic carbon fluxes in the Black Sea: sulfur isotopic evidence for origin of sulfur fluxes. Deep-Sea Res. 1991, 38 Suppl. 2, S1151-S1187. 26) Perry, K.A.; Pedersen, T.F. Sulphur speciation and pyrite formation in meromictic ex-fjords. Geochim. Cosmochim. Acta. 1993, 57, 4405-4418. 27) Wilkin, R.T.; Barnes, H.L. Pyrite formation in an anoxic estuarine basin. Amer. J. Sci. 297, 620-650. 28) Cutter, G.A.; Kluckhohn, R.S. The cycling of particulate carbon, nitrogen, sulfur and sulfur species (iron monosulfide, greigite, pyrite and organic sulfur) in the water columns of Framvaren Fjord and the Black Sea. Mar. Chem. 1999, 67,149-160. 29) Li, X.; Cutter, G.A.; Thunell, R.C.; Tappa, E.; Gilhooly, W.P.III; Lyons, T.W.; Astor, Y.; Scranton, M.I. Particulate sulfur species in the water column of the Cariaco Basin. Geochim. Cosmochim. Acta. 2011, 75,148-163. 30) Bura-Nakić, E.; Viollier E.; Ciglenečki, I. Electrochemical and Colorimetric Measurements Show the Dominant Role of FeS in a Permanently Anoxic Lake. Environ. Sci. Technol. 2013, 47, 741-749. 31) Lohmayer, R.; Reithmaier, G.M.S.; Bura-Nakić, E.; Planer-Friedrich, B. Ion-pair chromatography coupled to inductively coupled plasma-mass spectrometry (IPC-ICP-MS) as a method for thiomolybdate speciation in natural waters. Anal. Chem. 2015, 87, 33883395. 32) Vorlicek, T.P.; Chappaz, A.; Groskreutz, L.M.; Young, N.; Lyons, T.W. A new analytical approach to determining Mo and Re speciation in sulfidic waters. Chem. Geol. 2015, 403, 52-57. 33) Chandrasekaran, J., Ansari, M.A.; Sarkar, S. Aging of ammonium tetrathiomolybdate(VI) in air: An example of induced electron transfer by external oxidant, oxygen. J. Less Common Met. 1987, 134, L23-L25. 34) Erickson, B.E.; Helz, G.R. Molybdenum(VI) speciation in sulfidic waters: Stability and lability of thiomolybdates. Geochim. Cosmochim. Acta. 2000, 64, 1149-1158. 35) Erickson, B.E. The speciation of molybdenum in sulfidic natural waters. Ph.D. dissertation. University of Maryland, College Park, 1998.

23 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 604 3 605 4 606 5 6 607 7 608 8 609 9 610 10 11 611 12 612 13 613 14 614 15 615 16 17 616 18 617 19 618 20 619 21 22 620 23 621 24 622 25 623 26 27 624 28 625 29 626 30 627 31 32 628 33 629 34 630 35 631 36 37 632 38 633 39 634 40 635 41 42 636 43 637 44 638 45 639 46 47 640 48 641 49 642 50 643 51 644 52 53 645 54 646 55 647 56 648 57 58 649 59 650 60

Page 24 of 52

36) Vorlicek, T.P. Toward a better understanding of molybdenum fixation in sediments: The roles of mineral catalysis, zero-valent sulfur, and metal sulfides. Ph.D. dissertation. University of Maryland, College Park, 2002. 37) Ravel, Á.; Newville, M.. ATHENA, ARTEMIS, HEPHAESTUS: data analysis for X-ray absorption spectroscopy using IFEFFIT. J. Synchrotron Radiat, 2005, 12, 537–541. 38) Bura-Nakić, E.; Krznarić, D., Helz, G.R.; Ciglenečki, I. Characterization of iron sulfide species in model solutions by cyclic voltammetry. Revisiting an old problem. Electroanal. 2011, 23, 1376-1382. 39) Krznarić, D.; Ciglenečki, I. Voltammetric study of an FeS layer on a Hg electrode in supersaturated FeS chloride solution. Environ. Chem. 2015, 12, 123-129. 40) Davison, W.; Phillips, N.; Tabner, B.J. Soluble iron sulfide species in natural waters: Reappraisal of their stoichiometry and stability constants. Aquatic Sci. 1999, 61, 23-43. 41) King, D. W.; Farlow, R. Role of carbonate speciation on the oxidation of Fe(II) by H2O. Mar. Chem. 2000, 70, 201-209. 42) Rickard, D. Experimental concentration-time curves for the iron(II) sulphide precipitation process in aqueous solutions and their interpretation. Chem. Geol. 1989, 78, 315-324. 43) Bourdoiseau, J.-A.; Jeannin, M.; Sabot, R.; Rémazeilles C.; Refait, Ph. Characterisation of mackinawite by Raman spectroscopy: Effects of crystallization, drying and oxidation. Corrosion Sci. 2008, 50, 3247-3255. 44) Benning, L.G.; Wilkin, R. T.; Barnes, H.L. Reaction pathways in the Fe-S system below 100oC. Chem. Geol. 2000, 167, 25-51. 45) Wolthers, M.; Charlet, L.; Der Linde, P.R.; Rickard, D.; van der Weijden, C.H. Surface chemistry of disordered mackinawite (FeS). Geochim. Cosmochim. Acta. 2005, 69, 34693481. 46) Wagman, D.D.; Evans, W.H.; Parker, V.B.; Schumm, R.H.; Halow, I.; Bailey, S.M.; Chruney, K.L.; Nuttall, R.L. The NBS tables of chemical thermodynamic properties. Selected values for inorganic and C1 and C2 organic substances in SI units. J. Phys. Chem. Ref. Data. 1982, 11(Suppl. 2). 47) Wolthers, M; van der Gaast, SJ; Rickard, D. The structure of disordered mackinawite. Amer. Min. 2003, 88, 2007-2015. 48) Jeong, HY; Lee, JH; Hayes, KF. Characterization of synthetic nanocrystalline mackinawite: Crystal structure, particle size, and specific surface area. Geochim. Cosmochim. Acta. 2008, 72, 493-505. 49) Lennie, A.R.; Redfern, S.A.T.; Schofield, P.F.; Vaughan, D.J. Synthesis and Rietveld crystal structure refinement of mackinawite, tetragonal FeS. Min. Mag. 1995, 59, 677-683. 50) Nissenbaum, A.; Swaine, D.J. Organic matter-metal interactions in recent sediments: The role of humic substances. Geochim. Cosmochim. Acta 1975, 40, 809-816. 51) Quagraine, E.K.; Georgakaki, I.; Coucouvanis, D. Reactivity and kinetic studies of (NH4)2(MoS4) in acidic aqueous solution: Possible relevance to the angiostatic function of the MoS42- ligand. J. Inorg. Biochem. 2009, 103, 143-155. 52) Weber, Th.; Muijsers, J.C.; Niemantsverdriet, J.W. (1995) Structure of amorphous MoS3. J. Phys. Chem. 99, 9194-9200. 53) Walton, R.I.; Dent, A.J.; Hibble, S.J. In situ investigation of the thermal decomposition of ammonium tetrathiomolybdate using combined time-resolved x-ray absorption spectroscopy and x-ray diffraction. Chem. Mat. 1998, 10, 3737-3745. 54) Hibble, SJ; Wood G.B. Modeling the structure of amorphous MoS3: A neutron diffraction and reverse Monte Carlo Study. J. Amer. Chem. Soc. 2003, 126, 959-965. 24 ACS Paragon Plus Environment

Page 25 of 52 1 2 651 3 652 4 653 5 6 654 7 655 8 656 9 657 10 11 658 12 659 13 660 14 661 15 16 662 17 663 18 664 19 665 20 21 666 22 667 23 668 24 669 25 26 670 27 671 28 672 29 673 30 31 674 32 675 33 676 34 677 35 678 36 37 679 38 680 39 681 40 682 41 42 683 43 684 44 685 45 686 46 47 687 48 688 49 689 50 690 51 52 691 53 692 54 693 55 694 56 57 695 58 696 59 60

ACS Earth and Space Chemistry

55) Vorlicek, T.P.; Kahn, M.D.; Kasuya, Y.; Helz, G.R. Capture of molybdenum in pyriteforming sediments: Role of ligand-induced reduction by polysulfides. Geochim. Cosmochim. Acta. 2004, 68, 547-556. 56) Davison, W.) The solubility of iron sulphides in synthetic and natural waters at ambient temperature. Aquatic Sci. 1991, 53, 309-329. 57) Kowalski, N; Dellwig O.; Beck M.; Gräwe U.; Neubert N.; Nägler T.F.; Badewien T.H.; Brumsack H.-J.; van Beusekom J.E.E.; Böttcher M.E. Pelagic molybdenum concentration anomalies and the impact of sediment resuspension on the molybdenum budget in two tidal systems of the North Sea. Geochim. Cosmochim. Acta 2013, 119, 198-211. 58) Morford, J.L; Emerson, S.R.; Breckel, E.J.; Kim, S.H. Diagenesis of oxyanions (V, U, Re, and Mo) in pore waters and sediments from a continental margin. Geochim. Cosmochim. Acta, 2005, 69, 5021-5032. 59) McManus, J.; Berelson, W.M.; Severmann, S.; Poulson, R.L.; Hammond, D.E.; Klinkhammer, G.P.; Holm C. Molybdenum and uranium geochemistry in continental margin sediments: Paleoproxy potential. Geochim. Cosmochim. Acta, 2006, 70, 4643-4662. 60) Raiswell, R.; Canfield, D.E. Sources of iron for pyrite formation in marine sediments. Amer. J. Sci. 1998, 298, 219-245. 61) Xu, J.; Qian, J.; Wei, Q. Synthesis and structure of a novel molybdenum-iron-sufur cluster with Mo2Fe2 core and all-disulfide chelate ligands, {Mo2Fe2(μ3-S)4(S2CNEt2)5}CH3CN. Inorg. Chim. Acta. 1989, 164, 55-58. 62) Kawaguchi, H.; Yamada, K.; Ohnishi, S.; Tatsumi, K. Construction of a cyclic tricubane cluster {Cp*2Mo2Fe2S4}(μ-S4)3 from the Mo2Fe2S4 single cubane component. J. Amer. Chem. Soc. 1997, 119, 10871-10872. 63) Koutmos, M.; Georgakaki, I.P.; Tsiolis, P.; Coucouvanis, D. Synthesis and characterization of MoFe3S4 and Mo2Fe2S4 single cubanes. Z. Anorg. Allg. Chem. 2008, 634, 255-261. 64) Melnik, M.; Ondrejkovicova, I.; Sundberg, M. R.; Holloway, C.E. Crystallographic and structural analysis of iron heterometallic compounds. Part III. Heterotetranuclear compounds. Rev. Inorg. Chem. 2008, 28, 237-456. 65) Shibahara, T.; Sakane, G.; Naruse, Y.; Taya, K.; Akashi, H., Ichimura, A.; Adachi, H. Syntheses and characterization of cubane-type clusters with molybdenum-iron-sulfur (Mo3FeS4) or molybdenum-nickel-sulfur (Mo3NiS4) cores. X-ray structures of {Mo3FeS4(H2O)10}(CH3C6H4SO3)4.7H2O and {Mo3FeS4(H2O)(NH3)9}Cl4, and discrete variational (DV)-Xα calculation of the electronic structures of {Mo3FeS4(H2O)10}4+, {Mo3FeS4(H2O)(NH3)9}4+ and {Mo3NiS4(H2O)10}4+. Bull. Chem. Soc. Japan. 1995, 68, 2769-2782. 66) Müller, A.; Hildebrand, A.; Penk, M.; Bögge, H. {Fe2S2(MoS4)2}4-: The first complex with two MoS42- ions coordinated to a ferredoxin type central unit. Inorg. Chim. Acta. 1988, 148,11-13. 67) Teo, B.-K.; Antonio, M.R.; Averill, B.A. Molybdenum K-edge extended x-ray absorption fine structure (EXAFS) studies of synthetic Mo-Fe-S clusters containing the MoS4 unit: Development of a fine adjustment technique based on models. J. Amer. Chem. Soc. 1983, 105, 3751-3762. 68) Bostick, B.C.; Fendorf, S.; Helz, G.R. Differential adsorption of molybdate and tetrathiomolybdate on pyrite (FeS2). Environ. Sci. Technol. 2003, 37, 285-291. 69) Diemann, E. Zur Struktur des natürlichen amorphen Molybdänsulfids (Jordisite) Naturwissen. 1976, 63, 385-386.

25 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 697 3 698 4 699 5 6 700 7 701 8 702 9 703 10 11 704 12 705 13 706 14 707 15 708 16 17 709 18 710 19 711 20 712 21 22 713 23 714 24 715 25 716 26 27 717 28 718 29 719 30 720 31 32 721 33 722 34 723 35 724 36 37 725 38 726 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 52

70) Cunningham, C.G.; Rasmussen, J.D.; Steven, T.A.; Rye, R.O.; Rowley, P.D.; Romberger, S.B.; Selverstone, J. Hydrothermal uranium deposits containing molybdenum and fluorite in the Marysvale volcanic field, west-central Utah. Mineralium Deposita 1998, 33, 477-494. 71) Granger, H.C.; Ingram, B.L. Occurrence and indentification of jordisite at Ambrosia Lake, New Mexico. U.S. Geol. Survey Prof. Paper. 1966, 550-B, B120-B124. 72) Cheney, E.S.; Trammell, J.W. Isotopic evidence for inorganic precipitation of uranium roll ore bodies. Amer. Assoc. Petrol. Geol. Bull. 1973, 57, 1297-1304. 73) Brookins, D.G. Position of uraninite and/or coffinite accumulations to the hematite/pyrite interface in sandstone-type deposits. Econ. Geol. 1976, 71, 944-948. 74) Hofmann, B.A.; Helfer, M.; Diamond, L.W.; Villa, I.M.; Frei, R.; Eikenberg, J. Topography-driven hydrothermal breccia mineralization of Pliocene age at Grimsel Pass, Aar massif, Central Swiss Alps. Schweiz. Mineral. Petrograph. Mitteilungen. 2004, 84, 271302. 75) Saunders, J.A.; Pivetz, B.E.; Voorhies, N.; Wilkin, R.T. Potential aquifer vulnerability in regions down-gradient from uranium in situ recovery (ISR) sites. J. Environ. Monit. 2016, 183, 67-83. 76) Clark, A.H. Molybdenite 2H1, molybdenite 3R, and jordisite from Carrizal Alto, Atacama, Chile. Amer. Min. 1971, 56, 1832-1835. 77) LeAnderson, J.P.; Schrader, E.L.; Brake, S.; Kaback, D.S. Behavior of molybdenum during weathering of the Ceresco Ridge porphyry molybdenite deposit, Climax, Colorado and a comparison with the Hollister deposit, North Carolina. Appl. Geochem. 1987, 2, 399-416. 78) Staples, L.W. Ilsemannite and jordistie. Amer. Min. 1951, 36, 609-614. 79) Kao, L-S.; Peacor, D.R.; Coveney, R. M. Jr.; Zhao, G.; Dungey K.E.; Curtis M.D.; PennerHahn J.E. A C/MoS2 mixed-layer phase (MoSC) occurring in metalliferous black shales from southern China, and new data on jordisite. Amer. Min. 2001, 86, 852-861. 80) Orberger, B.; Vymazalova, A.; Wagner, C.; Fialin, M.; Gallien, J.P.; Wirth, R.; Pasava, J.; Montagnac, G. Biogenic origin of intergrown Mo-sulphide- and carbonaceous matter in Lower Cambrian black shales (Zunyi Formation, southern China). Chem. Geol. 2007, 238, 213-231.

26 ACS Paragon Plus Environment

Page 27 of 52 1 2 727 3 728 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 729 26 730 27 731 28 732 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

Table 1. Experimental data. All runs contained 20 mM borate buffer. After about 100 h of reaction, 0.1 M MgCl2 was added to promote flocculation. Run

Initial ΣMo (M)

Initial ΣFe (M)

React. time (h)

Final ΣMo (M)

Final ΣFe (M)

Final pH

Final ΣS-II (M)

Final MoS42(M)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18*

6.62E-04 4.63E-04 3.31E-04 6.78E-05 6.81E-05 1.99E-04 1.32E-04 6.81E-05 6.62E-05 3.31E-05 6.81E-05 2.18E-05 6.78E-05 6.78E-05 6.78E-05 6.78E-05 0 6.78E-05

6.62E-05 6.62E-05 6.62E-05 1.36E-05 2.27E-05 6.62E-05 6.62E-05 3.40E-05 6.62E-05 6.62E-05 1.36E-04 6.62E-05 2.03E-04 3.39E-04 4.75E-04 6.76E-04 6.62E-05 6.76E-04

404 385 2328 355 811 2592 2760 980 3816 3168 1149 1656 979 810 356 356 1610 314

5.20E-04 4.28E-04 2.68E-04 6.50E-05 6.43E-05 1.61E-04 9.16E-05 6.31E-05 5.67E-05 2.00E-05 2.77E-05 1.96E-05 5.85E-05 4.86E-05 3.09E-05 4.05E-06 0 2.39E-05

5.28E-05 5.00E-05 4.96E-05 5.24E-06 1.30E-05 5.39E-05 3.02E-05 1.49E-05 3.10E-05 8.88E-06 1.35E-05 1.54E-05 4.57E-05 3.30E-05 2.91E-05 8.90E-07 2.45E-06 2.91E-05

8.35 8.33 8.37 8.22 8.07 8.23 8.23 8.23 8.38 8.18 8.25 8.31 8.03 8.07 8.08 8.03 8.33 8.01

8.76E-03 9.22E-03 7.41E-03 7.65E-03 4.66E-03 3.07E-03 6.67E-03 6.91E-03 4.61E-03 5.44E-03 6.69E-03 3.61E-03 7.21E-03 8.72E-03 6.56E-03 5.96E-03 6.77E-03 5.86E-03

5.11E-04 3.67E-04 2.27E-04 5.14E-05 4.68E-05 1.04E-04 5.80E-05 4.32E-05 1.93E-05 1.11E-05 1.17E-05 8.38E-06 1.26E-05 8.97E-06 4.42E-06 5.00E-07 0 4.90E-06



NFe/Mo is the initial mole ratio of total Fe to total Mo. * Sulfur saturated.

27 ACS Paragon Plus Environment

NFe/Mo‡ Initial 0.100 0.143 0.200 0.201 0.333 0.333 0.502 0.499 1.000 2.000 1.997 3.037 2.994 5.000 7.006 9.971 ∞ 9.971

ACS Earth and Space Chemistry 1 2 733 3 4 734 5 735 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 736 28 29 737 30 31 738 32 739 33 34 35 740 36 741 37 742 38 39 40 41 42 43 44 45 46 47 48 49 50 743 51 744 52 53 54 55 56 57 58 59 60

Page 28 of 52

Table 2. Results from thermodynamic speciation calculations. Run

HS- (M)

Fe2+ (M)

Fe(HS)20 (M)

[(FeS)2(MoS4)2]4- (M)

ΣMo calc (M)

log QFeS

log QMoS3

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18

6.47E-03 7.22E-03 6.12E-03 7.08E-03 4.14E-03 2.27E-03 6.07E-03 6.38E-03 4.23E-03 5.12E-03 6.32E-03 3.41E-03 6.52E-03 8.02E-03 6.06E-03 5.58E-03 6.56E-03 4.58E-03

2.47E-11 3.16E-11 5.42E-11 9.40E-11 4.19E-10 4.54E-10 2.32E-10 2.07E-10 6.95E-10 8.66E-10 6.99E-10 1.62E-09 1.98E-09 1.74E-09 4.27E-09 2.55E-09 1.44E-07 6.00E-09

3.84E-10 6.09E-10 7.52E-10 1.75E-09 2.67E-09 8.68E-10 3.18E-09 3.12E-09 4.60E-09 8.40E-09 1.03E-08 6.96E-09 3.11E-08 4.14E-08 5.81E-08 2.95E-08 2.29E-06 4.66E-08

2.52E-05 2.40E-05 2.35E-05 2.43E-06 6.24E-06 2.50E-05 1.45E-05 7.07E-06 1.39E-05 4.15E-06 6.35E-06 6.70E-06 2.26E-05 1.62E-05 1.42E-05 3.89E-07 1.43E-05

5.98E-04 4.38E-04 2.93E-04 5.90E-05 6.00E-05 1.73E-04 9.13E-05 6.03E-05 5.25E-05 2.06E-05 2.58E-05 2.47E-05 5.86E-05 4.21E-05 3.36E-05 2.40E-06 3.42E-05

1.80 1.93 2.14 2.29 2.55 2.49 2.62 2.59 3.09 2.87 3.13 3.30 3.38 3.45 3.73 3.43 5.55 3.69

-16.87 -17.04 -17.22 -17.78 -17.46 -16.99 -17.67 -17.82 -18.14 -17.92 -18.4 -18.34 -18.17 -18.44 -18.64 -19.03 -18.4

Table 3: Mo average oxidation state for the precipitates derived from the relationship between Em and the Mo oxidation state. See Fig.7. Samples Reference materials Mo(IV)S2 Mo(VI)S42Runs 6 7 9 12

Em (eV)

Average Mo oxidation state

2.0030.0 20043.0

4.0 6.0

20030.7 20030.1 20030.2 20030.5

4.1 4.0 4.0 4.1

28 ACS Paragon Plus Environment

Page 29 of 52 1 2 745 3 4 5 6 7 8 9 10 11 12 13 14 746 15 16 17 747 18 19 748 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 749 39 750 40 751 41 42 752 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

Table 4. Molybdenum EXAFS data for selected precipitates. Runs

Mo--S1

Mo--S2

Mo--Fe

Amp

R factor

n

r (Å)

σ2

n

r (Å)

σ2

n

r (Å)

σ2

6

2

2.32±0.02

0.003

3

2.44±0.03

0.003

1

2.79±0.03

0.003

0.85

0.02

7

2

2.32±0.05

0.003

3

2.47±0.05

0.003

1

2.81±0.04

0.003

0.93

0.02

9

2

2.30±0.03

0.003

3

2.45±0.04

0.003

1

2.80±0.02

0.003

0.97

0.03

12

2

2.31±0.03

0.003

3

2.46±0.03

0.003

1

2.81±0.05

0.003

0.88

0.02

All σ2 values were set to 0.003

Table 5. Atomic proportions of Fe, Mo and S in some jordisites and Fe-Mo-S precipitates. Sample Type

Empirical Formula

Reference

Jordisites

Fe0.61MoS3.21 Fe0.28MoS2.24 Fe0.39MoS2.56 Fe0.48MoS2.86 Fe0.40MoS2.70 Fe0.31MoS2.27

78 71 71 71 71 74

Cambrian ore mineral*

(Fe,Ni)0.33MoS2.67C10*

79

Precipitates

Fe0.78MoS3.83 This work Fe0.83MoS3.75 This work Fe0.63MoS3.60 This work *An exotic graphite intercolation compound that occurs as an ore mineral in metalliferous black shales in China.79,80

29 ACS Paragon Plus Environment

ACS Earth and Space Chemistry

0.0006 0.0005 0.0004

MoCalculated (M)

1 2 753 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 754 22 23 755 24 756 25 26 757 27 758 28 759 29 760 30 761 31 762 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 763 51 52 764 53 54 55 56 57 58 59 60

Page 30 of 52

-6

E .5 -6

0.0003 0.0002 =

8 99 0.

M

o obs

c

M

0.0001 0.0000 0.0000

0.0001

o cal

0.0002

0.0003

0.0004

0.0005

0.0006

Moobserved (M)

Figure 1. Excellent agreement of calculated total dissolved Mo with measured total dissolved Mo. Mocalc is determined by summing the calculated concentrations of the dissolved Mo species (i.e. MoS42- + MoOS32- + [FeO(OH)MoS4]3- + 2[(FeS)2(MoS4)2]4-). The point represented by a triangle at the highest dissolved Mo concentration deviates from the trend of the remaining data and was excluded from the regression.

Figure 2. Typical time course for a run (here, Run 5; see Table 1).

30 ACS Paragon Plus Environment

15

10 (Run 16) 0.5 (Run 7)

Initial NFe/Mo = 10

2-

5

0 0

25

50

75

100

Hours precipitate aged before EDTA

Figure 3. Reversal of the precipitation process by 0.001 M EDTA after various aging times. Beyond the first day, most of the precipitated MoS42- is no longer released by this treatment.

g lo

-17

oS

QM

=

3

-l og 4.

-1

eS

-18

95 ± 0. 06

log QMoS3

QF

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 765 21 22 766 23 24 767 25 768 26 769 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 770 47 48 771 49 772 50 773 51 774 52 775 53 776 54 55 777 56 57 58 59 60

ACS Earth and Space Chemistry

MoS4 dissolved by EDTA (as %Mo precipitated)

Page 31 of 52

-19

1

2

3

4

5

log QFeS

Figure 4. Gibbs-Duhem plot relating the chemical potentials of MoS3 and FeS components (data in Table 2). Error bars display analytical uncertainties that have been propagated through the thermodynamic speciation calculations by the root-mean-square method. The line is an errorweighted regression line with slope forced to -1.0 (R2 0.95). The ±0.06 uncertainty in the intercept is the standard error.

31 ACS Paragon Plus Environment

ACS Earth and Space Chemistry

100

101

111

112

Counts (arbitrary units)

1 2 778 3 779 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 780 27 781 28 29 782 30 31 783 32 784 33 785 34 35 786 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 52

Fe/Mo 0.35/0

0.35/0.07 0.07/0.21

10

20

30

40

50

60

2O

Figure 5. X-ray diffraction patterns (Cu Kα radiation) for three freeze-dried products (polynomial background correction). Labels at the right side of the figure give initial concentrations (mM) of Fe2+ and MoS42- in the solutions that produced the precipitates. Dotted vertical lines show positions of major peaks for crystalline mackinawite.49

32 ACS Paragon Plus Environment

Page 33 of 52

30

a 0.70 mM Fe

20 10 0 10

b 0.70 mM Fe, 0.07 mM Mo

5 0 10

Number of cases

1 2 787 3 788 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 789 42 790 43 791 44 792 45 793 46 794 47 795 48 49 796 50 797 51 798 52 799 53 800 54 801 55 56 802 57 803 58 59 60

ACS Earth and Space Chemistry

c 0.35 mM Fe, 0.07 mM Mo

5 0 15

d 0.07 mM Fe, 0.07 mM Mo

10 5 0 30

e 0.07 mM Fe, 0.35 mM Mo

20 10 0 30

f 0.07 mM Fe, 0.70mM Mo

20 10 0 0.0

0.2

0.4

0.6

[Fe counts]/[Fe counts + (Mo+S) counts]

Figure 6. Histograms showing the ratio of EDS counts for Fe to the sum of counts for Fe plus overlapping counts for Mo + S in test precipitates. Counts for Mo and S have nearly the same energy and could not be measured independently. Test precipitates were taken from mother solutions initially containing 10 mM S2-and 20 mM borate buffer at pH 8.5 as well as (panel a) 0.70 mM Fe2+, (panel b) 0.70 mM Fe2+ and 0.07 mM MoS42-, (panel c) 0.35 mM Fe2+ and 0.07 mM MoS42-, (panel d) 0.07 mM Fe2+ and 0.07 mM MoS42-, (panel e) 0.07 mM Fe2+ and 0.35 mM MoS42-, and (panel f) 0.07 mM Fe2+ and 0.70 mM MoS42-.

33 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 804 36 805 37 806 38 807 39 808 40 809 41 42 810 43 811 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 52

Figure 7: XANES spectra of samples (solid lines) and Mo reference compounds (dashed lines). The black dots identify the maximum absorption energy (Em).

34 ACS Paragon Plus Environment

Page 35 of 52

0 25 50 75 100 125

C

B

A

ambient logQFeS

logQFeS= 4.87

90 m 150 m

500

Depth (m)

1 2 812 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 813 32 814 33 815 34 816 35 817 36 818 37 819 38 39 820 40 821 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

1000

1500

2-

{MoS4 } (1) (2) (3)

2000 2+

Total dissolved Mo

{H2S}

{Fe } 2500 -9

-8

-7

-6

-5

log activity

-4

-9

-8

-7

-6

-14

log concentration

-12

-10

-8

-6

log activity

Figure 8. Black Sea data. Panel A: Profiles of thermodynamic activities of Fe2+ and H2S in the sulfidic waters of the Black Sea calculated from water analyses. Panel B: Observed concentration of dissolved Mo from water analyses. Panel C: Curve 1, {MoS42-}BS calculated from water analyses. Curve 2, {MoS42─}eq at equilibrium with FeMoS4 if log QFeS at each depth is calculated from water analyses. Curve 3, {MoS42-}eq if log QFeS is fixed at 4.9. Activities and concentrations on molar scale. See Table S1 for sources of data.

35 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 822 3 4 823 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 824 28 29 825 30 826 31 827 32 828 33 829 34 35 830 36 831 37 832 38 833 39 834 40 835 41 42 836 43 837 44 838 45 839 46 840 47 841 48 842 49 50 843 51 844 52 845 53 846 54 847 55 848 56 849 57 58 850 59 851 60

Page 36 of 52

Figure 9. Left: Structure of an idealized aqueous cuboidal cluster cation, [FeMo3S4(H2O)10]4+;65 actual heterometallic cuboid structures are distorted from perfect cubic symmetry. Right: hypothetical monomer produced in this work. Mo-S interatomic distances in the [FeMo3S4(H2O)10]4+ cation range from 2.344 to 2.356 Å compared to averages of 2.31 and 2.46 Å in our precipitates; Mo-Fe distances in [FeMo3S4(H2O)10]4+ are 2.660 to 2.681 Å but average 2.80 Å in our precipitates. The cation has Mo-Mo distances of 2.77 Å, but a Mo-Mo shell was not resolved in our amorphous precipitates, possibly due to disorder.

36 ACS Paragon Plus Environment

Page 37 of 52 1 2 852 3 853 4 854 5 855 6 856 7 8 857 9 858 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 859 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

Table of Contents Art

37 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 38 of 52

Run

Initial ΣMo (M)

Initial ΣFe (M)

React. time (h)

Final ΣMo (M)

Final ΣFe (M)

Final pH

Final ΣS-II (M)

Final MoS42(M)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18*

6.62E-04 4.63E-04 3.31E-04 6.78E-05 6.81E-05 1.99E-04 1.32E-04 6.81E-05 6.62E-05 3.31E-05 6.81E-05 2.18E-05 6.78E-05 6.78E-05 6.78E-05 6.78E-05 0 6.78E-05

6.62E-05 6.62E-05 6.62E-05 1.36E-05 2.27E-05 6.62E-05 6.62E-05 3.40E-05 6.62E-05 6.62E-05 1.36E-04 6.62E-05 2.03E-04 3.39E-04 4.75E-04 6.76E-04 6.62E-05 6.76E-04

404 385 2328 355 811 2592 2760 980 3816 3168 1149 1656 979 810 356 356 1610 314

5.20E-04 4.28E-04 2.68E-04 6.50E-05 6.43E-05 1.61E-04 9.16E-05 6.31E-05 5.67E-05 2.00E-05 2.77E-05 1.96E-05 5.85E-05 4.86E-05 3.09E-05 4.05E-06 0 2.39E-05

5.28E-05 5.00E-05 4.96E-05 5.24E-06 1.30E-05 5.39E-05 3.02E-05 1.49E-05 3.10E-05 8.88E-06 1.35E-05 1.54E-05 4.57E-05 3.30E-05 2.91E-05 8.90E-07 2.45E-06 2.91E-05

8.35 8.33 8.37 8.22 8.07 8.23 8.23 8.23 8.38 8.18 8.25 8.31 8.03 8.07 8.08 8.03 8.33 8.01

8.76E-03 9.22E-03 7.41E-03 7.65E-03 4.66E-03 3.07E-03 6.67E-03 6.91E-03 4.61E-03 5.44E-03 6.69E-03 3.61E-03 7.21E-03 8.72E-03 6.56E-03 5.96E-03 6.77E-03 5.86E-03

5.11E-04 3.67E-04 2.27E-04 5.14E-05 4.68E-05 1.04E-04 5.80E-05 4.32E-05 1.93E-05 1.11E-05 1.17E-05 8.38E-06 1.26E-05 8.97E-06 4.42E-06 5.00E-07 0 4.90E-06

Table 1

ACS Paragon Plus Environment

NFe/Mo‡ Initial 0.100 0.143 0.200 0.201 0.333 0.333 0.502 0.499 1.000 2.000 1.997 3.037 2.994 5.000 7.006 9.971 ∞ 9.971

Page 39 of 52 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

Run

HS- (M)

Fe2+ (M)

Fe(HS)20 (M)

[(FeS)2(MoS4)2]4- (M)

ΣMo calc (M)

log QFeS

log QMoS3

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18

6.47E-03 7.22E-03 6.12E-03 7.08E-03 4.14E-03 2.27E-03 6.07E-03 6.38E-03 4.23E-03 5.12E-03 6.32E-03 3.41E-03 6.52E-03 8.02E-03 6.06E-03 5.58E-03 6.56E-03 4.58E-03

2.47E-11 3.16E-11 5.42E-11 9.40E-11 4.19E-10 4.54E-10 2.32E-10 2.07E-10 6.95E-10 8.66E-10 6.99E-10 1.62E-09 1.98E-09 1.74E-09 4.27E-09 2.55E-09 1.44E-07 6.00E-09

3.84E-10 6.09E-10 7.52E-10 1.75E-09 2.67E-09 8.68E-10 3.18E-09 3.12E-09 4.60E-09 8.40E-09 1.03E-08 6.96E-09 3.11E-08 4.14E-08 5.81E-08 2.95E-08 2.29E-06 4.66E-08

2.52E-05 2.40E-05 2.35E-05 2.43E-06 6.24E-06 2.50E-05 1.45E-05 7.07E-06 1.39E-05 4.15E-06 6.35E-06 6.70E-06 2.26E-05 1.62E-05 1.42E-05 3.89E-07 0.00E+00 1.43E-05

5.98E-04 4.38E-04 2.93E-04 5.90E-05 6.00E-05 1.73E-04 9.13E-05 6.03E-05 5.25E-05 2.06E-05 2.58E-05 2.47E-05 5.86E-05 4.21E-05 3.36E-05 2.40E-06 3.42E-05

1.80 1.93 2.14 2.29 2.55 2.49 2.62 2.59 3.09 2.87 3.13 3.30 3.38 3.45 3.73 3.43 5.55 3.69

-16.87 -17.04 -17.22 -17.78 -17.46 -16.99 -17.67 -17.82 -18.14 -17.92 -18.4 -18.34 -18.17 -18.44 -18.64 -19.03 -18.4

Table 2

ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Samples Reference materials Mo(IV)S2 Mo(VI)S42Runs 6 7 9 12

Page 40 of 52

Em (eV)

Average Mo oxidation state

20030.0 20043.0

4.0 6.0

20030.7 20030.1 20030.2 20030.5

4.1 4.0 4.0 4.1

Table 3

ACS Paragon Plus Environment

Page 41 of 52 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

Runs

Mo--S1

Mo--S2

Mo--Fe

Amp

R factor

n

r (Å)

σ2

n

r (Å)

σ2

n

r (Å)

σ2

6

2

2.32±0.02

0.003

3

2.44±0.03

0.003

1

2.79±0.03

0.003

0.85

0.02

7

2

2.32±0.05

0.003

3

2.47±0.05

0.003

1

2.81±0.04

0.003

0.93

0.02

9

2

2.30±0.03

0.003

3

2.45±0.04

0.003

1

2.80±0.02

0.003

0.97

0.03

12

2

2.31±0.03

0.003

3

2.46±0.03

0.003

1

2.81±0.05

0.003

0.88

0.02

All σ2 value were set to 0.003

Table 4

ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 42 of 52

Sample Type

Empirical Formula

Reference

Jordisites

Fe0.61MoS3.21 Fe0.28MoS2.24 Fe0.39MoS2.56 Fe0.48MoS2.86 Fe0.40MoS2.70 Fe0.31MoS2.27

73 66 66 66 66 69

Cambrian ore mineral*

(Fe,Ni)0.33MoS2.67C10*

74

Precipitates

Fe0.78MoS3.83 This work Fe0.83MoS3.75 This work Fe0.63MoS3.60 This work *An exotic graphite intercolation compound that occurs as an ore mineral in metalliferous black shales in China.74,75

Table 5

ACS Paragon Plus Environment

Page 43 of 52

0.0006 0.0005 0.0004

MoCalculated (M)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

0.0003 0.0002 =

99 0.

8M

o obs

5E 6.

-6

c

M

0.0001 0.0000 0.0000

0.0001

o c al

0.0002

0.0003

0.0004

Moobserved (M)

Figure 1

ACS Paragon Plus Environment

0.0005

0.0006

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2

ACS Paragon Plus Environment

Page 44 of 52

15

10 (Run 16) 0.5 (Run 7)

Initial NFe/Mo = 10

5

2-

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41

ACS Earth and Space Chemistry

MoS4 dissolved by EDTA (as %Mo precipitated)

Page 45 of 52

0 0

25

50

75

Hours precipitate aged before EDTA ACS Paragon Plus Environment

100

ACS Earth and Space Chemistry

g lo

-17

=

3 oS

QM

og

-l 9 4.

-1

S

-18

5 ± 6

0 0.

log QMoS3

Q Fe

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 46 of 52

-19

1

2

3

log QFeS

Figure 4

ACS Paragon Plus Environment

4

5

Page 47 of 52

100

101

111

112

Counts (arbitrary units)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

Fe/Mo 0.35/0

0.35/0.07 0.07/0.21

10

20

30

40

50

2O

Figure 5

ACS Paragon Plus Environment

60

ACS Earth and Space Chemistry

30

a 0.70 mM Fe

20 10 0 10

b 0.70 mM Fe, 0.07 mM Mo

5 0 10

Number of cases

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 48 of 52

c 0.35 mM Fe, 0.07 mM Mo

5 0 15

d 0.07 mM Fe, 0.07 mM Mo

10 5 0 30

e 0.07 mM Fe, 0.35 mM Mo

20 10 0 30

f 0.07 mM Fe, 0.70mM Mo

20 10 0 0.0

0.2

0.4

[Fe counts]/[Fe counts + (Mo+S) counts]

Figure 6

ACS Paragon Plus Environment

0.6

Page 49 of 52 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

Figure 7

ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 8

ACS Paragon Plus Environment

Page 50 of 52

Page 51 of 52 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

Figure 9

ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC Figure

ACS Paragon Plus Environment

Page 52 of 52