Adsorption and molecular fractionation of dissolved organic matter on

1. Adsorption and molecular fractionation of dissolved organic matter on. 1 iron-bearing mineral matrices of varying crystallinity. 2. 3. Elizabeth K...
0 downloads 0 Views 1MB Size
Subscriber access provided by READING UNIV

Article

Adsorption and molecular fractionation of dissolved organic matter on iron-bearing mineral matrices of varying crystallinity Elizabeth K Coward, Tsutomu Ohno, and Alain F. Plante Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b04953 • Publication Date (Web): 12 Jan 2018 Downloaded from http://pubs.acs.org on January 12, 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 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 33

Environmental Science & Technology

1

Adsorption and molecular fractionation of dissolved organic matter on

2

iron-bearing mineral matrices of varying crystallinity

3

4

Elizabeth K. Coward1,2*, Tsutomu Ohno3, Alain F. Plante2

5

1

Delaware Environmental Institute, University of Delaware, Newark DE 19716-7310, USA.

6

2

Department of Earth & Environmental Science, University of Pennsylvania, Philadelphia PA

7

19104-6316, USA.

8

3

School of Food & Agriculture, University of Maine, Orono ME 04469-5763, USA.

ACS Paragon Plus Environment

1

Environmental Science & Technology

Page 2 of 33

9

Abstract. Iron (Fe)-bearing mineral phases contribute disproportionately to adsorption of soil

10

organic matter (SOM), due to their elevated chemical reactivity and specific surface area (SSA).

11

However, the spectrum of Fe solid-phase speciation present in oxidation-reduction-active soils

12

challenges analysis of SOM-mineral interactions and may induce differential molecular

13

fractionation of dissolved organic matter (DOM). This work used paired selective dissolution

14

experiments and batch sorption of post-extraction residues to 1) quantify the contributions of Fe-

15

bearing minerals of varying crystallinity to DOM sorption, and 2) characterize molecular

16

fractionation using Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR-MS).

17

A substantial proportion of soil SSA was derived from extracted Fe-bearing phases, and FT-ICR-

18

MS analysis of extracted DOM revealed distinct chemical signatures across Fe-OM associations.

19

Sorbed carbon (C) was highly correlated with Fe concentrations, suggesting that Fe-bearing

20

phases are strong drivers of sorption in these soils. Molecular fractionation was observed across

21

treatments, particularly those dominated by short-range-order (SRO) mineral phases, which

22

preferentially adsorbed aromatic and lignin-like formulae, and higher-crystallinity phases,

23

associated with aliphatic DOM. These findings suggest Fe speciation-mediated complexation

24

acts as a physicochemical filter of DOM moving through the critical zone, an important

25

observation as predicted changes in precipitation may dynamically alter Fe crystallinity and C

26

stability.

27 ACS Paragon Plus Environment

2

Page 3 of 33

28

Environmental Science & Technology

INTRODUCTION

29

Humid tropical forest soils contain a substantial portion (~500 Pg) of the global terrestrial

30

carbon (C) pool1, yet the persistence of soil organic matter (SOM) in these systems remains

31

unclear, largely due to mechanistic uncertainty in the biogeochemical controls on soil C.

32

Organomineral complexation, driven by high concentrations of iron (Fe)-bearing mineral phases

33

accumulating in surface horizons, has been framed as a potential C stabilization mechanism2-4.

34

However, humid tropical soils often contain a diverse array of Fe-bearing minerals produced by

35

pedogenesis5,

36

confounding our understanding of Fe contributions to C accumulation and stabilization.

37

Particularly, the solid-phase speciation of Fe minerals in soils, such as their size, structure, and

38

surface area, can regulate the dynamics of C sorption reactions8-11, which have been linked to

39

long-term retention of C in soils12,

40

volcanic soils underlying many humid tropical forests14, are exceptionally reactive due to their

41

high specific surface area (SSA) and large numbers of hydroxyl groups present on mineral

42

surfaces15. As such, even small quantities of these minerals may provide an important sorbent for

43

organic compounds.

6

and frequent oxidation-reduction (redox) oscillations in surface horizons7,

13

. Short-range-order (SRO) minerals, abundant in the

44

Tropical systems are also characterized by large dissolved organic matter (DOM) input to

45

