Resolving eddies by local mesh refinement

Abstract Nesting in large-scale ocean modeling is used for local refinement to resolve eddy dynamics that would not be accessible otherwise. Unstructured meshes offer this functionality too by adjusting their resolution according to some goal function. However, by locally refining the mesh one does not necessarily achieve the goal resolution, because the eddy dynamics, in particular the ability of eddies to release the available potential energy, also depend on the dynamics on the upstream coarse mesh. It is shown through a suite of experiments with a zonally re-entrant channel that baroclinic turbulence can be out from equilibrium in wide (compared to a typical eddy size) zones downstream into the refined area. This effect depends on whether or not the coarse part is eddy resolving, being much stronger if it is not. Biharmonic viscosity scaled with the cube of grid spacing is generally sufficient to control the smoothness of solutions on the variable mesh. However, noise in the vertical velocity field may be present at locations where the mesh is varied if momentum advection is implemented in the vector invariant form. Smoothness of vertical velocity is recovered if the flux form of momentum advection is used, suggesting that the noise originates from a variant of the Hollingsworth instability.

[1]  I. Held,et al.  Parameterizing subgrid-scale eddy effects using energetically consistent backscatter , 2014 .

[2]  Sergey Danilov,et al.  Two finite-volume unstructured mesh models for large-scale ocean modeling , 2012 .

[3]  Baylor Fox-Kemper,et al.  Can Large Eddy Simulation Techniques Improve Mesoscale Rich Ocean Models , 2013 .

[4]  R. Pacanowski,et al.  Parameterization of Vertical Mixing in Numerical Models of Tropical Oceans , 1981 .

[5]  C. Mertens,et al.  Circulation and transports in the Newfoundland Basin, western subpolar North Atlantic , 2014 .

[6]  Jonathan V. Durgadoo,et al.  Agulhas Leakage Predominantly Responds to the Southern Hemisphere Westerlies , 2013 .

[7]  H. Hellmer,et al.  Southern Ocean warming and increased ice shelf basal melting in the twenty-first and twenty-second centuries based on coupled ice-ocean finite-element modelling , 2013, Ocean Dynamics.

[8]  Robert Hallberg,et al.  Using a resolution function to regulate parameterizations of oceanic mesoscale eddy effects , 2013 .

[9]  E. Blayo,et al.  Two-way embedding algorithms: a review , 2008 .

[10]  Eugene Kazantsev.  Parameterizing subgrid scale eddy effects in a shallow water model , 2017 .

[11]  Daniel R. Lynch,et al.  Resolution issues in numerical models of oceanic and coastal circulation , 2007 .

[12]  Thierry Penduff,et al.  Mesoscale Eddies in the Labrador Sea and Their Contribution to Convection and Restratification , 2008 .

[13]  A. Hollingsworth,et al.  An internal symmetric computational instability , 1983 .

[14]  Philip W. Jones,et al.  A multi-resolution approach to global ocean modeling , 2013 .

[15]  Almut Gassmann,et al.  A global hexagonal C‐grid non‐hydrostatic dynamical core (ICON‐IAP) designed for energetic consistency , 2013 .

[16]  Qiang Wang,et al.  The Canadian Arctic Archipelago throughflow in a multiresolution global model: Model assessment and the driving mechanism of interannual variability , 2013 .