Pair your accounts.

Export articles to Mendeley

Get article recommendations from ACS based on references in your Mendeley library.

Pair your accounts.

Export articles to Mendeley

Get article recommendations from ACS based on references in your Mendeley library.

You’ve supercharged your research process with ACS and Mendeley!

STEP 1:
Click to create an ACS ID

Please note: If you switch to a different device, you may be asked to login again with only your ACS ID.

Please note: If you switch to a different device, you may be asked to login again with only your ACS ID.

Please note: If you switch to a different device, you may be asked to login again with only your ACS ID.

MENDELEY PAIRING EXPIRED
Your Mendeley pairing has expired. Please reconnect
ACS Publications. Most Trusted. Most Cited. Most Read
Trophic Dynamics of Mercury in the Baltic Archipelago Sea Food Web: The Impact of Ecological and Ecophysiological Traits
My Activity

Figure 1Loading Img
  • Open Access
Contaminants in Aquatic and Terrestrial Environments

Trophic Dynamics of Mercury in the Baltic Archipelago Sea Food Web: The Impact of Ecological and Ecophysiological Traits
Click to copy article linkArticle link copied!

  • Riikka K. Vainio*
    Riikka K. Vainio
    Department of Biology, University of Turku, FI-20014 Turku, Finland
    *Email: [email protected]
  • Veijo Jormalainen
    Veijo Jormalainen
    Department of Biology, University of Turku, FI-20014 Turku, Finland
  • Rune Dietz
    Rune Dietz
    Department of Ecoscience, Aarhus University, Frederiksborgvej 399, Postbox 358, DK-4000 Roskilde, Denmark
    More by Rune Dietz
  • Toni Laaksonen
    Toni Laaksonen
    Department of Biology, University of Turku, FI-20014 Turku, Finland
  • Ralf Schulz
    Ralf Schulz
    iES Landau, Institute for Environmental Sciences, University of Koblenz-Landau, Fortstrasse 7, DE-76829 Landau, Germany
    More by Ralf Schulz
  • Christian Sonne
    Christian Sonne
    Department of Ecoscience, Aarhus University, Frederiksborgvej 399, Postbox 358, DK-4000 Roskilde, Denmark
  • Jens Søndergaard
    Jens Søndergaard
    Department of Ecoscience, Aarhus University, Frederiksborgvej 399, Postbox 358, DK-4000 Roskilde, Denmark
  • Jochen P. Zubrod
    Jochen P. Zubrod
    iES Landau, Institute for Environmental Sciences, University of Koblenz-Landau, Fortstrasse 7, DE-76829 Landau, Germany
    Zubrod Environmental Data Science, Friesenstrasse 20, DE-76829, Landau, Germany
  • Igor Eulaers*
    Igor Eulaers
    Department of Ecoscience, Aarhus University, Frederiksborgvej 399, Postbox 358, DK-4000 Roskilde, Denmark
    Norwegian Polar Institute, FRAM Centre, Postboks 6606 Stakkevollan, NO-9296 Tromsø, Norway
    *Email: [email protected]
    More by Igor Eulaers
Open PDFSupporting Information (1)

Environmental Science & Technology

Cite this: Environ. Sci. Technol. 2022, 56, 16, 11440–11448
Click to copy citationCitation copied!
https://doi.org/10.1021/acs.est.2c03846
Published August 3, 2022

Copyright © 2022 The Authors. Published by American Chemical Society. This publication is licensed under

CC-BY 4.0 .

Abstract

Click to copy section linkSection link copied!

We investigated trophic dynamics of Hg in the polluted Baltic Archipelago Sea using established trophic magnification (TMFs) and biomagnification factors (BMFs) on a comprehensive set of bird, fish, and invertebrate species. As different ecological and ecophysiological species traits may affect trophic dynamics, we explored the effect of food chain (benthic, pelagic, benthopelagic) and thermoregulatory strategy on trophic total Hg (THg) dynamics, using different approaches to accommodate benthopelagic species and normalize for trophic position (TP). We observed TMFs and most BMFs greater than 1, indicating overall THg biomagnification. We found significantly higher pelagic TMFs (3.58–4.02) compared to benthic ones (2.11–2.34) when the homeotherm bird species were excluded from models, but not when included. This difference between the benthic and pelagic TMFs remained regardless of how the TP of benthopelagic species was modeled, or whether TMFs were normalized for TP or not. TP-corrected BMFs showed a larger range (0.44–508) compared to BMFs representing predator–prey concentration ratios (0.05–82.2). Overall, the present study shows the importance of including and evaluating the effect of ecological and ecophysiological traits when investigating trophic contaminant dynamics.

This publication is licensed under

CC-BY 4.0 .
  • cc licence
  • by licence
Copyright © 2022 The Authors. Published by American Chemical Society

Synopsis

Ecophysiological and ecological traits impact the output of established models for trophic contaminant dynamics, potentially having consequences for regulatory decisions based thereon.

1. Introduction

Click to copy section linkSection link copied!

Trophic dynamics of contaminants that show biomagnification potential have traditionally been investigated using trophic magnification factors (TMFs) and biomagnification factors (BMFs). (1) TMFs measure the average increase in contaminant concentration per trophic level, often proxied by bulk stable nitrogen isotopes (δ15N) due to their characteristic trophic enrichment. (2) BMFs in turn measure the increase in contaminant concentration from prey to predator. However, species’ ecological and ecophysiological traits, such as food chain origin, trophic position, and thermoregulatory strategy, may cause variation in the underlying ecological and contaminant variability when studying the trophic dynamics of contaminants, highlighting the need for trait-based modeling of contaminant dynamics. (3) Nonetheless, risk assessment modeling, typically using such TMFs and BMFs, does not commonly take into account such species traits nor their potential impact on modeling outcomes.
Ecotoxicological studies on marine systems have indicated that contaminant loadings and pathways differ between benthic and pelagic food chains due to for example differences in bioavailability. (4−8) Moreover, often food chains are assumed linear and discrete, while generalist species often feed on prey residing at multiple trophic positions or in different food chains, which may in fact be closely interlinked (for example benthopelagic coupling). Aside from all ecological species characteristics complicating practical discretization of food chains and their composing trophic positions, species exhibit ecophysiological traits, such as thermoregulation and respiration, that influence contaminant bioaccumulation and biomagnification. Indeed, species-differences in thermoregulatory strategy, resulting in differences in metabolism and energy requirements, can cause homeotherm species to accumulate higher contaminant concentrations than poikilotherms. (9−11) In fact, BMFs are currently approached with multiple different formulas, (12) some formulas normalizing for trophic position differences for the predator–prey pairs, while some do not. Even so, such potential impacts of different species-traits on BMFs, or TMFs, have not been consequently reported or robustly investigated for their importance when studying food web dynamics of contaminants.
The Baltic Sea is particularly susceptible to pollution, due to its geophysical location and hydrological properties (13) and has a well-known legacy for persistent contaminants that typically show biomagnification. (14,15) Among a plethora of hazardous substances found in the Baltic Sea is mercury (Hg (14)), for which emissions have increased dramatically since the Industrial Revolution. (16) Mercury is environmentally persistent, and concentrations are found to be elevated in most Baltic Sea basins (17) and readily bioavailable to the lower functional groups of the food web especially as organic methylmercury (MeHg), resulting in Hg bioaccumulation in higher trophic species and ultimately Hg biomagnification through food chains. (18)
Here we investigate the trophic dynamics of total Hg (THg) in benthic and pelagic food chains of the Baltic Archipelago Sea, using established TMFs and BMFs on a comprehensive set of species (n = 30), across different functional groups, that are key in this low diversity ecosystem. To investigate the impact of different species traits on trophic magnification models, we explored the effect of two ecological traits, i.e., the trophic position food chain origin (benthic, pelagic, benthopelagic) and one ecophysiological trait, i.e., thermoregulatory strategy (homeotherm, poikilotherm) on trophic THg dynamics, using different approaches to accommodate benthopelagic species and normalizing for trophic position (TP).

2. Materials and Methods

Click to copy section linkSection link copied!

2.1. Sample Collection

We collected samples of birds, fish, and invertebrates during the period 2017–2019, except for the extended time frame (2013–2019) for white-tailed eagle (Haliaeetus albicilla), in the Archipelago Sea, a Baltic Sea basin at the southwestern coast of Finland (Figure 1).

Figure 1

Figure 1. Map of the sampling locations.

The Archipelago Sea is one of the largest archipelagos and is a shallow brackish water ecosystem with low salinity, inhabited by both marine and freshwater species. Information on the sampling methods, collection times, tissues used for chemical analysis, and use of composite versus individual samples are presented in Table S1. We measured body length for invertebrates and fish, and individual fish mass was recorded when possible (Table S2). All birds were dissected for muscle. From larger fish species the dorsal fillet was dissected, while from smaller species the intestinal tract was removed, except from lesser sand eel (Ammodytes tobianus) and sand goby (Pomatoschistus minutus), before homogenization of the entire individual using a shear mixer. We pooled individual invertebrates sampled at the same location to constitute a composite sample with sufficient mass for chemical analysis. We transported all samples under cool and dark conditions during field collection and stored them at −20 °C prior to sample preparation for chemical analysis, prior to which samples were freeze-dried at −50 °C, gravimetrically determined for its dry matter, and finally homogenized using pestle and mortar.

