The Structural Arrangement and Relative ... - ACS Publications

Oct 8, 2018 - Skolkovo Institute of Science and Technology, 143025 Skolkovo, Moscow region, Russia. ∥. Institute for Energy Problems of Chemical Phy...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of Sunderland

Environmental Measurements Methods

The Structural Arrangement and Relative Abundance of Aliphatic Units May Effect Long-Wave Absorbance of Natural Organic Matter as Revealed by 1H NMR Spectroscopy Irina V. Perminova, Evgeny A. Shirshin, Andrey I. Konstantinov, Alexander Ya. Zherebker, Ivan V. Dubinenkov, Vasily A. Lebedev, Natalia A. Kulikova, Eugene N. Nikolaev, Ekaterina Bulygina, and Robert Max Holmes Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b01029 • Publication Date (Web): 08 Oct 2018 Downloaded from http://pubs.acs.org on October 9, 2018

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

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

Page 1 of 33

Environmental Science & Technology

1

The Structural Arrangement and Relative Abundance of Aliphatic Units May

2

Effect Long-Wave Absorbance of Natural Organic Matter as Revealed by 1H

3

NMR Spectroscopy

4 5

I. V. Perminova1*,

E.A. Shirshin,2

A. I. Konstantinov,1

6

I. V. Dubinenkov,1 N. A. Kulikova, 1,6,7 E.N. Nikolaev,3,4 E. Bulygina,8 R. M. Holmes8

7

1

8

Moscow, Russia; 2Department of Physics, Lomonosov Moscow State University, Leninskie Gory 1-

9

2, 119991 Moscow, Russia; 3Skolkovo Institute of Science and Technology, 143025 Skolkovo,

10

Moscow region, Russia, 4Institute for Energy Problems of Chemical Physics of RAS, Leninskij pr.

11

38 - 2, 119334 Moscow, Russia 5Department of Materials Science, Lomonosov Moscow State

12

University, Leninskie Gory 1-73, 199991 Moscow, Russia;

13

Lomonosov Moscow State University, Leninskie Gory 1-12, 199991 Moscow, Russia; 7Bach

14

Institute of Biochemistry of RAS, Federal Research Center “Biotechnology”, Leninskij pr. 33–2,

15

119071 Moscow, Russia 8Woods Hole Research Center, 149 Woods Hole Rd, Falmouth,

16

Massachusetts 02540, United States

V. A. Lebedev,1,5

Department of Chemistry, Lomonosov Moscow State University, Leninskie Gory 1-3, 119991

17 18

A. Zherebker,1,3,4

*E-mail: [email protected]

19

ACS Paragon Plus Environment

6

Department of Soil Science,

Environmental Science & Technology

Page 2 of 33

20

Abstract

21

The objective of this study was to shed light on structural features which underlay intensity of long

22

wave absorbance of natural organic matter (NOM) using 1H NMR spectroscopy. For this purpose, a

23

set of the NOM samples was assembled from arctic and non-arctic sampling sites (the Kolyma river

24

basin and Moscow region, respectively). It was to assure a substantial difference in the humification

25

degree of the isolated organic matter - the biogeochemical proxy of the long-wave absorbance of

26

NOM. The assembled NOM set was analyzed using solution-state 1H NMR spectroscopy. The

27

distribution of both backbone and exchangeable protons was determined using acquisition of spectra

28

in three different solvents. The substantially higher contribution of non-functionalized aliphatic

29

moieties CHn (e.g., materials derived from linear terpenoids, MDLT) in the arctic NOM samples

30

was revealed as compared to the non-arctic ones. The latter were characterized with the higher

31

content of CHα protons adjacent to electron-withdrawing groups which belong to carboxyl rich

32

alicyclic moieties (CRAMs) or to aromatic constituents of NOM. We have calculated a ratio of CHn

33

to CHα protons as a structural descriptor which showed significant inverse correlation to intensity of

34

long wave absorbance assessed with a use of E4/E6 ratio and the slope of absorption spectrum. The

35

steric hindrance of aromatic chromophoric groups of the NOM ensemble by bulky non-

36

functionalized aliphatic moieties (e.g., MDLT) was set as a hypothesis for explanation of this

37

phenomenon. The bulky aliphatics might increase a distance between the interacting groups resulting

38

in inhibition of electronic (e.g., charge-transfer) interactions in the NOM ensemble. The obtained

39

relationships were further explored using Fourier transform mass spectrometry as complementary

40

technique to 1H NMR spectroscopy. The data obtained on correlation of molecular composition of

41

NOM with 1H NMR data and optical properties were very supportive of our hypothesis that

42

capabilities of NOM ensemble of charge transfer interactions can be dependent on structural

43

arrangement and relative abundance of non-absorbing aliphatic moieties. ACS Paragon Plus Environment

Page 3 of 33

Environmental Science & Technology

44 45

1. Introduction

46

Natural organic matter (NOM) undergoes intense abiotic and microbial transformation in the

47

environment. Its rate and depth depends on different environmental conditions including type and

48

quantity of biomass, temperature, biological activity, humidity, and others.1 The progressive

49

oxidation of peptides, lipids, and polysaccharides leads to accumulation of structurally diverse

50

aromatics-rich compounds bearing both donor and acceptor groups, whose interactions give rise to a

51

characteristic “long-wave tail” of the NOM absorption spectrum.2 The more extended system of

52

interacting chromophores provides for the higher long-wave absorption, which has been intensively

53

used in soil chemistry to characterize “humification degree” of soil organic matter with a ratio of

54

absorbance values at 465 and 665 nm (E4/E6 ratio).4-6 Two different photophysical models are

55

currently under discussion, which have been developed to describe characteristic optical properties

56

of NOM: one of them considers NOM as a system of non-interacting chromophores, it will be

57

referred to hereinafter as “superposition model”, whereas the other one implies charge-transfer

58

interactions between NOM chromophores – “charge-transfer (CT) model”.7-9 The both models are

59

widely used to predict optical properties of NOM, which largely impact aquatic photochemistry and

60

microbial activity.10 Recently, new arguments have been set forth against the charge-transfer model.9

61

However, they have not been backed up by the corresponding structural studies.

62

Given that important prerequisite for charge-transfer interactions is a close spatial proximity

63

of the donor and acceptor,11 it was of interest to find a structural parameter of NOM, which could be,

64

on one side, related to spatial arrangement of donor and acceptor centers in supramolecular ensemble

65

of NOM, and, on the other side, reasonably easy measurable. The major structural patterns of NOM

66

include oxidation products of lipids, lignins, tannins, cellulose, and peptides.12 Hence, the major

67

donating centers might be represented by lignin-derived phenolic units, whereas the major acceptor ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 33

68

centers might be tannin-derived quinones – or carboxyl enriched rings.13 Both types of groups are

69

aromatic, and their spatial proximity may be ruled out by the presence of bulky substituents, e.g.,

70

branched, multiple, or annealed aliphatic groups. This effect was studied in detail by Rathore et al.14

71

who demonstrated that a group of three (1,3,5) ethyl substituents was sufficient for the complete

72

steric inhibition of the face of a benzenoid (donor) chromophore for intermolecular association by a

73

π–acceptor. The authors concluded that the difference of 0.9 Å between the separation of 3.6 Å in

74

the CT complex [hexamethylbenzene, chloranil (CA)] and the predicted separation of 4.5 Å in

75

[hexaethylbenzene, CA] could represent a “gray” area in which very weak, but visually (color) and

76

spectrally (CT) observable association may be apparent.

77

The other type of steric hindrance was recently envisioned by Chen et al15 for theoretic

78

explanation of red shift absorption of eumelanin –skin pigment, which is also often considered as

79

humic substances (HS) precursor.16,17 The authors reported that the excitonic interaction between

80

molecules in aggregates lead to the appearance of new (shifted) spectral bands that is a classical

81

effect in optical spectroscopy.15 However, this fact alone could not explain the exponential tail of

82

absorption: simple superposition of melanin aggregates spectra did not yield NOM-like absorption

83

spectrum. Chen et al. could achieve this effect by introducing the disorder into the system (i.e.

84

Gaussian distribution of distance and angle between the aggregates of melanin molecules), which

85

lead to the effective averaging of all possible shifted spectra and resulted in an overall exponentially

