Relationships Between Dissolved Organic Matter Composition and

Jul 18, 2017 - Trace Element Removal in Distributed Drinking Water Treatment Systems by Cathodic H2O2 Production and UV Photolysis. James M. Barazesh ...
14 downloads 7 Views 11MB Size
Subscriber access provided by UNIV OF DURHAM

Article

Relationships Between Dissolved Organic Matter Composition and Photochemistry in Lakes of Diverse Trophic Status Andrew Chapin Maizel, Jing Li, and Christina K. Remucal Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b01270 • Publication Date (Web): 18 Jul 2017 Downloaded from http://pubs.acs.org on July 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 32

Environmental Science & Technology

1

Relationships Between Dissolved Organic

2

Matter Composition and Photochemistry in

3

Lakes of Diverse Trophic Status

4

Andrew C. Maizel,1 Jing Li,1 and Christina K. Remucal1, 2*

5

6 1

7

Department of Civil and Environmental Engineering

8

University of Wisconsin - Madison

9

Madison, Wisconsin 2

10

Environmental Chemistry and Technology Program

11

University of Wisconsin - Madison

12

Madison, Wisconsin

13 14 15

* Corresponding author address: 660 N. Park St., Madison, WI 53706; e-mail:

16

[email protected]; telephone: (608) 262-1820; fax: (608) 262-0454; Twitter: @remucal.

 

ACS Paragon Plus Environment

1

Environmental Science & Technology

17

Page 2 of 32

Abstract

18

The North Temperate Lakes-Long Term Ecological Research site includes seven

19

lakes in northern Wisconsin that vary in hydrology, trophic status, and landscape

20

position. We examine the molecular composition of dissolved organic matter (DOM)

21

within these lakes using Fourier transform-ion cyclotron resonance mass spectrometry

22

(FT-ICR MS) and quantify DOM photochemical activity using probe compounds.

23

Correlations between the relative intensity of individual molecular formulas and reactive

24

species production demonstrate the influence of DOM composition on photochemistry.

25

For example, highly aromatic, tannin-like formulas correlate positively with triplet

26

formation rates, but negatively with triplet quantum yields, as waters enriched in highly

27

aromatic formulas exhibit much higher rates of light absorption, but only slightly higher

28

rates of triplet production. While commonly utilized optical properties also correlate with

29

DOM composition, the ability of FT-ICR MS to characterize DOM subpopulations

30

provides unique insight into the mechanisms through which DOM source and

31

environmental processing determine composition and photochemical activity.

 

ACS Paragon Plus Environment

2

Page 3 of 32

32

Environmental Science & Technology

Introduction

33

Dissolved organic matter (DOM) is a compositionally diverse assembly of

34

molecules that is ubiquitous in natural waters. DOM contributes to numerous

35

environmental processes including carbon transport,1 redox cycling,2 and reactions with

36

environmental contaminants.3 Of special interest to lacustrine systems is the ability of

37

DOM to degrade xenobiotic compounds through the photochemical production of

38

reactive triplet states (3DOM)4 and reactive intermediates such as singlet oxygen (1O2)5

39

and hydroxyl radicals.6,7

40

DOM derives from dissimilar sources and, accordingly, is diverse in composition

41

and photochemical behavior.8-11 Broadly, DOM is considered allochthonous when derived

42

from degraded terrestrial plant material or autochthonous when derived from aquatic

43

microorganisms. Allochthonous DOM is typically higher in molecular weight,12 lower in

44

heteroatom content (e.g., N and S),13 and more aromatic than autochthonous DOM.14

45

DOM is further differentiated in natural systems by physical,15 chemical,16 and

46

biological17-19 processing. While it is challenging to apportion individual mechanisms to

47

DOM modification in environmental systems, DOM tends to become more aliphatic, less

48

oxidized, more diverse in elemental composition, and less chromophoric as it moves from

49

uplands to oceans.20-22 Similarly, formulas that are more aliphatic, less oxidized, and N-

50

rich are more persistent in lakes.23 DOM photochemistry also varies with source and environmental processing.

51 52

3

53

allochthonous DOM.24,25 Similarly, 1O2 quantum yields increase while 3DOM and 1O2

54

steady-state concentrations decrease as DOM moves from headwaters to the ocean.26,27

 

DOM and 1O2 quantum yields are generally higher in autochthonous DOM than

ACS Paragon Plus Environment

3

Environmental Science & Technology

Page 4 of 32

55

DOM photochemistry is commonly evaluated with the 3DOM probes sorbic acid (HDA)

56

and 2,4,6-trimethylphenol (TMP), and the 1O2 probe furfuryl alcohol (FFA).28 While

57

