Description: The fifth Workshop on Systematic Errors (WSE) in weather and climate models was hosted by Environment and Climate Change Canada (ECCC) on under the auspices of the Working Group on Numerical Experimentation (WGNE), jointly sponsored by the Commission of Atmospheric Sciences of the World Meteorological Organization (WMO) and the World Climate Research Programme (WCRP). This major event welcomed over 200 scientists from the weather and climate communities. The workshop primary goal was to increase the understanding of the nature and cause of systematic errors in numerical models across timescales. Out of 240 abstracts submitted to the workshop, 48 talks and 132 posters were presented.

Description: It is often unclear how to optimally choose horizontal resolution for the oceanic and atmospheric components of coupled climate models, which has implications for their ability to make best use of available computational resources. Here we investigate the effect of using different combinations of horizontal resolutions in atmosphere and ocean on the simulated climate in a global coupled climate model (Alfred Wegener Institute Climate Model [AWI‐CM]). Particular attention is given to the Atlantic Meridional Overturning Circulation (AMOC). Four experiments with different combinations of relatively high and low resolutions in the ocean and atmosphere are conducted. We show that increases in atmospheric and oceanic resolution have clear impacts on the simulated AMOC, which are largely independent. Increased atmospheric resolution leads to a weaker AMOC. It also improves the simulated Gulf Stream separation; however, this is only the case if the ocean is locally eddy resolving and reacts to the improved atmosphere. We argue that our results can be explained by reduced mean winds caused by higher cyclone activity. Increased resolution of the ocean affects the AMOC in several ways, thereby locally increasing or reducing the AMOC. The finer topography (and reduced dissipation) in the vicinity of the Caribbean basin tends to locally increase the AMOC. However, there is a reduction in the AMOC around 45°N, which relates to the reduced mixed layer depth in the Labrador Sea in simulations with refined ocean and changes in the North Atlantic current pathway. Furthermore, the eddy‐induced changes in the Southern Ocean increase the strength of the deep cell.

Description: Models from phase 5 of the Coupled Model Inter-comparison Project (CMIP5) show substantial biases in the deep ocean that are larger than the level of natural variability and the response to enhanced greenhouse gas concentrations. Here, we analyze the influence of horizontal resolution in a hierarchy of five multi-resolution simulations with the AWI Climate Model (AWI-CM), the climate model used at the Al-fred Wegener Institute, Helmholtz Centre for Polar and Ma-rine Research, which employs a sea ice–ocean model com-ponent formulated on unstructured meshes. The ocean grid sizes considered range from a nominal resolution of ∼ 1◦ (CMIP5 type) up to locally eddy resolving. We show that increasing ocean resolution locally to resolve ocean eddies leads to reductions in deep ocean biases, although these im-provements are not strictly monotonic for the five different ocean grids. A detailed diagnosis of the simulations allows to identify the origins of the biases. We find that two key re-gions at the surface are responsible for the development of the deep bias in the Atlantic Ocean: the northeastern North Atlantic and the region adjacent to the Strait of Gibraltar. Furthermore, the Southern Ocean density structure is equally improved with locally explicitly resolved eddies compared to parameterized eddies. Part of the bias reduction can be traced back towards improved surface biases over outcrop-ping regions, which are in contact with deeper ocean layers along isopycnal surfaces. Our prototype simulations provide guidance for the optimal choice of ocean grids for AWI-CM to be used in the final runs for phase 6 of CMIP (CMIP6) and for the related flagship simulations in the High Resolution Model Intercomparison Project (HighResMIP). Quite remarkably, retaining resolution only in areas of high eddy activity along with excellent scalability characteristics of the unstructured-mesh sea ice–ocean model enables us to per-form the multi-centennial climate simulations needed in a CMIP context at (locally) eddy-resolving resolution with a throughput of 5–6 simulated years per day.

Description: Mass loss around the Antarctic Ice Sheet is driven by basal melting and iceberg calving,which constitute the two dominant paths of freshwater flux into the Southern Ocean. Although of similarmagnitude, icebergs play an important and still not fully understood role in the balance of heat andfreshwater around Antarctica. This lack of understanding is partly due to operational difficulties inlarge-scale monitoring in polar regions, despite observational and remote sensing efforts. In this study, anovel machine learning approach, augmented by visual inspection, was applied to three SyntheticAperture Radar (SAR) mosaics of the whole Antarctic continent and its adjacent coastal zone. Althoughoriginally intended for a mapping of the Antarctic continent, the SAR mosaics allow us to document theevolution and distribution of the size (and mass) of icebergs in the pan-Antarctic near-coastal zone for theyears 1997, 2000, and 2008. Our novel algorithm identified 7,649 icebergs in 1997, 13,712 icebergs in 2000,and 7,246 icebergs in 2008 with surface areas between 0.1 and 4,567.82 km2and total masses of 4,641.53,6,862.81, and 5,263.69 Gt, respectively. Large regional variability was observed, although a zonal patterndistribution is present. This has implications for future climate modeling studies that try to estimate thefreshwater flux from melting icebergs, which demands a realistic representation of the interannuallyvarying near-coastal iceberg pattern to initialize the simulations.

