Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Experimental evidence for recovery of mercury-contaminated fish populations

## Abstract

Anthropogenic releases of mercury (Hg)1,2,3 are a human health issue4 because the potent toxicant methylmercury (MeHg), formed primarily by microbial methylation of inorganic Hg in aquatic ecosystems, bioaccumulates to high concentrations in fish consumed by humans5,6. Predicting the efficacy of Hg pollution controls on fish MeHg concentrations is complex because many factors influence the production and bioaccumulation of MeHg7,8,9. Here we conducted a 15-year whole-ecosystem, single-factor experiment to determine the magnitude and timing of reductions in fish MeHg concentrations following reductions in Hg additions to a boreal lake and its watershed. During the seven-year addition phase, we applied enriched Hg isotopes to increase local Hg wet deposition rates fivefold. The Hg isotopes became increasingly incorporated into the food web as MeHg, predominantly from additions to the lake because most of those in the watershed remained there. Thereafter, isotopic additions were stopped, resulting in an approximately 100% reduction in Hg loading to the lake. The concentration of labelled MeHg quickly decreased by up to 91% in lower trophic level organisms, initiating rapid decreases of 38–76% of MeHg concentration in large-bodied fish populations in eight years. Although Hg loading from watersheds may not decline in step with lowering deposition rates, this experiment clearly demonstrates that any reduction in Hg loadings to lakes, whether from direct deposition or runoff, will have immediate benefits to fish consumers.

## Main

The Minamata Convention on Mercury is an international treaty that aims to protect human health and the environment from adverse effects of MeHg by controlling Hg emissions, which should then decrease deposition and loading of anthropogenic Hg to aquatic environments10. Yet there is little direct evidence for how quickly fish MeHg concentrations will decline following reductions in current rates of Hg loading owing, in part, to a range of ecological factors that can influence both the microbial production and the bioaccumulation of MeHg in aquatic food webs9,11. Further complicating this relationship are human activities such as commercial fishing, introduction of exotic species and enhanced nutrient additions that trigger large-scale trophic disruptions12, which can in turn substantially alter fish tissue MeHg concentrations13,14,15, because fish acquire most of their MeHg through their diet16. Changing climatic conditions can also influence MeHg production8, as well as restructure food webs, alter the dominant pathways of energy flow and cause size-dependent changes in fish growth rates that shift population size structure14,17,18. In addition, until now there has been no way to evaluate the relative contribution of newly deposited Hg to contemporary MeHg production. Consequently, it is exceedingly difficult to unambiguously assess the recovery of contaminated fish populations due specifically to Hg control measures9.

The METAALICUS watershed is in an undisturbed remote region of Canada, such that our experimental addition rate increased wet Hg deposition approximately fivefold (from approximately 3.6 to 19 µg m−2 yr−1), to levels similar to more polluted regions of the world21. Most of the Hg added to the wetland and upland areas of the watershed either remained bound to vegetation and soils or evaded back to the atmosphere22,23. The wetland spike was below the detection level in all fish species. The Hg applied to the upland catchment accounted for only a small fraction (less than 1%) of all Hg in runoff to the lake19 and consequently contributed little (less than 2%) to the changes in MeHg concentrations of fish populations throughout the study (Extended Data Fig. 1). Hence, after six years of increased additions to the watershed (Fig. 1a), the observed increases in fish MeHg were due almost entirely to Hg added directly to the lake surface. This did not appear to be caused by preferential methylation of lake spike Hg. The evidence for this is that the seasonal production of lake spike MeHg, all of which had to be formed within the lake because it was added as inorganic Hg, varied synchronously with ambient MeHg in the lake and in biota19. This finding also indicates that in this headwater lake ambient MeHg is mainly derived from in-lake methylation of inorganic ambient Hg, most of which came from the upland catchment.

Delivery of lake spike Hg to the sediments and anoxic bottom waters, which are the dominant sites of methylation in the study lake19,24, resulted in formation of spike MeHg in these compartments, from where spike MeHg also migrated to surface waters. Lake spike MeHg rapidly accumulated in all lake biota (Fig. 1b), with concentrations in fish muscle increasing with continued loading for all species (Extended Data Table 1; linear regression, P < 0.05), apart from young-of-year (YOY) yellow perch, which showed high inter-annual variability in spike MeHg concentrations after an initial increase (Fig. 1b). By contrast, ambient MeHg concentrations in all fish species did not show any consistent trends during the addition phase (Fig. 2a–c), nor did they in a nearby reference lake (Extended Data Table 1; P > 0.05). Steady ambient MeHg concentrations in fish through time are indicative of relatively stable watershed inputs of Hg, which is the main source of ambient inorganic Hg for methylation in both the experimental and reference lakes19,25.

A critical question for the addition period was how much higher were MeHg concentrations than they would have been in the absence of the experimentally increased loading of Hg. The addition of lake spike Hg was roughly equivalent to all ambient Hg inputs (runoff plus direct deposition) to the lake, resulting in a doubling, or about 100% increase, in Hg loading to the lake19. In response to seven years of experimental additions to the lake, per cent increases in lake spike MeHg concentrations were highest in water (60%) and least in the upper 2 cm of sediments (30%) where large stores of ambient Hg existed (Fig. 1c). The response of food web organisms was intermediate to that of water and sediments, such that spike Hg additions to the lake raised MeHg concentrations by 45–57% in invertebrates and forage fishes and by more than 40% for large-bodied fish species (Fig. 1c).

