M. Ekrem
Karpuzcu‡
*ab,
Andrew J.
McCabe‡
a and
William A.
Arnold
a
aUniversity of Minnesota, Department of Civil, Environmental, and Geo-Engineering, Minneapolis, MN 55455, USA. E-mail: karpuzcu@itu.edu.tr
bIstanbul Technical University, Department of Environmental Engineering, Maslak 34469, Istanbul, Turkey
First published on 8th January 2016
Photochemical reactions involving a variety of photosensitizers contribute to the abiotic transformation of pesticides in prairie pothole lakes (PPLs). Despite the fact that triplet excited state dissolved organic matter (DOM) enhances phototransformation of pesticides by acting as a photosensitizer, it may also decrease the overall phototransformation rate through various mechanisms. In this study, the effect of DOM on the phototransformation of four commonly applied pesticides in four different PPL waters was investigated under simulated sunlight using photoexcited benzophenone-4-carboxylate as the oxidant with DOM serving as an anti-oxidant. For atrazine and mesotrione, a decrease in phototransformation rates was observed, while phototransformations of metolachlor and isoproturon were not affected by DOM inhibition. Phototransformation rates and the extent of inhibition/enhancement by DOM varied spatially and temporally across the wetlands studied. Characterization of DOM from the sites and different seasons suggested that the DOM type and variations in the DOM structure are important factors controlling phototransformation rates of pesticides in PPLs.
Environmental impactThe role of DOM as a photosensitizer in phototransformation of organic contaminants has been well known, but the potential of DOM as an antioxidant during photolysis has only recently been explored. By investigating whole prairie pothole waters and their seasonal variations, this study tries to directly address the environmental relevance of what so far has been a laboratory phenomenon. The inhibitory effect of DOM is especially important for the attenuation of pesticides that are recalcitrant to biotransformation and are inefficiently depleted by direct phototransformation. |
Previous studies reported that photochemical reactions play an important role in the abiotic transformation of a wide array of pesticides via direct and/or indirect pathways.16–19 Direct photolysis occurs when light energy is directly absorbed by a contaminant resulting in chemical transformation.20 Indirect photolysis occurs when another chemical species (photosensitizer) absorbs light and becomes electronically excited and consequently reacts directly with the contaminant of interest or produces “photochemically produced reactive intermediates” (PPRIs) which are capable of transforming the target compound. A variety of photosensitizers, including dissolved organic matter (DOM), nitrate/nitrite, and iron ions, as well as secondary species such as bicarbonate and halide ions are present in natural waters.20
DOM plays a critical role in aquatic health because it serves as an important component of the carbon cycle and plays a critical role in the transformation of aquatic contaminants, among several other roles. Inputs of water from agricultural fields,21 storm water,22,23 or wastewater effluent24 may affect the quantity and composition of the dissolved organic matter present. While DOM increases the phototransformation rates of pesticides by acting as a photosensitizer as mentioned above, it may also decrease the direct phototransformation rate by attenuating light as it proceeds through the water column (i.e., light screening), decrease indirect photolysis rates by quenching reactive species (PPRIs), or lead to reformation of the parent compound by acting as an antioxidant.25,26 The antioxidant pathways were first explored by Canonica and Laubscher.25 In their proposed model, the organic contaminant P reacts with an oxidizing radical denoted by Rad˙ or an excited triplet state denoted by 3Sens* to form a radical cation P˙+. Subsequently, P˙+ undergoes parallel reactions where it is either irreversibly oxidized to a product, Pox, or reduced back to its parent compound P by DOM (the antioxidant) forming an oxidized DOM radical, DOM˙+ (eqn (1)–(4)). The formed radical cation may further react with oxygen to form peroxy radicals.27
Rad˙ + P → Rad− + P˙+ | (1) |
3Sens* + P → Sens˙− + P˙+ | (2) |
P˙+ → Pox | (3) |
P˙+ + DOM → P + DOM˙+ | (4) |
Thus, reduction of P˙+ by DOM would decrease the overall rate of contaminant transformation. It has been suggested that the antioxidant properties of DOM manifest through a variety of organic moieties that DOM contains, such as phenolic groups.26,28
Most PPLs have large surface areas with relatively shallow depths (typically < 1.5 m), which greatly increases the potential of sunlight to penetrate the water column.10,29 Another feature that makes the surface waters of PPLs suitable for photosensitized reactions is that they often contain high levels of DOM.17 It has been previously shown that Suwannee River fulvic acid, used as a reference DOM, exhibited an inhibitory effect on the photolytic reactions of some organic contaminants such as anilines, cyanophenol and the antibiotic drug trimethoprim, while it enhanced the overall phototransformation of some other compounds such as isoproturon and 4-methylphenol.25 The main purposes of this study were to understand the potential of DOM as an antioxidant during phototransformation of pesticides in PPL waters and to determine the effect of season and DOM properties on the extent of inhibition or enhancement of phototransformation.
Fig. 1 Map of the prairie pothole region (PPR) of North America (map adapted from USGS). The inset map shows the Cottonwood Lake study area. The sampling sites are marked with triangles. |
Site/season | pH | DOC (mg L−1) | SUVA254 (L mg−1 m−1) | S 250–500 (nm−1) | S 275–295 (nm−1) | S 350–400 (nm−1) | S R | E 2:E3 | E 4:E6 | Fluorescence index (FI) |
---|---|---|---|---|---|---|---|---|---|---|
T9 summer | 9.70 | 36.86 | 7.70 | 0.0181 | 0.0184 | 0.0192 | 0.96 | 6.87 | 15.67 | 1.57 |
T9 fall | 8.10 | 20.70 | 6.53 | 0.0183 | 0.0187 | 0.0187 | 1.00 | 7.22 | 20.11 | 1.57 |
P7 summer | 8.94 | 32.59 | 4.17 | 0.0201 | 0.0244 | 0.0208 | 1.18 | 11.46 | 8.00 | 1.63 |
P7 fall | 8.72 | 29.78 | 4.31 | 0.0200 | 0.0241 | 0.0196 | 1.23 | 10.85 | 12.36 | 1.65 |
P1 summer | 8.75 | 38.41 | 4.36 | 0.0236 | 0.0276 | 0.0212 | 1.30 | 14.66 | 8.14 | 1.62 |
P1 fall | 8.99 | 35.26 | 4.29 | 0.0236 | 0.0284 | 0.0212 | 1.34 | 14.65 | 7.25 | 1.61 |
P8 summer | 8.94 | 23.45 | 5.96 | 0.0197 | 0.0217 | 0.0196 | 1.11 | 8.60 | 11.63 | 1.59 |
P8 fall | 8.24 | 21.85 | 5.39 | 0.0209 | 0.0227 | 0.0203 | 1.12 | 9.67 | 10.85 | 1.63 |
The photolytic reactions were conducted in an Atlas Suntest CPS+ solar simulator equipped with a 1500 W xenon arc lamp. The lamp was fitted with a UV-Suprax optical filter (passing wavelengths ranging from 290 to 800 nm), and the light intensity was set at 765 W m−2. The temperature of the reaction solutions was kept below 30 °C by circulating ambient air through the reaction chamber. The phototransformation experiments were carried out using 10 mL quartz test tubes. Test tubes were held at approximately a 30° angle from horizontal. All reactions were conducted in duplicate. Beside 10 μM of the target pesticide, aqueous solutions contained: treatment 1: borate buffer (blank); treatment 2: 40 μM CBBP and borate buffer; treatment 3: PPL water; treatment 4: PPL water and 40 μM CBBP; dark control: the same as treatment 4, but the test tubes were wrapped with aluminum foil to block light. These treatments follow the protocols outlined previously.25
Pesticide degradation was quantified with an Agilent 1100 HPLC equipped with a multiwavelength UV absorbance detector. A Discovery RP-amide C16 column (15 cm × 4.6 mm; 5 μm) was used for chromatographic separation. More technical details about the HPLC analysis are provided in the ESI.† Prior to irradiation experiments, UV-visible absorption spectra of the PPL water samples, target pesticides and CBBP were measured with a Shimadzu UV-1601PC spectrometer (Fig. S21, ESI†). Spectral slopes (S) were determined by fitting the spectra between 250 and 500 nm (corrected with absorbance at 700 nm) to an exponential trend line using Microsoft Excel®. Fluorescence spectra were collected using an Aqualog UV-800C with a 150 W ozone-free Xe source (Horiba, Jobin Yvon). Fluorescence indexes were computed from corrected spectra as the ratio of emission intensities at 470 nm to 520 nm at an excitation wavelength of 370 nm.31,32 Dissolved organic carbon (DOC) was measured on a Shimadzu TOC-L total organic carbon analyzer. Samples were acidified and purged with carbon-free air to remove any inorganic carbon prior to analysis. Iron(III) concentrations in the PPL samples were measured following a slightly modified version of the Ferrozine method of Viollier et al. (2000).33
Ct = C0e−kobst | (5) |
The rate constant for depletion of the target pesticide in the presence of both CBBP and DOM from the PPL water (treatment 4), kCB,PPL, was corrected by subtracting the rate constant for depletion in the presence of DOM from the PPL water only (treatment 3), kPPL, to compensate for pesticide loss due to DOM-induced phototransformation. When using kPPL in eqn (6), light screening by CBBP was taken into account by multiplying the kPPL value with a light screening factor, SDOM,CB (Table S2, ESI†). The rate constant for depletion of the target pesticide in the presence of CBBP (treatment 2), kCB, was corrected by subtracting the rate constant of the blank (treatment 1), kdirect, to compensate for the degradation caused by direct phototransformation rather than the presence of CBBP. When using kdirect in eqn (7), light screening by CBBP was taken into account and the kdirect value was corrected by multiplying it with a light screening factor, SCB (Table S2, ESI†). Screening factors were calculated as a ratio of the rates of light absorption in the presence and absence of the species responsible for screening. The details of light screening correction calculations are presented in the ESI.† Finally, the “Inhibition Factor”, IF,25 was defined as the ratio of the corrected rate constant in the presence of both CBBP and DOM, kcorr,CB,PPL, and the corrected rate constant in the presence of only CBBP, kcorr,CB, as described by eqn (6)–(8). When using kcorr,CB in eqn (8), light screening by DOM was taken into account by multiplying the kcorr,CB value with the light screening factor SCB,DOM (Table S2, ESI†). It should be noted that an IF value >1 indicates enhancement of phototransformation, while an IF value <1 indicates inhibition.
kcorr,CB,PPL = kCB,PPL − kPPL | (6) |
kcorr,CB = kCB − kdirect | (7) |
(8) |
Fig. 2 Phototransformation kinetics of atrazine in four photolysis treatments in (a) T9, (b) P7, (c) P1 and (d) P8 waters from the summer season. Error bars represent 95% confidence intervals. |
kindirect = kPPL − kdirect | (9) |
When using kdirect in eqn (9), light screening by DOM was taken into account by multiplying the kdirect value with the light screening factor SF,DOM (Table S2; ESI†). Pesticides showed seasonal and spatial variation in observed overall indirect phototransformation rates across the field sites (Fig. 3). Isoproturon had the fastest mean (±standard error of the mean) indirect phototransformation rate overall (kindirect, s−1 = 9.08(±0.70) × 10−5), followed by mesotrione (kindirect, s−1 = 9.46(±0.62) × 10−6), atrazine (kindirect, s−1 = 7.79(±1.79) × 10−6) and metolachlor (kindirect, s−1 = 3.67(±0.96) × 10−6). Zeng et al. (2013)17 observed similar results in P8 water using the same solar simulator where kindirect (s−1) was reported to be 7.56(±0.12) × 10−5 for isoproturon, 7.38(±0.15) × 10−6 for atrazine, 5.26(±0.10) × 10−6 for mesotrione, and 7.66(±0.15) × 10−6 for metolachlor. A potential explanation for the difference in the magnitude of the indirect photolysis reaction rates between isoproturon and the other three compounds is the difference in one electron oxidation potential (E1 or the Gibbs free energy change of electron transfer) of these compounds which has been suggested as a potential predictor of reaction rate constants with 3DOM* in photoactive systems.34 Using computational chemistry methods, E1 values for 70 compounds including the pesticides used in this study have been calculated by Arnold (2014).34 Isoproturon has the highest E1 value (−1.36 V vs. NHE) among the pesticides included in this study and therefore expected to have the highest indirect phototransformation rate. Although isoproturon reacts with carbonate radical35 and model triplet chromophores,36 it has been shown that the single electron oxidation observed in the presence of DOM is due to the reaction with 3DOM* rather than reactions with 1O2, CO3−˙ or ˙OH.19,29,37 Atrazine and metolachlor have similar E1 values (−2.41 V vs. NHE and −2.40 V vs. NHE, respectively) and these pesticides exhibited indirect phototransformation rates of the same magnitude. Although mesotrione had the lowest E1 value (−2.96 V vs. NHE), it exhibited a higher indirect phototransformation rate compared to atrazine and metolachlor. This might be due to other possible mechanisms described in previous studies.38 It has been reported that while H-donors present in DOM have an inhibitory effect on the photolysis of mesotrione, this effect is counterbalanced by the reaction of mesotrione with singlet oxygen produced by 3DOM*.38
A seasonal analysis of our results also indicates that indirect phototransformation rates were higher in summer compared to fall for atrazine and isoproturon at a significance level of 0.05, while there was no statistically significant seasonal difference for mesotrione and metolachlor (p = 0.40 and 0.22 for mesotrione and metolachlor, respectively). The seasonal differences may be explained by the loss of chromophores as a result of DOM photobleaching throughout the summer. Photobleaching depletes chromophoric dissolved organic matter (CDOM) which in turn may cause a reduction in indirect phototransformation rates of organic contaminants.39
The isoproturon phototransformation rate was enhanced in the presence of DOM for both seasons (fall IF = 1.40 ± 0.08, and summer IF = 1.20 ± 0.09). This result is in agreement with the expectation that the inhibitory effect of DOM would not play a significant role for electron-rich phenylureas.25 The metolachlor phototransformation rate was also enhanced in the presence of DOM for both seasons (fall IF = 1.24 ± 0.04, and summer IF = 1.67 ± 0.06). The observed seasonal and spatial variation in the inhibition (or enhancement) of phototransformation caused by DOM (Fig. 4) indicated possible differences in DOM properties and suggested the need for DOM characterization of the PPL waters at different seasons.
Previous studies have also documented temporal and spatial changes in the PPL water chemistry in terms of composition and reactivity of DOM.40–43 Both internal processes and external sources such as introduction of wastewater effluent, agricultural runoff, or stormwater flows may alter the DOM composition in a system. The exact chemical structure of aquatic DOM is unknown because DOM is comprised of a collection of numerous molecules originating from microbial (autochthonous) and terrestrial (allochtonous) precursor materials.44–46 Allochtonous DOM is characterized by its generally higher molecular weight, more hydrophobic nature, and greater aromaticity relative to autochthonous DOM.46 The two major classes of DOM also have quite different spectroscopic properties.46–48 Several spectral parameters have been defined to extract information from these spectra about CDOM properties. To track changes in the average molecular weight of DOM, an E2:E3 ratio (the ratio of absorption at 250 to 365 nm) was defined by de Haan and de Boer.49 They showed that the E2:E3 ratio was inversely proportional to the molecular size of DOM because of stronger light absorption by high molecular weight CDOM.49 Another absorption ratio, E4:E6 (the ratio of absorption at 465 to 665 nm), was reported to be inversely related to the aromaticity of CDOM.50 In many natural waters, absorption at 254 nm or 280 nm has been used as an indicator of aromaticity instead of the E4:E6 ratio, because there was often very little or no measurable absorption at 665 nm.46 Weishaar et al.48 showed that a parameter called specific UV absorbance or SUVA254 (UV absorption at 254 nm divided by dissolved organic carbon (DOC) concentration), correlated strongly with DOM aromaticity for a large number of humic substance isolates.
In addition to the parameters mentioned above, spectral slopes (S, nm−1) have been used to gain further insights into CDOM characterization such as determining the ratio of fulvic acids to humic acids.46S values have also been shown to correlate with the molecular weight of isolates of fulvic acids, however, their usefulness is limited by the fact that their values depend on the wavelength interval over which they are calculated.46 A wide range in S values is reported in the literature even for similar sample types and the lack of standardization has made comparisons of published S values difficult.46 To overcome the potential artifacts mentioned above, Helms et al. (2008)46 defined a dimensionless parameter called “slope ratio” or SR by calculating the ratio of the slope of the shorter wavelength region (275–295 nm) to that of the longer wavelength region (350–400 nm). The advantage of this approach is that it avoids the use of the spectral data near the instrumental detection limits and focuses on absorbance values that change significantly with the photochemical alteration of CDOM.46 Fluorescence spectroscopy has also been extensively used to characterize the source and properties of DOM.31,47 McKnight et al.47 defined the fluorescence index (FI) (the ratio of the emission intensity at 450 nm to 500 nm obtained with an excitation wavelength of 370 nm) and suggested that it can be used as a simple index for DOM characterization.
Based on the current literature, there is no consensus on which spectral parameter is the best for the characterization of DOM. In this study, all of the aforementioned spectral parameters were calculated for the PPL waters (Table 1, and Fig. S17–20 in the ESI†). It has been previously reported that iron is a source of interference in the spectroscopic analysis of DOM.51 Iron concentrations ([Fe2+] + [Fe3+]) in the PPL samples were below 3 μM, hence no iron related correction for the spectral parameters was necessary. It was determined that SR and the E2:E3 ratio exhibited the most reliable and consistent results in correlation with the Inhibition Factor (lowest p-values and highest correlation coefficients in Pearson correlation analyses, Table 2).
Compound | Season | Correlation with IF | S R | SUVA254 | E 2:E3 | E 4:E6 | S 250−500 | FI |
---|---|---|---|---|---|---|---|---|
Atrazine | Summer | Pearson correlation coefficient, r | 0.93* | −0.727 | 0.950* | 0.745 | 0.938 | 0.621 |
p | 0.049 | 0.273 | 0.050 | 0.255 | 0.062 | 0.379 | ||
Fall | Pearson correlation coefficient, r | −0.227 | 0.483 | −0.091 | 0.426 | −0.196 | −0.935 | |
p | 0.773 | 0.517 | 0.909 | 0.574 | 0.804 | 0.065 | ||
Isoproturon | Summer | Pearson correlation coefficient, r | −0.805 | 0.829 | −0647 | 0.844 | −0.685 | −0.679 |
p | 0.195 | 0.171 | 0.353 | 0.156 | 0.315 | 0.321 | ||
Fall | Pearson correlation coefficient, r | 0.989* | −0.950* | 0.949* | −0.817 | 0.837 | 0.511 | |
p | 0.011 | 0.050 | 0.050 | 0.183 | 0.163 | 0.489 | ||
Mesotrione | Summer | Pearson correlation coefficient, r | 0.123 | −0.500 | 0.099 | 0.470 | 0.140 | 0.635 |
p | 0.877 | 0.500 | 0.901 | 0.530 | 0.860 | 0.365 | ||
Fall | Pearson correlation coefficient, r | −0.311 | 0.565 | −0.166 | 0.471 | −0.246 | −0.967* | |
p | 0.689 | 0.435 | 0.834 | 0.529 | 0.754 | 0.033 | ||
Metolachlor | Summer | Pearson correlation coefficient, r | −0.454 | 0.252 | −0.272 | 0.289 | −0.289 | 0.001 |
p | 0.546 | 0.748 | 0.728 | 0.711 | 0.711 | 1.000 | ||
Fall | Pearson correlation coefficient, r | 0.985* | −0.899 | 0.996* | −0.933 | 0.964* | 0.463 | |
p | 0.015 | 0.101 | 0.004 | 0.067 | 0.036 | 0.537 |
There was a positive correlation between the inhibition factors and the SR and E2:E3 ratio values for atrazine in summer (Table 2). The inhibition factor for atrazine also had a relatively weak negative correlation (at a significance level of α = 0.10) with SUVA254 in summer samples (p = −0.73, r = −0.090). Inhibition factors for mesotrione did not exhibit any significant correlation with SR and E2:E3 ratio parameters, suggesting the role of other possible processes in its phototransformation as mentioned previously. Inhibition factors for both isoproturon and metolachlor were positively correlated with SR and the E2:E3 ratio in fall samples. Metolachlor also exhibited a positive correlation with the S250–500 parameter, while isoproturon exhibited a negative correlation with SUVA in fall (Table 2).
Previous studies reported evidence that photobleaching and shifts from high molecular weight (HMW) DOM to low molecular weight (LMW) DOM cause an increase in SR and S275–295 values.46 In this study, this phenomenon is supported by the fact that measured SR values increased in fall compared to summer at all PPL sites, suggesting the importance of photobleaching processes over the summer. The same relationship did not apply at all PPL sites for SUVA254 values, indicating that the observed seasonal differences in DOM behavior cannot be attributed to the changes in aromaticity alone. Together with the relationship of SUVA254 values to the aromaticity of DOM,48 our correlation analyses suggest that for the pesticide exhibiting the lowest IF values (atrazine), photobleaching and shift from HMW to LMW DOM were the dominant processes affecting the observed phototransformation inhibition by DOM followed by changes in aromaticity.
Thus, our results indicate that allochtonous, HMW and highly aromatic DOM (with higher SUVA and lower SR values) generally presents a relatively higher efficiency to inhibit excited triplet-induced oxidation than mainly autochthonous DOM. These results are in agreement with previous studies reporting that for the contaminants subject to inhibition of oxidation, their phototransformation induced by 3DOM* will vary greatly depending on the DOM type and will be inhibited more by highly aromatic DOM.26,44 These trends are largely consistent among the individual pesticides which are from different structural classes. More work is necessary to evaluate these correlations among more pesticides (across and within classes) as well as with additional water samples across a range of sites and seasons.
The results of this study are part of an effort to understand factors controlling the phototransformation of pesticides in PPL waters which still serve as collection points for non-point source agricultural runoff. Given the interconnectedness of water resources in the PPR, a better understanding of pesticide phototransformation in PPLs may shed light on future water and land use management strategies to mitigate the impact of pesticides on downstream surface water resources.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c5em00374a |
‡ Authors contributed equally. |
This journal is © The Royal Society of Chemistry 2016 |