Do we need non-Boussinesq effects in an ocean general circulation model for climate simulations?

The Boussinesq approximation is commonly made in ocean general circulation models (OGCMs). As a consequence, the model ocean is incompressible and conserves volume, but not mass. It has been argued that these consequence introduce errors at the noise level of coarse OGCMs, but that non-Boussinesq modeling is preferable simply for tidiness. Here, we use the height-pressure coordinate isomorphism implemented in the MITgcm to construct a non-Boussinesq OGCM and revisit the differences between Boussinesq and non-Boussinesq models at a resolution comparable to IPCC climate models. Subtleties such as the choice of a proper equation of state that includes the effect of pressure on heat capacity, but also the use of mass as a convenient alternative to pressure coordinates are discussed

Sea ice dynamics solvers in the MITgcm

Most dynamic sea ice models for climate type simulations are based on the viscous-plastic (VP) rheology. New rheologies such as the Maxwell-Elasto-Brittle (MEB) rheology are usually compared against traditional VP-schemes, but the new schemes also require revisiting the validity of VP-schemes. So far, comparisons between different schemes are confounded by factors unrelated to rheology, such as grid resolution, advection schemes, forcing by atmosphere and ocean, and last but not least, by differences in numerical details of different model codes. The sea ice component of the Massachusetts Institute of Technology general circulation model (MITgcm) offers an easy-to-use testbed for comparing different sea ice rheologies and implementation that avoids any confounders because all solvers share the same code and configuration environment. VP-rheologies with different flavors of Picard (or fixed point iterative) solvers, Newton methods, and different variants of the Elastic-Viscous-Plastic solver have been evaluated in this framework. With this framework, a new implementations such as an MEB solver may be compared to these traditional solvers in idealized geometries and in realistic Arctic configurations

Oceanic control of multidecadal variability in an idealized coupled GCM

Idealized ocean models are known to develop intrinsic multidecadal oscillations of the meridional overturning circulation (MOC). Here we explore the role of ocean–atmosphere interactions on this low-frequency variability. We use a coupled ocean–atmosphere model set up in a flat-bottom aquaplanet geometry with two meridional boundaries. The model is run at three different horizontal resolutions (4°, 2° and 1°) in both the ocean and atmosphere. At all resolutions, the MOC exhibits spontaneous variability on multidecadal timescales in the range 30–40 years, associated with the propagation of large-scale baroclinic Rossby waves across the Atlantic-like basin. The unstable region of growth of these waves through the long wave limit of baroclinic instability shifts from the eastern boundary at coarse resolution to the western boundary at higher resolution. Increasing the horizontal resolution enhances both intrinsic atmospheric variability and ocean–atmosphere interactions. In particular, the simulated atmospheric annular mode becomes significantly correlated to the MOC variability at 1° resolution. An ocean-only simulation conducted for this specific case underscores the disruptive but not essential influence of air–sea interactions on the low-frequency variability. This study demonstrates that an atmospheric annular mode leading MOC changes by about 2 years (as found at 1° resolution) does not imply that the low-frequency variability originates from air–sea interactions

Comparison of MERRA-2 and ECCO-V4 Ocean Surface Heat Fluxes: Consequences of Different Forcing Feedbacks on Ocean Circulation and Implications for Climate Data Assimilation

Forcing ocean models with reanalysis data is a common practice in ocean modeling. As part of this practice, prescribed atmospheric state variables and interactive ocean SST (Sea Surface Temperature) are used to calculate fluxes between the ocean and the atmosphere. When forcing an ocean model with reanalysis fields, errors in the reanalysis data, errors in the ocean model and errors in the forcing formulation will generate a different solution compared to other ocean reanalysis solutions (which also have their own errors). As a first step towards a consistent coupled ocean-atmosphere reanalysis, we compare surface heat fluxes from a state-of-the-art atmospheric reanalysis, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), to heat fluxes from a state-of-the-art oceanic reanalysis, the Estimating the Circulation and Climate of the Ocean Version 4, Release 2 (ECCO-v4). Then, we investigate the errors associated with the MITgcm (Massachusetts Institute of Technology general circulation model) ocean model in its ECCO-v4 ocean reanalysis configuration (1992-2011) when it is forced with MERRA- 2 atmospheric reanalysis fields instead of with the ECCO-v4 adjoint optimized ERA-interim state variables. This is done by forcing ECCO-v4 ocean with and without feedbacks from MERRA-2 related to turbulent fluxes of heat and moisture and the outgoing long wave radiation. In addition, we introduce an intermediate forcing method that includes only the feedback from the interactive outgoing long wave radiation. The resulting ocean circulation is compared with ECCO-v4 reanalysis and in-situ observations. We show that, without feedbacks, imbalances in the energy and the hydrological cycles of MERRA-2 (which are directly related to the fact it was created without interactive ocean) result in considerable SST drifts and a large reduction in sea level. The bulk formulae and interactive outgoing long wave radiation, although providing air-sea feedbacks and reducing model-data misfit, strongly relax the ocean to observed SST and may result in unwanted features such as large change in the water budget. These features have implications in a desired forcing recipe to be used. The results strongly and unambiguously argue for next generation data assimilation climate studies to involve fully coupled systems

