ALBERT

Description: This study assesses the impact of different sea ice thickness distribution (ITD) discretizations on the sea ice concentration (SIC) variability in ocean stand-alone NEMO3.6–LIM3 simulations. Three ITD discretizations with different numbers of sea ice thickness categories and boundaries are evaluated against three different satellite products (hereafter referred to as “data”). Typical model and data interannual SIC variability is characterized by K-means clustering both in the Arctic and Antarctica between 1979 and 2014. We focus on two seasons, winter (January–March) and summer (August–October), in which correlation coefficients across clusters in individual months are largest. In the Arctic, clusters are computed before and after detrending the series with a second-degree polynomial to separate interannual from longer-term variability. The analysis shows that, before detrending, winter clusters reflect the SIC response to large-scale atmospheric variability at both poles, while summer clusters capture the negative and positive trends in Arctic and Antarctic SIC, respectively. After detrending, Arctic clusters reflect the SIC response to interannual atmospheric variability predominantly. The cluster analysis is complemented with a model–data comparison of the sea ice extent and SIC anomaly patterns. The single-category discretization shows the worst model–data agreement in the Arctic summer before detrending, related to a misrepresentation of the long-term melting trend. Similarly, increasing the number of thin categories reduces model–data agreement in the Arctic, due to a poor representation of the summer melting trend and an overly large winter sea ice volume associated with a net increase in basal ice growth. In contrast, more thin categories improve model realism in Antarctica, and more thick ones improve it in central Arctic regions with very thick ice. In all the analyses we nonetheless identify no optimal discretization. Our results thus suggest that no clear benefit in the representation of SIC variability is obtained from increasing the number of sea ice thickness categories beyond the current standard with five categories in NEMO3.6–LIM3.

Description: The ice thickness distribution (ITD) is one of the core constituents of modern sea ice models. The ITD accounts for the unresolved spatial variability of sea ice thickness within each model grid cell. While there is a general consensus on the added physical realism brought by the ITD, how to discretize it remains an open question. Here, we use the ocean–sea ice general circulation model, Nucleus for European Modelling of the Ocean (NEMO) version 3.6 and Louvain-la-Neuve sea Ice Model (LIM) version 3 (NEMO3.6-LIM3), forced by atmospheric reanalyses to test how the ITD discretization (number of ice thickness categories, positions of the category boundaries) impacts the simulated mean Arctic and Antarctic sea ice states. We find that winter ice volumes in both hemispheres increase with the number of categories and attribute that increase to a net enhancement of basal ice growth rates. The range of simulated mean winter volumes in the various experiments amounts to ∼30 % and ∼10 % of the reference values (run with five categories) in the Arctic and Antarctic, respectively. This suggests that the way the ITD is discretized has a significant influence on the model mean state, all other things being equal. We also find that the existence of a thick category with lower bounds at ∼4 and ∼2 m for the Arctic and Antarctic, respectively, is a prerequisite for allowing the storage of deformed ice and therefore for fostering thermodynamic growth in thinner categories. Our analysis finally suggests that increasing the resolution of the ITD without changing the lower limit of the upper category results in small but not negligible variations of ice volume and extent. Our study proposes for the first time a bi-polar process-based explanation of the origin of mean sea ice state changes when the ITD discretization is modified. The sensitivity experiments conducted in this study, based on one model, emphasize that the choice of category positions, especially of thickest categories, has a primary influence on the simulated mean sea ice states while the number of categories and resolution have only a secondary influence. It is also found that the current default discretization of the NEMO3.6-LIM3 model is sufficient for large-scale present-day climate applications. In all cases, the role of the ITD discretization on the simulated mean sea ice state has to be appreciated relative to other influences (parameter uncertainty, forcing uncertainty, internal climate variability).

Description: The Ice Thickness Distribution (ITD) is one of the core constituents of modern sea ice models. The ITD accounts for the unresolved spatial variability of sea ice thickness within each model grid cell. While there is a general consensus on the added physical realism brought by the ITD, how to implement it remains an open question. Here, we use the ocean--sea ice general circulation model NEMO3.6-LIM3 forced by atmospheric reanalyses to test how the ITD discretization (number of ice thickness categories, positions of the category boundaries) impacts the simulated mean Arctic and Antarctic sea ice states. We find that winter ice volumes in both hemispheres increase with the number of categories, and attribute that increase to a net enhancement of basal ice growth rates. The range of simulated mean winter volumes in the various experiments amounts to ~ 30 % and ~ 10 % of the reference values (run with 5 categories) in the Arctic and Antarctic, respectively. This suggests that the way the ITD is discretized has a significant influence on the model mean state, all other things being equal. We also find that the existence of a thick category with lower bounds at ~ 4 m and ~ 2 m for the Arctic and Antarctic, respectively, is a prerequisite for allowing the storage of deformed ice, and therefore for fostering thermodynamic growth in thinner categories. Our analysis finally suggests that increasing the resolution of the ITD without changing the lower limit of the upper category results in small but not negligible variations of ice volume and extent. Our study proposes for the first time a bi-polar process-based explanation of the origin of mean state changes when the ITD discretization is modified. The sensitivity experiments conducted in this study, based on one model, emphasize that the choice of category positions, especially of thickest categories, has a primary influence on the simulated mean sea ice states while the number of categories and resolution have only a secondary influence.

Description: In this paper, we present and evaluate the skill of an EC-Earth3.3 decadal prediction system contributing to the Decadal Climate Prediction Project – Component A (DCPP-A). This prediction system is capable of skilfully simulating past global mean surface temperature variations at interannual and decadal forecast times as well as the local surface temperature in regions such as the tropical Atlantic, the Indian Ocean and most of the continental areas, although most of the skill comes from the representation of the external radiative forcings. A benefit of initialization in the predictive skill is evident in some areas of the tropical Pacific and North Atlantic oceans in the first forecast years, an added value that is mostly confined to the south-east tropical Pacific and the eastern subpolar North Atlantic at the longest forecast times (6–10 years). The central subpolar North Atlantic shows poor predictive skill and a detrimental effect of initialization that leads to a quick collapse in Labrador Sea convection, followed by a weakening of the Atlantic Meridional Overturning Circulation (AMOC) and excessive local sea ice growth. The shutdown in Labrador Sea convection responds to a gradual increase in the local density stratification in the first years of the forecast, ultimately related to the different paces at which surface and subsurface temperature and salinity drift towards their preferred mean state. This transition happens rapidly at the surface and more slowly in the subsurface, where, by the 10th forecast year, the model is still far from the typical mean states in the corresponding ensemble of historical simulations with EC-Earth3. Thus, our study highlights the Labrador Sea as a region that can be sensitive to full-field initialization and hamper the final prediction skill, a problem that can be alleviated by improving the regional model biases through model development and by identifying more optimal initialization strategies.