Exploring the potential in fl uence of climate change and particulate organic carbon scenarios on the fate of neutral organic contaminants in the Arctic environment †

The main objective of this study is to explore the potential influence of climate change and particulate organic carbon scenarios on the fate of organic chemicals in the Arctic marine environment using an evaluative modeling approach. Particulate organic carbon scenarios are included to represent changes such as enhanced primary production and terrestrial inputs. Simulations are conducted for a set of hypothetical chemicals covering a wide range of partitioning property combinations using a 40-year emission scenario. Differences in model output between the default simulations (i.e. contemporary conditions) and future scenarios during the primary emission phase are limited in magnitude (typically within a factor of two), consistent with other modeling studies. The changes to particulate organic carbon levels in the Arctic Ocean assumed in the simulations exert a relatively important influence for hydrophobic organic chemicals during the primary emission phase, mitigating the potential for exposure via the pelagic food web by reducing freely-dissolved concentrations in the water column. The changes to particulate organic carbon levels are also influential in the secondary emission/depuration phase. The model results illustrate the potential importance of changes to organic carbon levels in the Arctic Ocean and support efforts to improve the understanding of organic carbon cycling and links to climate change.


Introduction
6][17][18] The public availability of global and regional-scale projections of climatic conditions throughout the 21st century (e.g.][25][26][27][28][29] A common strategy adopted for model simulations aiming to explore the potential inuence of GCC on the fate and transport of organic chemicals is to dene a parameter set representing contemporary conditions ('before') and then a modied parameter set incorporating projected changes for a future climate ('aer').Simulations incorporating the temporal trend in projected climate changes over the 21st-century have also been conducted. 28,294][25][26][27][28][29][30] Modeling these legacy contaminants is useful because (i) POPs are of toxicological concern, (ii) monitoring data are oen available to evaluate model performance under contemporary conditions and hence build condence in the reliability of model output under GCC scenarios and (iii) their physical-chemical properties are also relatively well-known, reducing an important source of model input uncertainty.Excluding changes in emissions, model output under the assumed GCC scenarios is typically less than a factor of two different from model output using contemporary (i.e.20th century) conditions 30 and both positive (e.g.reduced contaminant levels) and negative (e.g.elevated contaminant levels) outcomes have been reported, depending on the model application and environmental medium being considered (e.g.atmosphere vs. surface water).
Climate projections indicate that the Arctic environment will experience the greatest warming over the 21st century [19][20][21][22] resulting in profound changes to key elements of the ecosystem, particularly the cryosphere (i.e.lake and sea-ice, permafrost, glaciers).6][37][38] A recent satellite-based study 33 conducted at regional-scale over the period 1998-2009 reported a statistically signicant increase (+20%) in total annual net primary productivity (NPP) in the Arctic Ocean, predominantly due to increases in both the extent of open water and the duration of the open water season.The main factors controlling primary production in the Arctic Ocean are light limitation (e.g. in areas presently covered by multi-year sea ice) and nutrient limitation (particularly in highly stratied waters). 222][33][34] All else being equal, enhanced primary production will lead to higher levels of organic carbon (OC) in the water column. 27][41] The link between primary production in aquatic environments, (particulate) organic carbon levels and exposure to hydrophobic organic contaminants is re-emerging as a topic of interest.For example, based on trends in Shexaand Shepta-PCB levels in burbot (Lota lota) and sedimentary organic matter over the period 1985-2009, it was recently hypothesized that increased primary production (related to climate change) corresponds to enhanced exposure potential to hydrophobic contaminants. 123][44][45][46] Concentrations in pelagic biota can be lower in eutrophic systems even if the total amount of chemical in the abiotic environment is higher compared to less productive systems. 45The implication of eutrophication for benthic infauna depends on how organic carbon normalized concentrations in the active sediment layer are inuenced as well as whether respiratory exchange occurs primarily with sediment pore water or overlying water.While organic carbon-normalized concentrations of SPCBs in surface sediments showed no relationship with trophic status in a survey of 19 Swedish lakes located in close proximity, 45 these ndings cannot be assumed to be universal.
To date, global-scale modeling studies seeking to explore the potential inuence of GCC on the fate and transport of organic chemicals have not included changes to organic carbon cycling in the Arctic environment.Accordingly, the main objective of the current study is to explore the potential inuence of particulate organic carbon (POC) scenarios on the fate of neutral organic contaminants in combination with other GCC projections.Long-range transport is assumed to be the only source of contaminants entering the Arctic environment.Model output is generated over a 40-year period for hypothetical chemicals covering a range of partitioning properties (i.e.airwater partition coefficient, K AW ; octanol-water partition coefficient, K OW ; octanol-air partition coefficient, K OA ) and susceptibility to degradation in the environment.The modeling approach is evaluative in nature but still provides useful insights into the potential inuence of GCC scenarios on the behaviour of organic contaminants in the environment.As exposure to organic contaminants in the Arctic via the marine food web remains a particular concern, 47,48 the results are discussed primarily in the context of exposure to neutral organic contaminants in the marine environment.