surface soils, a highly mobile and reactive pool of SOM 16. In many tropical forests, fine litterfall

46

accounts for ~60% of aboveground net primary productivity17 and rainfall is both frequent and

47

abundant. As such, the movement of litter-leached DOM represents an important, near-constant

48

input of C that fuels substantial heterotrophic respiration and could account for a large proportion

49

of soil CO2 fluxes18. This DOM is thought to be comprised of both largely plant-derived

50

compounds, such as lignin-derived phenols, mobilized through leaching, root tissue and

ACS Paragon Plus Environment

3

Environmental Science & Technology

Page 4 of 33

51

exudates19, 20, and char21. Recent evidence has shown DOM is also highly processed by microbial

52

assimilation and turnover21-23.

53

This abundance of DOM entering humid tropical soils may undergo molecular

54

fractionation at the mineral interface, driven by the molecular diversity of DOM and the

55

spectrum of reactive minerals present. DOM is comprised of numerous small molecular

56

components bound by hydrophobic interactions and hydrogen bonds, which may be broken when

57

external forces from mineral surfaces are greater than the intermolecular forces between the

58

molecular components of DOM24. Molecular components with greater affinity to mineral

59

surfaces, often hydrophobic-, aromatic-, and carboxylic-rich fractions, will likely be bound while

60

DOM characterize by lower affinity functional groups will remain in solution, leading to the

61

molecular fractionation of DOM25, 26. The study of such molecular fractionation has, however,

62

been historically limited to spectroscopic methods, which provide moiety-scale assessment of

63

DOM chemistry, but are not capable of resolving molecular formulae. The chemical

64

characterization of DOM composition has recently been enhanced by the ultrahigh resolution

65

capabilities of electrospray ionization Fourier transform ion cyclotron resonance mass

66

spectrometry (FT-ICR-MS), capable of resolving thousands of individual OM components27.

67

This work uses FT-ICR-MS analysis to investigate the role of Fe-bearing mineral

68

crystallinity in adsorbing and fractionating DOM in humid tropical soil matrices. Such chemical

69

fractionation of DOM by the physicochemical filter of Fe mineralogy may modulate the

70

composition of DOM reaching subsoils and the critical zone at large. Thus, characterizing the

71

contribution of mineral matrices with varying Fe-bearing solid-phase speciation may aid in

72

understanding the fate of DOM at the watershed or ecosystem level, in reconciling the wide

ACS Paragon Plus Environment

4

Page 5 of 33

Environmental Science & Technology

73

range of C turnover times observed in tropical systems, and in disentangling the complex

74

response of a substantial, but understudied, C reservoir to projected changes in climate.

75

We aim to assess the contribution of Fe content and crystallinity to these organomineral

76

complexes via two methodologies: 1) molecular characterization of DOM extracted with various

77

Fe-bearing mineral phases, and 2) assessment of the contribution of Fe-bearing mineral phases to

78

adsorptive molecular fractionation of DOM. Several recent studies have used FT-ICR-MS

79

analysis to provide strong evidence of adsorptive fractionation of organic acids and DOM by

80

synthetic metal oxides 24, 28, 29 and a temperate soil matrix 30. However, the role of Fe solid-phase

81

speciation in governing molecular adsorptive fractionation of complex, litter-derived DOM

82

remains uncertain Inorganic selective dissolution of distinct Fe-bearing mineral phases from

83

humid tropical Oxisol soils and associated OM were followed by FT-ICR-MS analysis of the

84

extracts to detect in-situ molecular fractionation of OM. Fe-induced fractionation was also

85

investigated in a laboratory batch sorption experiment. Soil residues post-dissolution, a

86

Ferrihydrite-coated (Fh+) soil and untreated control soil, comprising five distinct matrices

87

varying in iron content, crystallinity and SSA, were then reacted with DOM in a batch sorption

88

experiment. Unsorbed and post-sorption DOM solutions were analyzed for C concentrations to

89

assess bulk sorption effects and via FT-ICR-MS to determine preferential adsorption of OM

90

components.

91

MATERIALS AND METHODS

92

Soil Sampling and Selective Dissolution.

93

Soil samples originated from an extensive sampling campaign completed in 2010 31 spanning

94

the tropical montane forests of the Luquillo Critical Zone Observatory (LCZO). For the present

95

