Rapid Iron Reduction Rates Are Stimulated by High-Amplitude Redox

Feb 13, 2017 - Nicholas K. Masters-Waage , Katherine Morris , Jonathan R. Lloyd , Samuel Shaw , J. Frederick W. Mosselmans , Christopher Boothman , Pi...
0 downloads 20 Views 3MB Size
Subscriber access provided by University of Newcastle, Australia

Article

Rapid iron reduction rates are stimulated by highamplitude redox fluctuations in a tropical forest soil Brian R. Ginn, Christof Meile, Jared Wilmoth, Yuanzhi Tang, and Aaron Thompson Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b05709 • Publication Date (Web): 13 Feb 2017 Downloaded from http://pubs.acs.org on February 19, 2017

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 free 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 accessible to all readers and 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.

Environmental Science & Technology 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 29

Environmental Science & Technology

1 2 3

Rapid iron reduction rates are stimulated by high-amplitude redox fluctuations in a tropical

4

forest soil

5

Brian Ginn1, Christof Meile2, Jared Wilmoth1, Yuanzhi Tang3 and Aaron Thompson1,*

6 7

1

2

8 9

University of Georgia (Crop and Soil Sci.)

3

University of Georgia (Marine Sci.)

Georgia Institute of Technology (Earth and Atmospheric Sci)

10 11 12

For submission to Environmental Science and Technology

13 14 15

*Corresponding author E-mail: [email protected], (01) 706-410-1293

16

1

ACS Paragon Plus Environment

Environmental Science & Technology

Page 2 of 29

17 18 19

Abstract:

20

altered by redox-driven processes. This study examined the influence of temporal oxygen variations on

21

Fe speciation in soils from the Luquillo Critical Zone Observatory (Puerto Rico). We incubated soils

22

under cycles of oxic-anoxic conditions (τoxic:τanoxic = 1:6) at three frequencies with and without

23

phosphate addition. Fe(II) production, P availability, and Fe mineral composition were monitored using

24

batch analytical and spectroscopic techniques. The rate of FeII reduction increased from ~ 3 to > 45

25

mmol Fe(II) kg-1 soil d-1 over the experiment with a concomitant increase of an Fe(II) concentration

26

plateau within each anoxic period. The apparent maximum in reducible Fe is similar in all treatments,

27

but was hastened by P-amendment. Numerical modelling suggests the Fe(II) dynamics can be

28

explained by the formation of rapidly reducible Fe(III) phases derived from the progressive dissolution

29

of native Fe(III) oxides accompanied by minor increases in Fe reducer populations. The shift in Fe(III)

30

reactivity is evident from Fe-reducibility assays using Shewanella sp., however was undetectable by

31

chemical extractions, Mössbauer or X-ray Absorption spectroscopies. More broadly, our findings

32

suggest Fe reduction rates are strongly coupled to redox dynamics of the recent past, and that frequent

33

shifts in redox conditions can prime a soil for rapid Fe-reduction.

Iron oxides are important structural and biogeochemical components of soils that can be strongly

Rapidly&Reducible&Fe(III)& Reducible&Fe(III)&

Oxida1on&

34

Fe(II)&

Reduc1on&

Reduc1on&

TOC ART

2

ACS Paragon Plus Environment

Page 3 of 29

35

Environmental Science & Technology

Introduction

36 37

Although iron comprises only 1-3% of the total mass in soils1, it has a profound influence on

38

the properties of soils. The high reactivity and large surface area of Fe-(oxyhydr)oxide minerals

39

influence many biogeochemical processes through the sorption of key nutrients or contaminants2-10.

40

When exposed to low oxygen conditions, Fe minerals serve as important electron acceptors for

41

microbial respiration11, leading to the dissolution of solid phases and release of sorbed and

42

incorporated constituents12, 13. This reductive dissolution of Fe has significant implications for the

43

global C cycle in soils where a substantial fraction of the total soil organic C oxidation is coupled to the

44

microbial reduction of Fe-(oxyhydr)oxides (e.g., 44% in humid tropical forest soils14), which can occur

45

in upland soils with Fe(II) production rates ranging from 0.5 to 50 mmol kg-1 d-1 15, 16. The availability

46

of Fe minerals toward reduction is strongly correlated with their surface area and particle size17, 18. The

47

crystallinity of Fe-(oxyhydr)oxides also strongly influences the kinetics of microbial Fe(III)

48

reduction19. Microbes can rapidly reduce ferrihydrite and other short-range ordered (SRO) Fe minerals

49

(typically within hours), but reduce well-ordered hematite (α-Fe2O3), goethite (α-FeOOH), and

50

lepidocrocite (γ-FeOOH) at slower rates (e.g., several months)18, 20.

51

In upland soils, Fe reduction takes place within microsites where high-levels of organic matter

52

combine with ephemeral water saturation and low O2 diffusion rates to yield anoxic conditions21, 22.

53

Since O2 often re-infuses such microsites as the soils dry, Fe(II) produced during the anoxic period via

54

microbial Fe reduction is often re-oxidized to Fe(III) solid phases. This Fe redox cycling can influence

55

the rapid cycling of other elements (e.g., C, N, and metals23). Redox fluctuations can drive repeated

56

dissolution and precipitation of Fe (oxyhydr)oxide minerals with potential influence on the crystallinity

57

of those phases21, 24-27. During anoxic periods, the sorption of the Fe2+(aq) ions to Fe mineral surfaces

3

ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 29

58

can accelerate the process of Ostwald ripening from SRO to secondary crystalline mineral phases (e.g.,

59

goethite, magnetite) that are less chemically reactive and bioavailable than ferrihydrite26, 28, 29.

60

However, there are also studies illustrating SRO phases form or persist in fluctuating redox

61

environments30-34, indicating that the dynamics of the fluctuations are likely important in determining

62

the trajectory of mineral transformations.

63

Changes in the crystallinity of Fe (oxyhydr)oxides may also impact the fate of phosphorous (P)

64

in natural systems because Fe (oxyhydr)oxide surfaces have a great capacity for adsorbing P4, 5.

65

Ferrihydrite has a significantly higher sorption capacity than the crystalline (oxyhydr)oxides because of

66

its large surface area, but can form weaker surface complexes with aqueous species than the more

67

crystalline iron phases35-37. As a consequence, microbial Fe(III) reduction may simultaneously reduce

68

the overall capacity of soils to adsorb metal/metalloids while strengthening the remaining sorbate-

69

sorbent interactions.

70

Given the potential importance of Fe redox cycling for soil systems, the influence of redox

71

oscillations on Fe reduction rates needs to be quantified. A redox fluctuation can be characterized by

72

(a) the frequency of the redox oscillation; (b) the oscillation amplitude; and (c) the time spent under

73

oxic (τoxic) or anoxic (τanoxic) conditions. Here we hypothesize that shorter frequency oscillations will

74

lead to faster Fe reduction rates with associated changes in Fe-(oxyhydr)oxide mineral composition and

75

the distribution of phosphate. We test this by exposing a tropical soil to different redox oscillation

76

frequencies while holding the oscillation amplitude and ratio of τoxic: τanoxic constant. In order to

77

minimize spatial heterogeneity, soils were continuously mixed in slurries—separating the effects of

78

redox fluctuations from transport to the extent possible. We then measure changes in HCl-extractable

79

Fe(II), Fe mineral composition and P distribution, and simulate the experimental data in a numerical

80

model to calculate Fe reduction rates.

4

ACS Paragon Plus Environment

Page 5 of 29

Environmental Science & Technology

81 82

Methods

83 84 85

Field site and soil preparation Intact soil blocks were collected from the upper 10 cm of soil at a site similar to that described

86

by Peretyazhko and Sposito 38 in the Bisley watershed of the Luquillo Experimental Forest (LEF)

87

Puerto Rico (approximate GPS location 18.31553N,-65.74574W), which is part of the Luquillo Long-

88

term Ecological Research site and the Critical Zone Observatory. Samples were placed in 4-mil

89

polyethylene ziplock bags and transported without cooling inside an insulated container to the

90