Global-scale fate and transport model (GloboPOP)
All simulations were conducted with the GloboPOP model, 49,50 a fugacity-based chemical fate and transport model.Important features of this model are presented in Section S1 of the ESI.† The globe is divided into ten latitudinal bands, each with wellmixed compartments representing the atmosphere (four layers), surface ocean water (0-200 m depth) and the terrestrial environment which includes freshwater and underlying sediments, different soil types (e.g.agricultural, uncultivated, forest) and forest canopies (deciduous, coniferous).Seasonal snow packs are included in the North Temperature, North Boreal and North Polar zones (ESI, Section S1 †) whereas the South Polar zone has a permanent snow pack.Sea-ice cover (% of water surface area covered) is tracked for the North Polar, South Polar and South Sub-Polar zones.Details on the treatment of snow (as precipitation and snow pack) and other recent modications to the model are presented elsewhere. 51One advantage of using the current version of GloboPOP is that shelf sediments in the Arctic Ocean are explicitly included.This compartment is an important reservoir of hydrophobic organic contaminants in the Arctic environment and can act as major secondary source once primary emissions are exhausted. 51

Global climate change (GCC) scenarios
The GCC scenarios implemented for the current study reect changes in (i) air and ocean water temperatures, (ii) precipitation rates, and (iii) seasonal sea-ice cover (and ice-free periods for freshwater systems).General features of the GCC scenarios are summarized in Table 1 whereas specic details are presented in Section S2 of the ESI.† Given the coarse spatial resolution of the model and limited inuence reported in previous global-scale modeling exercises, 26,29 projected changes to atmospheric and oceanic circulation patterns are not considered.Note that all changes related to the GCC scenarios are introduced at the beginning of the simulation as opposed to being introduced incrementally over time.[21][22] Particulate organic carbon (POC) scenarios Levels of particulate organic carbon (POC) in surface ocean waters vary substantially in different regions of the Arctic Ocean due to factors such as the intensity and duration of primary production, distance from shore/intensity of inputs from terrestrial sources (e.g.riverine inow, coastal erosion) and the potential for resuspension from underlying sediment beds. 52his variability is reected in the ranges of values for POC levels in surface waters reported for a recent sampling campaign (ISSS-08 Cruise). 53For example, measured ranges were 6.5-17.8mM in the Barents Sea, 3.4-33.0mM in the Kara Sea, 1.2-151.9mM in the Laptev Sea, 1.0-19.4mM in the W-East Siberian Sea, 2.0-36.6 mM in the E-East Siberian Sea and 3.1-33.1 mM in the Chukchi Sea.Measured values from other studies range from approximately 10-30 mM in the Barents Sea, 54 1-10 mM in the Canada Basin 55,56 and <1 mM reported for the central Arctic. 57s the coarse spatial resolution of the GloboPOP model does not allow regional differences in POC levels to be represented, simulations for the default GloboPOP simulations and the GCC scenarios (Table 1) were conducted assuming average POC levels in Arctic surface waters of 20 mg L À1 ($1.5 mM) and 200 mg L À1 ($15 mM).
Changes in primary production (PP) projected for the Arctic Ocean described in the Arctic Climate Impact Assessment (ACIA) 22 are as follows; (i) light-limited regions: increased 2-5-fold, due to sea-ice disappearance, assuming no nutrient limitation; (ii) Barents Sea: increased $2-fold due to deeper wind-mixed layer and supply of nutrients from underlying Atlantic water; (iii) Arctic shelf seas: increased $2-fold due to retreat of permanent ice pack beyond shelf break and subsequent onset of enhanced upwelling and exchange with nutrientrich offshore waters.Recent studies incorporating model simulations suggest far more modest changes.For example, Ellingsen et al. 31 applied a coupled hydrodynamic-ecological model for the Barents Sea and reported only an 8% increase in PP over the next 65 years.Similarly, Lavoie et al. 32 employed a coupled sea ice-ocean-biological 1D model to project future PP  Table S1-S3 (ii) Altered precipitation rates (up to AE35%) Table S4 and  for the Mackenzie Shelf (Beaufort Sea); projected increases in average annual primary production are low (increased 6 and 9% for 2050 and 2090 respectively vs. present day), primarily due to the inuence of enhanced fresh water input (intensied strati-cation, reduced replenishment of nutrients).On the other hand, as noted above, a statistically signicant increase (+20%) in total annual net primary productivity (NPP) in the Arctic Ocean has already been reported for the period 1998-2009. 33hile primary productivity per unit area per unit time is similar, total annual NPP increased because both area for primary production and growth period have increased.
Assuming no nutrient limitations, the authors suggested that total annual NPP could increase nearly 2-fold by mid-century, when summer minimum ice coverage approaches zero.In agreement with this estimation, Slagstad et al. 34 reported increases in gross and net primary productivity in the range of approximately 1.1-4-fold based on model simulations of different regions of the Arctic Ocean assuming ice-free conditions in summer.
Based on the considerations presented above, the volume fraction of POC in the water column of the Arctic surface ocean was increased 2.5-fold in the GCC + OC scenario.In other words, the default values for the concentration of POC in the water column (C POC ) are multiplied by a factor of 2.5 (i.e.C POC ¼ 50 and 500 mg L À1 ).To be consistent with this change, the mass transport coefficient (MTC) characterizing deposition (i.e.settling) of particle-bound contaminant was increased by the same factor.Additionally, the organic carbon content of the shelf sediments (SOC) was also increased 2.5-fold.Changes in the parameterization of the freshwater compartment (see Table 1) are discussed briey in the ESI (Section S2 †).