study, one surface (0 - 20 cm) horizon from a quantitative soil pit of Oxisol soil order was

ACS Paragon Plus Environment

5

Environmental Science & Technology

Page 6 of 33

96

selected to assess the contribution of matrix Fe-bearing mineralogy to OM sorption and

97

molecular fractionation. Detailed site information and soil characteristics may be found in the SI.

98

Six soil subsamples from the source pit were subjected to three parallel inorganic chemical

99

extractions (n=18) for quantification of extractable Fe-containing secondary weathering phases

100

and associated C, as previously detailed in Coward et al.3.

101

Preparation of Sorbent Matrices and Dissolved Organic Matter.

102

The iron oxyhydroxide Ferrihydrite (Fh) was synthesized according to methods described 32

103

by Herbel and Fendorf

104

NaOH until stable at pH 7.5. The suspension was allowed to equilibrate for at least 2 h, during

105

which time the pH was adjusted to maintain a value of 7.5. Synthesized Fe oxyhydroxides were

106

repeatedly washed (4×) by centrifugation, re-suspended in Milli-Q water and immediately added

107

to substrates at 25 °C: one of 30 g soil and another of 30 g 50-mesh quartz sand. Each slurry was

108

shaken horizontally for 48 hours until dry; this process was repeated three times to maximize

109

surficial coating. Ferrihydrite coatings were isolated from quartz sand samples via suspended

110

sonication and centrifugation for mineralogical confirmation. Mineral identity was established by

111

powder X-ray diffraction (pXRD) methods (Fig. S1).

. Briefly, a ferric chloride solution was rapidly titrated with 1.0 N

112

DOM was extracted from O horizon material collected during the 2010 sampling

113

campaign. The composition of the original litter was 31.0 % C and 1.1 % N31. The extract was

114

obtained by shaking a 1:10 soil:Milli-Q water solution for 24 hours. The supernatant was

115

obtained by centrifugation at 5300 × g and vacuum filtration through 0.22 µm nylon filters. The

116

DOM solution was kept in the dark at 4 °C until TOC analysis via a persulfate digestion and

117

spectrophotometric method (Method 10173, Hach Company, Loveland CO). Initial DOC

118

concentration was 200 ± 1.48 mg C L-1, which is relatively high but within the range of

119

concentrations used in previous sorption studies (e.g., Feng et al.33). The high DOC

ACS Paragon Plus Environment

6

Page 7 of 33

Environmental Science & Technology

120

concentration was used in the aim of saturating the system and reaching maximum sorption

121

potential. pH of the final extract was 4.4 due to the acidic character of the water-soluble

122

compounds in the initial litter.

123

Sorption Experiments.

124

Batch sorption experiments were conducted by mixing 0.5 g aliquots of soil with 15 mL of

125

DOM in a 50 mL Erlenmeyer flask. Batch sorption experiments included four replicates of each

126

sorbent matrix and four control samples consisting solely of DOM solution (~200 mg C L-1) to

127

determine the extent of microbial DOM degradation during the experiment. Thus, a total of 24

128

Erlenmeyer flasks were shaken for 48 h on a reciprocal shaker (120 rpm). A biocide was not

129

used and experiments were performed at room temperature (25 °C). After sorption, suspensions

130

were centrifuged at 5,300 × g, filtered through 0.22 µm Whatman filters and stored in the dark at

131

4 °C. Supernatants from the batch sorption experiments were analyzed for DOC as described

132

above. The concentration of DOM adsorbed to the mineral matrices was calculated by the

133

difference between concentrations in the initial and post-sorption solutions. No detectable

134

differences in DOM control solutions were detected, and thus no correction was used in this

135

calculation.

136

FT-ICR-MS Analysis.

137

FT-ICR-MS analysis was used to characterize DOM of (i) selective dissolution extracts,

138

(ii) initial DOM, and (iii) post-sorption supernatants. Solid phase extraction (SPE) was

139

performed on all the samples to desalinate and concentrate the DOC solution according to

140

Dittmar et al.34. Briefly, extracts were prepared by diluting the extracts to 100 mg/L DOC with

141

mass spectrometry grade water and processing through Agilent PPL™ solid-phase extraction

142

cartridges. The adsorbed organic matter was eluted with mass spectrometry grade methanol for

ACS Paragon Plus Environment

7

Environmental Science & Technology