Temporal patterns of biomagnification for spike MeHg relative to ambient MeHg inform how quickly the different fish species came into equilibrium with their respective prey. For small-bodied yellow perch (1 year of age) feeding on zooplankton, it took three years before the biomagnification of lake spike resembled that of ambient MeHg (Fig. 2d), and a further two years for both the apex predator, northern pike, feeding on forage fishes (Fig. 2e), and lake whitefish feeding on Chaoborus (Fig. 2f). Relative to planktivorous yellow perch, final addition phase concentrations of spike MeHg were slightly higher for benthivorous lake whitefish (1.2×) and further increased for piscivorous northern pike (3.9×), similar to ambient MeHg (whitefish (1.5×) and pike (3.9×); relative to perch; Fig. 2a–c) and consistent with expectations of contaminant biomagnification among trophic guilds5,20. These findings imply that the key in-lake processes leading to the formation and trophic transfer of MeHg to the different fish populations became comparable for spike Hg and ambient Hg during the addition phase.

To then directly test the hypothesis that MeHg concentrations in fishes would decline following reductions in Hg loading to the lake, we ceased all experimental additions of enriched Hg isotopes (Fig. 1a). This resulted in a 100% reduction in loading of lake spike. Average concentrations of lake spike MeHg in fish populations rapidly declined (within less than ten years) in concert with the decline in the availability of spike MeHg through dietary and waterborne pathways (Fig. 1b). Within the first 3 years, the relative amount of lake spike MeHg declined by 81% in water, 35% in sediments, 66% in zooplankton and 67% in Chaoborus (Fig. 1c), leading to marked reductions (85–91%) in the concentration of spike MeHg in forage fish species by the end of the recovery phase (Fig. 1b). Eight years after addition, lake spike Hg contributed just a small fraction (approximately 6%) to MeHg concentrations in forage fishes and invertebrate prey (Fig. 1c). The more rapid decline in per cent spike MeHg in water compared with sediments, even in this relatively long water residence time lake (about 6 years), emphasizes that the magnitude and timing of responses by fish to Hg loading reductions could be influenced by their relative reliance upon pelagic versus benthic dietary pathways26.

The notably fast response of the lower food web to the cessation of lake spike Hg loadings initiated rapid recovery of large-bodied fish species. Within 8 years, lake spike MeHg concentrations declined by 76% in the northern pike population and by 38% in the lake whitefish population (Fig. 1b). During the recovery phase, spike MeHg concentrations for these large-bodied species initially increased for both populations before showing steady declines (Fig. 2b, c). The rate of decline in spike MeHg for northern pike, however, was roughly twice that of lake whitefish (Fig. 1c).

Differences in the lifespan of the fish populations had a key role in the rates of recovery following the reduction of Hg loadings to the lake. Lake whitefish were much older (median age = 17 years versus 3 years for pike) and larger (Extended Data Tables 2, 3) than northern pike, and more individuals in that population would have lived through some or all of the addition and recovery phases of the experiment. Lake whitefish had the coldest thermal preferences and greatest association with benthic habitats of any fish population, which probably also contributed to their delayed recovery.

Boreal fishes are known to eliminate MeHg very slowly once accumulated27. To further explain the recovery of the apex predator population, we tracked changes in the body burdens of spike MeHg in individual northern pike over time while also monitoring the population as a whole. As expected, individual responses were variable, but lake spike MeHg burdens in northern pike mostly increased during the early recovery phase with overall little to no loss of the spike MeHg 6–8 years after cessation of spike additions (Fig. 3). These findings parallel those observed for lake spike MeHg in individual northern pike moved from the study lake to a nearby reference lake28 and underscore how the prolonged retention of MeHg in fish muscle tissue can delay recovery of some fisheries7,29. Thus, it was the annual recruitment of new fish with low MeHg concentrations into the population, along with the loss of older fish (as evidenced by a stable population size structure; Extended Data Fig. 2), that enabled the swift recovery of the population from Hg contamination as a whole. Consequently, average burdens of spike MeHg in the northern pike population were reduced by 50% in less than five years, in spite of the efficient retention of spike MeHg by some older fish (Fig. 3, Extended Data Fig. 2).

Differentiating the relative importance of present day inputs versus previously deposited Hg to overall fish MeHg concentrations is a key uncertainty when predicting the efficacy of Hg pollution reduction29,30. Here we demonstrate that within a few years of abatement, the experimental Hg added previously to the lake was no longer an important source of MeHg to the lower food web or to forage fishes (Fig. 1c). Long-lived fish species of subsistence, commercial and recreational importance lagged behind their prey, but the contribution of recently deposited Hg to fish MeHg steadily diminished for these populations as well. There was a similar, rapid response for the upland spike when loading ceased, even though only a small amount appeared in the fish during the experiment (Extended Data Fig. 1). The small contribution of the terrestrial spike to fish MeHg supports our former conclusion19 that lakes with large watersheds will respond more slowly to changes in atmospheric deposition.

The most important outcome of this whole-ecosystem experiment is the demonstration that a decrease in a single factor (Hg loading to the lake) has a clear and timely effect on average MeHg concentrations in fish populations, even for long-lived species that eliminate MeHg slowly. The spike MeHg data show that fish populations will respond quickly to any change in loading rates—whether from direct deposition to the lake (Fig. 1) or runoff (Extended Data Fig. 1). Decreases in loading to the lake from these two sources will follow different time courses in response to lower atmospheric deposition19. However, as these two loads decrease, the fish populations in the receiving lake will soon afterwards have lower MeHg than they would have if nothing were done, thereby reducing human exposure.

