Environ. Sci. Technol. 2009, 43, 5692–5698
Importance of Ultraviolet Radiation in the Photodemethylation of Methylmercury in Freshwater Ecosystems IGOR LEHNHERR* AND VINCENT L. ST. LOUIS Department of Biological Sciences, University of Alberta, Edmonton, Alberta, Canada, T6G 2E9
Received January 27, 2009. Revised manuscript received May 29, 2009. Accepted June 9, 2009.
Photodemethylation (PD) is thought to be the most important biogeochemical sink of methylmercury (MeHg) in freshwater lakes. However, we possess little mechanistic knowledge of this important biogeochemical process with regard to, for instance, the role of ultraviolet (UV) radiation versus visible light in mediating MeHg PD. This information is critical to correctly model MeHg PD at the whole-lake level, since wavelengths in the UV and visible regions of the solar spectrum are attenuated at very different rates in the water column of lakes. Furthermore, the established methodology for quantifying MeHg PD requires the addition of a MeHg spike, which often increases the concentration of ambient MeHg by 1 to 2 orders of magnitude; however, the assumption that the MeHg spike behaves like ambient MeHg has never been verified. We quantified MeHg PD rates using an isotopically enriched Me199Hg tracer added to lake waters already containing high concentrations of ambient MeHg, allowing us to simultaneously monitor the decomposition rate of the spike and ambient MeHg. Experiments were conducted at the Experimental Lakes Area to quantify the first-order rate constant (kpd) of MeHg PD in samples exposed to (1) full solar radiation, (2) UV-A and visible light (i.e., with UV-B blocked), or (3) visible light only. We demonstrate for the first time that the use of a MeHg spike to quantify PD rates is appropriate since spike and ambient MeHgsboth in samples with and without a spike of Me199Hgsare photodemethylated at the same rate. We also show that rates of MeHg PD are reduced by an order of magnitude in the absence of UV radiation and that to correctly model MeHg PD at the whole-lake scale, both UV and visible light mediated MeHg PD rates must be independently calculated using the light-specific rate constants (kpd-UVB, kpd-UVA, kpd-VIS). By examining modeled areal MeHg PD fluxes, we observed that UV radiation accounts for 58% and 79% of MeHg PD activity in a clear and colored lake, respectively. Finally, we demonstrate that correcting kpd-overall for the attenuation of solar radiation by Teflon bottles, which are normally used for MeHg PD experiments, increases the measured value of 3.69 × 10-3 m2 E-1 to 4.41 × 10-3 m2 E-1.
Introduction Methylmercury (MeHg) is a potent neurotoxin (1) that has the ability to bioaccumulate and biomagnify in aquatic food * Corresponding author phone: (780) 492-0900; fax: (780) 4929234; e-mail:
[email protected]. 5692
9
ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 43, NO. 15, 2009
webs (2). As a result, fish in numerous lakes worldwide contain elevated MeHg concentrations and the consumption of fish is the main pathway for human exposure to MeHg (3). While many studies have focused on identifying sources of MeHg to aquatic systems and quantifying the biogeochemical processes responsible for the production of MeHg (see 4 for a review), decomposition pathways of MeHg have not been thoroughly examined, although these processes are likely to be equally important in controlling the availability of MeHg to biota. MeHg can be microbially demethylated in both the sediments (5, 6) and the water column (7) of lakes, as well as photolytically decomposed in surface waters (8). In fact, photodemethylation (PD) is the most important sink of MeHg in clear, oligotrophic lakes in both boreal (9) and tundra (10) ecoregions. Laboratory experiments using both artificial and natural light sources have shown that MeHg is degraded by ultraviolet (UV; 100-400 nm) radiation (11-14). However, because of the strong positive correlation between photosynthetically active radiation (PAR, 400-700 nm) incident at different depths in the water column and MeHg PD rates measured at those depths, it was recently suggested that visible light, and not UV radiation, is primarily responsible for MeHg PD in lake waters (15), despite evidence to the contrary (9, 16). Previous research has also demonstrated that rates of MeHg PD are similar in filtered and unfiltered waters taken from the same lake (8, 15), suggesting that particulate matter does not play a role in this process beyond attenuating and scattering incident light. PD is often one of the principal fluxes considered in building lake-scale mass balance budgets of MeHg (9, 10, 17), and MeHg PD rates are typically quantified using the experimental methodology originally described by Sellers et al. (8). However, this method requires the addition of a spike of MeHg to water samples prior to incubation in the study lake because the low ambient MeHg concentrations naturally present in lake waters usually limit MeHg PD to rates that are not detectable. The addition of the MeHg spike often increases the concentration of aqueous MeHg by 1 or 2 orders of magnitude, potentially disrupting the chemical speciation of MeHg in lake waters, including partitioning between the particulate and dissolved phases and binding of ligands, which could in turn yield measured MeHg PD rates which are not representative of in situ rates. To evaluate whether MeHg PD rates measured using the methodology described above are representative of true rates, we examined MeHg PD in a boreal lake at the Experimental Lakes Area (ELA) in northwestern Ontario, Canada, by adding an isotopically enriched Me199Hg tracer to lake waters already high in ambient MeHg. This approach allowed us to simultaneously quantify PD rates of the added Me199Hg tracer and ambient MeHg to determine whether the Me199Hg tracer behaves like ambient MeHg or otherwise affects the PD of ambient MeHg. To examine the effect of UV radiation on MeHg PD rates, we quantified MeHg PD rates in samples exposed to (1) full solar radiation including UV-B (280-320 nm), UV-A (320-400 nm), and visible light (400-700 nm), (2) UV-A and visible light (i.e., with UV-B blocked), and (3) visible light only. Experiments were performed with filtered and unfiltered water to elucidate whether particulate matter has a direct effect on measured MeHg PD rates. Finally, we present a model which incorporates light-specific rate constants of MeHg PD with UV and visible light attenuation to calculate whole-lake MeHg PD fluxes. 10.1021/es9002923 CCC: $40.75
2009 American Chemical Society
Published on Web 07/01/2009
Methods Experimental Design and Sample Collection. The MeHg PD experiments described here were conducted July 15-22, 2006, at the ELA located in the temperate boreal ecoregion, using water collected from Lake 979. Lake 979 was formed by the experimental flooding of a wetland complex in 1993 (18) and therefore had high concentrations of both DOC (12.8 mg L-1) and MeHg (0.81 ng L-1) due to the decomposition of flooded organic peatland soils and the resulting stimulation of Hg methylation. Unfiltered epilimnetic water samples were collected into a 1-L Teflon bottle using Teflon tubing connected to a battery-powered gear pump (Greylor, PQ-12DC). Filtered samples were pumped through an additional in-line Teflon cartridge equipped with a 2.2 µm quartz microfiber filter that was previously muffled overnight at 500 °C. Samples were handled using trace-metal sampling procedures such as the “clean hands, dirty hands” protocol (19), and all sampling equipment was acid-washed and rinsed with MilliPore water (resistivity >18.2 MΩ cm) prior to use. Bulk 1 L water samples were amended with a 1000 µL spike of a ∼1.2 µg L-1 solution of isotopically enriched Me199Hg tracer in MilliPore water. The Me199Hg tracer was synthesized from 199HgO (92% purity, Trace Sciences International) using methylcobalamin (20). Spiked bulk water (200 mL) was distributed into each of five separate 250mL FEP-Teflon bottles, and bottles were randomly assigned to one of five incubation times (0, 1, 2, 4, or 7 days). A rate of MeHg PD was obtained for each five-point time-series (see Data Analysis below) and every time-series was performed in duplicate. Samples were incubated on the surface of Lake 239, located near the main research facilities at the ELA, inside a raft with different compartments (Figure S1). Two of the raft compartments were covered with either Mylar or Lee (model 226) film to control the quality of light to which samples were exposed. Mylar film has a cutoff wavelength of 320 nm while Lee film does not transmit light of wavelengths shorter than 400 nm (Figure S2); furthermore, the UV-blocking properties of the Mylar and Lee films were assessed before and after the experiment using a UV-vis spectrophotometer and remained unchanged over the short 7-day duration of the incubations. Therefore, MeHg PD rates in filtered and unfiltered lake water were quantified under the following light conditions: (1) full solar radiation, including UV-B, UV-A, and visible light, (2) UV-A and visible light (i.e., with Mylar film blocking UV-B radiation), and (3) visible light only (i.e., with Lee film blocking all UV). Duplicate time-series were performed for each possible combination of the light quality and filtration factors, and within each time-series, PD rates of ambient MeHg and spike Me199Hg were simultaneously quantified (spike versus ambient MeHg factor). A number of treatments were also replicated in the absence of light, by wrapping aluminum foil around sample bottles, or without the addition of Me199Hg to verify that the spike did not interfere with the measurement of in situ MeHg PD rates. Bottles collected at the end of the desired incubation period were immediately frozen and stored in the dark at -20 °C until sample analysis. Sample Analysis. Water samples were distilled and analyzed for MeHg using isotope-dilution gas chromatography inductively coupled plasma mass spectroscopy (GC-ICP-MS) (21, 22). Additional analytical details and QA/QC are provided in the Supporting Information. Data Analysis. A rate constant of MeHg PD was obtained for each incubation time-series by modeling the MeHg concentration data using first-order chemical kinetics: ln(MeHg)t ) ln(MeHg)0 - kt
(1)
where (MeHg)t and (MeHg)0 are the MeHg concentrations at time t and time 0, respectively, and k is the first-order rate constant. However, for a photochemical process, incident photon flux density is a better predictor of the progress of reaction than reaction-time itself. Therefore, ln(MeHg)t was regressed against the cumulative solar radiation flux, measured using a LiCor PAR sensor (LI-190SA Quantum Sensor) over the duration of the incubation time, t. The rate constant of MeHg PD, kpd, was then obtained from the slope of the characteristic first-order kinetic plot (Figures S3-S5): ln(MeHg)t ) ln(MeHg)0 (kpd × Cumulative PAR photon flux)
(2)
The value of kpd, obtained for each time-series, was input into a repeated-measures ANOVA model with light quality and filtration as the between-time-series factors and ambient versus spiked MeHg as the within-time-series factor. All statistical analyses were performed at R ) 0.05, using Systat (v. 11).
Results and Discussion Effect of Light Quality on MeHg PD. We found that MeHg in lake water was decomposed when exposed to solar radiation (Figure 1), as has been shown in previous studies (8, 15, 17), and that UV radiation is the principal driver of MeHg PD in surface waters. Measured rate constants of MeHg PD (kpd, m2 E-1) varied significantly with incident light quality (p < 0.0001), and decreased noticeably in the absence of UV radiation (Figure 1). In samples exposed to visible light only, kpd was 7-12 times smaller than in samples exposed to the full solar spectrum (Table 1, Tukey’s posthoc multiple mean comparison p < 0.0001), and 5-8 times smaller than in samples exposed to UV-A and visible light radiation (Table 1, Tukey’s p ) 0.0003). Furthermore, blocking UV-B radiation resulted in 22-35% decreases in kpd, suggesting that UV-B radiation enhances MeHg PD above and beyond the effect of UV-A (Table 1, Tukey’s p ) 0.015 and 0.023 for ambient and spike MeHg, respectively). However, this effect could in part be due to the fact that the Mylar film used to block UV-B radiation also blocked ∼25% of the UV-A radiation in the range of 320-375 nm. Finally, kpds measured in filtered water samples incubated in the dark were ∼70% lower than in samples incubated under visible light only. The mean rate constant ((standard error) for the decomposition of MeHg in the absence of lightsexpressed with respect to time as opposed to cumulative photon fluxswas (5.36 ( 0.99) × 10-3 d-1; however, under these experimental conditions, such low rates of MeHg decomposition were at the limit of detection and are associated with a high degree of uncertainty (four regressionsofln(MeHg)dark versusTime,p-values)0.004-0.36). Effect of Filtration and Particulate Matter on MeHg PD. Overall, filtering the water prior to incubation had no effect on the measured value of kpd (p ) 0.105; see also Figure 2). However, kpd tended to be higher in unfiltered samples for the experimental treatments performed in the absence of UV-B and under visible light only (Table 1). This tendency is somewhat surprising because we expected that if particulate matter did have an effect on MeHg PD, it would be to lower kpd due to increased scattering of incident light and decreased photoreactivity of particulate-bound MeHg. However, particle content in Lake 979 water was low and particulate bound MeHg in Lake 979 comprised only ∼10% of the total aqueous MeHg; therefore, any potential differences in either light-scattering or photoreactivity of MeHg between filtered and unfiltered samples were likely too small to effect kpd. Our results are consistent with previous reports (8, 15) that there is no difference between MeHg PD rates quantified in filtered and unfiltered water. Therefore, parVOL. 43, NO. 15, 2009 / ENVIRONMENTAL SCIENCE & TECHNOLOGY
9
5693
FIGURE 1. Ambient MeHg (A, B) and spike Me199Hg (C, D) concentrations in filtered (A, C) and unfiltered (B, D) samples over the course of their incubation under various light conditions. Error bars represent the difference between duplicate samples. Data points represent incubation times of 0, 1, 2, 4, and 7 days.
TABLE 1. Rate Constants, kpd, of MeHg PD (Mean ± Standard Error; × 10-3 m2 E-1) Measured for the Various Experimental Treatments kpd (× 10-3 m2 E-1) full spectrum
UV-A + visible
visible only
dark
filtered water ambient MeHg in spiked samples spiked Me199Hg ambient MeHg in unspiked samples
3.87a 3.82a 3.93 ( 0.02
2.31 ( 0.19 2.26 ( 0.20 N/A
0.294 ( 0.010 0.269 ( 0.047 N/A
0.137 ( 0.008 0.101 ( 0.023 N/A
unfiltered water ambient MeHg in spiked samples spiked Me199Hg ambient MeHg in unspiked samples
3.59 ( 0.08 3.69 ( 0.19 3.84 ( 0.18
2.81 ( 0.005 2.85 ( 0.06 N/A
0.527 ( 0.002 0.505 ( 0.019 N/A
N/A N/A N/A
3.69 ( 0.07 (n ) 7)
2.44 ( 0.18 (n ) 4)
0.280 ( 0.069 (n ) 4)
0.119c ( 0.020 (n ) 2)
best estimateb
a Due to the presence of one outlier value (Q-test p < 0.05) in the duplicate five-point time-series performed with filtered water and incubated under full solar radiation, the kpd value obtained from this time-series was not used in the computation of mean kpds. b Best estimate values of kpds were corrected for the dark-decomposition of MeHg by subtracting the kpd measured in the absence of light from all other kpds. c The rate constant of dark MeHg decomposition expressed with respect to time is (5.36 ( 0.99) × 10-3 d-1.
ticulate matter does not appear to play a critical role in MeHg PD beyond increasing the scattering and absorption of incident solar radiation in the water column. However, because our results suggest that rates are slightly higher in unfiltered water, we recommend quantifying MeHg PD in unfiltered samples for the purpose of modeling this process at the whole-lake level. Effect of Spike Addition on MeHg PD. There was no difference between kpds obtained for ambient MeHg and spike Me199Hg (p ) 0.802). However, the interaction term Spike versus Ambient MeHg * Filtration was almost significant (p ) 0.081) suggesting that in unfiltered samples, the spike Me199Hg was photodemethylated at a slightly higher rate than the ambient MeHg, while the reverse was true in filtered samples. This could be due to the spike Me199Hg not partitioning between the particulate and dissolved phases 5694
9
ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 43, NO. 15, 2009
in exactly the same manner as the ambient MeHg such that in unfiltered samples, more of the spike Me199Hg remained in the dissolved, and more reactive, phase. However, the difference in kpd resulting from a potential Spike versus Ambient MeHg * Filtration interaction is quite small (Table 1), and is perhaps better described as a statistical effect that is not relevant on an ecological scale. Additional experiments performed using samples that did not receive a spike of Me199Hg showed that rates of ambient MeHg PD were no different than in samples spiked with Me199Hg (Figure 2) and that this holds true in both filtered and unfiltered water (t test for sample means with unequal variances, p ) 0.46 and p ) 0.38 for filtered and unfiltered samples, respectively, and p ) 0.16 when filtered and unfiltered data are pooled; Table 1). Together, these results validate the assumption that the added isotopically enriched MeHg tracer
FIGURE 2. Photodemethylation of ambient MeHg in samples either spiked with 1.2 ng L-1 of Me199Hg or not amended with spike. Data for both filtered and unfiltered water samples are shown and error bars represent the difference between duplicate samples. Samples were collected following incubation periods of 0, 1, 2, 4, and 7 days under full solar radiation.
TABLE 2. Relative Contribution of UV-B, UV-A, and Visible Light to the Whole-Lake Areal MeHg PD Flux, and Depth of MeHg PD Activitya light attenuation coefficient (m-1) % of overall areal MeHg PD flux areal PD fluxb (ng MeHg m-2 d-1) 10% depthc (m) -3
kpd (× 10
2
-1
m E )
UV-B UV-A all UV
0.871 (1.25) 3.22 (2.16)
visible
0.313 (0.280)
L979 35 11
2.9
L239
L979
L239
L979
L239
L979
L239
14 5.4
10 (18) 69 (58) 79 (76)
7 (12) 51 (42) 58 (54)
20 (16)
1.21 (0.98)
0.30 (0.30)
2.5 (2.5)
21 (24)
42 (46)
0.57
Rates were calculated using an average July PAR photon density flux at the ELA of 46.1 E m-2 d-1 and MeHg concentrations of 0.81 (this study) and 0.02 ng L-1 (9) for the colored-water Lake 979 and clear-water Lake 239, respectively. All values were corrected for the attenuation of solar radiation by Teflon bottles (25) and Mylar/Lee film (Figure S2) where applicable, and uncorrected values are shown in parentheses. b The areal MeHg PD flux was calculated assuming an equal volume of water in each 5 cm section of the water column within the zone of MeHg PD activity; this assumption is valid where lake depth is always greater than depth of MeHg PD activity, such as in the pelagic zone. c The 10% depth is defined as the depth below which MeHg PD activity is less than 10% of the overall areal MeHg PD flux. a
behaves similarly to ambient MeHg at the scale required to quantify MeHg PD; however, this does not preclude the possibility that fine-scale mass-dependent and massindependent isotope fractionation does take place, as has been recently demonstrated for the aqueous photoreduction of MeHg (23). Therefore, our results suggest that the standard methodology of quantifying MeHg PD using the addition of a spike of MeHg, whether isotopically enriched or not, is appropriate, at least in cases where the added spike does not increase the ambient MeHg concentration by more than an order of magnitude. Comparison with Previously Published kpds. As neither the addition of a spike of isotopically enriched MeHg nor filtration had a significant effect on kpd, the data were pooled, and corrected for the decomposition of MeHg in the absence of light, to obtain a value of kpd for each level of the light treatment (Table 1). The pooled value of kpd-overall ((standard error) measured in lake water exposed to full solar radiation was (3.69 ( 0.07) × 10-3 m2 E-1 which is very similar to the value of 3.72 × 10-3 m2 E-1 previously published for samples incubated under full solar radiation on the surface of an Alaskan tundra lake (15). However, if our value of kpd-overall is expressed as 0.170 d-1 (3.69 × 10-3 m2 E-1 × 46.1 E m-2 d-1), it no longer compares favorably with the value of 0.230 d-1 reported in ref 15 as a result of the higher mean daily PAR photon flux density
(62 E m-2 d-1) measured in Alaska in July. This underlines the importance of reporting kpd in units that are consistent with a photochemical process (i.e., related to photon flux density or irradiance). In a previous study conducted in Lake 240 at the ELA in summer-time, and therefore under light conditions similar to our experiments, kpd averaged 0.19 d-1 (9), comparable to our value of 0.170 d-1. The similarity of kpd values obtained in lakes from different regions and with different water-chemistry characteristics suggests that the kinetics of MeHg PD are not influenced by site characteristics apart from light intensity and MeHg concentration, and therefore, could be universal in freshwater bodies, as was previously suggested (8, 15). Furthermore, it also appears that MeHg PD proceeds via indirect photolysis and thus requires the presence of a photosensitizing species such as nitrate or dissolved organic matter ((14), Lehnherr, unpublished data), although theoretical considerations suggest that direct photolysis of MeHg by UV-B is possible (24). All previous field measurements of kpd, including the value we present here, are in good agreement but were all obtained using Teflon bottles (8, 9, 15, 17). However, Teflon bottles, which are more affordable and easier to work with in a field setting than quartz bottles, for example, are not completely transparent to UV radiation (9, 17). Therefore, values of kpd obtained using Teflon bottles are likely VOL. 43, NO. 15, 2009 / ENVIRONMENTAL SCIENCE & TECHNOLOGY
9
5695
FIGURE 3. Depth-specific (A, C) and depth-integrated (B, D) rates of UV and visible light mediated MeHg PD in the water column of a colored-water lake (Lake 979; A, B) and a clear-water lake (Lake 239; C, D). Depth-specific rates were summed over the specified water column depth to obtain depth-integrated rates. Note the different scales. underestimates of the true kpd. The transmission of visible and UV radiation through a lengthwise cross-section of a FEP-Teflon bottle was measured using a Shimadzu UV2401PC spectrophotometer (Figure S2). Similar to recently published results (28), we found that 79% of visible, 51% of UV-A and 30% of UV-B energy was transmitted by the Teflon bottle, which is less than previously reported (99%, 82%, and 66% transmittance for visible, UV-A, and UV-B radiation, respectively) by measuring solar irradiance at the surface of a lake with the sensor placed inside a Teflon bottle (25). Because the emission spectrum of the radiation source in the spectrophotometer (deuterium and halogen lamps) is different from solar radiation at the earth’s surface, and because having the sensor inside of the Teflon bottle floating at the surface of a lake more closely replicates the conditions to which samples were exposed, the results reported in ref 25 are probably more relevant in this case. To constrain the upper value of kpd-overall, we corrected for the attenuation of UV and visible radiation by Teflon (and Mylar/Lee films) using the spectra shown in Figure S2, resulting in a kpd-overall value of 5.57 × 10-3 m2 E-1 compared to 3.69 × 10-3 m2 E-1 without any correction. However, using the transmittance of solar radiation through Teflon bottles reported in ref 25, a likely more accurate kpd-overall value of 4.41 × 10-3 m2 E-1 is obtained. To improve the accuracy of future measurements of kpd, it would be helpful to measure UV and visible light irradiance inside a Teflon bottle throughout the experiment, using either a spectroradiometer or a solar actinometer. Alternatively, reaction vessels with better transmittance of solar radiation, such as PFA-Teflon bags (29) could also be used instead of FEPTeflon bottles. Modeling MeHg PD at the Whole-Lake Level. To scale up experimental results to the whole-lake level, we 5696
9
ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 43, NO. 15, 2009
calculated a rate of MeHg PD at 5-cm depth-intervals according to eq 3, for both a colored-water lake (Lake 979) and a clear, oligotrophic lake (Lake 239) at the ELA. MeHg PD Rate (ng L-1 d-1) ) kpd (m2 E-1) × MeHg (ng L-1) × PFD (E m-2d-1)
(3)
MeHg PD resulting from UV radiation and visible light was evaluated separately, such that three rate values were obtained for each depth (i.e., a separate rate value for UV-B, UV-A, and visible light mediated processes). kpd-UVB for UV-B mediated MeHg PD (1.25 × 10-3 m2 E-1) was calculated as the difference between the “full spectrum” kpd (UV-B + UV-A + visible) and kpd obtained with UV-B blocked (UV-A + visible) (Tables 1 and 2). Similarly, kpd-UVA for UV-A mediated MeHg PD (2.16 × 10-3 m2 E-1) was calculated as the difference between the “UV-A + visible” and “visible only” kpds (Tables 1 and 2). Photon flux density (PFDz) at each depth (z) was calculated according to eq 4, with PFD0 ) 46.1 E m-2 d-1, corresponding to the mean daily PAR photon flux density measured over the duration of our experiments. PFDz ) PFD0 × e-kz
(4)
Light attenuation coefficients (k) for UV-B and UV-A were estimated from DOC concentrations measured in July 2006 in each lake (Susan Kasian, Fisheries and Ocean Canada, Winnipeg, personal communication) using the equations provided in refs 26 and 27. Attenuation coefficients for PAR were calculated from duplicate light profiles conducted in July 2006 (Susan Kasian, personal communication). Finally, MeHg concentrations used to model MeHg PD rates for Lake
979 and Lake 239 were 0.81 ng L-1 (this study) and 0.02 ng L-1 (9), respectively. The depth-specific rates calculated according to eq 3 can then be multiplied by the appropriate volume of water in each 5-cm “slice” of the water column, which depends on lake bathymetry, and summed to obtain depth-integrated whole-lake fluxes of MeHg PD (ng d-1). For the purpose of our modeling exercise, however, we calculated areal fluxes of MeHg PD (ng m-2 d-1) assuming an equal volume of water in each 5-cm section of the water column within the zone of MeHg PD activity. This assumption is valid for the pelagic zone of a lake, where lake depth is always greater than the depth of MeHg PD activity. Areal fluxes were then obtained by summing depth-specific rates from the surface down to the depth at which the modeled MeHg PD rate was 1% of the surface rate. This was done to exclude contributions to the areal MeHg PD flux from depth-specific rates which are not significantly greater than zero but which could contribute significantly to the sum value in a deep and/or colored lake, thus artificially inflating the modeled areal MeHg PD flux. Modeled depth-specific and depth-integrated rates of MeHg PD (Figure 3) clearly demonstrate that for a colored, high DOC lake such as Lake 979, MeHg PD is primarily driven by UV radiation and is therefore largely restricted to the top 30 cm of the water column. Indeed, UV mediated MeHg PD accounts for over three-quarters of the overall areal flux, with UV-A radiation the single-most important contributor (58%) (Table 2). By comparison, in a clear, oligotrophic, low DOC lake such as Lake 239, visible light mediated MeHg PD becomes relatively more important as depth increases, such that UV-A and visible light contribute equally toward the overall areal MeHg PD flux. Furthermore, in Lake 239, MeHg PD activity encompassed the top 2.5 m of the water column, compared to only 30 cm for Lake 979. Despite the shallower zone of MeHg PD activity resulting from greater attenuation of solar radiation, the areal MeHg PD flux was much larger in Lake 979 than in Lake 239 (Table 2), demonstrating the important effect of aqueous MeHg concentration on MeHg PD rates and thus areal flux. These results indicate that UV radiation is the most important driver of MeHg PD in freshwaters, in contrast to previous reports suggesting that the contribution from UV to whole-lake fluxes is negligible (15). If the values of light-specific kpds are corrected for the attenuation of UV and visible radiation through Teflon bottles (25) and Mylar/Lee films (Figure S2), we obtain kpds of 8.71 × 10-4 m2 E-1, 3.22 × 10-3 m2 E-1, and 3.13 × 10-4 m2 E-1 for the UV-B, UV-A, and visible light mediated reactions, respectively (Table 2). As a result of this correction, the areal MeHg PD flux increases by 26% and 24% in Lakes 979 and 239, respectively, although the relative contributions of the UV and visible light mediated processes toward the overall flux remain almost unchanged. Using corrected kpds, we find that UV mediated MeHg PD contributes 79% of the areal flux (up from 76%) in Lake 979, and 58% of the areal flux (up from 54%) in Lake 239. Applying the modeling method described above to an Alaskan tundra lake (Toolik Lake; (15)) and using the values published in that study for mean daily PAR, PAR attenuation coefficient, and MeHg concentration, we obtained an areal MeHg PD flux of 0.57 µg MeHg m-2 y-1. This value is significantly smaller than the 1.31 µg MeHg m-2 y-1 reported by these authors, despite the close agreement in the value of kpd measured at the surface of Toolik Lake with the value we report here (see above). The discrepancy in the flux value results from the different approaches used to model MeHg PD; we broke the flux down into UV and visible light mediated components, whereas these authors attributed the PD of MeHg entirely to visible light wavelengths, which we found to degrade MeHg at a much slower rate than UV
radiation (Figure 1, Table 1). Because visible light is attenuated much less rapidly than UV radiation, attributing MeHg PD solely to visible light has the effect of deepening the zone of MeHg PD activity, resulting in a larger areal MeHg PD flux. In an effort to reconcile these differing results, we suggest that the role of UV radiation may have been previously underestimated due to the strong positive correlation observed between PAR intensity and MeHg PD rate measured at depths of 0.75, 1.5, 3, and 6 m in the water column of Toolik Lake (15). However, the rates measured at depths of 0.75 and 1.5 m, which were used to build the relationship between incident PAR and MeHg PD rate, may have been influenced by non-negligible levels of UV radiation. We estimate that in Toolik Lake, UV mediated MeHg PD still accounts for 67% and 40% of the overall rate at depths of 0.75 and 1.5 m, respectively. We therefore propose that further measurements of MeHg PD rates throughout the water column, coupled with measurements of incident UV radiation and PAR at each depth, are warranted to clarify this matter.
Acknowledgments We are grateful to Linnea Mowat for valuable assistance in the field and Po Yee Chan for assistance during sample analysis. We thank Jennifer Graydon for providing PAR data. Finally, we thank Rolf Vinebrooke and Rich Moses for contributing their knowledge of statistics to the experimental design and data analysis. This research was funded through a Natural Science and Engineering Research Council (NSERC) Discovery Grant to V.S.L., and a NSERC Post-Graduate Scholarship and Alberta Ingenuity Studentship to I.L.
Supporting Information Available Analytical details and five additional figures. This material is available free of charge via the Internet at http:// pubs.acs.org.
Literature Cited (1) Mergler, D.; Anderson, H. A.; Hing Man Chan, L.; Mahaffey, K. R.; Murray, M.; Sadamoto, M.; Stern, A. H. Methylmercury exposure and health effects in humans: a worldwide concern. Ambio 2007, 36, 3–11. (2) Morel, F. M. M.; Krepiel, A. M. L.; Amyot, M. The chemical cycle and bioaccumulation of mercury. Annu. Rev. Ecol. Syst. 1998, 29, 543–566. (3) Fitzgerald, W. F.; Clarkson, T. W. Mercury and monomethylmercury - present and future concerns. Environ. Health Perspect. 1991, 96, 159–166. (4) Ullrich, S. M.; Tanton, T. W.; Abdrashitova, S. A. Mercury in the Aquatic Environment: A Review of Factors Affecting Methylation. Crit. Rev. Environ. Sci. Technol. 2001, 31, 241–293. (5) Oremland, R. S.; Culbertson, C. W.; Winfrey, M. R. Methylmercury decomposition in sediments and bacterial cultures involvement of methanogens and sulfate reducers in oxidative demethylation. Appl. Environ. Microbiol. 1991, 57, 130–137. (6) Marvin-DiPasquale, M.; Agee, J.; McGowan, C.; Oremland, R. S.; Thomas, M.; Krabbenhoft, D.; Gilmour, C. C. Methyl-mercury degradation pathways: a comparison among three mercury impacted ecosystems. Environ. Sci. Technol. 2000, 34, 4908– 4916. (7) Schaefer, J. K.; Yagi, J.; Reinfelder, J. R.; Cardona, T.; Ellickson, K. M.; Tel-Or, S.; Barkay, T. Role of the bacterial organomercury lyase (MerB) in controlling methylmercury accumulation in mercury-contaminated natural waters. Environ. Sci. Technol. 2004, 38, 4304–4311. (8) Sellers, P.; Kelly, C. A.; Rudd, J. W. M.; MacHutchon, A. R. Photodegredation of methylmercury in lakes. Nature 1996, 380, 694–697. (9) Sellers, P.; Kelly, C. A.; Rudd, J. W. M. Fluxes of methylmercury to the water column of a drainage lake: the relative importance of internal and external sources. Limnol. Oceanogr. 2001, 46, 632–631. (10) Hammerschmidt, C. R.; Fitzgerald, W. F.; Lamborg, C. H.; Balcom, P. H.; Mao Tseng, C. Biogeochemical cycling of methylmercury in lakes and tundra watersheds of arctic Alaska. Environ. Sci. Technol. 2006, 40, 1204–1211. VOL. 43, NO. 15, 2009 / ENVIRONMENTAL SCIENCE & TECHNOLOGY
9
5697
(11) Inoko, M. Studies on the photochemical decomposition of organomercurials - methylmercury (II) chloride. Environ. Pollut. (Series B) 1981, 2, 3–10. (12) Suda, I.; Suda, M.; Hirayama, K. Degradation of methyl and ethyl mercury by singlet oxygen generated from sea water exposed to sunlight or ultraviolet light. Arch. Toxicol. 1993, 67, 365–368. (13) Gårdfeldt, K.; Sommar, J.; Stro¨mberg, D.; Feng, X. Oxidation of atomic mercury by hydroxyl radicals and photoinduced decomposition of methylmercury in the aqueous phase. Atmos. Environ. 2001, 35, 3039–3047. (14) Chen, J.; Pehkonen, S. O.; Lin, C.-J. Degradation of monomethylmercury chloride by hydroxyl radicals in simulated natural waters. Water Res. 2003, 37, 2496–2504. (15) Hammerschmidt, C. R.; Fitzgerald, W. F. Photodecomposition of methylmercury in an arctic Alaskan lake. Environ. Sci. Technol. 2006, 40, 1212–1216. (16) Sellers, P. Sediment Flux and Photodegradation of Methylmercury in Two Boreal Drainage Lakes (Ph.D. Thesis). University of Manitoba: Winnipeg, 1997. (17) Krabbenhoft, D. P.; Olson, M. L.; deWild, J. F.; Clow, D. W.; Striegl, R. G.; Dornblaser, M. M.; Vanmetre, P. Mercury loading and methylmercury production in high-altitude lakes from the western United States. Water, Air, Soil Pollut. Focus 2002, 2, 233–249. (18) St. Louis, V. L.; Rudd, J. W. M.; Kelly, C. A.; Bodaly, R. A.; Paterson, M. J.; Beaty, K. G.; Hesslein, R. H.; Heyes, A.; Majewski, A. R. The rise and fall of mercury methylation in an experimental reservoir. Environ. Sci. Technol. 2004, 38, 1348–1358. (19) St. Louis, V. L.; Rudd, J. W. M.; Kelly, C. A.; Beaty, K. G.; Bloom, N. S.; Flett, R. J. Importance of wetlands as sources of methyl mercury to boreal forest ecosystems. Can. J. Fish. Aquat. Sci. 1994, 51, 1065–1076. (20) Hintelmann, H. Comparison of different extraction techniques used for methylmercury analysis with respect to accidental formation of methylmercury during sample preparation. Chemosphere 1999, 39, 1093–1105.
5698
9
ENVIRONMENTAL SCIENCE & TECHNOLOGY / VOL. 43, NO. 15, 2009
(21) Hintelmann, H.; Evans, R. D.; Villeneuve, J. Y. Measurements of mercury methylation in sediments by using enriched stable mercury isotopes combined with methylmercury determination by gas-chromatography-inductively coupled plasma mass spectroscopy. J. Anal. At. Spectrom. 1995, 10, 619–624. (22) Hintelmann, H.; Evans, R. D. Application of stable-isotopes in environmental tracer studies - measurements of monomethylmercury (CH3Hg+) by isotope dilution ICP-MS and detection of species transformation. Fresnius J. Anal. Chem. 1997, 358, 378–385. (23) Bergquist, B. A.; Blum, J. D. Mass-dependent and -independent fractionation of Hg isotopes by photoreduction in aquatic systems. Science 2007, 318, 417–420. (24) Tossel, J. A. Theoretical study of the photodecomposition of methyl Hg complexes. J. Phys. Chem. A 1998, 102, 3587–3591. (25) Amyot, M.; Lean, D. R. S.; Mierle, G. Photochemical formation of volatile mercuy in high Arctic lakes. Environ. Toxicol. Chem. 1997, 16, 2054–2063. (26) Donahue, W. F.; Turner, M. A.; Findlay, D. L.; Leavitt, P. R. The role of solar radiation in structuring the shallow benthic communities of boreal forest lakes. Limnol. Oceanogr. 2003, 48, 31–47. (27) Scully, N. M.; Lean, D. R. S. The attenuation of ultraviolet radiation in temperate lakes. Arch. Hydrobiol. Beih. 1994, 43, 135–144. (28) Black, F. J.; Conaway, C. H.; Flega, A. R. Stability of dimethyl mercury in seawater and its conversion to monomethyl mercury. Enivron. Sci. Technol. 2009, 43 (11), 4056–4062. (29) Whalin, L. M.; Mason, R. P. A new method for the investigation of mercury redox chemistry in natural waters utilizing deflatable Teflon bags and addition of isotopically labeled mercury. Anal. Chim. Acta 2006, 558, 211–221.
ES9002923