Model application
A total of six scenarios were simulated for the current study (see Table 1 and ESI Section S2 †).Simulations were rst conducted using the default climate parameterization assuming a C POC in the Arctic Ocean surface water of 20 mg L À1 (see above).The simulations were then repeated assuming a C POC of 200 mg L À1 .Model outputs were then generated for the GCC scenario (Table 1) under the same two assumptions regarding POC (and SOC) in the Arctic environment.Next, model output was generated for the GCC + OC scenario (Table 1), where, for example, concentrations of POC in the Arctic Ocean surface water are 2.5-fold higher (i.e.C POC ¼ 50 mg L À1 and then 500 mg L À1 ).For each scenario, the model was run for a set of hypothetical chemicals covering a wide range of partitioning property values (log K OA ¼ 4 to 15, log K AW ¼ 3 to À5 where log K OW ¼ log K OA À log K AW and only log K OW values from 1 to 12 are included).All chemicals were assumed to be relatively persistent in the environment and degrade with the following characteristics; (i) a second-order reaction rate constant with atmospheric OH radicals (k OH ) at 25 C of 1 Â 10 À12 cm 3 per molecule per s and (ii) degradation half-lives at 25 C of 1, 2 and 10 years in water (fresh, marine), soils (all types) and sediments (fresh, marine) respectively.Degradation half-lives in snow pack and forest canopies were assumed to be equal to those in water.All All chemicals were emitted into the environment using the same 40-year emission scenario which species a constant (unit) emission rate for 20 years followed by zero emissions for the remaining 20 years.This approach allows model output to be generated and compared during the primary emission phase (Year 0-20) and as secondary sources (e.g.environmental reservoirs such as boreal forest soils) come to dominate (Year 20-40).Over the latter period, levels in the environment are declining due to degradation and export out of the model domain (e.g.sediment burial, export to deep oceans).The rate at which levels dissipate is specic to the environmental compartment and model zone however.Emissions were directed 100% to the atmosphere and the spatial distribution of emissions is based on estimates for PCBs. 59This pattern is assumed to represent a typical industrial chemical (i.e.emitted predominantly in the northern Hemisphere at mid-latitudes) (see Section S4 of the ESI †).Emissions to the North Polar zone of the model were purposely set to zero to ensure that all mass entering the Arctic environment is due to long-range transport.The model calculations are solved numerically with a time step ranging from 1-12 h, depending on the partitioning properties of the hypothetical chemical.