86

fading absorption.15 In case of NOM this would mean that aside from dependence on the amount and

87

quality of chromophoric groups (superposition model), the long-wave tailing of NOM spectra might

88

be governed by the quality and quantity of non-absorbing units which induce disorder in the system

89

of interacting chromophores (charge-transfer model). Of interest is that the both considered

90

phenomena imply a particular role of aliphatic units, which are both bulky and non-absorbing, in

91

imposing steric hindrance onto spatial arrangement of donor-acceptor pair which results in inhibition ACS Paragon Plus Environment

Page 5 of 33

Environmental Science & Technology

92

of their charge-transfer interactions. At the same time detailed assessment of aliphatic building

93

blocks was considered insignificant for describing UV/Vis behavior of NOM.

94

In case of the NOM system, these non-absorbing and bulky units can be easily attributed to

95

abundant aliphatic constituents such as materials derived from linear terpenoids (MDLT),18

96

carboxyl-rich alicyclic moieties (CRAM),19 as well as by branched sidechains of proteins, waxes,

97

suberins, and their associates.20-,22 Open chain terpenoids and proteins are much more amenable to

98

biodegradation and photooxidation as compared to alicyclic and branched structures.18 Hence,

99

progressive humification of biomass should be accompanied by an increase in contribution of

100

oxygenated alicyclic structures and branched methyl groups on the account of linear structures.21

101

This makes of particular diagnostic value the changes in aliphatic region of DOM.

102

The method which is most sensitive to the presence of saturated aliphatic moieties in the

103

NOM ensemble is 1H NMR spectroscopy. The main analytical challenge here comes from the

104

presence of both labile (exchangeable) and backbone (non-exchangeable) hydrogens in molecular

105

ensemble of NOM. The exchangeable protons are presented by O- and N-containing functional

106

groups: carboxylic (COOH), phenolic (ArO-H), amine (N-H), hydroxyl (O-H) etc., whereas the

107

backbone protons belong to non-substituted and functionalized methoxyl, methylene, methin groups,

108

and to aromatic rings.19,23-25 The major aliphatic compartments are presented by modified terpenes

109

such as the mentioned above MDLT18 and CRAM.19 In addition, side chains of peptides, O, and N-

110

acetyl polysaccharides, carbohydrates, ether bridges of lignins, etc. contribute to aliphatic region.21,22

111

As a result, the labile hydrogens of functionalized moieties can overlap with resonances of backbone

112

hydrogens when aprotic solvents are used for dissolution.26,27 At the same time, when NaOD/D2O is

113

used as a solvent, all exchangeable protons of NOM and HDO appear as a single coalescent peak in

114

the range of 4.5–4.8 ppm due to fast (on the NMR time scale) exchange.28

ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 33

115

For deriving quantitative data on proton distribution over NOM ensemble, a whole variety of

116

experimental techniques is described in the literature.23,25 This includes elimination of water from

117

the NOM sample, a use of aprotic solvents (e.g., DMSO-d6), a use of H-D exchange, and proton

118

downfield shift by addition of CD3COOD for enhancing resolution of 1H NMR spectrum with regard

119

to detection of non-exchangeable protons.26-28 Another approach uses suppression of the solvent

120

resonances in the NaOD/D2O solutions of NOM by selective presaturation of the HDO protons. 29

121

Improved solvent suppression techniques were applied for 1H NMR measurements of NOM in

122

D2O.30-33 In case of CD3OD (the solvent is widely used in 1H NMR studies of the solid phase

123

extraction (SPE) isolates of HS), presaturation is used to attenuate residual water.34,35 For excluding

124

residual water in case of DMSO-d6, thoroughly dried samples of NOM and HS were used.26,36

125

For elucidating a relationship between optical properties and structural features of both

126

aliphatic and aromatic constituents of the NOM ensemble, we have assembled a sample set

127

including the lesser humified materials from the high-latitude rivers in the Russian Arctic dominated

128

by the permafrost run-offs (the Kolyma River basin)37-39and the more humified materials from the

129

temperate-latitude rivers and lakes. We have measured proton distributions over exchangeable and

130

non-exchangeable structural moieties by acquiring triades of spectra in aprotic solvent (DMSO-d6)

131

alone, in DMSO-d6 in the presence of CF3COOD, and in D2O. We related the obtained 1H NMR

132

data to the optical properties measured with a use of UV-Vis spectroscopy. We backed up our

133

conclusions using complementary technique - Fourier transform ion cyclotron resonance mass

134

spectrometry (FTICR MS) for measuring molecular composition of the representative sample set.

135

2 Experimental

136

2.1. Sample set description

137

The set of high latitude permafrost-impacted NOM samples was isolated from the thaws of Yedoma

138

Ice Complex deposit and from the watersheds in the Kolyma River basin as a part of the “Polaris ACS Paragon Plus Environment

Page 7 of 33

Environmental Science & Technology

139

Project” expedition (www.thepolarisproject.org) during 2010-2013. The Yedoma deposit was

140

located in Duvanny Yar (DY),40 watersheds included different locations at the mainstream of the

141

Kolyma River, its four tributaries as well as lakes and streams. This “arctic” subset included 16

142

samples. The set of lower latitude samples (“non-arctic”) included four samples from the rivers and

143

lakes located in Moscow region (Russia) as well as the samples of fulvic acid (FA) and dissolved

144

organic matter (DOM) from the Suwannee River (SRFA and SRDOM, correspondingly) located in

145

the Georgia State, USA, – the standard samples of IHSS. The sampling locations and full sample

146

descriptions are shown in Fig. S1 and Table S1, respectively.

147

2.2. Isolation of NOM samples

148

For isolation of NOM from the Kolyma River basin watersheds, water was sampled by pre-cleaned

149

HDPE bottles. In-line 0.45 µm filters (AquaPrep 600, Gelman Sciences) were used for water

150

filtration. NOM isolation was carried out using Varian Bond Elut PPL cartridges (Mega Pack, 5 g),

151

Amberlite XAD-8 and XAD-16 resin (Rohm and Haas Ltd.). The Amberlite XAD-16 resin was used

152

as a feasible replacement to XAD-8 whose production is discontinued by the producer. Concentrated

153

HCl (HPLC grade) and solid NaOH (reagent grade) were used for samples acidification and for

154

preparing 0.1 M NaOH solution, respectively. Alkaline extracts were desalted on the cation

155

exchanging resin Amberlite 120R in H-form (Rohm and Haas Ltd.). HPLC-grade methanol was used

156

for elution of NOM from the PPL cartridges.

157

Water samples (20-40 liters) were filtered through 0.45 µm and acidified with concentrated

158

HCl until pH 2. Solid phase extraction (SPE) of NOM was conducted using Varian Bond Elut PPL

159

cartridges as described elsewhere.42,43 In brief, the PPL cartridge was rinsed with methanol; an

160

acidified water sample was discharged through the cartridge with a flow rate of 40 mL·min-1, the

161

cartridge was dried, and stored in the dark at 4ºC. The NOM was eluted using methanol at a flow

162

rate of 2 mL·min-1. The obtained extracts were dried on rotor evaporator. Column extraction of ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 33

163

NOM using XAD-8 was conducted according to the standard IHSS protocol.44 In brief, the filtered

164

water sample was acidified to pH 2 and discharged through the column packed with XAD-8 resin

165

(2x20 cm). The resin was eluted with 0.1 M NaOH. The alkaline extract was immediately passed

166

through the cation-exchanging resin in H-form. The desalted samples were rotor-evaporated to

167

dryness. The same protocol was applied to XAD-16 resin. All dry samples were stored at –20ºC in

168

the dark. All procedures are described in detail by Perminova et al.43 All sources of NOM samples,

169

the sorbents used for their isolation, and the attributions are listed in Table S1.

170

2.2. 1H NMR spectroscopy

171

Solvents and Reagents. DMSO-d6 had isotopic purity of 99.95 % (water content ≤ 0.02 %, Merck,

172

Darmstadt). Deuterated water (D2O, 99.95% D) and 40 % NaOD in D2O (99+% D) were purchased

173

from Aldrich (Milwaukee, WI). NaOD was used for facilitation of dissolution of some samples. All

174

