Description: We present four melt climatology estimates based on a simulation of Antarctic iceberg drift and melting that includes small, medium-sized, and giant tabular icebergs with a realistic size distribution. For the first time, an iceberg model is initialized with a set of nearly 7000 observed iceberg positions and sizes around Antarctica. We simulate drift and lateral melt using iceberg-draft averaged ocean currents, temperature, and salinity. A new basal melting scheme, originally applied in ice shelf melting studies, uses in situ temperature, salinity, and relative velocities at an iceberg's bottom. The climatology estimates based on simulations of small (SMA), 'small-to-medium'-sized (MED12 & MED123), and small-to-giant icebergs (ALL) exhibit differential characteristics: successive inclusion of larger icebergs leads to a reduced seasonality of the iceberg meltwater flux and a shift of the mass input to the area north of 58°S, while less meltwater is released into the coastal areas. This highlights the necessity to account for larger and giant icebergs in order to obtain accurate melt climatologies. The four monthly melt climatologies [mm/day] are available as netCDF files with 1°x1° spatial resolution and can be used, e.g., for sensitivity studies with uncoupled sea ice-ocean models, or as spatio-temporal templates for the redistribution of land ice from the Antarctic ice sheet over the Southern Ocean in climate models.

Description: We present a simulation of Antarctic iceberg drift and melting that includes small, medium‐sized, and giant tabular icebergs with a realistic size distribution. For the first time, an iceberg model is initialized with a set of nearly 7000 observed iceberg positions and sizes around Antarctica. The study highlights the necessity to account for larger and giant icebergs in order to obtain accurate melt climatologies. We simulate drift and lateral melt using iceberg‐draft averaged ocean currents, temperature, and salinity. A new basal melting scheme, originally applied in ice shelf melting studies, uses in situ temperature, salinity, and relative velocities at an iceberg's bottom. Climatology estimates of Antarctic iceberg melting based on simulations of small (≤2.2 km), “small‐to‐medium‐sized" (≤10 km), and small‐to‐giant icebergs (including icebergs 〉10 km) exhibit differential characteristics: successive inclusion of larger icebergs leads to a reduced seasonality of the iceberg meltwater flux and a shift of the mass input to the area north of 58°S, while less meltwater is released into the coastal areas. This suggests that estimates of meltwater input solely based on the simulation of small icebergs introduce a systematic meridional bias; they underestimate the northward mass transport and are, thus, closer to the rather crude treatment of iceberg melting as coastal runoff in models without an interactive iceberg model. Future ocean simulations will benefit from the improved meridional distribution of iceberg melt, especially in climate change scenarios where the impact of iceberg melt is likely to increase due to increased calving from the Antarctic ice sheet.

Description: We developed a new version of the Alfred Wegener Institute Climate Model (AWI-CM3), which has higher skills in representing the observed climatology and better computational efficiency than its predecessors. Its ocean component FESOM2 (Finite-volumE Sea ice-Ocean Model) has the multi-resolution functionality typical of unstructured-mesh models while still featuring a scalability and efficiency similar to regular-grid models. The atmospheric component OpenIFS (CY43R3) enables the use of the latest developments in the numerical-weather-prediction community in climate sciences. In this paper we describe the coupling of the model components and evaluate the model performance on a variable-resolution (25-125 km) ocean mesh and a 61 km atmosphere grid, which serves as a reference and starting point for other ongoing research activities with AWI-CM3. This includes the exploration of high and variable resolution and the development of a full Earth system model as well as the creation of a new sea ice prediction system. At this early development stage and with the given coarse to medium resolutions, the model already features above-CMIP6-average skills (where CMIP6 denotes Coupled Model Intercomparison Project phase 6) in representing the climatology and competitive model throughput. Finally we identify remaining biases and suggest further improvements to be made to the model.