## Methods

### Mercury additions to the study catchment

METAALICUS was conducted on the Lake 658 catchment at the Experimental Lakes Area (ELA; now IISD-ELA), a remote area in the Precambrian Shield of northwestern Ontario, Canada (49° 43′ 95″ N, 93° 44′ 20″ W) set aside for whole-ecosystem research31. The Lake 658 catchment includes upland (41.2 ha), wetland (1.7 ha) and lake surface (8.4 ha) areas. Lake 658 is a double basin (13 m depth), circumneutral, headwater lake, with a fish community consisting of forage (yellow perch (P. flavescens) and blacknose shiner (Notropis heterolepis)), benthivorous (lake whitefish (C. clupeaformis) and white sucker (Catostomus commersonii)), and piscivorous (northern pike (E. lucius)) fishes. The lake is closed to fishing.

Hg addition methods used in METAALICUS have been described in detail elsewhere19,32,33. In brief, three Hg spikes, each enriched with a different stable Hg isotope, were applied separately to the lake surface, upland and wetland areas. Upland and wetland spikes were applied once per year (when possible; Fig. 1a) by fixed-wing aircraft (Cessna 188 AGtruck). Mercury spikes (as HgNO3) were diluted in acidified water (pH 4) in a 500 l fiberglass tank and sprayed with a stainless-steel boom on upland (approximately 79.9% 200Hg) and wetland (approximately 90.1% 198Hg) areas. Spraying was completed during or immediately before a rain event, with wind speeds less than 15 km h−1 to minimize drift of spike Hg outside of target areas. Aerial spraying of upland and wetland areas left a 20-m buffer to the shoreline, which was sprayed by hand with a gas-powered pump and fire hose to within about 5 m of the lake32. Average net application rates of isotopically labelled Hg to the upland and wetland areas were 18.5 μg m−2 yr−1 and 17.8 μg m−2 yr−1, respectively.

The average net application rate for lake spike Hg was 22.0 μg m−2 yr−1. For each lake addition, inorganic Hg enriched with approximately 89.7% 202Hg was added as HgNO3 from four 20-l carboys filled with acidified lake water (pH 4). Nine lake additions were conducted bi-weekly at dusk over an 18-week (wk) period during the open-water season of each year (2001–2007) by injecting at 70-cm depth into the propeller wash of trolling electric motors of two boats crisscrossing each basin of the lake32,33. It was previously demonstrated with 14C additions to an ELA lake that this approach evenly distributed spike added in the evening by the next morning34.

We did not attempt to simulate Hg in rainfall for isotopic lake additions because it is impossible to simulate natural rainfall concentrations (about 10 ng l−1) in the 20-l carboys used for additions. Instead, our starting point for the experiment was to ensure that the spike was behaving as closely as possible to ambient surface water Hg very soon after it entered the lake. Several factors support this assertion. By the next morning each spike addition had increased epilimnetic Hg concentrations by only 1 ng l−1 202Hg. Average ambient concentrations were 2 ng l−1. Thus, while the Hg concentrations in the carboys were high (2.6 mg l−1), the receiving waters were soon at trace levels. Furthermore, we investigated if the additions altered the degree of bioavailability or photoreactivity of Hg(ii) in the receiving surface water. We examined the bioavailability of spike Hg(ii) as compared to ambient Hg in the lake itself using a genetically engineered bioreporter bacterium35. On seven occasions, epilimnetic samples were collected on the day before and within 12 h of spike additions. The spike was added to the lake as Hg(NO3)2, which is bioavailable to the bioreporter bacterium (detection limit = 0.1 ng Hg(ii) l−1), but we never saw bioavailable ambient or spike Hg(ii) in the lake, presumably because it was quickly bound to dissolved organic carbon (DOC). This indicates that, in terms of bioavailability, the spike Hg was behaving like ambient Hg soon after additions. Photoreactivity in the surface water was examined on seven occasions, by measuring the % of total Hg(ii) that was dissolved gaseous Hg for spike and ambient Hg, either 24 h or 48 h after the lake was spiked36. There was no significant difference (paired t-test, P > 0.05), demonstrating that by then the lake spike was behaving in the same way as ambient Hg during gaseous Hg production.

### Lake, food web and fish sampling

Water samples were collected from May to October every four weeks at the deepest point of Lake 658. Water was pumped from six depths through acid-cleaned Teflon tubing into acid-cleaned Teflon or glass bottles. Water samples were filtered in-line using pre-ashed quartz fibre filters (Whatman GFQ, 0.7 µm). Subsequently, Hg species were measured in the filtered water samples (dissolved Hg and MeHg) and in particles collected on the quartz fibre filter (particulate Hg and MeHg).

From 2001 to 2012, Lake 658 sediments were sampled at 4 fixed sites up to 5 times per year. Sampling frequency was highest in 2001, with monthly sampling from May to September, and declined over the course of the study. Fixed sites were located at depths of 0.5, 2, 3 and 7 m. A sediment survey of up to 12 additional sites was also conducted once or twice each year. Survey sites were selected to represent the full range of water depths in both basins. Cores were collected by hand by divers, or by subsampling sediments collected using a small box corer. Cores were capped and returned to the field station for processing within a few hours. For each site, three separate cores were sectioned and composited in zipper lock bags for a 0- to 2-cm depth sampling horizon, and then frozen at −20 °C.

