Description: Highlights: • Comparison of global NEMO and FESOM configurations with emphasis on the Agulhas system. • Both models simulate a reasonable and comparable large-scale circulation. • Both models have individual strengths and weaknesses to match the observations of the WBC system. • The numerical cost of FESOM is twice the one of NEMO. Abstract: Many questions in ocean and climate modelling require the combined use of high resolution, global coverage and multi-decadal integration length. For this combination, even modern resources limit the use of traditional structured-mesh grids. Here we compare two approaches: A high-resolution grid nested into a global model at coarser resolution (NEMO with AGRIF) and an unstructured-mesh grid (FESOM) which allows to variably enhance resolution where desired. The Agulhas system around South Africa is used as a testcase, providing an energetic interplay of a strong western boundary current and mesoscale dynamics. Its open setting into the horizontal and global overturning circulations also requires global coverage. Both model configurations simulate a reasonable large-scale circulation. Distribution and temporal variability of the wind-driven circulation are quite comparable due to the same atmospheric forcing. However, the overturning circulation differs, owing each model's ability to represent formation and spreading of deep water masses. In terms of regional, high-resolution dynamics, all elements of the Agulhas system are well represented. Owing to the strong nonlinearity in the system, Agulhas Current transports of both configurations and in comparison with observations differ in strength and temporal variability. Similar decadal trends in Agulhas Current transport and Agulhas leakage are linked to the trends in wind forcing. Although the number of 3D wet grid points used in FESOM is similar to that in the nested NEMO, FESOM uses about two times the number of CPUs to obtain the same model throughput (in terms of simulated model years per day). This is feasible due to the high scalability of the FESOM code.

Description: Mesoscale eddies and meanders have been shown to be one of the dominant sources of flow variability in the world's ocean. One example of an isolated eddy hotspot is the South-West Indian Ridge (SWIR). Several investigations have shown that the SWIR and the corresponding planetary potential vorticity field (f/H) exert a strong influence on the location and dynamics of the Antarctic Circumpolar Current (ACC), resulting in substantial fragmentation of the jets downstream of the ridge. The easterly extension of this eddy corridor appears to be restricted to the deep channel separating the Conrad Rise from the Del Cano and Crozet Plateau. However, while the fate of eddies formed at the SWIR has been widely investigated and the frontal character of this eastward extension is well known, the zone of diminishing variability that extends southwards to approximately 60°S remains poorly sampled. Using a combination of Argo, AVISO and NCEP/NCAR datasets, the character of this eddy corridor as a conduit for warm core eddies to move across the ACC into the Antarctic zone is investigated. In this study, we track a single warm-core eddy as it moves southwards from an original position of 31°E, 50°20'S to where it dissipates 10 months later in the Enderby Basin at 56°20'S. An Argo float entrained within the eddy confirms that its water masses are consistent with water found within the Antarctic Polar Frontal Zone north of the APF. Latent and sensible heat fluxes are on average 8W/m2 and 10W/m2 greater over the eddy than directly east of this feature. It is estimated that the eddy lost an average of 5W/m2 of latent heat and 5W/m2 of sensible heat over a 1-year period, an amount capable of melting approximately 0.92m of sea ice. In addition, using an eddy tracking algorithm a total of 28 eddies are identified propagating southwards, 25 of which are anti-cyclonic in rotation. Based on the new Argo float data, combined with AVISO and NCEP/NCAR datasets, these results suggest that the southward passage of warm-core eddies act as vehicles transporting heat, salt and biota southwards across the ACC and into the eastern boundary of the Weddell gyre.

Description: Ballast water treatment is required for vessels to prevent the introduction of potentially invasive neobiota. Some treatment methods use chemical disinfectants which produce a variety of halogenated compounds as disinfection by-products (DBPs). One of the most abundant DBP from oxidative ballast water treatment is bromoform (CHBr3) where we find an average concentration of 894±560nmolL-1 (226±142μgL-1) in the undiluted ballast water from measurements and literature. Bromoform is a relevant gas for atmospheric chemistry and ozone depletion, especially in the tropics where entrainment into the stratosphere is possible. The spread of DBPs in the tropics over months to years is assessed here for the first time. With Lagrangian trajectories based on the NEMO-ORCA12 model velocity field, we simulate DBP spread in the sea surface and try to quantify the oceanic bromoform concentration and emission to the atmosphere from ballast water discharge at major harbours in the tropical region of Southeast Asia. The exemplary simulations of two important regions, Singapore and the Pearl River Delta, reveal major transport pathways of the DBPs and the anthropogenic bromoform concentrations in the sea surface. Based on our simulations, we expect DBPs to spread into the open ocean, along the coast and also an advection with monsoon-driven currents into the North Pacific and Indian Ocean. Furthermore, anthropogenic bromoform concentrations and emissions are predicted to increase locally around large harbours. In the sea surface around Singapore we estimate an increase in bromoform concentration by 9% compared to recent measurement. In a moderate scenario where 70% of the ballast water is chemically treated bromoform emissions to the atmosphere can locally exceed 1000pmolm-2h-1 and double climatological emissions. In the Pearl River Delta all bromoform is directly outgassed which leads to an additional bromine (Br) input into the atmosphere of 495kmolBr (∼42tCHBr3) a-1. From Singapore ports the additional atmospheric Br input is calculated as 312kmolBr (∼26tCHBr3) a-1. We estimate the global anthropogenic Br input from ballast water into the atmosphere of up to 13Mmola-1. This is 0.1% global Br input from background bromoform emissions and thus probably not relevant for stratospheric ozone depletion.

Description: The Agulhas Current, the western boundary current of the South Indian Ocean, has been shown to play an important role in the connectivity between the Indian and Atlantic oceans. The greater Agulhas Current system is highly dominated by mesoscale dynamics. To investigate their influence on the regional and global circulations, a family of high-resolution ocean general circulation model configurations based on the NEMO code has been developed. Horizontal resolution refinement is achieved by embedding “nests” covering the South Atlantic and the western Indian oceans at 1/10∘ (INALT10) and 1/20∘ (INALT20) within global hosts with coarser resolutions. Nests and hosts are connected through two-way interaction, allowing the nests not only to receive boundary conditions from their respective host but also to feed back the impact of regional dynamics onto the global ocean. A double-nested configuration at 1/60∘ resolution (INALT60) has been developed to gain insights into submesoscale processes within the Agulhas Current system. Large-scale measures such as the Drake Passage transport and the strength of the Atlantic meridional overturning circulation are rather robust among the different configurations, indicating the important role of the hosts in providing a consistent embedment of the regionally refined grids into the global circulation. The dynamics of the Agulhas Current system strongly depend on the representation of mesoscale processes. Both the southward-flowing Agulhas Current and the northward-flowing Agulhas Undercurrent increase in strength with increasing resolution towards more realistic values, which suggests the importance of improving mesoscale dynamics as well as bathymetric slopes along this narrow western boundary current regime. The exploration of numerical choices such as lateral boundary conditions and details of the implementation of surface wind stress forcing demonstrates the range of solutions within any given configuration.