Description: Viscosity in the momentum equation is needed for numerical stability, as well as to arrest the direct cascade of enstrophy at grid scales. However, a viscous momentum closure tends to over-dissipate eddy kinetic energy. To return excessively dissipated energy to the system, the viscous closure is equipped with what is called dynamic kinetic energy backscatter. The amplitude of backscatter is based on the amount of unresolved kinetic energy (UKE). This energy is tracked through space and time via a prognostic equation. Our study proposes to add advection of UKE by the resolved flow to that equation to explicitly consider the effects of nonlocality on the subgrid energy budget. UKE can consequently be advected by the resolved flow before it is reinjected via backscatter. Furthermore, we suggest incorporating a stochastic element into the UKE equation to account for missing small-scale variability, which is not present in the purely deterministic approach. The implementations are tested on two intermediate complexity setups of the global ocean model FESOM2: an idealized channel setup and a double-gyre setup. The impacts of these additional terms are analyzed, highlighting increased eddy activity and improved flow characteristics when advection and carefully tuned, stochastic sources are incorporated into the UKE budget. Additionally, we provide diagnostics to gain further insights into the effects of scale separation between the viscous dissipation operator and the backscatter operator responsible for the energy injection. Oceanic swirls or "eddies" have a typical size of 10-100 km, which is close to the smallest scales that global ocean models commonly resolve. For physical and numerical reasons, these models require the addition of artificial terms that influence the flow near its smallest scales. Common approaches have the drawback of introducing systematic loss of kinetic energy contained in the eddies, which leads to errors that also affect the oceanic circulation on global scales. In our research, we compensate for this error by returning some of the missing energy back into the simulation, using a so-called kinetic energy backscatter scheme. In this work, we continue the development of an already existing and successful backscatter scheme, adding certain improvements to the way energy is budgeted and returned to the flow: we ensure that the local energy budget is attached to each fluid parcel as it is transported by the large-scale flow, and we also add a random forcing term that mimics unknown sources of such energy to bring its statistical properties closer to reality. We demonstrate that these modifications effectively improve the characteristics of the simulated flow. Extension of the subgrid energy equation of the kinetic energy backscatter parameterization by adding advection and a stochastic term Both additional terms improve several flow characteristics in two idealized test cases, a channel and a double-gyre Scale analysis reveals the necessity of sufficient scale separation between viscous energy dissipation and energy injection via backscatter. Key Points: - Extension of the subgrid energy equation of the kinetic energy backscatter parameterization by adding advection and a stochastic term - Both additional terms improve several flow characteristics in two idealized test cases, a channel and a double-gyre - Scale analysis reveals the necessity of sufficient scale separation between viscous energy dissipation and energy injection via backscatter

Description: The influence of a stochastic sea ice strength parametrization on the mean climate is investigated in a coupled atmosphere–sea ice–ocean model. The results are compared with an uncoupled simulation with a prescribed atmosphere. It is found that the stochastic sea ice parametrization causes an effective weakening of the sea ice. In the uncoupled model this leads to an Arctic sea ice volume increase of about 10–20% after an accumulation period of approximately 20–30 years. In the coupled model, no such increase is found. Rather, the stochastic perturbations lead to a spatial redistribution of the Arctic sea ice thickness field. A mechanism involving a slightly negative atmospheric feedback is proposed that can explain the different responses in the coupled and uncoupled system. Changes in integrated Antarctic sea ice quantities caused by the stochastic parametrization are generally small, as memory is lost during the melting season because of an almost complete loss of sea ice. However, stochastic sea ice perturbations affect regional sea ice characteristics in the Southern Hemisphere, both in the uncoupled and coupled model. Remote impacts of the stochastic sea ice parametrization on the mean climate of non-polar regions were found to be small.

Description: A new climate model has been developed that employs a multi-resolution dynamical core for the sea ice-ocean component. In principle, the multi-resolution approach allows one to use enhanced horizontal resolution in dynamically active regions while keeping a coarse-resolution setup otherwise. The coupled model consists of the atmospheric model ECHAM6 and the finite element sea ice-ocean model (FESOM). In this study only moderate refinement of the unstructured ocean grid is applied and the resolution varies from about 25 km in the northern North Atlantic and in the tropics to about 150 km in parts of the open ocean; the results serve as a benchmark upon which future versions that exploit the potential of variable resolution can be built. Details of the formulation of the model are given and its performance in simulating observed aspects of the mean climate is described. Overall, it is found that ECHAM6–FESOM realistically simulates many aspects of the observed climate. More specifically it is found that ECHAM6–FESOM performs at least as well as some of the most sophisticated climate models participating in the fifth phase of the Coupled Model Intercomparison Project. ECHAM6–FESOM shares substantial shortcomings with other climate models when it comes to simulating the North Atlantic circulation.

Description: State-of-the-art global ocean circulation models used in climate studies are only passing the edge of becoming “eddy-permitting” or barely eddy-resolving. Such models commonly suffer from overdissipation of mesoscale eddies by routinely used subgrid dissipation (viscosity) operators and a resulting depletion of energy in the large-scale structures which are crucial for draining available potential energy into kinetic energy. More broadly, subgrid momentum closures may lead to both overdissipation or pileup of eddy kinetic energy and enstrophy of the smallest resolvable scales. The aim of this chapter is twofold. First, it reviews the theory of two-dimensional and geostrophic turbulence. To a large part, this is textbook material with particular emphasis, however, on issues relevant to modeling the global ocean in the eddy- permitting regime. Second, we discuss several recent parameterizations of subgrid dynamics, including simplified backscatter schemes by Jansen and Held, stochastic superparameterizations by Grooms and Majda, and an empirical backscatter scheme by Mana and Zanna.