Bulk zooplankton and Chaoborus samples were collected from Lake 658 for MeHg analysis. Zooplankton were collected during the day from May to October (bi-weekly: 2001–2007; monthly: 2008–2015). A plankton net (150 μm, 0.5 m diameter) was towed vertically through the water column from 1 m above the lake bottom at the deepest point to the surface of the lake. Samples were frozen in plastic Whirl-Pak bags after removal of any Chaoborus using acid-washed tweezers. Dominant zooplankton taxa in Lake 658 included calanoid copepods (Diaptomus oregonensis) and Cladocera (Holopedium glacialis, Daphnia pulicaria and Daphnia mendotae). Chaoborus samples were collected monthly in the same manner at least 1 h after sunset. After collection, Chaoborus were picked from the sample using forceps and frozen in Whirl-Pak bags. Chaoborus were not separated by species for MeHg analyses, but both C. flavicans and C. punctipennis occur in the lake. Profundal chironomids were sampled at the deepest part of the lake using a standard Ekman grab sampler. Grab material was washed using water from a nearby lake and individual chironomids were picked by hand.

All work with vertebrate animals was approved by Animal Care Committees (ACC) through the Canadian Council on Animal Care (Freshwater Institute ACC for Fisheries and Oceans Canada, 2001–2013; University of Manitoba ACC for IISD-ELA, 2014–2015). Licenses to Collect Fish for Scientific Purposes were granted annually by the Ontario Ministry of Natural Resources and Forestry. Prior to any Hg additions, a small-mesh fence was installed at the outlet of Lake 658 to the downstream lake to prevent movement of fish between lakes. Sampling for determination of MeHg concentrations (measured as total mercury (THg), see below) occurred each autumn (August–October; that is, the end of the growing season in north temperate lakes) for all fish species in Lake 658, and for northern pike and yellow perch in nearby reference Lake 240 (Extended Data Tables 2, 3). Fish collections occurred randomly throughout the lakes. Forage fish (YOY and 1+ yellow perch, and blacknose shiner) were captured using small mesh gillnets (6–10 mm) set for <20 min, seine nets, and hoop nets. A small number of fish (up to n = 20) of each species and age class (determined by visual inspection) were euthanized immediately following capture in an overdose bath of 0.25 g l−1 tricaine methanesulfonate (TMS; Syndel Laboratories). After transport to the field station, fish were measured for fork length (FL; in mm) and mass (to 0.1 g), then immediately frozen (at −20 °C) in individual WhirlPak bags. A year class failure of yellow perch resulted in a single YOY collected in 2008 (data not presented) and no age 1+ fish in 2009 (Extended Data Table 3).

Large-bodied fish were captured by angling and multi-mesh gill nets (2.5–11.4 cm mesh) set for 20–30 min. Upon capture, each fish was anaesthetized with 0.06 g l−1 TMS, measured for FL (mm), weighed (to 1 g), tagged (Passive Integrated Transponder; Biomark), and a small biopsy of dorsal muscle (0.091 ± 0.002 g wet weight (mean ± s.e.m)) was collected using a dermal punch37. Only fish large enough for the biopsy procedure were sampled, such that our analyses include very few juveniles (pike: 317–850 mm FL; whitefish: 344–874 mm FL). Muscle samples were inserted into 0.6-ml polypropylene vials (Rose Scientific), immediately put on ice, and frozen within 4 h (−20 °C). This non-lethal method permitted repeated sampling of individual fish over time28,37. The first ray of either the pectoral or pelvic fin was collected for aging purposes upon first capture. Fish recovered from anaesthesia in a tub of fresh lake water (~15 min) before being released back into the lakes. From 2001–2015, we collected 690 biopsy muscle samples from 390 fish (238 northern pike, 114 lake whitefish and 38 white sucker) in Lake 658; 149 fish (90 northern pike, 38 lake whitefish and 21 white sucker) were biopsied more than once (2 to 6 per individual). Because of consistently low annual catches of white sucker (<10 individuals) across sampling years, we have excluded them from our analyses, but note here that their patterns of lake spike MeHg accumulation and recovery were similar to those of lake whitefish. We were unable to sex most fish because they were either immature or captured outside of their spawning season.

### Sample processing and analytical methods

Detailed methods on sample preparation and MeHg or THg analysis, as well as interlaboratory calibrations, have been reported elsewhere for the METAALICUS project19,38,39. In brief, MeHg was distilled from water samples and from sediment using atmospheric pressure water vapour distillation and measured after aqueous phase ethylation using sodium tetraethylborate (NaBEt4). Volatile Hg species were purged and trapped onto Tenax and MeHg was measured after thermodesorption and GC separation using inductively coupled plasma mass spectrometry (ICP-MS) detection (Micromass Platform or Perkin-Elmer Elan DRC II, respectively)39 and quantification by species specific isotope dilution mass spectrometry. The MeHg isotope dilution standards were synthesized and calibrated in-house. Isotope-dilution spikes were added prior to distillation, and MeHg external standards were routinely calibrated against degradation by measuring the standard against inorganic Hg before and after BrCl digestion. The QC strategy include the regular analysis of blanks, laboratory duplicates and certified reference materials (CRMs) IAEA 405 (International Atomic Energy Agency, Vienna, Austria) and NIST 1566b (National Institute of Standards and Technology, Gaithersburg, Maryland) for MeHg. No CRMs are commercially available for MeHg in water.

