Gsi6.0 Sea Ice Configuration Printer-friendly Version Interactive Discussion Development of Global Sea Ice 6.0 Cice Configuration for the Met Office Global Coupled Model Gmdd Gsi6.0 Sea Ice Configuration Printer-friendly Version Interactive Discussion

The new sea ice configuration GSI6.0, used in the Met Office global coupled configuration GC2.0, is described and the sea ice extent, thickness and volume are compared with the previous configuration and with observationally-based datasets. In the Arctic, the sea ice is thicker in all seasons than in the previous configuration, and there is now 5 better agreement of the modelled concentration and extent with the HadISST dataset. In the Antarctic, a warm bias in the ocean model has been exacerbated at the higher resolution of GC2.0, leading to a large reduction in ice extent and volume; further work is required to rectify this in future configurations.


Introduction
Within the Met Office's model development framework, there are four model components: atmosphere, using the Met Office Unified Model (MetUM, see Cullen and Davies, 1991;Davies et al., 2005); land surface, using the Joint UK Land Environment Simulator (JULES, see Best et al., 2011); ocean, using the Nucleus for European Modelling of the Ocean (NEMO, see Madec, 2008); and sea ice, using the Los Alamos Sea Ice Model, CICE (Hunke and Lipscomb, 2010).The UM and JULES run together as one executable, as do NEMO and CICE.UM-JULES and NEMO-CICE communicate via the OASIS coupler (Valcke, 2006).
The Met Office configurations of each component are known as Global Atmosphere (GA), Global Land (GL), Global Ocean (GO) and Global Sea Ice (GSI), and the combined system is known as the Global Coupled (GC) configuration.These terms are suffixed by a version number (e.g., "GA6.0","GC2.0").The second coupled configuration, GC2.0 (Williams et al., 2014), includes GA6.0 and GL6.0 (both described by Walters et al., 2014), GO5.0 (Megann et al., 2014) and GSI6.0.GC2.0 will be used on a range of spatial scales (regional and global), and on a range of temporal scales, from ocean forecasting (FOAM; see Blockley et al., 2014) prediction (GloSea4; see MacLachlan et al., 2014), to centennial-scale climate projections (HadGEM3; see Hewitt et al., 2011).In the present paper, we consider only the climate configuration, HadGEM3.Sea ice is a key component of the earth system because of its role in the energy balance of the polar regions.An accurate simulation of sea ice is therefore essential in fully-coupled atmosphere-ocean-ice models run on any timescale.Here, we describe the model setup and parameterisations used in GSI6.0 as part of GC2.0, and discuss how the change from the previous configuration (GSI4.0) to GSI6.0 has affected simulated sea ice extent, thickness and volume.
2 Description of GSI6.0 Thorndike et al. (1975) defined the ice thickness distribution (ITD), g, as a dimensionless function such that g(h)dh is the fraction of ice in thickness range h to h + dh; g is described by where ∇ • (v g) is the rate of change of g due to dynamical processes (v is the ice velocity), f is the rate of change of ice thickness due to thermodynamic growth and melt, and ψ gives the contribution from mechanical redistribution (ridging).A full explanation is given by Thorndike et al. (1975).

Horizontal, temporal and vertical discretisation
The GSI configurations discussed here use code revision 430 of CICE version 4.1, which allows a tripolar grid to be employed.These configurations use essentially the same family of ORCA grids as the NEMO model (see Appendix C of Hewitt et al., 2011), although CICE uses an Arakawa B grid rather than a C grid and so the CICE velocity grid points are not coincident with the NEMO velocity points.The grid and land-mask definitions required by CICE are read in directly from a file, as are the initial conditions.The sub-grid-scale ITD is modelled by dividing the ice pack at each grid point into a number of thickness categories.GSI uses five categories, plus an openwater category, which has been shown to be sufficient for climate modelling (Bitz et al., 2001).The lower bounds for the five thickness categories are 0, 0.6, 1.4, 2.4 and 3.6 m.GSI uses the zero-layer thermodynamic model of Semtner (1976) to calculate the growth and melt of the sea ice, with one layer of snow and one layer of ice in the vertical.This is not the standard scheme implemented in CICE, which has a multilayer ice model (Bitz and Lipscomb, 1999).It was not possible to use the CICE multilayer thermodynamics in GSI because the surface temperature at sea ice points, and the conductive heat flux into the ice, are currently calculated by the JULES land-surface model (which also models surface exchange over the ocean and sea ice).This would not be consistent with the CICE multilayer thermodynamics scheme, which calculates these quantities itself, so for GSI CICE has been adapted to use the zero-layer surface fluxes received from the UM atmosphere.