University of Georgia. Within 24 h of sampling, all moist soils were processed by air-drying at 30°C

91

for 24 h followed by 2-mm sieving/homogenization (to remove large roots and gravel) under field soil

92

moisture conditions within a 95%:5% N2:H2 glovebox (Coy). All soils were stored at 22°C under oxic

93

conditions in the dark. Ginn et al. 39 showed these soils maintained the highest levels of microbial Fe-

94

reduction when air-dried prior to long-term storage. Soils in the LEF are classified predominately as

95

Ultisols formed from volcanic parent material with quartz diorite intrusions40, 41. The major primary

96

minerals consist of quartz and plagioclase feldspars with lesser amounts of biotite, hornblende, K-

97

feldspar, and accessory magnetite, sphene, apatite, and zircon42. Prior work by Silver et al.21 has

98

identified periodic fluctuations in soil O2 concentrations sufficient to generate anoxic conditions within

99

the soil. Further information on this soil collection and characterization are described elsewhere39, 43.

100 101

Redox oscillation experiments

102

Three oscillation treatments were performed simultaneously. In each experiment, 4.5 g (dry-

103

weight) of live soil (not sterilized) was suspended in a buffer solution containing 2 mM KCl and 10

104

mM solution of 2-(N-morpholino)ethanesulfonic acid (MES) buffer with a 45 g final suspension mass. 5

ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 29

105

The suspension was supplemented with 200 mg kg-1 of P (i.e. 6.46 mmol kg-1 soil and 0.58 mmol L-1

106

suspension). The slurries were incubated for 56 d in a glovebox (95% N2, 5% H2) on an end-over-end

107

shaker at 250 rpm. Although these experiments likely had sufficient labile organic carbon to support Fe

108

reduction rates (see below), the presence of H2 in the glovebox may have served as an additional

109

electron donor for Fe reduction44. In each of the experiments, the soil slurries were exposed to

110

alternating oxic-anoxic conditions at frequencies of 3.5, 7, and 14 d with the ratio of time under

111

oxic/anoxic conditions (τoxic:τanoxic) maintained at 1:6 for all treatments (i.e., the 3.5 d oscillation had 3

112

d of anoxic time and 0.5 d of oxic time; the 7 d period had 6 d anoxic and 1 d of oxic; and the 14 d

113

period had 12 d anoxic and 2 d oxic time). Samples were transferred between shaker tables inside and

114

outside of the glovebox to create anoxic and oxic conditions, respectively. Control soil slurries were

115

maintained at constant anoxic or oxic conditions for 56 d.

116

Samples were extracted at the beginning and end of each oscillation cycle using wide-orifice

117

pipette tips that allowed complete collection of soil particles in the slurry. Aqueous Fe(II) was

118

extracted from the soil slurries by centrifuging the samples at 11,000 rcf for 30 min, which will remove

119

particles > 40 nm based on Stokes’ law with an assumed mean particle density of 2.72 g cm-3 for

120

basalt-influenced soils45. Acid-extractable Fe(II) was solubilized by suspending the remaining pellet in

121

1.0 mL of 0.5 M HCl and shaking it for 2 h on a horizontal shaker at 150 rpm. The extracts were then

122

centrifuged at 11,000 rcf for 30 min and the supernatants analyzed for Fe(II) using a modified ferrozine

123

protocol24 on a Shimadzu-1700 UV-vis spectrophotometer at λ=562 nm. The dry weight of each pellet

124

from the centrifugation step was measured to correct for minor differences in suspension density

125

between treatments and allow accurate dry-weight normalization. All values are reported as averages of

126

triplicate measurements with error bars indicating standard deviation.

127 128

FeIII-reduction rate experiment 6

ACS Paragon Plus Environment

Page 7 of 29

129

Environmental Science & Technology

To quantify Fe(III) reduction rates midway through the experiment, additional samples were

130

oscillated at the three frequencies given above for 28 d and then exposed to 7 d of continuously anoxic

131

conditions during which they were sampled for aqueous and acid-extractable Fe(II) every 0.5 d as

132

described above. This data captured the iron reduction dynamics during a reducing period and was used

133

to parameterize the reduction kinetics in the numerical model described below.

134 135 136

Citrate-ascorbate extractable Fe To quantify the short-range-ordered solid Fe phases in the samples, we conducted citrate-

137

ascorbate extractions46 on freeze-dried soils from the unreacted samples or at the end of the oscillation

138

experiment (freeze-dried immediately upon completion of the experiment). Briefly, the samples were

139

suspended in a solution containing 0.2 M Na-citrate and 0.05 M ascorbic acid (pH = 6) at a 60:1

140

reagent:sample weight ratio. The suspensions were shaken for 16 h on a horizontal shaker and

141

centrifuged at 11,000 rcf for 30 min. The supernatants were analyzed for total Fe on an atomic

142

absorption (AA) spectrometer (Perkin Elmer, AAnalyst 200).

143 144

Mössbauer spectroscopy and X-ray absorption spectroscopy

145

Fe speciation of the soil samples was characterized via 57Fe Mössbauer spectroscopy and X-ray

146

absorption spectroscopy (XAS). Details on data collection, analysis, and interpretation are given in the

147

Supporting Information (SI) Sections 2 and 3.

148 149 150

Fe reduction assay: Shewanella incubations We quantified the availability of Fe(III) for microbial reduction by incubating the soil with a

151

culture of Shewanella oneidensis MR-1. The purpose of this inoculation experiment was to determine

152

whether redox fluctuations altered the iron oxide availability for microbial Fe-reduction. S. oneidensis 7

ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 29

153

MR-1 was used because it is a good model Fe-reducers, capable of fast Fe-reduction and can live in

154

both oxic and anoxic conditions. This experiment is performed solely to assess Fe availability across

155

treatments, and is not designed to quantify the in situ iron reduction rate. The assay was performed on

156

each replicate of soils freeze-dried from the end of the three different oscillation frequencies, the

157

anaerobic and aerobic controls, and unreacted soil. Post-oxidation, freeze-dried soil was used to

158

maintain consistency across treatments. Shewanella oneidensis MR-1 was grown to stationary phase in

159

media containing (per liter): 0.5 g KH2PO4, 1.0 g Na2SO4, 2.0 g NH4Cl, 1.0 g yeast extract, 0.5 mM

160

CaCl2, 0.1 mM MgSO4, 10 mM Na-lactate, and 50 mM Fe-citrate. The cells were washed by

161

centrifugation at 3,000 rcf for 30 min and resuspended in fresh media without yeast extract. The

162

washing step was repeated twice to thoroughly remove spent media, after which 0.9 g of washed cell

163

culture was added to 2-mL centrifuge tubes containing 0.1 g of freeze-dried soil so each sample

164

contained the same ratio of bacterial culture to soil mass. The cultures were incubated for a week under

165

anoxic conditions. Freeze-dried soils with no S. oneidensis inoculation were used as controls. The

166

cultures were sampled for acid-extractable Fe(II) by pipetting 0.1 mL of sample from each culture and

167

centrifuging at 11,000 rcf for 30 min. The supernatants were poured off and the pellets were re-

168

suspended in 0.1 mL of 0.5 M HCl. The suspensions were shaken for 2 h on a horizontal shaker. The

169

extracts were then centrifuged at 11,000 rcf for 30 min. Samples were analyzed for acid-extractible

170

Fe(II) using the ferrozine protocol described above.

171 172 173

Phosphate extractions Sequential extractions were performed according to Hedley et al.47. A soil mass of 0.01 g was

174

sequentially extracted for 16 h (each) on a horizontal shaker with 0.6 mL of ultrapure water, 0.5 M

175

NaHCO3, 0.1 M NaOH, and 1 M HCl. After each extraction, the tubes were centrifuged at 11,000 rcf

176

for 30 min. In addition, we quantified microbial biomass P following Hedley and Stewert48. Briefly, 1 8

ACS Paragon Plus Environment

Page 9 of 29

Environmental Science & Technology

177

mL of CHCl3 was added to 0.01 g of soil in a centrifuge tube and capped. The solution was shaken

178