Assessment metrics
The assessment of the potential inuence of GCC scenarios on the behavior of neutral organic contaminants in the Arctic environment during the primary emission phase is based on the values of the following metrics at Year 20, (i) Arctic Contamination Potential (eACP 20 ), 60 (ii) average concentration in the surface layer of the atmosphere (C A ) and (iii) average freelydissolved concentration in the surface ocean (C O ). eACP 20 is calculated from the model output as: where M T is the total mass of chemical in the North Polar zone of the model, M A is the total mass of chemical in the atmosphere of the North Polar zone and E G is the total mass of chemical emitted to the global environment.Model output from all GCC scenarios (Table 1) is compared to default model output for each partitioning property combination and then the ratios (i.e., GCC:Default, GCC + OC:Default) displayed as a function of log K OA and log K AW .Ratios greater than one mean that the GCC scenario has resulted in elevated contamination of the Arctic environment compared to the default scenario.We chose to focus on eACP 20 , C A and C O because these outputs represent long-range transport potential, a frequently sampled abiotic medium in the Arctic and the most relevant consideration for exposure to contaminants via the marine pelagic food web (e.g.water / phytoplankton / zooplankton / Arctic cod / ringed seal / human/polar bear) respectively.Patterns of response in the atmosphere are also of interest because atmospheric concentrations of some POPs in the Arctic (e.g.HCHs) have already been suggested to be responding to climate forcings. 61Priority is given to the marine pelagic food web because consumption of seal blubber is associated with elevated body burdens of POPs in top predators (i.e.humans and polar bears). 15,16 different metric is used to assess the potential inuence of GCC scenarios on the fate of neutral organic contaminants in the Arctic environment aer primary emissions have ceased.In this case, the ratio between model output at Year 40 and Year 20 for each scenario is calculated (R) and then the assessment metric (M 40 ) is based on the comparison between the Year 40 : 20 ratios for the GCC and default scenarios, e.g.
where C A,40,DEF or GCC is the average surface air concentration at Year 40 and C A,20,DEF or GCC is the average surface air concentration at Year 20.This approach was adopted because the current study does not sequentially introduce changes in model parameterization during the simulation and the modeled concentrations (and masses) at Year 20 are not identical between the different scenarios.The M 40 values are surrogates for apparent dissipation half-life and values greater than one mean that the GCC scenario has resulted in a slower rate of depuration compared to the default scenario (e.g. the time required to reduce ambient concentrations in the Arctic environment by a factor of two is longer).Other interpretations of this assessment metric are discussed in the following section.An illustration of the temporal trend in emissions, representative model outputs under default and GCC scenarios and how model output at Year 20 and 40 is used to generate the assessment metrics is presented in Fig. 1.Note that the metrics at Year 20 and Year 40 are not directly comparable in terms of the magnitude of change.All results of the present set of model simulations are considered evaluative in nature (i.e.best suited to explore sensitivities and trends at a screening level).