Thermodynamics
The sea ice albedo in GSI is calculated in the JULES land-surface model, and is a function of temperature and snow cover, including a parametrisation to represent the impact of melt ponds, and -in the zero-layer model -a parameterisation to account for the ef- 2011).The total albedo is calculated from the ice albedo α i and the snow albedo α s , weighted by the fraction of the gridbox area that is covered by snow.Bare ice albedo α b is set as a single value.The ice albedo α i is then calculated by applying corrections to α b to account for the presence of melt ponds, and for scattering within the ice pack.Melt ponds are assumed to form on bare ice when the ice temperature reaches a threshold temperature T p .As the temperature increases between T p and the melting temperature T m , melt ponds are assumed to reduce the ice albedo α i linearly, where T m is fixed at 0 • C for all simulations while the values of T p and dα i dT can be set as parameters for each simulation.
Because the ice model configuration uses a zero-layer approximation, an additional parametrization is required to account for the effects of internal scattering (e.g. from brine pockets) on the albedo.Following the suggestion of Semtner (1987), a correction ∆α i is applied to the ice albedo, where f is the fraction of incident radiation which penetrates the ice pack, and β is an attenuation factor to take account of backscatter.
Snow albedo α s is assumed to vary linearly with temperature between that of cold, dry snow (α c ) at a threshold temperature T c , and that of melting snow (α m ) at the melting point, T m , where T m is fixed at 0 • C while T c , α c and α m can be varied.Figures

Back Close
Full The JULES land-surface scheme also calculates the sea-ice surface temperature and the atmosphere-to-ice fluxes as in HadGEM1 (see McLaren et al., 2006, for details).Within CICE these fluxes (downward latent heat flux, conductive flux through the ice, and surface heat flux), along with the ocean-ice heat flux (McPhee, 1992), determine the rate at which the ice grows or melts in each thickness category.This calculation also uses the enthalpy of each snow and ice layer.The enthalpy is defined as the negative of the energy required to melt a unit volume of ice/snow and raise its temperature to 0 • C. For the zero-layer thermodynamics used in the GSI configurations, the enthalpy is simply the negative product of the density and the latent heat of fusion.The calculated thermodynamic growth or melt rates are then used in the linear remapping scheme of Lipscomb (2001) to exchange the ice between thickness categories.