HDA and TMP directly measure 3DOM and 1O2 is formed from the reaction of 3DOM

58

and O2, there is evidence that these probes measure distinct 3DOM subpopulations.15,29

59

Recent reviews have proposed using a combination of multiple probes to better describe

60

3

DOM photoreactivity.30,31

61

The molecular composition of DOM is increasingly assessed with Fourier

62

transform-ion cyclotron resonance mass spectrometry (FT-ICR MS).32,33 The ability of

63

FT-ICR MS to identify individual molecular formulas in DOM has increased our

64

understanding of how physical,34 chemical,35 and biological processes36 modify DOM and

65

how DOM changes across environmental transects.20,21,37 For example, FT-ICR MS can

66

determine the composition of heteroatom-containing formulas and identify the source of

67

specific DOM subpopulations.16 Additionally, FT-ICR MS can detect bimodal

68

distributions in aromaticity and heteroatom composition that are not apparent using bulk

69

measurements.20 However, as with other methods for determining DOM composition,

70

FT-ICR MS is subject to analytical biases, such as the preferential detection of readily

71

ionized formulas.38

72

There is growing interest in relating DOM photochemistry and composition, as

73

xenobiotic compounds are increasingly identified in remote waters.39,40 However, 3DOM

74

production has not been previously related to molecular composition of DOM. Therefore,

75

we combine photochemical analyses with FT-ICR MS to investigate how DOM

76

composition controls photochemical activity in seven diverse lakes. Samples were

77

collected from the North Temperate Lakes-Long Term Environmental Research (NTL-

 

ACS Paragon Plus Environment

4

Page 5 of 32

Environmental Science & Technology

78

LTER) site in northern Wisconsin. Little is known about the DOM composition and

79

photoreactivity of NTL-LTER lakes, and the diversity of DOM sources and distinct UV-

80

vis absorbance spectra make them an ideal site to relate these attributes.

81 82

Methods

83

Sample Location. The NTL-LTER site includes two dystrophic (Crystal and

84

Trout Bogs), one mesotrophic (Allequash Lake), and four oligotrophic lakes (Big

85

Muskellunge, Crystal, Sparkling, and Trout Lakes). The NTL-LTER lakes vary in surface

86

area and water source, as well as in DOM source and concentration (Table S1).41-44 The

87

site is in a forested region of northern Wisconsin that is dominated by lakes.45

88

Samples were collected from the center of the NTL-LTER lakes between 02/1990

89

and 11/2014 and analyzed within 2-3 weeks for the concentration of dissolved organic

90

carbon ([DOC]) and by UV-visible spectroscopy (UV-vis). Additional samples (1–4 L)

91

were collected from the center of each lake in June 2015 and the edge of each lake in

92

August 2016 for photochemical analysis. Detailed information about sample collection,

93

analysis, and materials is available in Sections S1-S2.

94

UV-Vis Spectroscopy. Absorbance of samples collected in 2015 and 2016 was

95

determined using a Shimadzu UV-2401 PC, relative to a Milli-Q reference. SUVA254 is

96

the ratio of absorbance at 254 nm (A254) to [DOC].14 E2:E3 is the ratio of absorbance at

97

250 nm to 365 nm.22

98

Mass Spectrometry. DOM was extracted from the NTL-LTER samples collected

99

in August 2016 by solid phase extraction (SPE) and analyzed by FT-ICR MS. Details

100

about the SPE protocol are included in Section S3.46 Sample extracts were diluted 1:10 in

 

ACS Paragon Plus Environment

5

Environmental Science & Technology

Page 6 of 32

101

1:1 acetonitrile:Milli-Q and aspirated with 0.3 psi pressure into an electrospray source

102

with an applied voltage of -1.4 V. Analysis was performed with a SolariX XR 12T FT-

103

ICR MS (Bruker), coupled to a Triversa NanoMate sample delivery system (Advion), as

104

described previously.47 Details about instrument settings and formula identification are

105

available in Section S3. Formula relative intensities were determined by dividing the

106

intensity of each formula by the sum of intensities of all identified formulas in a sample.

107

Weighted average compositional values (e.g., H:Cw.avg) are the average of that

108

compositional value in all formulas, weighted by the relative intensity of each formula. It

109

should be noted that these average compositional values are determined only from

110

formulas identified by FT-ICR MS, rather than bulk elemental analysis. Average

111

compositional values exhibit similar qualitative trends to bulk analysis (e.g., elevated H:C

112

and N:C in autochthonous fulvic acid isolates)13 and provide insight into compositional

113

differences between natural DOM samples.21 Double bond equivalents per carbon

114