Page 8 of 33

143

FT-ICR-MS analysis, which was conducted with a 12-T Bruker Daltonics Apex Qe FT-ICR-MS

144

housed at the College of Sciences Major Instrumentation Cluster at Old Dominion University.

145

The solid-phase extraction processed samples (eluted in methanol) were diluted by a factor of 5

146

with LCMS grade methanol and water to give a final sample composition of 50:50 (v/v)

147

methanol:water. Further details on SPE and FT-ICR-MS parameters are included in the

148

Supplementary Information.

149

RESULTS & DISCUSSION

150

Selective Dissolution of Fe-associated OM.

151

Chemical properties of soil samples and selective extractants post-dissolution (Table 1)

152

were consistent with previous applications of this selective dissolution method3. Briefly,

153

inorganic dithionite-HCl (DH) extracted crystalline and SRO Fe oxides, hydroxides, and

154

oxyhydroxides derived largely from pedogenesis. HCl-hydroxylamine (HH) extracted SRO Fe-

155

bearing mineral phases, and sodium pyrophosphate (PP) extracted chelated and dispersable Fe

156

phases. Hydroxylamine- and dithionite-extracted Fe were found to be the dominant extractable

157

Fe phases, suggesting crystalline and SRO Fe phases are prevalent in these soils. The highest

158

concentrations of C were solubilized by pyrophosphate extraction, however this extractant

159

solubilized the lowest observed concentrations of Fe. A maximum molar ratio of adsorption C to

160

Fe oxides is 1.0 for the co-extracted species35-37. Co-precipitation, chelation or aggregation of

161

organic ligands generates low- density, organic-rich structures with higher (> 1) C:Fe values Fe-

162

bearing mineral phases, largely via available surface area.36,

163

indicate that hydroxylamine (C:Fe = 0.88 ± 0.054) and dithionite (C:Fe = 0.54 ± 0.038) extracts

164

likely consist of Fe-C bound via inner-sphere adsorptive mechanisms, while pyrophosphate-

38

Resultant molar C:Fe ratios

ACS Paragon Plus Environment

8

Page 9 of 33

Environmental Science & Technology

165

extractable C (C:Fe = 7.44 ± 0.32) may be associated via OM-rich mechanisms, such as co-

166

precipitation or aggregation.

167

Ultrahigh resolution mass spectrometry of the extracts was used to probe differences in

168

OM chemical composition liberated by each extract (Fig. 1). After post-spectral processing to

169

identify peaks present in at least two out of three analytical replicates, > 3,800 mass spectral

170

peaks were present in all three extractants, and elemental compositions were assigned in the

171

range 300 < m/z < 800. van Krevelen diagrams and relative distributions (Fig. 1A-C) of chemical

172

classifications reveal substantial differences in classification class abundances among extracts,

173

while highlighting the predominance of lignin-like compounds across all extracts. The strong

174

adsorption of lignin-like compounds to mineral surfaces has been well documented, and is likely

175

a function of both ligand exchange reactions through prevalent carboxyl groups20,

176

efficiency of microbial degradation of lignin and lignin-like compounds40,

177

persistence of such compounds in DOM.

41

39

and low

, leading to the

178

Beyond the abundance of lignin-like molecular formulae in all three extractant types, we

179

observed strong chemical fractionation among the three extractants, indicating OM associated

180

with Fe phases of varying crystallinity differs at the molecular level. Most notably,

181

pyrophosphate-solubilized DOM (Fig. 1A) differed substantially from DOM solubilized in other

182

extracts (Figs. 1B, 1C, S2). The relative abundances of compounds in each group reveal

183

pyrophosphate-extracted OM contained greater unsaturated aliphatic and saturated carbohydrate-

184

like contributions and fewer condensed lignin and aromatic contributions. DOM molecular

185

formulae in the aliphatic region (H/C > 1.5) of van Krevelen space are likely to be microbially-

186

derived, but may also may be sourced from root necromass or other root inputs to soils, such as

187

root exudates42. Carbohydrates in DOM are derived from decaying microorganisms, and are

ACS Paragon Plus Environment

9

Environmental Science & Technology

Page 10 of 33

188

mainly composed of large aggregates of carbohydrates such as polysaccharides as well as

189

proteins and peptides43. The predominance of carbohydrate-based compounds associated to

