We introduce PARASO, a novel five-component fully coupled regional climate model over an Antarctic circumpolar domain covering the full Southern Ocean. The state-of-the-art models used are the fast Elementary Thermomechanical Ice Sheet model (f.ETISh) v1.7 (ice sheet), the Nucleus for European Modelling of the Ocean (NEMO) v3.6 (ocean), the Louvain-la-Neuve sea-ice model (LIM) v3.6 (sea ice), the COnsortium for Small-scale MOdeling (COSMO) model v5.0 (atmosphere) and its CLimate Mode (CLM) v4.5 (land), which are here run at a horizontal resolution close to $\mathrm{1}/\mathrm{4}{}^{\circ }$ . One key feature of this tool resides in a novel two-way coupling interface for representing ocean–ice-sheet interactions, through explicitly resolved ice-shelf cavities. The impact of atmospheric processes on the Antarctic ice sheet is also conveyed through computed COSMO-CLM–f.ETISh surface mass exchange. In this technical paper, we briefly introduce each model's configuration and document the developments that were carried out in order to establish PARASO. The new offline-based NEMO–f.ETISh coupling interface is thoroughly described. Our developments also include a new surface tiling approach to combine open-ocean and sea-ice-covered cells within COSMO, which was required to make this model relevant in the context of coupled simulations in polar regions. We present results from a 2000–2001 coupled 2-year experiment. PARASO is numerically stable and fully operational. The 2-year simulation conducted without fine tuning of the model reproduced the main expected features, although remaining systematic biases provide perspectives for further adjustment and development.

The Southern Hemisphere cryosphere has recently shown regionally-contrasted responses to climate change, in particular to the positive phases of the Southern Annular Mode. However, the understanding of the impacts of this mode on ice-shelf basal melt at a circum-Antarctic scale is still limited. Here, we performed idealized experiments with a pan-Antarctic regional ice-shelf cavity-resolving ocean—sea-ice model for different phases of the Southern Annular Mode. We show that positive phases lead to increased upwelling and subsurface ocean temperature and salinity close to ice shelves, while the opposite occurs for negative phases. A one-standard-deviation increase of the Southern Annular Mode leads to a net basal mass loss of 40 Gt yr⁻¹, with strong regional contrasts: increased ice-shelf basal melt in the Bellingshausen and Western Pacific sectors and the opposite response in the Amundsen sector. Estimates of 1000–1200 and 2090–2100 ice-shelf basal melt changes due to the Southern Annular Mode are −86.6 Gt yr⁻¹and 55.0 to 164.9 Gt yr⁻¹, respectively, compared to the present.

The seasonal cycle of the Antarctic sea ice extent is strongly asymmetric, with a relatively slow increase after the summer minimum followed by a more rapid decrease after the winter maximum. This cycle is intimately linked to the seasonal cycle of the insolation received at the top of the atmosphere, but sea ice processes as well as the exchanges with the atmosphere and ocean may also play a role. To quantify these contributions, a series of idealized sensitivity experiments have been performed with an eddy-permitting ($\mathrm{1}/\mathrm{4}$ ^∘) NEMO-LIM3 (Nucleus for European Modelling of the Ocean–Louvain-la-Neuve sea ice model version 3) Southern Ocean configuration, including a representation of ice shelf cavities, in which the model was either driven by an atmospheric reanalysis or coupled to the COSMO-CLM² regional atmospheric model. In those experiments, sea ice thermodynamics and dynamics as well as the exchanges with the ocean and atmosphere are strongly perturbed. This perturbation is achieved by modifying snow and ice thermal conductivities, the vertical mixing in the ocean top layers, the effect of freshwater uptake and release upon sea ice growth and melt, ice dynamics, and surface albedo. We find that the evolution of sea ice extent during the ice advance season is largely independent of the direct effect of the perturbation and appears thus mainly controlled by initial state in summer and subsequent insolation changes. In contrast, the melting rate varies strongly between the experiments during the retreat, in particular if the surface albedo or sea ice transport are modified, demonstrating a strong contribution of those elements to the evolution of ice coverage through spring and summer. As with the advance phase, the retreat is also influenced by conditions at the beginning of the melt season in September. Atmospheric feedbacks enhance the model winter ice extent response to any of the perturbed processes, and the enhancement is strongest when the albedo is modified. The response of sea ice volume and extent to changes in entrainment of subsurface warm waters to the ocean surface is also greatly amplified by the coupling with the atmosphere.

We coupled together high-resolution versions of the ocean–sea ice model NEMO and the ice sheet model BISICLES configured to the Totten Glacier area and ran a series of simulations over the recent past (1995–2014) and under warming conditions (2081–2100; SSP4-4.5) with NEMO in stand-alone mode and with the coupled model to assess the effects of the coupling. During the recent past, the ocean–ice sheet coupling has increased the time-averaged value of the basal melt rate in both the Totten and Moscow University ice shelf cavities by $6.7\%$ and $14.2\%$, respectively. The relationship between the changes in ice shelf thickness and ice shelf basal melt rate suggests that the effect of the coupling is not a linear response to the melt rate but rather a more complex response, driven partly by the dynamical component of the ice sheet model. The response of the ice sheet–ocean coupling due to the ocean warming is a 10% and 3% basal melt rate decrease in the Totten and Moscow University ice shelf cavities, respectively. This indicates that the ocean–ice sheet coupling under climate warming conditions dampens the basal melt rates. Our study highlights the importance of incorporating ocean–ice sheet coupling in climate simulations, even over short time periods.