2.2. Mercury Analysis

The analysis for total mercury was performed using a direct mercury analyzer (DMA-80, Milestone, Italy) at the accredited Trace Element Laboratory at Department of Ecoscience, Aarhus University, following the US-EPA Method 7473 (19) during the period 2019–2021. The analytical quality was controlled by concurrently analyzing aqueous standard solutions (10 ng and 100 ng of Hg, prepared from 1000 ± 4 mg L–1 stock solution, Sigma-Aldrich, Switzerland), procedural blanks, duplicates, and Certified Reference Material (DORM-4; fish protein from the National Research Council, Canada; 0.412 ± 0.036 μg Hg g–1 dry weight, dw). Measurements were subsequently corrected for instrumental drift using results from the aqueous 10 and 100 ng standards applied to the low- and high-concentration cell, respectively (drift was always < 10%). All samples were corrected for concurrent blanks (<0.1 ng). The measured recovery long-term percentage of the DORM-4 was 102 ± 6% (mean ± SD; n = 284 during 2019–2021). This shows satisfactory accuracy and a precision (extended analytical uncertainty) (2 SD) of ca. 10% for the analyses. The laboratory is accredited by the Danish Accreditation Fond DANAK following ISO 17025 to DMA-80 analyses of Hg in biota with a detection limit of 0.001 μg g–1 dw. As part of the laboratory QA/QC, the lab participates twice a year in the international laboratory performance study program QUASIMEME (www.quasimeme.org). The participation has proved excellent long-term measurement accuracy and precision (n = 15 during 2017–2021; assigned concentration from 0.009 to 0.931 μg g–1 dw; Z-scores from −1.0 to 0.7 with a mean of −0.1).

2.3. Stable Isotope Analysis

The stable isotope analysis was carried out at the Stable Isotope Lab of the University of Koblenz-Landau (Germany). Ratios of stable nitrogen isotopes (15N:14N) were determined on a freeze-dried and homogenized subsample (mean ± SD: 1.48 ± 0.06 mg) using a Flash 2000 HT elemental analyzer coupled via a ConFlo IV interface to a Delta V Advantage isotope ratio mass spectrometer (all Thermo Fisher Scientific, Bremen, Germany). The stable nitrogen isotope values are reported following best practices suggested by Bond and Hobson (20) and are relative to their respective international measurement standard of atmospheric N2. Internal reference material (casein) was measured concurrently in duplicate every 10 samples, revealing an imprecision (±SD) ≤ 0.06 ‰.

2.4. Data Analysis

To investigate the effect of ecological and ecophysiological traits on trophic magnification of THg, we categorized the study species, based on their diet, for their food chain (18 benthic, five pelagic, and seven benthopelagic species) and thermoregulatory strategy (three homeotherm, 27 poikilotherm). Rather than relying on the empirical stable isotope values we used literature on the studied species (see Supporting Information for full reference list) to define species to be pelagic or benthic, when mainly feeding in a food chain depending on pelagic or benthic basal energy sources, respectively, or to be benthopelagic when feeding in both.
We calculated the TP of an individual of a benthic and pelagic species following Post: (21)
(1)
where the trophic enrichment factor (TEF) Δ15N was set at 3.40, as recommended by Borgå et al., (1) and where δ15Nbaseline is the δ15N value for the chosen benthic or pelagic primary consumer species (assumed to feed at TP = 2), that is, Gammarus spp. and zooplankton, respectively for benthic and pelagic species. Primary consumer species are preferred as food chain baselines as primary producers’ δ15N values can show large spatiotemporal variability due to seasonal variation in growth rate and different nutrient inputs in different locations. (21) Due to this variability, we omitted primary producers from the trophic magnification models (See Supporting Information for more on the primary producers).
The calculation of TPs for benthopelagic species can however not follow this one-source model, as mixed feeding on both benthic and pelagic prey comes along with dependency on two δ15N baselines. Because of this, we calculated the TP for individuals of benthopelagic species following two approaches. In the first approach, we assumed that the benthopelagic species feeds entirely in either benthic or pelagic habitat, that is, from one basal energy source, and therefore the above-described one-source TP calculation was used for both food chains separately. In the second approach, we used a two-source model, assuming that the benthopelagic species feeds in both benthic and pelagic habitats, calculating the TP following the two-source formula outlined by Post: (21)
(2)
where δ15Nbaseline_B and δ15Nbaseline_P are the δ15N values of benthic and pelagic baselines, respectively. As the actual proportions in which benthopelagic species feed from benthic and pelagic food chains are unknown, we set the proportion of feeding from the first food chain to α = 0.50.
For birds, a TEF of Δ15N = 2.4 is more appropriate to represent trophic enrichment between bird tissue and that of their prey. (22) Hence, we calculated TPs for benthic common eiders in the one-source model following Hop et al.: (9)
(3)
while TPs for the benthopelagic white-tailed eagle and great cormorant were calculated after adaptation from Hop et al. (9) and Post: (21)
(4)
with α = 0.50. As the δ15N signal of whole fish and muscle tissue can show small differences, the TP estimated for whole fish might differ 0–0.6 trophic levels compared to TP determined from the muscle tissue. (23,24)
We calculated TMFs, also referred to as food web magnification factors (FWMFs), separately for the benthic and pelagic food chains, using both one-source and two-source models due to the included benthopelagic species. We used a linear mixed model with log10[THg dw] as dependent variable, TP or δ15N as independent variable, and species as a random variable due to unbalanced sample sizes across the studied species. Using the slope b of each fitted model we calculated the TMF following
(5)
The slope of the regression measures the increase in THg concentrations with the calculated TPs across the food chain, while the intercept represents the biological and environmental factors that determine the input of THg into the food chain at its base. (1) We used dw concentrations due to smaller variation of wet weight concentrations negatively effecting the reliability of the models. To assess the impact of homeotherm species on the trophic magnification models, we repeated all models excluding the three bird species. To compare the goodness of fit for the models, we calculated coefficients of determination for each model according to Nagawa et al.: (25) marginal Rm2, which represents the total variance explained by the fixed effects, and conditional Rc2, representing proportion of variance explained by both fixed and random effects. Differences in the slopes and intercepts between the benthic and pelagic food chain models were tested using a Welch’s t test.
Finally, to assess the effect of normalizing BMFs for TP, we calculated BMFs for all likely predator–prey pairs, based on literature recordings of the diet of the species (see Supporting Information), with normalization for TP:
(6)
and without:
(7)
The difference between the values produced by the two BMF approaches was calculated for each predator–prey pair by taking the absolute value of the subtraction of the BMFTP and BMFR. Due to consisting of only one individual, we omitted Atlantic salmon from the BMF calculations. Full BMF results with Atlantic salmon can be found in Tables S4–S7.

3. Results

Click to copy section linkSection link copied!

Summary statistics for δ15N, TP, and THg for all studied species are given in Table 1, and the extended summary statistics are given in Table S3. All TMFs exceeded one, indicating food web magnification (Figure 2; Table 2).

Figure 2

Figure 2. Linear regressions between log10[THg dw] (μg g–1) and trophic position for (A) benthic habitat (all species), (B) benthic habitat (no birds), (C) pelagic habitat (all species), and (D) pelagic habitat (no birds). The solid lines are for regression of the one-source model (eqs 1 and 3), while the dashed lines are for the two-source model (eqs 2 and 4).