Results and discussion
Primary emissions phase (Year 0-20) A comparison of model output from the GCC scenarios (GCC, GCC + OC) to the default scenario in Year 20 for the selected model outputs (eACP 20 , C A , C O ) is presented in Fig. 2.These results are for the scenarios assuming default POC concentrations in Arctic Ocean surface water of 200 mg L À1 .A comparison of the default and GCC + OC scenario assuming POC levels of 20 mg L À1 is presented in Section S5 of the ESI.† All model outputs presented in Fig. 2 were generated assuming the longer degradation half-lives (i.e.slower degradation rates, see above).
Results obtained using the shorter degradation half-lives are broadly similar to those shown in Fig. 2 and hence are not discussed further.The general patterns in long-range transport potential and accumulation in the Arctic environment for the hypothetical chemicals considered are summarized in Section S6 of the ESI.† Key features of the results presented in Fig. 2 are as follows.First, regardless of the model output considered, the changes caused by implementation of the GCC scenarios (Table 1) typically result in values within a factor of two of the default scenario.For example, the Arctic Contamination Potential (eACP 20 ) is relatively insensitive to both GCC scenarios across a wide range of partitioning property combinations, exhibiting changes in the range of AE20%.Such results stem from the fact that changes implemented for the GCC scenarios can act antagonistically (i.e.inuences are competing and so tend to balance out).Furthermore, some changes implemented in the GCC scenarios (e.g.precipitation rates) are limited in magnitude.The eACP 20 s of volatile chemicals (log K AW $ 1) with log K OA # 7 are more sensitive to the GCC scenarios (i.e.values reduced to 50% or less of default model output) but such chemicals tend to distribute into the gas phase in the atmosphere with only a relatively small fraction present in surface media (ESI, Section S6 †).Higher temperatures under the GCC scenarios are particularly effective in reducing net deposition for these chemicals, partly because of the non-linear relationship between temperature and partitioning property values.
Secondly, increases in average air concentration in Year 20 (up to +40%) are observed under both GCC scenarios for chemicals with log K AW $ À4.5 and log K OA between approximately 7 and 10.5.This behaviour reects alterations in longrange transport efficiency of such contaminants to the Arctic

Paper
Environmental Science: Processes & Impacts environment and changes in the net deposition uxes and accumulation in surface media (e.g.soils and shelf sediments).
Note that this range of partitioning properties overlaps with known POPs such as PCBs.The suggestion that atmospheric concentrations of chemicals with such partitioning properties may increase in direct response to higher air temperatures is consistent with previous global-scale modeling studies. 26,29owever, it is also clear that the responses to the GCC scenarios seen in the surface compartments are different from the responses seen in the atmosphere (e.g.compare the eACP 20 and C O plots for the GCC scenario to the output for C A in Fig. 2).Accordingly, the direction (and magnitude) of the response to GCC scenarios predicted for the atmosphere cannot be used to predict responses in other compartments or speculate directly on the implications for exposure to organic contaminants via dietary uptake.
The most important feature of the results presented in Fig. 2 in the context of exposure to neutral organic contaminants is that C O is reduced over a much wider range of partitioning property combinations in the GCC + OC scenario compared to the GCC scenario.For chemicals with log K OW $ 5.5 and log K OA > 7, the response is roughly in proportion to the extent that the volume fraction of POC in the water column is assumed to increase in the GCC + OC scenario (i.e., a factor of 2.5; Table 1).The model output for C O is driven by a mass dilution effect (i.e. as the volume of particulate OC is increased, less chemical mass is present in the dissolved phase) as well as enhanced deposition of particulate-bound chemical out of the surface water to the shelf sediments and deep ocean.While atmospheric deposition becomes more favorable thermodynamically, such inputs are not sufficient to counter the enhanced losses associated with particle settling and consequently, OC-normalized concentrations in suspended solids are also reduced.These results are consistent with past studies on the inuence of eutrophication and contaminant dynamics in aquatic systems [42][43][44][45] and also with recent empirical data characterizing the temporal trends in bioaccumulation associated with the rapid changes in primary production observed in an alpine lake. 46The decrease in C O is expected to lead to a concomitant reduction in exposure at the base of the pelagic food web, which then propagates upwards. 27It is also worth noting that organic carbon-normalized concentrations in the shelf sediments of the Arctic Ocean (active layer only, i.e. top 5 cm) are reduced in the GCC + OC scenario across a similar range of partitioning properties that exhibit reduced freely-dissolved concentrations (results not shown).Lower organic carbon-normalized concentrations in the shelf sediments imply reduced exposure potential in benthic organisms via respiration of pore-water and ingestion of sediment OC.As shown in the ESI (Section S5 †), the comparison between the default model output assuming C POC ¼ 20 mg L À1 and the GCC + OC scenario is broadly similar to the results shown in Fig. 2 (default C POC ¼ 200 mg L À1 ).However, a shi occurs in the threshold of response for chemicals with log K OW between 5 and 8 and log K AW < À1.Essentially, when the default POC concentrations are reduced 10-fold, K OW must be increased 10-fold to achieve the same reduction in freely-dissolved water concentration.Because of the heterogeneity in POC levels in the Arctic Ocean, the implication of these ndings is that changes in the exposure potential for hydrophobic chemicals may not be uniform across the region, even given the same proportional change in POC levels.

