Katrine
Borgå‡
*a,
Melissa A.
McKinney‡
*b,
Heli
Routti
c,
Kim J.
Fernie
d,
Julia
Giebichenstein
a,
Ingeborg
Hallanger
c and
Derek C. G.
Muir
e
aDepartment of Biosciences, University of Oslo, NO-0316 Oslo, Norway. E-mail: katrine.borga@ibv.uio.no
bDepartment of Natural Resource Sciences, McGill University, Sainte-Anne-de-Bellevue, QC H9X 3 V9, Canada. E-mail: melissa.mckinney@mcgill.ca
cNorwegian Polar Institute, Fram Centre, NO-9296 Tromsø, Norway
dEcotoxicology & Wildlife Health, Environment and Climate Change Canada, Burlington, ON L7S 1A1, Canada
eAquatic Contaminants Research Division, Environment and Climate Change Canada, Burlington, ON L7S 1A1, Canada
First published on 18th February 2022
This review summarizes current understanding of how climate change-driven physical and ecological processes influence the levels of persistent organic pollutants (POPs) and contaminants of emerging Arctic concern (CEACs) in Arctic biota and food webs. The review also highlights how climate change may interact with other stressors to impact contaminant toxicity, and the utility of modeling and newer research tools in closing knowledge gaps on climate change-contaminant interactions. Permafrost thaw is influencing the concentrations of POPs in freshwater ecosystems. Physical climate parameters, including climate oscillation indices, precipitation, water salinity, sea ice age, and sea ice quality show statistical associations with POPs concentrations in multiple Arctic biota. Northward range-shifting species can act as biovectors for POPs and CEACs into Arctic marine food webs. Shifts in trophic position can alter POPs concentrations in populations of Arctic species. Reductions in body condition are associated with increases in levels of POPs in some biota. Although collectively understudied, multiple stressors, including contaminants and climate change, may act to cumulatively impact some populations of Arctic biota. Models are useful for predicting the net result of various contrasting climate-driven processes on POP and CEAC exposures; however, for some parameters, especially food web changes, insufficient data exists with which to populate such models. In addition to the impact of global regulations on POP levels in Arctic biota, this review demonstrates that there are various direct and indirect mechanisms by which climate change can influence contaminant exposure, accumulation, and effects; therefore, it is important to attribute POP variations to the actual contributing factors to inform future regulations and policies. To do so, a broad range of habitats, species, and processes must be considered for a thorough understanding and interpretation of the consequences to the distribution, accumulation, and effects of environmental contaminants. Given the complex interactions between climate change, contaminants, and ecosystems, it is important to plan for long-term, integrated pan-Arctic monitoring of key biota and ecosystems, and to collect ancillary data, including information on climate-related parameters, local meteorology, ecology, and physiology, and when possible, behavior, when carrying out research on POPs and CEACs in biota and food webs of the Arctic.
Environmental significanceAs it affects physical, biological, and ecological processes in the environment, global climate change has the potential to influence the uptake and fate of persistent organic pollutants (POPs) and contaminants of emerging Arctic concern (CEACs) in biota and food webs through multiple mechanisms. This review paper summarizes the current state of knowledge on the influence of both climate change-modulated physical environmental processes and ecological processes on the accumulation of POPs and CEACs in biota and food webs within Arctic marine, terrestrial, and freshwater ecosystems. It also covers climate change impacts on contaminant toxicity, and modeling and research tools to improve understanding of interactions between climate change and contaminants in biota and food webs. This study reveals that climate change is impacting the long-range transport of pollutants to the Arctic and within the Arctic by altering atmospheric, environmental, and ecological processes and, thus, is influencing the exposure and accumulation of POPs in Arctic wildlife. |
The potential effects of climate change on the food web accumulation of contaminants in polar regions was first identified just after the turn of the 21st century.1 As summarized in UNEP/AMAP (2011),3 the predicted effects of climate change on contaminant bioavailability and subsequent uptake and bioaccumulation were investigated based on various modeling studies constrained by Intergovernmental Panel on Climate Change (IPCC) forecasts. Prior to that, the focus was predominantly on understanding how the unique characteristics of polar ecosystems influence pollutant uptake and accumulation.4,5 In the decade that followed, national and international efforts, including the International Polar Year, obtained empirical data on the effects of climate change on POP food web accumulation in the Arctic, as summarized for marine ecosystems by McKinney et al.6
The bioavailability, uptake, bioaccumulation, and fate of POPs in Arctic ecosystems are affected by a multitude of complex and interacting factors, including environmental conditions, such as temperature, precipitation, and presence of sea ice; POP physicochemical properties, such as hydrophobicity and recalcitrance; and a number of biological factors, including physiological and ecological variables such as phenology, energy allocation and lipid dynamics, reproductive strategy, body size, age, sex, life-cycle stage, biotransformation capacity, habitat use, migration and feeding ecology (reviewed in Borgå et al.).7 POPs taken up by Arctic biota may therefore be directly impacted by physical changes in the environment, such as changes in sea ice extent, thickness, age (i.e., multi-year vs. first-year ice), length of the ice season in seasonal ice regions, presence of glaciers, snow cover and permafrost, as well as increased freshwater runoff and altered nutrient availability.8 Additionally, biological and ecological changes, including increased primary production, reduced population sizes of some ice-dependent species, northward range shifts of sub-Arctic and temperate marine and terrestrial species, and altered trophic structuring can influence POP dynamics in Arctic food webs.9–12 Climate change is projected to concurrently alter many of these processes, potentially resulting in decreased or increased contaminant exposure at the base of the food web and changes to food web accumulation (Fig. 1 and S1†).
Here, we review the current state of knowledge regarding climate change-driven physical and ecological factors impacting contaminant fate and accumulation in Arctic species and food webs (Table S1†), with a focus on the POPs and CEACs as defined in de Wit et al.13 We summarize findings from available studies published over the past decade since the previous UNEP/AMAP report3 to address two overarching questions: (1) how do climate change-driven physical processes influence POPs in Arctic wildlife? and (2) How do climate change-driven ecological processes influence the levels of POPs in Arctic biota and food webs? We also highlight how climate change may interact with other stressors to impact contaminant toxicity, as well as the utility of modeling and newer research tools in closing gaps in our understanding of climate change-contaminant interactions. The processes described herein are connected to those reviewed in the other papers in this issue on the abiotic environment14 and temporal trends in biota.15 Finally, we provide conclusions on our current understanding of climate change influences on POPs in Arctic species and food webs, identify knowledge gaps, and provide recommendations for future research and monitoring. Research findings for each topic are reviewed in order, from low- to high-trophic levels for marine, terrestrial, and freshwater ecosystems, when studies were available, and include results from empirical and model-based studies.
To date, research on the direct relationship between temperature and POP uptake and elimination rates in Arctic biota is still scarce.17 Uptake and elimination are both likely to increase with warming temperatures due increases in metabolism (e.g., ref. 18); however, this relationship may be limited by species-specific biotransformation capacities related to physiological differences (inclusive of thermal tolerances), as well as variations in life history strategies (i.e., body size, longevity, age at maturity, and fecundity and age). In one study, the bioaccumulation potential of POP-like and hypothetical chemicals modeled in the food chain of a non-Arctic fish species, the round goby (Apollonia melanostomus), projected only a minimal impact of climate change on temperature-dependent uptake and elimination rates.19 Similarly, a modeling study of the Barents Sea food web concluded that annual temperature increases would elevate some processes and decrease others, to have only limited effects on net bioaccumulation.20
In addition to its effects on chemical uptake and elimination rates, temperature also affects the bioavailability of, and thus exposure to, contaminants in the Arctic. Warmer air temperatures can increase mobilization and long-range transport of POPs from primary emissions21,22 and re-mobilization from secondary sources, such as soil, melting permafrost, seawater, glacial ice, and sea ice.14,23,24 For example, sediment from Resolute Lake in the high Canadian Arctic might act as a net secondary source of polychlorinated biphenyls (PCBs) to the atmosphere.25 Additionally, long-term temporal trends of PCBs in Arctic air measured at Alert (Canada), Stórhöfði (Iceland), Zeppelin (Svalbard, Norway) and Pallas (Finland), between 1993 and 2012 showed a lower rate of decline in recent years and an increase of PCBs in air after the year 2000.26 This increase may be connected to re-emissions of PCBs from secondary sources due to higher ambient temperatures. The changing of environmental contaminant sinks to contaminant sources could have major effects on the exposure and thus bioaccumulation potential of POPs in wildlife. Nonetheless, it is important to note that both observed and predicted changes in the environmental fate of chemicals strongly depend on compound-specific physical–chemical properties, such as volatility and water solubility.27
Modeling and correlative studies strongly suggest that climate change affects the uptake and accumulation of POPs in Arctic wildlife through altered wind and precipitation impacts on contaminant transport or wet deposition.22 Correlations have been documented between the AO or NAO indices and contaminant trends in glaucous gulls (Larus hyperboreus) in the Norwegian Arctic and ringed seals (Pusa hispida) in the Greenlandic and Canadian Arctic.27,29,30 Glaucous gulls from Bjørnøya, Svalbard, showed positive relationships between POP concentrations and AO conditions in the preceding summer and preceding winter.29 Concentrations of POPs in ringed seals from West Greenland were also positively associated with winter AO conditions.31 Similarly, POP concentrations in ringed seals from the Canadian Arctic Archipelago and Hudson Bay were positively related to AO and NAO states, respectively, in the preceding year.30 These findings indicate that POP concentrations were higher in gulls and seals following years with greater transport of air masses from North America and Europe toward the Arctic. In glaucous gulls, negative relationships between concentrations of POPs and AO in the current winter were suggested to be related to changes in diet or overwintering areas.29 The cold winters associated with a negative AO phase (AO−) may decrease the availability of high trophic level prey during the following summer; alternatively, cold winters may force glaucous gulls to migrate further south, exposing them to higher levels of POPs.
Climate-related increases in precipitation were previously predicted to increase the deposition of POPs through scavenging.1 Consistent with this prediction, organochlorine contaminant concentrations in thick-billed murre (Uria lomvia) and northern fulmar (Fulmarus glacialis) eggs in the Canadian Arctic were related to the NAO and/or rainfall amounts between 1975 and 2014 after controlling for differences in chemical partitioning characteristics, feeding ecology, and chemical emissions, which explained the majority (70% in murres; 77% in fulmars) of variability in egg organochlorine contaminant concentrations.32 More specifically, years of higher rainfall were followed by years with thick-billed murre eggs having higher concentrations of chlorobenzenes, octachlorostyrene, dieldrin, DDE and most PCBs, and lower concentrations of heptachlor epoxide, oxychlordane, cis-nonachlor, PCB-170 and PCB-180. In fulmars, years experiencing positive NAO phase (NAO+) conditions, when industrial chemicals are transported to the Arctic from eastern North America and Eurasia, were followed by years when eggs contained higher concentrations of chlorobenzenes, trans- and cis-nonachlor, dieldrin, photomirex and mirex. In contrast, negative NAO phase (NAO−) conditions were associated with birds laying eggs with lower concentrations of these chemicals. A four- to nine-year time lag occurred between the occurrence of the climate pattern and the apparent influence on POP concentrations deposited in the bird eggs.32 POP concentrations in landlocked Arctic char (Salvelinus alpinus) from the Canadian Arctic were also positively correlated with interannual variations of the NAO.25
Changes in the movement of water masses (e.g., oceanic currents) can also influence contaminant transport to, and within, the Arctic, and thus affect exposure to biota. From 1994 until 2010, some POPs, including PCB-52, PCB-153, p,p′-dichlorodiphenyldichloroethylene (p,p′-DDE), hexachlorobenzene (HCB), and α- and β-HCHs decreased in the blubber of ringed seals from central West Greenland.31 β-HCH, a chemical known to undergo long-range oceanic transport, showed the slowest decline, which was correlated to changes in salinity resulting from the dominance of the certain currents around Greenland. That is, years dominated by warm, saline water experienced a greater influx of the Atlantic water branch, around West Greenland, and potentially with it, a higher influx of ocean-transported POPs. Similarly, an increased flux of warm, saline Atlantic water into the European Arctic33 will likely increase oceanic contaminant transport into the high Arctic. Passive samplers deployed at different depths in the Fram Strait, an important deep-water channel to the Arctic Ocean, revealed a net influx of 0.16 megagrams of PCBs per year from the Atlantic into the Arctic, and a net export of HCB and HCH from the Arctic into the Atlantic Ocean, highlighting the continuous omnipresence of these chemicals.34
Several studies have shown relationships between sea ice parameters and POP concentrations in marine mammals (Table S1†), which may be related to various processes. Studies on ringed seals from several areas of the Canadian Arctic have reported significant relationships between concentrations of PCBs and/or chlorinated pesticides and sea ice coverage.30 The results mostly suggest that contaminants accumulate in ringed seals to a higher extent in years with greater total sea ice coverage and/or years with more multi-year sea ice (MYI). This may be related to the ability of sea ice to facilitate the delivery of organic contaminants to Arctic marine food webs,36 as sea ice, along with snowpack and glaciers, are reservoirs for organic contaminants. In addition, re-mobilization of contaminants from sea ice and water reservoirs is increasing due to increasing temperatures and decline of sea ice.23 However, changes in sea ice phenology also affect the prey composition of ringed seals.30 In contrast to the positive correlations between contaminants and total and MYI sea ice coverage, first-year sea ice (FYI) coverage was negatively correlated with POP levels in ringed seal blubber; therefore, the authors suggested that accumulation of contaminants in the cryosphere may also depend on the age of sea ice.30 Higher levels of POPs have also been associated with shorter sea ice seasons in ringed seals from the Beaufort Sea and Greenland, possibly related to variations in feeding opportunities.31,37 Although changes in prey composition were not evaluated in these studies, in both cases, the authors suggested POP changes could be related to prey type.
Studies from the Barents Sea have reported higher levels of lipophilic POPs and per- and polyfluoroalkyl substances (PFASs) in polar bears (Ursus maritimus) that use high-quality sea ice habitats (i.e., the marginal ice zone) in eastern Svalbard/Barents Sea, compared to bears using habitats with less or no sea ice in western Svalbard.38–41 These differences may be influenced by the presence/absence of sea ice,14 in addition to other biotic factors, which are discussed in Section C. The presence of sea ice is likely to increase uptake of atmospherically-deposited POPs into the marine food web. For example, PFASs deposited on surface snow are released and concentrated into the meltwater at the base of the snowpack by the end of the melting period,42 which is followed by phytoplankton bloom and subsequently large increases of zooplankton biomass.43 In contrast, PFASs deposited on open water surfaces are diluted by sea water and found at lower concentrations than in meltwater.44 Concentrations of lipophilic POPs in Barents Sea polar bears were also related to variations in the extent of sea ice on a seasonal scale.45 Plasma and adipose tissue concentrations of POPs in these polar bears were higher during seasons and in areas where sea ice conditions were poor. These relationships were, however, mainly related to changes in body fat, and are further discussed in Sections C.3 and C.4.
For Arctic foxes (Vulpes lagopus) in Svalbard, Norway, concentrations of β-HCH and several PFASs increased with increasing sea ice availability.38,46 However, these relationships were related to an increased availability of seals as further discussed in Section C.2.
In polar regions, runoff from terrestrial ecosystems includes meltwater from snow and glaciers, which can result in seasonal pulses of freshwater with contaminants to terrestrial, freshwater, and coastal ecosystems. Snow melt releases contaminants recently deposited from the atmosphere,53 as found for polycyclic aromatic hydrocarbons (PAHs) in Antarctica,54,55 whereas glacial melt transfers contaminants stored over the long-term to the terrestrial system, coastal waters, and lake sediments.56–58 However, in the Adventfjorden on Svalbard (Norway), which receives riverine inputs containing runoff from both a large tundra valley and glacier, benthic marine amphipods (Gammarus setosus) showed decreasing total PCB concentrations seasonally from April to August.59 The decrease in total PCB concentrations in amphipods coincided with an increase in terrestrial organic matter from riverine runoff, which may have both diluted and lowered the bioavailable fraction of contaminants at the base of the food web. Nonetheless, less hydrophobic PCBs peaked in amphipods in May/June, coinciding with snow melt. May and June dietary data (i.e., fatty acids) suggested pelagic diatoms made up a high proportion of the amphipod diets, and thus may be the link between the PCBs in snow melt and those detected in the amphipods.59 The net result of terrestrial runoff on contaminant accumulation in Arctic ecosystems is not known, and requires further investigation.60
A general decline in the duration of ice cover has been observed in Arctic lakes greater than 1 km2 in area. Satellite imagery surveyed for 13300 lakes across the circumpolar Arctic from 2000–2013 showed significant trends towards earlier ice break-up, ranging from −0.6 days per year in northern Alaska, to −0.1 days per year in northern Scandinavia.65 Lakes of this size generally have fish populations, with salmonids, and especially Arctic char, and are found largely above the treeline.66–68 These changes are likely to affect water temperatures and within-lake biogeochemical cycles, leading to increases in primary productivity and, consequently, changes in lake trophic relationships. Shifts in the diatom community structures of high Arctic lakes have been rationalized primarily by changing climatic factors such as shorter durations of ice coverage, longer growing seasons, and more favorable thermal stratification patterns.69
Climate warming is also thawing permafrost in lake catchments and shorelines resulting in changes to lake and river water chemistry. The thawing of permafrost results in inputs of dissolved organic carbon (DOC), inorganic solutes, including nutrients such as ammonium, nitrate, and phosphate and major ions such as Ca2+, Mg2+ and SO42−, and suspended solids to streams and lakes. Nonetheless, limnological properties are highly influenced by regional surroundings.64 Paleolimnological and water chemistry studies of retrogressive thaw slump-(RTS70) affected lakes located east of the Mackenzie River delta area of the Canadian Arctic have shown increases in the abundance and diversity of periphytic diatoms consistent with increased water clarity and subsequent development of aquatic macrophyte communities.71 There is broad agreement from several studies of lakes in this region that higher concentrations of Ca2+, Mg2+, SO42−, and other ions, lower DOC, and increased water clarity occur in RTS-impacted lakes compared to unimpacted lakes.72–74 The decline in DOC was attributed to the sedimentation of organics and particulates from the water column.74 In related work, reduced DOC in slump-affected lake water resulted in higher concentrations of POPs in sedimentary organic matter.75 Although the impacts of these changes on the bioaccumulation of POPs in these lake food webs have not been investigated in detail, amphipods (Gammarus sp.) from slump-affected lakes had higher mean ΣPCBs and ΣDDTs (27.5 ng g−1 lipid weight (lw) and 18.5 ng g−1 lw, respectively) than those from unaffected lakes (17.0 ng g−1 lw and 10.9 ng g−1 lw, respectively).76 Additionally, amphipod POP concentrations increased with the percentage of the catchment slumped.76 Mean ΣHCHs also differed between amphipods from slump-affected lakes (7.0 ng g−1 lw) and unaffected lakes (5.1 ng g−1 lw), although less so than for the more hydrophobic POPs. These water chemistry changes resulting from permafrost thawing and runoff might also affect the physiological condition of freshwater fish (Case study S1†). Other studies in northern Alaska (USA), northern Quebec (Canada), Siberia (Russia) and northern Sweden focused on tundra streams have shown elevated suspended sediments and DOC downstream of thermokarst activity.64,77
Climate change might affect exposure, accumulation, and effects of POPs on wildlife through shifts in the seasonal timing of their biological responses (e.g., migration, availability of prey, reproduction). In multiple seabird species in the European Arctic, plasma or body POP levels measured in parents and their offspring at summer Arctic breeding grounds reflected their exposures to higher contaminant concentrations on the overwintering areas in industrialized and agricultural regions.79–81 Changes in migration timing and the relative time spent at winter and breeding grounds could therefore affect POP exposure. Although links to POPs have not been examined, climate-associated changes in the timing of migrations by Pacific Arctic beluga and bowhead whales have also been documented.82,83 Due to climate change, the Arctic will experience increasingly earlier springs, with earlier availability of breeding grounds and food.9 If migrating birds do not respond to this earlier onset of spring, they may experience an increasing mismatch between the timing of resource availability and their arrival at breeding grounds, eventually resulting in poorer condition and reduced fitness (e.g., ref. 84). Birds may respond to this shift in phenology by migrating north earlier, which could also affect their migration route, timing, and energy used for reproduction at breeding sites (i.e., income vs. capital breeders), and thus, the POPs accumulated and transferred to offspring. Animals that rely mostly on stored reserves for breeding purposes (e.g., egg formation) are considered “capital breeders”, and species relying mostly on concurrent intake for reproduction are “income breeders”; there is a continuum across species in using stored nutrients for breeding.85 In the terrestrial barnacle goose (Branta leucopsis), PCB levels were slightly increased in later-hatching eggs, suggesting females arriving later to Arctic breeding grounds spend more time consuming vegetation at more-contaminated overwintering grounds, compared to earlier-arriving females that predominantly rely on less-contaminated vegetation at the breeding grounds to accumulate the energy reserves for reproduction and egg-laying.86 Based on this example, animals arriving at breeding grounds earlier and relying on the resources consumed there for reproduction, might be expected to have offspring with lower contaminant concentrations.
In the marine ecosystem, the seasonal melting of marginal ice zones influences the timing of algal blooms. The earlier reduction in sea ice extent and volume can lead to a mismatch in the timing of energy availability for secondary producers such as copepods, a key species in the Arctic marine food web.43 This may in turn alter exposures to POPs in these consumers and their predators. Seasonal changes in a pelagic food web can offer insight into climate change impacts in this regard. In Svalbard (Norway), an effort to elucidate the effects of climate change on contaminant bioaccumulation and biomagnification in a pelagic food web was conducted between 2007–2008 using a ‘space-for-time’ approach, in which two or more contrasting sites embodying characteristics representative of present-day and expected future conditions are used to model potential changes over time. Contaminant transfer in the pelagic food web of two fjord systems, Kongsfjorden (Atlantic fjord, representative of future Arctic climate conditions) and Liefdefjorden (Arctic fjord, representative of present-day Arctic climate conditions), were studied with regards to seasonal changes in POP food web bioaccumulation (only Kongsfjorden) and differences between POP bioaccumulation in the food webs of the two fjord systems.87–89 At Kongsfjorden, seasonal differences in POP levels were observed in zooplankton (decrease from May to October), various fishes (different trends depending on species), little auks (Alle alle; similar trends as for zooplankton), and kittiwakes (Rissa tridactyla; different trends depending on contaminant).88,89 Given the diversity of seasonal changes in POPs experienced by different trophic levels and species of the fjord food web, knowledge of seasonal variation will be essential for differentiating climate change-induced alterations from seasonal changes in POP bioaccumulation and biomagnification.
In those studies, differences between the Atlantic- and Arctic-influenced fjords were observed both with regards to POP transfer and fjord-specific characteristics. Higher POP concentrations and bioaccumulation factors (BAFs) were generally observed in zooplankton species from Kongsfjorden. However, there were compound-specific differences, with Liefdefjorden having higher concentrations and proportions of the more volatile compounds, such as α-HCH and chlordanes.87 These differences were attributed to different water masses, ice cover, and freshwater input: all factors which have clear climate relevance but may represent confounding factors associated with the fjords. The seasonal input of freshwater due to the melting of glaciers and snow was earlier in southern Kongsfjorden than northern Liefdefjorden. Further, Kongsfjorden did not have sea ice during winter, while the sea ice in Liefdefjorden broke up shortly before sampling. These differences made the two fjords out of phase with each other with regard to the release of volatile POPs from sea to air. Delayed volatilization of compounds in Liefdefjorden due to the later melting of snow, glacial ice, and sea ice resulted in higher concentrations/exposure to volatile POPs in seawater in Liefdefjorden compared to Kongsfjorden.88,90,91 Therefore, the differences observed in zooplankton and the food webs of the two fjords could be due to seasonal variation and not climate change.88 Overall, there seemed to be a seasonal lag in the peak concentrations detected in biota beginning with an early season peak in zooplankton and progressing to mid- and late-season peaks in fish and birds. Environmental changes in contaminant exposure appear to be reflected faster in zooplankton than in longer-lived, higher trophic-level organisms.92 Thus, extrapolating the results of this study to approximate future climate change effects should be done with caution. Space-for-time approaches require two contrasting sites that are comparable with regards to other potential confounding factors within- or between-sites. This was not achieved within these two fjord systems. Additional investigations with sampling taking place over several years, and/or using Arctic and Atlantic locations outside the fjords that are not influenced by sea ice and meltwater, could be more suitable for elucidating climate change influenced alterations in bioaccumulation and biomagnification.
In consumers at lower trophic levels of marine and freshwater food webs, POP exposure occurs both through direct uptake from water and through dietary exposure.87,96–98 Concentrations of POPs in Arctic marine zooplankton from Lancaster Sound in the Canadian Arctic decreased in the open water period, being lowest in August–September.96 This was due in part to reduced water and algae organochlorine contaminant concentrations, but may also have been related to increasing zooplankton lipid content during this time. Similar trends were observed in Arctic zooplankton from Kongsfjorden, Svalbard, with overall POP levels decreasing from May to October.87 In the western Canadian Arctic, burbot (Lota lota) from the McKenzie River showed increasing PCB concentrations from 1985–2009, despite decreasing or stable atmospheric concentrations over the same time period.99 The authors showed evidence that the rising PCB concentrations in these freshwater fish were likely related to increases in primary production and ‘algal contaminant scavenging’ into the food web.99
Increased temperature and ocean current velocities due to climate change are leading to changes in the northern boundary of various species, with Arctic species generally following the receding marginal ice zone northward, and boreal species increasing their presence and abundance in the Arctic, a phenomenon referred to as ‘borealization’.11,105 Borealization will lead to changes in the species composition of Arctic food webs, and thus changes in trophic interactions, and ultimately alterations in biomagnification. This concept is referred to as ‘Atlantification’ in the marine environment of the European Arctic and ‘Pacification’ in the western North American Arctic. Altered community composition gives rise to novel interactions among species, such as new competitive and predatory interactions. Disturbed predator-prey interactions, in particular, are expected to have outsized consequences for ecological communities through trophic cascades.106 These changing vertical interactions can be anticipated to impact POP burdens in biota, since most exposures in consumers occur largely via diet.5,87,107 Moreover, given that the greatest rates of species invasions, at least for marine ecosystems, are projected to occur in polar regions,102 and rapid Arctic borealization has already been documented (e.g., in the Barents Sea11), we can speculate that polar marine wildlife may be more at risk from climate-driven ecological changes and POP exposures, than other regions.108
With the current lack of understanding of how climate affects food web structure and the strengths of interspecific interactions for most species, the ability to forecast shifting biotic interactions under climate change remains limited.109 Hence, models exploring the future impacts of shifting food webs on POPs are non-existent. However, in the past decade, several studies have empirically examined the role of changing trophic relationships on tissue POP concentrations in Arctic species and food webs. The current state of knowledge is based on research from Canadian, Alaskan (USA), Greenlandic (Denmark), and Norwegian Arctic ecosystems, however, data from the Russian Arctic are lacking. Thus far, most published research on the impacts of shifting species interactions on POPs has focused on marine environments, with marine mammal and seabird studies predominating, but a few studies on marine fishes and a terrestrial mammal (Arctic foxes) have also been performed; findings from these studies are summarized below.
In three regions of varying latitudes within the eastern Canadian Arctic, concentrations of PCBs, organochlorine pesticides (OCPs), and various flame retardants were compared between marine prey fish communities, including native and redistributing non-native species.110 The comparisons particularly focused on the keystone species, Arctic cod (Boreogadus saida), and the sub-Arctic ‘replacement’ species, capelin (Mallotus villosus) and sand lance (Ammodytes sp.). Although not all species were collected in all regions, native benthic species, specifically sculpin (Cottoidea sp.) and northern shrimp (Pandalus borealis), showed the highest concentrations of ∑PCBs and ∑OCPs of all species. Nonetheless, capelin, and to a lesser extent, sand lance, showed higher lipid-normalized concentrations of many POPs than Arctic cod (Fig. 2). Furthermore, neither feeding habitat (based on δ13C values), trophic position (based on δ15N values), or fish length accounted for this difference. This finding, along with congener/compound pattern differences among the species (e.g., higher proportion of DDT vs. DDE in capelin relative to Arctic cod), suggested that the higher POP concentrations in invading sub-Arctic capelin may, at least in part, be a consequence of the fishes acquiring temperate contaminant signatures during seasonal migrations. Thus, these species may act as climate-driven contaminant biovectors into Arctic marine food webs.
Fig. 2 Mean (± SE) lipid weight (lw) concentrations of PCBs (left) and DDTs (right) in muscle of marine forage fish and invertebrates of the low-(LA), mid-(MA), and high-(HA) Canadian Arctic from 2012–2014. Arctic cod (green), the two sub-Arctic species, capelin, and sand lance (orange), sculpin (blue), and northern shrimp (red) are shown. Significant differences in contaminant concentrations among species are indicated by different letters above each bar. Note that the y-axes scales differ. Adapted from Pedro et al. (2017).110 |
Similarly, both transient marine fish and mammals showed higher POP concentrations than similar trophic-position Arctic-resident species in a study of the food web of Cumberland Sound, Canada.111 Elevated concentrations of several current-use pesticides (chlorothalonil, chlorpyrifos, dacthal, and endosulfan) were also found in the capelin of Cumberland Sound relative to other species at a similar trophic position.112 Nonetheless, the differences in POP concentrations between these native and non-native replacement fishes were small (two-fold or less), suggesting changes in these forage fish communities may have a limited impact on POP concentrations in piscivorous predators (i.e., predator fish, marine mammals, and seabirds). Direct impacts of these changes in prey fish on predator POP concentrations have been investigated for a few seabird species (as discussed further below), but not for any piscivorous fish or marine mammals to date.
Killer whales (Orcinus orca), perhaps the most contaminated species on Earth with respect to POPs,113 are undergoing a northward range-shift into Arctic and sub-Arctic waters in multiple circumpolar locations.114–116 Within the eastern Canadian Arctic and around southeast Greenland, killer whales appear to be targeting marine mammals for food, which may differ from their presumably more typical diet of fish in the North Atlantic.80,114,117–119 However, recent studies have suggested that some North Atlantic killer whale populations may also eat prey other than fish. Killer whales sampled between 2012–2014 around southeast Greenland showed mean blubber concentrations of PCBs (40 mg kg−1 lw) and OCPs (70 mg kg−1 lw) which were an order of magnitude higher than those found in killer whales sampled over similar years in the Faroe Islands.120 This substantial difference in contaminant burdens was suggested to be due, at least in part, to North Atlantic whales increasingly feeding on marine mammals in new Arctic habitats (Fig. 3).
Fig. 3 Mean concentrations (±SE) of ΣPCBs, ΣCHLs, and ΣDDTs (mg kg−1 lw) in blubber of eighteen killer whales sampled in 2012–2014 in southeast Greenland (left) and two killer whales sampled in 2008 in the Faroe Islands (right). Red horizontal lines indicate thresholds for risk (10 μg g−1 lw), and high risk (100 μg g−1 lw) of health effects associated with PCB exposure as reported by Dietz et al.190 Adapted from Pedro et al.120 |
Recent studies using stable isotopes as dietary markers in a killer whale population off the Norwegian coast also indicate that some whales may feed on marine mammals,121 resulting in elevated contaminant levels.122 Analysis of PCBs indicates 100% of the seal-eating killer whales were above the threshold for a risk of effects on immune and hormonal systems, compared to 43% of fish-eating killer whales (Fig. 4).123
Fig. 4 ΣPCB levels in blubber of seal-eating (n = 7) and fish-eating (n = 24) killer whales from Norway. (Left) Box plot of ΣPCB levels in killer whales by dietary group. Horizontal lines represent the median concentration, whiskers represent the lower (first) and upper (fourth) quartiles, and dots represent individual whales, with points outside the box and whiskers being outliers. (Right) ΣPCB levels of individual whales, in decreasing order by dietary group. Red horizontal lines indicate thresholds for risk (10 μg g−1 lw), and high risk (100 μg g−1 lw) of health effects associated with PCB exposure as reported by Dietz et al.190 Adapted from Andvik et al.122 |
With climate change has come a shift in the dietary composition of seabirds in the Canadian Arctic124 and in East Greenland,125 that in turn, has likely contributed to changes in their exposure and accumulation of POPs. In the Canadian Arctic, and particularly the Hudson Bay-Hudson Strait region, the diet of thick-billed murres has been shifting since the 1970s due to changing sea ice conditions. The diets of these low-Arctic breeding murres have changed to contain less ice-associated, cold-water Arctic cod and more of the sub-Arctic species, capelin, and sand lance.126–128 Conversely, high-Arctic breeding murres have continued to consume predominantly Arctic cod.129,130 These dietary patterns may explain inter- and intra-species differences in POPs of Canadian Arctic thick-billed murres and northern fulmars measured in 2008.124 Similarly, temporal increases in hepatic organochlorine concentrations in murres from the high Arctic measured in 1998, 2003 and 2008, may be partially related to concurrent changes in diet and trophic position.124 In East Greenland between 2004 and 2015, the diet of little auks (Alle alle) contained less lipid-rich, ice-associated prey when sea ice coverage decreased, with adult birds consuming zooplankton at a higher trophic position, as reflected by an increase in blood δ15N values.125 Concurrently, feather mercury (Hg) concentrations of these auks increased by 3.4% each year, notably occurring during the summer and not when overwintering.125,131 Surprisingly, additional research suggests that interannual variation in little auk feather Hg concentrations reflect changes in food contaminant levels, notably increased Hg levels in zooplankton during summertime, and not the reorganization of the food web or modification of seabird trophic ecology, based on δ15N and δ13C values in little auk whole blood.131 Similar patterns may occur for lipid-soluble POPs. For further studies of Hg, we refer the reader to the recent AMAP report on Hg in the Arctic in relation to climate change.132
Climate change is altering the structure of some Arctic food webs and influencing trophic positions of predatory seabirds and potentially influencing long-term trends in environmental contaminants. After adjusting for the reduced trophic position of thick-bil led murres in Hudson Bay, Canada (based on egg δ15N values), declining trends of organochlorine concentrations (HCB, heptachlor epoxide, oxychlordane, dieldrin, p,p′-DDE and ∑PCB69) measured in eggs between 1993–2003 were more moderate than originally thought.133 In contrast, adjusting for the increased trophic position of high-Arctic breeding murres, resulted in increased rates of organochlorine concentration declines, although p,p′-DDE and ∑PCB69 remained relatively unchanged. A direct link between climate change and contaminant trends was not investigated in this study.
Also, in the Canadian sub-Arctic region of western Hudson Bay, dietary and/or food web shifts were linked to changes in POP trends in polar bears. The diet change was associated with faster increases in adipose concentrations of total polybrominated diphenyl ethers (∑PBDEs) over time, and with a switch from slowly decreasing to slowly increasing concentrations of ∑PCBs from 1991–2007.134 Based on the decline of δ13C values and particular shift in fatty acid profiles observed, the authors speculated that the change in contaminant trends may have been related to changes in the proportion of bearded seal (Erignathus barbatus), sub-Arctic harbor seal (Phoca vitulina) and harp seal (Pagophilus groenlandicus) in the bears' diet, although other food web changes could also have contributed. A similar change was observed in East Greenland polar bears between 1984–2011; a shift towards feeding on sub-Arctic hooded seal (Cystophora cristata) and/or harp seal was directly estimated based on quantitative fatty acid signature analysis, and was associated with NAO conditions.135 However, at least over the time period studied, the shift did not significantly alter trends in adipose concentrations of any of the POPs examined (PCBs, OCPs, and PBDEs) (Fig. 5).135 A re-evaluation of this relationship including recent years through 2016 in East Greenland polar bears, suggested that hooded and/or harp seals now represent the majority of the polar bears' diet, and that a concurrent rise in concentrations of certain POPs, including PCBs, occurred at least until 2013.136 The effect of this ecological change on POPs concentrations may be due to harp seals, and particularly, hooded seals occupying a higher trophic position than the bears' historical main prey, ringed seals, and/or due to these migratory sub-Arctic seals acting as biovectors transporting contaminants from further south into Arctic marine food webs.111
Fig. 5 Dietary trends and POP concentrations in East Greenland polar bears. (Top) Contribution of ringed seal and hooded seal to polar bear diets in spring-summer from 1984–2011, based on estimates from quantitative fatty acid signature analysis. Open circles represent individual bears, while black circles represent annual means. Non-significant trends for bearded seal, harp seal, narwhal, and walrus not shown. (Bottom) Annual change of POP concentrations (mean ± SE) in these polar bears, both unadjusted and adjusted (controlled) for changes in percentage of ringed seal in the diet over this same time period. Note that controlling for dietary shift did not significantly alter trends Adapted from McKinney et al.135 |
In the southern Beaufort Sea, some individual polar bears are spending extended periods of time onshore during the reduced ice season and consuming onshore foods, namely the remains of subsistence-harvested bowhead whales (Balaena mysticetus).137,138 These onshore polar bears have shown lower concentrations of ∑CHLs, but not other legacy POPs, compared to individuals remaining on the sea ice year-round, likely because bowhead whales occupy a lower trophic position than the bears' main ice-associated prey, ringed seals.139 Recent observations of polar bears foraging on terrestrial species, including colonial nesting seabirds, have also been reported in the Beaufort Sea and for other subpopulations of polar bears,140–142 which all else being equal, would serve to lower their POP exposures.
Temporal trends of lipophilic POPs and PFASs in relation to changes in feeding habits and body condition have also been investigated in adult female polar bears sampled in spring from the Barents Sea.38,143 Temporal changes in δ15N and δ13C values indicate that the winter diet of Barents Sea female polar bears has shifted towards less marine, less ice-associated and lower trophic level prey items,38,143 which is likely related to the rapid decline of sea ice in the area.144 Although not statistically significant, the same polar bears also tended to get thinner over the period 1997–2006, whereas from 2006–2017 they got fatter.143 Contaminant concentrations were strongly related to polar bear diet (lipophilic POPs and PFASs) and body condition (lipophilic POPs only), however, temporal changes in diet and body condition did not significantly affect trends of the measured pollutants.38,143
Climate change responses also affect POPs in Arctic terrestrial communities. The Arctic fox is an opportunistic predator and scavenger that is among the most contaminated of Arctic mammals. Although a terrestrial species, it also feeds on marine foods at coastal tundra areas, such as Svalbard (Norway). Studies investigating the influence of feeding habits and food availability on contaminant temporal trends in Arctic foxes from Svalbard reported that concentrations of all measured lipophilic POPs and PFASs increased with a higher intake of marine food items. Additionally, concentrations of β-HCH, perfluorooctane sulfonic acid (PFOS), perfluoroheptanesulfonic acid (PFHpS) and perfluorotridecanoic acid (PFTrDA) were positively associated with sea ice availability.38,46 A similar, but non-significant tendency was also observed for p,p′-DDE, HCB, CHLs, PBDEs and for perfluorohexane sulfonic acid (PFHxS) and the C9–C11 perfluorocarboxylates (PFCAs). Arctic foxes use sea ice to scavenge the remains of seals killed by polar bears and to hunt newborn ringed seal pups.145 The increased availability of sea ice thus increased Arctic foxes' access to seals, which have higher levels of pollutants than terrestrial prey. Additionally, decreasing concentrations of HCB and PFASs in fox tissues were related to the increasing mortality of reindeer (Rangifer tarandus platyrhynchus).38,46 Reindeer mortality has been connected to rain-on-snow events,146 which entirely encapsulate short-growing vegetation in ice across large areas of high Arctic tundra.147 The frequency of these extreme weather events, such as warm spells and rain-on-snow events, has been increasing in Arctic terrestrial regions.148 Icing of the tundra prevents winter foraging opportunities for Arctic tundra herbivores, which affects their vital rates and has further consequences for the abundance of Arctic carnivores.146 In conclusion, climate change may decrease POP exposures in Arctic foxes as (1) decreases in sea ice will reduce opportunities for Arctic foxes to eat more-contaminated seals and (2) increased numbers of reindeer carcasses become available due to more frequent rain-on-snow events, providing foxes with less contaminated reindeer as a food source.
Behavioral changes and alternative uses of landscapes and food resources may, however, buffer the impacts of environmental change in some species.149,150 For example, the proportion of Svalbard reindeer feeding along the shoreline has increased along with icier winters, and stable isotope values indicate that these reindeer are likely feeding on washed-ashore kelp.150 Although not yet investigated, this may have further implications for contaminant intake of reindeer and their predators. In addition to increasing reindeer mortality, rain-on-snow events also promote the expansion of boreal predators into Arctic tundra. In Yamal, Russia, a high mortality event of domestic reindeer due to icing of the snow layer was followed by the increased presence of red foxes (Vulpes vulpes) and hooded crows (Corvus cornix).151 Other changes, such as the ‘shrubification’ of the Arctic tundra, are also underway,152,153 and like their marine counterparts, terrestrial fauna are also undergoing northward range shifts.9 These shifts in terrestrial plant and animal communities have the potential to alter both contaminant movement and fate within terrestrial Arctic food webs.
Neutral storage lipids represent the main organic phase that non-ionic chemicals, such as most POPs, partition into.154 The relationship between lipid content and POPs has rarely been studied over time in organisms with very high seasonal variation in lipid content. Arctic marine zooplankton adapted to fasting during winter hibernation, such as Calanoid copepods, store energy as energy-dense wax esters,155–157 and not only as triacylglycerols, which act as the dominant neutral storage lipid class in vertebrates and winter-active zooplankton. As POPs partition into wax esters more than into triacylglycerols,158 a change in the overwintering strategy of Arctic zooplankton to include more active feeding, combined with increased abundances of winter-active southern zooplankton in the Arctic, will change (and likely decrease) POP accumulation at the base of the aquatic food web.
Borealization is also expected to contribute to changes in lipid-dynamics within food webs by altering low- and mid-trophic species compositions. As the climate warms, the Arctic will likely experience an increasing abundance of species with lower lipid content, smaller size, faster growth, and an adaptation for income breeding rather than capital reproduction based on stored energy (Fig. 6).11,159,160 As a result, the lipid mass at the base of the Arctic food web is predicted to decrease with a greater influx of southern secondary producers such as Calanus finmarchicus, and reduced influence of large, lipid-rich Arctic species, such as Calanus glacialis and Calanus hyperboreus.160,161 However, at an ecosystem level, the overall lipid mass is expected to increase due to a more rapid turnover of copepods and a higher general biomass of zooplankton.160
Fig. 6 Projected effects of borealization on Arctic marine zooplankton and associated changes in organism POP concentrations. Modified from Renaud et al.160 |
A shift from capital to income breeding (i.e., being less dependent on stored lipid reserves) may in turn result in lower maternal transfer of POPs to offspring due to growth dilution, and lower POP storage capacity due to lower lipid content (see Section B.6). However, a meta-analysis of maternal transfer in marine mammals and seabirds did not show any relationship between reproductive strategy and maternal transfer of POPs to offspring for existing mother-pup data in pinnipeds162 and maternal investment in birds.86 These studies also identified significant knowledge gaps as most results to date come from capital breeders and not income breeders (Fig. 7).
Fig. 7 Relationship between reproductive strategy and maternal transfer of mercury (Hg) and POPs, including lipid-soluble PCBs and protein-associated PFOS, in marine mammals (upper panel) and seabirds (lower panel). Upper panel: Association of the estimated capital breeding index value and mother-to-pup contaminant transfer for various pinniped species. Response ratios indicate the ratio between concentration in offspring (Cpup) and mother (Cmother), with values above one indicating higher concentrations in offspring, and values below one indicating higher levels in mother. Gray areas depict confidence intervals. Capital breeding index values indicate the species’ relative reproductive strategy (i.e., the degree to which the pup is dependent on income (newly acquired) or capital (stored body lipids)) energy from the mother, with higher values indicating greater reliance on capital breeding strategy. Lower panel: Association of the seabirds’ estimated clutch-to-female mass ratio and the ratio between concentration in adult males and females, with values above one indicating higher concentrations in males, and values below one indicating higher levels in females (mothers). Clutch-to-female mass ratio represents the energy invested by the mother into the offspring. Figure from Hitchcock et al.86,162 |
The common eider (Somateria mollissima) is a large-bodied marine duck which undergoes large seasonal changes in lipid mass. Prior to breeding, females accumulate large lipid reserves and may subsequently lose up to 45% of their body mass during the egg laying and incubation fasting period due to lipid metabolism.163 One study investigated blood concentrations of PCB-153, p,p′-DDE and HCB in incubating common eiders in Svalbard and the sub-Arctic northern Norwegian mainland.164 The authors found that as a result of their higher metabolism, Arctic birds experienced higher blood-level exposures to POPs during the incubation period than sub-Arctic birds, as their use of energy reserves led to a greater re-mobilization of lipid stores, and thus movement of associated lipophilic POPs from fat into the blood. Moreover, increases in blood POP levels were more pronounced in females with low body condition compared to females with high body condition. Based on these findings, increasing temperature in the Arctic can be expected to lead to reduced metabolism and mobilization of lipids in breeding birds and consequently lower blood-level POP exposure during breeding.164
Arctic marine mammals also undergo seasonal changes in body fatness. For example, polar bears generally feed extensively in spring and early summer when ringed seals are pupping and molting on the sea ice.165 During sea ice-free periods, polar bears have reduced, if any, access to food, and therefore largely fast. Pregnant females spend winters in dens, fasting up to eight months in a row.165 The dynamics of lipophilic pollutants in polar bears are tightly connected to annual changes in lipid accumulation, as well as overall body condition, which has been documented to have decreased over time in several polar bear management areas.166–169 For example, studies from the Barents Sea have shown that concentrations of PCBs are several times higher in both plasma and adipose tissue of thin female polar bears compared to fat ones, and that body condition is a more important predictor for lipophilic POP concentrations than diet as inferred from stable isotopes.45 In seasons and areas with reduced sea ice extent, polar bears were thinner, and consequently had higher tissue concentrations of lipophilic contaminants. Although inverse relationships between sea ice and body condition with levels of pollutants have been shown at seasonal and inter-annual scales, long-term changes in springtime polar bear body condition, which are like ly related to sea ice conditions, did not seem to affect POP trends in polar bears.143 Plasma concentrations of proteinophilic POPs (i.e., PFASs) were not related to changes in body fatness of polar bears from the Barents Sea.38,39 However, PFAS concentrations were higher in fasting polar bears than in non-fasting bears, which may be related to higher protein concentrations or reduced metabolic rates, and thus lower contaminant excretion in fasting animals.39 Concentrations of the proteinophilic Hg (measured as total Hg) were inversely related to the body mass index (BMI) of polar bears from the southern Beaufort Sea.170
In the freshwater environment, wet weight PCB concentrations and lipid content in Arctic char from a lake in southwest Greenland decreased from 1994–2008 and were negatively correlated with summer air temperatures; however, the relationship between lipid-normalized PCB concentrations and air temperatures was less clear and did not show a significant trend.171 As there was no indication of changes in trophic morph or trophic level, air temperature was suggested as the influential factor impacting PCB levels in the fish by indirectly altering processes that resulted in reduced condition (i.e., lower lipid content) of the fish.
Substantial changes in the sea ice habitat of polar bears have forced them to spend longer portions of the year in lower-quality habitats with reduced access to high-quality prey176–178 or to move longer distances with greater energy expenditure.177,179,180 A comparison of selected polar bears with similar body conditions in the Barents Sea, showed that offshore polar bears were exposed to higher concentrations of pollutants than coastal bears and that this was related to differences in feeding habits, energy expenditure, and geographical distribution.41 Nonetheless, since offshore bears were, on average, fatter than coastal ones,181,182 plasma concentrations of lipophilic POPs were overall similar in bears with different strategies, and only proteinophilic PFASs were higher in pelagic bears.181 Compound-specific and bulk stable isotopes (δ15N and δ13C), home range, field metabolic rates (based on telemetry), as well as contaminant levels in harp seal prey from different locations, cumulatively indicated that higher POP concentrations in offshore Barents Sea polar bears were related to a combination of factors, including the consumption of greater proportions of high-trophic level and marine-based prey, higher levels of POPs in prey species, larger energy requirements, and distribution in marginal ice zones (Fig. 8).181,183 Conclusively, these studies indicate that the fate of pollutants in polar bears from the Barents Sea is the net effect of multiple factors that are driven by climatic conditions. If the costs of migration override the benefits of the energy-rich, high trophic-level prey in the future, offshore bears may exhibit higher concentrations of both lipophilic POPs and PFASs than coastal bears. Also, any change in movement strategies is likely to affect the fate of pollutants in polar bears.
Fig. 8 Combination of factors influencing differences in contaminant exposure and tissue concentrations of coastal and offshore polar bears from the Barents Sea. (A) Concentrations of lipid-soluble POPs are similar in fatter, offshore polar bears and in thinner, coastal bears from the Barents Sea, whereas PFAS concentrations are more elevated in the offshore bears.181 (B) When comparing offshore and coastal bears of similar condition, however, offshore polar bears are exposed to higher levels of POPs than coastal polar bears due to differences in feeding habits, energy expenditure, and geographical distribution.183 Tracks of offshore polar bears shown in blue, and tracks of coastal polar bears shown in orange. |
Novel infectious diseases, changes in pathogen distribution, and an increasing presence of zoonotic pathogens are likely to occur with climate change.190 The physiological stress induced by increased temperature might impair the immune system in cold-blooded organisms such as fish, rendering them more susceptible to infection and disease.191 The rate of phocine distemper virus (PDV) infection in pinnipeds from the North Pacific Ocean was recently found to be significantly higher in years following the presence of an open water route along the northern coast of Russia.192 The authors concluded that reduced sea ice removed a barrier to pinniped movement, thus allowing viral transmission between the North Pacific and eastern North Atlantic Oceans. Although the role of contaminants on disease transmission or animal health was not investigated in this study, some POPs have been shown to suppress immune responses, specifically in marine mammals.193 Taken together, it is conceivable that the combined effects of altered environmental conditions (i.e., temperature-stress), contaminant exposure (i.e., immunosuppression) and ecological changes (e.g., ice loss leading to greater disease transmission), could have greater impacts to animal health than any one factor alone.
Adding another layer of complexity, organisms are simultaneously exposed to multiple types of environmental contaminants with varying toxicities, including but not limited to, the POPs, CEACs, Hg and other potentially unrecognized chemicals or substances. More recently, plastics, including microplastics and nanoplastics, have been identified as substances that may interact with contaminants, or regardless, affect organisms directly.194 Contaminants have been shown to bind to plastic particles with high affinity, and thus plastics may act as vectors transporting chemicals to the Arctic, within the Arctic, and facilitate their uptake by organisms. Moreover, plastics contain chemical additives, and therefore may serve a source of new contaminants to Arctic ecosystems. Plastics are also substrates for microbiota which may play a role in ecosystem functioning. However, as there are currently no studies addressing how climate change directly or indirectly interacts with microplastics to affect contaminants within Arctic food webs, this topic is not addressed further here. See Halsall et al.14 for additional details on microplastics in the abiotic environment.
The effects of climate change may also affect an individual's sensitivity to toxicants (i.e., climate-induced toxicant sensitivity)18,195,196 by altering the toxicokinetics (i.e., uptake, distribution, metabolism, and elimination of toxicants) and toxicodynamics (i.e., toxicant interactions with biological receptors and enzymes) of contaminants (see Section B.1). Such changes may also alter a species' sensitivity and thresholds for effects, including lowering effect thresholds such that impairments are experienced with lower levels of contaminant exposure, as seen in developing children.197
Most studies investigating environmental pollution in combination with other stressors usually focus on the impacts of climate-related environmental parameters (e.g., temperature increase or drought). Most, if not all species have an optimal range of temperatures for which they can maintain homeostasis, outside of which, their ability to function may be impaired. However, the effects of environmental temperature changes are especially pronounced in ectothermic organisms which heavily rely on external sources of heat to control body temperature. Ectotherms display a reaction norm (general performance curve) in response to temperature, whereby they have an optimal temperature window, as well as upper- and lower-critical limits. The response variable (i.e., performance) reflects important determinants of fitness, such as life history traits (e.g., age at maturity), that directly affect the vital rates of a life table (e.g., survival, fecundity, growth). The ability of an organism to cope with other stressors, such as toxicants, depends on where within this temperature range they are.198
Understanding the effects of multiple stressors, compared to a single stressor such as temperature, can be more challenging. Usually, exposures to multiple stressors are predicted to have synergistic outcomes, in other words, a combined effect that is greater than the sum of the two. However, most laboratory-based studies are designed such that synergistic effects are also the most likely outcome; for example, many studies investigate the effects resulting from two stressors under continuous exposure at the same time,199 when real exposure regimes most likely differ from this scenario. Organisms do not always experience constant and simultaneous exposure to all stressors; rather, exposures can be pulsed, or experienced out of phase with one another, depending on habitat, season, or other variables. Exposure to stressors such as climate change and toxicants may be on a small or large spatiotemporal scale, static or dynamic, and the timing of responses and their duration might differ depending on the stressor. Taking these more realistic exposure scenarios into account by recreating pulsed and/or staggered exposures and including an organisms' time for a compensatory physiological response to maintain homeostasis, might result in combined effects that are additive or antagonistic, rather than synergistic, and provide more ecologically-relevant outcomes. For example, recent non-Arctic studies suggest that pesticide toxicity increases with daily temperature fluctuations as illustrated by the increased mortality and growth rate of aquatic damselfly (Ischnura elegans) larvae after chlorpyrifos exposure under conditions of fluctuating daily temperatures, but not under constant temperatures.200 Chlorpyrifos exposure under fluctuating daily temperatures also resulted in reduced energy storage (fat content) and increased enzyme activity compared to exposures under constant temperature conditions. These findings demonstrate that natural variations and fluctuations can influence responses to multiple stressors, and therefore are important considerations for both laboratory- and field-based studies.
Correspondingly, the health effects of multiple stressors, including those related to climate change, are increasingly being identified as knowledge gaps for the Arctic.3,190 Multi-stressor laboratory experiments can aid in understanding the combined effects of climate change and pollution, however, experiments specific to Arctic biota are still scarce, as organisms can be difficult to obtain due to the inaccessibility, harsh conditions, and high costs of research in the region. Although relevant POP multi-stressor laboratory studies are currently lacking, studies on PAHs can provide insight, as these petroleum-related substances share some properties of POPs and therefore behave similarly in some species.201–203
Modeling was used to evaluate the effect of multiple stressors, including climate change, pollution, and egg predation, on population dynamics and viability in the sub-Arctic common eider duck.204 Eider ducks feed low in the food web, and generally have low contaminant levels, apart from periods in the breeding season, such as late in the egg incubation period, when the female has fasted for approximately 20 days. During incubation the use of lipids for energy results in the re-mobilization of contaminants from lipid stores and a corresponding increase of contaminants in blood.205 Thus, when body condition is poor (i.e., low lipid stores), circulating contaminant levels are high. Projections of eider duck population growth and abundances using the Leslie matrix population model showed that egg predation alone was a strong enough factor to cause population extinctions, whereas an increase in sea surface temperature and reduced clutch size due to elevated contaminant exposure alone did not.204 However, any combination of the three stressors reduced the population growth rate, and the effect of pollution on clutch size was increasingly negative when co-occurring with a warming climate and increased egg predation. Population viability was lowest when all stressors occurred simultaneously. For pollution alone to be a significant stressor on the population of eider ducks through effects on clutch size, the level of exposure needed to be unrealistically high, however, when co-occurring with increased predation pressure, the threshold level needed for pollution to cause extinctions was lower.204
Field-collected data on wildlife populations have also been used to investigate the impacts of multiple stressors. In the great skua (Stercorarius skua), a seabird that breeds in the North Atlantic from Iceland and the Shetland Islands, Scotland to Svalbard, Norway, organochlorine and organobromine contaminants had negative effects on their return rate (i.e., survival) and different sensitivities to contaminant exposure were observed among breeding colonies.206 Effects occurred at lower contaminant exposure levels in breeding colonies with birds in poorer body condition, probably because of food limitations. Complex correlations among circulating PFAS congeners, total triiodothyronine, adult body mass and hatch dates, were evident in thick-billed murres breeding in Hudson Bay (Canada), and the authors concluded that the interaction of PFAAs on thyroid activity, in conjunction with several indirect effects of climate change, may cause additional stress to the murres.207 The impact of multiple stressors on contaminant levels in free-ranging wildlife has also been shown for polar bears in the marginal ice zone of the Barents Sea.183
Several studies suggest that POP exposure has adverse effects on lipid metabolism, immune function, circulating thyroid hormones and neurochemistry of polar bears,208 and climate change is likely to affect the sensitivity of polar bears to these adverse effects of pollutants. A wide range of biomarkers for energy metabolism, specifically gene transcript levels in fat, physiological parameters, and metabolomic and lipidomic markers in plasma, were investigated for associations with pollutant exposure in female polar bears from Svalbard, Norway.209 Several biomarkers involved in lipid metabolism were related to POP exposure. Furthermore, the differences between the biomarker responses of low- and highly-polluted bears were more contrasted during a period with poor sea ice conditions (Fig. 9), suggesting that contaminant exposure and sea ice decline may have compound adverse effects on polar bears.
Fig. 9 Transcript levels of a key regulator of lipid metabolism, peroxisome proliferator-activated receptor gamma (PPARG), in adipose tissue of low- and highly-polluted female polar bears from the Barents Sea are more contrasted during a period with poor sea ice conditions. Data from Tartu et al.209 |
The decline of sea ice and increased pollutant exposure may also have additive or synergistic effects on polar bear susceptibility to disease. In vitro and correlative field studies indicate that the polar bear immune system is compromised by pollutant exposure, as well as by climate change. For example, polar bears from the southern Beaufort Sea staying for extended time periods on shore had heightened immune responses and were more exposed to pathogens than bears remaining on the sea ice year-round.139,210 Also, concentrations of a stress hormone, cortisol, were higher in polar bears with low body condition, indicating that thinner polar bears were more stressed, and thus potentially more prone to adverse effects of pollutants, than fatter bears.211
The many overlapping changes occurring in the Arctic will challenge the ability of ecological communities to adapt and maintain resilience.212–214 The interaction between toxicants, climate change and other environmental changes is complex and adds to the difficulty in predicting effects in wildlife; however, the impacts observed in the Arctic environment may be able to function as early warning signals for other regions.108
The published work on mechanistic modeling of food web POP accumulation in response to climate change over the past decade is scarce and does not always apply to Arctic or polar regions. Nevertheless, the task of evaluating how climate change will affect food web POP bioaccumulation can be divided into the efforts of (i) combining abiotic contaminant transport models with biotic energetics and bioaccumulation models, and (ii) isolating the climate dimension and its contribution to bioaccumulation. The combined environmental fate and transport model, CoZMo-POP,217 and the bioaccumulation model, ACC-Human to CoZMoMAN,218,219 are good examples of models combining abiotic transport, biotic energetics, and bioaccumulation models that are built on similar fugacity-based frameworks. To identify the influence of climate change on POPs using these or similar models, a climatic mode is included in the abiotic transport model to estimate changes to contaminant transport and exposure at the base of the food web (e.g., ref. 220). This modeling framework has been used to assess the potential implications of various climate change scenarios for long-term human exposure to POPs in the Arctic.221 When modeling the most important processes contributing to future Arctic human exposure to POPs (i.e., emissions, transport, abiotic vs. biotic changes), Armitage et al.221 concluded that climate-induced changes to food web biomagnification would be most influential. Details on climate-related impacts to the physical environment that would influence contaminant exposure are found in Halsall et al.14
The AQUAWEB bioaccumulation food web model222 was used to explore how future climate scenarios of increased ambient temperatures (+2 °C and +4 °C), increased primary production (+50% and +100%), and decreased organism lipid content (−10%), would affect contaminant levels through the Arctic marine food web from copepods to seabirds.223 Due to higher predicted levels of primary production, and therefore an increase in organic particles in the water, the net overall bioaccumulation in the food web was reduced. This reduction was greatest for the most hydrophobic compounds subjected to the largest temperature increase used in the model. The reduced bioavailability of POPs predicted by the model resulted from the projected increases in primary production. However, the actual effect of climate change on primary production is uncertain, and the potential for both increased and decreased levels of primary production under future climate change scenarios have been reported.
To approximate a more realistic climate scenario, one study modeled how seasonal changes in abiotic exposure, water temperature, lipid content, and food web structure affected contaminant biomagnification within an Arctic marine food web (Fig. 10).215 The annual cycle of temperature and baseline water contaminant exposure was estimated using the Danish Eulierian Hemispheric Model (DEHM) in climate mode to compare two years a century apart (2007 and 2107).215 The AQUAWEB bioaccumulation food web model,222 parameterized for an Arctic marine food web,223 was used to model seasonally-dynamic, temperature-sensitive parameters affecting contaminant uptake and elimination, seasonal lipid content, and food web structure. The best fit between present day (2007) empirical and modeled seasonal food web biomagnification, as assessed by trophic magnification factors (TMFs), resulted primarily from the inclusion of seasonal variations in lipid content; seasonal changes in temperature, abiotic exposure, and dietary relations had only minor effects on food web biomagnification. Since food web structure did not contribute to explaining variation in seasonal food web biomagnification, it was not altered in the model scenario of future Arctic climate. Modeled food web biomagnification for the year 2107 differed only marginally from that of 2007, with magnitude and direction depending on the contaminant modeled; the mean change from 2007 to 2107 was negative for PCB-52 (−3.1%) and PCB-153 (−1.7%), and positive for γ-HCH (0.4%). Thus, the seasonal variation in food web biomagnification observed over a single year is greater than the change expected in response to a predicted future climate state, at least for the three POPs modeled. These results also emphasize the importance of understanding the different dimensions of environmental and physiological variation for interpreting the future impacts of climate change. Additionally, the model's best-fit seasonal scenario was that which included annual variations in lipid content, indicating this physiological factor is also crucial for making sound predictions of contaminant food web bioaccumulation for the future. Changes in abiotic drivers alone are not sufficient to explain food web bioaccumulation on a temporal scale. As the focus of this study was the effect of climate change on food web accumulation, the DEHM model scenarios for each year were run using the same contaminant emission levels, however, it is expected that changes in contaminant emissions will occur in the future and will have a greater influence than climate change on exposure.14,22
Models have also been used to identify the effect of climate change on food web accumulation by comparing the outcomes of different ecosystem conditions under IPCC-based future climate scenarios. In the Laurentian Great Lakes of Lake Erie and Lake Superior, four climate scenarios were contrasted by combining model parameters on bioenergetics (i.e., temperature and body-size effects on consumption and respiration) and bioaccumulation of PCB-153 for various fish species with different temperature regimes.222,224,225 The authors predicted that future climate conditions would lead to higher growth and feeding rates, and thus increased bioaccumulation. However, the specific conditions leading to higher growth differed by species; high temperatures hampered growth in cold-water species such as the native forage fish, mottled sculpin (Cottus bairdii) and the predator species, lake trout (Salvelinus namaycush), whereas low temperatures lowered the growth of warm-water species, such as the invasive forage fish, the round goby.224 If we speculate on how this translates to Arctic ecosystems, expected temperature increases would result in reduced growth and lower bioaccumulation in Arctic-resident, cold-water species, and increased growth and higher bioaccumulation in warm-water species moving north due to borealization. Note that changes in diet due to climate change-induced alterations in predator-prey interactions were not considered, and thus the effects of feeding and dietary uptake on growth rate and accumulation due to climate change are not known.
More commonly used indicators of risk, such as biomarkers of exposure or effects, are usually tied to a specific toxicant or specific modes and mechanisms of action. There has been a shift, however, towards realizing the need for a broader perspective to address the combined effect of multiple stressors, where integrated responses at a higher level of biological organization are addressed.227 One such response is the overall allocation and use of energy by an organism, which can be assessed using dynamic energy budgets.228–230 Typically, such integrated responses can be linked to the fitness of organisms and are easier to extrapolate to population-level effects using, for example, individual-based population models, as shown for eider ducks in Section D.204 Ideally, climate change and exposure models (see Section E) should be combined with effect models that assess the combined effects of multiple stressors on population levels. These may include models investigating changes in vital rates in response to pollutants and environmental parameters, as suggested for polar bears.208 Still, empirical studies of climate change-contaminant effects on population relevant parameters, such as survival and reproduction, are needed to support and ground-truth model projections.
With the large and growing number of contaminants present in the Arctic, there is still a need to investigate their interrelationships and the underlying structure in their variance. Although not new, statistical approaches (such as principal components analysis) are still very useful tools for condensing, and understanding the multidimensional space represented by contaminants and their drivers in the ecosystem.
Various omics methods are increasingly being used to identify changes in physiological and cellular processes due to contaminant exposure and other factors affected by climate change (e.g., sea ice extent, body condition). Combined with adverse outcome pathways that identify the sequence of cellular events leading to an adverse organism-level effect (e.g., ref. 231), these molecular-based methods hold promise for a more holistic understanding of individual and ecosystem health under multiple stressors by reflecting approaches used in human toxicology and having an integrated ‘one health’ view as goals on the horizon.
In terms of the influence of large-scale climate patterns on POPs in Arctic biota, currently, only a small fraction of the vast collections of climate data and POPs monitoring data that are available have been studied collectively. Of those that have, findings show that the exposure and accumulation of POPs is higher in Arctic biota following time periods with enhanced transport of air and water masses from North America and Europe toward the Arctic, but other underlying mechanisms may also be involved. Additional studies that include more species, time points, and locations (including freshwater and terrestrial ecosystems), are required for examining the relationship of meteorological parameters (i.e., wind and precipitation patterns) to biota POP concentrations; however, the use of climate indices to investigate the effects of climate change on pollutant exposure in Arctic ecosystems is not straightforward.58,232,233 The use of meta-analysis in ecotoxicology is still in its infancy, but could be applied in these studies if both monitoring programmes and individual research studies make raw data publicly available, or at least report metadata beyond its specific use in a particular study.234,235 For future studies, high-resolution spatial and temporal data and variability (e.g., local sea ice conditions, daily and seasonal temperature variations) should be considered, instead of climate variables averaged over large spatial and temporal scales, if possible. Considering the potential effects of intermittent and/or extreme weather events (e.g., rain-on-snow events) is also important as these anomalies are predicted to increase in response to a warming climate.
The melting of multi-year sea ice, just like the melting of permafrost and glaciers, may release stored contaminants into the Arctic marine food web. Although this could cause short-term increases in POP levels in biota, the eventual complete replacement of multi-year ice with annual ice, which may be less contaminated due to more frequent freeze–thaw cycles, should lead to reduced POP exposures via this pathway. Research findings available to date show generally, but not consistently, associations between sea ice condition and contaminant levels in biota, likely resulting from the influences of other species-, food web-, and ecosystem-specific factors. Currently, it is difficult to separate the direct effects of declining sea ice cover on POP levels in biota from the indirect effects resulting from changes in food web structure and function. As making robust conclusions by comparing single climate- and ecological-parameters (e.g., sea ice extent, primary productivity, dietary composition) to POP concentrations in biota can be challenging or impossible, on-going studies should strive to measure both abiotic and biotic parameters to simultaneously compare the strength of their associations with POP levels in Arctic biota. Studies should include ancillary biological and ecological data (e.g., body condition, lipid content, fatty acid profiles, stable carbon, nitrogen, and sulphur isotopes, and compound-specific isotopes (e.g., amino acid nitrogen isotopes)) to disentangle the impacts of climate-driven physical changes from those resulting from ecological or physiological effects on exposures in biota. Use of models and other research tools that account for multiple environmental and biotic variables could be useful in determining the contribution of single variables to changes in contaminant levels observed in wildlife.
Increased terrestrial runoff is predicted to affect both exposure and food web biomagnification of POPs, but there are few data to test or support this prediction. Terrestrial runoff affects nutrients, organic matter quality and quantity, and contaminants, in a complex process that influences both food web ecology and the availability and uptake of POPs and CEACs in biota. Additional experimental evidence for effects of terrestrial runoff on biota POP levels in aquatic and coastal systems is required.
Recent studies suggest that permafrost degradation and thaw can increase the amount of POPs bound to sedimentary organic matter and increase levels in amphipods in small lakes. There are knowledge gaps requiring further research to better understand the impacts of changing water quality (i.e., increasing nutrient concentrations and turbidity) on contaminant dynamics in freshwater food webs. These gaps include impacts on fish physiological conditions in lakes and rivers impacted by permafrost RTS, links between water quality and POP trends in fish from Arctic freshwater environments, and changes in water quality due to permafrost degradation and RTS. Additional studies are also needed from a wider range of lakes across the circumpolar Arctic.
Changes in seasonal timing and species phenology, e.g., migration, reproduction, and food availability, may have potential impacts on levels of POPs in Arctic biota. Many of these changes are dependent on species- and ecosystem-specific characteristics and may result in specific windows of elevated vulnerability during the year and during the organisms' life span. As the effects of phenological changes on contaminant bioaccumulation are not known, future empirical studies should assess whether further investigation of these factors is warranted. The cascading effects of changes in seasonal timing and potential mismatches in Arctic communities are unknown, warranting interdisciplinary collaborations with biologists/ecologists in future contaminant studies. At the local- and regional-levels, indigenous knowledge from hunters and fishers may provide useful data in this regard.
The net effect of increasing primary production on contaminant dynamics in Arctic food webs is unknown as predictions of future primary production and downstream impacts on dissolved versus particulate POP concentrations are uncertain. Yet, increased primary production in combination with temporal mismatches with grazers could lead to increased sedimentation and contaminant delivery to the benthic ecosystem, although changes in benthic-pelagic coupling due to the northward-range shifts of pelagic fishes may influence this. Given the potentially counteracting effects of increased primary production on contaminant dynamics (i.e., lower food web exposures via biodilution versus higher food web exposures via increased dietary intake of particulate-bound POPs), additional research on these processes should be undertaken.
With many unknowns and variable findings from different species and ecosystems, the net effect of climate change-driven alterations in species interactions on contaminant accumulation in food webs is not known. Sea ice decline may reduce the dietary importance of marine foods leading to lower intakes of POPs for some species (e.g., polar bears). Conversely, a foraging switch from Arctic- to boreal-prey species may lead to higher POP concentrations in predators if their new prey are seasonal migrants that act as contaminant biovectors into Arctic marine food webs, although the importance of this phenomenon is not fully understood. Dietary tracers and other ancillary biological and ecological data should be included in national biomonitoring programs to track changes and interpret ongoing contaminant datasets. Such data may assist in evaluating the role of changing species interactions on contaminant levels if/when trends deviate from expected patterns. Adjusting for some of these parameters will aid the interpretation of trend data. Modeling future scenarios of shifting food webs and their impacts on POP exposures in top predators should be pursued in subsequent studies once a better understanding of how climate change affects food web structure and function is developed.
Lipid amount, quality, and dynamics within organisms and ecosystems are expected to change with climate change, and thus have potential consequences for contaminants, especially lipophilic POPs. Decreased body condition can lead to increased blood levels of POPs in biota. When feasible, lipid content and body condition metrics should be routinely measured in contaminant trend studies and incorporated into data interpretations as changing lipid dynamics could introduce variability into POP temporal trends and cross-comparisons between individuals or populations. Investigation of the behavior of proteinophilic contaminants (e.g., PFASs and Hg) with changes in lean body mass may also be warranted.
Few studies have examined the consequences of behavioral changes resulting from climate change on exposure and accumulation of POPs in biota. Additional research examining the impact of behavioral changes on POP exposure is warranted. This should be done by studying relationships between pollutant exposure, movement ecology, diet, food web structure and energetics, and would benefit from multidisciplinary research teams and approaches to data collection. As behavioral and dietary changes may affect contaminant levels in biota, they may also introduce variation into long-term temporal trends that should be identified and accounted for.
Few studies have examined how climate change may interact with other stressors to affect POP toxicity in biota in the Arctic or elsewhere. An improved understanding of the effects of multiple stressors, inclusive of climate change impacts, on contaminant toxicity in Arctic biota is crucial. Knowledge of the effects of contaminants on the abundance, development, health, and behavior of wildlife populations, and the mechanisms and processes underlying these effects, is needed for toxicants individually and in combination with other stressors (i.e., climate change, pathogen and parasite exposure, etc.), as effect thresholds under these more realistic scenarios may be different from those based on exposure to a single contaminant or stressor in isolation. Protecting Arctic biota and people from the effects of POPs requires that contaminant researchers expand their scope of research and collaboration to consider POP exposures and effects within the true context of many co-occurring and potentially interacting anthropogenic and environmental stressors facing populations, communities, and ecosystems in the Arctic.
Models are important tools for evaluating the net result of contrasting processes, as climate change can lead to both increased and decreased exposures and enhanced or reduced bioaccumulation processes, depending on the processes, areas, chemicals, and species interactions affected. Mechanistic process-oriented models combining physicochemical transport models with bioenergetic and food web models, can help to understand how the variation in different parameters, such emission sources, transport processes, temperature, lipid dynamics, and food web structure, affect exposure, uptake, elimination, and overall bioaccumulation of contaminants. Nonetheless, for some of these parameters (e.g., food web changes), there remains insufficient understanding or empirical data to populate models. Models combining climate change impacts on abiotic and biotic processes affecting POPs and CEACs should be generated, as possible, and evaluated with measured data. In some cases, and related to other recommendations, additional research is required to constrain model parameters. In addition to including processes and physicochemical parameters tied to chemical emissions, transport and partitioning, more ecological realism, inclusive of variation, is needed for parameters related to ecosystem structure and functioning that may affect bioaccumulation and magnification.
We recommend that future research identifies and characterizes the direct linkages between climate-induced changes in LRT of contaminants to the Arctic, and the exposure, accumulation, and toxicity of POPs and metals to wildlife. It will be challenging to separately identify the effects of climate change and LRT on Arctic wildlife that are also forced to cope with many other stressors that also influence the accumulation and toxicity of these chemicals and metals. However, forward progress may be made with comprehensive monitoring of Arctic ecosystems and wildlife in combination with lab-based mechanistic studies, new research methods and models that reflect the complexity of natural variation and species interactions, and risk assessments that holistically consider the effects of multiple exposures and stressors on individual and population health. As changes in the structure and function of Arctic ecosystems are not currently predictable, to project the effects of climate-driven changes in food webs on POPs in Arctic wildlife, we recommend that the contaminants community continue to work closely with other research disciplines.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1em00469g |
‡ These authors contributed equally to the work. |
This journal is © The Royal Society of Chemistry 2022 |