Table 1. Mean ± SD Values for δ15N, TP, and THg Concentrations (μg g–1 ww and dw) for the Studied Species, Grouped According to Their Food Chaina
scientific namecommon namenδ15NTPHg dwHg ww
Pelagic
Ammodytes tobianuslesser sand eel1*+10.732.950.0140.0027
Clupea harengusherring10+13.00 ± 2.393.61 ± 0.700.15 ± 0.170.035 ± 0.041
Coregonus albulavendace2+12.29 ± 0.223.41 ± 0.060.087 ± 0.0130.018 ± 0.003
Salmo salarAtlantic salmon1+13.003.610.150.055
zooplankton 3*+7.52 ± 0.142.00 ± 0.040.0059 ± 0.00150.0006 ± 0.0002
Benthic
Somateria mollissimacommon eider13+10.32 ± 0.794.11 ± 0.230.67 ± 0.310.17 ± 0.09
Abramis bramacommon bream6+13.80 ± 1.374.84 ± 0.400.17 ± 0.060.042 ± 0.018
Gobius nigerblack goby5+11.80 ± 0.274.25 ± 0.080.028 ± 0.00870.0060 ± 0.0019
Gymnocephalus cernuaEurasian ruffe10+14.40 ± 0.695.01 ± 0.200.13 ± 0.080.026 ± 0.017
Myoxocphalus quadricornisfourhorn sculpin5+16.659 ± 0.755.66 ± 0.220.45 ± 0.380.083 ± 0.055
Neogobius melanostomusround goby10+10.45 ± 0.643.85 ± 0.190.018 ± 0.0050.0044 ± 0.0020
Pomatoschistus minutussand goby1*+13.734.820.0190.0045
Rutilus rutiluscommon roach10+11.67 ± 1.584.21 ± 0.460.15 ± 0.060.035 ± 0.014
Zoarces viviparusviviparous eelpout10+12.78 ± 0.494.54 ± 0.140.066 ± 0.0360.015 ± 0.008
Gammarus spp. 6*+4.15 ± 0.842.00 ± 0.250.018 ± 0.0260.0037 ± 0.0050
Idotea spp. 4*+4.25 ± 1.222.03 ± 0.360.0060 ± 0.00350.0016 ± 0.0010
Macoma balthicaBaltic clam4*+8.70 ± 0.293.34 ± 0.090.037 ± 0.0060.0074 ± 0.0013
Mytilus edulisblue mussel9*+7.16 ± 1.402.88 ± 0.410.037 ± 0.0130.0059 ± 0.0023
Palaemon adpersus 1+9.903.690.00860.0019
Palaemon elegans 1*+8.443.260.00910.0022
Rhithropanopeus harrisii 1*+8.243.200.0280.0071
Saduria entomon 2*+10.55 ± 0.793.88 ± 0.230.05 ± 0.030.016 ± 0.014
Theodoxus fluviatilis 3*+4.29 ± 0.532.0 ± 0.160.0048 ± 0.00500.0021 ± 0.0022
Benthopelagic
Haliaeetus albicillawhite-tailed eagle7+12.09 ± 1.833.64 ± 0.54 (P)1.9 ± 1.70.51 ± 0.34
4.63 ± 0.54 (B)
4.13 ± 0.54 (T)
Phalacrocorax carbogreat cormorant6+14.48 ± 1.864.34 ± 0.55 (P)1.3 ± 1.00.35 ± 0.27
5.33 ± 0.55 (B)
4.84 ± 0.55 (T)
Coregonus lavaretuscommon whitefish2+11.85 ± 0.163.27 ± 0.05 (P)0.084 ± 0.0490.019 ± 0.011
4.26 ± 0.05 (B)
3.77 ± 0.05 (T)
Esox luciusnorthern pike9+15.84 ± 0.464.45 ± 0.13 (P)0.68 ± 0.430.15 ± 0.09
5.44 ± 0.13 (B)
4.9 ± 0.13 (T)
Gasterosteus aculeatusthree-spined stickleback6+11.52 ± 0.473.18 ± 0.14 (P)0.039 ± 0.0180.010 ± 0.004
4.17 ± 0.14 (B)
3.67 ± 0.14 (T)
Osmerus eperlanusEuropean smelt10+12.78 ± 0.653.55 ± 0.22 (P)0.12 ± 0.080.024 ± 0.017
4.54 ± 0.19 (B)
4.04 ± 0.19 (T)
Perca fluviatilisEuropean perch10+13.73 ± 1.683.83 ± 0.49 (P)0.48 ± 0.380.095 ± 0.070
4.82 ± 0.49 (B)
4.32 ± 0.49 (T)
a

The TPs for benthopelagic species are calculated to comply with either the pelagic (P) or benthic (B) food chain baseline species both using a one-source model (eqs 1 and 3), or to accommodate both baseline species, using a two-source model (T) (eqs 2 and 4). Pooled samples are marked with an asterisk (∗).

Table 2. Statistical Output of Linear Models and the Resulting TMFs for THg (μg g–1 dw) in the Archipelago Sea, Accommodating Different Configurations to Test the Impact of the Species Traits Trophic Position, Food Chain Origin, and Thermoregulatory Strategya
modelnRm2Rc2pintercept (SE)slope (SE)TMF
Pelagic Food Chain
one-source700.140.78<0.01–2.15 (0.34)0.36 (0.08)2.27
one-source (no birds)540.550.70<0.01–3.15 (0.31)0.60 (0.09)4.02
two-source700.260.76<0.01–2.33 (0.32)0.38 (0.08)2.38
two-source (no birds)540.590.73<0.01–3.13 (0.30)0.55 (0.08)3.58
δ15Nb700.110.80<0.01–2.14 (0.36)0.10 (0.03)1.26
δ15N (no birds)b540.550.70<0.01–3.28 (0.33)0.18 (0.06)1.51
Benthic Food Chain
One-source1540.110.80<0.01–2.12 (0.29)0.23 (0.07)1.69
One-source (no birds)1250.420.71<0.01–2.81 (0.25)0.37 (0.06)2.34
Two-source1540.060.82<0.01–1.96 (0.30)0.19 (0.07)1.56
Two-source (no birds)1250.290.71<0.01–2.61 (0.28)0.32 (0.07)2.11
δ15Nb1540.080.81<0.01–1.87 (0.25)0.06 (0.02)1.15
δ15N (no birds)b1250.460.71<0.01–2.52 (0.21)0.11 (0.02)1.28
a

The TP of benthopelagic species was estimated using only pelagic or benthic baseline species (one-source model, eqs 1 and 3) or both (two-source model, eqs 2 and 4). n = sample size, Rm2 = marginal coefficient of determination, Rc2 = conditional coefficient of determination, p = significance for the slope, SE = standard error.

b

Note that the TMFs for the models using δ15N represent THg biomagnification per unit increase of δ15N (‰) instead of per trophic level.

We did not observe statistically significant differences between the pelagic one- and two-source model regardless of whether homeotherms were included or not (−0.38 ≤ t ≤ 0.44; 0.66 ≤ p ≤ 0.97). Similarly, we did not observe differences between benthic one- and two-source models with or without the homeotherms (−0.54 ≤ t ≤ 0.47; 0.59 ≤ p ≤ 0.71). The slopes or intercepts did not differ between the respective benthic models or the pelagic two-source models with or without the homeotherms (−0.21 ≤ t ≤ 1.76; 0.08 ≤ p ≤ 0.83; pelagic two-source: −1.57 ≤ t ≤ 1.76; 0.08 ≤ p ≤ 0.12). However, there was a difference between both slopes (t120 = −2.06; p = 0.04) and intercepts (t120 = 2.19; p = 0.03) for the pelagic one-source models with and without the homeotherms. We did not find differences between the benthic and pelagic models, irrespective of the one- or two-source approach, when homeotherms were included (−1.73 ≤ t ≤ 0.82; 0.09 ≤ p ≤ 0.96). There were also no differences between the model intercepts for the benthic and pelagic food chains, neither using a one-source model (t175 = 0.85; p = 0.40) nor a two-source one (t175 = 1.28; p = 0.20) when homeotherms were excluded.
However, in the models excluding homeotherm bird species, there was a significant difference in the slopes between the pelagic and benthic one-source (t175 = −2.24; p = 0.03) and two-source models (t175 = −2.20; p = 0.03). Regression slopes showed to be steeper in the pelagic food chain in both cases (Table 2). There was also a significant difference between benthic and pelagic slopes (t175 = −2.24; p = 0.03) as well as between the intercepts (t175 = −1.94; p = 0.05) of the TMFδ15N when excluding birds. Again, the slope showed to be steeper in the pelagic food chain, while intercept was higher in the benthic food chain (Table 2). Omitting birds from the models increased the explanatory power (Rm2; Table 2) compared to the respective models including these homeotherm species.
All obtained BMFs calculated using eqs 6 and 7 are given in Tables S4–S7. We observed 95.8% of BMFTP (eq 6) and 91.5% of BMFR (eq 7) to be greater than 1, indicating overall biomagnification. For homeotherms, we observed a range of BMFTP = 4.40–508 and BMFR = 1.48–56.8, while for poikilotherms BMFTP = 0.44–33.5, and BMFR = 0.19–82.2, values of BMFTP, but not BMFR being higher for homeotherms than poikilotherms. For pelagic predators, BMFTP = 1.61–9.35 and BMFR = 2.11–14.4, for benthic predators BMFTP = 0.44–30.9 and BMFR = 0.19–77.6, and for benthopelagic predators BMFTP = 1.61–508 and BMFR = 1.48–82.2, with highest BMFs in both cases found in the benthopelagic species, followed by benthic and pelagic species. Higher variation was observed for the obtained BMFTP (eq 6; mean ± SD = 14.5 ± 45.4; range, 0.44–508; Figure 3A), than for BMFR (eq 7; 11.2 ± 13.4; 0.19–82.2; Figure 3B). There was high variability in the difference between BMFTP and BMFR for the same predator–prey pair (mean ± SD = 11.5 ± 44.0; range, 0.004–494.2), the difference being greater than 10 in 19.9% of the pairs. The denominator in eq 6, that is, the normalization for TP, also was a cause for 9.6% of the predator–prey pairs negative BMFTP (removed from Figure 3A) due to the predator’s TP being lower than that of its prey (see Discussion in the Supporting Information).

