Assessment of sea ice thickness simulations in the CMIP6 models with CICE components

Arctic sea ice plays a critical role in modulating our global climate system and the exchange of heat fluxes in the polar region, but its impact on climate varies across different sea ice thickness (SIT) categories. Compared to sea ice cover, the performance of ice models in simulating SIT has been less evaluated, particularly in the sixth Coupled Model Intercomparison Project Phase (CMIP6). Here, we chose 12 CMIP6 models with the Community Ice Code model (CICE) components and compared their SIT simulations with the satellite observations and the Pan-Arctic Ice Ocean Modeling and Assimilation System (PIOMAS) model between 1980 and 2014. Our results show that the seasonal cycle of the PIOMAS SIT is consistent with satellite observations. Compared to the PIOMAS reanalysis, the multi-model ensemble mean (MME) well represents the sea ice extent in both the thin ice (<0.6 m) and thick ice (> 3.6 m). However, the MME SIT has larger biases in the Chukchi Sea, the Beaufort Sea, the central Arctic, and the Greenland Sea during winter and mainly in the central Arctic during summer. Both the MME and PIOMAS show decreasing trends in SIT over the entire Arctic Ocean in all seasons, but the interannual variability of SIT in MME is smaller than that in PIOMAS. Among the 12 CMIP6 models, the FIO-ESM-2.0 model shows the best simulation of the annual mean SIT, but the SAM0-UNICON and NESM3 models have the largest biases in the climatological mean SIT over the Arctic Ocean. We also demonstrate that the FIO-ESM-2.0 performs the best in the seasonal cycles of SIT. Our study suggests that more attention should be paid to the coupling of the CICE model with ocean and atmosphere models, which is vital to improving the SIT simulation in CMIP6 models and to better understanding the impact of Arctic sea ice on our climate system.