periodically over 30 min and then uncapped so that the CHCl3 could evaporate overnight. In all cases,

179

supernatants were analyzed for P using a molybdate blue protocol49 on a Shimadzu-1700

180

spectrophotometer at λ = 880 nm.

181 182

Phosphate adsorption isotherms

183

To assess the adsorption capacity of the soil, we constructed phosphate adsorption isotherms on

184

unamended soils. Briefly, 1 mL of 0, 80, 160, 320, and 640 ppm phosphate was added to 0.1 g of each

185

freeze-dried soil in a centrifuge tube. The suspensions were shaken on a horizontal shaker for 2 h and

186

then centrifuged for 30 min at 11,000 rcf. The supernatants were collected and measured as described

187

above.

188 189

Numerical model

190

Experiments were simulated with MATLAB (version R2015b) using a simple, parsimonious,

191

process-based kinetic model to describe and test mechanisms of microbial Fe(III) reduction, growth of

192

the Fe-reducing microbial community, abiotic Fe(II) oxidation by O2, and solid phase partitioning of Fe

193

minerals into distinct pools. Redox oscillations are imposed by setting O2 concentrations to 200 μM

194

during oxic periods (slightly below saturation) and 0 during anoxic periods, with Fe(II), Fe- reducing

195

microbial community, and a fast and a slowly reacting Fe(III) solid phase (FeOxfast and FeOxslow,

196

respectively) as state variables. During anoxic periods, the rate of microbial Fe(III) reduction is given

197

as:

198

Rred = ∑ ki [ BM ] i

199

[ FeOxi ] K m, i + [ FeOxi ]

Eq. 1

where ki is a cell specific rate constant, [BM] is the concentration of Fe reducers, and Km(i) is a half 9

ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 29

200

saturation constants and i represents the fast and slowly reacting Fe(III) solid phases, respectively. The

201

net growth rate of Fe-reducing microbes (RBM) is set proportional to the Fe reduction rate, RBM = γ Rred

202

, with a growth efficiency γ set to 0.06 (based on growth yields from Kostka et al.50). As the oxic

203

periods are relatively short in our experiments, no net growth or death of the Fe reducers is considered

204

during oxic periods. Under oxic conditions, Fe reduction is set to 0, and the oxidation of Fe(II) is

205

implemented following Stumm and Morgan51: 𝑅!" = 𝑘!" 𝐹𝑒 !! 𝑂!

206

Eq. 2

207

where the Fe(II) oxidation rate constant kox is set to 3 × 102 M-1 min-1, 1000 fold higher than the value

208

reported for homogeneous oxidation kinetics at the pH of 5.5 in our the experiments to accommodate

209

heterogeneous reaction rates52. Fast oxidation can promote the production of less ordered and more

210

reactive Fe oxides53, 54, thus we partition all the oxidized Fe(II) to the fast pool given the rapid

211

oxidation rate observed in our experiment. During our experiments, we routinely observed that a

212

portion of the HCl-extractable Fe(II) formed during the anoxic half-cycle is resistant to oxidation,

213

resulting in an apparent accumulation of Fe(II) over the course of the experiment. We represented this

214

in the numerical model by removing a small portion of the Fe(II) from the reactive pool, which we

215

hypothesize results from physical protection (encasement) against oxidation. The rate of Fe(II) encased

216

was set to 2% (q) of the rate of iron oxidation to fit the experimental data. Thus, the dynamics of fast

217

and slowly reacting Fe(III) solid phases are given by:

218 219

!!"#$!"#$ !" !!"#$!"#$ !"

!"!"

= 1 − 𝑞 𝑅!" − 𝑅!"#

Eq. 3

!"#$ = −𝑅!"#

Eq. 4

220

Initial reducible Fe(III) solid concentrations were based on the observed values of 10 and 140 mmol L-1

221

for fast and slowly reacting Fe oxides, respectively, with zero Fe(II). The fast pool represents Fe that

222

can be reduced in less than 3-d39; while the slow pool Fe is set to reproduce the Fe(II) plateau levels 10

ACS Paragon Plus Environment

Page 11 of 29

Environmental Science & Technology

223

observed in the oscillation experiments. To match the Fe reduction rates across all treatments, the two

224

model parameters representing kfast *[BM]t=0 and kslow*[BM]t=0 were set to 6 and 1.5 mmol kg-1 d-1 in

225

experiments without P additions, and 10 and 2.3 mmol kg-1 d-1 in experiments with P addition,

226

reflecting alleviation of P limitation. The fast and slow rate parameters were determined by fitting the

227

data from FeIII-reduction rate experiment conducted after 28-days of incubation (see above). All other

228

model parameters were identical across all simulations.

229 230

Results and Discussion

231 232 233

Fe reduction dynamics Fe(II) concentrations responded to changing redox conditions (Figure 1) with peak values at the

234

end of the anoxic periods and the lowest values after an oxic period. Peak Fe(II) concentrations

235

increased from 10 to 148-165 mmol kg-1 soil in the P-amended treatments as the experiment progressed

236

through successive redox cycles (Figure 1). The peak Fe(II) concentrations plateaued near 150 mmol

237

FeII kg-1 soil after about 1 month and increased only slowly for the remainder of the anoxic half-cycles.

238

The phosphate-amended treatments reached this ~150 mmol FeII kg-1 soil plateau at earlier oscillation

239

cycles than the unamended treatments, but by the end of the experiment all treatments exhibited similar

240

amplitudes of Fe(II) oscillation, except the 3.5 d unamended control soil (Figure 1).

241

The increase in peak Fe(II) concentrations was similar across most oscillation frequency

242

treatments, despite different anoxic incubation lengths (Figure SI-1). Thus, the time-averaged iron

243

reduction rate is lower in treatments with shorter oscillation frequencies, but the similar peak Fe(II)

244

concentrations between treatments suggests total Fe reduction extent is independent of oscillation

245

frequency. To assess how the frequency of redox oscillations affected the reduction kinetics, we

246

repeated our oscillation experiments (for the phosphorus addition treatments only) initiating a 11

ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 29

247

continuous 2-week anoxic incubation after one month of redox oscillations (Figure 2). By tracking the

248

accumulation of Fe(II) over time (Figure 2), we surmise there are two phases in the production of

249

Fe(II) during anoxic incubations: (1) an initial, rapid Fe(II) production interval followed by (2) an

250

interval with much slower Fe(II) production that brings the Fe(II) concentrations into proximity of the

251

peak Fe(II) concentrations at the end of each anoxic cycle (Figure 1). Irrespective of oscillation

252

frequency, a slow Fe-reduction phase (or plateau) occurs after 4 days of anoxic incubation. This plateau

253

could explain why all of the phosphate-amended soil slurries had similar peak Fe(II) concentrations

254

(even the 3.5-day treatment), as 90% of the plateau FeII concentration is achieved within 3 days of

255

incubation (Figure 1). Although the second experiment was only conducted with the P-amended soils,

256

we suspect this plateau may not manifest as prominently in the 3.5-d unamended treatment because the

257

overall slower Fe reduction rates in the unamended controls.

258

Our experiments were conducted on re-wetted soils that were oxidized without subsequent

259

drying during the redox oscillations. This experimental design minimized heterogeneity during the

260

experiment and tightened reproducibility between replicates, however it may have impacted the iron

261

dynamics in two important ways. First, the re-wetting of soils is well known to produce a flush of

262

organic carbon55, and thus our treatments likely contain more soluble, labile organic carbon than if we

263

had initiated the experiments from a field moist state. This likely provided an excess of electron donor

264

in the form of soluble organic carbon and it may have increased iron availability through complexation

265

or by facilitating electron shuttling relative to the field moist state56-60. Second, oxidizing the samples

266

under saturated conditions represents the common situation in moist soils or near-surface sediments61

267

in which pulses of O2 accompany the infiltration of rainwater—which has high dissolved O2 levels;

268

whereas is our design does not represent as well soil-drying events, where oxidation fronts follow the

269

major fluid flow paths and create a heterogeneous distribution of FeIII phases62.

270 12

ACS Paragon Plus Environment

Page 13 of 29