Figure 3

Figure 3. Heatmap of biomagnification factors (BMFs) for relevant predator–prey pairs (A) normalized for TP (eq 6) or (B) not (eq 7). Negative BMF values were removed.

4. Discussion

Click to copy section linkSection link copied!

We investigated the impact of two ecological traits, trophic position and food chain origin, and an ecophysiological trait, thermoregulatory strategy, on established risk assessment modeling approaches, that is, TMFs and BMFs, to study trophic contaminant dynamics of mercury in the Baltic Archipelago Sea. Overall, TMFs and BMFs showed THg biomagnification to be affected by distinguishing food chains and the presence of homeotherm species. We also found most BMFs to be highly variable between different predator–prey pairs and formulas of calculation.
For both the benthic and pelagic food chains large positive TMF values showed strong THg food web magnification. The observed TMFδ15N values for the pelagic food chain (1.28–1.51) were similar to the TMFδ15N value of 1.50 reported for a pelagic food chain of zooplankton, Mysis spp., and herring in the Baltic Sea. (26) The trophic magnification slopes on δ15N in our study (0.06–0.18) also matched the average of all sites (mean ± SD: 0.16 ± 0.11), temperate sites (0.17 ± 0.10), and marine coastal sites (0.19 ± 0.08) evaluated in a worldwide meta-analysis by Lavoie et al. (27) The number of trophic levels (TPtop predator – TPbaseline) in both benthic (3.7) and pelagic (2.9) food chains in our study included more trophic levels than the studies (mean ± SD 1.7 ± 0.7) included in the same meta-analysis. (27) However, we found the TMFs derived separately for the benthic and pelagic food chains of the Archipelago Sea to reveal higher trophic magnification of THg in the pelagic food chain when excluding birds. This result is robust as it remained independently of how TPs of benthopelagic species were modeled. Similar to our study, higher TMFs in the pelagic food chain compared to those in the benthic one have been reported in the Gulf of St. Lawrence, Canada, (4) and Santos continental shelf, Brazil. (6) Pelagic pathways have been suggested to be more efficient in transferring highly bioaccumulative and biomagnifying MeHg compared to transfer from sediments, MeHg being more readily bioavailable in the pelagic food chain. (28,29) Although studies have shown both inorganic and organic species of Hg being capable of being transferred through the diet and assimilated into organisms’ tissue, both the rate of absorption and assimilation are higher for MeHg than inorganic species of Hg, (30) leading to higher rates of both bioaccumulation and biomagnification for MeHg, possibly contributing to the higher TMF observed for the pelagic food chain. Also, although the THg concentrations at the base of both benthic and pelagic food chains of our study were similar and dietary pathways dominate the Hg uptake in higher trophic levels (for example refs (30and31)), possible additional uptake of Hg through bioconcentration may partly explain the difference in TMFs between the benthic and pelagic food chains. These possibly different pathways of Hg into the benthic and pelagic food chain and proportions of MeHg should be further investigated to unravel their effect on the trophic dynamics of Hg. However, trophic contaminant dynamics between food chains are system-specific due to differences in for example physiochemical properties of the environment, and thus inherent variability in the resulting contamination pathways. Our results demonstrate the importance of investigating differences in trophic transfer of contaminants between different food chain pathways to better understand the contaminant dynamics.
In the present study, we modeled the TPs of benthopelagic species in two ways: assuming linearity of the food chain with its composing species feeding from only one basal energy source and allowing for the possibility of cross-food chain feeding depending on two basal energy sources. Often TP and trophic magnification models based thereon accommodate only one baseline species. However, most generalist and mobile species forage in different environments and food chains. Prey-switching can happen due to availability and ontogeny, and individuals can move between environments between seasons and stages of life cycle. (32−34) Recently, studies have investigated the coupling of benthic and pelagic systems (35) suggesting that benthic and pelagic habitats are coupled more closely in various ways than previously assumed. (36) The modeling approach (one- or two-source model) for the benthopelagic species did not affect the TMFs. While our two-source model approach, assuming benthopelagic species feeding 50% in either food chain, is perhaps still not entirely accurate as actual consumer-specific prey proportions would show, it offers a more ecologically realistic solution to accommodate species that show cross-food chain feeding strategies when using bulk stable isotopes for TP estimation. Moreover, applying an analysis for compound-specific stable isotopes (37,38) may provide additional insights and, perhaps, an even more robust assessment of the trophic dynamics of Hg in the studied food web.
Homeotherm species may accumulate higher concentrations of recalcitrant contaminants, such as Hg, than poikilotherms due to higher energy requirements demanding higher food consumption and thus intake of contaminants. (1,9) Moreover, as long-lived species, birds bioaccumulate recalcitrant contaminants during their lifetime. (39) Trophic magnification models including homeotherm birds did not show different trophic magnification between the benthic and pelagic food chain, in sharp contrast to the models excluding homeotherms. Moreover, all trophic magnification models without homeotherm species showed better statistical fit (Rm2) compared to models including homeotherms. However, omitting the homeotherm birds affected TMFs only for the pelagic one-source model, with TMF being higher without the birds. We found THg concentrations in homeotherms to be relatively elevated in comparison to species from other functional groups exhibiting similar TP. Similarly, Wang et al. (10) modeled MeHg biomagnification in floodplain food webs and found that the regression line for homeotherms had a higher intercept than the one for poikilotherms. While such a comparison should be made with precaution, as the present study cannot provide insight on MeHg dynamics, our findings also seem to suggest that the impact of thermoregulatory strategy on both within food chain and between food chain TMF estimates merits future consideration of designing and comparing trophic contaminant magnification modeling attempts.
More than 90% of all BMFs exceeded 1, further demonstrating THg biomagnification in the Archipelago Sea food web. For BMFTP, homeotherms had a larger range and higher maximum values than poikilotherms, while for BMFR ranges of homeotherms and poikilotherms were similar. When comparing BMFs across food chains, the highest BMFTP and BMFR were found for benthopelagic predators, and the lowest ones were found for pelagic predators. The total range of BMFs was large, especially BMFTP. This was not only due to differences in THg concentrations between predator and its prey, but also due to a small TP difference: small TP differences are therefore observed to inflate BMFTP. Conversely, large differences in δ15N-derived TPs may also lead to underestimation of the biomagnification potential. Large disparities in BMFs with and without the normalization for TPs may indicate that the actual diet of the predator includes also other sources of contamination than the prey used for calculating BMFs, (40) while also other factors, such as disparity in temporality, may affect BMFs. Currently, there seem to be no unbiased BMF estimations as each model (see Franklin (12)) has trade-offs between ecological reality and statistical applicability. Different ecological and physiological factors, such as uncertainties in species feeding ecology and metabolism, also affect BMFs, and cause variability in BMFs of the same contaminants between different studies. (12) Indeed, our results show that different ecological and ecophysiological species traits, such as the here studied food chain origin, TP, and thermoregulatory strategy, introduce uncertainty in BMF comparisons making comparisons within and across studies more difficult. Our finding of highly variable BMFs supports the recommendation of using TMFs rather than BMFs, as TMFs represent the average biomagnification across the food chain and are less likely to be affected by individual species traits than the BMFs. (41)
A recent risk assessment of THg concentrations on the Baltic Sea (42) showed the fish and bivalves being at low risk of Hg mediated health effects, while in white-tailed eagles the risk of negative health effects was assessed to be moderate to high. However, the results of the risk assessment should be taken as indicative, as variation in the proportions of highly toxic MeHg between both species and individuals can occur. (42) Also, climate change mediated increases in runoff combined with other factors, such as shorter icing period, may increase Hg fluxes into the Baltic Sea as well as Hg bioavailability, (43) altogether potentially increasing the loadings and changing the dynamics of Hg in the food web. Because of these concerns, it is important to understand the role of ecological and ecophysiological traits on trophic contaminant dynamics, as well as the vulnerability of current models used for regulatory decisions. Moreover, the novel functional ecological understanding and comparative model framework of the present study should be further investigated for implications due to Hg speciation dynamics.

Supporting Information

Click to copy section linkSection link copied!

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.est.2c03846.

  • Methods and discussion of primary producer and trophic positions of the species, literature references on the diet of the study species, Tables S1–S7 with additional information on the sampling, summary statistics, and full biomagnification factors (PDF)

Terms & Conditions

Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.

Author Information