Secondary source/depuration phase (Year 20-40)
A comparison of model output from the GCC scenarios (GCC, GCC + OC) to the default scenario in Year 40 versus Year 20 (see Fig. 2 Comparison of model output from the GCC scenarios (GCC, GCC + OC) to the default scenario in Year 20 for (i) Arctic Contamination Potential (eACP 20 ), (ii) average concentration in surface air (C A ) and (iii) average freely-dissolved concentration in surface ocean water (C O ).Ratios greater than 1.0 (+) mean that model output under the GCC scenario is higher than output from the default scenario.The dotted line indicates a transition from ratios less than 1.0 to those greater than 1.0.The approximate regions of chemical partitioning space occupied by PCBs and HCHs are also indicated (see GCC, C A panel).eqn (2)) for average concentration in surface air (C A ) and average freely-dissolved concentration in surface ocean water (C O ) is presented in Fig. 3.These comparisons represent the changes in the apparent dissipation half-life of the chemical from the Arctic environment as a whole and in the atmospheric boundary layer and surface ocean water.As in Fig. 2, all model outputs presented in Fig. 3 were generated assuming the longer degradation half-lives.Also, as before, results obtained using the shorter degradation half-lives are broadly similar to those shown in Fig. 3 and are not explicitly discussed.
The most obvious feature of the results presented in Fig. 3 is that the apparent dissipation half-life in the atmospheric boundary layer (C A ) is greater under the GCC scenarios compared to the default conditions (i.e.ratio greater than one; see eqn ( 2)) across a wide range of partitioning property combinations.For all such chemicals, the reduction in air concentrations is lower in relative terms over the last 20 years of the simulation.The highest sensitivity in the atmospheric boundary layer (C A ) under the GCC scenarios is seen for chemicals with log K OA > 10.5.Enhanced revolatilization from secondary sources is an important factor driving the results for such chemicals.This phenomenon is also important for other partitioning property combinations, including those overlapping known POPs.The results follow intuitively from the higher temperatures and reduced period of time that transport from soils and the surface ocean water to the atmosphere is impeded by snow pack and sea-ice respectively in the GCC scenarios.There are no important differences between the GCC and GCC + OC scenarios.The inuence of the GCC scenarios on temporal trends in atmospheric levels during the depuration phase (Year 20-40) is consistent with a recent analysis of POP air monitoring data from the Arctic (i.e.slowing decline in atmospheric levels due to enhanced inputs from secondary reservoirs). 61However, the changes implemented in the GCC scenarios are more extreme than the changes observed in the Arctic environment to date (e.g.sea-ice cover reduced to negligible fraction of Arctic Ocean in the GCC scenarios) and hence the model results should not be over-interpreted when discussing trends in contemporary monitoring data.
The other feature of the results presented in Fig. 3 worth discussing is that response patterns observed in the atmosphere (C A ) are different from the pattern of those observed in the freely-dissolved concentrations in surface ocean water (C O ).Furthermore, the responses in C O are substantially different between the GCC and GCC + OC scenarios for chemicals with log K OW $ 5. Generally speaking, the increases in the apparent dissipation half-lives in GCC + OC scenario are more pronounced (i.e. higher sensitivity).The direction of the response in the C O metric in the depuration phase (Fig. 3) is also opposite to what was observed during the primary emission phase (Fig. 2).These patterns can be explained by considering the relative importance of sediment resuspension as a source term for the surface ocean over the simulation period.During the primary emission phase (Year 0-20), atmospheric deposition is more important than sediment resuspension in terms of driving levels in the surface ocean and hence alterations in sediment-water exchange (i.e.enhanced mass transfer coefficient, MTC, for particulate deposition) are not inuential.However, resuspension from shelf sediments increasingly dominates aer Year 20 and hence assumptions regarding the MTCs for sediment deposition, resuspension and burial become important.In these simulations, sediment resuspension is assumed to be a constant fraction of the sediment deposition ux and therefore any increase in the sediment deposition MTC (Table 1) is matched by a proportional increase in the resuspension MTC (i.e. shelf sediments are more efficient at resupplying the water column in the GCC + OC scenario).The greater apparent dissipation half-lives seen in the surface ocean water are primarily due to this model parameterization approach.Such a response is plausible in reality though, if, for example, enhanced inputs in the water column gradually lead to enhanced OC content in sediments and resuspended particles (i.e. more chemical could be delivered to the overlying water column per unit time).Although these results imply that declines from peak exposure levels in the pelagic food web are delayed as a consequence of GCC, no specic comments on the toxicological implications can be made given the generic nature of the current model application.