271 272

Environmental Science & Technology

Fe mineral composition To quantify the changes in the availability of Fe for microbial reduction63, we inoculated the

273

same freeze-dried soils (before and after treatment) with S. oneidensis MR-1 in a defined media to

274

quantify the rapidly reducible Fe(III) content (Figure 3). The purpose of this experiment was not to

275

simulate Fe-reduction that occurs in the field, but to determine the bioavailability of Fe after exposure

276

to fluctuating redox conditions. The S. oneidensis cultures produced FeII more rapidly following the

277

3.5-day, 7-day, 14-day, and anoxic control slurry incubations than in the oxic control and slurries with

278

soils prior to the oscillation experiments. We found Fe(III) reduction in the non-inoculated soils was

279

negligible compared to that occurring in the inoculated soils. Thus, outside of a possible stimulatory

280

effect of Shewanella on the growth of the native soil bacteria32, we assume the Fe(III)-reduction assays

281

were not biased by differences in soil microbial population or activity.

282

The 3.5-day, 7-day, 14-day, and anoxic control soils exhibited similarly higher rates of Fe(III)

283

reduction in the S. oneidensis assays than the oxic and initial soils (Figure 3, SI section 4). Although

284

each treatment had different redox oscillation periods, all are likely to have similar abundances of

285

recently oxidized Fe(III), as each oscillation treatment reached similar Fe(II) concentration plateaus in

286

their last anoxic cycle, whereas the oxic control and the starting soil had much lower HCl-extractable

287

Fe(II) concentrations. Since we terminated each treatment following an oxidation half cycle (including

288

the anoxic control), the rapid introduction of O2 oxidizes nearly all of the recently produced Fe2+(aq),

289

likely forming a very reactive, short-range-ordered (SRO) phase54.

290

However, we did not find significant changes in Fe mineral composition between treatments

291

using citrate-ascorbate chemical extractions (Figure S2), 57Fe Mössbauer spectroscopy (MBS, SI

292

section 2), or X-ray adsorption spectroscopy (XAS, SI section 3). This suggests one or more of the

293

following: (1) the proportional change in the SRO pool is too small relative to the total background soil

294

Fe in the soil for these chemical and spectroscopic measures; (2) the re-oxidized Fe is more proximal to 13

ACS Paragon Plus Environment

Environmental Science & Technology

295

the microbial iron reducers—perhaps precipitating as thin surface coatings on the original SRO-Fe

296

phases, while retaining similar local structure or by forming complexes with extracellular polymeric

297

substances on the bacterial surface64-66; (3) the re-oxidized Fe co-precipitated with organic electron

298

shuttling agents that can increase the availability for microbial reduction without yielding detectable

299

changes in mineral composition.

300

Page 14 of 29

Prior work from redox fluctuation experiments illustrates that Fe mineral crystallinity can

301

increase or decrease24, 32, 33, 67, or remain undetected68, although the specific conditions driving the

302

trajectory of mineral change are unclear. Documented increases in Fe mineral crystallinity might be

303

associated with slow or incomplete Fe oxidation. Liptzin and Silver67 found increases in Fe

304

crystallinity when exposing intact, un-homogenized soil samples to redox fluctuations and Thompson

305

et al.24 similarly found Fe crystallinity increases when redox fluctuating soil slurries were exposed to

306

progressive oxidation intended to approximate conditions within an aggregate. Likewise, if the initial

307

phase has very low crystallinity, then exposure to Fe2+(aq) during the anoxic periods may promote

308

increases in crystallinity as Tomaszewski et al.26 showed when exposing ferrihydrite to a redox

309

fluctuating environment. However, decreases in Fe crystallinity have also observed during redox

310

fluctuations, especially when accompanied with soil drying69, 70, or when the starting material is of

311

higher crystallinity32. Furthermore, there are other reports, akin to our study, where changes in Fe

312

crystallinity cannot be detected following redox fluctuations68.

313 314

Process identification

315

We identified process dynamics consistent with the observational data by implementing

316

different processes in a numerical model. Our implementation is able to reproduce the observed

317

increase in Fe(II) concentrations (aqueous plus 0.5M HCl extractable) during the anoxic cycles, as well

318

as the plateau phase revealed in our 2nd experiment (Figure 1). The basic structure of the model 14

ACS Paragon Plus Environment

Page 15 of 29

Environmental Science & Technology

319

includes two pools of iron with different cell-specific Fe reduction rates (kfast and kslow). This was

320

critical for reproducing the high Fe(II) production rates at the beginning of each anoxic period,

321

followed by slower Fe(II) production rates associated with reduction of the less reactive iron pool.

322

A key feature of the whole data set (P-amended and non-amended treatments) is that net Fe(II)

323

production rates increase over the course of the experiment (Figure 1). The initial rates of Fe reduction

324

were ~ 3 mmol Fe(II) kg-1 soil d-1, which are similar to rates measured previously on this soil39.

325

However, during the later half of the experiment, the Fe reduction rates increased to 17-66 mmol FeII

326

kg-1 soil d-1 across all treatments (Figure 4). This increase in peak Fe(II) concentrations over the course

327

of the experiment can be explained by the combined effect of an increased pool of reactive iron

328

oxide—formed during rapid oxidation of Fe(II) accumulated during the anoxic periods—and a minor

329

increase in the activity of the iron reducing community. We measured Fe reduction population densities

330

at the beginning and end of the experiment (Table S1) and found no difference within the log unit

331

precision of our methods, which—consistent with the model simulations—implies the Fe reducer

332

population increased less than tenfold over the course of the experiment. Fe(II) concentrations also

333

increase slightly during the oxic half-cycles over the course of the experiment, amounting to what

334

appears as a upward drifting Fe(II) baseline (Figure 1). This may reflect the formation of a mixed-

335

valance iron mineral, a stable Fe(II)-organic matter complex or an Fe sulfide phase, represented in the

336

model by the partitioning a portion of Fe(II) into a non-oxidizable form that is still extractable via 0.5M

337

HCl (see methods). Thus, in our final model, we assumed that the “active Fe”—defined as Fe that

338

could be mobilized via Fe reduction or other processes within the timescale of our experiment—moved

339

between four reservoirs: (1) a rapidly reducible Fe(III) solid; (2) a slowly reducible Fe(III) solid; (3)

340

dissolved or adsorbed Fe(II); and (4) a non-oxidizable Fe(II) phase. Both Fe(III) solid phases

341

contributed to the trend of increasing FeII peaks over the course of the oscillation experiment. During

342

each oxidation cycle, we partitioned all of the Fe(II) oxidized into the rapidly reducing Fe(III) solid 15

ACS Paragon Plus Environment

Environmental Science & Technology

343

phase, hence providing each subsequent simulated reduction cycle with a larger pool of rapidly

344

reducible Fe(III).

345

Page 16 of 29

An alternative explanation for the plateau in Fe(II) concentrations is that sorption of Fe(II) to Fe

346

mineral surfaces may pacify the Fe(III) toward further reductive dissolution. Roden and Zachara2

347

attributed a microbial Fe(III)-reduction plateau they observed in pure culture experiments to this effect.

348

However, in our experiment, peak Fe(II) concentrations clearly increase with subsequent oscillation

349

cycles, suggesting Fe reduction rates across all treatments increased as the experiment progressed

350

(Figure 4).

351

Imposing the above constraints and processes does not yield a model fit to all frequency

352

treatments unless the model parameters are tuned to each treatment. However, by imposing a short

353

suppression of simulated microbial growth at the beginning of each anoxic cycle—representing a

354

predominance of intracellular metabolic adjustments over the formation of new biomass—we can

355

achieve a consistent match with the experimental Fe(II) data. A single global parameter set matches the

356

observations across all oscillation frequencies (Figure 1), except that the P-amended treatments

357

required a higher cell-specific Fe reduction rate than the unamended treatments.

358

Considering first the non-amended treatments, the simulated peak rates of Fe reduction are

359

similar for the 3.5-d and 7-d treatments (peak rates of Fe(II) production ~ 30 mmol kg-1 soil d-1), but

360

