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