Modelling PCB-153 in northern ecosystems across time, space, and species using the nested exposure model †

There is concern over possible e ﬀ ects on ecosystems and humans from exposure to persistent organic pollutants (POPs) and chemicals with similar properties. The main objective of this study was to develop, evaluate, and apply the Nested Exposure Model (NEM) designed to simulate the link between global emissions and resulting ecosystem exposure while accounting for variation in time and space. NEM, using environmental and biological data, global emissions, and physicochemical properties as input, was used to estimate PCB-153 concentrations in seawater and biota of the Norwegian marine environment from 1930 to 2020. These concentrations were compared to measured concentrations in (i) seawater, (ii) an Arctic marine food web comprising zooplankton, ﬁ sh and marine mammals, and (iii) Atlantic herring ( Clupea harengus ) and Atlantic cod ( Gadus morhua ) from large baseline studies and monitoring programs. NEM reproduced PCB-153 concentrations in seawater, the Arctic food web, and Norwegian ﬁ sh within a factor of 0.1 – 31, 0.14 – 3.1, and 0.09 – 21, respectively. The model also successfully reproduced measured trophic magni ﬁ cation factors for PCB-153 at Svalbard as well as geographical variations in PCB-153 burden in Atlantic cod between the Skagerrak, North Sea, Norwegian Sea, and Barents Sea, but estimated a steeper decline in PCB-153 concentration in herring and cod during the last decades than observed. Using the evaluated model with various emission scenarios showed the important contribution of European and global primary emissions for the PCB-153 load in ﬁ sh from Norwegian marine o ﬀ shore areas.


Introduction
Legacy persistent organic pollutants (POPs) are regulated on a global scale due to their persistence, toxicity, and potential for bioaccumulation and long-range transport. 1As a result of these properties, POPs are ubiquitous in the environment, also in Arctic ecosystems far away from emission sources. 2,3For instance, polychlorinated biphenyls (PCBs) reaching high concentrations in Arctic top predators such as seabirds and polar bears, [4][5][6] are still of ecotoxicological concern decades aer they were prohibited. 7,8][11][12][13][14][15][16][17] Marine ecosystems are a key source of nutrients and contaminants for both marine top predators [18][19][20] and humans, 21,22 with sh playing a key role as vectors of contaminant exposure.Hence, understanding contaminant sources, pathways and exposure of marine sh is also vital for understanding and managing contaminant exposure of humans and top predators.The Norwegian marine areas are important commercial shing grounds, 23 as well as feeding areas for numerous seabird, seal, and whale species.These areas range from the North Sea and Skagerrak impacted by local sources in the south, to the remote Barents Sea in the north where long-range environmental transport is the main source.7][28] In addition, eld-based research contributes important insights into POPs exposure across environmental media 29,30 and trophic levels. 31,32owever, environmental sampling and chemical analysis for POPs are expensive and time-consuming, limiting the availability of eld-based measurements across time, space, and species.Moreover, measured concentrations alone are oen not enough to provide information on contaminant sources, particularly when exposure may arise from global emissions.As a result, a comprehensive understanding of the full link between global emissions and biotic exposure is difficult to achieve based on measurements alone.
Multimedia physical fate and bioaccumulation models have been developed to further increase our understanding of the environmental fate and bioaccumulation behavior of organic contaminants. 33Objectives of earlier modeling studies on contaminants in sh and biota in Norwegian marine areas included increasing the understanding of bioaccumulation in Arctic marine species and food webs 32,34,35 and how this may be impacted by climate change. 36Models have also been used to identify chemicals of concern, 37 and to elucidate the role of plastic in organic contaminant bioaccumulation 38 in this region.Most of these studies relied on steady-state models not accounting for temporal or spatial variability, even though high latitude ecosystems are highly dynamic due to a strong seasonality in light regimes, temperatures, primary production, etc.An exception is the OMEGA model, which has been used to investigate the seasonality of PCB bioaccumulation in the Barents Sea 34 as well as time-trends of a range of contaminants. 37However, this model is not spatially resolved.Environmental fate and bioaccumulation models can greatly benet from taking into account that emissions and environmental characteristics vary in space and time. 39In fact, increased capability for estimating temporally and spatially variable concentrations has been identied as a key research need to improve exposure estimation in ecological risk assessment in the 21st century. 40,41urthermore, none of the earlier studies modeling bioaccumulation in Norwegian marine areas included a simulation of chemical fate and transport in the physical environment.They instead relied either on aqueous concentrations as input, taken from measurements in water 32,35,38 or derived from measured concentrations in biota, 32,35,37 or on hypothetical scenarios. 34,36This approach does not facilitate an understanding of the whole link between chemical emissions to the environment and biotic exposure, which is crucial for enabling scientically-sound management strategies.Combined models of environmental fate and bioaccumulation have proven fruitful to increase our understanding of contaminant behavior and exposure.An example is the CoZMoMAN model 42 which has been used both to evaluate the environmental fate of shortchain chlorinated paraffins in the Nordic environment 43 and to increase our understanding of human exposure to PCBs. 44,45he CoZMoMAN model is limited by its lack of spatial resolution and because it cannot realistically model contaminant transport into the model domain from the outside world.The bioaccumulation model ACC-Human 46 has been combined with global dynamic and spatially resolved environmental fate models like Globo-POP 2,47,48 and BETR-Global 49 to quantitatively describe the whole link between global contaminant emissions and exposure to humans and wildlife.For example, BETR-Global and ACC-Human have been used to estimate human exposure to PCB-153 on a global scale. 50,51Also, Globo-POP and an Arctic version of ACC-Human have been used to distinguish the physico-chemical properties of chemicals likely to be Arctic contaminants. 3][55] Recently, the Nested Exposure Model (NEM) 56 was introduced for the physical environment as a global dynamic multimedia environmental fate and transport model that can zoom into a user-dened region of interest with increasing spatial resolution.The nested feature of NEM reduces computational demands, while allowing for simulations that are global in scale yet with a high spatial resolution in the region of interest.NEM has been shown to satisfactorily reproduce fate and concentrations of PCB-153 in the European atmosphere, with differences between observations and model estimates comparable to sampling and analytical variability across involved laboratories. 56However, NEM has not yet been evaluated for environmental compartments other than air.Moreover, while Breivik et al. 56 noted that "NEM is also integrated with the ACC-Human bioaccumulation model", this was not further evaluated or explored.Hence, the main objectives of this study are to (i) further develop NEM for northern marine ecosystems, (ii) evaluate NEM's ability to estimate PCB-153 concentrations across time, space, and species, and (iii) apply the model to explore the sources of PCB-153 to northern marine ecosystems.The focus of the evaluation is on an Arctic marine food web in Svalbard, Norway, and on two ecologically and commercially important sh species (Atlantic herring (Clupea harengus) and Atlantic cod (Gadus morhua)).These case-studies were selected because of the availability of consistent measured data across time, space, and species and their relevance to management for protecting northern ecosystems and human health.