Click to copy section linkSection link copied!

  • Corresponding Authors
    • Riikka K. Vainio - Department of Biology, University of Turku, FI-20014 Turku, FinlandOrcidhttps://orcid.org/0000-0003-2649-6799 Email: [email protected]
    • Igor Eulaers - Department of Ecoscience, Aarhus University, Frederiksborgvej 399, Postbox 358, DK-4000 Roskilde, DenmarkNorwegian Polar Institute, FRAM Centre, Postboks 6606 Stakkevollan, NO-9296 Tromsø, Norway Email: [email protected]
  • Authors
    • Veijo Jormalainen - Department of Biology, University of Turku, FI-20014 Turku, Finland
    • Rune Dietz - Department of Ecoscience, Aarhus University, Frederiksborgvej 399, Postbox 358, DK-4000 Roskilde, Denmark
    • Toni Laaksonen - Department of Biology, University of Turku, FI-20014 Turku, Finland
    • Ralf Schulz - iES Landau, Institute for Environmental Sciences, University of Koblenz-Landau, Fortstrasse 7, DE-76829 Landau, GermanyOrcidhttps://orcid.org/0000-0002-6348-6971
    • Christian Sonne - Department of Ecoscience, Aarhus University, Frederiksborgvej 399, Postbox 358, DK-4000 Roskilde, Denmark
    • Jens Søndergaard - Department of Ecoscience, Aarhus University, Frederiksborgvej 399, Postbox 358, DK-4000 Roskilde, Denmark
    • Jochen P. Zubrod - iES Landau, Institute for Environmental Sciences, University of Koblenz-Landau, Fortstrasse 7, DE-76829 Landau, GermanyZubrod Environmental Data Science, Friesenstrasse 20, DE-76829, Landau, Germany
  • Notes
    The authors declare no competing financial interest.

Acknowledgments

Click to copy section linkSection link copied!

The study was performed under the project BONUS BALTHEALTH. BONUS BALTHEALTH has received funding from BONUS (Art. 185), funded jointly by the EU, Innovation Fund Denmark (Grants 6180-00001B and 6180-00002B), Forschungszentrum Jülich GmbH, German Federal Ministry of Education and Research (Grant FKZ 03F0767A), Academy of Finland (Grant 311966) and Swedish Foundation for Strategic Environmental Research (MISTRA). Funding was also received from the Turku University Foundation (Grant 080552) and the Finnish Foundation for Nature Conservation and Finnish Cultural Foundation, Varsinais-Suomi Regional fund (Grant 85201728). The authors thank Finnish Natural History Museum, Minna Kivimäki, Olavi Sahlsten, Johannes Sahlsten, Luca Rugiu, Jarmo Väisänen, Aki Perälä, and Lasse Heimo for their help with the sample collection.

References

Click to copy section linkSection link copied!