All biota samples were handled using clean techniques with Teflon or stainless steel tools cleaned with 95% ethanol19,38. Zooplankton and Chaoborus were freeze dried, ground with an acid-washed mortar and pestle, subsampled, and weighed to the nearest 0.00001 g. For determination of MeHg concentrations (ambient, lake spike, upland spike and wetland spike) in invertebrate samples, MeHg was solubilized by treatment with a solution of KOH in ethanol (20 % w/v), ethylated by additions of NaBEt4, and the resulting volatile Hg species were purged and trapped on carbotrap39. Samples were thermally desorbed and separated by gas chromatography before quantification by ICP-MS as above39. Samples of CRMs (TORT2 (2001–2013), IAEA452 (2014–2015); National Research Council of Canada, Ottawa, Ontario) were subjected to the same procedures; measured MeHg concentrations in the reference materials were not statistically different from certified values (P > 0.05).

Prey fish were kept frozen to maintain consistent wet weights. Approximately 0.2 g of skinless dorsal muscle was removed from each fish, weighed (0.0001 g), and placed in an acid-washed glass vial with a Teflon-lined cap (National Scientific Company). Muscle biopsy samples were weighed to the nearest 0.00001 g (Sartorius BP211D, Data Weighing Systems) before and after freeze-drying (Lyph-lock 12-l freeze dry system Model 77545, Labconco) to obtain wet and dry sample masses, and dry weight proportion28. Fish samples were analysed for THg, which is the sum of organic and inorganic Hg. Because we had previously determined that >90% of the Hg in muscle tissue from yellow perch in Lake 658 is MeHg40,41, here we report fish mercury data as MeHg.

THg concentrations (ambient, lake spike, upland spike and wetland spike) in fish muscle samples were quantified by ICP-MS39. Samples were digested with HNO3/H2SO4 (7:3 v/v) and heated at 80 °C until brown NOx gases no longer formed. The THg in sample digests was reduced by SnCl2 to Hg0 which was then quantified by ICP-MS (Thermo-Finnigan Element2) using a continuous flow cold vapour generation technique41. To correct for procedural recoveries, all samples were spiked with 201HgCl2 prior to sample analysis. Samples of CRMs (DORM2 (2001–2011), DORM3 (2012–2013), DORM4 (2014–2015); National Research Council of Canada) were submitted to the same procedures; measured THg concentrations in the reference materials were not statistically different from certified values (P > 0.05). Detection limit for each of the spikes was 0.5% of ambient Hg.

### Calculations and statistical methods

Analyses were completed with Statistica (6.1, Statsoft) and Sigmaplot (11.0, Systat Software). We present wet weight (w.w.) MeHg concentrations for all samples, except sediments which are dry weight (d.w.) concentrations. For zooplankton, Chaoborus, and profundal chironomids, d.w. MeHg concentrations were multiplied by a standard proportion (0.15) to yield w.w. concentrations for each sample42. The resulting w.w. concentrations were averaged over each open water season to determine annual means. For fish muscle biopsies, d.w. MeHg concentrations were multiplied by individual d.w. proportions to yield w.w. MeHg concentrations for each sample. To avoid any size-related biases, we calculated standardized annual MeHg concentrations (ambient and lake spike) for northern pike and lake whitefish by determining best-fit relationships between FL and MeHg concentrations for each year (quadratic polynomial, except for a linear fit for lake whitefish in 2004), and using the resulting regression equations to estimate MeHg concentrations at a standard FL43 (the mean FL of all fish sampled for each species: northern pike, 475 mm; lake whitefish, 530 mm). Square root transformation of raw northern pike data was required to satisfy assumptions of normality and homoscedasticity prior to standardization. The resulting data represent standardized concentrations of lake spike and ambient MeHg for each species each year.

We used the ratio of lake spike and ambient Hg in each sample as a measure of the amount by which Hg concentrations were changed with the addition of isotopically enriched Hg:

$${\rm{P}}{\rm{e}}{\rm{r}}{\rm{c}}{\rm{e}}{\rm{n}}{\rm{t}}\,{\rm{i}}{\rm{n}}{\rm{c}}{\rm{r}}{\rm{e}}{\rm{a}}{\rm{s}}{\rm{e}}={[{\rm{l}}{\rm{a}}{\rm{k}}{\rm{e}}{\rm{s}}{\rm{p}}{\rm{i}}{\rm{k}}{\rm{e}}{\rm{H}}{\rm{g}}]}_{i}/{[{\rm{a}}{\rm{m}}{\rm{b}}{\rm{i}}{\rm{e}}{\rm{n}}{\rm{t}}{\rm{H}}{\rm{g}}]}_{i}\times 100$$
(1)

where [lake spike Hg]i is the concentration of lake spike MeHg in sample i, and [ambient Hg]i is the concentration of ambient MeHg in sample i. For northern pike and lake whitefish, we calculated the mean annual relative increase from all individuals (not the size-standardized concentration data).

Biomagnification factors (BMF) were calculated to describe differences in Hg concentrations between predator and prey5:

$${\rm{BMF}}={\log }_{10}({[{\rm{MeHg}}]}_{{\rm{p}}{\rm{r}}{\rm{e}}{\rm{d}}{\rm{a}}{\rm{t}}{\rm{o}}{\rm{r}}}/{[{\rm{MeHg}}]}_{{\rm{p}}{\rm{r}}{\rm{e}}{\rm{y}}})$$
(2)