Materials and methods
The NEM model structure The dynamic and spatially resolved fugacity-based NEM model (Fig. 1) consists of two parts: (i) a module that simulates environmental fate and behavior of organic contaminants in the physical environment 56 and (ii) a bioaccumulation module based on the ACC-Human, 46 Arctic ACC-Human, 3 and expanded Arctic ACC-Human models. 52,54NEM uses global emission scenarios to estimate concentrations in air, water, soil, and sediment in a user-dened domain (in time and space) with a user-dened spatial resolution ranging from 30°× 30°l atitude/longitude (lat/long) to 0.5°× 0.5°lat/long.The nested feature of NEM means that it is possible to model the entire globe with a coarse resolution, then zoom into an area of interest with increased spatial resolution while accounting for chemical inow from the outside world into the nested domain.In doing this, the model uses the output from a larger domain with coarser resolution as boundary conditions for a smaller domain with ner resolution. 56Next, estimated concentrations or fugacities in the physical environment are used to simulate bioaccumulation behavior of organic contaminants in a pelagic food web typical for Norwegian marine ecosystems (Fig. 1).The bioaccumulation module can be run for a user-dened spatial domain that is embedded in the domain of the physical module with the same spatial resolution.Currently, organisms cannot move between grid cells, but are assumed to remain within a single grid cell throughout their life.This simplication can lead to errors when assessing contaminant exposure of migratory species.The temporally resolved concentrations or fugacities in the physical environment driving the bioaccumulation module in NEM need not be output of the physical fate module, but can be user-generated les, e.g. based on environmental measurements.

Species included and biota parameterization
Nine species of key importance to Norwegian sub-Arctic and Arctic pelagic marine ecosystems make up the food web in the current version of NEM (Fig. 1): 57 three pelagic zooplankton groups with different feeding strategies 58,59 (copepods, krill, and amphipods, ESI S1.1 †), four sh species (Atlantic herring, 46 capelin (Mallotus villosus), 52 polar cod (Boreogadus saida), 3 and Atlantic cod 46 ), and two Arctic marine mammals (ringed seal (Pusa hispida) 3,52,60 and white whale (Delphinapterus leucas) 52 ).Existing modules for lichen, caribou, narwhal, and humans are also integrated in NEM but are not included in this study. 46,52,54ontaminant concentrations in all three zooplankton groups are assumed to be in equilibrium with seawater based on bioaccumulation factors measured in King's Bay, Svalbard. 59Agedependent lipid-normalized concentrations of contaminants in sh and mammals are calculated from mass balance equations comprising key uptake and elimination mechanisms, as described previously. 3,46,52,54,60Details of the parameterization of the organisms, especially updates to earlier versions of the Arctic ACC-HUMAN models, 3,52,54,60 can be found in ESI S1. † In particular, seasonal lipid-dynamics have been included for all plankton and sh species, except for Atlantic cod which has a relatively constant lipid content (S1.1 and S1.2 †).Parameterization of the two mammals has not been changed. 3,52,60her input parameters All environmental input parameters to the physical fate module of NEM have been described in detail previously. 56Most spatially resolved environmental parameters have a resolution of 0.5°× 0.5°, and the model automatically calculates the spatially resolved environmental parameters at the selected resolution based on these parameters at 0.5°× 0.5°. 56The bioaccumulation module uses the same spatially and temporally resolved seawater and air temperatures as direct input as used in the physical fate module of NEM. 56or evaluation and application of the model, we selected PCB-153 (2,2 ′ ,4,4 ′ ,5,5 ′ -hexachlorobiphenyl) because (i) NEM has already been shown to reproduce PCB-153 concentrations in the European atmosphere within a factor of 1.5 of observations, 56 (ii) it is still of toxicological importance in wildlife, 7 (iii) empirical data for PCB-153 in marine and Arctic biota is plentiful, and (iv) a global historical emission scenario 61,62 is available.The most up-to-date default global emission scenario for PCB-153 from 1930 to 2100 (ref.63) was used as input, assuming emissions solely to air, as in Breivik et al. 56 Physicalchemical property data used as input for PCB-153, including equilibrium partition ratios between air and water (K AW ), octanol and water (K OW ), and octanol and air (K OA ) 64 and their temperature dependence were also harmonized with the NEM module for the physical environment and are described in Breivik et al. 56 Biotransformation of PCB-153 was assumed negligible in zooplankton and sh.Biotransformation rates and feces-blood partition coefficients (K FB ) of PCB-153 in the mammals were used as in Binnington et al. 52 Simulations for model evaluation Simulations for PCB-153 in the physical environment were carried out for the years 1930 to 2020, starting with global simulations at 5°× 5°resolution, followed by additional simulations for the Norwegian marine areas and bordering regions (55°N -80°N, 10°W -45°E) at 1°× 1°resolution (Fig. 1).Based on the calculated fugacities of PCB-153 in the atmospheric gas phase and in the dissolved seawater phase, PCB-153 fugacities and concentrations were estimated for zooplankton, cod, and herring with spatial resolutions of both 5°× 5°and 1°× 1°(Fig.S2 †).Simulations for all species, including marine mammals, were run for the Svalbard region only (75°N -80°N, 10°E − 20°E) (Fig. S2 †), also at both resolutions.We ran simulations for all species for the Svalbard region only because measurement data for model evaluation are available for this region and to save computing time.All simulations were run with a 12 h time-step, with results stored aer every 1752 h (73 days), and results based on both spatial resolutions are presented.