substantially slower in the 14-d treatment (peak rates of Fe(II) production ~ 17 mmol kg-1 soil d-1).

361

Because the iron reduction rates are only high during the first ~ 4 days of each anoxic period (Figure

362

4), the 14-d treatment has longer period in which the microbes are largely inactive. With microbial

363

growth linked to metabolic activity, this leads to less growth. Over the course of the experiment, the

364

simulated biomass increased 4.3× in both the 3.5-d and 7-d treatments, compared with an increase of

365

only 2.7× in the 14-d treatment (Figure S3), which in turn leads to lower Fe(II) peak concentrations

16

ACS Paragon Plus Environment

Page 17 of 29

366

Environmental Science & Technology

(Figure 1).

367 368 369

Phosphorus considerations Augmenting the soil slurries with phosphate enhanced the production of Fe(II) during the

370

anoxic cycles (Figure 1), suggesting that microbial activity was stimulated. Consistent with this

371

hypothesis, the amount of phosphorus released following exposure to CHCl3 increased over the course

372

of the experiment (Figure S5), suggesting P accumulated in the microbial biomass. No other changes in

373

the phosphorus distributions (based on the Hedley protocol47) were observed in the treatments, with

374

nearly all of the P residing in the NaOH extractable pool (negligible P was found in the aqueous,

375

NaHCO3, and HCl extracts and thus were not reported). NaOH extractable P is considered to be

376

associated with Fe and Al oxide surfaces in the standard interpretation of the Hedley protocol47.

377

Previous measurements of Luquillo soils also found P to be predominately associated with mineral

378

surfaces38, 67. We found there was substantial potential for P uptake in these soils, but the amount and

379

shape of the phosphate adsorption isotherms (which followed typical Langmuir curves) did not differ

380

between treatments and controls (Figure S4, Table S2).

381

The phosphorus addition treatments had higher simulated rates of Fe reduction than the non-

382

amended controls, with peak rates of Fe(II) production ~66, ~52, and ~31 mmol kg-1 soil d-1 for the 3.5-

383

d, 7-d, and 14-d treatments, respectively. Higher cell-specific Fe reduction rate constants are required

384

to match the temporal evolution of Fe(II) oscillations in the P-amended treatments than are required in

385

the non-amended controls (an increase of 66.7% and 50.3% for fast and slowly reacting iron oxides,

386

respectively). Similar to the non-amended treatments, our simulations increase microbial biomass more

387

in the shorter frequency treatments as a direct consequence of their proportionally longer time with

388

access to rapidly reducible iron (Fe(III)RR) than the longer frequency treatments. Over the course of the

389

experiment, the simulated biomass in the P-amended treatments increased 7.0×, 5.1× and 3.2× in the 17

ACS Paragon Plus Environment

Environmental Science & Technology

390

Page 18 of 29

3.5-d, 7-d, and 14-d treatments, respectively.

391

This may be the first study illustrating that additions of P promote an increase in Fe(II)

392

production rates in native soils or sediments however, there is a large body of work based on synthetic

393

minerals. In these pure-system studies, addition of phosphate often alters the rate of Fe(II) production,

394

although the direction and magnitude of change is inconsistent, depending on P concentrations, Fe

395

mineral composition and other factors71-78. Our P-amended soils had 6.46 mmol kg-1 soil (or 0.58 mmol

396

L-1 suspension) of P added, which is lower than most of the P amendments used in prior synthetic iron

397

mineral studies—typical P amendments are 2 - 20 mmol L-1 in the suspension74, 75, with 0.4 mmol L-1

398

often considered a low P treatment73. Thus, given the low amount of P added in our treatments, we

399

suspect the increase in Fe(II) production rates results from an alleviation a P limitation on microbial

400

growth, which is consistent with P uptake to the biomass (Figure S5).

401 402

Ecosystem Implications

403

Upland ecosystems are subjected to short periods of anoxia that often lead to some ephemeral

404

iron reduction, which can influence C cycling and the fate of many nutrients and chemicals. We have

405

shown that successive periods of redox fluctuations over a month timescale with ample organic C can

406

increase soil Fe reduction rates, regardless of fluctuation frequency. Thus, Fe reduction rates during the

407

onset of anoxia are strongly coupled to physiochemical dynamics of the recent past. Terrestrial systems

408

exposed to rain events that yield iron reducing conditions and a build-up of Fe(II) thus might be primed

409

for more rapid iron reduction rates than systems emerging from a long drought or with continuously

410

wet soils.

411

The amount of Fe(II) produced following a prolonged anoxic incubation is often predicted

412

through selective extractions (oxalate, citrate-ascorbate, hydroxylamine, etc.) targeting the short-range-

413

ordered (SRO) phases. Our work is consistent with this in that the total amount of Fe(II) produced 18

ACS Paragon Plus Environment

Page 19 of 29

Environmental Science & Technology

414

during an anoxic cycle later in the experiment (or in the anoxic control) was similar to the citrate-

415

ascorbate extractable iron (~ 30% of total Fe). However, SRO abundance is a poor predictor of iron

416

reduction rate, because the SRO abundance does not change throughout the experiment, while iron

417

reduction rate changes dramatically. Instead, we propose the availability of SRO phases for Fe

418

reduction is enhanced through rapid oxidation events by atmospheric air (~21% O2), creating a sub-set

419

of rapidly reducible SRO-Fe [Fe(III)RR]. We found Shewanella cultures could reduce Fe(III) much

420

faster from soils exposed to redox fluctuations than from the oxic controls, however we could not

421

detect any corresponding shifts in mineral composition via selective extractions, Mössbauer

422

spectroscopy, or X-ray adsorption spectroscopy. Regardless, by the later half of the experiment nearly

423

the entire SRO Fe pool undergoes reduction within four days.

424

Soils that are frequently exposed to redox fluctuations have a large pool of rapidly reducible

425

iron to transfer electrons to during periods of oxygen depletion—a biogeochemical battery 79 that is

426

recharged by rapid oxidation of previously anoxic pore waters. Our work highlights the importance of

427

identifying the reaction timescales and limiting reactants/catalysts (microbes) in predicting iron

428

reduction rates in soils. Prior work suggests this soil does not have significant C- limitations 39, but we

429

did find the soil responded to P-amendments with more rapid increases in iron reduction rate than

430

unamended soils. The P-amended treatments required a faster cell-specific Fe reduction rate when P

431

was added and suggests iron reducer activity (population growth, CO2 production, etc.) is higher when

432

redox fluctuations occur on timescales similar to the process timescales (length of time to oxidize most

433

of the Fe(II) or the length of time to a Fe(II) concentration pseudo-plateau). Soils that are exposed to

434

near-weekly redox pulsing (due to variable rain events or water-table shifts) may therefore be more

435

dominated by iron biogeochemical dynamics than soils that experience redox shifts only seasonally.

436

19

ACS Paragon Plus Environment

Environmental Science & Technology

437

ASSOCIATED CONTENT

438

Supporting Information

439

Four supporting information sections are included: Section 1 contains two additional tables and five

440

additional figures cited in the main text; Section 2 contains a detailed description of the Mössbauer

441

(MBS) characterization of soil Fe; Section 3 contains a detailed description of the XAS

442

characterization of soil Fe and a comparison of the XAS and MBS findings; Section 4 contains a

443

description of the statistical test of Shewanella sp. assay. This material is available free of charge via

444

the Internet at http://pubs.acs.org.

445 446

AUTHOR INFORMATION

447

Corresponding Author

448

*E-mail: [email protected].

449

Notes

450

The authors declare no competing financial interest.

451 452 453 454

ACKNOWLEDGEMENTS We thank Whendee Silver for providing access to the Bisley site; Michelle Scherer, Drew Latta, and

455

Tim Pasakarnis for Mössbauer measurements; Chris Gorski and Prachi Joshi for preliminary

456

electrochemical analyses of these soils; Qing Ma (APS 5-BM-D) and Ryan Davis (SSRL 4-1) for

457

assistance with XAS analysis; and Nehru Mantripragada for laboratory assistance. This work was

458