reagents and solvents were used as received from the suppliers.

175

1

176

(Bruker, Germany) operating at 400 MHz 1H frequency. For measurements in DMSO-d6, the

177

samples were prepared as described by Kovalevskii et al.26 In brief, a weight of 15 mg of NOM

178

sample and 0.75 mL of anhydrous DMSO-d6 sealed in an ampoule were placed into a 5 mm NMR

179

tube and 13 mm pyrex tube, respectively, which were adjusted to each other as a Schlenk-type

180

apparatus. The apparatus was attached to a vacuum line and evacuated under high vacuum (< 5×10-2

181

mm Hg) for 6 hours to remove hygroscopic water from the NOM sample. The apparatus was

182

detached from the vacuum line and shaken to break DMSO-d6 ampoule and to drain the solvent into

183

an NMR tube with the dried NOM sample. The NMR tube was sealed and used for spectrum

184

acquisition. Then the tube was re-opened, added with 20 to 40 µL of CF3COOD and the spectrum

185

was reacquired. For D2O measurements, a weight of 15 mg NOM was dissolved in 0.6 mL of D2O.

186

All 1H NMR spectra were acquired in a 5 mm tube using 90 excitation pulses (acquisition time 3.2 s,

H NMR measurements. 1H NMR spectra were acquired using an Avance-400 NMR spectrometer

ACS Paragon Plus Environment

Page 9 of 33

Environmental Science & Technology

187

relaxation delay 2 s, ~100 scans). Fourier transformation, phase correction and integration were

188

performed using MestReC software (Mestrelab Research, USA).

189

The integration of spectra was conducted using the following assignments (ppm): 0-1.15 ppm

190

(CH3 – protons), 0.–2.05 – protons of alkyl chains (CHn – protons); 2.05–3.2 – protons of carbon

191

located in α-position to carboxyl group or to aromatic ring (CHα – protons) (X-CαH, X: COOH,

192

COOR, CAr), it might also carry resonances of amines and other functional groups; 3.2–6.0 – protons

193

of alkoxy and aliphatic hydroxyl groups, amines and amides (CHnO(N), CHn-O(N)H); 6.0–7.0 –

194

olephins, peptides, aminoacids, 7.0 -10.0 – aromatic and phenolic protons (CArH and CArOH; where

195

CAr corresponds to aromatic carbon atoms); 10.0–18.0 – protons of carboxyls (COOH).26,31,45

196

The spectral intensity of residual H atoms of DMSO-d6 at 2.5 ppm was subtracted from the

197

spectral integral of CHα- protons by using the NOM signature at δ H ~ 2.5 ppm as a continual

198

background. The spectral intensity of HDO peak was calculated as an integral from 4 to 5.5 ppm and

199

subtracted from the total integral in the range from 0 to 18 ppm. The values of CHn and CHα

200

integrals were used for calculation of CHn/CHα ratio.

201

2.4. UV-Vis spectroscopic measurements and data treatment

202

A weight of the NOM sample was dissolved in 0.03 M phosphate buffer at pH 6.8 to yield

203

concentration of 20 mg·L-1. UV-Vis spectra were recorded using Cary-50 spectrophotometer (Varian

204

Inc.) in the range from 200 to 800 nm. The recorded data were used to calculate E4/E6 ratio as a ratio

205

of absorbance values at 465 and 665 nm. This ratio is traditionally used to estimate humification

206

degree of organic matter.3,4 Another parameter calculated from the UV-Vis data is similar to the

207

chromaticity coefficient of HS.1 For this purpose the long wave (>500 nm) part of a UV-Vis

208

spectrum was approximated by eq (1):

209

A(λ) = A0·exp(-λ/Λ) + a,

ACS Paragon Plus Environment

(1)

Environmental Science & Technology

Page 10 of 33

210

where: A is the optical density; λ is the wavelength (nm); Λ (nm) is the slope of absorption in

211

semilogarithmic scale, and a is the offset connected to nonideal measurements of absorption spectra.

212

Eq 1 is based on logarithmic approximation of the long-wave tail of the absorbance spectra

213

described by the Urbach rule.46 The calculated slope value (Λ ) was used for evaluation of long-wave

214

tail of NOM.

215

2.5 Fourier transform ion cyclotron resonance mass spectrometry (FTICR MS)

216

FTICR mass-spectra were recorded using a hybrid LTQ FT Ultra (Thermo Electron Corp., Bremen,

217

Germany) mass-spectrometer equipped with a 7T superconducting magnet. Negative ions were

218

produced by an IonMax Electrospray ion source (Thermo Electron Corp., Bremen, Germany).

219

Needle voltage was -3kV. The direct infusion rate was 1 µL·min-1 Temperature of the desolvating

220

capillary was 200 oC. Broadband MS spectra (m/z 200-1000) were acquired in the FT ICR with a

221

resolution (R) of 400000 at m/z 400 by a consecutive summation of 300 scans. Internal calibration

222

was systematically conducted using the known peak series of NOM and the internal carboxylated

223

styrene calibrant.47,48 The NOM samples were dissolved in methanol (HPLC grade) at concentration

224

of 1 mg·mL-1. The data treatment is described in the Supporting Information,

225

3. Results and discussion

226

3.1. Solution state 1H NMR spectroscopy and UV-Vis on the NOM samples used in this study

227

1

228

of CF3COOD (TFA-d), and D2O. This was to observe full distribution of both backbone and

229

exchangeable protons in the NOM samples used in this study. The spectrum in DMSO-d6 is a

230

superposition of both exchangeable and backbone protons, whereas the spectra in DMSO-d6/TFA-d

231

and D2O contain only backbone protons and differ in the position of residual solvent protons: at 2.5

232

ppm (DMSO-d6) and 4.7 ppm (HDO). Acquisition of all three spectra enables normalization of

233

spectral intensity and calculation of full proton distribution. The “triads” of the obtained 1H NMR

H NMR spectra were acquired using different NMR solvents – DMSO-d6, DMSO-d6 in the presence

ACS Paragon Plus Environment

Page 11 of 33

Environmental Science & Technology

234

spectra typical for the “arctic” versus “non-arctic” NOM samples used in this study as well as for the

235

SPE-PPL versus the XAD-isolates from the same source are shown in Fig. 1. The “triads” were

236

acquired for 17 out of 27 samples used in this study. They are shown in full in Fig. S2 of the SI.

237

A

238

C

B

D

239

Figure 1. The typical triads of 1H NMR spectra acquired in DMSO-d6, DMSO-d6/TFA-d, and D2O

240

(highlighted in dark-purple, green, and blue, respectively) for the arctic NOM (Yedoma thaw, AHF-

241

DYP 11) (A) versus the non-arctic NOM (the Suwannee River DOM, SRDOM) (B), and the 1H

242

NMR spectra of the PPL isolate – AHF-LSP-12 (highlighted in blue) versus the XAD-8-isolate –

243

AHF-LSX12 (highlighted in orange) from the same source – the Sion Lake (Moscow Region)

244

acquired in D2O (C) and DMSO-d6 (D).

ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 33

245

The major difference among the spectra in the three solvents used can be seen in the region

246

of backbone aliphatics (0-3.2 ppm) occupied by the most “rigid” or sterically encumbered building

247

blocks of NOM, nominally, by MDLT and CRAM.18,19 In case of the typical “arctic” sample (Fig.

248

1A) it is largely dominated by non-functionalized alkyl protons of MDLT and other lipids signaling

249

from 0 to 2.05 ppm (CHn). Their contribution exceeds substantially intensity of downshifted CHα

250

protons usually attributed to CRAM (2.05 -3.2 ppm). In the shown representative of the “non-arctic”

251

samples (Fig. 1B), the intensity of CHn protons is drastically reduced (it becomes comparable with

252

CHα) – mostly, on the account of protons of carbohydrates and peptides signaling from 3.2 to 6

253

ppm. These trends might be indicative of the higher content of the least altered aliphatic units into

254

the structure of arctic versus non-arctic NOM samples, which is consistent with the differences in

255

their genesis: the low temperatures of the permafrost-dominated Kolyma region preclude rapid

256

decomposition of organic matter, whereas the elevated temperatures of temperate and other non-

257

arctic regions accelerate transformation processes of organic matter. However, except for the genetic