Measurement data for model evaluation
6][67] When a sample had been collected along ship transects spanning up to several degrees latitude/ longitude, the median and range of concentrations estimated at the time of sampling in the grid cells covering the transect were used for comparison.
Next, we used the two case-studies ("Svalbard" and "Norwegian sh") for more detailed model evaluation.For evaluation of model performance across compartments and species in the Svalbard area (case-study "Svalbard"), we compared estimated concentrations with PCB-153 concentrations measured in bulk air at the Zeppelin observatory (2007-2016), 68 in dissolved seawater phase (July 2008) 59 and whole-body zooplankton, herring, polar cod and Atlantic cod (May-Oct 2007) 31 from King's Bay, Svalbard (Fig. S4 †), as well as in blubber of ringed seal (May-Sept 2014) and white whale (July-August 2014-2016) from the larger Svalbard area (Fig. S5 †). 69King's Bay falls entirely within one grid cell based on NEM's spatial resolution of 5°× 5°, but is located at the intersection of four cells in a 1°× 1°grid (Fig. S4 †).Ringed seal and white whale were sampled across Svalbard, spanning two 5°× 5°grid cells (Fig. S5 †).For estimates based on 1°× 1°, the median and range of concentrations estimated for the three grid cells covering most of King's Bay (Fig. S4 †) were used.For estimates based on 5°× 5°, median and range of estimated concentrations in the one (air, seawater, plankton, and sh) or two cells (marine mammals) covering King's Bay and the wider Svalbard region, respectively, were used.In all cases, concentrations estimated for the same year and season as when the samples were collected were used for comparison.For sh and white whale, estimated concentrations for all age groups were used for comparison.For ringed seal whose age had been determined, 69 only estimates for seal aged 5-20 years (females) and 0-30 years (males) were used to increase comparability to the measured data.For marine mammals, the median estimated concentration in total blubber was used. 69or evaluation of model performance across space and time, we relied on an extensive dataset of PCB-153 concentrations measured in muscle of Norwegian spring-spawning (NSS) herring and liver of Atlantic cod from large baseline and monitoring studies in Norwegian marine areas performed by the Institute of Marine Research (case-study "Norwegian sh") (Fig. S6 and Table S1 †). 26,28,70This dataset includes measured data from the Skagerrak, North Sea, Norwegian Sea and Barents Sea, as well as Norwegian coastal and ord areas (Fig. S6-S8 †) covering the years 1995-2017 for herring (n = 1254) and 2002-2018 for cod (n = 3651).Some of the cod have been sampled in areas impacted by local sources, mostly from historically PCBcontaminated sediment acting as a secondary source to the local marine environment.These samples were agged as "may be impacted", "slightly impacted", or "heavily impacted" (Table S1 †), based on the relative vicinity of the sampling location to local sources.Heavily impacted samples (e.g., from the inner Osloord and Bergen harbor) (n = 229) were removed from the dataset, as the spatial resolution of NEM is too coarse to resolve concentration gradients at the local scale.Estimated lipidnormalized concentrations (ng per g lipid weight (lw)) of PCB-153 in individual sh were used for comparison, if the measured data had accompanying information about (i) lipid content and (ii) age of the sh, (iii) the time of sampling (date and year), and (iv) the sampling position (in lat/long) (Table S1 †).This was the case for more than 3000 sh (n = 783 for herring, n = 2291 for cod).For these, the measured lipidnormalized PCB-153 concentration in individual sh was compared to the estimated lipid-normalized PCB-153 concentration for a cod or herring of the same age at the time and location closest to the time and place of sampling.For sampled sh older than 10 years, the estimated concentration for a 10 year-old sh was used (n = 38 for herring and n = 4 for cod).Because such auxiliary data are missing for sh sampled prior to 2006, earlier data (n = 310 for herring, n = 29 for cod) were only included to evaluate the overall time-trends on a ng per g wet weight (ww) basis.The NEM model was evaluated for its ability to reproduce (i) PCB-153 concentrations in individual sh, (ii) spatial trends of PCB-153 in sh, and (iii) temporal trends of PCB-153 in sh.

Exploration of emission scenarios
Three country-or region-specic emission scenarios for PCB-153 for 1930-2020 were explored to evaluate the impact on concentrations in Norwegian sh: (1) primary emissions within