funded by USDA - NIFA Soil Processes Program AFRI- NIFA Grant no. 2009-65107-05830 to AT

459

and CM; NSF grants EAR-1331841, EAR-1053470, EAR- 1451508, and DEB-1457761 to AT; and

460

NSF OCE-1559087 to YT. Portions of this research were conducted at the Stanford Synchrotron 20

ACS Paragon Plus Environment

Page 20 of 29

Page 21 of 29

461

Environmental Science & Technology

Radiation Lightsource (SSRL) and the Advanced Photon Source (APS).

462 463

21

ACS Paragon Plus Environment

Environmental Science & Technology

Page 22 of 29

464

References

465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517

(1) Loeppert, R. H.; Inskeep, W. P., Iron. In Methods of soil analysis: Part 3-Chemical methods, Sparks, D. L., Ed. Soil Science Society of America: 1996; pp 639-664. (2) Roden, E. E.; Zachara, J. M., Microbial reduction of crystalline iron(III) oxides: Influenceof oxide surface area and potential for cell growth. Env. Sci. Technol. 1996, 30, 1618-1628. (3) Frenzel, P.; Bosse, U.; Janssen, P. H., Rice roots and methanogenesis in a paddy soil: Ferric iron as an alternative electron acceptor in the rooted soil. Soil Biol. Biochem. 1999, 31, 421-430. (4) Ann, Y.; Reddy, K. R.; Delfino, J. J., Influence of redox potential on phosphorus solubility in chemically amended wetland organic soils. Ecol. Eng. 2000, 14, 169-180. (5) Miller, A. J.; Schuur, E. A. G.; Chadwick, O. A., Redox control of phosphorus pools in hawaiian montane forest soils. Geoderma 2001, 102, 219-237. (6) Fimmen, R.; Richter, D.; Vasudevan, D.; Williams, M.; West, L., Rhizogenic Fe–C redox cycling: A hypothetical biogeochemical mechanism that drives crustal weathering in upland soils. Biogeochemistry 2008, 87, 127. (7) Sulzberger, B.; Suter, D.; Siffert, C.; Banwart, S.; Stumm, W., Dissolution of Fe(III) (hydr)oxides in natural waters; laboratory assessment on the kinetics controlled by surface coordination. Mar. Chem. 1989, 28, 127-144. (8) Bruemmer, G. W.; Gerth, J.; Tiller, K. G., Reaction kinetics of the adsorption and desorption of nickel, zinc and cadmium by goethite. I. Adsorption and diffusion of metals. J. Soil Sci. 1988, 39, 37-52. (9) Sparks, D. L., Kinetics of soil chemical processes. Academic Press: 2013. (10) Spadini, L.; Manceau, A.; Schindler, P. W.; Charlet, L., Structure and stability of Cd2+ surface complexes on ferric oxides: 1. Results from EXAFS spectroscopy. J. Colloid Interf. Sci. 1994, 168, 73-86. (11) Canfield, D. E.; Thamdrup, B.; Hansen, J. W., The anaerobic degradation of organic-matter in danish coastal sediments: Iron reduction, manganese reduction, and sulfate reduction. Geochim. Cosmochim. Acta 1993, 57, 3867-3883. (12) Lovley, D. R., Organic-matter mineralization with the reduction of ferric iron - a review. Geomicrobiol. J. 1987, 5, 375-399. (13) McBride, M. B., Environmental chemistry of soils. Oxford University Press: NY, 1994. (14) Dubinsky, E. A.; Silver, W. L.; Firestone, M. K., Tropical forest soil microbial communities couple iron and carbon biogeochemistry. Ecology 2010, 91, 2604-2612. (15) Yang, W. H.; Liptzin, D., High potential for iron reduction in upland soils. Ecology 2015, 96, 2015-2020. (16) Hall, S. J.; Liptzin, D.; Buss, H. L.; DeAngelis, K.; Silver, W. L., Drivers and patterns of iron redox cycling from surface to bedrock in a deep tropical forest soil: A new conceptual model. Biogeochemistry 2016, 130, 177-190. (17) Bonneville, S.; Behrends, T.; Van Cappellen, P., Solubility and dissimilatory reduction kinetics of iron(III) oxyhydroxides: A linear free energy relationship. Geochim. Cosmochim. Acta 2009, 73, 5273-5282. (18) Bonneville, S.; Van Cappellen, P.; Behrends, T., Microbial reduction of iron(III) oxyhydroxides: Effects of mineral solubility and availability. Chem. Geol. 2004, 212, 255-268. (19) Roden, E. E.; Wetzel, R. G., Organic carbon oxidation and suppression of methane production by microbial Fe(III) oxide reduction in vegetated and unvegetated freshwater wetland sediments. Limnol. Oceanogr. 1996, 41, 1733-1748. (20) Roden, E. E., Geochemical and microbiological controls on dissimilatory iron reduction. C. R. Geosci. 2006, 338, 456-467. (21) Silver, W. L.; Lugo, A. E.; Keller, M., Soil oxygen availability and biogeochemistry along rainfall and topographic gradients in upland wet tropical forest soils. Biogeochemistry 1999, 44, 301-328. (22) Sexstone, A. J.; Revsbech, N. P.; Parkin, T. B.; Tiedje, J. M., Direct measurement of oxygen profiles and denitrification rates in soil aggregates. Soil Sci. Soc. Am. J. 1985, 49, 645-651. (23) Roden, E. E., Analysis of long-term bacterial vs. Chemical Fe(III) oxide reduction kinetics. Geochim. Cosmochim. Acta 2004, 68, 3205-3216. (24) Thompson, A.; Chadwick, O. A.; Rancourt, D. G.; Chorover, J., Iron-oxide crystallinity increases during soil redox oscillations. Geochim. Cosmochim. Acta 2006, 70, 1710-1727. (25) Mejia, J.; Roden, E. E.; Ginder-Vogel, M., Influence of oxygen and nitrate on Fe (hydr) oxide mineral transformation and soil microbial communities during redox cycling. Env. Sci. Technol. 2016, 50, 3580-3588. (26) Tomaszewski, E. J.; Cronk, S. S.; Gorski, C. A.; Ginder-Vogel, M., The role of dissolved Fe (II) concentration in the mineralogical evolution of Fe (hydr) oxides during redox cycling. Chem. Geol. 2016, 438, 163-170. (27) Vogelsang, V.; Kaiser, K.; Wagner, F. E.; Jahn, R.; Fiedler, S., Transformation of clay-sized minerals in soils exposed to prolonged regular alternation of redox conditions. Geoderma 2016, 278, 40-48. (28) Zachara, J. M.; Kukkadapu, R. K.; Fredrickson, J. K.; Gorby, Y. A.; Smith, S. C., Biomineralization of poorly

22

ACS Paragon Plus Environment

Page 23 of 29

518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573

Environmental Science & Technology

