scholarly journals On the discretization of the ice thickness distribution in the NEMO3.6-LIM3 global ocean–sea ice model

2019 ◽  
Author(s):  
François Massonnet ◽  
Antoine Barthélemy ◽  
Koffi Worou ◽  
Thierry Fichefet ◽  
Martin Vancoppenolle ◽  
...  

Abstract. 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.

2019 ◽  
Vol 12 (8) ◽  
pp. 3745-3758 ◽  
Author(s):  
François Massonnet ◽  
Antoine Barthélemy ◽  
Koffi Worou ◽  
Thierry Fichefet ◽  
Martin Vancoppenolle ◽  
...  

Abstract. 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).


2020 ◽  
Vol 13 (10) ◽  
pp. 4773-4787
Author(s):  
Eduardo Moreno-Chamarro ◽  
Pablo Ortega ◽  
François Massonnet

Abstract. 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.


2021 ◽  
Author(s):  
Yushi Morioka ◽  
Doroteaciro Iovino ◽  
Andrea Cipollone ◽  
Simona Masina ◽  
Swadhin Behera

<p>Skillful sea-ice prediction in the Antarctic Ocean remains a big challenge due to paucity of sea-ice observations and insufficient representation of sea-ice processes in climate models. This study demonstrates that the Antarctic sea-ice concentration (SIC) prediction is significantly improved using a coupled general circulation model (SINTEX-F2) in which the model’s SIC and sea-ice thickness (SIT) are initialized with the ocean/sea-ice reanalysis product (C-GLORSv7). It is found that the wintertime SIT initialization adds positive values to the prediction skills of the summertime SIC, most effectively in the Weddell Sea where the SIT climatology and variability are the largest among the Antarctic Seas. Examination of the SIT balance during low sea-ice years of the Weddell Sea shows that negative SIT anomalies initialized in June retain the memory throughout austral winter (July-September) owing to horizontal advection of the SIT anomalies by sea-ice velocities. The negative SIT anomalies continue to develop in austral spring (October-December) owing to more incoming solar radiation via ice-albedo feedback and the associated warming of mixed layer. This results in further sea-ice decrease during austral summer (January-March). Concomitantly, the model reasonably reproduces atmospheric circulation anomalies in the Amundsen-Bellingshausen Seas as well as the Weddell Sea during the development of the negative sea-ice anomalies. These results provide solid evidence that the wintertime SIT initialization benefits skillful summertime sea-ice prediction in the Antarctic Seas.</p>


2013 ◽  
Vol 7 (1) ◽  
pp. 441-473 ◽  
Author(s):  
L. Rabenstein ◽  
T. Krumpen ◽  
S. Hendricks ◽  
C. Koeberle ◽  
C. Haas ◽  
...  

Abstract. A combined interpretation of synthetic aperture radar (SAR) satellite images and helicopter electromagnetic (HEM) sea-ice thickness data has provided an estimate of sea-ice volume formed in Laptev Sea polynyas during the winter of 2007/08. The evolution of the surveyed sea-ice areas, which were formed between late December 2007 and middle April 2008, was tracked using a series of SAR images with a sampling interval of 2–3 days. Approximately 160 km of HEM data recorded in April 2008 provided sea-ice thicknesses along profiles that transected sea-ice varying in age from 1–116 days. For the volume estimates, thickness information along the HEM profiles was extrapolated to zones of the same age. The error of areal mean thickness information was estimated to be between 0.2 m for younger ice and up to 1.55 m for older ice, with the primary error source being the spatially limited HEM coverage. Our results have demonstrated that the modal thicknesses and mean thicknesses of level ice correlated with the sea-ice age, but that varying dynamic and thermodynamic sea-ice growth conditions resulted in a rather heterogeneous sea-ice thickness distribution on scales of tens of kilometers. Taking all uncertainties into account, total sea-ice area and volume produced within the entire surveyed area were 52 650 km2 and 93.6 ± 26.6 km3. The surveyed polynya contributed 2.0 ± 0.5% of the sea-ice produced throughout the Arctic during the 2007/08 winter. The SAR-HEM volume estimate compares well with the 112 km3 ice production calculated with a high resolution ocean sea-ice model. Measured modal and mean-level ice thicknesses correlate with calculated freezing-degree-day thicknesses with a factor of 0.87–0.89, which was too low to justify the assumption of homogeneous thermodynamic growth conditions in the area, or indicates a strong dynamic thickening of level ice by rafting of even thicker ice.