View Article Online
Norway only, (2) primary emissions within the countries that were European Union (EU) members as of 1973 only (Germany, France, Italy, Netherlands, Belgium, Luxembourg, Denmark, Ireland, United Kingdom), and (3) primary emissions within Russia only.The results from these scenarios were compared to the default scenario based on global emissions to investigate the relative contribution of historical primary emissions in Norway, EU, and Russia to PCB-153 contamination in Norwegian sh.For this purpose, estimated lipid-normalized concentrations in a 5 year-old herring were compared between scenarios.For example, the contribution from Norwegian primary emissions to the total body burden of PCB-153 concentrations in a 5 year old herring was calculated as the estimated concentrations in the herring based on primary emissions within Norway only (scenario 1) divided by the estimated concentrations in herring based on global emissions (default scenario).As the model responds linearly to emission changes, results will be identical for cod and herring in the same grid cell.

Evaluation of model performance
To evaluate model performance, several quantitative metrics were used.The ratio between estimated and measured concentrations (prediction-measurement ratio, PMR) can be calculated for individual pairs of estimated and measured concentrations and provides information on both the extent of agreement between estimates and measurements (the closer to 1, the better agreement) and direction of agreement (underprediction (PMR < 1) or overprediction (PMR > 1)).
It has been argued that an agreement between model estimates and eld measurements of an order of magnitude can be used as a rule of thumb when evaluating the performance of multimedia models. 71This compares to a PMR between 0.1 and 10.
To evaluate overall model performance, the model bias (MB) 72 was calculated as Hence, MB is the geometric mean of the PMRs for e.g.all species (Svalbard) or all individual cod or herring (Norwegian sh), and provides information on the overall performance of the model, including the extent (the closer to 1, the better agreement) and direction of agreement (underprediction (MB < 1) or overprediction (MB > 1)). 72imilarly, the root mean square error of the log-transformed concentrations 39 was calculated as and summarizes both random error and systematic bias. 39The smaller RMSE log is, the better model performance.
In addition, Spearman's coefficient (r s ) was calculated for correlations between estimated and measured PCB-153 concentrations in individual organisms.The higher r s is, the better the correlation.

Results and discussion
Model evaluation for seawater NEM has already been shown to successfully estimate PCB-153 concentrations in air. 56Prior to applying NEM to marine biota, we evaluated its ability to estimate PCB-153 concentrations in seawater.Measured concentrations included in the evaluation ranged from 0.02 pg L −1 (Barents Sea, 1999) 32 to 12.6 pg L −1 (Skagerrak, 1989), 66 spanning more than two orders of magnitude.This range represents both a temporal decline due to the global phase-out of PCBs as well as spatial variation due to variable distance from source areas (see map in Fig. S3 †).Measuring concentrations of PCB-153 dissolved in seawater is difficult due to its low water solubility, and the measurements have relatively high inherent uncertainty. 32Estimated concentrations agreed better with measurements when based on a spatial resolution of 1°× 1°(MB = 1.9, RMSE log = 5.9) than 5°× 5°(MB = 3.2, RMSE log = 6.9) (Fig. S9 †).This mirrors the results of the model evaluation for the atmosphere and illustrates the effect of numerical diffusion. 56Based on 1°× 1°, the PMR in seawater ranged from 0.1 to 31, with 64% of estimates within a factor of 4 of measurements (n = 33) (Fig. S9 †).