(DBE/C) was calculated as: {C – 0.5(H+Cl) + 0.5(N+S) +1}/C.48 Optical properties and

115

photochemical measurements are correlated by Pearson’s coefficients with the relative

116

intensities of individual formulas.

117

Photochemistry. Samples collected in 2015 and 2016 were irradiated in a

118

Rayonet photoreactor with UV-A bulbs (365 ± 9 nm).25,49 The initial photon flux of the

119

experimental apparatus was determined to be 7.35 x 10-8 E cm-2 s-1 with p-nitroanisole

120

(PNA)-pyridine actinomtery.25 Irradiation durations varied according to probe and sample

121

(10 – 360 minutes), and all irradiations were conducted in triplicate. HDA was added at

122

initial concentrations of 10, 100, 250, 500, and 1000 µM. 3DOM quantum yields (Φ3DOM),

123

formation rates (F3DOM), steady-state concentrations ([3DOM]ss) and first-order loss rate

 

ACS Paragon Plus Environment

6

Page 7 of 32

Environmental Science & Technology

124

constants (kd) were calculated from the isomerization rate of t,t-HDA and previously

125

estimated rate constants.25,50 Additionally, the observed loss rates of the 3DOM probe

126

TMP (initial concentration = 10 µM) were used to calculate 3DOM quantum yield

127

coefficients (fTMP) and steady-state concentrations of

128

estimated rate constants as described previously.51,52 The quantum yields (Φ1O2) and

129

steady-state concentrations of 1O2 ([1O2]ss) were calculated from the observed loss rates of

130

FFA (initial concentration = 10 µM) and previously measured rate constants.53,54

131

Quantum yields and fTMP were calculated with PNA-pyridine actinometry. All

132

calculations are described in detail in Section S4.55 FFA, PNA, TMP, and HDA isomer

133

concentrations were quantified by high-performance liquid chromatography as described

134

previously.25

3

DOM ([3DOM]ss,TMP) using

135 136

Results and Discussion

137

[DOC] and Optical Properties. [DOC] and UV-vis measurements taken from

138

1990 through 2014 distinguish the seven NTL-LTER lakes according to trophic status

139

and reflect the source and transformation of DOM within each lake. [DOC] represents the

140

concentration of bulk DOM, while A254 quantifies light absorption by chromophoric

141

DOM. DOM composition is described by SUVA254, which correlates with aromaticity,14

142

and E2:E3, which is inversely related to molecular weight.25,56 The bogs have the highest

143

mean [DOC], A254, and SUVA254, and lowest mean E2:E3 (Figure 1; Table S2). The high

144

SUVA254 and low E2:E3 values are typical of terrestrially-derived DOM,57-59 and are

145

reflective of the high fraction of DOM from adjacent wetlands (i.e., ~65%), high DOC

146

loading rates, shallow photic zones, and short hydraulic residence times (HRTs).44

 

ACS Paragon Plus Environment

7

Environmental Science & Technology

Page 8 of 32

147

In contrast, optical properties cannot be used to distinguish the source of DOM in

148

the oligotrophic lakes. These lakes have the lowest [DOC] and A254 due to longer HRTs

149

and lower carbon loading rates (Table S1),44 but are highly variable in SUVA254 and

150

E2:E3. The major sources of DOM to the oligotrophic lakes are precipitation, aerial

151

deposition, surface waters, and groundwater, rather than terrestrially-derived DOM.44

152

Additionally, more DOM is lost to mineralization compared with export and

153

sedimentation, indicating higher rates of DOM processing.44 The optical properties of the

154

oligotrophic lakes could indicate the presence of autochthonous DOM, which typically

155

has higher E2:E3 and lower SUVA254 than allochthonous DOM.14,22,59-61 However,

156

photobleaching of terrestrially-derived DOM decreases SUVA254 and increases E2:E3,15

157

and the presence of measureable autochthonous DOC in these lakes is debated.44,62

158

Collectively, the lower SUVA254 and higher E2:E3 values in oligotrophic lakes could

159

reflect DOM from autochthonous sources or allochthonous DOM that has undergone

160

extensive processing.

161

Mesotrophic Allequash Lake is unique among the seven NTL-LTER lakes.

162

Although the DOC loading rate is similar to the bogs, Allequash Lake has the shortest

163

HRT and, uniquely, exports more DOC than is lost to mineralization and sedimentation.44

164

Additionally, it receives more DOM from adjacent wetlands than the oligotrophic lakes,

165

but less than the bogs (i.e., ~10%).44 These factors produce DOM that shares some

166

properties with the bogs (high SUVA254, low E2:E3) and some properties with the

167

oligotrophic lakes (low [DOC] and A254).

168