2011 ◽  
Vol 5 (1) ◽  
pp. 131-167
Author(s):  
A. Oikkonen ◽  
J. Haapala

Abstract. Changes of the mean sea ice thickness and concentration in the Arctic are well known. However, comparable little is known about the ice thickness distribution and the composition of ice pack in quantity. In this paper we determine the ice thickness distributions, mean and modal thicknesses, and their regional and seasonal variability in the Arctic under different large scale atmospheric circulation modes. We compare characteristics of the Arctic ice pack during the periods 1975–1987 and 1988–2000, which have a different distribution in the AO/DA space. The study is based on submarine measurements of sea ice draft. The prevalent feature is that the peak of sea ice thickness distributions has generally taken a narrower form and shifted toward thinner ice. Also, both mean and modal ice thickness have generally decreased. These noticeable changes result from a loss of thick, mostly deformed, ice. In the spring the loss of the volume of ice thicker than 5 m exceeds 35% in all regions except the Nansen Basin, and the reduction is as much as over 45% at the North Pole and in the Eastern Arctic. In the autumn the volume of thick, mostly deformed ice has decreased by more than 40% in the Canada Basin only, but the reduction is more than 30% also in the Beaufort Sea and in the Chukchi Sea. In the Beaufort Sea region the decrease of the modal draft has been so strong that the peak has shifted from multiyear ice to first-year type ice. Also, the regional and seasonal variability of the sea ice thickness has decreased, since the thinning has been the most pronounced in the regions with the thickest pack ice (the Western Arctic), and during the spring (0.6–0.8 m per decade).


2021 ◽  
Author(s):  
Nicholas Williams ◽  
Nicholas Byrne ◽  
Daniel Feltham ◽  
Peter Jan Van Leeuwen ◽  
Ross Bannister ◽  
...  

<div><span>A modified, standalone version of the Los Alamos Sea Ice Model (CICE) has been coupled to the Parallelized Data Assimilation Framework (PDAF) to produce a new Arctic sea ice data assimilation system CICE-PDAF, with routines for assimilating many types of recently developed sea ice observations. In this study we explore the effects of assimilating a sub-grid scale sea ice thickness distribution derived from Cryosat-2 Arctic sea ice estimates into CICE-PDAF. The true state of the sub-grid scale ice thickness distribution is not well established, and yet it plays a key role in large scale sea ice models and is vital to the dynamical and thermodynamical processes necessary to produce a good representation of the Arctic sea ice state. We examine how assimilating sub-grid scale sea ice thickness distributions can affect the evolution of the sea ice state in CICE-PDAF and better our understanding of the Arctic sea ice system.</span></div>


2016 ◽  
Author(s):  
Petteri Uotila ◽  
Dorotea Iovino ◽  
Martin Vancoppenolle ◽  
Mikko Lensu ◽  
Clement Rousset

Abstract. A set of hindcast simulations with the new NEMO3.6 ocean-ice model in the ORCA1 grid and forced by the DFS5.2 atmospheric data was performed from 1958–2012. We focussed on simulations that differ only in their sea-ice component: the old standard version LIM2 and its successor LIM3. Main differences between these sea-ice models are the parameterisations of sub-grid-scale sea-ice thickness distribution, ice deformation, thermodynamic processes, and sea-ice salinity. Our main objective was to diagnose the ocean-ice sensitivity to the updated NEMO-LIM sea-ice physics. Results of such analysis have not been published for this new NEMO version. In the polar regions, NEMO-LIM3 compares better with observations, while NEMO-LIM2 deviates more, producing too much ice in the Arctic, for example. Differences between NEMO-LIM2 and NEMO-LIM3 do not change in simulations even when the freshwater adjustments are turned off. In the extra-polar regions, the oceanographic conditions of the two NEMO-LIM versions remain relatively similar, although they slowly drift apart over decades. A simplified NEMO-LIM3 configuration, having a virtual, single-category sea-ice thickness distribution, produced sea ice with a skill sufficient for ocean-ice hindcasts that target oceanographic studies. We conclude that NEMO3.6 is ready to be used as a stand-alone ocean-ice model and as a component of coupled atmosphere-ocean models.


Sign in / Sign up

Export Citation Format

Share Document