Model evaluation across species (Svalbard)
Estimated concentrations (1°× 1°resolution) of PCB-153 in air and seawater for King's Bay displayed PMRs of 0.89 and 1.1, respectively (Table S2 †).Estimated concentrations in plankton, sh and marine mammals from King's Bay and Svalbard are compared to measurements in Fig. 2 (1°× 1°resolution) and S10 (5°× 5°resolution).† As discussed previously for the atmosphere 56 and seawater (above), model estimates for Svalbard biota agreed best with measurements based on the spatial resolution of 1°× 1°(Table S2 and Fig. 2, S10 †), demonstrating again improved model performance at increased spatial resolution.This is partly due to the good agreement between measured and estimated seawater concentrations in King's Bay (PMR = 1.1), because the estimated concentrations in waterrespiring biota scale linearly with the seawater concentrations (i.e., if water concentrations are doubled, so do concentrations in biota).Seen in the context of the overall evaluation of seawater concentrations where PMRs range from 0.1 to 31 (1°× 1°, Fig. S9 †), this may be partly fortuitous.
Model estimates (1°× 1°) across species displayed PMRs between 0.14 and 3.1, an overall MB of 1.3 and RMSE log of 3.7 (Table S2 †).Estimated concentrations in female white whale displayed the greatest deviation from measurements and were underpredicted compared to measurements (PMR = 0.14).This was also observed for PCB-153 estimates in female white whales in the Canadian Arctic using the Expanded Arctic ACC-Human model by Binnington et al., 52 who suggested that the model overestimates reproductive loss of PCBs through calving and/or lactation.If female white whale is excluded, the PMR varied between 0.89-3.1 across plankton, sh, and marine mammals.This is deemed acceptable considering all inherent uncertainties in input properties, model assumptions, as well as measurement data.Most importantly, this means that NEM succeeds in simulating the whole continuum of processes linking global historical emissions of PCB-153 with marine ecosystem exposure in the Norwegian Arctic, including bioaccumulation across trophic levels.
To investigate NEMs ability to reproduce bioaccumulation across trophic levels in more detail, we also compared measured and estimated trophic magnication factors (TMF) using the data in Fig. 2. For zooplankton groups and sh species, average trophic positions (TP) from Hallanger et al. 31 were used.For mammals, trophic positions were calculated from stable nitrogen ratios (d 15 N) in skin, 73 measured in the same individuals of ringed seal and white whale as the PCB-153 concentrations. 69The TP of the marine mammals and the TMFs were calculated following the method by Hallanger et al. 31 for consistency.Separate TMFs were calculated for poikilotherms only (zooplankton and sh) and for all organisms (including ringed seal and white whale), as including homeotherms can signicantly impact calculated TMFs.Linear regressions between log-normalized PCB-153 concentrations and TP showed signicant positive relationships for both measurements and estimates (Fig. S11 †).The TMFs for poikilotherms were 3.89 and 4.21 for measurements and estimates, respectively.This is within, or close to, the range of the seasonal TMFs for PCB-153 (1.53-4.07)calculated by Hallanger et al. 31 based on the same measurement data, but including some more species of sh (haddock and pollock) and zooplankton (arrow worms (Chaetognatha)).When we included the homeotherms in the regressions, the TMFs increased to 8.82 and 6.08 for measurements and estimates, respectively.Increased TMFs when homeotherms are included have also been seen previously. 31,74here was a larger discrepancy between measured and estimated TMF when homeotherms were included.This is partly due to the underprediction of PCB-153 concentrations in female white whale.If female white whale is excluded from the regression of the estimated concentrations, the estimated TMF increase to 7.51.
Herring and cod display migratory behavior to varying degrees and have not spent their entire lives at the sampling position.Hence, the measured concentrations in sh reect exposure from a wider geographical area than their sampling position.Northeast Atlantic cod migrates seasonally between the Barents Sea and the Norwegian Sea to spawn, 75 while coastal cod and North Sea cod are more stationary. 76The NSS herring undertake seasonal migrations between coastal and offshore areas in the Norwegian Sea. 77This is not yet accounted for in NEM, where sh are assumed to reside in individual grid cells for their entire lifecycle.Hence, using a coarser spatial resolution (5°× 5°) may be more representative for their real-life exposure than the ner spatial resolution (1°× 1°).This agrees with the results obtained herein, where all metrics of model performance (PMR, MB, RMSE log and r s ) for estimates of PCB-153 in individual sh were better when NEM was run with the coarse resolution (5°× 5°), except for the average PMR and model bias for herring (Table S4 †).The improved MB and average PMR for herring at 1°× 1°are because estimates are more spread out with concentrations both lower and higher than measurements (Fig. S13 †), whereas there is a clearer tendency for overprediction but less spread at 5°× 5°resolution (Fig. 4).To explore the effect of migration, two simplied migration scenarios for herring and cod, respectively, were constructed to explore how spatial migration impacts estimated concentrations of PCB-153 as well as model performance (5°× 5°r esolution, ESI Section S8 †).Overall, accounting for migration did not improve estimates for herring and cod (ESI Section S8 †).Hence, the remainder of the results and discussion is based on 5°× 5°spatial resolution without accounting for migration between grid cells.
Model performance (5°× 5°) was better for individual herring than for individual cod, with a lower RMSE log (4.71 compared to 5.66), higher spearman r s (0.65 compared to 0.22), and comparable MB (2.32 compared to 0.45) (Fig. 4 and Table S4 †).There was a tendency of overpredicting concentrations in herring (MB = 2.32, i.e. 2.32-fold overprediction), and underpredicting concentrations in cod (MB = 0.45, i.e. 2.22-fold underprediction), with 77% and 70% of estimates for herring and cod, respectively, within a factor 4 of measurements.NSS herring is a schooling, pelagic sh stock migrating throughout the Norwegian Sea, and is relatively homogenous in terms of diet, habitat, and geographical distribution. 77Consequently, the dataset used for model evaluation of PCB-153 in herring is also relatively homogenous.On the other hand, cod has a much more diverse diet, and may feed on both pelagic and benthic species. 78,79A benthic diet of cod is not yet accounted for in NEM.9][80] The model does not account for substantial differences in the biology of these cod stocks, including differences in migration, life history, and diet.Furthermore, some of the cod have been sampled in areas suspected, or known, to be inuenced by local (secondary) sources of pollution (Tables S1 and  S5 †).The larger diversity of the cod dataset may explain the poorer model performance when compared to herring, as NEM View Article Online is currently not able to capture this diversity in diet, ecology, and space-use strategies.
To explore whether there was any spatial variation in the performance of NEM, the logarithm of the PMR was plotted against number, latitude, and longitude of the grid cell for both herring and cod (Fig. S14 †).Even though there was variation in the PMR, there was no systematic bias with either grid cell number, latitude, or longitude.This means that there were no particular cells or spatial regions that could explain the overprediction of herring or the underprediction of cod.This also indicates that NEM performs equally well for the whole Norwegian marine area, ranging from the North Sea in the South to the Barents Sea in the North.To explore this in more detail, the Norwegian marine areas were divided into four regions (Skagerrak, North Sea, Norwegian Sea, Barents Sea) as in Nøstbakken et al., 81 and the cod measurements were allocated to these regions depending on their sampling location.Measured concentrations of PCB-153 in cod liver decreased in the sequence: Skagerrak > North Sea > Norwegian Sea > Barents Sea (Table S6, Fig. 5 and S15 †), in agreement with spatial patterns observed by Ho et al. 28 and Everaert et al. 25 NEM reproduced this trend of decreasing concentrations from south to north (Table S6 † and Fig. 5), with the exception of slightly higher estimated concentrations in the Norwegian Sea than in the North Sea, contrary to observations.Lower median estimated PCB-153 concentrations in North Sea cod compared to measurements (PMR = 0.25) than in the three other marine regions (PMRs 0.53-0.89)(Table S6 and Fig. S15, S16 †) could have several explanations.Broadly, it can be caused by uncertainties in (i) the global emission estimates used as input to the model, (ii) estimated fate and persistence in the physical environment, and/or (iii) estimated bioaccumulation.NEM could reproduce measured temporal trends in air at the Birkenes observatory in southern Norway, close to the North Sea. 56This suggests that the discrepancy by NEM for cod in the North Sea is not likely to be due to the emission estimates.However, seawater and seawater sediments are expected to respond more slowly to emission changes than the atmosphere.Thus, it cannot be excluded that PCB-153 concentrations in the North Sea have responded more slowly to declining primary emissions the last decades than estimated by the model, e.g., due to underestimating secondary sources.Moreover, there could be aspects of the diet or ecology of North Sea cod that is not adequately represented in NEM and underestimates their reallife exposure.For example, North Sea cod grows faster than cod in the Barents Sea. 80As a result, North Sea cod is larger at the same age than cod in the Norwegian Sea, and could possibly be feeding at a higher trophic level.A strong benthic link, particularly in combination with historically contaminated sediment, could also underestimate their real-life exposure as the prey species for cod in NEM is currently only pelagic.
The estimated temporal trends of PCB-153 concentrations in ng per g ww in herring and cod in Norwegian marine areas for 1930-2020 are shown in Fig. 6.Concentrations were not lipidnormalized to allow for the use of data for sh collected before 2006 when lipid content was not recorded.2][63] The large range in estimated concentrations in Fig. 6, particularly for cod, is caused by combining results for several grid cells to look at the overall results for Norwegian marine areas (see Fig. S17-S19 † for selected individual grid cells).The estimated time-trends overlapped with the measured time-trends in both herring and cod (Fig. 6).To investigate the agreement between the estimated and measured time-trends in more detail, we calculated the environmental half-lives (t 1/2 ) of PCB-153 in herring and cod for both measured and estimated concentrations assuming rst-order kinetics using eqn (1) and (2): where C is the PCB-153 concentration in herring or cod t years aer the rst year of the measured time-series (t 0 ), C 0 is the median PCB-153 concentration in herring or cod at t 0 , and k is the rate constant of decline.For cod, we calculated separate half-lives for the North Sea, Norwegian Sea, and the Barents Sea, as the number of sh sampled every year in each region varies, and this can impact the calculated half-life since the PCB-153 concentrations vary between the different sea regions (Fig. 5).Skagerrak was not included in this analysis, as the dataset contained cod sampled in Skagerrak for only two separate years.
The measured concentrations showed a signicant decline with time for PCB-153 only in cod liver from the Barents Sea (Fig. S20 and Table S7 †), with a calculated half-life of 32.6 years.In contrast, NEM estimated a consistent, signicant decline in PCB-153 concentrations for both species and all sea regions with half-lives ranging between 7.5-9.5 years.Measured PCB-153 concentrations showed non-signicant declines for herring and cod in the Norwegian Sea, and a signicant increase with time for North Sea cod (Fig. S20 and Table S7 †).It should be noted that this is not an in-depth time-trend analysis of the measurement data, and any differences between years in size, age, lipid-content, or sampling season for the sh have not been accounted for.Hence, these results should be interpreted with caution.The results nonetheless agree with monitoring for Norwegian marine sh, showing a decline in contaminant concentrations in the 1990s followed by a stabilization of concentrations in the 2000s. 82Independent of the exact halflives calculated for the measurement data herein, NEM estimates a steeper decline in PCB-153 concentrations in Norwegian marine sh than observed.This was also observed for certain stations in the evaluation of NEM for the atmosphere, including Stórhöfdi (Iceland) and Aspvreten (Sweden). 56ossible explanations for this include (i) uncertainties in the emission estimates, namely an overestimation of the success of efforts to reduce primary PCB emissions, (ii) uncertainties in the estimated fate and persistence of PCB-153 in the physical environment, (iii) uncertainties in estimated bioaccumulation, and (iv) environmental and/or biotic parameters that experienced a temporal trend during the twenty years of sampling, but which the model assumes to be constant.It has been suggested that climate change may perturbate observed time-trends of POPs in the physical environment and biota in the Arctic, 83,84 which may also be the case for Norwegian marine areas in general.As the properties of the physical environment and biota in NEM only varies with seasons, and not over longer time scales, NEM is currently not able to capture this.
The model performance of NEM is comparable to, or better than, the performance of other bioaccumulation models that have been applied to Norwegian marine and Arctic areas.Borgå and Di Guardo 32 used a mechanistic steady-state bioaccumulation model (without spatial resolution) and measured concentrations in water to estimate PCB concentrations in copepods, krill, amphipods and polar cod in the Barents Sea.Estimated concentrations were 20 to 200 times lower than measured, which was attributed to uncertainty in the measured concentrations in water. 32However, observed food chain biomagnication could be reproduced, also in a follow-up study with a model that had been expanded to include the seabird black-legged kittiwake (Rissa tridactyla). 36ikewise, observed bioaccumulation factors (BAF) of PCBs in capelin and cod in the Barents Sea could be reproduced with the dynamic OMEGA model based on concentrations in water as input. 34This study was followed up by two broader studies using the OMEGA model for a range of POPs and marine species in the Svalbard area that estimated 61% and 87% of concentrations within a factor of 5 of measured concentrations, respectively. 35,37lso here, concentrations in water that had been measured, or derived from measured concentrations in air or biota, were used as input.The main strength of NEM compared to these models, is that NEM enables simulation of the complete link from global emissions to biotic exposure, including temporal and spatial variation in both emissions and environmental factors.To our knowledge, this is the rst study to estimate PCB-153 concentrations in marine biota in Norwegian marine areas based on global emissions.Our results demonstrate that NEM performs better in reproducing average measured concentrations, than for individual organisms.This has been observed previously for individual estimates of PCBs in humans made with the ACC-Human model. 45,54,85Individual estimates of PCB-153 in Norwegian women, mothers in the Canadian Arctic, and in American citizens correlated with measured levels with r s of 0.13-0.67,0.17-0.43,and 0.44, respectively.This is in the same range as our results for herring (r s = 0.65) and cod (r s = 0.22).Unreliable dietary intake data at the individual level was identied as a main contributor to poor model performance for estimates of PCB-153 in individual humans. 54,85Similarly, we believe that failing to represent the diversity in the diet of cod is one of the main reasons for the poorer model performance for cod than for herring in this study.Also, one of the main reasons for the good model performance for individual estimates of PCB-153 in herring is likely attributed to the homogeneity of the NSS herring including its diet.

