Impact of partial steps and momentum advection schemes in a global ocean circulation model at eddy-permitting resolution

Series of sensitivity tests were performed with a z-coordinate, global eddy-permitting (1/4°) ocean/sea-ice model (the ORCA-R025 model configuration developed for the DRAKKAR project) to carefully evaluate the impact of recent state-of-the-art numerical schemes on model solutions. The combination of an energy–enstrophy conserving (EEN) scheme for momentum advection with a partial step (PS) representation of the bottom topography yields significant improvements in the mean circulation. Well known biases in the representation of western boundary currents, such as in the Atlantic the detachment of the Gulf Stream, the path of the North Atlantic Current, the location of the Confluence, and the strength of the Zapiola Eddy in the south Atlantic, are partly corrected. Similar improvements are found in the Pacific, Indian, and Southern Oceans, and characteristics of the mean flow are generally much closer to observations. Comparisons with other state-of-the-art models show that the ORCA-R025 configuration generally performs better at similar resolution. In addition, the model solution is often comparable to solutions obtained at 1/6 or 1/10° resolution in some aspects concerning mean flow patterns and distribution of eddy kinetic energy. Although the reasons for these improvements are not analyzed in detail in this paper, evidence is shown that the combination of EEN with PS reduces numerical noise near the bottom, which is likely to affect current–topography interactions in a systematic way. We conclude that significant corrections of the mean biases presently seen in general circulation model solutions at eddy-permitting resolution can still be expected from the development of numerical methods, which represent an alternative to increasing resolution.

[1]  Ross J. Murray,et al.  Explicit Generation of Orthogonal Grids for Ocean Models , 1996 .

[2]  David G. Vaughan,et al.  BEDMAP: a new ice thickness and subglacial topographic model of Antarctica , 2001 .

[3]  P. Xie,et al.  Global Precipitation: A 17-Year Monthly Analysis Based on Gauge Observations, Satellite Estimates, and Numerical Model Outputs , 1997 .

[4]  W. Dewar.  Topography and barotropic transport control by bottom friction , 1998 .

[5]  Gurvan Madec,et al.  A global ocean mesh to overcome the North Pole singularity , 1996 .

[6]  P. M. Saunders,et al.  Bottom Currents Derived from a Shipborne ADCP on WOCE Cruise A11 in the South Atlantic , 1995 .

[7]  D. Stammer,et al.  On the mid‐depth circulation in the Labrador and Irminger Seas , 2001 .

[8]  E. Balopoulos,et al.  Objective analysis of temperature and salinity historical data set over the Mediterranean basin , 1998, IEEE Oceanic Engineering Society. OCEANS'98. Conference Proceedings (Cat. No.98CH36259).

[9]  K. Trenberth,et al.  Estimates of Meridional Atmosphere and Ocean Heat Transports , 2001 .

[10]  Audrey Estublier,et al.  Choice of an advection scheme for biogeochemical models , 2001 .

[11]  A. Arakawa,et al.  A Potential Enstrophy and Energy Conserving Scheme for the Shallow Water Equations , 1981 .

[12]  Gurvan Madec,et al.  Agulhas eddy fluxes in a 1/6° Atlantic model , 2002 .

[13]  H. Bryden,et al.  Thermohaline circulation at three key sections in the North Atlantic over 1985–2002 , 2005 .

[14]  C. Provost,et al.  A numerical study of quasi-geostrophic flow over isolated topography , 1985, Journal of Fluid Mechanics.

[15]  Gurvan Madec,et al.  Salt conservation, free surface, and varying levels: A new formulation for ocean general circulation models , 2000 .

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

[17]  Y. Masumoto.  A fifty-year eddy-resolving simulation of the world ocean : Preliminary outcomes of OFES (OGCM for the Earth Simulator) , 2004 .

[18]  J. Verron,et al.  How Topographic Smoothing Contributes to Differences between the Eddy Flows Simulated by Sigma- and Geopotential-Coordinate Models , 2002 .

[19]  B. Bourlès,et al.  Diagnosing the Annual Cycle of the Equatorial Undercurrent in the Atlantic Ocean from a General Circulation Model , 2006 .

[20]  R. Reynolds,et al.  The NCEP/NCAR 40-Year Reanalysis Project , 1996, Renewable Energy.

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