190

chelated or colloidal Fe is in accordance to previous findings43, 44 that these inherently polar and

191

highly soluble compounds may not be adsorbed directly to mineral surfaces. More broadly,

192

chelated or colloidal OM appears to be largely microbial in origin, and, due to the weaker

193

electrostatic interactions stabilizing these structures, may represent a more rapidly cycling pool

194

of C. In accordance with the “cytoplasmic solute” hypothesis proposed by Fierer and Schimel45,

195

such hydrophilic compounds may be easily dissolved, and subsequently mineralized, during the

196

frequent precipitation (i.e., Fe-reducing) events in humid forest soils.

197

Hydroxylamine-extracted (Fig. 1B) and dithionite-extracted (Fig. 1C) OM occupied less

198

aliphatic space (H/C > 1.5), and were instead comprised of formulae with higher O/C and higher

199

double-bond equivalence (DBE), a measure of the number of double bonds and rings in a

200

molecule (Fig. S2, modified from analysis presented by Riedel et al., 201246), comprised of

201

polycyclic aromatic, aromatic and lignin components. Preferential adsorption of aromatic

202

structures to hydrous Fe phases is well-established25,

203

solubilized by these two extractants were expected, due to the likely extraction of SRO Fe phases

204

by both, differences between the Hydroxylamine-extracted (Fig. 1B) and dithionite-extracted

205

(Fig. 1C) OM were observed. Hydroxylamine-solubilized OM contained relatively lower

206

contributions of lignin-like components than those observed in dithionite-solubilized OM,

207

accompanied by greater polycyclic aromatic contributions. The increase in dithionite-extracted

208

lignin-like components is likely due to overall increased mineral dissolution and subsequent

209

breakage of ligand-exchange bonds. We suggest that polycyclic aromatic compounds may be

210

preferentially bound to those SRO mineral phases solubilized by hydroxylamine extraction.

47

. While similarities in the formulae

ACS Paragon Plus Environment

10

Page 11 of 33

211

Environmental Science & Technology

DOM sorption by mineral matrices.

212

Post-extraction solid-phases residues, Fh+ soil, and untreated soils, were exposed to O

213

horizon-derived DOM in a batch sorption experiment to assess the contribution of Fe mineral

214

crystallinity and reactivity to matrix adsorption capacity. All sorbent matrices were characterized

215

for total C and Fe content and SSA prior to sorption (Table 1, Figure 2). The mineral matrix

216

treatments represent a wide range of Fe concentrations, which are strongly correlated (R2 = 0.90)

217

with measured SSA (m2 g-1) (Fig. 2A). The strength of the observed correlation confirms the

218

hypothesis that Fe-bearing mineral phases are a critical driver of available surface area in these

219

highly-weathered tropical soils.

220

Measurements of total C and SSA were obtained post-sorption to determine the

221

concentration of C sorbed and C loadings pre- and post-sorption. Notably, sorbed C

222

concentrations, derived from initial and post-sorption DOC measurements, were well-correlated

223

(R2 = 0.76) with Fe content (Fig. 2B). Such robust correlations suggest that Fe-bearing mineral

224

phases, largely via available surface area, are strong controls on C accumulation in these soils. It

225

is important to note untreated soils adsorbed anomalously low C concentrations, which we

226

suggest may be due to desorption or exchange of native SOM or to the existing occlusion of Fe-

227

bearing minerals by OM, as the other treatments either removed Fe but also C (selective

228

dissolution), or added Fe phases without increasing OM content (Fh coating). Comparison of C

229

loadings (mg C m-2) pre- and post-sorption (Fig. 2C) revealed an expected increase in loadings

230

due to sorption across all treatments, consistently correlated with Fe content (R2 = 0.79, R2 =

231

0.80).

232

Molecular fractionation of matrix-sorbed OM.

ACS Paragon Plus Environment

11

Environmental Science & Technology

Page 12 of 33

233

FT-ICR-MS characterization of pre- and post-sorption DOM were used to identify

234

compounds that were selectively removed from the solution due to adsorptive molecular

235

fractionation of the bulk DOM (Fig. 3). After post-spectral processing to identify peaks present

236

in at least two out of three analytical replicates, each treatment contained > 2,000 mass spectral

237

peaks that were assigned elemental compositions in the range 220