258

features, there is a purely analytical factor, which could also impact greatly the structural

259

composition of the NOM isolates: this is sorbent selectivity, which was considered in detail in our

260

previous publication.43 Its impact on structural features of the NOM isolates extracted from the same

261

source with a use of different sorbents is shown in Figs. 1C,D on the example of the NOM samples

262

extracted from the same source – the Sion Lake located in the temperate latitudes of the Moscow

263

region – using the SPE PPL cartridge and the Amberlite XAD-8 resin. The 1H NMR spectra for the

264

both samples were acquired in two different solvents – D2O and DMSO-d6, which are shown in Figs.

265

1C and 1D, respectively. It can be seen that the PPL isolate is characterized with the higher content

266

of aliphatic units as compared to the XAD-8 isolate: the latter, on contrary, is enriched with the

267

aromatic units. This is in line with the trends which were revealed in our previous studies.43 Hence,

268

the particular caution should be exercised upon comparing the NOM isolates, which were obtained ACS Paragon Plus Environment

Page 13 of 33

Environmental Science & Technology

269

with a use of different sorbents. This topic was studied and described in detail by Li et al.35 The

270

quantitative assessment of the proton distribution in all NOM samples used in this study is shown in

271

Fig. 2 for DMSO-d6 and D2O, the corresponding numeric values are given in Tables S2-S4 (SI). The

272

full size 1H NMR spectra in all three NMR solvents are shown in Figs. S3-S5 in the SI.

273

A

274

Figure 2. 1H NMR spectra acquired in DMSO-d6 (left panel, A) and D2O (right panel, B) with the

275

