Assessment of the sea-ice carbon pump : Insights from a three-dimensional ocean-sea-ice-biogeochemical model ( MPIOM / HAMOCC )

1Max Planck Institute for Meteorology, Hamburg, Germany 2Department of Biology, University of Southern Denmark, Odense, Denmark 3Nordic Center for Earth Evolution, University of Southern Denmark, Odense, Denmark 4Scottish Association for Marine Science, Oban, United Kingdom 5Arctic Research Centre, Aarhus University, Aarhus, Denmark 6Centre for Earth Observation Science, Department of Environment and Geography, University of Manitoba, Winnipeg, Canada 7Department of Geological Sciences, University of Manitoba, Winnipeg, Canada 8Greenland Climate Research Centre, Greenland Institute of Natural Resources, Nuuk, Greenland


Introduction
The magnitude of CO 2 exchange between the atmosphere and the ocean is determined by the difference in partial pressure of CO 2 (pCO 2 ) between the two reservoirs.This difference is affected by at least three oceanic pumps that carry carbon from the ocean surface to greater depths: the solubility pump, the biological pump, and the carbonate counter pump (e.g., Sarmiento and Gruber, 2006).In addition to these three carbon pumps, recent studies have suggested that processes related to sea-ice growth and melt may also affect the air-sea CO 2 exchange (Rysgaard et al., 2007(Rysgaard et al., , 2009(Rysgaard et al., , 2011;;Dieckmann et al., 2010;Loose et al., 2011).However, a robust quantification of this so-called sea-ice carbon pump (SICP) has been lacking so far.This study quantifies the contribution of the SICP to air-sea CO 2 exchange in the state-of-the-art global ocean-sea-ice biogeochemistry model MPIOM/HAMOCC.
The basic functioning of the SICP can be understood as follows.During sea-ice growth, all dissolved tracers, including dissolved inorganic carbon (DIC) and the total alkalinity (TA, the acid neutralization capacity of seawater), are concentrated in the brine.Hence, when drainage of brine into the ocean occurs, it reduces the content of tracers in the bulk sea ice and increases their content in the underlying ocean.Simultaneously, the precipitation of the carbonate crystal ikaite (CaCO 3 (s) • 6H 2 O), can occur within the sea-ice brine, as observed in both natural and artificial sea ice (Dieckmann et al., 2008(Dieckmann et al., , 2010;;Rysgaard et al., 2012Rysgaard et al., , 2013Rysgaard et al., , 2014;;Geilfus et al., 2013).Because one mole of DIC is consumed and TA decreases by two moles for each mole of ikaite formed, ikaite crystals have a TA:DIC ratio of 2:1 (for details on the carbonate chemistry see, e.g., Sarmiento and Gruber, 2006).If the formed ikaite crystals remain in the ice matrix while some of the CO 2 (aq) produced during ikaite formation is rejected via brine drainage to the underlying ocean, then a TA:DIC ratio between 1 and 2 remains in the bulk sea ice.This ratio can be further modified through other biogeochemical processes such as the growth of sea-ice algae.During seaice melt, the release of melt-water to the ocean dilutes the oceanic TA and DIC concentrations, and the dissolution of ikaite potentially increases the oceanic TA:DIC ratio which is usually 1.1.These changes in the surface ocean TA and DIC affect the oceanic pCO 2 and determine the direction and strength of the SICP-induced air-sea CO 2 flux (Figure 1).If, on average, the polar mixed layer becomes enriched in brine-associated TA and DIC of potentially low TA:DIC ratio, then a SICP-induced flux of CO 2 to the atmosphere can be expected.Such CO 2 flux to the atmosphere would occur, for example, in regions with a net export of sea ice or where the oceanic circulation causes a net import of brine-associated TA and DIC to the mixed layer.On the other hand, if melt-water dilutes the DIC and TA concentrations in the mixed layer potentially accompanied by an increase in the mixed layer TA:DIC ratio due to the dissolution of ikaite, then a SICP-induced oceanic uptake of atmospheric CO 2 can be expected.Such CO 2 uptake by the ocean would occur, for example, in regions with a net import of sea ice or where the oceanic circulation transports brine-associated TA and DIC away from the mixed layer.Hence, the total SICP-induced oceanic CO 2 uptake will increase with the fraction of brine being exported from the polar mixed layer (Rysgaard et al., 2007).
Budget calculations of the potential size of the SICP-induced annual mean CO 2 flux have suggested a possible maximum oceanic uptake of 9 to 20% (14 to 31 Tg C yr −1 ) of the current net uptake in the region north of 62°N, and 42 to 116% (19 to 52 Tg C yr −1 ) of the current net uptake in the region south of 50°S (Rysgaard et al., 2011).Including an idealized representation of the SICP in our prognostic three-dimensional global ocean-sea-ice-biogeochemical model allows us to evaluate the robustness of these estimates in a more comprehensive setup.
To do so, we analyzed the impact of sea-ice thermodynamics and dynamics and of oceanic circulation on the magnitude of the SICP-induced air-sea CO 2 exchange.We considered two idealized equilibrium climate scenarios with fixed atmospheric CO 2 concentrations: (i) a climate state representing the present-day oceanic and sea-ice state; and (ii) a climate with reduced Arctic sea-ice extent and volume and reduced strength of  the oceanic overturning circulation.These scenarios allowed us to disentangle the various contributors to the overall magnitude of the SICP.Because of the high spatial and temporal variability in the observed sea-ice TA and DIC concentrations, we ran several simulations for each climate scenario using different prescribed ratios and concentrations of TA and DIC in the sea ice that cover the observed range.This study complements a related study that is based on the model NEMO-LIM-PISCES and quantifies the SICP under pre-industrial climate conditions (Moreau et al., 2016).

Ocean-sea-ice-biogeochemical model
Our model simulations were performed using the Max Planck Institute global ocean-sea-ice general circulation model MPIOM (Marsland et al., 2003;Jungclaus et al., 2013) coupled to the Hamburg ocean carbon cycle model HAMOCC (Maier-Reimer et al., 2005;Ilyina et al., 2013).MPIOM is a free surface model with hydrostratic and Boussinesq assumptions.The model computes isopycnal diffusion of the thermohaline fields.Tracer transport by unresolved eddies is parameterized following Gent et al. (1995).MPIOM includes a sea-ice model (Marsland et al., 2003;Notz et al., 2013) formulated with a viscous-plastic rheology following Hibler (1979).MPIOM is formulated on an orthogonal C-grid in the horizontal and z-grid in the vertical.In the setup used in this study one pole is located over Greenland and the other over Antarctica.The horizontal resolution (GR15) is 1.5°, which corresponds to about 12 km around the poles and 150 km in the tropics.The time-step of the model is set to 1.2 hours.The water column is divided into 40 levels in the vertical, 9 of which are in the uppermost 100 m, and 18 in the uppermost 500 m.
HAMOCC is a state-of-the-art ocean carbon cycle model resolving oceanic and sedimentary organic and inorganic carbon processes and therefore includes the biological pump, the carbonate counter pump, and the solubility pump (Six and Maier-Reimer, 1996;Heinze and Maier-Reimer, 1999;Heinze et al., 1999;Ilyina et al., 2013).The temporal evolution of TA and DIC is computed prognostically.All biogeochemical tracers of HAMOCC are transported with the same advection scheme as the thermohaline fields of MPIOM.
The air-sea CO 2 flux, F, in HAMOCC is derived at each time step by computing the pCO 2 difference between the atmosphere and the uppermost oceanic layer, ΔpCO 2, multiplied by a gas transfer velocity, k, the solubility of CO 2, α, according to Weiss (1974), and the open-water fraction, A op , as gas exchange is assumed to occur only in the open-water part of the grid cell.Hence, where k depends on the Schmidt number (Gröger and Mikolajewicz, 2011) and wind speed at the surface (Wanninkhof, 1992).The oceanic pCO 2 in the model is computed as a function of temperature, salinity, DIC, TA, and sea-level pressure.More details on the biogeochemical processes simulated in HAMOCC are given by Six and Maier-Reimer (1996) and Ilyina et al. (2013).

Treatment of brine and tracer release
To simulate the SICP including the formation of ikaite, we prescribed the mass concentration of dissolved tracers (solute) in sea ice (e.g., salinity, TA, DIC) as a constant value C sol,ice [g kg −1 ].The solute concentration change in the surface ocean ΔC sol,oce /Δt [g kg −1 s −1 ] from changes in ice mass was computed as Here Δm sol,oce /Δt denotes the rate of change of the mass of a dissolved tracer in seawater [g s −1 ], m seawater denotes the mass of seawater [kg], ρ seawater denotes seawater density [kg m −3 ], z denotes the thickness of the surface ocean layer [m], and Δm ice /Δt denotes the gross sea-ice growth rate (positive values) and melt rate (negative values) [kg m −2 s −1 ].Hence, with a constant tracer concentration prescribed in the sea ice, growth of sea ice leads to an increase of the oceanic tracer concentration, while melt leads to a decrease of oceanic tracer concentration.Within each time step, Δt (1.2 h), the model computes the thermodynamics of sea ice (i.e., growth and melt) and the dynamics of sea ice (i.e., convergence and divergence).In the model, sea-ice salinity is set to a constant value of 5 g kg −1 , which results in a distribution of modeled water masses and large-scale circulation features that are in good agreement with observations ( Jungclaus et al., 2013), even though open ocean convection in the Southern Ocean might be too strong (Stössel et al., 2015).

Model initialization and boundary conditions
The model was initialized from the oceanic, sea-ice and biogeochemical fields obtained from the year 1979 of the fifth phase of the Coupled Model Intercomparison Project (CMIP5) using the Max-Planck-Institute for Meteorology's Earth System Model (MPI-ESM) in GR15/L40 (1.5°) resolution (Ilyina et al., 2013).
Elementa: Science of the Anthropocene • 4: 000136 • doi: 10.12952/journal.elementa.000136 We ran a spin-up simulation of 920 years by forcing the model each year with daily heat, freshwater and momentum fluxes computed from the OMIP climatological atmospheric forcing (Röske, 2006).To correct for the unbalanced globally integrated surface freshwater flux and other errors in the forcing fields of precipitation and river runoff, surface salinity is restored with a time-linear relaxation coefficient of 3.3•10 −7 s −1 to the conditions derived from the Polar Science Center Hydrographic Climatology (PHC2, which is based on the climatology provided by Levitus et al. (1998), with modifications in the Arctic based on Steele et al. (2001)).The atmospheric CO 2 concentration was kept constant at the level at the end of the year 1979 (337.4 ppm) throughout the spin-up simulation.
The spin-up simulation reaches quasi-equilibrium concerning deep-water salinity and temperature, the meridional overturning circulation, the sea-ice state, and the air-sea CO 2 flux after 700 years of simulation, and was continued until the year 920.

Experimental design
To identify the influence of climatic changes on the strength of the SICP, we simulated two idealized climate scenarios (Section 2.3.1).For both idealized climate scenarios, we performed several simulations with different prescribed TA and DIC concentrations in the sea ice, spanning the range of observed values (Section 2.3.2).All simulations were initialized at the end of the spin-up simulation, and were run for 500 years of simulation.

Idealized climate scenarios
In the equilibrium (EQ) scenario, which roughly represents today's physical forcing on the ocean and atmosphere with fixed atmospheric CO 2 concentrations, we forced the model with the OMIP atmospheric forcing as used in the spin-up simulation.
In the global warming (WARM) scenario, we forced the model with a warming climate, while the atmospheric CO 2 concentration is kept constant at 337.4 ppm as in the EQ scenario.As we used an ocean-only model, we needed to provide forcing fields for the WARM scenario.To obtain these, we first computed the climatological anomaly of each field between the last 30 years of the 1000 year long CMIP5 control experiment and the climatological mean of the years 66 to 95 of the CMIP5 experiment with an increase in atmospheric CO 2 concentration of 1% yr −1 , calculated both using MPI-ESM (Ilyina et al., 2013).By adding these anomaly fields to the climatological fields of the OMIP forcing, we generate a 'new' forcing.This 'new' forcing was applied in the WARM scenario by linearly increasing the added anomaly over the first 65 years of the simulation from zero to the full anomaly, at which level it is then kept constant until the end of the 500 year model run.The same anomaly forcing was applied to sea-surface salinity relaxation in the WARM scenario.We further note that an ocean-only model implies no feedback from ocean processes to the atmospheric climate.

SICP parameter setting
For each idealized climate scenario, we used a constant sea-ice TA:DIC of 2000:2000 mmol m −3 as our control parameter setting (Table 1).This control parameter setting is very similar to standard oceanic conditions, and thus minimizes fractionation of DIC and TA, and consequently oceanic pCO 2 changes, upon sea-ice growth and melt.The impact of the SICP can then readily be evaluated by assessing changes that occur for more realistic settings for TA and DIC relative to this control setting.The individual simulations were designed to cover the range of observed concentrations and ratios of TA:DIC in the sea ice (Table 1; Rysgaard et al., 2007Rysgaard et al., , 2009Rysgaard et al., , 2013;;Fransson et al. 2011;Geilfus et al. 2013).In addition to these observed parameter settings, we also performed simulations with a parameter setting of both TA and DIC at zero in the sea ice, as this setting is the default setup in most ocean-sea-ice biogeochemical models, including MPIOM/HAMOCC (e.g., Steiner et al., 2014, and  not change the oceanic TA:DIC ratio but simulates full drainage of brine-associated TA and DIC into the ocean upon sea-ice growth and full dilution of oceanic TA and DIC by melt water.For our experiments we follow a naming scheme that is detailed in Table 1.

Quantification of the SICP
Within each simulated, idealized climate scenario the circulation fields of the SICP experiments are identical, and the changes due to biological production are insignificant (Text S1, Figure S1).Correspondingly, in each idealized climate scenario the effect of the SICP-induced changes in ocean properties, e.g., the air-sea CO 2 flux, can be separated by examining differences between simulations with different SICP parameter settings.We quantified the changes in biogeochemical properties induced by the SICP by subtracting within each idealized climate scenario the simulation results based on the control parameter setting (CTRL) from the simulation results based on the SICP parameter settings.The differences are denoted by "Δ"; e.g., ΔEQ_R2h denotes EQ_R2h minus EQ_CTRL (where R2h refers to a TA:DIC ratio of 2 with high tracer concentrations; Table 1).

Determination of oceanic brine and melt-water distribution
SICP-induced changes in TA can be used to illustrate the footprints of brine and melt-water in the ocean, as TA is not influenced by gas-exchange processes and changes to the SICP due to biological production are not significant (Text S1, Figure S1).We analyzed SICP-induced TA changes between the CTRL and R0 experiments in both the EQ scenario and the WARM scenario.In this analysis, a positive SICP-induced change in oceanic TA indicates that the increase of TA from local brine release and/or from oceanic transport of brine is larger than the decrease of TA from local sea-ice melt and/or from the oceanic transport of seaice melt-water.In the following, we refer to such a positive SICP-induced change in oceanic TA as "brine footprint", and to a negative SICP-induced change in oceanic TA as "melt-water footprint".

Sea-ice volume
In the Arctic, the simulated sea-ice extent in the EQ scenario is in good agreement with the observational data of the National Oceanic and Atmospheric Administration (NOAA)/National Snow and Ice Data Center (NSIDC) (Meier et al., 2011).The Arctic sea-ice volume in the EQ scenario approximates the years 1963 to 2000 of the CMIP5 historical MPI-ESM runs (Notz et al., 2013; Figure 2a).In the Antarctic, the model underestimates the observed sea-ice extent, in particular during summer, resulting in a minimum simulated sea-ice volume that is close to zero in both idealized climate scenarios (Figure 2b).This underestimation is a known bias of MPIOM (Marsland et al. 2003;Jungclaus et al. 2013).

Figure 2
Simulated annual sea-ice volume and production over 500 years.
Annual maximum (solid lines) and minimum (dashed lines) seaice volume (1000 km 3 ) derived from monthly mean values of the EQ scenario (black lines) and WARM scenario (red lines) for the Arctic (a) and Antarctic region (b); annual sea-ice production (1000 km 3 ) for the Arctic (c) and Antarctic region (d) derived by subtracting the annual minimum from the annual maximum sea-ice volume.
doi: 10.12952/journal.elementa.000136.f002 Elementa: Science of the Anthropocene • 4: 000136 • doi: 10.12952/journal.elementa.000136 Compared to the EQ scenario, the WARM scenario exhibits a continuous reduction of both the summer and winter sea-ice volume during the first 65 years of simulation in both hemispheres, which subsequently equilibrate at lower volumes (red lines in Figure 2a,2b).This evolution of sea-ice volume in both hemispheres is a direct consequence of the warming applied in the WARM scenario, which after 65 years of simulation is kept constant.In the WARM scenario, the annual sea-ice production is enhanced by 8% in the Arctic and reduced by 20% in the Antarctic compared to the EQ scenario (Figure 2c,2d).The Arctic sea-ice volume of the WARM scenario stabilizes at similar levels to those obtained in MPI-ESM by the end of the century in the CMIP5 RCP4.5 scenario (Notz et al., 2013).

Annual net sea-ice growth and melt
In both idealized climate scenarios on an annual mean, sea ice melts along the sea-ice margins and grows in the ice interior in both hemispheres (Figure 3).In the WARM scenario, the position of the Arctic sea-ice margin retreats northward relative to the EQ scenario (Figure 3a, 3c).In addition, in the WARM scenario over larger areas of the central Arctic Ocean, the actual annual sea-ice production is increased but the annual net sea-ice growth is reduced compared to the EQ scenario (Figure 1c).This reduction is because in the WARM scenario the sea ice grown in the central Arctic is also melted within this basin, and less efficiently exported southward than in the EQ scenario.In contrast, the Antarctic annual mean position of the sea-ice margin and regions of annual net sea-ice growth and melt barely change between the EQ scenario and the WARM scenario (Figure 3b and 3d).

Mixed layer depth
Deep mixing occurs in the Labrador Sea and also to some extent in the Greenland Sea in the EQ scenario (Figure 4a).In the WARM scenario, the warmer climate reduces the mixed layer depth in the North Atlantic compared to the EQ scenario (Figure 4c).Similarly, the winter mixed layer depth is reduced in the central Arctic Ocean in the WARM scenario compared to the EQ scenario.In the Southern Ocean, very deep mixing occurs in the Weddell and Ross Seas in both the EQ scenario and the WARM scenario (Figure 4b and 4d).This open ocean deep mixing in the Southern Ocean is likely overestimated by the model, while dense-water formation on the shelves around Antarctica is underestimated (e.g., Marsland et al., 2004;Jungclaus et al., 2013;Stössel et al., 2015).Any uncertainty in the modeled location and strength of dense-water formation might be relevant for the strength of the SICP-induced air-sea CO 2 flux.As for the northern hemisphere, in the WARM scenario deep mixing is reduced compared to the EQ scenario, particularly in the Weddell Sea.

Meridional overturning stream function
The meridional overturning circulation is characterized by a clockwise rotating cell above 3000 m depth and north of 35°N, which is driven by the sinking of dense water in the North Atlantic Ocean (Figure 5b).Dense-water sinking in the Southern Ocean drives a counterclockwise rotating cell adjacent to Antarctica, and at depth below 3000 m (Figure 5a and 5b).The model simulates a maximum southward transport in the Atlantic Ocean of about 20 Sv (1 Sv = 10 6 m³ s −1 ) at 40°N, and 12 Sv at 60°N for the EQ scenario (Figure 6).In the WARM scenario, the strength of the southward transport continuously weakens during the first 65 years of changing forcing to a minimum of 13 Sv at 40°N and 6 Sv at 60°N (Figure 6).While at 40°N the southward transport recovers during the period of constant forcing (years 65 to 500) to the initial strength of about 20 Sv after 120 years of simulation, at 60°N the recovery of the strength of the southward transport is much slower (Figure 6), and the transport is still around 3 to 4 Sv weaker than in the EQ scenario after 500 years of simulation (Figure 5d, Figure 6).Similarly, the strength of the counterclockwise rotating cell adjacent to Antarctica, and at depth below 3000 m, is reduced in the WARM compared to the EQ scenario (Figure 5c).Any change in the meridional overturning circulation might influence the brine export into the deep ocean, and thus the strength of the SICP-induced air-sea CO 2 flux.

Distribution of brine and sea-ice melt-water in the ocean
To analyze the distribution of brine and melt-water footprints, we examined the last 50 years of simulation, when the SICP-induced change in TA is stable in both the EQ scenario and the WARM scenario (Figures S2 and S3).

EQ scenario
A prevailing brine footprint in the surface ocean (positive SICP-induced change in TA) is simulated in areas of net sea-ice growth in the interior of the ice pack, and a melt-water footprint (negative SICP-induced change in TA) is simulated in net sea-ice melt areas along and partly beyond the sea-ice margins in both hemispheres (Figure 3a, 3b; Figure 7a, 7b).Generally the strength of the simulated brine footprint increases with net sea-ice growth, and the melt-water footprint increases with net sea-ice melt, although ocean circulation can modify these regional patterns (Figure 3a, 3b; Figure 7a, 7b).One prominent regional modification is simulated in the Beaufort Gyre; here the brine footprint is larger in the center of the gyre than along the coast of the Beaufort Sea (Figure 7a), although net sea-ice growth is larger along the coast (Figure 3a).This difference is because the Beaufort Gyre causes convergence towards its center, and hence a trapping of brine in the center of the gyre.At depth, the brine footprint is traceable to intermediate depth (about 350 m) in the Arctic, locally reaching as deep as 900 m in the Beaufort Sea and 2200 m in the Baffin Bay (not shown in detail; Figure 8a).The North Atlantic Ocean to a depth of 3500 m is characterized by a prevailing melt-water footprint, indicating that brine contribution to dense-water formation and its spreading into the North Atlantic Ocean is low in the model.In the Southern Ocean, on a zonal mean the brine footprint is transported from the upper ocean (south of 60°S) into the deep Pacific and Indian Ocean (below 2500 m depth) and into the deep Atlantic Ocean (below 3500 m depth) as far north as 55°N (Figure 8a, 8b), suggesting some brine contribution to dense-water formation in the Southern Ocean.

WARM scenario
As in the EQ scenario, the regional distribution of the annual mean surface ocean brine and melt-water footprints in the WARM scenario are mainly determined by the annual net sea-ice growth and melt (Figure 3c, 3d; Figure 7c, 7d).Considerably reduced annual net sea-ice growth in the Beaufort Sea strongly reduces the brine footprint in the WARM scenario compared to the EQ scenario (Figure 7c).The retreat of the sea-ice margin towards the north and reduced sea-ice import and melt in the Greenland and Labrador Seas (Figure 3c) in the WARM scenario (Figure 3c) lead to a considerably lower melt-water footprint compared to the EQ scenario in these seas (Figure 7c).In the Antarctic, both the brine and melt-water footprints in the WARM scenario are reduced compared to the EQ scenario (Figure 7d), resulting from the reduced annual sea-ice production (Figure 2d).
On a zonal mean, at depth the brine footprint in the Arctic Ocean is reduced and reaches to shallower depths (about 200 m) in the WARM compared to the EQ scenario (Figure 8c).This reduction is caused by reduced net sea-ice growth in the central Arctic Ocean (Figure 3c), and downwelling of the melt-water footprint in the center of the Beaufort Gyre (Figure 7c).The water spreading southward in the North Atlantic Ocean (south of 63°N and to a depth of 3000 m) is characterized by a melt-water footprint similar to the EQ scenario (Figure 8c).This similarity indicates that the reduced sea-ice melt in the Greenland and Labrador Seas (Figure 3c) is in the WARM scenario roughly compensated by the reduced southward transport of water and brine (Figure 6) in the WARM scenario.In the Antarctic, the reduced annual production of sea ice (Figure 2d), and hence of brine, and the reduced strength of the global meridional overturning circulation adjacent to Antarctica and below 3000 m depth (Figure 5c) result in a reduced brine footprint in the deep Southern, Pacific, and Indian Oceans in the WARM compared to the EQ scenario (Figure 8d).

The SICP-induced air-sea CO 2 flux
We define a positive sign for a net CO 2 flux from the atmosphere to the ocean.SICP-induced air-sea CO 2 flux denotes the difference between the SICP experiments and the CTRL experiments (e.g., EQ_R0 minus EQ_CTRL), with a positive sign denoting that the SICP experiment shows more CO 2 uptake or less CO 2 out-gassing than its corresponding CTRL experiment, and vice versa for a negative sign.As one would expect from the temporally rather stable SICP-induced changes in surface ocean TA in both idealized climate scenarios (Figures S2 and S3), the simulated SICP-induced air-sea CO 2 fluxes are quasi-stationary after 150 years of simulation in the Arctic (Figures S4c, S4d and S5c, S5d), and after 400 to 450 years in the Antarctic (Figures S4b and S5b).The SICP-induced air-sea CO 2 flux is quantified by averaging the results over the simulation years 451 to 500.

The EQ climate scenario
The regional distribution of the SICP-induced air-sea CO 2 exchange in both hemispheres is broadly determined by the regional distribution of surface brine and melt-water footprints on both annual (Figures 7a, 7b and 9c-9f,) and seasonal scales (Figures S6a-S6d and S7e-S7l).Over most areas of the Arctic and Antarctic, an increased brine footprint in the surface ocean leads to a SICP-induced CO 2 out-gassing, whereas an increased melt-water footprint results in a SICP-induced oceanic CO 2 uptake.A large TA:DIC ratio of 2 prescribed in the sea ice (EQ_R2h versus EQ_R0) strengthens the SICP-induced air-sea CO 2 exchange, with the regional flux patterns remaining the same (Figures 9c-9f and S7e-S7l).
The strength of the SICP-induced air-sea CO 2 exchange, however, is not always just linearly related to the brine or melt-water footprint at the ocean surface.For example, in the Beaufort Gyre the annual mean SICP-induced oceanic CO 2 out-gassing is relatively small (Figure 9c, 9e), but the annual mean brine content is large (Figure 7a).This difference is because in the Beaufort Sea the almost complete winter sea-ice cover inhibits air-sea CO 2 exchange (Figure S7c, S7g, S7k), while the reduced summer sea-ice cover combined with a high surface brine footprint (Figure S6a) leads to an enhanced SICP-induced CO 2 out-gassing during summer (Figure S7e, S7i).This difference demonstrates that the annual SICP-induced CO 2 out-gassing in the Beaufort Sea (Figure 9c, 9e) is dominated by the summer conditions with reduced brine footprint and reduced sea-ice coverage compared to winter (Figure S6a, S6c).
In both the Arctic and Antarctic sea-ice zones, the SICP-induced oceanic CO 2 uptake increases in the experimental order of R1.1h < R1.1l < R0 < R2l < R2h (Table 2).This ordering indicates that during seaice melt a maximum dilution of the surface TA and DIC concentrations without changing the oceanic ratio (R0) increases the SICP-induced oceanic CO 2 uptake more strongly than the melt-water release of some concentration of TA and DIC with a ratio of 1.1 (R1.1h,R1.1l).Because the oceanic ratio is already 1.1, the R1.1h and R1.1l experiments do not change the oceanic TA:DIC ratio.In contrast, the dilution effect on the SICP-induced oceanic CO 2 uptake is overwhelmed by the melt-water release of a TA:DIC ratio of 2 (R0 < R2).For this sea-ice ratio of 2, the increase in the oceanic TA:DIC ratio is larger for higher concentrations of TA and DIC released upon melt (R2l < R2h).Independent of the concentrations and ratios of TA:DIC prescribed in the sea ice, the SICP-induced oceanic CO 2 uptake in the Antarctic is larger than in the Arctic (Table 2).This difference is a consequence of the more efficient brine export away from the Antarctic mixed layer (Figure 8a, 8b; Section 3.7).

The WARM climate scenario
Similar to the EQ scenario, the magnitude of the SICP-induced air-sea CO 2 exchange in the WARM scenario is determined by the strength of the surface brine and melt-water footprints, with the regional patterns being enhanced for a larger concentration and ratio of TA and DIC prescribed in the sea ice (WARM_R0 versus WARM_R2h) (Figures 7c, 7d and 10c-10f ).The larger summer melt-water footprint in the central Arctic Ocean compared to the EQ scenario (Figure S6e) causes an enhanced SICP-induced oceanic CO 2 uptake in the WARM scenario in summer (Figure S8e, 8i) and annually (Figure 10c, 10e; Table 2).The strength of the SICP-induced enhanced oceanic CO 2 uptake is reduced in the Greenland and Labrador Sea in the WARM compared to the EQ scenario (Figure 10c, 10e).This loss in strength is caused by the reduced sea-ice import and melt (Figure 3c), which reduces the melt-water footprint (Figure 7c).Similar to the EQ scenario, an Annual mean air-sea CO 2 flux (g C m -2 yr -1 ) in the EQ_CTRL experiment (a,b), annual mean SICP-induced air-sea CO 2 flux (g C m -2 yr -1 ) in the EQ_R0 experiment (c,d), and annual mean SICP-induced air-sea CO 2 flux in the EQ_R2h experiment (e,f ) for the Arctic (a,c,e) and Antarctic (b,d,f ) regions.The airsea CO 2 flux in the EQ_CTRL experiment is positive for a net CO 2 flux from the atmosphere to the ocean (oceanic CO 2 uptake), and negative for the other direction (oceanic CO 2 out-gassing).SICPinduced air-sea CO 2 flux denotes the difference between the SICP experiment (EQ_R0, EQ_R2h) and the EQ_CTRL experiment, with a positive sign indicating that the SICP experiment shows more CO 2 uptake or less CO 2 out-gassing than the CTRL experiment; negative values indicate the opposite.Black lines show the annual mean sea-ice margins.Latitudinal bands are shown for 80°, 60° and 40°.Results are averaged over the simulation years 451 to 500.Note that the color scales are nonlinear.Regional extent is defined as the annual mean sea-ice extent averaged over the years 451 to 500 (Figure 3).doi: 10.12952/journal.elementa.000136.t002 enhanced SICP-induced oceanic out-gassing, caused by intensive mixing (Figure 4c) and brine footprint entrainment to the surface, is simulated in the North Atlantic Ocean south of the average sea-ice edge at 45°W-5°W/40°N-61°N (Figure 10c, 10e).The strength of this North Atlantic out-gassing is weaker in the WARM scenario than in the EQ scenario, which is related to the reduced North Atlantic southward water (and brine footprint) transport (Figure 6).
In contrast to the EQ scenario, in the WARM scenario the relative order of the Arctic SICP-induced oceanic CO 2 uptake rate is R0 < R1.1h < R2h (Table 2).This ordering shows that in the WARM scenario a low ratio of TA:DIC released upon sea-ice melt (R1.1h) increases the SICP-induced oceanic CO 2 uptake more strongly than a maximum surface ocean TA and DIC dilution upon sea-ice melt without changing the oceanic TA:DIC ratio (R0).This outcome is because in the WARM scenario, the Arctic surface ocean TA:DIC ratio is considerably lower than 1.1, which is mainly caused by the summer ice-free central Arctic Ocean allowing for an enhanced net oceanic CO 2 uptake compared to the EQ scenario (Table 2; Figure S9c).Consequently, the release of a low TA:DIC ratio of 1.1 upon sea-ice melt increases the oceanic TA:DIC ratio and induces an enhanced CO 2 uptake in the WARM scenario compared to the EQ scenario.In the Southern Ocean, TA:DIC ratios are around 1.1 both in the WARM scenario and in the EQ scenario.Hence, both scenarios have a similar order of the SICP-induced oceanic CO 2 uptake strength (R1.1h < R0 < R2h) (Table 2).As in the EQ scenarios, in the WARM scenario the SICP-induced oceanic CO 2 uptake is larger in the Antarctic than in the Arctic (Table 2), which is related to a more efficient brine export from the Antarctic than from the Arctic (Figure 8c, 8d; Section 3.7).

Estimation of the fraction of brine drainage below the polar mixed layer
In a conceptual view of the SICP, one would expect that in full equilibrium the additional CO 2 uptake caused by the SICP is out-gassed in non-polar regions in order to sustain the global flux balance.Thus, in equilibrium the mass of additional SICP-induced non-polar out-gassed CO 2 equals the mass of brine-associated DIC leaving the polar oceanic mixed layer.In our simulations, the global SICP-induced oceanic CO 2 flux balance (zero flux) is not established yet (Table 2).Nevertheless, because our simulations show a stable SICP-induced air-sea CO 2 exchange and brine or melt-water footprints in the sea-ice zones of both hemispheres (Figures S2,S3,S4b,S4d,and S5b,S5d), we can use this conceptual understanding to estimate the fraction of DIC that is transported away from the polar mixed layer.
In the EQ scenario, the model simulates an annual sea-ice production of 14.1 • 10 12 m 3 in the Arctic (Figure 2c).Based on this annual sea-ice production, we infer that the EQ_R0 experiment releases 300 Tg C yr −1 more brine-associated DIC to the surface ocean in the Arctic than the EQ_CTRL experiment ( 300  Annual mean air-sea CO 2 flux (g C m -2 yr -1 ) in the WARM_ CTRL experiment (a,b), annual mean SICP-induced air-sea CO 2 flux (g C m -2 yr -1 ) in the WARM_ R0 experiment (c,d), and annual mean SICP-induced air-sea CO 2 flux in the WARM_R2h experiment (e,f ) for the Arctic (a,c,e) and Antarctic (b,d,f ) regions.The air-sea CO 2 flux in the WARM_CTRL experiment is positive for a net CO 2 flux from the atmosphere to the ocean (oceanic CO 2 uptake), and negative for the other direction (oceanic CO 2 out-gassing).SICPinduced air-sea CO 2 flux denotes the difference between the SICP experiment (WARM_R0, WARM_ R2h) and the WARM_CTRL experiment, with a positive sign denoting that the SICP experiment shows more CO 2 uptake or less CO 2 out-gassing than the CTRL experiment; negative values indicate the opposite.Black lines show the annual mean sea-ice margins.Latitudinal bands are shown for 80°, 60° and 40°.Results are averaged over the simulation years 451 to 500.Note that the color scales are nonlinear.doi: 10.12952/journal.elementa.000136.f010 Table 2) equals the export flux of brine-associated DIC, we infer that 0.6% of the 300 Tg C of the brineassociated DIC produced every year leaves the Arctic mixed layer.In the Antarctic, about 3.1% of the annually produced brine-associated DIC is exported from the mixed layer.Performing this calculation for the SICP-induced air-sea CO 2 flux in the EQ_R2h experiment, we find a slight increase in the fraction of brine-associated DIC exported from the mixed layer, with values of 2.0% in the Arctic and 5.9% in the Antarctic.All other SICP experiments in the EQ scenario range between the two endmembers.
In the WARM climate scenario, we infer a brine-associated DIC export of −0.2% to 0.6% in the Arctic, with the lower estimate indicating net brine accumulation in the Arctic mixed layer.In the Antarctic, the brine-associated DIC export ranges between 1.9 and 4.7%.

Discussion
In this study we separate the mechanisms that influence the SICP-induced air-sea CO 2 flux.To do so, we have used a three-dimensional ocean-sea-ice-biogeochemical model and run a set of idealized SICP experiments combined with idealized climate scenarios.The SICP scenarios cover the range of observed TA and DIC compositions in the sea ice.The idealized climate scenarios (both forced with constant atmospheric CO 2 concentration) are designed to represent (i) the present-day ocean and sea-ice state (EQ scenario); and (ii) a future reduced oceanic overturning and Arctic sea-ice state (WARM scenario) that is comparable to the RCP4.5 scenario.Rysgaard et al. (2007) emphasize that the strength of the SICP-induced air-sea CO 2 flux is particularly sensitive to the fraction of brine being transported away from the polar mixed layer.Under the present-day strength of oceanic overturning (EQ scenario), the model-derived fractional export of brine-associated DIC ranges between 0.6 and 2% in the Arctic, and between 3.1 and 5.9% in the Antarctic.A reduction in the strength of the oceanic overturning circulation (WARM scenario) decreases the brine export to a maximum of 0.6% in the Arctic and 4.7% in the Antarctic.These differences in the fractions of brine-associated DIC export between the Arctic and Antarctic explain the larger SICP-induced oceanic CO 2 uptake in the Antarctic compared to the Arctic in both idealized climate scenarios.An earlier study using MPIOM (Stössel et al., 2002) has shown that sinking of dense water in the North Atlantic is rather insensitive to brine density.This insensitivity is confirmed by our model results showing that a relatively large brine footprint remains above the Arctic cold halocline, whereas the water spreading into the intermediate to deep North Atlantic contains more sea-ice melt-water than brine (Figure 8a, 8c).In contrast, dense-water sinking to depth in the Southern Ocean (in MPIOM) is to a larger extent driven by surface density changes induced by sea-ice brine drainage (Stössel et al., 2002).This effect is shown by the increased brine footprint in waters spreading from the Southern Ocean along the seafloor into the global ocean (Figure 8).In general, dense-water formation in global ocean general circulation models depends critically on the representation of mixing and of open water in the sea-ice cover (leads).Over leads, intense cooling and new-ice formation induce brine rejection that potentially contributes to the transport of dense water to depth.Hence, a model with higher resolution that more accurately resolves atmosphere-ocean interactions in leads and/or uses a different mixing parameterization might result in a stronger contribution of highly saline brine to dense-water formation, and hence a stronger SICP-induced oceanic CO 2 uptake.On the regional scale, the SICP-induced air-sea CO 2 flux is largely determined by the cycle of sea-ice growth and melt, the oceanic circulation redistributing brine and sea-ice melt-water throughout the ocean, and the fraction of sea-ice coverage blocking the air-sea CO 2 exchange.Areas that are characterized by a distinct surface ocean brine footprint cause a SICP-induced oceanic CO 2 out-gassing.In contrast, surface ocean areas that are characterized by a melt-water footprint show an enhanced SICP-induced oceanic CO 2 uptake.

Mechanisms determining the strength of the SICP-induced air-sea CO 2 flux
Both the concentrations of TA and DIC and the TA:DIC ratio prescribed in the sea ice relative to the basic oceanic TA and DIC setting further determine the strength of the SICP-induced CO 2 uptake, which is higher for high ratios of TA:DIC prescribed in the sea ice (R2l, R2h) than for low ratios (R1.1l,R1.1h).This general grouping of the experiments highlights the need for refined observations of bulk sea-ice TA and DIC concentrations both on the regional and seasonal scale.Such refined observations will then also allow for a more realistic representation of the SICP in global and regional models beyond the commonly used implementation of setting both TA and DIC to zero in the simulated sea ice.
Our results suggest that a retreat of the summer Arctic sea ice, as simulated in the WARM scenario, increases the SICP-induced oceanic CO 2 uptake over large areas of the central Arctic Ocean.Because the sensitivity towards changes in the TA and DIC concentrations and the TA:DIC ratio is increased under higher oceanic pCO 2 conditions (e.g., Sarmiento and Gruber, 2006), we expect that this SICP-induced oceanic CO 2 uptake in the central Arctic will be enhanced under the current rising atmospheric CO 2 concentrations.

Comparison of the simulated SICP-induced air-sea CO 2 flux to related studies
In a specific sea-ice melt region of the Greenland Sea (74°N-80°N, 15°W-0°E), the SICP-induced oceanic CO 2 uptake in summer (1 May to 1 September) reaches values of up to 12 g C m −2 yr −1 in the EQ_R2h scenario and 9 g C m −2 yr −1 in the WARM_R2h scenario.These fluxes are lower than the mean summer oceanic CO 2 uptake of 18 g C m −2 yr −1 of the SICP estimated from a simple model of melted sea ice in the mixed layer for the same season and region and a comparable SICP scenario (Rysgaard at al., 2009).This SICP-induced CO 2 flux in the EQ scenario constitutes around 30% of the present-day mean uptake in spring (45 ± 21 g C m −2 yr −1 ) and summer (37 ± 14 g C m −2 yr −1 ) in the Greenland Sea (Nakaoka et al., 2006).The simulated annual maximum for modeled SICP-induced oceanic CO 2 uptake in this region (74°N-80°N, 15°W-0°E) ranges between 13 and 28% of the observed uptake in the EQ scenario, and between 7 and 16% of the observed uptake in the WARM scenario.The observed uptake values used for these estimates are based on measurements of the annual present-day CO 2 flux from Nakaoka et al. (2006; mean of 52 ± 20 g C m −2 yr −1 ) and Yasunaka et al. (2016; mean of 48 ± 13 g C m −2 yr −1 ) derived from different regions in the Greenland Sea.Similarly, in other net sea-ice melt or growth regions, we find a considerable influence of the SICP-induced air-sea CO 2 flux on both the seasonal and annual scale.
Although the SICP considerably influences the high-latitude air-sea CO 2 flux, when averaged over the polar regions the SICP-induced CO 2 uptake only drives a small fraction of the net oceanic CO 2 uptake.This fraction is up to 4% for the Arctic and 12% for the Antarctic in the EQ scenario, dropping to 1% for the Arctic and 5% for the Antarctic sea-ice zone in the WARM scenario.In absolute numbers, the SICP drives a maximum polar oceanic CO 2 uptake of 14 Tg C yr −1 in the EQ scenario, and 8 Tg C yr −1 in the WARM scenario (Table 2), which is up to 7% of the simulated net polar oceanic CO 2 uptake.These values are lower than the derived estimates of 33 to 83 Tg C yr −1 assuming full brine export (Rysgaard et al., 2011), because considerable amounts of brine remain in the mixed layer in our experiments.The rather small SICP-induced variations between the climatic scenarios further demonstrate that changes in the physical ocean-sea-ice state barely influence the SICP-induced air-sea CO 2 flux when averaged over the polar regions.However, as discussed above, climatic changes induce considerable variations in the strength of the SICP on regional scales.For example, in the future summer sea-ice free central Arctic Ocean, the annual mean SICP-induced CO 2 uptake is likely to increase significantly.
An accompanying study simulating the SICP for pre-industrial climate conditions using a different global-ocean-sea-ice-biogeochemical model (NEMO-LIM-PISCES, Moreau et al., 2016) shows a similarly small fraction of brine transport away from the polar mixed layer, and similar regional patterns of the SICPinduced air-sea CO 2 flux that are also mainly determined by net sea-ice growth and melt areas.Compared to the current global oceanic net CO 2 uptake (mean of 2600 ± 500 Tg C yr −1 ) which is forced by the rising atmospheric CO 2 concentrations (Le Quéré et al., 2015), the SICP-induced CO 2 uptake is small in the simulated idealized scenarios.However, for long term climatic changes, such as the transition between the last interglacial to glacial, which was accompanied by an increase in the sea-ice extent (Sarnthein et al., 2003;de Vernal et al., 2005;Gersonde et al., 2005), the SICP might have contributed to a significant drawdown of the atmospheric CO 2 concentrations.Indeed, using a model of intermediate complexity, Bouttes et al. (2010) show that a direct transport of brine to depth can explain a large part of the interglacial to glacial oceanic CO 2 uptake.In this case, the brine-associated DIC is transported directly from the surface to the deep ocean and the dense brine transported to depth stratifies the ocean, which are crucial for the long term trapping and accumulation of DIC in the deep ocean.To further constrain the influence of the SICP on the interglacial-glacial atmospheric CO 2 drawdown, simulations with an ocean-sea-ice-biogeochemical general circulation model resolving both regional circulation features and sea-ice dynamical processes are required.

Summary
Using a three-dimensional ocean-sea-ice-biogeochemical model for a set of SICP experiments within two idealized climate scenarios, we have separated the mechanisms that determine the strength of the SICPinduced air-sea CO 2 flux.The efficiency of the SICP itself is determined by the concentrations of TA and DIC and the TA:DIC ratio prescribed in the sea ice relative to the polar oceanic TA and DIC conditions; the efficiency of the SICP-induced oceanic CO 2 uptake is determined by the export of brine-associated DIC from the polar mixed layer to depth.Because that export to depth is larger in the Antarctic than in the Arctic, the SICP-induced oceanic CO 2 uptake is higher in the Antarctic.Particular uncertainties in the quantification of the SICP hence arise from modeling issues related to the sites and strength of densewater formation and its transport throughout the ocean.Within the sea-ice zone, the strength and sign of the SICP-induced air-sea CO 2 exchange is controlled mainly by the relative distribution of brine produced during sea-ice formation (SICP-induced oceanic CO 2 out-gassing), and melt-water released during sea-ice melt (SICP-induced oceanic CO 2 uptake).These general air-sea CO 2 exchange features are modified further by the influence of sea-ice coverage on the air-sea CO 2 exchange.Surface layer is defined as the upper 12 m of the ocean.Results are annually averaged.SICP-induced change in TA (mmol m −3 ) denotes the difference in TA between the SICP experiment (WARM_R0, WARM_R2h) and the WARM_CTRL experiment, with a positive sign denoting that the SICP causes an increase in the oceanic TA concentration compared to the WARM_CTRL experiment (brine footprint), and a negative sign denoting that the SICP causes a decrease in the oceanic TA concentration compared to the WARM_CTRL experiment (melt-water footprint).

•
Figure S4.Simulated SICP-induced air-sea CO 2 flux in the EQ scenario.
Annual mean SICP-induced air-sea CO 2 flux (Tg C yr −1 ) averaged over the global ocean (a), the Antarctic sea-ice zone (b), the central Arctic Ocean (c), and the Arctic sea-ice zone (d).SICP-induced air-sea CO 2 flux denotes the difference between the SICP experiment (EQ_R0, EQ_R2h) and the EQ_CTRL experiment, with a positive sign denoting that the SICP experiment shows more CO 2 uptake or less CO 2 out-gassing than the CTRL experiment.
Regional extent of the central Arctic Ocean as shown in Figure S10, regional extent of the Arctic and Antarctic sea-ice zone as defined by the annual mean sea-ice extent (Figure 2).Seasonal mean air-sea CO 2 flux (g C m -2 yr −1 ) in the EQ_CTRL experiment (a-d), SICP-induced air-sea CO 2 flux (g C m −2 yr −1 ) in the EQ_R0 experiment (e-h), and SICP-induced air-sea CO 2 flux in the EQ_R2h experiment (i-l).The air-sea CO 2 flux in the EQ_CTRL experiment is positive for a net CO 2 flux from the atmosphere to the ocean (oceanic CO 2 uptake), and negative for the other direction (oceanic CO 2 out-gassing).SICP-induced air-sea CO 2 flux denotes the difference between the SICP experiment (EQ_R0, EQ_R2h) and the EQ_CTRL experiment, with a positive sign denoting that the SICP experiment shows more CO 2 uptake or less CO 2 out-gassing than the CTRL experiment.Seasonal mean air-sea CO 2 flux (g C m -2 yr −1 ) in the WARM_CTRL experiment (a-d), SICP-induced air-sea CO 2 flux (g C m −2 yr −1 ) in the WARM_R0 experiment (e-h), and SICP-induced air-sea CO 2 flux in the WARM_R2h experiment (i-l).The air-sea CO 2 flux in the WARM_CTRL experiment is positive for a net CO 2 flux from the atmosphere to the ocean (oceanic CO 2 uptake), and negative for the other direction (oceanic CO 2 out-gassing).SICP-induced air-sea CO 2 flux denotes the difference between the SICP experiment (WARM_R0, WARM_R2h) and the WARM_CTRL experiment, with a positive sign denoting that the SICP experiment shows more CO 2 uptake or less CO 2 out-gassing than the CTRL experiment.

Figure 3
Figure 3 Simulated annual net sea-ice growth and melt.Annual net sea-ice growth (positive values, red) and melt (negative values, blue) (m) in the EQ climate scenario (a,b), and the WARM climate scenario (c,d).The black lines indicate the annual mean position of the sea-ice margins.Latitudinal bands are shown for 80°, 60° and 40°.Results are averaged over the simulation years 451 to 500.Note that the color scale is nonlinear.doi: 10.12952/journal.elementa.000136.f003

Figure 4
Figure 4 Simulated mixed layer depth.Results for the mean mixed layer depth (m) in the EQ climate scenario (a,b) and the WARM climate scenario (c,d) for the Arctic region in March (a,c) and the Antarctic region in September (b,d).The mixed layer depth is diagnosed in the model as the depth where the potential density is 0.125 kg m -3 larger than at the surface.Latitudinal bands are shown for 80°, 60° and 40°.Results are averaged over the simulation years 451 to 500.Note that the color scale is nonlinear.doi: 10.12952/journal.elementa.000136.f004 Figure 5 Simulated meridional overturning stream function.Annual mean meridional overturning stream function (Sv) in the EQ climate scenario (a,b), where positive values indicate clockwise circulation, and change in the meridional overturning stream function in the WARM climate scenario relative to the EQ scenario (c,d).The panels show results for the global ocean (a,c) and the Atlantic Ocean (b,d).All results are averaged over the simulation years 451 to 500.doi: 10.12952/journal.elementa.000136.f005 Figure 7Simulated SICP-induced change in surface oceanic TA.Annual mean SICP-induced change in TA (mmol m -3 ) in the EQ_R0 experiment (a,b), and the WARM_R0 experiment (c,d).SICP-induced change in TA denotes the difference in TA between the SICP experiment (EQ_R0, WARM_R0) and the corresponding CTRL experiment (EQ_CTRL, WARM_CTRL), with a positive sign denoting that the SICP causes an increase in the oceanic TA concentration compared to the corresponding CTRL experiment (brine footprint), and a negative sign denoting a decrease in the oceanic TA concentration compared to the corresponding CTRL experiment (melt-water footprint).Surface layer is defined as the upper 12 m of the ocean.Black lines show the annual mean sea-ice margins.Latitudinal bands are shown for 80°, 60° and 40°.Results are averaged over the simulation years 451 to 500.Note that the color scale is nonlinear.doi: 10.12952/journal.elementa.000136.f007 Figure 8 Simulated SICP-induced change in TA at depth.Zonally averaged annual mean SICP-induced change in TA (mmol m -3 ) in the EQ_R0 experiment (a,b), and the WARM_ R0 experiment (c,d).SICPinduced change in TA denotes the difference in TA between the SICP experiment (EQ_R0, WARM_ R0) and the corresponding CTRL experiment (EQ_CTRL, WARM_CTRL), with a positive sign denoting that the SICP causes an increase in the oceanic TA concentration (brine footprint), and a negative sign denoting a decrease in the CTRL oceanic TA concentration (melt-water footprint).Results are averaged over the Atlantic and Arctic Oceans (a,c), and zonally averaged over the Pacific, Indian and Southern Oceans (b,d), over the simulation years 451 to 500.Note that the color scale is nonlinear.doi: 10.12952/journal.elementa.000136.f008 Figure 9 Simulated air-sea CO 2 flux and SICP-induced air-sea CO 2 flux in the EQ climate scenario.
The dotted black horizontal line indicates the zero-line for SICP-induced air-sea CO 2 flux.Note the different y-axis scales.• Figure S5.Simulated SICP-induced air-sea CO 2 flux in the WARM scenario.Annual mean SICP-induced air-sea CO 2 flux (Tg C yr −1 ) averaged over the global ocean (a), the Antarctic sea-ice zone (b), the central Arctic Ocean (c), and Arctic sea-ice zone (d).SICP-induced air-sea CO 2 flux denotes the difference between the SICP experiment (WARM_R0, WARM_R2h) and the WARM_CTRL experiment, with a positive sign denoting that the SICP experiment shows more CO 2 uptake or less CO 2 out-gassing than the CTRL experiment.Regional extent of the central Arctic Ocean as shown in Figure S10, regional extent of the Arctic and Antarctic sea-ice zone as defined by the annual mean sea-ice extent (Figure 2).The dotted black horizontal line indicates the zero-line for SICP-induced air-sea CO 2 flux.Note the different y-axis scales.• Figure S6.Simulated SICP-induced change in surface oceanic TA.Seasonal mean SICP-induced change in TA (mmol m −3 ) in the EQ_R0 experiment (a-d), and the WARM_R0 experiment (e-h).SICP-induced change in TA implies the difference in TA between the SICP experiment (EQ_ R0, WARM_R0) and the corresponding CTRL experiment (EQ_CTRL, WARM_CTRL), with a positive sign denoting that the SICP causes an increase in the oceanic TA concentration compared to its corresponding CTRL experiment (brine footprint), and a negative sign denoting a decrease in the oceanic TA concentration compared to its corresponding CTRL experiment (melt-water footprint).Results are averaged over the summer months MJJAS (a,e) and NDJFM (b,f ), and averaged over the winter months NDJFM (c,g) and MJJAS (d,h).Superimposed is the sea-ice concentration with cross-shaded areas showing > 0.95, squared areas > 0.5, and dotted areas > 0. Surface layer is defined as the upper 12 m of the ocean.Latitudinal bands are shown for 80°, 60° and 40°.Results are averaged over the simulation years 451 to 500.Note that the color scales are nonlinear.• Figure S7.Simulated air-sea CO 2 flux and SICP-induced air-sea CO 2 flux in the EQ scenario.
First two rows show mean summer air-sea CO 2 flux and SICP-induced air-sea CO 2 flux (MJJAS (a,e,i), NDJFM (b,f,j)), and last two rows mean winter air-sea CO 2 flux and SICP-induced air-sea CO 2 flux (NDJFM (c,g,k), MJJAS (d,h,l)).Superimposed is the sea-ice fraction with cross-shaded areas showing > 0.95, squared areas > 0.5, and dotted areas > 0. Latitudinal bands are shown for 80°, 60° and 40°.Results are averaged over the simulation years 451 to 500.Note that the color scales are nonlinear.• Figure S8.Simulated air-sea CO 2 flux and SICP-induced air-sea CO 2 flux in the WARM scenario.
First two rows show mean summer air-sea CO 2 flux and SICP-induced air-sea CO 2 flux (MJJAS (a,e,i), NDJFM (b,f,j)), and last two rows mean winter air-sea CO 2 flux and SICP-induced air-sea CO 2 flux (NDJFM (c,g,k), MJJAS (d,h,l)).Superimposed is the sea-ice fraction with cross-shaded areas showing > 0.95, squared areas > 0.5, and dotted areas > 0. Latitudinal bands are shown for 80°, 60° and 40°.Results are averaged over the simulation years 451 to 500.Note that the color scales are nonlinear.• Figure S9.Simulated air-sea CO 2 flux.Time series of the annual mean air-sea CO 2 flux (Tg C yr −1 ) averaged over the global ocean (a), the area south of 40°S (b), the central Arctic Ocean (c), and the area north of 40°N (d).The air-sea CO 2 flux is positive for a net CO 2 flux from the atmosphere to the ocean (oceanic CO 2 uptake), and negative for the other direction (oceanic CO 2 out-gassing).Regional extent of the central Arctic Ocean as shown in Figure S10.Note the different y-axis scales.• Figure S10.Extent of the central Arctic Ocean.Latitudinal bands are shown for 80°, 60° and 40°.