Description: We present a simulation of Antarctic iceberg drift and melting that includes small, medium-sized, and giant tabular icebergs with a realistic size distribution. For the first time, an iceberg model is initialized with a set of nearly 7000 observed iceberg positions and sizes around Antarctica. The study highlights the necessity to account for larger and giant icebergs in order to obtain accurate melt climatologies. We simulate drift and lateral melt using iceberg-draft averaged ocean currents, temperature, and salinity. A new basal melting scheme, originally applied in ice shelf melting studies, uses in situ temperature, salinity, and relative velocities at an iceberg's bottom. Climatology estimates of Antarctic iceberg melting based on simulations of small (≤ 2.2 km), 'small-to-medium'-sized (≤ 10 km), and small-to-giant icebergs (including icebergs 〉 10 km) exhibit differential characteristics: successive inclusion of larger icebergs leads to a reduced seasonality of the iceberg meltwater flux and a shift of the mass input to the area north of 58 °S, while less meltwater is released into the coastal areas. This suggests that estimates of meltwater input solely based on the simulation of small icebergs introduce a systematic meridional bias; they underestimate the northward mass transport and are, thus, closer to the rather crude treatment of iceberg melting as coastal runoff in models without an interactive iceberg model. Future ocean simulations will benefit from the improved meridional distribution of iceberg melt, especially in climate change scenarios where the impact of iceberg melt is likely to increase due to increased calving from the Antarctic ice sheet.

Description: This study introduces a new flow‐dependent distribution sampling (FDDS) scheme for air–sea coupling. The FDDS scheme is implemented in a climate model and used to improve the simulated mean and variability of atmospheric and oceanic surface fields and thus air–sea fluxes. Most coupled circulation models use higher resolutions in the sea ice and ocean compared to the atmospheric model component, thereby explicitly simulating the atmospheric subgrid‐scale at the interface. However, the commonly applied averaging of surface fields and air–sea fluxes tends to smooth fine‐scale structures, such as oceanic fronts. The stochastic FDDS scheme samples the resolved spatial ocean (and sea ice) subgrid distribution that is usually not visible to a coarser‐resolution atmospheric model. Randomly drawn nodal ocean values are passed to the corresponding atmospheric boxes for the calculation of surface fluxes, aiming to enhance surface flux variability. The resulting surface field perturbations of the FDDS scheme are based on resolved dynamics, displaying pronounced seasonality with realistic magnitude. The AWI Climate Model is used to test the scheme on interannual time‐scales. Our set‐up features a high ocean‐to‐atmosphere resolution ratio in the Tropics, with grid‐point ratios of about 60:1. Compared to the default deterministic averaging, changes are largest in the Tropics leading to an improved spatial distribution of precipitation with bias reductions of up to 50%. Enhanced sea‐surface temperature variability in boreal winter further improves the seasonal phase locking of temperature anomalies associated with the El Niño–Southern Oscillation. Mean 2m temperature, sea ice thickness and concentration react with a contrasting dipole pattern between hemispheres but a joint increase of monthly and interannual variability. This first approach to implement a flow‐dependent stochastic coupling scheme shows considerable benefits for simulations of global climate, and various extensions and modifications of the scheme are possible.

Description: Atmosphere-ocean coupling is of major importance for numerical climate models. Systematic errors in surface flux computations affect both components simultaneously. As an example, surface flux variations caused by (sub-)mesoscale oceanic features and their impact on the atmospheric circulation are generally underrepresented in current models. In the case of coupled climate models, atmospheric sub-grid scale information is in fact available at the surface, as part of it is typically resolved by the oceanic model component. We present a stochastic scheme that couples resolved spatial ocean (and ice) variability, previously not visible to the atmospheric model. Consequently, the resulting SST perturbations of the stochastic scheme are based on resolved dynamics, displaying a pronounced seasonality and realistic magnitude. We apply the stochastic coupling method in the AWI Climate Model, a new multi-resolution climate model with an ocean component (FESOM) supporting unstructured triangular grids. Our specific setup features a high ocean-to-atmosphere resolution ratio in the tropics, with grid point ratios of about 60:1 (less than 25km to about 200km). Compared to the default deterministic coupling, changes are largest in the tropics, leading to an improved distribution of convective precipitation, and reductions of prominent biases of up to 50%. The scheme does not rely on a simultaneous increase of resolution in the atmospheric component, saving computational resources. The results are also of interest to other modeling centres employing high ocean-to-atmosphere resolution ratios, as the coupling scheme could be implemented relatively easily as an additional coupling option, without a change to the model code. (Oral in AS5.6/BG4.14/CL5.09/OS1.14 - Recent Developments in Numerical Earth System Modelling)