Limitations and future considerations
The results presented in Fig. 2 and 3 are broadly consistent with expectations based on previous global-scale modeling characterizing the potential direct inuence of GCC on chemical fate and transport 26,29,30 in terms of magnitude of response.These results lend support to the assertion that indirect effects and behavioural responses to GCC (e.g.dietary shis, northward shis in industrial or agricultural activity) may be more important considerations overall compared to direct effects of GCC. 10,30However, it is also important to recognize the limitations inherent to the evaluative modeling approach employed and the scope of the GCC scenarios considered in the present study.For example, potential GCC-related changes in emission strength and spatial distribution over time and the response to GCC scenarios of chemicals exhibiting different modes of entry into the environment (e.g. for a pesticide emitted 95% emission to soil, 5% to freshwater as opposed to 100% to the air) were not simulated.The simulations in the present study were also conducted using degradation half-lives and activation energies typical of most known POPs and hence are not broadly representative of the entire realm of chemical space.Some of these aspects were explored in a recent global-scale modeling exercise. 29As the global fate of chemicals with lower environmental persistence than assumed here increasingly becomes dominated by degradation kinetics, higher temperatures under GCC scenarios increasingly result in reduced exposure potential relative to the default scenario (i.e. as T increases, mass of chemical in environment decreases). 29As the selected emission pattern (mode of entry, spatial distribution) and temperaturedependencies of partitioning properties and reactivity are most representative of industrial chemicals like PCBs, caution is necessary when extrapolating the results to chemicals that deviate substantially from these assumptions.
It is also worth reiterating that changes in atmospheric and oceanic circulation patterns were not considered in the present study.While changes in long-term (average) atmospheric circulation patterns do not appear to have an important inuence on long-range transport potential, 26,29 there are still valid reasons to explore these aspects using more appropriate (i.e.spatially-resolved) tools.For example, the potential implications of changes in the frequency and intensity of positive NAO events (NAO+), conditions which favour long-range transport of chemicals to the Arctic from European and North American sources, 24,62 could be simulated.The potential short-term and long-term inuence of altered frequencies of episodic transport events and extreme weather conditions could also be assessed.Finally, although the available empirical data do not generally support the hypothesis that sea-ice is an important mass reservoir and transport pathway for the organic compounds studied to date or that food webs intimately linked with sea-ice (i.e.including epontic organisms) have higher exposure potential than other pelagic food webs, [63][64][65] a more detailed treatment of this compartment in the model is required to fully assess the potential implications of GCC for exposure in the Arctic environment.The role of glacier melt in the Arctic (e.g.input of stored chemical vs. dilution through melt of pre-industrial ice) also requires further consideration. 2,30

