Attribution of horizontal and vertical contributions to spurious mixing in an Arbitrary Lagrangian-Eulerian ocean model

Abstract We examine the separate contributions to spurious mixing from horizontal and vertical processes in an ALE ocean model, MOM6, using reference potential energy (RPE). The RPE is a global diagnostic which changes only due to mixing between density classes. We extend this diagnostic to a sub-timestep timescale in order to individually separate contributions to spurious mixing through horizontal (tracer advection) and vertical (regridding/remapping) processes within the model. We both evaluate the overall spurious mixing in MOM6 against previously published output from other models (MOM5, MITGCM and MPAS-O), and investigate impacts on the components of spurious mixing in MOM6 across a suite of test cases: a lock exchange, internal wave propagation, and a baroclinically-unstable eddying channel. The split RPE diagnostic demonstrates that the spurious mixing in a lock exchange test case is dominated by horizontal tracer advection, due to the spatial variability in the velocity field. In contrast, the vertical component of spurious mixing dominates in an internal waves test case. MOM6 performs well in this test case owing to its quasi-Lagrangian implementation of ALE. Finally, the effects of model resolution are examined in a baroclinic eddies test case. In particular, the vertical component of spurious mixing dominates as horizontal resolution increases, an important consideration as global models evolve towards higher horizontal resolutions.

[1]  Knut Klingbeil,et al.  Quantification of spurious dissipation and mixing – Discrete variance decay in a Finite-Volume framework , 2014 .

[2]  Jeffery R. Scott,et al.  Southern Ocean warming delayed by circumpolar upwelling and equatorward transport , 2016 .

[3]  Alistair Adcroft,et al.  Rescaled height coordinates for accurate representation of free-surface flows in ocean circulation models , 2004 .

[4]  A. Hogg,et al.  Adjustment of the Meridional Overturning Circulation and Its Dependence on Depth of Mixing , 2016 .

[5]  Mehmet Ilicak,et al.  Spurious dianeutral mixing and the role of momentum closure , 2012 .

[6]  Robert Pinkel,et al.  Global Patterns of Diapycnal Mixing from Measurements of the Turbulent Dissipation Rate , 2014 .

[7]  R. Ferrari,et al.  Influence of Enhanced Abyssal Diapycnal Mixing on Stratification and the Ocean Overturning Circulation , 2015 .

[8]  Stephen M. Griffies,et al.  Spurious Diapycnal Mixing Associated with Advection in a z-Coordinate Ocean Model , 2000 .

[9]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[10]  Christian Tenaud,et al.  High order one-step monotonicity-preserving schemes for unsteady compressible flow calculations , 2004 .

[11]  Todd D. Ringler,et al.  Evaluation of the arbitrary Lagrangian–Eulerian vertical coordinate method in the MPAS-Ocean model , 2015 .

[12]  M. Prather.  Numerical advection by conservation of second-order moments. [for trace element spatial distribution and chemical interaction in atmosphere] , 1986 .

[13]  Eric A. D'Asaro,et al.  Available potential energy and mixing in density-stratified fluids , 1995, Journal of Fluid Mechanics.

[14]  A. Hogg,et al.  Effect of topographic barriers on the rates of available potential energy conversion of the oceans , 2014 .

[15]  H. T. Huynh.  Schemes and constraints for advection , 1997 .

[16]  C. W. Hirt,et al.  An Arbitrary Lagrangian-Eulerian Computing Method for All Flow Speeds , 1997 .

[17]  Andrew T. Wittenberg,et al.  Impacts on Ocean Heat from Transient Mesoscale Eddies in a Hierarchy of Climate Models , 2015 .

[18]  Gurvan Madec,et al.  z˜-Coordinate, an Arbitrary Lagrangian–Eulerian coordinate separating high and low frequency motions , 2011 .

[19]  G. Hughes,et al.  Estimating Lorenz’s Reference State in an Ocean with a Nonlinear Equation of State for Seawater , 2015 .

[20]  Richard A. Gentry,et al.  An Eulerian differencing method for unsteady compressible flow problems , 1966 .

[21]  Laurent White,et al.  A high-order finite volume remapping scheme for nonuniform grids: The piecewise quartic method (PQM) , 2008, J. Comput. Phys..

[22]  Mehmet Ilıcak.  Quantifying spatial distribution of spurious mixing in ocean models , 2016, Ocean modelling.

[23]  Rainer Bleck,et al.  An oceanic general circulation model framed in hybrid isopycnic-Cartesian coordinates , 2002 .

[24]  Stephen M. Griffies,et al.  Fundamentals of Ocean Climate Models , 2004 .

[25]  Hartmut Peters,et al.  Improving Oceanic Overflow Representation in Climate Models: The Gravity Current Entrainment Climate Process Team , 2009 .

[26]  Laurent White,et al.  High-order regridding-remapping schemes for continuous isopycnal and generalized coordinates in ocean models , 2009, J. Comput. Phys..

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

[28]  H. Burchard,et al.  Comparative quantification of physically and numerically induced mixing in ocean models , 2008 .

[29]  F. Gouillon.  Internal wave propagation and numerically induced diapycnal mixing in oceanic general circulation models , 2010 .

[30]  P. Gent,et al.  Isopycnal mixing in ocean circulation models , 1990 .

[31]  A. Kasahara.  Various Vertical Coordinate Systems Used for Numerical Weather Prediction , 1974 .

[32]  M. A. Morales Maqueda,et al.  Second-order moment advection scheme applied to Arctic Ocean simulation , 2006 .

[33]  R. Hallberg,et al.  Parameterization of eddy fluxes based on a mesoscale energy budget , 2015 .

[34]  Len G. Margolin,et al.  Second-order sign-preserving conservative interpolation (remapping) on general grids , 2003 .