color-coded integral intervals, in ppm: 0-1.15 (yellow), 0-2.05 (light brown), 2.05-3.2 (violet, 3.2 -6

276

(purple), 6-7 (turquoise), 7-10 (dark green), 10-18 (light green).

277

As follows from Fig. 2, the largest amount of hydrogen in the NOM samples under study belongs to

278

non-functionalized aliphatic building blocks (e.g., MDLT),51 which are highlighted in brown. Their

279

content reaches up to 40% of the total (labile and backbone) hydrogen (Fig. 2A), and the consistently

280

higher values are observed in the arctic (30-40%) versus non-arctic (20 to 30%) isolates. The

281

respective values account for 40-50% and 30-35% of non-exchangeable hydrogen as follows from

ACS Paragon Plus Environment

B

Environmental Science & Technology

Page 14 of 33

282

D2O and TFA spectra (Fig, 2B and Table S4, respectively). Given the high relevance of branched

283

aliphatic units for the steric hindrance effect, we estimated the amount of methyl groups (shown in

284

yellow in Fig. 2) in the samples under study. It accounted for about 15% of the total hydrogen for the

285

arctic samples and did not exceed 10 % in the non-arctic ones. The corresponding ratios of CH3 to

286

(CH2+CH) and (CH2+CH+CHα) protons varied from 1:(1.6 – 1.8) and 1:(3.5-3.7), respectively, for

287

arctic NOM, to 1:(2.1-2.3) and 1:(4.5-5.5) for SRDOM and other non-arctic NOM. The obtained

288

values are indicative of rather extended methin-methylene core of the samples under study, in

289

particular, in those of arctic origin. Its branching degree is lower than in linear terpenoids (e.g.

290

squalene) or bacterial hopanoids, where the corresponding ratios could be estimated from 1:1.5 to

291

1:1.49 This could be connected to substantial contribution of non-branched lipids such as waxes and

292

fatty acids into terpenoid-like material, maybe, in the form of associates of hopanoids and cell-wall

293

lipids49 or with other methylene-rich bacterial compounds, e.g. branched glycerol dialkyl glycerol

294

tetraethers (brGDGT), which are important biomarkers of the high latitude OM.50

295

The content of CHα - protons varied not so remarkably – between 20 and 25 % of the total

296

hydrogen for all isolates and reached its minimum value (15%) for the Yedoma thaws (AHF-DYP10

297

and –DYP11). These protons are usually assigned to CRAM,19 which are produced by photo- and

298

microbial oxidation of terpenoids.18,21 The obtained trend looks rather consistent with the samples

299

origin: the Yedoma NOM was depleted with these highly processed NOM constituents. It is also in

300

line with the higher content of carbohydrate/peptide units (from 3.2. to 6 ppm) in these samples as

301

compared to the other arctic samples confirming their least degraded character.39-41 The CHα protons

302

also reflect amount of CH-substituted aromatic rings in the system, which might explain the highest

303

intensity of this interval for the non-arctic samples enriched with aromatics. This corroborates the

304

higher content of aromatic and COOH protons in these samples signaling from 6 to 14 ppm.

ACS Paragon Plus Environment

Page 15 of 33

Environmental Science & Technology

305

Summarizing the obtained results, the conclusion can be made that the “arctic” part of the

306

assembled NOM sample set was characterized with much higher contribution of methylene-enriched

307

non-functionalized aliphatic units usually assigned to MDLT compartments as compared to the non-

308

arctic isolates. The content of these CHn or “MDLT” protons in the arctic samples was up to a factor

309

of two higher than the content of CHα protons of CRAM- and aromatic constituents; in the non-

310

arctic samples it dropped to one. This motivated us for using a ratio of CHn to CHα protons (it can

311

be also considered as a proxy for MDLT to CRAM ratio) as a parameter reflecting the prevalence of

312

rigid non-oxidized hydrophobic core over the more polar, C-H acidic aliphatic units. Given that

313

these are non-functionalized rigid alicyclic moieties and branched chain proteins which contribute

314

the most into the amount of CHn protons, these very moieties might induce an effect of steric

315

hindrance by causing distortions in planar orientation needed for charge-transfer interactions of the

316

aromatic rings. In this case, a ratio of CHn to CHα protons might serve as a structural parameter

317

indicative of the size and “bulkiness” of the aliphatic part of the supramolecular ensemble of NOM

318

reflecting its feasibility for charge-transfer interactions and their contribution into the optical

319

properties of NOM. Hence, the NOM samples with the higher CHn/CHα ratio (e.g., the arctic subset

320

in case of this study) will be characterized with least pronounced charge-transfer interactions,

321

whereas the samples with the low CHn/CHα ratio (e.g., the non-arctic subset) will be characterized

322

with much higher contribution of charge-transfer interactions into their optical properties. On the

323

other side, this very ratio might serve as an NMR equivalent of humification degree of the NOM

324

ensemble while it reflects contribution of least transformed extended aliphatic structures versus

325

residual aliphatic structures adjacent to aromatic group of carboxyl resulting from progressive

326

oxidation of biomass during its transformation. For proving our working hypothesis on the

327

relationship between structural arrangement of aliphatics and its feasibility to interfere with charge-

ACS Paragon Plus Environment

Environmental Science & Technology

Page 16 of 33

328

transfer interactions in the NOM ensemble, we have measured the optic properties of the NOM set

329

used in this study.

330

3.2. Structure – optical properties relationships as deduced from 1H NMR and UV-Vis data

331

UV-Vis spectra of the samples under study were typical for NOM showing monotonous decrease of

332

absorbance along with an increase in a wavelength (Fig. 3A). They were characterized using the

333

E4/E6 value,3,4 which reflects steepness of the “long-wave tail” of the NOM spectrum, and the Λ

334

parameter, which is a reverse value of the slope of red edge absorption plotted in semilogarithmic

335

scale (Fig. 3A). The obtained E4/E6 and Λ values along with CHn/CHα ratio are given in Table 1.

336 337

Figure 3. UV-Vis spectra of NOM in semilogarithmic scale show good linearity in the range from

338

450 up to 700 nm, the red lines correspond to fits of absorption spectra using eq (1) (A);

339

visual representation of the Spearman correlation coefficient for the 1H NMR data and both optical

340

parameters used in this study: the E4/E6 ratio and the Λ parameter (n=32, p = 0.95) (B)

341

ACS Paragon Plus Environment

Page 17 of 33

Environmental Science & Technology

342

Table 1. The calculated CHn/CHα ratio and optical parameters characterizing long wave absorption

343

of the NOM samples used in this study Sample

Description of NOM source

Sorbent

the Yedoma deposit thaw – “arctic” NOM

CHn/CHα D2O

DMSOd6

E4/E6

Λ nm

AHF-DYP-10

Duvanny Yar (2010)

PPL

2.0

2.4

15.2

75

AHF-DYP-11

Duvanny Yar (2011)

PPL

2.1

2.0

13.7

70

AHF-DYP-13

Duvanny Yar (2013)

PPL

1.6

1.6

16.1

75

the Kolyma River basin watersheds (‘arctic” NOM) AHF-AOP-10 KR* estuary (2010) PPL

1.8

1.8

19.0

69

AHF-FPP-10

Flood plane stream (2010)

PPL

1.8

2.4

23.1

66

AHF-LTdP-10

Lake Tubdispanser (2010)

PPL

1.8

1.8

20.1

67

AHF-RK5P-10

KR near Chersky (2010)

PPL

-

1.8

26.1

70

AHF-RK6P-10

KR at Omolon River (2010)

PPL

1.9

2.0

16.3

73

AHF-RKP-12

KR mainstream (2012)

PPL

1.9

-

17.1

66

AHF-RKX8-11

KR near Chersky (2011)

XAD-8

-

1.4

15.1

83

AHF-RMAP-10

Malyi Anyui River (2010)

PPL

1.9

-

17.4

73

AHF-RMAX16-12

Malyi Anyui River (2012)

XAD-16

1.3

-

10.4

86

AHF-RMdP-10

Medvedka River (2010)

PPL

1.9

1.6

17.7

72

AHF-RPP-10

Panteleikha River (2010)

PPL

1.8

1.7

14.7

76

AHF-RPP-11

Panteleikha River (2011)

PPL

1.8

1.8

12.0

72

AHF-RPX8-11

Panteleikha River (2011)

XAD-8

1.8

1.7

10.0

83

AHF-RSkP-13

Sukharnaya River (2013)

PPL

1.8

-

16.2

66

the lower latitude watersheds (“non-arctic” NOM)

344

AHF-RIsX8-04

Istra River, Moscow, Russia

XAD-8

1.6

-

12.6

81

AHF-RMX8-97

Moscow River, Russia

XAD-8

1.6

-

9.6

88

AHF-LSP-12

Sion Lake, Moscow, Russia

PPL

1.5

1.5

11.0

87

AHF-LSX8-12

Sion Lake, Moscow, Russia

XAD-8

1.0

1.0

10.0

104

SRDOM

Suwannee River, IHSS

RO**

1.4

0.9

11.8

84

SRFA

Suwannee River, IHSS

XAD-8

1.1

1.1

12.7

79

*KR stands for the Kolyma River, **RO stands for reverse osmosis.

345 ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 33

346

As it can be deduced from Table 1, the highest E4/E6 values were characteristic for all PPL isolates

347

of arctic NOM reaching the maximum value of 26.1, whereas much lower values ranging from 9.6 to

348

12.7 were observed for the non-arctic samples (Table 1). In general, the obtained values of E4/E6

349

estimates corroborated well the reported ones for the NOM samples from different watersheds.

350

O’Driscoll et al. studied NOM from lakes and obtained E4/E6 ratios in the range from 4.5 to 16.52

351

Chin et al. reported the E4/E6 values for fulvic acids isolated from different rivers in the range from

352

9.7 (the Yakima River) to 21 (the Suwannee River).53 At the same time, the authors stated that the

353

optical density at 665 nm for the NOM samples could be quite low bringing about significant errors

354

in E4/E6 estimation.54 This was the same reason why we have used the Λ parameter, whose

355

measurement is based on approximation of the whole long-wave range of absorption spectra of

356

NOM. The corresponding Λ values varied from 66 to 91 for the arctic NOM used in this study, and

357

from 81 to 104 for the NOM samples isolated from the lower latitude watersheds (Table 1).

358

The revealed trends in optical properties indicated the smaller long-wave absorption of the

359

arctic versus non arctic NOM isolates used in this study. The obtained data set was used for

360

correlation analysis with the three sets of 1H NMR data calculated from the spectra measured in

361

D2O, DMSO-d6, and DMSO-d6/TFA-d. The results of correlation analysis are shown in Fig. 3B. The

362

Spearman correlation was applied due to the different sizes of the sample subsets. For all three 1H

363

NMR sets, the strongest statistically significant correlation with the Λ parameter (P=95%) was

364

observed for CHn/CHα (D2O –0.65, DMSO-d6 –0.82, DMSO-d6/TFA-d –0.81); CHn (D2O –0,6;

365

DMSO-d6 –0.83, DMSO-d6/TFA-d – 0.87), and HAr (D2O + 0,61; DMSO-d6 +0.79, DMSO-d6/TFA-

366

d +0.93). There was no statistically significant correlation observed for COOH. The similar

367

correlation trends were observed for the E4/E6 parameter. Given the specific feature of the NOM

368

sample set used in our study that the “arctic” samples were isolated mostly with a use of the SPE

369

PPL cartridges, whereas the”non-arctic” samples were isolated mostly using the Amberlite XAD ACS Paragon Plus Environment

Page 19 of 33

Environmental Science & Technology

370

resins, it was of interest to validate the obtained correlation relationship for the subset of samples

371

isolated with a use of the same isolation technique. For this purpose we conducted the same

372

correlation analysis for the PPL subset composed of 16 samples for all three 1H NMR data sets:

373

D2O, DMSO-d6, DMSO-d6/TFA-d (given in Tables S5-S6). The results of correlation analysis

374

(shown in Tables S7-S8) confirmed the trends obtained for the full data set: the closest correlation

375

with the Λ parameter (P=95%) was observed for CHn/CHα (D2O –-0.652, DMSO-d6 –0.639,

376

DMSO-d6/TFA-d –0.777).

377

In terms of charge-transfer interactions, the obtained relationships demonstrate that

378

enrichment of the NOM samples with the non-functionalized aliphatics weakens charge-transfer

379

interactions in the NOM ensemble, which is reflected in less pronounced long-wave tailing. The

380

obtained results support our initial hypothesis that bulky (or extended) aliphatic moieties might

381

prevent proximity of aromatic donor and acceptor centers which is necessary for formation of

382

charge-transfer complexes.

383

3.3. Molecular composition of the NOM samples under study and its relation to 1H NMR

384

structures and optical properties

385

For further exploration of the revealed structure-optics relationships, we have undertaken in depth

386

studies on molecular composition of the NOM samples used in this study using complementary to

387

1

388

FTICR MS including both arctic and non-arctic subsets. Fig 4 represents Van Krevelen (VK)

389

diagrams for the arctic (left panel) and non-arctic (right panel) NOM samples used in this study. All

390

four samples shown in Fig. 4 were characterized by the dominant contribution of unsaturated CHO

391

compounds with O/C≈0.5, which are considered as the most conservative part of NOM8 and could

392

be related to CRAM.19

H NMR spectroscopy technique – FTICR MS. Nine out of 22 NOM samples were analyzed using

ACS Paragon Plus Environment

Environmental Science & Technology

Page 20 of 33

393 394

Fig. 4. Van Krevelen diagrams for A) AHF-DYP10, B) AHF-AOP10, C) AHF-LSP12, D) SRDOM.

395

Dot size represents relative abundance of the corresponding peak in FTICR mass-spectra.

396

Attribution of the VK regions is according to55,56

397

The Yedoma-derived sample (AHF-DYP10) was characterized with the highest contribution

398

of aliphatic and peptide components (MDLT), whereas non-arctic sample SRDOM was enriched

399

with the most oxidized species (O/C>0.6) characteristic of hydrolyzable tannins. The latter were

400

much less abundant in the arctic NOM samples. These trends become even more obvious upon

401

segregation of the VK diagram field into the molecular classes using aromaticity index:55,56 the

402

major contribution was observed for the unsaturated species with H/C1.5) was higher in the arctic versus and non arctic NOM samples. For

404

relating the found trends in molecular composition to 1H NMR data and optical parameters, we have

405

identified common formulas present in all nine NOM samples (~650 out of ~4000 total) and used

406

their intensities for further correlation analysis57 (Fig. 5). ACS Paragon Plus Environment

Page 21 of 33

Environmental Science & Technology

407 408