Conclusions
The different model results seen in the atmospheric boundary layer (C A ) compared to the surface ocean (C O ) demonstrate that responses to GCC are best evaluated from a multimedia perspective, ideally using the most exposurerelevant metrics (e.g.freely-dissolved as opposed to total water concentrations, organic carbon-normalized as opposed to total sediment concentrations).The differences in model output between the GCC and GCC + OC scenarios illustrate that enhanced particulate OC in the Arctic marine environment may be one of the more important considerations for hydrophobic chemicals.During the primary emission phase, enhanced POC in the Arctic marine environment can exert a mitigating inuence on exposure to hydrophobic chemicals.Potential changes in POC cycling and sediment-water exchange dynamics linked to GCC are therefore key processes to consider further, particularly in the context of determining if responses in pelagic food webs are representative of responses in benthos (e.g.benthic infauna inhabiting shelf sediments).Unfortunately, projecting changes to inputs of organic carbon to the Arctic environment due to enhanced primary production, permafrost melt, coastal erosion and other alterations to organic carbon cycling is challenging and subject to many data gaps. 66][69][70] For example, the extent to which mineralization of POC occurs in the water column was hypothesized to exert an important inuence on how chemicals distribute between suspended solids, sediments and water. 71For legacy pollutants such as PCBs, the contaminant levels associated with terrestrial OC sources in the Arctic may also be important consideration.Examples of research that could provide useful insights in this context include (i) contaminant concentration vs. depth proles and sorption/desorption experiments with permafrost and eroded coastline soils/ sediments, (ii) characterization of the susceptibility to degradation of various organic carbon sources in the Arctic environment and (iii) development and implementation of an improved parameterization of the organic carbon mass balance in the Arctic Ocean for contemporary and potential conditions in the future.More explicit treatment of the fate of organic carbon in shelf sediments is also worth considering.Given the many uncertainties related to the potential for enhanced levels of OC in the Arctic marine environment to inuence exposure to neutral organic chemicals, further empirical studies and model development aimed to clarify these aspects are warranted.View Article Online Altered temperatures in atmosphere/ terrestrial environment and surface ocean (up to + 8-10 C in Arctic, +3-4 C elsewhere)

5 ( 5 a
Fig. S2 (v) Y snow pack melt period in Spring À10 days GCC + OC GCC (see above)+ ESI, Section S2 (i) [ volume fraction of POC in surface ocean water b Â2.5 simulations were repeated assuming 10-fold slower reaction kinetics, i.e. a k OH at 25 C of 1 Â 10 À13 cm 3 per molecule per s and degradation half-lives at 25 C of 10, 20 and 100 years in water (marine), soils (all types) and sediments.Taken together, this range of degradation half-lives is representative of many POPs.Note that the two k OH values listed above correspond to a degradation half-life in the atmosphere of approximately 8 and 80 days respectively at 25 C, assuming a global average atmospheric OH radical concentration 58 of 9.8 Â 10 5 molecules cm À3 .However, atmospheric degradation half-lives in each model zone are determined on a daily basis following the seasonal patterns in temperature and OH radical concentrations.Assumptions regarding the temperature-dependence of the partitioning constants and degradation kinetics are presented in Section S3 of the ESI.†

Fig. 1
Fig. 1 Plot illustrating the temporal trend of the emission scenario, representative model outputs under default and GCC scenarios and how model output at Year 20 and 40 is used to generate the assessment metrics.The solid line represents model output under the default parameterization.The dotted lines represent model output under the GCC scenarios (which can result in higher or lower values compared to the default scenario, depending on the chemical properties).

Fig. 3
Fig. 3 Comparison of model output from the GCC scenarios (GCC, GCC + OC) to the default scenario in Year 40 versus Year 20 for (i) average concentration in surface air (C A ) and (ii) average freely-dissolved concentration in surface ocean water (C O ).Ratios greater than 1.0 (+) mean that the apparent dissipation half-life is greater in the GCC scenario than the default scenario.The dotted line indicates the transition from ratios less than 1.0 to those greater than 1.0.The approximate regions of chemical partitioning space occupied by PCBs and HCHs are also indicated (see GCC, C A panel).
Sciences and Engineering Research Council of Canada are acknowledged for funding support.Paper Environmental Science: Processes & Impacts Open Access Article.Published on 18 October 2013.Downloaded on 7/19/2018 7:27:30 PM.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.

Table 1
Summary of the GCC scenario parameterization implemented for the current study compared to the default GloboPOP parameterization a