This article references 43 other publications.

  1. 1
    Borgå, K.; Kidd, K. A.; Muir, D. C. G.; Berglund, O.; Conder, J. M.; Gobas, F.; Kucklick, J.; Malm, O.; Powell, D. E. Trophic magnification factors: Considerations of ecology, ecosystems, and study design. Integr Environ. Assess Manag 2012, 8, 6484,  DOI: 10.1002/ieam.244
  2. 2
    Cabana, G.; Rasmussen, J. B. Modelling food chain structure and contaminant bioaccumulation using stable nitrogen isotopes. Nature 1994, 372, 255257,  DOI: 10.1038/372255a0
  3. 3
    Van den Brink, P. J.; Baird, D. J.; Baveco, H. J. M.; Focks, A. The use of traits based approaches and eco(toxico)logical models to advance ecological risk assessment framework for chemicals. Integr Environ. Assess Manag 2013, 9, e47e57,  DOI: 10.1002/ieam.1443
  4. 4
    Lavoie, R. A.; Hebert, C. E.; Rail, J. F.; Braune, B. M.; Yumvihoze, E.; Hill, L. G.; Lean, D. R. S. Trophic structure and mercury distribution in a Gulf of St. Lawrence (Canada) food web using stable isotope analysis. Sci. Total Environ. 2010, 408, 55295539,  DOI: 10.1016/j.scitotenv.2010.07.053
  5. 5
    Chen, C. Y.; Borsuk, M. E.; Bugge, D. M.; Hollweg, T.; Balcom, P. H.; Ward, D. M.; Williams, J.; Mason, R. P. Benthic and pelagic pathways of methylmercury bioaccumulation in estuarine food webs of the Northeast United States. PLoS One 2014, 9, e89305,  DOI: 10.1371/journal.pone.0089305
  6. 6
    Muto, E. Y.; Soares, L. S. H.; Sarkis, J. E. S.; Hortellani, M. A.; Petti, M. A. V.; Corbisier, T. N. Biomagnification of mercury through the food web of the Santos continental shelf, subtropical Brazil. Mar. Ecol.: Prog. Ser. 2014, 512, 5569,  DOI: 10.3354/meps10892
  7. 7
    Ciesielski, T. M.; Pastukhov, M. V.; Leeves, S. A.; Farkas, J.; Lierhagen, S.; Poletaeva, V. I.; Jenssen, B. M. Differential bioaccumulation of potentially toxic elements in benthic and pelagic food chains in Lake Baikal. Environ. Sci. Pollut Res. 2016, 23, 1559315604,  DOI: 10.1007/s11356-016-6634-0
  8. 8
    Fan, S.; Wang, B.; Liu, H.; Gao, S.; Li, T.; Wang, S.; Liu, Y.; Liu, X.; Wan, Y. Trophodynamics of organic pollutants in pelagic and benthic food webs of Lake Dianchi: Importance of ingested sediment as uptake route. Environ. Sci. Technol. 2017, 51, 1413514143,  DOI: 10.1021/acs.est.7b03681
  9. 9
    Hop, H.; Borga, K.; Gabrielsen, G. W.; Kleivane, L.; Skaare, J. U. Food web magnification of persistent organic pollutants in poikilotherms and homeotherms from the Barents Sea. Environ. Sci. Technol. 2002, 36, 25892597,  DOI: 10.1021/es010231l
  10. 10
    Wang, J.; Newman, M. C.; Xu, X.; Liang, L. Higher and more variable methylmercury biomagnification factors for floodplain than the contiguous river (South River, Virginia USA). Ecotoxicol. Environ. Saf. 2013, 92, 191198,  DOI: 10.1016/j.ecoenv.2012.04.023
  11. 11
    Romero-Romero, S.; Herrero, L.; Fernandez, M.; Gomara, B.; Acuna, J. L. Biomagnification of persistent organic pollutants in a deep-sea, temperate food web. Sci. Total Environ. 2017, 605–606, 589597,  DOI: 10.1016/j.scitotenv.2017.06.148
  12. 12
    Franklin, J. How reliable are field-derived biomagnification factors and trophic magnification factors as indicators of bioaccumulation potential? Conclusions from a case study on per- and polyfluoroalkyl substances. Integr. Environ. Assess. Manag. 2016, 12, 620,  DOI: 10.1002/ieam.1642
  13. 13
    Reusch, T. B. H.; Dierking, J.; Andersson, H. C.; Bonsdorff, E.; Carstensen, J.; Casini, M.; Czajkowski, M.; Hasler, B.; Hinsby, K.; Hyytiäinen, K.; Johannesson, K.; Jomaa, S.; Jormalainen, V.; Kuosa, H.; Kurland, S.; Laikre, L.; MacKenzie, B. R.; Margonski, P.; Melzner, F.; Oesterwind, D.; Ojaveer, H.; Refsgaard, J. C.; Sandström, A.; Schwarz, G.; Tonderski, K.; Winder, M.; Zandersen, M. The Baltic Sea as a time machine for the future coastal ocean. Sci. Adv. 2018, 4, 8195,  DOI: 10.1126/sciadv.aar8195
  14. 14
    HELCOM Thematic Assessment of Hazardous Substances 2011–2016; Baltic Marine Environment Protection Commission, 2018; p 157.
  15. 15
    de Wit, C. A.; Bossi, R.; Dietz, R.; Dreyer, A.; Faxneld, S.; Garbus, S. E.; Hellström, P.; Koschorreck, J.; lohrmann, N.; Roos, A.; Sellström, U.; Sonne, C.; Treu, G.; Vorkamp, K.; Yuan, B.; Eulaers, I. Organohalogen compounds of emerging concern in Baltic Sea biota: Levels, biomagnification potential and comparisons with legacy contaminants. Environ. Int. 2020, 144, 106037,  DOI: 10.1016/j.envint.2020.106037
  16. 16
    Driscoll, C. T.; Mason, R. P.; Man Chan, H.; Jacob, D. J.; Pirrone, N. Mercury as a Global Pollutant: Sources, Pathways, and Effects Terms of Use. Environ. Sci. Technol. 2013, 47, 39,  DOI: 10.1021/es305071v
  17. 17
    HELCOM (2018b) Metals (lead, cadmium and mercury). HELCOM core indicator report.https://helcom.fi/media/core%20indicators/Metals-HELCOM-core-indicator-2018.pdf (accessed 2021-08-20).
  18. 18
    Morel, F. M. M.; Kraepiel, A. M. L.; Amyot, M. The chemical cycle and bioaccumulation of mercury. Annu. Rev. Ecol Syst 1998, 29, 543566,  DOI: 10.1146/annurev.ecolsys.29.1.543
  19. 19
    US-EPA. Method 7473 (SW-846): Mercury in Solids and Solutions by Thermal Decomposition, Amalgamation, and Atomic Absorption Spectrophotometry, rev. 0; US-EPA: Washington, DC, 1998; 17 pp (https://www.epa.gov/sites/production/files/2015-07/documents/epa-7473.pdf).
  20. 20
    Bond, A. L.; Hobson, K. A. Reporting stable-isotope ratios in ecology: recommended terminology, guidelines, and best practices. Waterbirds 2012, 35, 324331,  DOI: 10.1675/063.035.0213
  21. 21
    Post, D. M. Using stable isotopes to estimate trophic position: Models, methods, and assumptions. Ecology 2002, 83, 703718,  DOI: 10.1890/0012-9658(2002)083[0703:USITET]2.0.CO;2
  22. 22
    Mizutani, H.; Kabaya, Y.; Wada, E. Nitrogen and Carbon Isotope Compositions relate linearly in Cormorant Tissues and its Diet. Isot Environ. Heal Stud 1991, 27, 166168
  23. 23
    Kurle, C. M.; Sinclair, E. H.; Edwards, A. E.; Gudmundson, C. J. Temporal and spatial variation in the d15N and d13C values of fish and squid from Alaskan waters. Marine Biology 2011, 158, 23892404,  DOI: 10.1007/s00227-011-1741-4
  24. 24
    Fincel, M. J.; VanDeHey, J. A.; Wuestewald, A.; Chipps, S. R. Comparing isotope signatures of prey fish: does gut removal affect d13C or d15N?. Journal of Freshwater Ecology 2012, 27, 5562
  25. 25
    Nakagawa, S.; Johnson, P. C. D.; Schielzeth, H. The coefficient of determination and intra-class correlation coefficient from generalized linear mixed-effects models revisited and expanded. J. R Soc. Interface 2017, 14, 20170213,  DOI: 10.1098/rsif.2017.0213
  26. 26
    Nfon, E.; Cousins, I. T.; Järvinen, O.; Mukherjee, A. B.; Verta, M.; Broman, D. Trophodynamics of mercury and other trace elements in a pelagic food chain from the Baltic Sea. Sci. Total Environ. 2009, 407, 62676274,  DOI: 10.1016/j.scitotenv.2009.08.032
  27. 27
    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. 2013, 47, 1338513394,  DOI: 10.1021/es403103t
  28. 28
    Chen, C. Y.; Dionne, M.; Mayes, B. M.; Ward, D. M.; Sturup, S.; Jackson, B. P. Mercury bioavailability and bioaccumulation in estuarine food webs in the Gulf of Maine. Environ. Sci. Technol. 2009, 43, 18041810,  DOI: 10.1021/es8017122
  29. 29
    Chen, C. Y.; Borsuk, M. E.; Bugge, D. M.; Hollweg, T.; Balcom, P. H.; Ward, D. M.; Williams, J.; Mason, R. P. Benthic and pelagic pathways of methylmercury bioaccumulation in estuarine food webs of the northeast Unitd States. PLoS One 2014, 9, e89305,  DOI: 10.1371/journal.pone.0089305
  30. 30
    Bradley, M. A.; Barst, B. D.; Basu, N. A review of mercury bioavailability in humans and fish. International Journal of Environmental Research and Public Health 2017, 14, 169,  DOI: 10.3390/ijerph14020169
  31. 31
    Hall, B. D.; Bodaly, R. A.; Fudge, R. J.; Rudd, J. W. M.; Rosenberg, D. M. Food as the dominant pathway of methylmercury uptake by fish. Water, Air and Soil Pollution 1997, 100, 1324,  DOI: 10.1023/A:1018071406537
  32. 32
    Lankov, A.; Ojaveer, H.; Simm, M.; Põllupüü, M.; Möllmann, C. Feeding ecology of pelagic fish species in the Gulf of Riga (Baltic Sea): the importance of changes in the zooplankton community. J. Fish Biol. 2010, 77, 22682284,  DOI: 10.1111/j.1095-8649.2010.02805.x
  33. 33
    Jacobson, P.; Bengström, U.; Eklöf, J. Size-dependent diet composition and feeding of Eurasian perch (Perca fluviatilis) and northern pike (Esox lucius) in the Baltic Sea. Boreal. Environ. Res. 2019, 24, 137153
  34. 34
    Bergström, U.; Olsson, J.; Casini, M.; Eriksson, B. K.; Fredriksson, R.; Weenhage, H.; Appelberg, M. Stickleback increase in the Baltic Sea - A thorny issue for coastal predatory fish. Estuar Coast Shelf Sci. 2015, 163, 134142,  DOI: 10.1016/j.ecss.2015.06.017
  35. 35
    Griffiths, J. R.; Kadin, M.; Nascimento, F. J. A.; Tamelander, T.; Törnroos, A.; Bonaglia, S.; Bonsdorff, E.; Brüchert, V.; Gårdmark, A.; Järnström, M.; Kotta, J.; Lindegren, M.; Nordström, M. C.; Norkko, A.; Olsson, J.; Weigel, B.; Zydelis, R.; Blenckner, T.; Niiranen, S.; Winder, M. The importance of benthic–pelagic coupling for marine ecosystem functioning in a changing world. Glob. Chang. Biol. 2017, 23, 21792196,  DOI: 10.1111/gcb.13642
  36. 36
    Kiljunen, M.; Peltonen, H.; Lehtiniemi, M.; Uusitalo, L.; Sinisalo, T.; Norkko, J.; Kunnasranta, M.; Torniainen, J.; Rissanen, A. J.; Karjalainen, J. Benthic-pelagic coupling and trophic relationships in northern Baltic Sea food webs. Limnol Oceanogr 2020, 65, 17061722,  DOI: 10.1002/lno.11413
  37. 37
    Won, E. J.; Choi, B.; Hong, S.; Khim, J. S.; Shin, K. H. Importance of accurate trophic level determination by nitrogen isotope of amino acids for trophic magnification studies: A review. Environ. Pollut. 2018, 238, 677690,  DOI: 10.1016/j.envpol.2018.03.045
  38. 38
    Elliott, H.; Braune, B. M.; Elliott, J. E. Beyond bulk δ15N: Combining a suite of stable isotopic measures improves the resolution of the food webs mediating contaminant signals across space, time and communities. Environ. Int. 2021, 148, 106370,  DOI: 10.1016/j.envint.2020.106370
  39. 39
    Isomursu, M.; Koivusaari, J.; Stjernberg, T.; Hirvelä-Koski, V.; Venäläinen, E.-R. Lead poisoning and other human-related factors cause significant mortality in white-tailed eagles. Ambio 2018, 111
  40. 40
    Borgå, K.; Fisk, A. T.; Hoekstra, P. F.; Muir, D. C. G. Biological and chemical factors of importance in the bioaccumulation and trophic transfer of persistent organochlorine contaminants in arctic marine food webs. Environ. Toxicol. Chem. 2004, 23, 23672385,  DOI: 10.1897/03-518
  41. 41
    Gobas, F. A.P.C.; de Wolf, W.; Burkhard, L. P.; Verbruggen, E.; Plotzke, K. Revisiting bioaccumulation criteria for POPs and PBT assessments. Integ Environ. Assess Manag 2009, 5, 624637,  DOI: 10.1897/IEAM_2008-089.1
  42. 42
    Dietz, R.; Fort, J.; Sonne, C.; Albert, C.; Bustnes, J. O.; Christensen, T. K.; Ciesielski, T. M.; Danielsen, J.; Dastnai, S.; Eens, M.; Erikstad, K. E.; Galatius, A.; Garbus, S. E.; Gilg, O.; Hanssen, S. E.; Helander, B.; Helberg, M.; Jaspers, V. L. B.; Jenssen, B. M.; Jonsson, J. E.; Kauhala, K.; Kolbeinsson, Y.; Kyhn, L. A.; Labansen, A. L.; Larsen, M. M.; Linstrøm, U.; Reiertsen, T. K.; Rigét, F. F.; Roos, A.; Strand, J.; Strøm, H.; Sveegaard, S.; Søndergaard, J.; Sun, J.; Teilmann, J.; Therkildsen, O. R.; Thórarinsson, T. L.; Tjørnløv, R. S.; Wilson, S.; Eulaers, I. A risk assessment of the effects of mercury on Baltic Sea, Greater North Sea and North Atlantic wildlife, fish and bivalves. Environ. Int. 2021, 146, 106178,  DOI: 10.1016/j.envint.2020.106178
  43. 43
    HELCOM. Inputs of hazardous substances to the Baltic Sea. Baltic Sea Environmental Proceedings No. 162; Baltic Marine Environment Protection Commission, 2018a.

Cited By

Click to copy section linkSection link copied!

This article is cited by 2 publications.

  1. Kang Luo, Wei Yuan, Zhiyun Lu, Zichun Xiong, Che-Jen Lin, Xun Wang, Xinbin Feng. Unveiling the Sources and Transfer of Mercury in Forest Bird Food Chains Using Techniques of Vivo-Nest Video Recording and Stable Isotopes. Environmental Science & Technology 2024, 58 (13) , 6007-6018. https://doi.org/10.1021/acs.est.3c10972
  2. Xiaobo Zheng, Xiaodan Wu, Ruifeng Lu, Xingpei Cao, Bi-Xian Mai. Identification of Species-Specific Prey Uptake and Biotransformation of Chiral Polychlorinated Biphenyls (PCBs) in Riparian and Aquatic Food Webs. Environmental Science & Technology 2023, 57 (48) , 20282-20291. https://doi.org/10.1021/acs.est.3c07377
Open PDF

Environmental Science & Technology

Cite this: Environ. Sci. Technol. 2022, 56, 16, 11440–11448
Click to copy citationCitation copied!
https://doi.org/10.1021/acs.est.2c03846
Published August 3, 2022

Copyright © 2022 The Authors. Published by American Chemical Society. This publication is licensed under

CC-BY 4.0 .

Article Views

1502

Altmetric

-

Citations

Learn about these metrics

Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.

Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.

The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.

  • Abstract

    Figure 1

    Figure 1. Map of the sampling locations.

    Figure 2

    Figure 2. Linear regressions between log10[THg dw] (μg g–1) and trophic position for (A) benthic habitat (all species), (B) benthic habitat (no birds), (C) pelagic habitat (all species), and (D) pelagic habitat (no birds). The solid lines are for regression of the one-source model (eqs 1 and 3), while the dashed lines are for the two-source model (eqs 2 and 4).

    Figure 3

    Figure 3. Heatmap of biomagnification factors (BMFs) for relevant predator–prey pairs (A) normalized for TP (eq 6) or (B) not (eq 7). Negative BMF values were removed.

  • References


    This article references 43 other publications.

    1. 1
      Borgå, K.; Kidd, K. A.; Muir, D. C. G.; Berglund, O.; Conder, J. M.; Gobas, F.; Kucklick, J.; Malm, O.; Powell, D. E. Trophic magnification factors: Considerations of ecology, ecosystems, and study design. Integr Environ. Assess Manag 2012, 8, 6484,  DOI: 10.1002/ieam.244
    2. 2
      Cabana, G.; Rasmussen, J. B. Modelling food chain structure and contaminant bioaccumulation using stable nitrogen isotopes. Nature 1994, 372, 255257,  DOI: 10.1038/372255a0
    3. 3
      Van den Brink, P. J.; Baird, D. J.; Baveco, H. J. M.; Focks, A. The use of traits based approaches and eco(toxico)logical models to advance ecological risk assessment framework for chemicals. Integr Environ. Assess Manag 2013, 9, e47e57,  DOI: 10.1002/ieam.1443
    4. 4
      Lavoie, R. A.; Hebert, C. E.; Rail, J. F.; Braune, B. M.; Yumvihoze, E.; Hill, L. G.; Lean, D. R. S. Trophic structure and mercury distribution in a Gulf of St. Lawrence (Canada) food web using stable isotope analysis. Sci. Total Environ. 2010, 408, 55295539,  DOI: 10.1016/j.scitotenv.2010.07.053
    5. 5
      Chen, C. Y.; Borsuk, M. E.; Bugge, D. M.; Hollweg, T.; Balcom, P. H.; Ward, D. M.; Williams, J.; Mason, R. P. Benthic and pelagic pathways of methylmercury bioaccumulation in estuarine food webs of the Northeast United States. PLoS One 2014, 9, e89305,  DOI: 10.1371/journal.pone.0089305
    6. 6
      Muto, E. Y.; Soares, L. S. H.; Sarkis, J. E. S.; Hortellani, M. A.; Petti, M. A. V.; Corbisier, T. N. Biomagnification of mercury through the food web of the Santos continental shelf, subtropical Brazil. Mar. Ecol.: Prog. Ser. 2014, 512, 5569,  DOI: 10.3354/meps10892
    7. 7
      Ciesielski, T. M.; Pastukhov, M. V.; Leeves, S. A.; Farkas, J.; Lierhagen, S.; Poletaeva, V. I.; Jenssen, B. M. Differential bioaccumulation of potentially toxic elements in benthic and pelagic food chains in Lake Baikal. Environ. Sci. Pollut Res. 2016, 23, 1559315604,  DOI: 10.1007/s11356-016-6634-0
    8. 8
      Fan, S.; Wang, B.; Liu, H.; Gao, S.; Li, T.; Wang, S.; Liu, Y.; Liu, X.; Wan, Y. Trophodynamics of organic pollutants in pelagic and benthic food webs of Lake Dianchi: Importance of ingested sediment as uptake route. Environ. Sci. Technol. 2017, 51, 1413514143,  DOI: 10.1021/acs.est.7b03681
    9. 9
      Hop, H.; Borga, K.; Gabrielsen, G. W.; Kleivane, L.; Skaare, J. U. Food web magnification of persistent organic pollutants in poikilotherms and homeotherms from the Barents Sea. Environ. Sci. Technol. 2002, 36, 25892597,  DOI: 10.1021/es010231l
    10. 10
      Wang, J.; Newman, M. C.; Xu, X.; Liang, L. Higher and more variable methylmercury biomagnification factors for floodplain than the contiguous river (South River, Virginia USA). Ecotoxicol. Environ. Saf. 2013, 92, 191198,  DOI: 10.1016/j.ecoenv.2012.04.023
    11. 11
      Romero-Romero, S.; Herrero, L.; Fernandez, M.; Gomara, B.; Acuna, J. L. Biomagnification of persistent organic pollutants in a deep-sea, temperate food web. Sci. Total Environ. 2017, 605–606, 589597,  DOI: 10.1016/j.scitotenv.2017.06.148
    12. 12
      Franklin, J. How reliable are field-derived biomagnification factors and trophic magnification factors as indicators of bioaccumulation potential? Conclusions from a case study on per- and polyfluoroalkyl substances. Integr. Environ. Assess. Manag. 2016, 12, 620,  DOI: 10.1002/ieam.1642
    13. 13
      Reusch, T. B. H.; Dierking, J.; Andersson, H. C.; Bonsdorff, E.; Carstensen, J.; Casini, M.; Czajkowski, M.; Hasler, B.; Hinsby, K.; Hyytiäinen, K.; Johannesson, K.; Jomaa, S.; Jormalainen, V.; Kuosa, H.; Kurland, S.; Laikre, L.; MacKenzie, B. R.; Margonski, P.; Melzner, F.; Oesterwind, D.; Ojaveer, H.; Refsgaard, J. C.; Sandström, A.; Schwarz, G.; Tonderski, K.; Winder, M.; Zandersen, M. The Baltic Sea as a time machine for the future coastal ocean. Sci. Adv. 2018, 4, 8195,  DOI: 10.1126/sciadv.aar8195
    14. 14
      HELCOM Thematic Assessment of Hazardous Substances 2011–2016; Baltic Marine Environment Protection Commission, 2018; p 157.
    15. 15
      de Wit, C. A.; Bossi, R.; Dietz, R.; Dreyer, A.; Faxneld, S.; Garbus, S. E.; Hellström, P.; Koschorreck, J.; lohrmann, N.; Roos, A.; Sellström, U.; Sonne, C.; Treu, G.; Vorkamp, K.; Yuan, B.; Eulaers, I. Organohalogen compounds of emerging concern in Baltic Sea biota: Levels, biomagnification potential and comparisons with legacy contaminants. Environ. Int. 2020, 144, 106037,  DOI: 10.1016/j.envint.2020.106037
    16. 16
      Driscoll, C. T.; Mason, R. P.; Man Chan, H.; Jacob, D. J.; Pirrone, N. Mercury as a Global Pollutant: Sources, Pathways, and Effects Terms of Use. Environ. Sci. Technol. 2013, 47, 39,  DOI: 10.1021/es305071v
    17. 17
      HELCOM (2018b) Metals (lead, cadmium and mercury). HELCOM core indicator report.https://helcom.fi/media/core%20indicators/Metals-HELCOM-core-indicator-2018.pdf (accessed 2021-08-20).
    18. 18
      Morel, F. M. M.; Kraepiel, A. M. L.; Amyot, M. The chemical cycle and bioaccumulation of mercury. Annu. Rev. Ecol Syst 1998, 29, 543566,  DOI: 10.1146/annurev.ecolsys.29.1.543
    19. 19
      US-EPA. Method 7473 (SW-846): Mercury in Solids and Solutions by Thermal Decomposition, Amalgamation, and Atomic Absorption Spectrophotometry, rev. 0; US-EPA: Washington, DC, 1998; 17 pp (https://www.epa.gov/sites/production/files/2015-07/documents/epa-7473.pdf).
    20. 20
      Bond, A. L.; Hobson, K. A. Reporting stable-isotope ratios in ecology: recommended terminology, guidelines, and best practices. Waterbirds 2012, 35, 324331,  DOI: 10.1675/063.035.0213
    21. 21
      Post, D. M. Using stable isotopes to estimate trophic position: Models, methods, and assumptions. Ecology 2002, 83, 703718,  DOI: 10.1890/0012-9658(2002)083[0703:USITET]2.0.CO;2
    22. 22
      Mizutani, H.; Kabaya, Y.; Wada, E. Nitrogen and Carbon Isotope Compositions relate linearly in Cormorant Tissues and its Diet. Isot Environ. Heal Stud 1991, 27, 166168
    23. 23
      Kurle, C. M.; Sinclair, E. H.; Edwards, A. E.; Gudmundson, C. J. Temporal and spatial variation in the d15N and d13C values of fish and squid from Alaskan waters. Marine Biology 2011, 158, 23892404,  DOI: 10.1007/s00227-011-1741-4
    24. 24
      Fincel, M. J.; VanDeHey, J. A.; Wuestewald, A.; Chipps, S. R. Comparing isotope signatures of prey fish: does gut removal affect d13C or d15N?. Journal of Freshwater Ecology 2012, 27, 5562
    25. 25
      Nakagawa, S.; Johnson, P. C. D.; Schielzeth, H. The coefficient of determination and intra-class correlation coefficient from generalized linear mixed-effects models revisited and expanded. J. R Soc. Interface 2017, 14, 20170213,  DOI: 10.1098/rsif.2017.0213
    26. 26
      Nfon, E.; Cousins, I. T.; Järvinen, O.; Mukherjee, A. B.; Verta, M.; Broman, D. Trophodynamics of mercury and other trace elements in a pelagic food chain from the Baltic Sea. Sci. Total Environ. 2009, 407, 62676274,  DOI: 10.1016/j.scitotenv.2009.08.032
    27. 27
      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. 2013, 47, 1338513394,  DOI: 10.1021/es403103t
    28. 28
      Chen, C. Y.; Dionne, M.; Mayes, B. M.; Ward, D. M.; Sturup, S.; Jackson, B. P. Mercury bioavailability and bioaccumulation in estuarine food webs in the Gulf of Maine. Environ. Sci. Technol. 2009, 43, 18041810,  DOI: 10.1021/es8017122
    29. 29
      Chen, C. Y.; Borsuk, M. E.; Bugge, D. M.; Hollweg, T.; Balcom, P. H.; Ward, D. M.; Williams, J.; Mason, R. P. Benthic and pelagic pathways of methylmercury bioaccumulation in estuarine food webs of the northeast Unitd States. PLoS One 2014, 9, e89305,  DOI: 10.1371/journal.pone.0089305
    30. 30
      Bradley, M. A.; Barst, B. D.; Basu, N. A review of mercury bioavailability in humans and fish. International Journal of Environmental Research and Public Health 2017, 14, 169,  DOI: 10.3390/ijerph14020169
    31. 31
      Hall, B. D.; Bodaly, R. A.; Fudge, R. J.; Rudd, J. W. M.; Rosenberg, D. M. Food as the dominant pathway of methylmercury uptake by fish. Water, Air and Soil Pollution 1997, 100, 1324,  DOI: 10.1023/A:1018071406537
    32. 32
      Lankov, A.; Ojaveer, H.; Simm, M.; Põllupüü, M.; Möllmann, C. Feeding ecology of pelagic fish species in the Gulf of Riga (Baltic Sea): the importance of changes in the zooplankton community. J. Fish Biol. 2010, 77, 22682284,  DOI: 10.1111/j.1095-8649.2010.02805.x
    33. 33
      Jacobson, P.; Bengström, U.; Eklöf, J. Size-dependent diet composition and feeding of Eurasian perch (Perca fluviatilis) and northern pike (Esox lucius) in the Baltic Sea. Boreal. Environ. Res. 2019, 24, 137153
    34. 34
      Bergström, U.; Olsson, J.; Casini, M.; Eriksson, B. K.; Fredriksson, R.; Weenhage, H.; Appelberg, M. Stickleback increase in the Baltic Sea - A thorny issue for coastal predatory fish. Estuar Coast Shelf Sci. 2015, 163, 134142,  DOI: 10.1016/j.ecss.2015.06.017
    35. 35
      Griffiths, J. R.; Kadin, M.; Nascimento, F. J. A.; Tamelander, T.; Törnroos, A.; Bonaglia, S.; Bonsdorff, E.; Brüchert, V.; Gårdmark, A.; Järnström, M.; Kotta, J.; Lindegren, M.; Nordström, M. C.; Norkko, A.; Olsson, J.; Weigel, B.; Zydelis, R.; Blenckner, T.; Niiranen, S.; Winder, M. The importance of benthic–pelagic coupling for marine ecosystem functioning in a changing world. Glob. Chang. Biol. 2017, 23, 21792196,  DOI: 10.1111/gcb.13642
    36. 36
      Kiljunen, M.; Peltonen, H.; Lehtiniemi, M.; Uusitalo, L.; Sinisalo, T.; Norkko, J.; Kunnasranta, M.; Torniainen, J.; Rissanen, A. J.; Karjalainen, J. Benthic-pelagic coupling and trophic relationships in northern Baltic Sea food webs. Limnol Oceanogr 2020, 65, 17061722,  DOI: 10.1002/lno.11413
    37. 37
      Won, E. J.; Choi, B.; Hong, S.; Khim, J. S.; Shin, K. H. Importance of accurate trophic level determination by nitrogen isotope of amino acids for trophic magnification studies: A review. Environ. Pollut. 2018, 238, 677690,  DOI: 10.1016/j.envpol.2018.03.045
    38. 38
      Elliott, H.; Braune, B. M.; Elliott, J. E. Beyond bulk δ15N: Combining a suite of stable isotopic measures improves the resolution of the food webs mediating contaminant signals across space, time and communities. Environ. Int. 2021, 148, 106370,  DOI: 10.1016/j.envint.2020.106370
    39. 39
      Isomursu, M.; Koivusaari, J.; Stjernberg, T.; Hirvelä-Koski, V.; Venäläinen, E.-R. Lead poisoning and other human-related factors cause significant mortality in white-tailed eagles. Ambio 2018, 111
    40. 40
      Borgå, K.; Fisk, A. T.; Hoekstra, P. F.; Muir, D. C. G. Biological and chemical factors of importance in the bioaccumulation and trophic transfer of persistent organochlorine contaminants in arctic marine food webs. Environ. Toxicol. Chem. 2004, 23, 23672385,  DOI: 10.1897/03-518
    41. 41
      Gobas, F. A.P.C.; de Wolf, W.; Burkhard, L. P.; Verbruggen, E.; Plotzke, K. Revisiting bioaccumulation criteria for POPs and PBT assessments. Integ Environ. Assess Manag 2009, 5, 624637,  DOI: 10.1897/IEAM_2008-089.1
    42. 42
      Dietz, R.; Fort, J.; Sonne, C.; Albert, C.; Bustnes, J. O.; Christensen, T. K.; Ciesielski, T. M.; Danielsen, J.; Dastnai, S.; Eens, M.; Erikstad, K. E.; Galatius, A.; Garbus, S. E.; Gilg, O.; Hanssen, S. E.; Helander, B.; Helberg, M.; Jaspers, V. L. B.; Jenssen, B. M.; Jonsson, J. E.; Kauhala, K.; Kolbeinsson, Y.; Kyhn, L. A.; Labansen, A. L.; Larsen, M. M.; Linstrøm, U.; Reiertsen, T. K.; Rigét, F. F.; Roos, A.; Strand, J.; Strøm, H.; Sveegaard, S.; Søndergaard, J.; Sun, J.; Teilmann, J.; Therkildsen, O. R.; Thórarinsson, T. L.; Tjørnløv, R. S.; Wilson, S.; Eulaers, I. A risk assessment of the effects of mercury on Baltic Sea, Greater North Sea and North Atlantic wildlife, fish and bivalves. Environ. Int. 2021, 146, 106178,  DOI: 10.1016/j.envint.2020.106178
    43. 43
      HELCOM. Inputs of hazardous substances to the Baltic Sea. Baltic Sea Environmental Proceedings No. 162; Baltic Marine Environment Protection Commission, 2018a.
  • Supporting Information

    Supporting Information


    The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.est.2c03846.

    • Methods and discussion of primary producer and trophic positions of the species, literature references on the diet of the study species, Tables S1–S7 with additional information on the sampling, summary statistics, and full biomagnification factors (PDF)


    Terms & Conditions

    Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.