Fig. 5. Correlation of the relative abundance of the common formulas of nine NOM samples

409

measured using FTICR MS with 1H NMR (DMSO-d6) and optical parameters. The colors designate

410

a positive (red) or negative (blue) value of Spearman correlation coefficient (P>95%).

411

It can be seen that 1H NMR data produce rather meaningful correlations with the abundance trends

412

in the molecular constituents of the NOM set used in this study. Positive significant correlations are

413

observed both for the CHn and CHα parameters with the molecular constituents located in the

414

[unsaturated and CRAM] region for the molecules with O/C0.5 region [aromatic and polyphenols] and negatively - to the aliphatic region. The content

418

of COOH was directly related to the population density of carboxyl-rich hydrolyzable tannin region.

419

Of particular interest are correlations with the structural and optical parameters proposed in our

420

study as characteristics of long-wave absorption: CHn/CHα and Λ, respectively. CHn/CHα directly

421

correlated with rather narrow suite of formulas located in the unsaturated region with O/C~0.5, and

422

was inversely proportional to abundance of a broad suite of formulas in the polyphenol region. The

ACS Paragon Plus Environment

Environmental Science & Technology

Page 22 of 33

423

similar trends were observed both for Λ and E4/E6 parameters: they were most sensitive to the

424

abundance of MDLT and aromatic regions.

425

The observed inverse relationship between the Λ parameter and the abundance of aliphatic

426

and low-oxidized unsaturated compounds might be seen as supportive of our working hypothesis on

427

the steric hindrance of charge-transfer interactions. Indeed, unsaturated compounds with H/C> 1 and

428

low O/C ratios can be assigned both to aromatic structures with aliphatic side-chains and to CRAM,

429

which are bulky alicyclic moieties and cause an increase in the distance between interacting aromatic

430

fluorophores. The direct correlation between Λ and the content of highly oxidized unsaturated

431

compounds, in turn, might indicate importance of the relative abundance of donors and acceptors in

432

the molecular ensemble of NOM for its optical properties. This is in line with the reported data on

433

reduction of carbonyl groups of DOM by sodium borohydride causing a loss of visible absorption

434

and inducing blue-shift of fluorescence spectra.58 The observed effects were explained as a

435

significant weakening of charge-transfer interactions. The similar trends might be deduced from the

436

data obtained in our study: firstly, the content of oxidized species– the major electron acceptors in

437

the NOM ensemble was significantly correlated with the Λ value (Table S9), which is proportional

438

to long-wave tailing caused by charge-transfer interactions; secondly, an increase in the CHn/CHα

439

ratio was directly proportional to the content of aliphatic and low-oxidized unsaturated molecules,

440

which cause a decrease in the Λ value. Of interest is also a direct relationship between the value of

441

CHα integrals in the 1H NMR spectra and the abundance of highly-oxidized unsaturated molecules

442

in the FTICR MS spectra (Table S9, Fig. S6), which might be indicative of oxidative transformation

443

of methyl groups during aromatization of terpene residues. In general, the obtained data provides

444

additional arguments in favor of steric mechanism of the impact which aliphatic compounds

445

(terpenoids and lipids) and substituents might impose on the visible absorbance of NOM.

ACS Paragon Plus Environment

Page 23 of 33

Environmental Science & Technology

446

3.4. Rationale on the role of aliphatic moieties in charge-transfer interactions in the

447

supramolecular ensemble of NOM

448

In the context of the obtained results and in line with the both previously discussed mechanisms of

449

the steric hindrance,14,15 it can be suggested that the prevalence of non-functionalized aliphatic

450

molecules/fragments (e.g. MDLTs) in the molecular ensemble of NOM might change distribution of

451

the relative positions of major chromophoric groups by shifting the respective intermolecular

452

distances to the values higher than 3.6 Å (Fig. 6). This might lead to a decrease in the strength of

453

charge-transfer interactions resulting in the lower absorption in the red edge of the spectrum. Increased fraction of aliphatic moieties

Aliphatic moieties

r r NOM NOM Chromophoric groups

454

Distance between chromophoric groups (r)

455

Figure 6. Distribution of distances between the chromophoric groups depends on the ratio of

456

aliphatic/aromatic moieties (CHn/CHα ratio) in the NOM supramolecular ensemble.

457 458

This looks overall plausible given that the MDLTs are mostly represented by sterically

459

encumbered modified linear and alicyclic terpenoids (e.g., squalene, sterans, hopanoids), other low-

460

functionalized lipids (e.g. waxes, fatty acids, branched glycerol dialkyl glycerol tetraethers),

461

branched side-chain peptides. Due to extended sp3 backbone these compounds are essentially

462

hydrophobic and CH2-chains might produce crystalline-like domains as it was previously envisioned

ACS Paragon Plus Environment

Environmental Science & Technology

Page 24 of 33

463

for soil HS.59 Mere hydrophobic interactions of these large aliphatic domains with aromatic rings

464

might increase the distance between the chromophoric groups to more than 3.6 Å leading to partial

465

inhibition of the charge-transfer interactions or to their complete extinction if the separation exceeds

466

4.5 Å.14 Progressive oxidation of the constituents of these hydrophobic domains (e.g., CRAM

467

formation) will increase solubility and reduce their size leading to more polar, functionalized and

468

labile (much less rigid) molecular ensemble of NOM, which will facilitate spatial proximity needed

469

for stronger charge-transfer interactions. The proposed conceptual model is in good agreement with

470

the existing data on NOM interactions with hydrophobic compounds with extended π-system such as

471

polycyclic aromatic hydrocarbons.60,61 It might also have direct implications for predicting

472

antioxidant and “pro-oxidant” activities of NOM and HS from the proton distribution in their

473

molecular ensemble.

474 475

Acknowledgements

476

Authors would like to acknowledge the outstanding contribution of the anonymous ES&T reviewer

477

whose insightful comments and suggestions enabled substantial improvement of the initial version

478

of the manuscript. AK, AZ, VL, NK and IP appreciate support of the Russian Science Foundation:

479

grant # 16-14-00167 (1H NMR spectroscopic analysis of the NOM samples and data interpretation),

480

AZ and EN – RSF grant #14-24-00114 (high-resolution mass spectrometry studies), AK, AZ, NK,

481

ID and IP appreciate support of the Russian Foundation for Basic Research (grant # 16-04-01753 -

482

isolation and primary characterization of the NOM samples), ES and IP – grant RFBR 15-05-09284

483

(optical studies of the DOM samples), ID, VL, EB and MH appreciate support of the NSF Polaris

484

project (sampling campaign in the Kolyma River basin). ES was supported by the Young

485

Investigator Research grant from the International Humic Substances Society. Authors are grateful to

486

developers of the open source Matplotlib software. ACS Paragon Plus Environment

Page 25 of 33

Environmental Science & Technology

487 488

The authors declare no conflict of interests

489

Supporting Information (SI) available for this manuscript contains extended information on the

490

samples (Fig. S1,Table S1), 1H NMR spectra and spectral integrals (Figs. 2-5), Tables S2-8,

491

description of FTICR MS data treatment and results (Table S9, Fig. S6). This material is available

492

free of charge via the Internet at http://pubs.acs.org.

493

References

494

1. G.R. Aiken, D.M. McKnight, R.L. Wershaw and P. McCarthy, Humic substances in soil,

495

sediment, and water: Geochemistry, Isolation, and Characterization, John Wiley & Sons, New

496

York, 1985.

497 498

2. F.J. Stevenson, Humus Chemistry: Genesis, Composition, Reactions, John Wiley & Sons, New York, 1994.

499

3. M.M. Kononova, Soil organic matter, Pergamon Press, Oxford, 1966.

500

4. Y. Chen, N. Senesi and M. Schnitzer, Information provided on humic substances by E4/E6

501

ratios, Soil Sci. Soc. Am. J., 1977, 41, 352-358.

502

5. L.T. Shirshova, E.A. Ghabbour and G. Davies, Spectroscopic characterization of humic

503

acid fractions isolated from soil using different extraction procedures, Geoderma, 2006, 133(3-4),

504

204–216.

505

6. R. Albrecht, J.L. Petit, G. Terrom and C. Périssol, Comparison between UV spectroscopy

506