crystalline Fe(III) oxides by dissimilatory metal reducing bacteria (DMRB). Geomicrobiol. J. 2002, 19, 179-207. (29) Hansel, C. M.; Benner, S. G.; Neiss, J.; Dohnalkova, A.; Kukkadapu, R. K.; Fendorf, S., Secondary mineralization pathways induced by dissimilatory iron reduction of ferrihydrite under advective flow. Geochim. Cosmochim. Acta 2003, 67, 2977-2992. (30) Chadwick, O. A.; Gavenda, R. T.; Kelly, E. F.; Ziegler, K.; Olson, C. G.; Elliott, W. C.; Hendricks, D. M., The impact of climate on the biogeochemical functioning of volcanic soils. Chem. Geol. 2003, 202, 195-223. (31) Chorover, J.; Amistadi, M. K.; Chadwick, O. A., Surface charge evolution of mineral-organic complexes during pedogenesis in hawaiian basalt. Geochim. Cosmochim. Acta 2004, 68, 4859-4876. (32) Coby, A. J.; Picardal, F.; Shelobolina, E.; Xu, H.; Roden, E. E., Repeated anaerobic microbial redox cycling of iron. Appl. Environ. Microbiol. 2011, 77, 6036-6042. (33) Li, Y.; Yu, S.; Strong, J.; Wang, H., Are the biogeochemical cycles of carbon, nitrogen, sulfur, and phosphorus driven by the “FeIII–FeII redox wheel” in dynamic redox environments? J. Soil Sed. 2012, 12, 683-693. (34) Thompson, A.; Rancourt, D. G.; Chadwick, O. A.; Chorover, J., Iron solid-phase differentiation along a redox gradient in basaltic soils. Geochim. Cosmochim. Acta 2011, 75, 119-133. (35) Dzombak, D. A.; Morel, F. M. M., Surface complexation modeling: Hydrous ferric oxide. John Wiley & Sons: 1990. (36) Axe, L.; Anderson, P. R., Experimental and theoretical diffusivities of Cd and sr in hydrous ferric oxide. J. Colloid Interf. Sci. 1997, 185, 436-448. (37) Trivedi, P.; Axe, L., Modeling Cd and Zn sorption to hydrous metal oxides. Env. Sci. Technol. 2000, 34, 2215-2223. (38) Peretyazhko, T.; Sposito, G., Iron(III) reduction and phosphorous solubilization in humid tropical forest soils. Geochim. Cosmochim. Acta 2005, 69, 3643. (39) Ginn, B. R.; Habteselassie, M. Y.; Meile, C.; Thompson, A., Effects of sample storage on microbial Fe-reduction in tropical rainforest soils. Soil Biol. Biochem. 2014, 68, 44-51. (40) Beinroth, F. H., Some highly weathered soils of Puerto Rico, 1. Morphology, formation and classification. Geoderma 1982, 27, 1-73. (41) Frizano, J.; Johnson, A. H.; Vann, D. R.; Scatena, F. N., Soil phosphorus fractionation during forest development on landslide scars in the luquillo mountains, Puerto Rico. Biotropica 2002, 34, 17-26. (42) White, A. F.; Blum, A. E.; Schulz, M. S.; Vivit, D. V.; Stonestrom, D. A.; Larsen, M.; Murphy, S. F.; Eberl, D., Chemical weathering in a tropical watershed, luquillo mountains, Puerto Rico: I. Long-term versus short-term weathering fluxes. Geochim. Cosmochim. Acta 1998, 62, 209-226. (43) Tishchenko, V.; Meile, C.; Scherer, M. M.; Pasakarnis, T. S.; Thompson, A., Fe2+ catalyzed iron atom exchange and re-crystallization in a tropical soil. Geochim. Cosmochim. Acta 2015, 148, 191-202. (44) Achtnich, C.; Bak, F.; Conrad, R., Competition for electron donors among nitrate reducers, ferric iron reducers, sulfate reducers, and methanogens in anoxic paddy soil. Biol. Fertil. Soils 1995, 19, 65-72. (45) Thompson, A.; Chadwick, O. A.; Boman, S.; Chorover, J., Colloid mobilization during soil iron redox oscillations. Env. Sci. Technol. 2006, 40, 5743-5749. (46) Reyes, I.; Torrent, J., Citrate-ascorbate as a highly selective extractant for poorly crystalline iron oxides. Soil Sci. Soc. Am. J. 1997, 61, 1647-1654. (47) Hedley, M. J.; Stewart, J. W. B.; Chauhan, B., Changes in inorganic and organic soil phosphorus fractions induced by cultivation practices and by laboratory incubations. Soil Sci. Soc. Am. J. 1982, 46, 970-976. (48) Hedley, M. J.; Stewart, J. W. B., Method to measure microbial phosphate in soils. Soil Biol. Biochem. 1982, 14, 377385. (49) Murphy, J.; Riley, J. P., A modified single solution method for the determination of phosphate in natural waters. Analytica chimica acta 1962, 27, 31-36. (50) Kostka, J. E.; Dalton, D. D.; Skelton, H.; Dollhopf, S.; Stucki, J. W., Growth of iron(III)-reducing bacteria on clay minerals as the sole electron acceptor and comparison of growth yields on a variety of oxidized iron forms. Appl. Environ. Microbiol. 2002, 68, 6256-6262. (51) Stumm, W.; Morgan, J. J., Aquatic chemistry: Chemical equilibrium and rates in natural waters. 3rd ed.; John Wiley & Sons, Inc.: New York, NY, 1996; p 1022. (52) Van Beek, C.; Hiemstra, T.; Hofs, B.; Nederlof, M.; Van Paassen, J.; Reijnen, G., Homogeneous, heterogeneous and biological oxidation of iron (II) in rapid sand filtration. J. Water Supply Res. T. 2012, 61, 1-13. (53) Senko, J. M.; Dewers, T. A.; Krumholz, L. R., Effect of oxidation rate and Fe (II) state on microbial nitratedependent Fe (III) mineral formation. Appl. Environ. Microbiol. 2005, 71, 7172-7177. (54) Steefel, C. I.; Van Cappellen, P., A new kinetic approach to modeling water-rock interaction: The role of nucleation, precursors, and ostwald ripening. Geochim. Cosmochim. Acta 1990, 54, 2657-2677. (55) Lundquist, E. J.; Jackson, L. E.; Scow, K. M., Wet–dry cycles affect dissolved organic carbon in two california agricultural soils. Soil Biol. Biochem. 1999, 31, 1031-1038.

23

ACS Paragon Plus Environment

Environmental Science & Technology

574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629

Page 24 of 29

