The impact of assimilating weekly CryoSat-2 sea ice thickness data together with daily SMOS sea ice thickness and daily SSMIS sea ice concentration data on the sea ice fields of a coupled sea ice-ocean model of the Arctic Ocean is investigated. The sea-ice model is based on the Massachusetts Institute of Technology general circulation model (MITgcm) and the assimilation is performed by a localized Singular Evolutive Interpolated Kalman (LSEIK) filter coded in the Parallel Data Assimilation Framework (PDAF). A period of three months from 1 November 2011 to 30 January 2012 is selected to assess the skill of the assimilation system in the cold season. Compared to the unassimilated solution and a solution where only sea ice concentration is assimilated, the model-data misfits are substantially reduced in areas of both thick and thin ice. The sea ice thickness estimates agree significantly better with in situ observations in the central Arctic Ocean than the sea ice thickness obtained from assimilating SMOS data alone, while the sea ice concentration shows very small improvements. The sea ice fields obtained by the joint assimilation of SMOS and CryoSat-2 data also have lower errors in thickness and concentration than those obtained from directly assimilating a statistically merged SMOS and CryoSat-2 sea ice thickness product. These lower errors suggest that model dynamics play a significant role in data blending.

Exploiting the complementary character of CryoSat‐2 and Soil Moisture and Ocean Salinity satellite sea ice thickness products, daily Arctic sea ice thickness estimates from October 2010 to December 2016 are generated by an Arctic regional ice‐ocean model with satellite thickness assimilated. The assimilation is performed by a Local Error Subspace Transform Kalman filter coded in the Parallel Data Assimilation Framework. The new estimates can be generally thought of as combined model and satellite thickness (CMST). It combines the skill of satellite thickness assimilation in the freezing season with the model skill in the melting season, when neither CryoSat‐2 nor Soil Moisture and Ocean Salinity sea ice thickness is available. Comparisons with in situ observations from the Beaufort Gyre Exploration Project, Ice Mass Balance Buoys, and the NASA Operation IceBridge demonstrate that CMST reproduces most of the observed temporal and spatial variations. Results also show that CMST compares favorably to the Pan‐Arctic Ice‐Ocean Modeling and Assimilation System product and even appears to correct known thickness biases in the Pan‐Arctic Ice‐Ocean Modeling and Assimilation System. Due to imperfect parameterizations in the sea ice model and satellite thickness retrievals, CMST does not reproduce the heavily deformed and ridged sea ice along the northern coast of the Canadian Arctic Archipelago and Greenland. With the new Arctic sea ice thickness estimates sea ice volume changes in recent years can be further assessed.

A new global climate model setup using FESOM2.0 for the sea ice-ocean component and ECHAM6.3 for the atmosphere and land surface has been developed. Replacing FESOM1.4 by FESOM2.0 promises a higher efficiency of the new climate setup compared to its predecessor. The new setup allows for long-term climate integrations using a locally eddy-resolving ocean. Here it is evaluated in terms of (1) the mean state and long-term drift under preindustrial climate conditions, (2) the fidelity in simulating the historical warming, and (3) differences between coarse and eddy-resolving ocean configurations. The results show that the realism of the new climate setup is overall within the range of existing models. In terms of oceanic temperatures, the historical warming signal is of smaller amplitude than the model drift in case of a relatively short spin-up. However, it is argued that the strategy of "de-drifting" climate runs after the short spin-up, proposed by the HighResMIP protocol, allows one to isolate the warming signal. Moreover, the eddy-permitting/resolving ocean setup shows notable improvements regarding the simulation of oceanic surface temperatures, in particular in the Southern Ocean.

This paper describes and evaluates the assimilation component of a seamless sea ice prediction system, which is developed based on the fully coupled Alfred Wegener Institute, Helmholtz Center for Polar and Marine Research Climate Model (AWI-CM, v1.1). Its ocean/ice component with unstructured-mesh discretization and smoothly varying spatial resolution enables seamless sea ice prediction across a wide range of space and time scales. The model is complemented with the Parallel Data Assimilation Framework to assimilate observations in the ocean/ice component with an Ensemble Kalman Filter. The focus here is on the data assimilation of the prediction system. First, the performance of the system is tested in a perfect-model setting with synthetic observations. The system exhibits no drift for multivariate assimilation, which is a prerequisite for the robustness of the system. Second, real observational data for sea ice concentration, thickness, drift, and sea surface temperature are assimilated. The analysis results are evaluated against independent in situ observations and reanalysis data. Further experiments that assimilate different combinations of variables are conducted to understand their individual impacts on the model state. In particular, assimilating sea ice drift improves the sea ice thickness estimate, and assimilating sea surface temperature is able to avert a circulation bias of the free-running model in the Arctic Ocean at middepth. Finally, we present preliminary results obtained with an extended system where the atmosphere is constrained by nudging toward reanalysis data, revealing challenges that still need to be overcome to adapt the ocean/ice assimilation. We consider this system a prototype on the way toward strongly coupled data assimilation across all model components.Plain Language Summary Sea ice prediction over seasonal time scale has attracted the focus from both the scientific and socioeconomic communities recently. We develop a system aiming for seamless sea ice prediction using a coupled model that is equipped with an unstructured ocean/sea ice component. The high-resolution mesh over the polar regions allows us to explore its possible benefits on the prediction across a wide range of time scales. To this end, the sea surface temperature, sea ice concentration, sea ice thickness, and sea ice drift observations are assimilated in the current system to diagnose the performance of the model initialization for future forecasts.

