Localized general vertical coordinates for quasi‐Eulerian ocean models: The Nordic overflows test‐case

A generalized methodology to deploy different types of vertical coordinate system in arbitrarily defined time-invariant local areas of quasi-Eulerian numerical ocean models is presented. After detailing its characteristics, we show how the general localization method can be used to improve the representation of the Nordic Seas overflows in the UK Met Office NEMO-based eddy-permitting global ocean configuration. Three z*-levels with partial steps configurations localizing different types of hybrid geopotential/terrain-following vertical coordinates in the proximity of the Greenland-Scotland ridge are implemented and compared against a control configuration. Experiments include a series of idealized and realistic numerical simulations where the skill of the models in computing pressure forces, reducing spurious diapycnal mixing and reproducing observed properties of the Nordic Seas overflows are assessed. Numerical results prove that the localization approach proposed here can be successfully used to embed terrain-following levels in a global geopotential levels-based configuration, provided that the localized vertical coordinate chosen is flexible enough to allow a smooth transition between the two. In addition, our experiments show that deploying localized terrain-following levels via the multi-envelope method allows the crucial reduction of spurious cross-isopycnal mixing when modeling bottom intensified buoyancy driven currents, significantly improving the realism of the Nordic Seas overflows simulations in comparison to the other configurations. Important hydrographic biases are found to similarly affect all the realistic experiments and a discussion on how their interaction with the type of localized vertical coordinate affects the realism of the simulated overflows is provided.