and NIRS to assess humification process during sewage sludge and green wastes co- composting,

507

Bioresour. Technol., 2011, 102(6), 4495–4500.

508

7. C.M. Sharpless, N.V. Blough, The importance of charge-transfer interactions in

509

determining chromophoric dissolved organic matter (CDOM) optical and photochemical properties.

510

Environ. Sci.: Processes Impacts 2014, 16 (4), 654−671. ACS Paragon Plus Environment

Environmental Science & Technology

511 512

Page 26 of 33

8. R. Del Vecchio, N. V. Blough, On the origin of the optical properties of humic substances. Environ. Sci. Technol. 2004, 38 (14), 3885−3891.

513

9. G. McKay, J.A. Korak, P.R. Erickson, D. E. Latch, K. McNeill, and F. L. Rosario-Ortiz,

514

The case against charge transfer interactions in dissolved organic matter photophysics. Environ. Sci.

515

Technol. 2018. DOI: 10.1021/acs.est.7b03589

516

10. E. M. Thurman, Organic geochemistry of natural waters; Kluwer, 1985.

517

11. R. Foster, Organic Charge-Transfer Complexes; Academic Press: New York, 1969.

518

12. E. Kujawinski, M. Behn, Automated analysis of electrospray ionization fourier transform

519

ion cyclotron resonance mass spectra of natural organic matter. Anal. Chem. 2006, 78, 4363–4373.

520

13. D. T. Scott, D. M. McKnight, E. L. Blunt-Harris, S. E.Kolesar, D. R. Lovley, Quinone

521

moieties act as electron acceptors in the reduction of humic substances by humics-reducing

522

microorganisms. Environ. Sci. Tech. 1998, 32, 2984–2989.

523

14. R. Rathore, S.V. Lindeman and J.K. Kochi, Charge-Transfer Probes for Molecular

524

Recognition Via Steric Hindrance in Donor-Acceptor Pairs, J. Am. Chem. Soc., 1997, 119,

525

9393–9404.

526

15. C.T. Chen, C. Chuang, J. Cao, V. Ball, D. Ruch and M.J. Buehler, Excitonic effects from

527

geometric order and disorder explain broadband optical absorption in eumelanin, Nat. Commun.,

528

2014, 5, 3859.

529

16. G. Zonios, A. Dimou, I. Bassukas, D. Galaris, A. Tsolakidis and E. Kaxiras, Melanin

530

absorption spectroscopy: new method for noninvasive skin investigation and melanoma detection, J.

531

Biomed. Opt., 2008, 13(1), 014017.

532 533

17. M. Schnitzer, Y. K. Chan, Structural Characteristics of a Fungal Melanin and a Soil Humic Acid. Soil Sci. Soc. Am. J., 1985, 50(1), 67-71.

ACS Paragon Plus Environment

Page 27 of 33

Environmental Science & Technology

534

18. B. Lam, A. Baer, M. Alafe, B. Lefebvre, A. Moser, A. Williams and A.J. Simpson, Major

535

structural components in freshwater dissolved organic matter, Environ, Sci. Technol., 2007,

536

41(24), 8240-8247.

537

19. N. Hertkorn, R. Benner, M. Frommberger, P. Schmitt-Kopplin, M. Witt, K. Kaiser, A.

538

Kettrup and J.I. Hedges, Characterization of a Major Refractory Component of Marine Dissolved

539

Organic Matter, Geochim. Cosmochim. Acta, 2006, 70(12), 2990–3010

540 541

20. J.A. Leenheer, Progression from model structures to molecular structures of natural organic matter components, Annals of Environmental Science, 2007, 1, 57-68.

542

21. J.A. Leenheer, M.A. Nanny and C. McIntyre, Terpenoids as major precursors of

543

dissolved organic matter in landfill leachates, surface water, and groundwater, Environ. Sci.

544

Technol., 2003, 37, 2323–2331.

545 546

22. P.G. Hatcher, Aliphatic Structure of Humic Acids; a Clue to Their Origin. Org. Geochem, 1981, 3, 43-48.

547

23. N. Hertkorn and A. Kettrup, Molecular level structural analysis of natural organic matter

548

and of humic substances by multinuclear and higher dimensional NMR spectroscopy, in Use of

549

Humic Substances to Remediate Polluted Environments: From Theory to Practice, ed. I.V.

550

Perminova, K. Hatfield, Kirk, N. Hertkorn, Springer-Verlag, Berlin/Heidelberg, 2005, 391–435.

551

24. N. Hertkorn, A.B. Permin, I.V. Perminova, D.V. Kovalevskii, M.V. Yudov, A, Kettrup,

552

Comparative Analysis of Partial Structures of a Peat Humic and Fulvic Acid using One and Two

553

Dimensional Nuclear Magnetic Resonance Spectroscopy. J. Environ. Qual. 2002, 31, 375-387.

554

25. A.J. Simpson, D.J. McNally and M.J. Simpson, NMR spectroscopy in environmental

555

research: From molecular interactions to global processes, Prog. Nucl. Magn. Reson. Spectrosc.,

556

2011, 58, 97–175.

ACS Paragon Plus Environment

Environmental Science & Technology

Page 28 of 33

557

26. D.V. Kovalevskii, A.B. Permin, I.V. Perminova, D.V. Konnov, and V.S. Petrosyan,

558

Quantitative 1H NMR spectroscopic determination of exchangeable and backbone protons of humic

559

substances, Moscow University Chemistry Bulletin, Series 2 (Chemistry), 1999. 40(6), 375–380.

560

27. P. Ruggiero, F.S. Interesse and O. Sciacovelli, [1H] and [13C] NMR studies on the

561

importance of aromatic structures in fulvic and humic acids, Geochim. Cosmochim. Acta, 1979, 43,

562

1771–1775.

563

28. P. Ruggiero, F.S. Interesse, L. Cassidei and O. Sciacovelli, 1H NMR spectra of humic

564

and fulvic acids and their peracetic oxidation products, Geochim. Cosmochim. Acta, 1980, 44, 603–

565

609.

566

29. G.S.H. Lee, M.A. Wilson and B.R. Young, The application of the “WATERGATE”

567

suppression technique for analyzing humic substances by nuclear magnetic resonance, Organic

568

Geochemistry, 1998, 28, 549–559.

569 570 571 572

30. A.J. Simpson and S.A. Brown, Purge NMR: Effective and easy solvent suppression, J. Magn. Reson., 2005, 175(2), pp 340–346. 31. B. Lam and A.J. Simpson, Direct 1H NMR spectroscopy of dissolved organic matter in natural waters, Analyst, 2008, 133, 263–269.

573

32. A.S. Wozniak, R.U. Shelley, R.L. Sleighter, H.A.N. Abdulla, P.L. Morton, W.M.

574

Landing and P.G. Hatcher, Relationships among aerosol water soluble organic matter, iron and

575

aluminum in European, North African, and Marine air masses from the 2010 US GEOTRACES

576

cruise, Mar. Chem., 2013, 154, 24–33.

577

33. V. Kovacevic, A.J. Simpson and M.J. Simpson, 1H NMR-based metabolomics of

578

Daphnia magna responses after sub-lethal exposure to triclosan, carbamazepine and ibuprofen,

579

Comp. Biochem. Physiol., Part D: Genomics Proteomics, 2016, 19, 199–210.

ACS Paragon Plus Environment

Page 29 of 33

Environmental Science & Technology

580

26. F. Zhang, M. Harir, F. Moritz, J. Zhang, M. Witting, Y. Wu, P. Schmitt-Kopplin, A.

581

Fekete, A. Gaspar and N. Hertkorn, Molecular and structural characterization of dissolved organic

582

matter during and post cyanobacterial bloom in Taihu by combination of NMR spectroscopy and

583

FTICR mass spectrometry, Water Res., 2014, 57, 280–294.

584

35. Y. Li, M. Harir, J. Uhl, B. Kanawati, M. Lucio, K.S. Smirnov, B.P. Koch, P. Schmitt-

585

Kopplin and N. Hertkorn, How representative are dissolved organic matter (DOM) extracts? A

586

comprehensive study of sorbent selectivity for DOM isolation, Water Res., 2017, 116, 316–323.

587

