Evaluation of the arbitrary Lagrangian–Eulerian vertical coordinate method in the MPAS-Ocean model

Abstract The vertical coordinate of the Model for Prediction Across Scales-Ocean (MPAS-Ocean) uses the Arbitrary Lagrangian–Eulerian (ALE) method, which offers a variety of configurations. When fully Eulerian, the vertical coordinate is fixed like a z-level ocean model; when fully Lagrangian there is no vertical transport through the interfaces so that the mesh moves with the fluid; additional options for vertical coordinates exist between these two extremes, including z-star, z-tilde, sigma, and isopycnal coordinates. Here we evaluate spurious diapycnal mixing in MPAS-Ocean in several idealized test cases as well as real-world domains with full bathymetry. Mixing data is compared to several other ocean models, including the Parallel Ocean Program (POP) z-level and z-star formulations. In three-dimensional domains, MPAS-Ocean has lower spurious mixing that other ocean models. A series of simulations show that this is likely due to MPAS-Ocean’s hexagon-type horizontal grid cells combined with a flux-corrected transport tracer advection scheme designed for these unstructured meshes. The frequency-filtered vertical coordinate of Leclair and Madec (2011) (also called z-tilde) has been implemented and analyzed in MPAS-Ocean. This addition allows low-frequency vertical transport to pass through the vertical interface in an Eulerian manner, while high-frequency vertical oscillations, such as internal gravity waves, are treated in a Lagrangian manner. Z-tilde leads to a substantial reduction in vertical transport across layer interfaces, and a reduction in spurious diapycnal mixing.

[1]  McCarthyEd,et al.  A Unified Approach , 2005 .

[2]  William Skamarock,et al.  Conservative Transport Schemes for Spherical Geodesic Grids: High-Order Reconstructions for Forward-in-Time Schemes , 2010 .

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

[4]  Hailong Wang,et al.  Evaluation of Scalar Advection Schemes in the Advanced Research WRF Model Using Large-Eddy Simulations of Aerosol–Cloud Interactions , 2009 .

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

[6]  M. Alford,et al.  Redistribution of energy available for ocean mixing by long-range propagation of internal waves , 2003, Nature.

[7]  Todd D. Ringler,et al.  A Multiscale Nonhydrostatic Atmospheric Model Using Centroidal Voronoi Tesselations and C-Grid Staggering , 2012 .

[8]  Darryl D. Holm,et al.  Implementation of the LANS-alpha turbulence model in a primitive equation ocean model , 2007 .

[9]  T. Benjamin.  Gravity currents and related phenomena , 1968, Journal of Fluid Mechanics.

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

[11]  Rainer Bleck,et al.  On the Use of Hybrid Vertical Coordinates in Numerical Weather Prediction Models , 1978 .

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

[13]  Darryl D. Holm,et al.  Implementation of the LANS-α turbulence model in a primitive equation ocean model , 2008, J. Comput. Phys..

[14]  William C. Skamarock,et al.  Conservative Transport Schemes for Spherical Geodesic Grids: High-Order Flux Operators for ODE-Based Time Integration , 2011 .

[15]  Y. Song,et al.  A General Pressure Gradient Formulation for Ocean Models. Part I: Scheme Design and Diagnostic Analysis , 1998 .

[16]  Todd D. Ringler,et al.  An eddy closure for potential vorticity , 2009 .

[17]  A. Huerta,et al.  Arbitrary Lagrangian–Eulerian Methods , 2004 .

[18]  Alexander F. Shchepetkin,et al.  Accurate Boussinesq oceanic modeling with a practical, “Stiffened” Equation of State , 2011 .

[19]  Sonya Legg,et al.  Comparison of entrainment in overflows simulated by z-coordinate, isopycnal and non-hydrostatic models☆ , 2006 .

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

[21]  S. Manabe,et al.  A Global Ocean-Atmosphere Climate Model. Part II. The Oceanic Circulation , 1975 .

[22]  R. C. Malone,et al.  Parallel ocean general circulation modeling , 1992 .

[23]  Stephen G. Yeager,et al.  Diurnal to decadal global forcing for ocean and sea-ice models: The data sets and flux climatologies , 2004 .

[24]  Jean-Marc Molines,et al.  Circulation characteristics in three eddy-permitting models of the North Atlantic , 2001 .

[25]  William C. Skamarock,et al.  Numerical representation of geostrophic modes on arbitrarily structured C-grids , 2009, J. Comput. Phys..

[26]  Frank O. Bryan,et al.  Coordinated Ocean-ice Reference Experiments (COREs) , 2009 .

[27]  Alexander F. Shchepetkin,et al.  The regional oceanic modeling system (ROMS): a split-explicit, free-surface, topography-following-coordinate oceanic model , 2005 .

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

[29]  James C. McWilliams,et al.  A method for computing horizontal pressure‐gradient force in an oceanic model with a nonaligned vertical coordinate , 2003 .

[30]  Todd D. Ringler,et al.  A multiresolution method for climate system modeling: application of spherical centroidal Voronoi tessellations , 2008 .

[31]  W. Large,et al.  Oceanic vertical mixing: a review and a model with a nonlocal boundary layer parameterization , 1994 .

[32]  Rüdiger Gerdes,et al.  Formulation of an ocean model for global climate simulations , 2005 .

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

[34]  Mariana Vertenstein,et al.  The Parallel Ocean Program (POP) reference manual: Ocean component of the Community Climate System Model (CCSM) , 2010 .

[35]  K. Koltermann,et al.  WOCE Global Hydrographic Climatology , 2004 .

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

[37]  C. W. Hirt.  An arbitrary Lagrangian-Eulerian computing technique , 1971 .

[38]  Stephen M. Griffies,et al.  Biharmonic Friction with a Smagorinsky-Like Viscosity for Use in Large-Scale Eddy-Permitting Ocean Models , 2000 .

[39]  William C. Skamarock,et al.  A unified approach to energy conservation and potential vorticity dynamics for arbitrarily-structured C-grids , 2010, J. Comput. Phys..

[40]  J. Marshall,et al.  The Role of Eddy Transfer in Setting the Stratification and Transport of a Circumpolar Current , 2002 .

[41]  M. Maltrud,et al.  An eddy resolving global 1/10° ocean simulation , 2005 .

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

[43]  A. Adcroft,et al.  Representation of Topography by Shaved Cells in a Height Coordinate Ocean Model , 1997 .

[44]  T. McDougall,et al.  Minimal Adjustment of Hydrographic Profiles to Achieve Static Stability , 1995 .

[45]  Alistair Adcroft,et al.  On methods for solving the oceanic equations of motion in generalized vertical coordinates , 2006 .

[46]  Anand Gnanadesikan,et al.  Transient Response in a Z-Level Ocean Model That Resolves Topography with Partial Cells , 1998 .

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

[48]  S. Zalesak.  Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[49]  M. Gerritsen,et al.  An unstructured-grid, finite-volume, nonhydrostatic, parallel coastal ocean simulator , 2006 .