where [MeHg]predator is the mean (forage fish) or standardized (large-bodied fish) concentration of MeHg in the predator (ng g−1 w.w.) and [MeHg]prey is the mean concentration of MeHg in the prey (ng g−1 w.w.). MeHg concentration of prey items were averaged from samples collected throughout the open-water season immediately prior to autumn sampling of fish species to represent an integrated exposure for calculation of BMF. We used a dominant prey item to represent the diet of each fish species. For age 1+ yellow perch, northern pike, and lake whitefish, dominant prey items were zooplankton, forage fishes (YOY and 1+ yellow perch, and blacknose shiner) and Chaoborus, respectively.

To assess loss of lake spike MeHg by northern pike during the recovery period (2008–2015), we calculated28 whole body burdens (in μg) of lake spike MeHg for the standardized population and for individuals that had been sampled in autumn 2007 (t0 is the final time spike Hg was added to the lake) and again in at least one subsequent year during annual autumn sampling (n = 16 fish, of which 1–9 individuals were recaptured annually from 2008–2015). This calculation of MeHg burden is a relative measure of whole fish Hg content because MeHg is higher in muscle tissue than in other tissue types28,40. For the standardized population data, we used best-fit relationships between FL (in mm) and body weight (in g; quadratic polynomial) to determine body weight at the standard FL. We multiplied this body weight by standard ambient and spike MeHg concentrations (in ng g−1 w.w.) in muscle tissue for each year to determine body burdens over time (in ng). For individual fish, we multiplied spike MeHg concentration (in ng g−1 w.w.) by body weight (in g) to yield individual body burdens (in ng). To account for differences among individuals and between individuals and the population, we normalized the data to examine the mean proportion of original (t0) lake spike MeHg burden present in northern pike each year of the recovery period (2008–2015).

$${\rm{change}}\,{\rm{in}}\,{\rm{burden}}\,{\rm{from}}\,{t}_{0}={{\rm{burden}}}_{{\rm{tx}}}/{{\rm{burden}}}_{{\rm{t}}0}$$
(3)

We used a best fit regression (exponential decay, beginning in the second year of recovery) to estimate the half-life (50% of original burden) of lake spike MeHg for the population.

Northern pike and lake whitefish ages were determined by cleithra and otoliths, respectively, if mortality had occurred, but most ages were quantified using fin rays collected from live fish44 (K. H. Mills, DFO or North/South Consultants). Northern pike of the sizes selected for biopsy sampling had a median age of 3 years (range: 2–12 years; n = 305); the median age of lake whitefish was 17 years (range: 3–38 years; n = 86).

### Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this paper.

## Data availability

Datasets generated in this study are available at https://doi.org/10.5061/dryad.nzs7h44sfSource data are provided with this paper.

## References

1. Streets, D. G. et al. Five hundred years of anthropogenic mercury: spatial and temporal release profiles. Environ. Res. Lett. 14, 084044 (2019).

2. Fitzgerald, W. F., Engstrom, D. R., Mason, R. P. & Nater, E. A. The case for atmospheric mercury contamination in remote areas. Environ. Sci. Technol. 32, 1–7 (1998).

3. Driscoll, C. T., Mason, R. P., Chan, H. M., Jacob, D. J. & Pirrone, N. Mercury as a global pollutant: sources, pathways, and effect. Environ. Sci. Technol. 47, 4967–4983 (2013).

4. Mergler, D. et al. Methylmercury exposure and health effects in humans: a worldwide concern. Ambio 36, 3–11 (2007).

5. Lavoie, R. A., Jardine, T. D., Chumchal, M. M., Kidd, K. A. & Campbell, L. M. Biomagnification of mercury in aquatic food webs: a worldwide meta-analysis. Environ. Sci. Technol. 47, 13385–13394 (2013).

6. Sunderland, E. M., Li, M. & Bullard, K. Decadal changes in the edible supply of seafood and methylmercury exposures in the United States. Environ. Health Perspect. 126, 017006 (2018).

7. Munthe, J. et al. Recovery of mercury-contaminated fisheries. Ambio 36, 33–44 (2007).

8. Krabbenhoft, D. P. & Sunderland, E. M. Global change and mercury. Science 341, 1457–1458 (2013).

9. Eagles-Smith, C. A. et al. Modulators of mercury risk to wildlife and humans in the context of rapid global change. Ambio 47, 170–197 (2018).

10. Global Mercury Assessment 2013: Sources, Emissions, Releases and Environmental Transport, http://www.unep.org/hazardoussubstances/Mercury/Informationmaterials/ReportsandPublications/tabid/3593/Default.aspx (UNEP Chemicals Branch, 2013).

11. Clayden, M. G. et al. Mercury biomagnification through food webs is affected by physical and chemical characteristics of lakes. Environ. Sci. Technol. 47, 12047–12053 (2013).

12. Vander Zanden, M. J., Casselman, J. M. & Rasmussen, J. B. Stable isotope evidence for the food web consequences of species invasions in lakes. Nature 401, 464–467 (1999).

13. Kelly, E. N., Schindler, D. W., St. Louis, V. L., Donald, D. B. & Vladicka, K. E. Forest fire increases mercury accumulation by fishes via food web restructuring and increased mercury inputs. Proc. Natl Acad. Sci. USA 103, 19380–19385 (2006).