Abstract. Data assimilation integrates information from observational measurements with numerical models. When used with coupled models of Earth system compartments, e.g., the atmosphere and the ocean, consistent joint states can be estimated. A common approach for data assimilation is ensemble-based methods which utilize an ensemble of state realizations to estimate the state and its uncertainty. These methods are far more costly to compute than a single coupled model because of the required integration of the ensemble. However, with uncoupled models, the ensemble methods also have been shown to exhibit a particularly good scaling behavior. This study discusses an approach to augment a coupled model with data assimilation functionality provided by the Parallel Data Assimilation Framework (PDAF). Using only minimal changes in the codes of the different compartment models, a particularly efficient data assimilation system is generated that utilizes parallelization and in-memory data transfers between the models and the data assimilation functions and hence avoids most of the file reading and writing, as well as model restarts during the data assimilation process. This study explains the required modifications to the programs with the example of the coupled atmosphere–sea-ice–ocean model AWI-CM (AWI Climate Model). Using the case of the assimilation of oceanic observations shows that the data assimilation leads only to small overheads in computing time of about 15 % compared to the model without data assimilation and a very good parallel scalability. The model-agnostic structure of the assimilation software ensures a separation of concerns in which the development of data assimilation methods can be separated from the model application.

The Alfred Wegener Institute Climate Model (AWI‐CM) participates for the first time in the Coupled Model Intercomparison Project (CMIP), CMIP6. The sea ice‐ocean component, FESOM, runs on an unstructured mesh with horizontal resolutions ranging from 8 to 80 km. FESOM is coupled to the Max Planck Institute atmospheric model ECHAM 6.3 at a horizontal resolution of about 100 km. Using objective performance indices, it is shown that AWI‐CM performs better than the average of CMIP5 models. AWI‐CM shows an equilibrium climate sensitivity of 3.2°C, which is similar to the CMIP5 average, and a transient climate response of 2.1°C which is slightly higher than the CMIP5 average. The negative trend of Arctic sea‐ice extent in September over the past 30 years is 20–30% weaker in our simulations compared to observations. With the strongest emission scenario, the AMOC decreases by 25% until the end of the century which is less than the CMIP5 average of 40%. Patterns and even magnitude of simulated temperature and precipitation changes at the end of this century compared to present‐day climate under the strong emission scenario SSP585 are similar to the multi‐model CMIP5 mean. The simulations show a 11°C warming north of the Barents Sea and around 2°C to 3°C over most parts of the ocean as well as a wetting of the Arctic, subpolar, tropical, and Southern Ocean. Furthermore, in the northern middle latitudes in boreal summer and autumn as well as in the southern middle latitudes, a more zonal atmospheric flow is projected throughout the year.

An ensemble-based data assimilation framework for a coupled oceanatmosphere model is applied to investigate the influence of assimilating different types of ocean observations on the ocean and atmosphere simulation. The data assimilation is performed with the parallel data assimilation framework (PDAF) for the climate model AWI-CM. Observations of the ocean, namely satellite sea-surface temperature (SST) and temperature and salinity profiles, are assimilated into the ocean component. The atmospheric state is only influenced by the model dynamics. Different assimilation scenarios were carried out with different combinations of observations to investigate to what extent the assimilation into the coupled model leads to a better estimation of the state of the ocean as well as the atmosphere. The influence of the data assimilation is assessed by comparing the ocean prediction with dependent and independent ocean observations. For the atmosphere, the assimilation result is compared with the ERA-Interim atmospheric reanalysis data. The ocean temperature and salinity are improved by all the assimilation scenarios in the coupled system. The assimilation leads to a response of the atmosphere throughout the troposphere and impacts the global atmospheric circulation. Globally the temperature and wind speed are improved in the atmosphere on average. K E Y W O R D S coupled model, data assimilation, sea-surface temperature, temperature and salinity profiles 1 INTRODUCTION Traditionally, different components of the Earth system such as the ocean and the atmosphere are simulated by separate models with influences of other components being modelled as boundary conditions or forcings. However, the oceans and the atmosphere are connected and interact with each other. A consistent initial condition for these different components is required and is expected to provide a better forecast for both the ocean and the atmosphere. Earth system models simulate different components like the ocean, atmosphere, sea ice and land This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

No abstract