Description: Digital drawings of the sensors and instruments that have been deployed on the sea ice as part of the "Distributed Network" around Polarstern - during leg 1 of the year-long MOSAiC expedition to the Arctic Ocean. The drawings can be freely used in presentations and for science communication, e.g. if a wide audience needs to be addressed. We would appreciate a proper citation (e.g. Krueger and Rackow, 2020). The drawings have all been done on board of the supporting Russian research vessel Akademik Fedorov in the Arctic Ocean.

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: In September 2019, the research icebreaker Polarstern started the largest multidisciplinary Arctic expedition to date, the MOSAiC (Multidisciplinary drifting Observatory for the Study of Arctic Climate) drift experiment. Being moored to an ice floe for a whole year, thus including the winter season, the declared goal of the expedition is to better understand and quantify relevant processes within the atmosphere–ice–ocean system that impact the sea ice mass and energy budget, ultimately leading to much improved climate models. Satellite observations, atmospheric reanalysis data, and readings from a nearby meteorological station indicate that the interplay of high ice export in late winter and exceptionally high air temperatures resulted in the longest ice-free summer period since reliable instrumental records began. We show, using a Lagrangian tracking tool and a thermodynamic sea ice model, that the MOSAiC floe carrying the Central Observatory (CO) formed in a polynya event north of the New Siberian Islands at the beginning of December 2018. The results further indicate that sea ice in the vicinity of the CO (〈40 km distance) was younger and 36 % thinner than the surrounding ice with potential consequences for ice dynamics and momentum and heat transfer between ocean and atmosphere. Sea ice surveys carried out on various reference floes in autumn 2019 verify this gradient in ice thickness, and sediments discovered in ice cores (so-called dirty sea ice) around the CO confirm contact with shallow waters in an early phase of growth, consistent with the tracking analysis. Since less and less ice from the Siberian shelves survives its first summer (Krumpen et al., 2019), the MOSAiC experiment provides the unique opportunity to study the role of sea ice as a transport medium for gases, macronutrients, iron, organic matter, sediments and pollutants from shelf areas to the central Arctic Ocean and beyond. Compared to data for the past 26 years, the sea ice encountered at the end of September 2019 can already be classified as exceptionally thin, and further predicted changes towards a seasonally ice-free ocean will likely cut off the long-range transport of ice-rafted materials by the Transpolar Drift in the future. A reduced long-range transport of sea ice would have strong implications for the redistribution of biogeochemical matter in the central Arctic Ocean, with consequences for the balance of climate-relevant trace gases, primary production and biodiversity in the Arctic Ocean.

Description: 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.

Description: Drifting icebergs represent a significant hazard for polar navigation and are able to impact the ocean environment around them. Freshwater flux and the associated cooling from melting icebergs can locally decrease salinity and temperature and thus affect ocean circulation, biological activity, sea ice, and –on larger spatial scales– the whole climate system. However, despite their potential impact, the large-scale operational monitoring of drifting icebergs in sea ice-covered regions is as of today typically restricted to giant icebergs, larger than 18.5 km in length. This is due to difficulties in accurately identifying and following the motion of much smaller features in the polar ocean from space. So far, tracking of smaller icebergs from satellite imagery thus has been limited to open-ocean regions not covered by sea ice. In this study, a novel automated iceberg tracking method, based on a machine learning-approach for automatic iceberg detection, is presented. To demonstrate the applicability of the method, a case study was performed for the Weddell Sea region, Antarctica, using 1213 Advanced Synthetic Aperture Radar (ASAR) satellite images acquired between 2002 and 2011. Overall, a subset of 414 icebergs (3134 re-detections in total) with surface areas between 3.4 km² and 3612 km² were investigated with respect to their prevalent drift patterns, size variability, and average disintegration. The majority of the tracked icebergs drifted between 1.3 km and 2679.2 km westward around the Antarctic continent, following the Antarctic Coastal Current (ACoC) and the Weddell Gyre, at an average drift speed of 3.6 ± 7.4 km day⁻¹. The method also allowed us to estimate an average daily disintegration (i.e. iceberg area decrease) rate of ~0.13% (~37% year⁻¹) for all icebergs. Using the sum of all detected individual surface area reductions, we estimate a total iceberg mass decrease of ~683 Gt year⁻¹, which can be freshwater input and/or new ‘child’ icebergs calved from larger icebergs. The extension to an automated long-term tracking method for icebergs is challenging as the iceberg shape can vary significantly due to abrupt disintegration or calving of bergy bits. However, our machine learning approach extended by automatic shape-based tracking capabilities proved to be a reliable alternative for automatic detection and tracking of icebergs, even under the ambiguous SAR background signatures often found in the Southern Ocean. In particular, the method works in the challenging near-coastal environment where the presence of sea ice and coastal ocean dynamics such as surface waves usually pose major obstacles for other approaches.

Description: The dominant feature of large-scale mass transfer in the modern ocean is the Atlantic meridional overturning circulation (AMOC). The geometry and vigour of this circulation influences global climate on various timescales. Palaeoceanographic evidence suggests that during glacial periods of the past 1.5 million years the AMOC had markedly different features from today; in the Atlantic basin, deep waters of Southern Ocean origin increased in volume while above them the core of the North Atlantic Deep Water (NADW) shoaled. An absence of evidence on the origin of this phenomenon means that the sequence of events leading to global glacial conditions remains unclear. Here we present multi-proxy evidence showing that northward shifts in Antarctic iceberg melt in the Indian–Atlantic Southern Ocean (0–50°E) systematically preceded deep-water mass reorganizations by one to two thousand years during Pleistocene-era glaciations. With the aid of iceberg-trajectory model experiments, we demonstrate that such a shift in iceberg trajectories during glacial periods can result in a considerable redistribution of freshwater in the Southern Ocean. We suggest that this, in concert with increased sea-ice cover, enabled positive buoyancy anomalies to ‘escape’ into the upper limb of the AMOC, providing a teleconnection between surface Southern Ocean conditions and the formation of NADW. The magnitude and pacing of this mechanism evolved substantially across the mid-Pleistocene transition, and the coeval increase in magnitude of the ‘southern escape’ and deep circulation perturbations implicate this mechanism as a key feedback in the transition to the ‘100-kyr world’, in which glacial–interglacial cycles occur at roughly 100,000-year periods.