Representing grounding line migration in synchronous coupling between a marine ice sheet model and a z-coordinate ocean model

Synchronous coupling is developed between an ice sheet model and a z-coordinate ocean model (the MITgcm). A previously-developed scheme to allow continuous vertical movement of the ice-ocean interface of a floating ice shelf (“vertical coupling”) is built upon to allow continuous movement of the grounding line, or point of floatation of the ice sheet (“horizontal coupling”). Horizontal coupling is implemented through the maintenance of a thin layer of ocean ( ∼ 1 m) under grounded ice, which is inflated into the real ocean as the ice ungrounds. This is accomplished through a modification of the ocean model’s nonlinear free surface evolution in a manner akin to a hydrological model in the presence of steep bathymetry. The coupled model is applied to a number of idealized geometries and shown to successfully represent ocean-forced marine ice sheet retreat while maintaining a continuous ocean circulation

Dynamic adjustment of the ocean circulation to self-attraction and loading effects

The oceanic response to surface loading, such as that related to atmospheric pressure, freshwater exchange, and changes in the gravity field, is essential to our understanding of sea level variability. In particular, so-called self-attraction and loading (SAL) effects caused by the redistribution of mass within the land–atmosphere–ocean system can have a measurable impact on sea level. In this study, the nature of SAL-induced variability in sea level is examined in terms of its equilibrium (static) and nonequilibrium (dynamic) components, using a general circulation model that implicitly includes the physics of SAL. The additional SAL forcing is derived by decomposing ocean mass anomalies into spherical harmonics and then applying Love numbers to infer associated crustal displacements and gravitational shifts. This implementation of SAL physics incurs only a relatively small computational cost. Effects of SAL on sea level amount to about 10% of the applied surface loading on average but depend strongly on location. The dynamic component exhibits large-scale basinwide patterns, with considerable contributions from subweekly time scales. Departures from equilibrium decrease toward longer time scales but are not totally negligible in many places. Ocean modeling studies should benefit from using a dynamical implementation of SAL as used here

ECCO Version 4: Second Release

The purpose of this note is twofold: (1) document the second release of ECCO version 4 state estimates (ECCO v4-r2); (2) provide a citable identifier to distinguish it from the previous release (ECCO v4-r1; Forget et al 2015).Major support for this work was provided through NASA’s Physical Oceanography Program. The bulk of the calculations was performed on the NASA Advanced Supercomputing (NAS) division’s Pleiades supercomputer at NASA/ARC

CATKE: a turbulent-kinetic-energy-based parameterization for ocean microturbulence with dynamic convective adjustment

We describe CATKE, a parameterization for ocean microturbulence with scales between 1 and 100 meters. CATKE is a one-equation model that predicts diffusive turbulent vertical fluxes a prognostic turbulent kinetic energy (TKE) and a diagnostic mixing length that features a dynamic model for convective adjustment (CA). With its convective mixing length, CATKE predicts not just the depth range where microturbulence acts but also the timescale over which mixing occurs, an important aspect of turbulent convection not captured by convective adjustment schemes. As a result, CATKE can describe the competition between convection and other processes such as baroclinic restractification or biogeochemical production-destruction. We estimate CATKE's free parameters with a posteriori calibration to eighteen large eddy simulations of the ocean surface boundary layer, and validate CATKE against twelve additional large eddy simulations with stronger and weaker forcing than used during calibration. We find that a CATKE-parameterized single column model accurately predicts the depth structure of buoyancy and momentum at vertical resolutions between 2 and 16 meters and with time steps of 10-20 minutes. We propose directions for future model development, and future efforts to recalibrate CATKE's parameters against more comprehensive and realistic datasets.Comment: submitted to J. Adv. Model. Earth Sy., 24 pages, 8 figure

Air-Sea Interactions in a High-Resolution Ocean-Atmosphere Simulation

During the past few years the Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have produced, respectively, global atmosphere-only and ocean-only simulations with km-scale grid spacing. These simulations have proved invaluable for process studies and for the development of satellite and in-situ sampling strategies. Nevertheless, a key limitation of these "nature" simulations is the lack of interaction between the ocean and the atmosphere, which limits their usefulness for studying air-sea interactions and for designing observing missions to study these interactions. We present here results from a coupled GEOS-MIT "nature run" simulation, wherein we have coupled a cubed-sphere-720 (~ 1/8) configuration of the GEOS atmosphere to a lat-lon-cap-1080 (~ 1/12) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere simulation to equivalent atmosphere-only and ocean-only simulations. A particular focus of the comparisons is the coupled versus uncoupled differences in interactions between Sea Surface Temperature (SST) and ocean surface wind. We discuss, in particular, a several-day mode of temporal variability in the SST-wind cycle and how it is represented in the different model simulations and in observationally-based products. A mechanism for the cycle, which is driven by SST-wind feedback, is proposed