36. J.-J. Wang, O. Pisani, L.H. Lin, O.O.Y. Lun, R.D. Bowden, K. Lajtha, A.J. Simpson and

588

M.J. Simpson, Long-term litter manipulation alters soil organic matter turnover in a temperate

589

deciduous forest, Sci. Total Environ., 2017, 607–608, 865–875.

590

37. B.G. Pautler, A.J. Simpson, D.J. McNally, S. F. Lamoureux and M.J. Simpson, Arctic

591

permafrost active layer detachments stimulate microbial activity and degradation of soil organic

592

matter, Environ. Sci. Technol., 2010, 44(11), 4076–4082.

593 594

38. K. Kaiser, M. Canedo-Oropeza, R. McMahon, R.M. W. Amon Origins and transformations of dissolved organic matter in large Arctic rivers Sci. Rep. 2017, 7, 13064.

595

39. J. Vonk, S. Tank, P. Mann,R. Spencer, C. Treat, R. Striegl, B. Abbott, K.Wickland

596

Biodegradability of dissolved organic carbon in permafrost soils and aquatic systems: a meta-

597

analysis. Biogeosciences. 2015, 12. 6915-6930.

598

40. J. Strauss, L. Schirrmeister, G. Grosse, D. Fortier, G. Hugelius, C. Knoblauch, V.

599

Romanovsky, C. Schädel, T. Schneider von Deimling, E. Schuur, D.G. Shmelev, M.Ulrich, A.

600

Veremeeva, Deep Yedoma Permafrost: A Synthesis of Depositional Characteristics and Carbon

601

Vulnerability. Earth-Science Reviews. 2017, .172, 75-86.

ACS Paragon Plus Environment

Environmental Science & Technology

Page 30 of 33

602

41. J E Vonk, P J Mann, K L Dowdy, A Davydova, S P Davydov, N Zimov, R G M Spencer,

603

E B Bulygina, T I Eglinton, R M Holmes, Dissolved organic carbon loss from Yedoma permafrost

604

amplified by ice wedge thaw. Environ. Res. Lett. 2013 8, 035023.

605

42. T. Dittmar, B.P. Koch, N. Hertkorn and G. Kattner, A simple and efficient method for the

606

solid-phase extraction of dissolved organic matter (SPE-DOM) from seawater, Limnol. Oceanogr.:

607

Methods, 2008, 6, 230–235.

608

43. I.V. Perminova, I.V. Dubinenkov, A.S. Kononikhin, A.I. Konstantinov, A.Ya. Zherebker,

609

M.A. Andzhushev, V.A. Lebedev, E. Bulygina, R.M. Holmes, Yu.I. Kostyukevich, I.A. Popov and

610

E.N. Nikolaev, Molecular mapping of sorbent selectivities with respect to isolation of arctic

611

dissolved organic matter as measured by Fourier transform mass spectrometry, Environ. Sci.

612

Technol., 2014, 48(13), 7461–7468.

613

44. G. R. Aiken, D. M. M cKnight, K. A. Thorn, E. M. Thurman, Isolation of hydrophilic

614

organic acids from water using nonionic macroporous resins. Org. Geochem. 1992, 18 (4), 567−573.

615

45. S.E.-M. Dvorski, M. Gonsior, N. Hertkorn, J. Uhl, H. Müller, C. Griebler and P. Schmitt-

616

Kopplin, Geochemistry of Dissolved Organic Matter in a Spatially Highly Resolved Groundwater

617

Petroleum Hydrocarbon Plume Cross-Section, Environ. Sci. Technol., 2016, 50(11), 5536–5546.

618 619

46. F. Urbach, The long-wavelength edge of photographic sensitivity and of the electronic absorption of solids, Phys. Rev., 1953, 92, 1324.

620

47. R.L. Sleighter, G.A. McKee, Z. Liu, P.G. Hatcher, Naturally Present Fatty Acids as

621

Internal Calibrants for Fourier Transform Mass Spectra of Dissolved Organic Matter, Limnol.

622

Ocean. Methods, 2008, 6(6), 246–253.

623

48. A. Zherebker, A.V. Turkova, Y. Kostyukevich, A. Kononikhin, K.V. Zaitsev, I.A. Popov,

624

E. Nikolaev and I.V. Perminova, Synthesis of Carboxylated Styrene Polymer for Internal Calibration

ACS Paragon Plus Environment

Page 31 of 33

Environmental Science & Technology

625

of Fourier Transform Ion Cyclotron Resonance Mass-Spectrometry of Humic Substances, Eur. J.

626

Mass Spectrom., 2017, 23(4), 156–161.

627

49. B.J. Belin, N. Busset, E. Giraud, A. Molinaro, A. Silipo and D.K. Newman, Hopanoid

628

lipids:

from

membranes

to

629

DOI:10.1038/nrmicro.2017.173.

plant–bacteria

interactions,

Nat.

Microbiol.,

2018,

630

50. F. Peterse, J. Vonk, R. Holmes, L. Giosan, N. Zimov, T. Eglinton, Branched glycerol

631

dialkyl glycerol tetraethers in Arctic lake sediments: Sources and implications for paleothermometry

632

at high latitudes. JGR: Biogeosciences, 2014, 119. 10.1002/2014JG002639.

633 634

51. B. Lam, A.J. Simpson, Direct H-1 NMR spectroscopy of dissolved organic matter in natural waters, Analyst, 2008, 133, 263–269.

635

52. N.J. O'Driscoll, S.D. Siciliano, D. Peak, R. Carignan and D.R.S. Lean, The influence of

636

forestry activity on the structure of dissolved organic matter in lakes: implications for mercury

637

photoreactions, Sci. Total Environ., 2006, 366(2), 880–893.

638 639

53. Y.P. Chin, G.R. Aiken and E. O'Loughlin, Molecular weight, polydispersity, and spectroscopic properties of aquatic humic substances, Environ. Sci. Technol., 1994, 28, 1853–1858.

640

54. N. Senesi, T.M. Miano, M.R. Provenzano, and G. Brunetti, Spectroscopic and

641

compositional comparative characterization of IHSS reference and standard fulvic and humic acids

642

of various origin, Sci. Total Environ., 1989, 81, 143–156.

643

55. B.P. Koch, T. Dittmar, From Mass to Structure: An Aromaticity Index for High-

644

Resolution Mass Data of Natural Organic Matter. Rapid Commun. Mass Spectrom. 2006, 20 (5),

645

926–932.

646

.56. A.M. Kellerman, T. Dittmar, D.N. Kothawala, L.J. Tranvik, Chemodiversity of

647

Dissolved Organic Matter in Lakes Driven by Climate and Hydrology. Nat. Commun. 2014, 5, 3804.

ACS Paragon Plus Environment

Environmental Science & Technology

Page 32 of 33

648

57. A.M. Kellerman, D.N. Kothawala, T.Dittmar and L.J. Tranvik, Persistence of dissolved

649

organic matter in lakes related to its molecular characteristics, Nat. Geosci., Published online: 25

650

May 2015, DOI: 10.1038/NGEO2440

651

58. J. Ma, R. Del Vecchio, K. S.Golanoski, E.S. Boyle, N.V. Blough. Optical properties of

652

humic substances and CDOM: Effects of borohydride reduction. Environ. Sci. Technol. 2010, 44,

653

5395–5402.

654

59. W.-G. Hu, J. Mao, B. Xing and K. Schmidt-Rohr, Poly(methylene) Crystallites in Humic

655

Substances Detected by Nuclear Magnetic Resonances. Environ. Sci. Technol., 2000, 34(3), 530–

656

534.

657

60. I.V. Perminova, N.Yu. Grechishcheva and V.S. Petrosyan, Relationships between

658

structure and binding affinity of humic substances for polycyclic aromatic hydrocarbons: relevance

659

of molecular descriptors, Environ. Sci. Technol., 1999, 33, 3781–3787.

660 661

61. B. Chefetz, A.P. Deshmukh, P. G. Hatcher, E. A. Guthrie, Pyrene Sorption by Natural Organic Matter, Environ. Sci. Technol., 2000, 34 (14), pp 2925–2930

ACS Paragon Plus Environment

Page 33 of 33

Environmental Science & Technology

271x114mm (96 x 96 DPI)

ACS Paragon Plus Environment