(56) Kleber, M.; Eusterhues, K.; Keiluweit, M.; Mikutta, C.; Mikutta, R.; Nico, P. S., Mineral–organic associations: Formation, properties, and relevance in soil environments. In Advances in agronomy, Donald, L. S., Ed. Academic Press: 2015; Vol. 130, pp 1-140. (57) Kappler, A.; Newman, D. K., Formation of Fe(III)-minerals by Fe(II)-oxidizing photoautotrophic bacteria. Geochim. Cosmochim. Acta 2004, 68, 1217-1226. (58) Taillefert, M.; Beckler, J. S.; Carey, E.; Burns, J. L.; Fennessey, C. M.; DiChristina, T. J., Shewanella putrefaciens produces an Fe(III)-solubilizing organic ligand during anaerobic respiration on insoluble Fe(III) oxides. J. Inorg. Biochem. 2007, 101, 1760-1767. (59) Royer, R. A.; Burgos, W. D.; Fisher, A. S.; Unz, R. F.; Dempsey, B. A., Enhancement of biological reduction of hematite by electron shuttling and Fe(II) complexation. Env. Sci. Technol. 2002, 36, 1939-1946. (60) Royer, R. A.; Burgos, W. D.; Fisher, A. S.; Jeon, B.-H.; Unz, R. F.; Dempsey, B. A., Enhancement of hematite bioreduction by natural organic matter. Env. Sci. Technol. 2002, 36, 2897-2904. (61) Couture, R.-M.; Charlet, L.; Markelova, E.; Madé, B. t.; Parsons, C. T., On–off mobilization of contaminants in soils during redox oscillations. Env. Sci. Technol. 2015, 49, 3015-3023. (62) Schulz, M.; Stonestrom, D.; Lawrence, C.; Bullen, T.; Fitzpatrick, J.; Kyker-Snowman, E.; Manning, J.; Mnich, M., Structured heterogeneity in a marine terrace chronosequence: Upland mottling. Vadose Zone J. 2016, 15. (63) Eusterhues, K.; Neidhardt, J.; Hädrich, A.; Küsel, K.; Totsche, K., Biodegradation of ferrihydrite-associated organic matter. Biogeochemistry 2014, 119, 45-50. (64) Warren, L. A.; Ferris, F. G., Continuum between sorption and precipitation of Fe (III) on microbial surfaces. Env. Sci. Technol. 1998, 32, 2331-2337. (65) Frankel, R. B.; Bazylinski, D. A., Biologically induced mineralization by bacteria. Rev. Mineral. Geochem. 2003, 54, 95-114. (66) Mikutta, C.; Mikutta, R.; Bonneville, S.; Wagner, F.; Voegelin, A.; Christl, I.; Kretzschmar, R., Synthetic coprecipitates of exopolysaccharides and ferrihydrite. Part I: Characterization. Geochim. Cosmochim. Acta 2008, 72, 11111127. (67) Liptzin, D.; Silver, W. L., Effects of carbon additions on iron reduction and phosphorus availability in a humid tropical forest soil. Soil Biol. Biochem. 2009, 41, 1696-1702. (68) Parsons, C. T.; Couture, R. M.; Omoregie, E. O.; Bardelli, F.; Greneche, J. M.; Roman-Ross, G.; Charlet, L., The impact of oscillating redox conditions: Arsenic immobilisation in contaminated calcareous floodplain soils. Environ. Pollut. 2013, 178, 254-263. (69) McGeehan, S. L.; Fendorf, S. E.; Naylor, D. V., Alteration of arsenic sorption in flooded-dried soils. Soil Sci. Soc. Am. J. 1998, 62, 828-833. (70) Mansfeldt, T.; Schuth, S.; Häusler, W.; Wagner, F.; Kaufhold, S.; Overesch, M., Iron oxide mineralogy and stable iron isotope composition in a gleysol with petrogleyic properties. J. Soil Sed. 2012, 12, 97-114. (71) Fredrickson, J. K.; Zachara, J. M.; Kennedy, D. W.; Dong, H. L.; Onstott, T. C.; Hinman, N. W.; Li, S. M., Biogenic iron mineralization accompanying the dissimilatoryreduction of hydrous ferric oxide by a groundwater bacterium. Geochim. Cosmochim. Acta 1998, 62, 3239-3257. (72) Zachara, J. M.; Fredrickson, J. K.; Li, S. M.; Kennedy, D. W.; Smith, S. C.; Gassman, P. L., Bacterial reduction of crystalline Fe3+ oxides in single phasesuspensions and subsurface materials. Am. Mineral. 1998, 83, 1426-1443. (73) Glasauer, S.; Weidler, P. G.; Langley, S.; Beveridge, T. J., Controls on Fe reduction and mineral formation by a subsurface bacterium. Geochim. Cosmochim. Acta 2003, 67, 1277-1288. (74) Kukkadapu, R. K.; Zachara, J. M.; Fredrickson, J. K.; Kennedy, D. W., Biotransformation of two-line silicaferrihydrite by a dissimilatory Fe(III)-reducing bacterium: Formation of carbonate green rust in the presence of phosphate. Geochim. Cosmochim. Acta 2004, 68, 2799-2814. (75) Borch, T.; Masue, Y.; Kukkadapu, R. K.; Fendorf, S., Phosphate imposed limitations on biological reduction and alteration of ferrihydrite. Env. Sci. Technol. 2007, 41, 166-172. (76) O’Loughlin, E. J.; Boyanov, M. I.; Flynn, T. M.; Gorski, C. A.; Hofmann, S. M.; McCormick, M. L.; Scherer, M. M.; Kemner, K. M., Effects of bound phosphate on the bioreduction of lepidocrocite (γ-FeOOH) and maghemite (γ-Fe2O3) and formation of secondary minerals. Env. Sci. Technol. 2013, 47, 9157-9166. (77) O’Loughlin, E. J.; Gorski, C. A.; Scherer, M. M.; Boyanov, M. I.; Kemner, K. M., Effects of oxyanions, natural organic matter, and bacterial cell numbers on the bioreduction of lepidocrocite (γ-FeOOH) and the formation of secondary mineralization products. Env. Sci. Technol. 2010, 44, 4570-4576. (78) Zachara, J. M.; Kukkadapu, R. K.; Peretyazhko, T.; Bowden, M.; Wang, C. M.; Kennedy, D. W.; Moore, D.; Arey, B., The mineralogic transformation of ferrihydrite induced by heterogeneous reaction with bioreduced anthraquinone disulfonate (AQDS) and the role of phosphate. Geochim. Cosmochim. Acta 2011, 75, 6330-6349. (79) Byrne, J. M.; Klueglein, N.; Pearce, C.; Rosso, K. M.; Appel, E.; Kappler, A., Redox cycling of Fe (II) and Fe (III) in

24

ACS Paragon Plus Environment

Page 25 of 29

630 631

Environmental Science & Technology

magnetite by Fe-metabolizing bacteria. Science 2015, 347, 1473-1476.

632

25

ACS Paragon Plus Environment

Environmental Science & Technology

633

Page 26 of 29

Figure 1

634

Unamended

635 636 637 638 639 640 641 642 643 644

40

0

14

7d

160

28

Days

120 80 40

200

0

7

Days

120 80 40 0

0

7

200

Phosphate amended 14 d

160

Figure 1: Fluctuations and model

120 80 40 200 56

0

14

7d

160

28

42

56

Days

120 80 40

200 14 21 28 35 42 49 56

3.5 d

160

42

2+

200

2+

80

0.5 M HCl extractable Fe -1 mmol kg

120

0.5 M HCl extractable Fe -1 mmol kg

2+

14 d

160

0.5 M HCl extractable Fe -1 mmol kg

0.5 M HCl extractable Fe

0.5 M HCl extractable Fe -1 mmol kg

2+

2+

(mmol kg )

0.5 M HCl extractable Fe 0.5 M HCl extractable Fe 2+ -1 -1 -1 mmol kg mmol kg

2+

200

0

160

7

14 21 28 35 42 49 56

3.5 d

Days

120 80 40

0 14 21 28 35 42 49 56

0

7

Days

14 21 28 35 42 49 56

Days

Figure 1. Fluctuations in FeII concentrations due to redox oscillations with periods of 3.5 d (bottom), 7 d (middle) or 14 d (top). All treatments have the same oxic:anoxic exposure time ratio of 6:1 and the oxic cycles involve exposure to 21% O2. Graphs on the left side are for un-amended slurries and those on the right are for treatments with added phosphorous. Symbols indicate measured values with 1 s.d. error bars based on three replicates. The solid lines represent the numerical model.

26

ACS Paragon Plus Environment

Page 27 of 29

645 646

Environmental Science & Technology

Figure 2

647

648 649 650 651 652 653 654

Figure 2. FeII-concentrations under anoxic conditions following 28 d of redox oscillations with periods of 3.5 d, 7 d, or 14 d in treatments amended with phosphorous. Error bars are 1 s.d. This data was used to fit the dynamics during the reduction period in the numerical model across all treatments.

27

ACS Paragon Plus Environment

Environmental Science & Technology

655

Page 28 of 29

Figure 3

Oscillations & anoxic control after oxidation

Unreacted & oxic control after oxidation

656 657 658 659 660 661 662 663 664 665

Figure 3 Figure 3. FeII concentrations in a liquid culture of S. oneidensis MR-1 incubated with the sole source of ferric iron deriving from freeze-dried soil samples from the indicated redox oscillation treatments and controls (all of which ended with oxic exposure, including the anoxic control). Each data point is the difference in acid-extractable FeII concentrations between inoculated soils and the non-inoculated controls. Error bars are 1 s.d. See SI-Section 4 for more details.

28

ACS Paragon Plus Environment

Page 29 of 29

666

Environmental Science & Technology

Figure 4

667

Phosphate Amended 80

Net Fe(II) Production Rate -1 mmol kg

40

14 d

30 20 10

Net Fe(II) Production Rate -1 -1 mmol kg d

668 669 670 671 672 673 674 675 676

Net Fe(II) Production Rate -1 -1 mmol kg d

Net Fe(II) Production Rate -1 -1 mmol kg d

40

7d

30 20 10 40

3.5 d

30 20 10 0

0

7

14 21 28 35 42 49 56

Net Fe(II) production rate -1 -1 mmol kg d

Net Fe(II) Production Rate -1 -1 mmol kg d

Unamended

14 d

60 40 20 80

7d

60 40 20

80

3.5 d

60 40 20 0

0

7

14 21 28 35 42 49 56

Days

Days

Figure 4 Figure 4. Net Fe(II) production rates calculated at each time point from the numerical model for the unamended (left graphs) and P-amended (right graphs) treatments. Values shown only for anoxic periods.

29

ACS Paragon Plus Environment