Dynamics and ridging
The ice velocities are calculated by solving the 2-D momentum equation for the force balance per unit area in the ice pack (Hibler, 1979), including terms for wind stress, ocean stress, internal ice stress, and stresses due to Coriolis effects.The internal ice stress is calculated using the elastic viscous plastic (EVP) scheme (Hunke and Dukowicz, 2002), which assumes the ice has a viscous plastic rheology, and incorporates an elastic wave modification to improve the computational efficiency.The GSI configurations use the Rothrock et al. (1975) formulation for ice strength.The sea ice is advected using the CICE incremental remapping scheme (Lipscomb and Hunke, 2004).2014) investigated the sensitivity of Arctic and Antarctic sea ice extent, thickness and volume in GSI4.0 to changes in several sea ice physical parameters, as well as to changes in the resolutions of the atmosphere and ocean models.By testing each of these sensitivities in isolation, they identified an optimum set of sea ice parameters for use in the Met Office coupled configuration.They found the Arctic sea ice to be most sensitive to changes in the albedos and thermal conductivities of ice and snow, while the Antarctic sea ice was most sensitive to changes in ice salinity, atmospheric and oceanic forcing, and ice-ocean model resolution.This forms the basis for the set of parameters used in GSI6.0, with some adjustments to account for the effect of changes in the atmosphere model made at the same time (see Walters et al., 2014).Parameter values are given in Table 1 and the CICE namelist used for GSI6.0 is included in Appendix A. The albedo parameters α m , f and β, were set in such a way as to increase the surface albedo, thereby reducing summer melt; the other albedo parameters were left unchanged.The values of the thermal conductivities of ice and snow, κ ice and κ snow , were chosen to increase the heat flux through the ice in autumn and winter, thereby increasing ice growth.The ice salinity, S, was increased, because Rae et al. (2014) found that this led to greater Antarctic ice growth due to a colder ocean mixed layer through the effect of salinity on ocean mixing.does not calculate any surface fluxes or the surface ice temperature.Note that setting calc_Tsfc to false also means that the albedo settings in the CICE namelist are irrelevant as the albedo is not calculated by CICE.Wind stresses are passed from the UM atmosphere rather than being calculated in CICE, so calc_strair is set to false.A constant value for the freezing point of sea water is used (1.8 • C), by setting Tfrzpt='constant'.This is required for consistency with the UM atmosphere-ice thermodynamics.The variable ns_boundary_type is set to tripole for the ORCA1 grid (i.e. in GSI4.0), indicating a tripolar grid with the "north fold" occurring along velocity points.The alternative setting tripoleT is used for the ORCA025 grid (i.e. in GSI6.0)where the north fold occurs along temperature points.The CICE cpp keys used in HadGEM3 at GC2.0 are shown in Table 2.The CICE namelist used in GSI6.0, which has been edited to detail the scientific options only, is given in Appendix A.

Experimental setup
We compare sea ice simulations from GSI6.0 (within GC2.0) with those from the previous configuration, GSI4.0 (within an earlier configuration of the coupled model).Both simulations were performed with a fully-coupled configuration of the Met Office's modelling system.The atmosphere and land-surface models were run on an N96 grid (equivalent to a resolution of 1.875 • in longitude and 1.25 • in latitude); the ocean and sea-ice models were on an ORCA1 grid (nominal 1 • resolution) for GSI4.0, and an ORCA025 grid (nominal 0.25 • resolution) for GSI6.0.The model setups and parameter values used are given in Table 1.Both simulations used initial conditions, greenhouse gas concentrations, and emissions of aerosols and their precursors appropriate for the present day (equivalent to year 2000).In both cases, we consider 50 years of output following an 80 year spin-up.Figures

Back Close
Full

Model evaluation
In GSI4.0, the Arctic ice volume (Fig. 2c; Table 3) was too low relative to that from the Pan Arctic Ice Ocean Modeling and Assimilation System (PIOMAS, see Schweiger et al., 2011), a coupled ice-ocean model that includes assimilation of observations.The findings of Rae et al. (2014), and the poor agreement of GSI4.0 Arctic sea ice with observational datasets, informed the choice of parameter values for GSI6.0 (see Table 1).In this section, the differences between GSI6.0 and GSI4.0 will be discussed, and put in the context of the findings of Rae et al. (2014).

Arctic
In GSI6.0, we see thickening of the Arctic ice pack at the end of winter relative to GSI4.0 (Fig. 1a and b), resulting in improved agreement with observations (see Fig. 1 of Laxon et al., 2013).Rae et al. (2014) found that increased conductivities of both ice and snow led to an increased upward conductive heat flux through the Arctic sea ice in late summer and early autumn, leading to reduced basal melt in July and August, and increased basal growth in winter.This in turn led to thicker ice in the Arctic in winter.
We also see an increase in summer ice extent, thickness and volume in GSI6.0 compared to GSI4.0 (Figs. 1c and d, 3d and e, 2a and c; Table 3).This mirrors the behaviour seen by Rae et al. (2014) with increased ice and snow thermal conductivities, where the increased ice thickness seen in winter persisted through the following melt season.In addition to this, Rae et al. ( 2014) also found that in the Arctic increased snow albedo led to reduced surface melt in summer, and thus to increased summer ice extent, thickness and volume.It is likely that similar effects are occurring here in GSI6.0.The summer ice concentration and extent are now more in agreement with the HadISST dataset of Rayner et al. (2003) (Figs. 3f, 2a), and the agreement of the volume with PIOMAS has also improved (Fig. 2c; Table 3).
In winter, there are also overall improvements in the total extent relative to HadISST (Fig. 2a; Table 3), largely due to reduced ice cover in the Labrador Sea (Fig. 3a-c).Figures