14. Schartup, A. T. et al. Climate change and overfishing increase neurotoxicant in marine predators. Nature 572, 648–650 (2019).

15. Lepak, R. F. et al. Mercury source changes and food web shifts alter contamination signatures of predatory fish from Lake Michigan. Proc. Natl Acad. Sci. USA 116, 23600–23608 (2019).

16. Hall, B. D., Bodaly, R. A., Fudge, R. J. P., Rudd, J. W. M. & Rosenberg, D. M. Food as the dominant pathway of methylmercury uptake by fish. Water Air Soil Pollut. 100, 13–24 (1997).

17. Guzzo, M. M., Blanchfield, P. J. & Rennie, M. D. Behavioral responses to annual temperature variation alter the dominant energy pathway, growth, and condition of a cold-water predator. Proc. Natl Acad. Sci. USA, 114, 9912–9917 (2017).

18. Gårdmark, A. & Huss, M. Individual variation and interactions explain food web responses to global warming. Phil. Trans. R. Soc. B 375, 20190449 (2020).

19. Harris, R. C. et al. Whole-ecosystem study shows rapid fish-mercury response to changes in mercury deposition. Proc. Natl Acad. Sci. USA 104, 16586–16591 (2007).

20. Depew, D. C. et al. An overview of mercury concentrations in freshwater fish species: a national fish mercury dataset for Canada. Can. J. Fish. Aquat. Sci. 70, 436–451 (2013).

21. St. Louis, V. L. et al. Atmospheric concentrations and wet/dry loadings of mercury at the remote Experimental Lakes Area, Northwestern Ontario, Canada. Environ. Sci. Technol. 53, 8017–8026 (2019).

22. Graydon, J. A. et al. The role of terrestrial vegetation in atmospheric Hg deposition: Pools and fluxes of spike and ambient Hg from the METAALICUS experiment. Global Biogeochem. Cycles 26, GB1022 (2012).

23. Oswald, C. J., Heyes, A. & Branfireun, B. A. Fate and transport of ambient mercury and applied mercury isotope in terrestrial upland soils: insights from the METAALICUS watershed. Environ. Sci. Technol. 48, 1023–1031 (2014).

24. Eckley, C. S. & Hintelmann, H. Determination of mercury methylation potentials in the water column of lakes across Canada. Sci. Total Environ. 368, 111–125 (2006).

25. 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. 46, 623–631 (2001).

26. Vander Zanden, M. J. & Vadeboncoeur, Y. Fishes as integrators of benthic and pelagic food webs in lakes. Ecology 83, 2152–2161 (2002).

27. Madenjian, C. P., Chipps, S. R. & Blanchfield, P. J. Time to refine mercury mass balance models for fish. FACETS 6, 272–286 (2021).

28. Van Walleghem, J. L. A., Blanchfield, P. J., Hrenchuk, L. E. & Hintelmann, H. Mercury elimination by a top predator, Esox lucius. Environ. Sci. Technol. 47, 4147–4154 (2013).

29. Knightes, C. D., Sunderland, E. M., Barber, M. C., Johnston, J. M. & Ambrose, R. B. Application of ecosystem‐scale fate and bioaccumulation models to predict fish mercury response times to changes in atmospheric deposition. Environ. Toxicol. Chem. 28, 881–893 (2009).

30. Chen, C. Y. et al. A critical time for mercury science to inform global policy. Environ. Sci. Technol. 52, 9556–9561 (2018).

31. Blanchfield, P. J., Paterson, M. J., Shearer, J. A. & Schindler, D. W. Johnson and Vallentyne’s legacy: 40 years of aquatic research at the Experimental Lakes Area. Can. J. Fish. Aquat. Sci. 66, 1831–1836 (2009).

32. Sandilands, K. A. et al. Application of Enriched Stable Mercury Isotopes to the Lake 658 Watershed for the METAALICUS Project, at the Experimental Lakes Area, Northwestern Ontario, Canada. Can. Tech. Rep. Fish. Aquat. Sci. vol. 2597 (Fisheries and Oceans Canada, 2005).

33. Sandilands, K. A. et al. Application of Enriched Stable Mercury Isotopes to the Lake 658 Watershed for the METAALICUS Project, at the Experimental Lakes Area, Northwestern Ontario, Canada, 2001−2007. Can. Tech. Rep. Fish. Aquat. Sci. vol. 2813 (Fisheries and Oceans Canada, 2008).

34. Bower, P. M. et al. Simultaneous measurement of primary production by whole-lake and bottle radiocarbon additions. Limnol. Oceanogr. 32, 299–312 (1987).

35. Golding, G. R. et al. Evidence for facilitated uptake of Hg(ii) by Vibrio anguillarum and Escherichia coli under anaerobic and aerobic conditions. Limnol. Oceanogr. 47, 967–975 (2002).

36. Southworth, G. et al. Evasion of added isotopic mercury from a northern temperate lake. Environ. Technol. Chem. 26, 53–60 (2007).

37. Baker, R. F., Blanchfield, P. J., Paterson, M. J., Flett, R. J. & Wesson, L. Evaluation of non-destructive biopsy methods for the analysis of mercury in fish tissue. Trans. Amer. Fish. Soc. 133, 568–576 (2004).

38. Paterson, M. J. et al. Bioaccumulation of newly-deposited mercury by fish and invertebrates: an enclosure study using stable mercury isotopes. Can. J. Fish. Aquat. Sci. 63, 2213–2224 (2006).