Exploration of emission scenarios
The estimated contribution from Norwegian primary emissions to the total body burden of PCB-153 concentrations in 5 year-old herring was small for the whole time-period and all grid cells For the whole period 1930-2020, the highest contribution from Norwegian primary emissions was estimated to be in the 1960s (Fig. S23 †).Even then, the estimated maximum contribution from Norwegian primary emissions to PCB-153 burden in Norwegian herring was not larger than 3.4% (1963, grid cell 48).This compares to a Norwegian fraction of the total global cumulative PCB-153 emissions (1930-2020) of 0.11% (Fig. S24 †), with a yearly maximum of 0.18% in 1965 (Fig. S25 †).Hence, even if the contribution from Norwegian primary emissions to the body burden of PCB-153 in Norwegian herring is estimated to be small, it is still larger than the Norwegian fraction of global PCB-153 emissions.Nevertheless, this shows that PCB-153 concentrations in sh in Norwegian off-shore marine areas are mainly a result of historical and/or ongoing primary emissions outside of Norway, and not from historical and/or ongoing primary emissions within Norway.This is an important result for Norwegian environmental management, as this means that the regulatory focus needs to be placed on international measures (e.g. the Stockholm Convention) more than national efforts to control and reduce PCB-153 emissions.This conclusion is not valid for locations known to be impacted by local sources including historically contaminated sediment, such as those excluded from this model evaluation (Tables S1 and S5 †).
The estimated contribution of European and Russian primary emissions to the total body burden of PCB-153 concentrations in 5 year-old herring in January 2020 is shown in Fig. 7 and S22, † respectively.The relative contribution of EU emissions to the total PCB-153 load in herring varied with both space (Fig. 7 and S21 †) and time (Fig. S23 †).The contribution of EU primary emissions decreased further north and further east within the model domain (Fig. 7), ranging in 2020 from 31% (cell 55; White Sea) to 74% (cell 60; Skagerrak/Kattegat area) (Fig. S21 †).This compares to a European fraction of the total global cumulative PCB-153 emissions (1930-2020) of 26% (Fig. S24 †), with a yearly maximum of 38% in 1977 (Fig. S25 †).The temporal trend of the relative contribution of European primary emissions for 1930-2020 was similar for all cells (shown for four selected cells in Fig. S23 †).This trend is inuenced by the time-trend of European production of PCBs, which increased from the 1950s (Fig. S25 †). 61Decreasing contribution from European emissions from the 1980s (Fig. S23 †) indicates that emissions ceased faster in the included EU countries than in other countries (Fig. S25 †).This is reected in the temporal trend of the relative contribution of Russian primary emissions, which increases somewhat from the 1980s onward (Fig. S23 and S25 †).Russia was the last country to stop production of PCBs. 61he estimated contribution of Russian primary emissions to PCB-153 body burden in herring is for obvious reasons highest along the Russian coast (White Sea area) (max in 2020: 31%, grid cell 55) and into the Barents Sea, while it is small in the southwestern part of the model domain (min in 2020: 3% in grid cell 60).This compares to a Russian fraction of the total global cumulative PCB-153 emissions (1930-2020) of 4.8% (Fig. S24 †), with a yearly maximum of 13% in 1953 (Fig. S25 †).This means that Russian emissions may be important for e.g., PCB-153 exposure in Northeast Atlantic cod residing in the Barents Sea for large parts of the year, while it is small for sh in the North Sea and Norwegian Sea areas.The contributions of both European and Russian primary emission are smallest in the northwestern part of the model domain, where the total body burden of PCB-153 in 5 year-old herring is mainly explained by emissions from other countries not included in the explored emission scenarios (Fig. S22 †).This makes sense considering that this area is furthest away from both the EU countries, Russia, and Norway, and may be more impacted by primary emissions originating from elsewhere.In grid cell 16, which covers the western part of Svalbard, including King's Bay, PCB-153 concentrations in herring in 2020 are estimated to be the result of primary emissions from "other countries" (49%, i.e. from all countries except Norway, Russia, and EU member states as of 1973), EU (39%), Russia (12%), and Norway (0.4%) (Fig. S21 and S23 †).This clearly shows the importance of global chemical regulations to protect Arctic ecosystems from persistent organic pollutants like PCB-153.