[22]  Bruno Blanke,et al.  Variability of the Tropical Atlantic Ocean Simulated by a General Circulation Model with Two Different Mixed-Layer Physics , 1993 .

[23]  A. Biastoch,et al.  The Role of Mesoscale Eddies in the Source Regions of the Agulhas Current , 1999 .

[24]  Dirk Olbers,et al.  The Antarctic Circumpolar Current System , 1999 .

[25]  David J. Webb,et al.  The first main run of the OCCAM global ocean model , 1997 .

[26]  Walter H. F. Smith,et al.  Global Sea Floor Topography from Satellite Altimetry and Ship Depth Soundings , 1997 .

[27]  C. Wunsch,et al.  Large-Scale Ocean Heat and Freshwater Transports during the World Ocean Circulation Experiment , 2003 .

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

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

[30]  R. Davis.  Intermediate-Depth Circulation of the Indian and South Pacific Oceans Measured by Autonomous Floats , 2005 .

[31]  W. Dewar,et al.  On the dynamics of the Zapiola Anticyclone , 1999 .

[32]  Valérie Dulière,et al.  On the representation of high latitude processes in the ORCA-LIM global coupled sea ice–ocean model , 2005 .

[33]  John Woodward,et al.  New grid of Arctic bathymetry aids scientists and mapmakers , 2000 .

[34]  P. Myers.  SPOM: A regional model of the sub‐polar north Atlantic , 2002 .

[35]  P. Killworth,et al.  The distribution of kinetic energy in the Southern Ocean: a comparison between observations and an eddy resolving general circulation model , 1992 .

[36]  A. Coward,et al.  Eddy mass transport for the Southern Ocean in an eddy-permitting global ocean model , 2003 .

[37]  M. Rhein,et al.  Modeling CFC inventories and formation rates of Labrador Sea Water , 2003 .

[38]  W. Large,et al.  A Global Ocean Wind Stress Climatology Based on ECMWF Analyses , 1989 .

[39]  F. Bryan,et al.  Large-scale transport processes in high-resolution circulation models , 1996 .

[40]  R. Sadourny.  The Dynamics of Finite-Difference Models of the Shallow-Water Equations , 1975 .

[41]  J. McWilliams,et al.  Dynamically balanced absolute sea level of the global ocean derived from near‐surface velocity observations , 2003 .

[42]  Andrew C. Coward,et al.  The OCCAM 66 level model: model description, physics, initial conditions and external forcing , 2005 .

[43]  Paul J. Martin,et al.  Validation of interannual simulations from the 1/8° global Navy Coastal Ocean Model (NCOM) , 2006 .

[44]  C. Provost,et al.  Influence of bottom topography roughness on the jet and inertial recirculation of a mid-latitude gyre , 1993 .

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

[46]  W. R. Holland,et al.  On the effects of a midocean ridge on the general circulation: numerical simulations with an eddy-resolved ocean model , 1987 .

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

[48]  M. Maltrud,et al.  Numerical simulation of the North Atlantic Ocean at 1/10 degrees , 2000 .

[49]  R. C. Malone,et al.  Global eddy-resolving ocean simulations driven by 1985-1995 atmospheric winds , 1998 .

[50]  Michael Steele,et al.  PHC: A Global Ocean Hydrography with a High-Quality Arctic Ocean , 2001 .

[51]  Gilles Reverdin,et al.  Global high-resolution mapping of ocean circulation from TOPEX/Poseidon and ERS-1 and -2 , 2000 .

[52]  M. Maqueda,et al.  Sensitivity of a global sea ice model to the treatment of ice thermodynamics and dynamics , 1997 .

[53]  G. Madec,et al.  An eddy-permitting model of the Atlantic circulation: evaluating open boundary conditions , 2001 .

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

[55]  R. Davis,et al.  The intermediate depth circulation of the western South Atlantic , 1999 .

[56]  W. Dewar,et al.  Mode waters and subduction rates in a high-resolution South Atlantic simulation , 1999 .

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

[58]  B. Cuevas,et al.  Dynamical budgets of the Antarctic Circumpolar Current using ocean general‐circulation models , 2005 .

[59]  R. Döscher,et al.  A Method for Improved Representation of Dense Water Spreading over Topography in Geopotential-Coordinate Models , 1997 .

[60]  Thierry Penduff,et al.  On the use of current meter data to assess the realism of ocean model simulations , 2005 .