Samples were collected in 2016 for photochemistry and FT-ICR MS

169

measurements, and have [DOC], A254, and E2:E3 slightly above historical means (Figure

 

ACS Paragon Plus Environment

8

Page 9 of 32

Environmental Science & Technology

170

1; Table S2). The elevated [DOC] and A254 of samples collected in 2016 may reflect their

171

collection from the lake edges, rather than centers. However, UV-vis compositional

172

measurements are generally within a standard deviation of historical means, indicating

173

that the photochemistry and FT-ICR MS results are reflective of representative DOM for

174

NTL-LTER lakes. Long-term trends in the optical properties of the NTL-LTER lakes

175

have been discussed previously.63

176

Molecular Composition of DOM. [DOC] and UV-vis measurements differ

177

between lakes according to trophic status, but the cause of these variations is unclear.

178

Analysis of DOM by FT-ICR MS provides more detailed compositional information.

179

3037 unique CHON0-1P0-1S0-1 formulas are identified in the NTL-LTER lakes (Table 1),

180

with 975 - 1791 unique formulas in individual lakes. In each lake, 60 – 81% of identified

181

formulas contain only CHO. The predominance of CHO formulas is typical of

182

allochthonous DOM.13,17 Similarly, the bulk compositional measurements of CHO

183

formulas are characteristic of lignin-like, terrestrially-derived molecules (H:Cw.avg = 1.05

184

– 1.36; O:Cw.avg = 0.46 – 0.56; Figures 2 and S2).21,64-66 While the terms “lignin-like” and

185

“tannin-like” are not definitive (i.e., a lignin-like formula may not necessarily be derived

186

from lignin), these elemental ratio cutoffs are useful in visualizing the FT-ICR MS data.  

187

H:Cw.avg increases and DBE/Cw.avg decreases from the bogs to Allequash Lake to the

188

oligotrophic lakes, confirming that the DOM in the bogs is more aromatic (Table 1).

189

O:Cw.avg is highest in the bogs and lower in Allequash Lake and the oligotrophic lakes.

190

CHON formulas comprise 12 – 31% of all identified formulas (Table 1) and

191

occupy similar van Krevelen space as CHO formulas from the same water body (Figures

192

2 and S3; Table S5). Higher fractions of formulas containing N are observed in the

 

ACS Paragon Plus Environment

9

Environmental Science & Technology

Page 10 of 32

193

oligotrophic lakes (21 – 31%) than in the bogs (12 – 15%) or Allequash Lake (17%). The

194

enrichment of CHON formulas in oligotrophic lakes could derive from the extensive

195

irradiation of allochthonous DOM or the presence of autochthonous DOM.17 However,

196

the compositional similarity of CHON and CHO formulas suggests a terrestrial source.

197

Few CHOP and CHOS formulas are detected (Figures S4 and S5; Tables 1 and

198

S6). Crystal Bog and Allequash Lake have the highest abundance of P-containing

199

formulas (9.4 and 6.1%, respectively), which could indicate preferential uptake of P-

200

containing formulas in the nutrient-poor lakes. Phospholipid formulas have been

201

identified in offshore coastal waters, but the CHOP formulas observed here have higher

202

O:C than expected for phospholipids.21 While few CHOS formulas are identified in NTL-

203

LTER lakes, some are observed at very high intensity, likely reflecting their ease of

204

ionization.38

205

Molecular Composition Differences Among Lakes. Bray-Curtis dissimilarity

206

analysis according to the log-weighted intensity of all identified formulas separates the

207

lakes by trophic status (Figure S6). Crystal and Trout Bogs are highly similar, with

208

Allequash Lake being slightly more dissimilar, while the oligotrophic lakes are highly

209

dissimilar to the bogs. This trend reflects the differences in the 25-year record of optical

210

properties. For example, SUVA254 measurements were highest in bogs and lowest in

211

Crystal Lake (Figure 1).

212

DOM compositional variation with trophic status is apparent in comparisons of

213

CHON0-1 formulas that are commonly identified in 5 or more lakes (n = 844; Figure S7).

214

The 481 CHO and 14 CHON formulas with higher relative intensity in the bogs and

215

Allequash Lake have lower average H:C (1.01 ± 0.24) and higher average O:C (0.54 ±

 

ACS Paragon Plus Environment

10

Page 11 of 32

Environmental Science & Technology

216

0.15) than 250 CHO and 99 CHON formulas enhanced in the oligotrophic lakes (H:C =

217

1.37 ± 0.02; O:C = 0.47 ± 0.14). 25% of formulas enhanced in oligotrophic lakes have

218

H:C values ≥1.5, which are typical of microbially-derived compounds,64,67 compared with

219