Conclusions
NEM allows for an integrated, quantitative analysis of the whole continuum of processes linking global emissions of contaminants with exposure of organisms in northern ecosystems.We have shown that NEM is able to reproduce medians and ranges of measured PCB-153 concentrations across compartments, including observed variation across space and trophic levels.Further research is needed to investigate why NEM estimates a steeper decline in PCB-153 concentrations in Norwegian sh with time than those observed in the measurement data.We have used the model to explore how the current PCB-153 burden in Atlantic herring from Norwegian marine areas can largely be rationalized by historical primary emissions from outside Norway.This is just an example of how NEM can be applied to explore the relative role of emissions originating from different countries/regions, and to tease out the most important sources for a given contaminant in a specic species or region over time.This is of key importance for policymakers, as information about emission sources is a prerequisite for scientically sound regulatory decision making.
A key advantage of NEM is its mechanistic and integrated nature that allows for scenario testing of not only various emission scenarios, but also evaluation of the estimated effects of changes in other input parameters on estimated fate, bioaccumulation, and concentrations.The model can be used to explore how changes in both the physical environment and the food web are estimated to impact contaminant concentrations, and to tease out to which parameters the model estimates are most sensitive.This is particularly relevant in the context of climate change.Climate change has already been identied as a driver of future changes in the distribution and viability of Norwegian and Arctic marine ecosystems. 23,86,87Climate changeinduced alterations in the physical environment and food webs are also expected to impact contaminant dynamics in Arctic ecosystems, but the extent and direction of it is still largely unknown. 83,84,88,89As NEM is based on both temporally and spatially resolved emission inventories and environmental and food-chain properties (e.g., temperatures, air and water circulation, organic carbon and lipid contents, dietary relationships etc.), it is well-equipped to evaluate concurrent changes in primary emissions and climate-induced changes in environmental and ecosystem characteristics on contaminant dynamics and exposure.Even if our understanding of future climate-induced changes in the physical and biological environment may still be limited, 88 the model can be used for scenario-based hypothesis testing to shed light on processes and parameters that could be important.The model's capability to describe complex climate change scenarios will be restricted to key aspects, such as temperatures or ice cover varying between years/decades, and/or shiing diets and lipid dynamics in the food web.
The NEM bioaccumulation module currently does not allow for animal movement between grid cells, and biota properties are only seasonally and not spatially resolved.Further development is needed to include spatial migrations, differentiated for different species or for different populations within species.It should also be evaluated whether spatial and long-term temporal variation in biota properties should be incorporated in the model.As the physical module of NEM is global in scope, it is possible to apply NEM to other regions of interest.While the NEM bioaccumulation module is currently parameterized for species found in Norwegian marine waters, it is possible to further develop, parameterize, and evaluate it for other species and ecosystems of interest.Moreover, NEM has so far only been evaluated for PCB-153.Further research is needed to evaluate NEM for a range of chemicals with diverse physico-chemical properties, emission scenarios and regulatory status.Based on the good performance of NEM for PCB-153, we are optimistic that the model is also applicable to other POPs and POP-like chemicals.Hence, NEM could represent a future tool for scientists and decision-makers to evaluate the link between emissions and exposure also for chemicals considered for regulation.View Article Online