39. Hintelmann, H. & Ogrinc, N. in Biogeochemistry of Environmentally Important Trace Elements (eds Cai, Y. & Braids, O. C.) 321–338 (ACS, 2003).

40. Van Walleghem, J. L. A., Blanchfield, P. J. & Hintelmann, H. Elimination of mercury by yellow perch in the wild. Environ. Sci. Technol. 41, 5895–5901 (2007).

41. Hrenchuk, L. E., Blanchfield, P. J., Paterson, M. J. & Hintelmann, H. H. Dietary and waterborne mercury accumulation by yellow perch: a field experiment. Environ. Sci. Technol. 46, 509–516 (2012).

42. Howmiller, R. P. Effects of preservatives on weight of some common macrobenthic invertebrates. Trans. Amer. Fish. Soc. 101, 743–746 (1972).

43. Tremblay, G., Legendre, P., Doyon, J.-F., Verdon, R. & Schetagne, R. The use of polynomial regression analysis with indicator variables for interpretation of mercury in fish data. Biogeochemistry 40, 189–201 (1998).

44. Mills, K. H. & Beamish, R. J. Comparison of fin-ray and scale age determinations for lake whitefish (Coregonus clupeaformis) and their implications for estimates of growth and annual survival. Can. J. Fish. Aquat. Sci. 37, 534–544 (1980).

## Acknowledgements

We thank A. Robinson for aircraft support, M. Stainton and R. Hesslein for the preparation of mercury isotopes, M. Lyng, J. T. Bell, R. Burrell, B. Dimock, M. Dobrin, K. Fetterly, C. Miller, M. Rearick, J. Reid, G. Riedel, J. Shead, J. Zhou and staff and students at the Experimental Lakes Area for assistance with field collections and analysis of samples. K. Kidd, D. Orihel and J. Smol provided reviews of this manuscript. We thank the following agencies for their financial support for this study: Electric Power Research Institute, Environment and Climate Change Canada, Fisheries and Oceans Canada, Natural Sciences and Engineering Research Council of Canada, National Science Foundation Grant DEB 0451345 (C.C.G.) and DEB 0351050 (A.H. and C.C.G.), Southern Company, University of Alberta, U.S. Department of Energy, U.S. Environmental Protection Agency, and Wisconsin Focus on Energy Program Project Grant no. 4900-02-03 (J.P.H.). Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

## Author information

Authors

### Contributions

P.J.B., J.W.M.R., M.A., K.G.B., R.A.B., B.A.B., C.C.G., R.C.H., A.H., H.H., J.P.H., C.A.K., D.P.K., S.E.L., R.P.M., M.J.P., C.L.P. and V.L.S.L. contributed to the design of the whole-ecosystem experiment. The METAALICUS project was overseen by J.W.M.R. and R.C.H., who along with C.L.B., R.A.B., J.A.G., H.H., J.P.H., C.A.K., D.P.K., K.A.S., V.L.S.L. and M.T.T., applied mercury to the lake and watershed. P.J.B., L.E.H. and L.S.T. conducted the fish sampling. Water, sediment and lower trophic level data were collected and prepared by C.C.G., B.D.H., H.H., C.A.K., M.J.P., C.L.P. and K.A.S. Field samples in this study were analysed for mercury by C.C.G., H.H., D.P.K. and M.T.T. Data analyses were performed by P.J.B. and L.E.H. All authors collected and discussed project-level data that contributed to the interpretation of the data presented in this study. L.E.H. produced the figures and wrote the methods. P.J.B., J.W.M.R., C.A.K., and V.L.S.L. wrote the manuscript with input from all authors.

### Corresponding author

Correspondence to Paul J. Blanchfield.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Peer review information

Nature thanks Bridget Bergquist and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Extended data figures and tables

### Extended Data Fig. 1 Temporal dynamics of upland mercury in fish.

Isotopic 198Hg added to the upland area of Lake 658 was above the detection limit (0.5% of ambient MeHg; dashed line) in all fish species, but contributed little (<2%) to overall MeHg concentrations (percent increase = [upland spike MeHg]/[ambient MeHg] × 100). Mean annual MeHg concentration data for each species or age-class is presented and based on fish collected during fall population sampling (n = 1,052; sample size details in Extended Data Tables 2 and 3); dotted lines indicate missing data

.

### Extended Data Fig. 2 Comparison of individual and population body sizes of northern pike.

Mean (± s.e.m.) body size of all northern pike sampled in the fall of each year (population, black circles; n = 442) for muscle MeHg concentration using a biopsy method was stable over time. Individual northern pike (grey triangles) captured in 2007 and again in at least one subsequent year (n = 16 fish with 1–9 individuals recaptured each year 2008–2015) were used to determine individual losses of lake spike MeHg during recovery (see Fig. 3). These individual fish, which were also captured prior to 2007, showed an increase in body size over time (linear regression: y = 17.43x−34470.0, R2 = 0.55, F1,51 = 61.1, P < 0.0001)

.

## Rights and permissions

Reprints and Permissions

Blanchfield, P.J., Rudd, J.W.M., Hrenchuk, L.E. et al. Experimental evidence for recovery of mercury-contaminated fish populations. Nature 601, 74–78 (2022). https://doi.org/10.1038/s41586-021-04222-7

• Accepted:

• Published:

• Issue Date:

• DOI: https://doi.org/10.1038/s41586-021-04222-7