Back Close
Full The investigations of Rae et al. (2014) suggest that this is attributable to the increased ice-ocean model resolution.They found that the increased resolution led to warmer sea surface temperatures in the Labrador sea, leading in turn to a reduced sea ice concentration there, and thus to a lower total Arctic winter sea ice extent.Despite this reduced winter ice extent, the increased ice thickness has led to an increased ice volume, with the result that it is now more in agreement with that from PIOMAS (Fig. 2c; Table 3).

Antarctic
The GC simulations have been found to display a warm bias in sea-surface temperatures (SST) in the Southern Ocean (Megann et al., 2014), due to a positive bias in downward heat flux from the atmosphere into the ocean (Williams et al., 2014).In GSI4.0, this led to a low Antarctic sea ice extent in Austral summer, although the winter ice extent compared favourably with HadISST (Fig. 2b; Table 3).Rae et al. (2014) found that the Antarctic ice extent and volume were generally insensitive to perturbations in the ice physics parameters (other than salinity), but that the effects of the warm SST bias were exacerbated at higher ice-ocean resolution.They attributed this to the removal of the Gent-McWilliams eddy parameterisation at the eddy-permitting resolution of ORCA025.It is thought that this parameterisation helps to mask the warm bias at lower resolution, but that its removal in the higher-resolution runs leads to increased southward heat transport in the ocean.
As discussed in Sect.3, GSI6.0 is run at the higher resolution of ORCA025 (see Table 1).The exacerbation of the warm bias in the Southern Ocean therefore has an impact on the Antarctic sea ice in GSI6.0, and there is a substantial reduction in ice extent and volume in all seasons (Fig. 2b, d; Table 3).Thus, while the transition from GSI4.0 to GSI6.0 leads to some improvements in the Arctic, the same is not true in the Antarctic.Work is ongoing to resolve the warm bias in the Southern Ocean, and it is anticipated that this will lead to improved simulations of Antarctic sea ice in future configurations.Figures

Back Close
Full

Conclusions
We have described and evaluated the new Global Sea Ice configuration, GSI6.0, run within the Met Office Global Coupled model configuration GC2.0.The choice of parameters for GSI6.0 was informed by the work of Rae et al. (2014), who conducted an extensive sea ice parameter sensitivity study within the Met Office coupled modelling system and in addition isolated the impact of ice physics changes from that of forcing and resolution changes.In the new configurations, the values of several sea ice parameters have been changed, and the ice-ocean model resolution has been increased from ORCA1 (nominal 1 • resolution) to ORCA025 (nominal 0.25 • resolution).This has resulted in thicker Arctic ice in all seasons, and Arctic ice concentration and extent that agree better with the HadISST observational dataset (Rayner et al., 2003).In the Antarctic, the higher ice-ocean model resolution has resulted in the exacerbation of an existing warm bias in the Southern Ocean.This has in turn led to a large reduction in ice extent and volume.Rectification of this bias will require further development work on atmosphere-ocean heat transfer in the coupled model.
While the sea ice simulation in GSI6.0 represents an improvement over that in GSI4.0 -at least in the Arctic -there are still several areas in which there is potential for further model enhancement.First, while the GSI configurations use five ice thickness categories in the CICE model, the sea ice surface fluxes of latent and sensible heat are calculated in the JULES land-surface scheme as gridbox means.In the next configuration, these calculations will be performed on all five thickness categories.Second, the sea ice surface albedo scheme used in GSI4.0 and GSI6.0 is the same broadband scheme used in HadGEM1 (McLaren et al., 2006).The next configuration will include separate calculations for four radiation bands -direct and diffuse radiation for both visible and near-infrared bands -as well as for each ice thickness category.
It is anticipated that future configurations will also include an explicit representation of the effect of melt ponds on surface albedo.As mentioned in Sect.ture configurations, this freezing temperature will be calculated as a function of ocean salinity.Finally, as mentioned in Sect.2.1, the current GSI configurations use the zerolayer thermodynamics of Semtner (1976Semtner ( , 1987)), rather than the full multi-layer CICE scheme.Planned modifications to CICE, the UM, and JULES will enable the CICE multilayer model to be used with the UM atmosphere in the future.The versions and revisions of each model used in this paper are given in Table 1.A number of branches are applied to these codes.Please contact the authors for more information on these branches and how to obtain them.Introduction