Fig. 1
Fig. 1 Illustration of the NEM version used in this study.A global environmental fate module is used to run dynamic simulations with 5°×5°lat/long resolution, with subdomains used for the Norwegian marine areas (red) and the Svalbard region (blue) with increased spatial resolution (1°× 1°) highlighted.Concentrations are estimated in three (red) and nine (blue) selected species for the subdomain in the same color.Arrows illustrate feeding relationships.

Fig. 2
Fig. 2 Model evaluation across species in Svalbard for PCB-153 based on 1°× 1°resolution.Median (point) and range (variability bars) of estimated concentrations are plotted against median (point) and range (variability bars) of measured concentrations.M = male, F = female.The solid diagonal line represents a perfect match between estimated and measured concentrations, while dashed diagonal lines show deviations of one order of magnitude.The ranges for the measured concentrations do not include concentrations below the limit of detection.

Fig. 3
Fig. 3 Maps of estimated concentrations of PCB-153 (ng per g lw) in 5-year-old herring and cod in northern European marine areas in January 2020 based on 5°× 5°resolution.The data has been interpolated with the MATLAB 2023a provided script contourf, because it gives a better impression of the concentration distribution (see Fig. S12 †).

Fig. 4
Fig. 4 Model evaluation for individual estimates of PCB-153 in Norwegian spring-spawning herring and Atlantic cod in Norwegian marine areas based on 5°× 5°resolution.The solid diagonal line represents a perfect match between estimated and measured concentrations, while dashed diagonal lines show deviations of one order of magnitude.

Fig. 5
Fig. 5 Comparison of median and range of measured (orange) and estimated (blue) PCB-153 concentrations (ng per g lw) in cod liver of the four different sea regions of the Norwegian marine areas included in the model evaluation.

Fig. 6
Fig. 6 Evaluation of time-trends for concentrations of PCB-153 (ng per g ww) in herring and cod based on 5°× 5°resolution.The black line and the grey area represent the median and range (min-max), respectively, of estimated concentrations for all ages of fish in all grid cells where samples have been collected after 2006 (n = 10 cells for herring, n = 19 cells for cod).The orange markers and their variability bars are median and range (min-max) of the measured PCB-153 concentrations.

Fig. 7
Fig.7Fraction (in percent (%)) of total estimated PCB-153 concentrations in herring muscle in January 2020 estimated to originate from historical or ongoing primary emissions in the EU (member countries as of 1973).The data has been interpolated with the MATLAB 2023a provided script contourf, to give a better impression of the spatial distribution (see Fig.3and S12 †).