Conclusions References
Tables Figures

Back Close
Full  Full  Full Discussion Paper | Discussion Paper | Discussion Paper | The CICE sea ice model solves this equation to determine the evolution of g in time and space.Full details of the model are available in the CICE user manual (Hunke and Lipscomb, 2010); here we summarise the main features of the model used in GSI, and detail the specific settings and choices for the previous configuration (GSI4.0)and the new configuration (GSI6.0).Much of the basic model description is as in Appendix D of Hewitt et al. (2011), but it is reproduced here for completenessDiscussion Paper | Discussion Paper | Discussion Paper | fects of scattering.This is the same scheme used in HadGEM1 (McLaren et al., 2006), HadGEM2 (HadGEM2 Development Team et al., 2011) and HadGEM3 (Hewitt et al., Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | The mechanical redistribution (or ridging) scheme in CICE converts thinner ice to thicker ice and open water, and is applied after the advection of ice.When the ice is converging, enough ridging takes place such that the ice area does not exceed the grid-cell area.The scheme is based on work byThorndike et al. (1975);Hibler (1980);Flato andHibler (1995), andRothrock et al. (1975).It favours the closing of open water and ridging of the thinnest ice over the ridging of thicker ice.In GSI the ridging participation function suggested byLipscomb et al. (2007) is used.The ridged ice is then distributed between thickness categories assuming an exponential ITD(Lipscomb et al., 2007).Discussion Paper | Discussion Paper | Discussion Paper | 2.4 CICE settings used for GSI6.0 Rae et al. ( Rae et al. (2014) found the Arctic and Antarctic sea ice extent and volume to be relatively insensitive to the value of the ridging parameter µ rdg(Hunke, 2010); however, the value was reduced from 4 to 3 m 1/2 as this is now the recommended value.The roughness lengths of pack ice and the marginal ice zone, z 0 (ice) and z 0 (MIZ), previously had different values in the climate and Numerical Weather Prediction (NWP) configurations of the model.In GSI6, the values in the climate configuration have been increased to make them consistent with those in the NWP configuration.For coupling with the UM atmosphere, heat_capacity and calc_Tsfc are both set to false.This means that zero-layer thermodynamics are used and that CICE Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | 2.4, GSI currently uses a fixed reference value of −1.8 • C for the freezing temperature of sea water.In fu-Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Global Land 3.0/3.1 configurations, Geosci.Model Dev., 4, 919-941, doi:10.5194/gmd-4-919-2011,2011.2530, 2535 Wettlaufer, J. S.: Heat flux at the ice-ocean interface, J. Geophys.Res., 96, 7215-7236, 1991.Williams, K. D., Harris, C. M., Bodas-Salcedo, A., Camp, J., Comer, R. E., Copsey, D., Fereday, D., Graham, T., Hill, R., Hinton, T., Hyder, P., Ineson, S., Masato, G., Milton, S. F., Roberts, Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Table 1 .
Model setup and values of sea ice parameters in GSI4.0 and GSI6.0 simulations.

Table 2 .
The preprocessor keys used for CICE in GC1.0-GSI6.0.available for input and output files CICE_IN_NEMO CICE is run in the NEMO environment.CICE is called from the NEMO surface module which also exchanges the coupling fields between NEMO and CICE ORCA_GRID Controls reading in grid, land masks and forcing data on the ORCA family of grids key_oasis3 Coupling uses OASIS3 REPRODUCIBLE Ensures global sums bit compare for parallel model runs with different grid decompositions Introduction

Table 3 .
Modelled and observationally-based sea ice extent and volume.