Articles | Volume 13, issue 1
Model description paper
09 Jan 2020
Model description paper |  | 09 Jan 2020

MOMSO 1.0 – an eddying Southern Ocean model configuration with fairly equilibrated natural carbon

Heiner Dietze, Ulrike Löptien, and Julia Getzlaff

We present a new near-global coupled biogeochemical ocean-circulation model configuration. The configuration features a horizontal discretization with a grid spacing of less than 11 km in the Southern Ocean and gradually coarsens in meridional direction to more than 200 km at 64 N, where the model is bounded by a solid wall. The underlying code framework is the Geophysical Fluid Dynamics Laboratory (GFDL)'s Modular Ocean Model coupled to the Biogeochemistry with Light, Iron, Nutrients and Gases (BLING) ecosystem model of Galbraith et al. (2010). The configuration is unique in that it features both a relatively equilibrated oceanic carbon inventory and an eddying ocean circulation based on a realistic model geometry/bathymetry – a combination that has been precluded by prohibitive computational cost in the past. Results from a simulation with climatological forcing and a sensitivity experiment with increasing winds suggest that the configuration is sufficiently equilibrated to explore Southern Ocean carbon uptake dynamics on decadal timescales. The configuration is dubbed MOMSO, a Modular Ocean Model Southern Ocean configuration.

1 Introduction

The Southern Ocean, also known as the Antarctic Ocean, comprises the southernmost waters of the world oceans. Its southern boundary is set by the Antarctic continent. With regard to its northern boundary to the Atlantic, Pacific and Indian oceans, there is no consensus. One common definition is the location of the Subtropical Front (STF), which separates the relatively saline subtropical waters from the fresher Subantarctic waters. This definition is straightforward because it characterizes the transition from one oceanographic regime to another distinctly different one. However, because the position of the front varies with time, the definition adds complexity to certain analyses such as oceanic inventories. A pragmatic solution is to define a fixed latitude as the northern boundary. We choose 40 S in order to facilitate a comparison with, e.g., Lovenduski et al. (2013) and Dietze et al. (2017).

The atmospheric conditions in the Southern Ocean are characterized by frequent cyclonic storms that travel eastward around the Antarctic continent and high-pressure areas over the poles. This averages to a strong westerly wind belt between 40 and 60 S, and rather weak and irregular polar easterlies from 60 S on southwards. The dominant low-frequency mode of atmospheric variability is summarized by the Southern Annular Mode (Limpasuvan and Hartmann1999; Marshall2003) which essentially describes the strength and position of the westerly wind belt.

The strong westerly wind belt drives the Antarctic Circumpolar Current (ACC), the “mightiest current in the oceans” (Pickard and Emery1990), which is unobstructed by continents. The ACC circumnavigates the globe zonally and thus links the Atlantic, Pacific and Indian oceans. Northward Ekman transport (also driven by the westerlies) drives up- and downwelling (Marshall and Speer2012) and is part of a meridional overturning circulation. Among the major processes controlling the net overturning is the so-called “eddy compensation” which largely cancels changes in wind-driven overturning (see Marshall and Radko2003; Viebahn and Eden2010; Abernathey et al.2011, 2016; Hallberg and Gnanadesikan2006; Thompson et al.2014; Tamsitt et al.2017). Further complexity is added by sea ice which modulates air–sea buoyancy fluxes and thus exerts control on convection, a prerequisite for Antarctic Bottom Water (AABW) formation.

Being one of the few places worldwide where deep water (such as AABW) is formed, the Southern Ocean plays a major role in the global carbon budget. Convection events tap into the abyssal ocean and modulate the difference between atmospheric and oceanic CO2 concentrations which drive net air–sea fluxes. A comprehensive quantitative understanding has not been reached yet, but there is consensus that the variability in the extent to which deep-water masses in the Southern Ocean are isolated from the atmosphere is among the major drivers regulating atmospheric CO2 variability (e.g., Anderson et al.2009; van Heuven et al.2014; Ritter et al.2017). Consequently, the Southern Ocean has shifted into the limelight of climate research (DeVries et al.2017; Tamsitt et al.2017; Langlais et al.2017, and many more).

To date, the Southern Ocean accounts for almost half of the global oceanic CO2 uptake from the atmosphere (Takahashi et al.2012). But there is concern that anticipated climate change (e.g., via changes of atmospheric circulation and sea ice cover) may trigger substantial changes in the Southern Ocean carbon budget (e.g., Heinze et al.2015; Abernathey et al.2016) such that the current rate of uptake may well decline in decades to come. Indications for the existence of such potential triggers have been revealed by observation-based atmospheric reanalysis products which show an ongoing strengthening and a poleward shift of the southern westerly winds starting from the 1970s (Thompson and Solomon2002).

This observed trend is projected by climate scenarios to intensify (e.g., Simpkins and Karpechko2012), and it is straightforward to assume that the associated wind-driven circulation affects biogeochemical dynamics and, eventually, the oceanic carbon budget. A comprehensive understanding of the link between changing winds and oceanic upwelling of carbon-rich deep water (which, in turn, affects surface saturation and net air–sea CO2 exchange) has, however, not been achieved yet. To this end, the role of mesoscale ocean eddies is especially uncertain: the current generation of coarse-resolution (non-mesoscale-resolving) models suggests that a poleward shift and an intensification of the Southern Ocean westerlies results in a strengthening of the subpolar meridional overturning cell (e.g., Saenko et al.2005; Hall and Visbeck2002; Getzlaff et al.2016) and, consequently, in increased upwelling of deep water south of the circumpolar flow, which is rich in dissolved inorganic carbon (e.g., Zickfeld et al.2007; Lenton and Matear2007; Lovenduski et al.2008; Verdy et al.2007). The net effect here is that changes of the atmospheric circulation reduce the capability of the Southern Ocean to sequester carbon away from the atmosphere. Early on, however, there were indications (e.g., Böning et al.2008; Hallberg and Gnanadesikan2006; Hogg et al.2008; Screen et al.2009; Thompson and Solomon2002) that this model behavior is a spurious consequence of the underlying eddy parameterization in coarse-resolution models which cannot afford to resolve mesoscale dynamics explicitly. Supporting evidence came from Munday et al. (2014), who showed that mesoscale eddies indeed reduce the sensitivity of oceanic carbon sequestration towards changing wind stress in an idealized model.

To date, we know that very different state-of-the-art approaches to parameterize eddies yield surprisingly similar sensitivities of oceanic carbon inventories to changing winds (Dietze et al.2017). Further, it has been demonstrated that “the conversion of dense waters back to light waters is sensitive to the degree of parameterization of mesoscale eddies in models” (Spence et al.2010, 2009). The question, however, as to how a realistic (as opposed to one featuring an idealized geometry) high-resolution coupled biogeochemical ocean-circulation model that actually resolves eddies and carbon dynamics explicitly compares to a coarse-resolution model has not been answered yet. Among the reasons is the prohibitive computational cost that is associated with equilibrating simulated dissolved inorganic carbon concentrations specifically in the deep ocean.

In this model description paper, we present the realistic high-resolution Modular Ocean Model Southern Ocean model configuration (MOMSO 1.0). The configuration features a realistic topography and simulated levels of eddy kinetic energy which do not undercut observed values. The latter suggests that the configuration explicitly resolves a substantial part of mesoscale-related variability rather than relying on parameterizing their effect. MOMSO is designed to explore the sensitivity of the Southern Ocean carbon uptake to atmospheric changes on decadal scales. The configuration is rendered feasible by recent advances in compute hardware and by its similarity to a spun-up coarse-resolution model which delivered the initial conditions for the biogeochemical module. (Note that module refers to the algorithmic entity that calculates the local biogeochemical sources and sinks of the prognostic biogeochemical variables such as nutrients and carbon. The algorithmic entity comprises both the underlying partial differential equations and their approximated representation within the numerical solver.) More specifically, we will showcase that the “level of equilibration” of simulated dissolved inorganic carbon allows us to test the sensitivity of the Southern Ocean carbon budget to anticipated climate change patterns.

In summary, this model description paper aims to (1) describe a new eddying coupled ocean-circulation biogeochemical model configuration of the Southern Ocean and (2) outline research questions for which MOMSO may serve as the base of a tool bench.

The project is dubbed MOMSO, a configuration of the Geophysical Fluid Dynamics Laboratory (GFDL)'s Modular Ocean Model version 4p1 with enhanced resolution in the Southern Ocean. The naming is a homage to the underlying framework, the MOM4p1 release of NOAA's GFDL Modular Ocean Model (Griffies2009). The ocean-circulation model is coupled to a sea ice model and the Biogeochemistry with Light, Iron, Nutrients and Gases (BLING) model from Galbraith et al. (2010).

2 Model setup

This model description paper describes simulations with the Modular Ocean Model (MOM), version MOM4p1 coupled to the Sea Ice Simulator (SIS) (Griffies2009). The configuration is near global, bounded by Antarctica and 64 N. In the Southern Ocean, the horizontal resolution is higher than 11 km. From 40 S, the meridional grid resolution coarsens towards the north (Fig. 1). There are no open boundaries and there is no tidal forcing.

Figure 1MOMSO model domain and spatial (finite difference) discretization. Panel (a) shows the model bathymetry. Panels (b), (c) and (d) show the vertical, meridional and zonal resolutions, respectively. In total, there are 2400×482×55 grid boxes in zonal, meridional and vertical directions, respectively.

The biogeochemical model BLING (Galbraith et al.2010) is coupled online to the ocean–sea ice model. Atmospheric CO2 concentrations are prescribed to a pre-industrial level of 278 ppmv. The respective carbon inventories and fluxes are referred to as “natural carbon”.

2.1 Grid and bathymetry

The underlying bathymetry is ETOPO5 (see Data Announcement 88-MGG-02, Digital relief of the Surface of the Earth. NOAA, National Geophysical Data Center, Boulder, Colorado, 1988). Using a bilinear scheme, the bathymetry is interpolated onto an Arakawa B grid (Arakawa and Lamp1977) with 2400×482 tracer grid boxes in the horizontal. The ocean-circulation and sea ice models share the same horizontal grid.

The vertical discretization comprises a total of 55 levels. Figure 1 shows the nominal depth and thickness of each level. The model bathymetry is smoothed with a filter similar to the Shapiro filter (Shapiro1970). The filter weights are 0.25, 0.5 and 0.25. The filtering procedure can only decrease the bottom depth; i.e., essentially, it fills rough holes. The filter is applied three times consecutively because we found this to be a good compromise between unnecessary smoothing on the one hand and numerical instabilities introduced by overly steep topography on the other hand in other high-resolution model configurations (Dietze and Kriest2012; Dietze et al.2014). The resulting bathymetry contained lakes which we filled after visual inspection. In addition, we filled narrow inlets which had a width of less than three grid boxes. In total, MOMSO has 42 429 759 wet tracer grid boxes.

Figure 2Close-up of MOMSO bathymetry in the Drake Passage. The color denotes depths in meters. The red rectangle denotes the area where, from 3250 m down to the bottom, an additional Laplacian horizontal viscosity of 1000 m2 s−2 is applied for stability reasons.

We use the zstar coordinate in the vertical (Stacey et al.1995; Adcroft and Campin2004) which is essentially an extension of the non-linear free-surface method of Campin et al. (2004) to all model levels. The algorithm is well-established for its very accurate conservation of tracers. zstar (z) is calculated as a function of nominal depth (z), respective water depth (H) and the free sea surface height (η) which varies with time:

(1) z = H z - η H + η

(Eq. 6.6 in Griffies2009). The approach overcomes the problem with vanishing surface grid boxes which appears in generic z-level discretization when sea surface height variations are of similar magnitude to the thickness of the uppermost grid box.

2.2 Ocean circulation

MOM4p1 is a z coordinate, free-surface ocean general circulation model which discretizes the ocean's hydrostatic primitive equations on a fixed Eulerian grid. The vertical mixing of momentum and scalars is parameterized by the k-profile parameterization approach of Large et al. (1994) with the same parameters applied in eddy-permitting global configurations of Dietze and Kriest (2012), Dietze and Loeptien (2013), Dietze et al. (2014), and Liu et al. (2010). The relevant parameters are (1) a critical bulk Richardson number of 0.3 and (2) a constant vertical background diffusivity and viscosity of 10−5 m2 s−1. These background values apply also below the surface mixed layer throughout the water column. Both parameterizations of the non-local and the double diffusive (vertical) scalar tracer fluxes are applied.

We apply a state-dependent horizontal Smagorinsky viscosity scheme (Griffies and Hallberg2000; Smagorinski1963, 1993) to keep friction at the minimal level necessitated by numerical stability. We apply a respective coefficient that sets the scale of the Smagorinsky isotropic viscosity to 0.01. Note that this value is much smaller than the range between 2 and 4 recommended by Griffies and Hallberg (2000) for stability reasons. The rationale behind our choice is to minimize friction which has little physical justification and at the same time has been shown to degrade the performance of ocean-circulation models (e.g., Jochum et al.2008).

We use the PPM (piecewise parabolic method) advection scheme (Colella and Woodward1984) for active tracers because it has been shown to perform well in treating sharp gradients in small-scale circulation (e.g., Carpenter Jr. et al.1989). For the biogeochemical tracers, we use a flux-limited scheme following Sweby (1984) because we experienced problems in the past with other configurations where negative values occurred for biogeochemical tracers when using non-flux-limited schemes. The choice of two different advection schemes may have repercussions. For example, Lévy et al. (2001) have shown that the choice spuriously affects biogeochemistry in regions of sharp vertical macronutrient gradients such as in oligotrophic regions. We expect that this is not so much an issue in MOMSO because the Southern Ocean vertical nutrient gradients are smoother. Further, we argue that all advection schemes introduce spurious errors (either of diffusive or dispersive nature). These are directly related to the respective property gradient in flow direction. Since all properties (such as temperature, salinity and biogeochemical variables) feature differing gradients, they are affected to a differing degree by the spurious behavior of the advection scheme. Hence, even when using the same advection scheme for all properties, the actual transport of properties is still inconsistent in that they are all affected by differing degrees by spurious errors. This inconsistency can be larger than the inconsistency introduced by switching from one advection scheme to another. Note in this context that we do not apply an explicit horizontal background diffusivity other than the contribution that is implicit to the advection scheme.

Several decades into the spin-up, our configuration became unstable in coarsely resolved places where strong currents met rough topography. This may well have been the result of our very low choice of the Smagorinsky isotropic viscosity (as reviewer 1 suspected). Our fix to the instability problem was to set an additional horizontal isotropic Laplacian viscosity of 600 m2 s−2 from 10 S to 50 N, of 1200 m2 s−2 above 50 N and 1800 m2 s−2 above 60 N outside the Southern Ocean towards the northern boundary of the model domain. This kept the respective oscillations in check. In addition, we added Laplacian viscosity at the exit of the Drake Passage (Fig. 2).

2.3 Sea ice

The ocean component is coupled to a dynamical sea ice module, the GFDL SIS. SIS uses elastic–viscous–plastic rheology adapted from Hunke and Dukowicz (1997). In the standard version, the simulated sea ice impacts sea surface height. This led to a vicious cycle at some places where sea ice attracted ever more sea ice, resulting in unrealistic anomalies in sea surface height and finally in a numerical blow-up of the simulation. Because of the associated computational cost, we did not investigate thoroughly into this but – instead – switched to levitating sea ice by deleting line 1353 in the file ice_model.F90, which reads Ice%p_surf(i,j) = Ice%p_surf(i,j) + grav*x(i,j). Note that this is now also used, e.g., in Hordoir et al. (2019), which is based on the Nucleus for European Modelling of the Ocean (NEMO) framework (Robinson Hordoir, personal communication, 2019).

2.4 Biogeochemistry

In our setup, the ocean–sea ice module is coupled to the BLING ecosystem model of Galbraith et al. (2010). BLING is a prognostic model that, in the basic version, explicitly resolves only four biogeochemical tracers: dissolved inorganic phosphorous, dissolved organic phosphorous, dissolved iron and dissolved oxygen. Here, we use BLING in conjunction with a carbon module that explicitly resolves dissolved inorganic carbon and alkalinity as described, e.g., in Bernadello et al. (2014).

The design idea behind the “reduced-tracer” model BLING is, on the one hand, to minimize the number of prognostic tracers that are actually advected by the ocean-circulation model and, on the other hand, to explicitly resolve the most influential environmental processes that control the net biotic carbon uptake (macro- and micronutrients phosphorous and iron; oxygen which controls heterotrophic respiration). The rationale behind minimizing the number of prognostic tracers is that even simple models (resolving only macronutrients, two phytoplankton species, zooplankton and detritus) are hard to constrain with available observations already (cf. Löptien and Dietze2015, 2017, 2019) – a difficulty that increases with the number of tracers. As a favorable side effect, the reduced number of prognostic tracers reduces the computational cost.

2.5 Initial conditions and spin-up procedure

The circulation model starts from rest (i.e., initial velocities are nil) with initial values for temperature and salinity taken from the World Ocean Atlas (WOA) 2009 (Locarnini et al.2010; Antonov et al.2009, respectively,). After 20 years of physics-only spin-up, the biogeochemical model is hooked on. In contrast to the physics, we do not use an observational product to initialize the prognostic variables of the biogeochemical model – also because the observations of dissolved iron are so sparse that Orr et al. (2017) suggest not to use them for model initialization. Instead, the initial conditions for the biogeochemical tracers are interpolated from the fully spun-up coarse-resolution configuration used by Dietze et al. (2017) (their “FMCD” simulation). This coarse-resolution configuration is, apart from the spatial discretization (and related parameters in the physical parameterizations of unresolved processes), identical to the MOMSO configuration. A beneficial aspect of this spin-up is that it accelerates the equilibration of the carbon dynamics substantially compared to initializing with an observational products. After a subsequent 60-year spin-up with online biogeochemistry, the model allows already (as we will put forward in Sect. 4) for an investigation of circulation-driven decadal changes of the Southern Ocean carbon budget.

2.6 Boundary conditions and sponges

The boundaries towards the Arctic (i.e., the northern end of the model domain shown in Fig. 1) are closed (i.e., they are represented by solid and flat walls). Temperature and salinity are restored to climatological annual mean estimates (Locarnini et al.2010; Antonov et al.2009) in so-called “sponge zones” located in the coarse-resolution domain. The sponge zones along with restoring timescales are shown in Fig. 3. The purpose of these sponges is to ensure realistic deep-water characteristics even though Northern Hemisphere deep-water formation processes are handicapped by the combination of coarse resolution with the absence of eddy parameterizations.

Figure 3MOMSO domains where temperature and salinity are restored to observed climatological values (Locarnini et al.2010; Antonov et al.2009, respectively) throughout the water column. The red and blue patches denote e-folding restoring timescales of 10 and 5 years, respectively.

At the air–sea boundary, we apply climatological atmospheric conditions taken from the Corrected Normal Year Forcing (COREv2; Large and Yeager2004). In addition, we apply a surface salinity restoring to climatological values (Antonov et al.2009) with a timescale of half a year throughout the model domain.

Atmospheric CO2 concentrations are prescribed to a pre-industrial level of 278 ppmv. Thus, the simulated oceanic carbon is also referred to as “natural carbon”. Biogeochemical air–sea fluxes (of iron) are identical to the ones applied in Galbraith et al. (2010) and Dietze et al. (2017).

This paper refers to output from two simulations dubbed REF and WIND. Both simulations share the same 80-year spin-up described in Sect. 2.5. WIND branches off from REF during the nominal year 1980 and is exposed to ever-increasing wind speeds south of 40 S. The increase is linear at a rate of 14 % in 50 years, consistent with results from a reanalysis of the period 1958 to 2007 (Lovenduski et al.2013). The rationale behind presenting results from the sensitivity experiment WIND in this model description paper is that they serve as a reference point against which the remaining model drift in REF (which is there since we cannot afford thousands of years of spin-up) can be compared. Based on this comparison, Sect. 4 sketches – very briefly – research questions that may be tackled with MOMSO.

IO-related hardware problems caused data loss such that REF covers the period from 1980 to 2024, while WIND covers 1980 to 2022 only.

Please note that a comprehensive analysis of the sensitivity experiment WIND is beyond the scope of this model description paper. Our major aim here is to show that REF is sufficiently equilibrated so that trends affected by decadal-scale changes in winds can clearly be distinguished against the background trend that is still persistent.

3 Results

In the following, we evaluate our model (simulation REF) by comparing our climatological results from the nominal years 1980–2024 to observational data. One problem is the tradeoff between data coverage and the length of the period the data are representative of. For any given year, the data coverage is insufficient to compile a comprehensive three-dimensional gridded data product. Binning the observational data of several years or even decades into one product closes spatial data gaps, but then this blurs the referencing to an ever-changing (anthropogenically driven) system state. This problem is especially pronounced in the Southern Ocean where in situ data acquisition is complicated by hostile environmental conditions.

The climatological atmospheric boundary conditions which drive our ocean model are representative of the period 1958–2000. Climatological data products are typically biased in that they contain more recent observations, being the result of recent technological advances (such as the development of autonomous platforms). Hence, a model evaluation is not straightforward and it is difficult to define meaningful model–data misfit metrics. Further complexity comes from the fact that misfit metrics must be tailored to the actual application of the model in order to be meaningful.

A potential application for the model configuration presented here is to explore the role of mesoscale features, or eddies, in determining the CO2 uptake of the Southern Ocean. One hypothesis this model is set up to test is whether spatially unresolved dynamics in IPCC-type coarse-resolution models bias their oceanic carbon uptake sensitivity. In order to come to a meaningful conclusion on this, our high-resolution model has to perform with a fidelity similar or superior to that of IPCC-type coarse-resolution models. In the following, we list and explain our choice of model assessments which we deem relevant in this respect. Wherever applicable, we will follow the approach taken by Sallée et al. (2013a, b), who assessed the performance of CMIP5 models, and Russell et al. (2018), who suggested evaluation metrics for the Southern Ocean in coupled climate models and Earth system models.

Please note that we refrain from using state estimates based on data assimilation such as the Biogeochemical Southern Ocean State Estimate (B-SOSE) (Verdy and Mazloff2017) for biogeochemical parameters. On the one hand, these state estimates are especially useful in the Southern Ocean, where observational data are particularly sparse. On the other hand, these estimates rely – especially when observational data are sparse – on the realism of the underlying ocean-circulation biogeochemical model framework. In the case of B-SOSE, the underlying model framework is so similar to our model configuration (high-resolution ocean-circulation model coupled to the biogeochemical model of Galbraith et al.2010) that we essentially would assess the benefit of the assimilation scheme in fitting observational data rather than assessing our model configuration.

The model–data comparison is divided into the following subsections:

  • Ocean circulation (Sect. 3.1), e.g., affects the transport of carbon-rich deepwater to the surface, shapes the locations of fronts and constitutes a major pathway for nutrients that are essential for phytoplankton growth. For the evaluation of surface currents, we use exemplary snapshots showcasing main circulation paths and spatial variability, along with climatological sea surface height which is indicative of the barotropic circulation. More quantitative measures of transport characteristics are provided for the ACC (Drake Passage transport) and the Southern Ocean meridional overturning circulation. Further, we assess the strength of the cyclonic Ross and Weddell polar gyres which are formed by interactions between the ACC and the Antarctic continental shelf. Special emphasis is here on the larger Weddell Gyre. This gyre entrains heat and salt from the ACC and carries them to the Antarctic continental shelves, where deep and bottom waters are produced and thus establish an intermittent connection between the atmosphere (surface ocean) and the deep oceanic carbon (nutrient) pool.

  • Eddy kinetic energy (EKE; Sect. 3.2) is an important measure for the mesoscale activity and thus a key proxy for realistically reproduced eddy dynamics. At the surface, the EKE can be derived from the variability of the sea surface height (SSH), a measure that can be directly observed from space by satellite altimetry.

  • Surface mixed layer depth (MLD; Sect. 3.3) is determined by (1) the stability of the water column (i.e., the vertical density gradient), (2) wind-induced turbulence which provides energy for eroding the stability and thereby deepens the MLD and (3) air–sea buoyancy fluxes which can, depending on their sign, increase or decrease the stability of the water column and thus shallow or deepen the MLD. The MLD is an important concept or metric because it locates that fraction of the upper water column that is in direct contact with the atmosphere and sets air–sea gradients of temperature and partial pressure of gases (such as carbon dioxide or oxygen). Further, changes in MLD modulate the level of average light levels experienced by phytoplankton: deep MLDs define conditions where phytoplankton spend more time mixed downwards away from the sunlit surface. This reduces phytoplankton growth and associated biotic carbon uptake. Antagonistic to this light deprivation, however, deepening MLDs are typically associated with transports of nutrients (essential for phytoplankton) from depth to the sunlit nutrient-depleted surface.

  • Temperatures (Sect. 3.4), which influence the rate of biological turnover, affect the solubility pump of carbon and the density of seawater. Cooling typically results in a reduction of buoyancy and causes convection. Spatial temperature gradients are associated with geostrophic circulation (if they are not salinity compensated). Vertical temperature gradients affect the stability of the water and precondition vertical mixing. The sea surface temperature (SST) can be directly observed from space and therefore is available in an unrivaled (compared to in situ measured properties) spatial and temporal resolution. In contrast, temperature observations at depth are unfortunately much sparser even though deep temperatures determine important processes, e.g., for basal melting, and thus are related to the formation rate of Antarctic Bottom Water (AABW). AABW formation, in turn, affects the solubility and biotic pump of carbon.

  • Salinity (Sect. 3.5) is related to the density of seawater. Saltening by brine rejection can cause convection; meltwater, on the other hand, can form buoyant lenses, thus increasing the local stability of the water column which prevents vertical mixing. Spatial salinity gradients are associated with geostrophic circulation (if they are not temperature compensated). Vertical salinity gradients affect the stability of the water column and precondition vertical mixing.

  • Sea ice (Sect. 3.6) caps the direct exchange between atmosphere and ocean and thus controls the air–sea gas exchange of CO2. It also modulates the air–sea buoyancy forcing by, e.g., insulating the surface from heat loss or by brine rejection during ice formation. Further, it shields the surface water from solar irradiance and hampers the assimilation of CO2 by autotrophic plankton.

  • Nutrients (Sect. 3.7) are essential to the growth of autotrophic plankton, and their availability exerts a major control on the biological pump. The most important macronutrient is bioavailable phosphorous (such as phosphate, PO4), because its availability is essential to all phytoplankton (and cyanobacteria). The distribution of PO4 is determined by the interplay between the ocean circulation (transporting PO4 dissolved in seawater) and marine biota which utilize phosphorous to build biomass (typically at the surface) and release PO4 in the course of degradation of organic material (typically at depth). In addition, we assess simulated iron concentrations since the Southern Ocean is well known for being a site where iron is limiting the growth of autotrophs (e.g., Boyd and Ellwood2010). Further, in order to guide the interpretation of simulated PO4 concentrations at depth, we include a meridional section of dissolved oxygen in this section: oxygen is reset to values close to saturation at the surface, while its sinks in the interior are typically assumed to be linearly linked to interior sources of PO4 (which is not reset at the surface). Hence, cases where biases in simulated PO4 concentrations are inconsistent with biases in simulated oxygen concentrations contain information on what exactly is the cause for the simulated biases within the interplay of ocean circulation and biogeochemistry.

Figure 4Climatological mean sea surface height in units of centimeters. Panel (a) shows a 1993 to 1998 climatological average observed from space (Maps of Absolute Geostrophic Current (MADT) Aviso data). Panel (b) shows a 6-year average (nominal years 1993 to 1998) simulated with the reference simulation. The white patches in panels (a) and (b) indicate missing data and (spurious) model land mask, respectively.

3.1 Ocean circulation

Figure 4 shows simulated climatological SSH, along with an estimate based on observations from space (Maps of Absolute Geostrophic Current (MADT) Aviso). The model captures the main features of the global barotropic circulation, even in the coarsely resolved domain north of 40 S, such as in the subtropical gyres in both hemispheres north and south of the Equator (positive SSH anomalies) and in the subpolar gyres around 50 and 60 N in the Pacific and Atlantic (negative SSH anomalies), respectively. Further, as indicated by the transition from yellow to blue in the Southern Ocean in Fig. 4, the model reproduces the northern boundary of the Antarctic Circumpolar Current (ACC), also referred to as the northern edge of the Subantarctic Front (SAF).

A quantitative numerical comparison of observed (Aviso) and simulated climatologic surface speeds yields a spatial correlation of 0.6 in the high-resolution domain south of 40 N. The spatial variance in the speed of simulated surface currents is typically 30 % higher than the observational estimate and the respective simulated mean is 20 % higher than in the observations. Apparently, a (yet-to-be-quantified) fraction of this misfit is attributable to errors in the estimates from space (cf. Fratantoni2001).

Figure 5 shows exemplary snapshots of the surface velocities as observed from space (Fig. 5a, MADT Aviso) and simulated (Fig. 5b, MOMSO). The boundary conditions of the model are not identical to the atmospheric boundary conditions on that specific day, and the highly non-linear characteristics of eddy dynamics render an “eddy-to-eddy” similarity without data assimilation impossible. So the purpose of Fig. 5 is to demonstrate the similarity of patterns and the major transport pathways in the high-resolution domain south of 40 N. The close-up into the Agulhas retroflection zone (Fig. 6), which is located right in the transition zone from very high to very coarse resolution (see Fig. 1c), highlights that the realism of simulated patterns is sustained right up into the transition to coarser, non-eddy-resolving resolution.

Figure 5Magnitude of surface velocities. Panels (a) and (b) show typical snapshots as observed from space (Aviso, 3 May 1995) and as simulated (MOMSO, 3 May of nominal year 2020), respectively.

Figure 6Magnitude of surface velocities. Panels (a) and (b) show typical snapshots as observed from space (Aviso, 3 May 1995) and as simulated (MOMSO, 3 May of nominal year 2020), respectively. This is a close-up from Fig. 5, focusing on those latitudes where the meridional resolution transitions from eddy permitting in the south to eddy prohibiting in the north.


Figure 7 features the simulated transport through the Drake Passage. (Here, we refer to the reference simulation only; simulation WIND is discussed in Sect. 4.) The simulated Drake Passage transport averages to 99 Sv (nominal years 1980–2023). This is biased low compared to available observational estimates that range from 110 to 170 Sv (Withworth1983 for 1979; Cunningham et al.2003 for 1993–2000; Chidichima et al.2014 for 2007–2011). Note that similar biases have been reported in other high-resolution configurations and may, according to Dufour et al. (2015), be related to a deficient representation of the overflow of dense waters formed along the Antarctic coasts. If so, the ACC bias may well be endemic to z-level models which struggle to represent complex topography (in comparison to more elaborate numerical approaches, e.g., finite elements). A comprehensive investigation is beyond the scope of this paper. But, still, the problem is an intriguing one – especially since, historically, (coarse-resolution) models started out from an opposing bias dubbed “Hidaka's dilemma” (Hidaka and Tsuchiya1953), where an excessive ACC transport could only by application of unrealistically high friction be fenced into realistic bounds.

Figure 7Simulated volume transport of the Antarctic Circumpolar Current through the Drake Passage in units of 106 m3 s−1. The black (red) line refers to the reference (increasing wind) simulation.


In terms of Eulerian meridional overturning in the Southern Ocean, our model values are consistent with the Southern Ocean state estimate of Mazloff et al. (2010): Mazloff et al. (2010) find a surface meridional overturning cell across 32 S of 12±12 Sv and an abyssal cell of 13±6 Sv. We find a climatological mean value for the upper cell of 12±4 and 8±3 Sv for the lower cell. Please note that more elaborate measures such as the physically more meaningful overturning calculated on density coordinates or its approximation given by the residual mean streamfunction (cf. McIntosh and McDougall1996; Viebahn and Eden2010) are beyond the scope of this paper (which is intended to serve as a model description only).

In terms of transports of the Weddell and Ross gyres, our simulation is slightly biased high. In the Weddell Gyre, we simulate 70 Sv, while published estimates range from 40±8 Sv (Southern Ocean state estimate; Mazloff et al.2010) to 56±10 Sv (recent SODA estimate Yongliang et al.2016), 55.9±9.8 and 61–66 Sv (Schröder and Fahrbach1999). In the Ross Gyre, we simulate 35 Sv, while published estimates range from 20±5 Sv (Mazloff et al.2010) to 15–30 Sv (Chu and Fan2007) and 37±6 Sv (Yongliang et al.2016).

Figure 8Eddy kinetic energy. Panel (a) is calculated from 1993–1998 satellite altimetry (Maps of Sea Level Anomalies and Geostrophic Velocity Anomalies (MSLA) Aviso data) observations. Panel (b) corresponds to a 6-year average (nominal years 1993 to 1998) calculated from daily averaged surface velocities from the reference simulation. The white patches in the left and right panels indicate missing data and (spurious) model land mask, respectively.

3.2 Eddy kinetic energy and sea surface height

Figure 8 shows that the simulated climatological EKE reproduces the amplitudes and spatial pattern observed from space (Maps of Sea Level Anomalies and Geostrophic Velocity Anomalies (MSLA) Aviso data) as can be expected from high-resolution configurations (e.g., Delworth et al.2012; Barnier et al.2006). This applies especially in the high-resolution region south of 40 S, where typical deviations in simulated EKE are rather small (<0.03 m2 s−2) compared to the maximum values found in the ACC (>0.3 m2 s−2; cf. Fig. 9). The correlation coefficient between simulated and observed EKE patterns is 0.7. The simulated EKE features 7 % more spatial variance and an average of 8 % more energy within the high-resolution domain south of 40 S. By averaging all simulated values with their immediate spatial neighbors two times consecutively, the simulated variance levels become more similar to observed values and the respective correlation increases to 0.71. We interpret this as an indication that the model resolution is higher than the effective spatial resolution of the MSLA Aviso data. This may also explain why simulated EKE amplitudes are generally higher than estimates from space (cf. Fratantoni2001).

Figure 9Difference between observed and simulated eddy kinetic energy. Blue colors denote regions where the observed levels calculated from 1993–1998 satellite altimetry (MSLA Aviso data) observations are lower than simulated levels. The region outside the high-resolution nest (north of 40 S) is masked by a white patch.

We conclude that we found no evidence of an underestimation of simulated mesoscale activity. This, in turn, implies that our combination of spatial resolution and (parameterization of) friction in the Southern Ocean is suitable to allow for investigations into eddy-driven processes.

Figure 10Surface mixed layer depth in units of meters. The upper (lower) line refers to austral winter (summer). The columns refer to (observed) CARS 2009 climatology based on Condie and Dunn (2006), Argo Al. and Argo Thresh. computed from (observed) profiles of Argo floats using two different methods by Holte et al. (2017) and our simulation (MOMSO). White patches in Argo Al. and Argo Thresh. denote missing data.

3.3 Surface mixed layer depth

MLD is a concept that locates that fraction of the upper water column that is well mixed. MLD is typically derived from vertical profiles of temperature and salinity rather than from direct measurements of mixing intensity (i.e., turbulence). As a consequence, algorithms calculating MLDs struggle with distinguishing between identifying actively mixed layers and those which have been actively mixed in the past but are no longer fed by energy fueling actual turbulence and associated mixing. The first three columns of Fig. 10 provide a comparison between contemporary (in the sense that they are still used in current peer-reviewed literature) databases and algorithms to compute MLD: CARS 2009 distributed by CSIRO (last access: 1 November 2019) and introduced by Condie and Dunn (2006) is predominantly based on historical conductivity–temperature–depth (CTD) observations, while both Argo Al. and Argo Thresh. are based on observations from Argo floats. The difference between Argo Al. and Argo Thresh. is the algorithm used to derive MLDs in Holte et al. (2017). We find a high correlation between the different algorithms Argo Al. and Argo Thresh. in the region south of 40 S (0.98 and 0.91 in summer and winter, respectively). The correlation from one database to another (i.e., Argo vs. CARS 2009) is much lower (≈0.7 for both winter and summer – irrespective of algorithm). Argo Thresh. is closest to our simulated pattern (fourth row in Fig. 10), and we find correlations of 0.5 and 0.6 in summer and winter, respectively.

In summary (cf. winter statistics in Fig. 11), we find that the average of variance in simulated MLDs is higher than in the observations and that the differences in correlation from one MLD product to another are comparable to our model–observation misfit. In terms of bias, we find that our simulated winter MLD is, on average, 18 m deeper than Argo Thresh. This is a major improvement compared to (coarse-resolution and fully coupled) CMIP5 models which feature similar correlations but winter mixed layer depths biased low by typically 50–100 m (Sallée et al.2013a, their Fig. 4a).

Figure 11Taylor diagram of surface mixed layer depth in austral winter referred to as Argo Thresh. observational data compiled by Holte et al. (2017). The units for the standard deviation and root mean square deviations are meters. Argo Al. refers also to observational data compiled by Holte et al. (2017) but based on a different algorithm. CARS 2009 is a climatology based on observations compiled by Condie and Dunn (2006), and MOMSO refers to the high-resolution reference simulation presented in this study.


3.4 Temperature

Figure 12 shows a comparison of the simulated SST with WOA09 (Locarnini et al.2010) and Argo-float-based (Roemmich and Gilson2009) observations. The overall pattern is well reproduced – an exception being a local high bias of about 3 C close to the Antarctic coast (between 120 and 160 E). Averaged over the Southern Ocean (Fig. 13), the simulated SSTs are within the observational range (although at the lower edge) of the years 1960 to 2010 (Fig. 14 calculated from HadISST; Rayner et al.2003).

Figure 12Climatological mean sea surface temperature. WOA09 and Argo (2004–2017 period) refer to observations compiled by Locarnini et al. (2010) and Roemmich and Gilson (2009), respectively. MOMSO refers to an average over the nominal 1993–1998 period of the reference simulation. The upper (lower) panels show sea surface temperature (differences).

Figure 15 summarizes the fidelity of simulated temperatures in the interior. Following Sallée et al. (2013b) and their density-based water mass criteria, we discuss the five major water masses (subtropical water, mode water, intermediate water, circumpolar deep water, bottom water) separately:

  • With subtropical water (TW), MOMSO features a slight cold bias at the surface (less than 0.5 C) and a warm bias of ≈1C in the interior (Fig. 15). This is apparently very good when compared against the 2.7 warm bias in the CMIP5 models (Sallée et al.2013b). There are, however, indications that the 2.7 bias is a consequence of biased air–sea heat fluxes in the Southern Hemisphere in the coupled models rather than the consequence of an unrealistic oceanic circulation.

  • The subduction of surface waters is dominated by mode water (MW) in the Southern Ocean. Hence, MW is an important agent in sequestering anthropogenic carbon away from the atmosphere. MOMSO's MW is too warm by ≈1C, which compares favorably to the much larger bias of an ensemble mean bias of 3.5C in CMIP5 models (Sallée et al.2013b).

  • Intermediate water (IW) is a denser type of MW which ventilates the thermocline and is sequestered away from the atmosphere longer than MW because it protrudes deeper into the water column. MOMSO's IW is biased warm by ≈1.1C (Fig. 15), which compares favorably against the 3.5 C warm bias in CMIP5 models (Sallée et al.2013b).

  • During its formation at the surface, circumpolar deep water (CDW) taps into the abyssal waters rich in natural carbon. Hence, it plays a key role for the oceanic sequestration of natural carbon dioxide (e.g., Le Quéré et al.2009). MOMSO overestimates the CDW temperature by ≈1C. This is substantially more than the 0.4 C overestimation of the CMIP5 ensemble mean but still within the envelope of the ensemble, which peaks at an overestimation of 1.5 C (Sallée et al.2013b).

  • MOMSO overestimates Antarctic Bottom Water (BW) temperatures by 1C, while the CMIP5 ensemble mean compares more favorably with an underestimation of only 0.4 C. These biases are consistent with MOMSO featuring too much sea ice coverage which caps the ocean from the cooling atmosphere (cf. Sect. 3.6). Likewise, CMIP5 models tend to underestimate sea ice coverage (Turner et al.2013), which may overly expose the ocean to the cooling atmosphere and thereby cause the respective cold bias in BW temperature.

Figure 13Simulated sea surface temperature averaged over the Southern Ocean (i.e., south of 40 S). The black (red) line refers to the reference (increasing wind) simulation.


Figure 14Observational estimate of annual mean sea surface temperatures south of 40 S calculated from HadISST (Rayner et al.2003).


Figure 15Zonally averaged meridional section of temperature in units of C. Panel (a) refers to simulated concentrations averaged over the nominal 1993–1998 period of the reference simulation; (b) refers to observed climatological values (WOA09); (c) refers to the difference between simulated and observed values, with blue colors denoting simulated temperatures that are biased low. The thick black, red, cyan, blue and green contours refer to mean densities (Sallée et al.2013b, their Table 2) of subtropical water, mode water, intermediate water, circumpolar deep water and bottom water, respectively.


3.5 Salinity

Figure 16 shows a comparison of the simulated climatological mean sea surface salinity (SSS) with three observation-based products (WOA09, Antonov et al.2009; Argo, Roemmich and Gilson2009; and SMOS, Köhler et al.2015). Within the spatially highly resolved Southern Ocean, the simulated sea surface salinity is in good agreement with the observations. Towards the north, where the resolution coarsens, the model fidelity disintegrates; particularly in the Atlantic and Indian sectors, the model is biased low.

Figure 16Climatological mean sea surface salinity in psu. WOA09, Argo (2004–2017 period) and SMOS in the upper panels refer to observations compiled by Antonov et al. (2009), Roemmich and Gilson (2009) and Köhler et al. (2015), respectively. MOMSO in the lower panel refers to an average over the nominal 1993–1998 period of the reference simulation. MOMSO-WOA09 and MOMSO-Argo refer to sea surface salinity differences between the reference simulation and respective observations.

Averaged over the Southern Ocean, the simulated SSS is biased low (compare reference in Fig. 17) compared to recent observational estimates during 2005–2017 (Fig. 18). Figure 17 suggests that increasing winds can increase simulated SSSs up to observed values. This could be interpreted as a mismatch between the climatological forcing (representing the years 1958–2000) and the observation period. (Please note that a comprehensive analysis is beyond the scope of this model description paper.) Figure 17 suggests further that the surface salinity restoration (cf. Sect. 2.6) is weak enough to allow for substantial SSS dynamics in response to, e.g., increasing winds.

Figure 17Simulated sea surface salinity averaged over the Southern Ocean (i.e., south of 40 S) for the nominal years 1980–2024. The black (red) line refers to the reference (increasing wind) simulation.


Figure 19 summarizes the fidelity of simulated salinities in the interior. Following Sallée et al. (2013b) and their density-based water mass criteria, we discuss the five major water masses (subtropical water, mode water, intermediate water, circumpolar deep water, bottom water) separately:

  • For TW, MOMSO is biased low by 0.5 at the surface, and this makes the simulated TW too light in comparison with WOA09 observations. Sallée et al. (2013b) find a smaller ensemble mean bias of 0.2 in the CMIP5 models which is – in contrast to MOMSO – enhanced, in terms of biasing the density towards smaller values, by a 2.7 C warm bias.

  • For MW, MOMSO is biased low by 0.1 at the surface and high by 0.1 in the interior. This is in line with the overall poor performance of CMIP5 models in representing the MW (Sallée et al.2013b). Coincidentally, the ensemble mean over all the CMIP5 models averages to a very small fresh bias of 0.02.

  • For IW, MOMSO is, again, biased low by 0.1 at the surface and high by 0.1 in the interior. This dampens the equatorward protruding of IW's salinity-minimum signature. Compared against the high bias of only 0.008 in the ensemble mean over all the CMIP5 models, MOMSO's bias appears substantial. Compared, however, against individual CMIP5 models whose characteristics vary widely (Sallée et al.2013b), MOMSO is comparable (i.e., within the envelope spanned by the CMIP5 models).

  • For CDW, MOMSO is biased high by 0.07 south of 75 S at the surface and less than 0.02 equatorwards. These small biases are well within the range set by the CMIP5 models and of the same order of magnitude as the bias of 0.02 of the CMIP5 ensemble.

  • For BW, MOMSO's bias is identical to the high bias of 0.06 in the CMIP5 ensemble mean (Sallée et al.2013b) and much better than the up-to-0.5 biases of many individual CMIP5 members. The bias is consistent with the bias in sea ice: an overestimated ice production releases too much brine and shields the ocean from atmospheric cooling. This drives spuriously elevated salinities and temperatures.

Figure 18Observed sea surface salinity averaged over the Southern Ocean (i.e., south of 40 S) based on Argo data (Roemmich and Gilson2009).


Figure 19Zonally averaged meridional section of salinity. Panel (a) refers to simulated salinities averaged over the nominal 1993–1998 period of the reference simulation. Panel (b) refers to observed climatological salinities (WOA09). Panel (c) refers to the difference between simulated and observed values, with blue colors denoting simulated salinities that are biased low. The thick black, red, cyan, blue and green contours refer to mean densities (Sallée et al.2013b, their Table 2) of subtropical water, mode water, intermediate water, circumpolar deep water and bottom water, respectively.


3.6 Sea ice

Figure 20 shows a comparison of the simulated number of sea-ice-covered months per year with an observation estimate (HadISST; Rayner et al.2003). Overall, the agreement is good, with the following two exceptions: (1) in the Weddell and Ross seas, the period of ice coverage is underestimated by 2 months. We speculate that this triggers elevated air–sea momentum fluxes and, eventually, biases the respective gyre strengths high (cf. Sect. 3.1). (2) Overall, the simulated ice extent is biased high (compare Fig. 21, black line, to Fig. 22). Following Russell et al. (2018) we compare the data, in addition, against the fractional sea ice coverage obtained from the National Snow and Ice Data Center (NSIDC; pub/DATASETS/NOAA/G02202_v2/, last access: 1 September 2019). The comparison in Fig. 23 confirms the impression that the model overestimates sea ice coverage. Figure 21 suggests that increasing the wind speeds to levels that are more representative of the time period of observations shown in Fig. 22 alleviates this model bias. In addition, Fig. 23 reveals that the model bias features a seasonality with underestimation in January until March and an overestimation of ice coverage from May to December. Ranked against the current generation of IPCC models, MOMSO fits well into the envelope of the ensemble reported by Turner et al. (2013).

Figure 20Ice-covered months in a year. Panel (a) refers to a 1980–2000 average based on the HadISST observational estimate (Rayner et al.2003). Panels (b) and (c) refer to the nominal year 2022 of the reference and increasing wind simulations, respectively.

Figure 21Simulated annual mean sea ice cover south of 40 S in units of 107 km2. The black (red) line refers to the reference (increasing wind) simulation.


Figure 22Observational estimate of annual mean sea ice cover south of 40 S in units of 107 km2 calculated from HadISST (Rayner et al.2003).


Figure 23Simulated (dashed line, reference simulation) and observed (solid line, reanalysis from NSIDC) climatological annual cycle of sea ice cover south of 40 S in units of 107 km2.


3.7 Nutrients

Simulated Southern Ocean PO4 surface concentrations are biased low, locally down to 0.6 mmol P m−3 (Fig. 24a, c). The reason is not straightforward to identify because it could be associated with a deficient physical module, a deficient biogeochemical module or both. In the following, we will present an indication that the problem is associated with a deficient formulation of iron limitation, argue why the formulation of light limitation is unlikely to be the main problem and put the model–data misfit into perspective.

Figure 24Comparison of simulated and observed phosphate concentrations. Panel (a) refers to simulated surface phosphate concentrations (reference simulation, averaged over nominal years 1993 to 1998). The color bar denotes phosphate concentrations in units of mmol P m−3. Panel (b) shows simulated (blue line) and observed (black line) phosphate concentrations averaged horizontally in the Southern Ocean (i.e., south of 40 S), along with their respective spatial standard deviations (grey and blue horizontal bars). Panel (c) shows the difference between simulated and observed (Garcia et al.2010) surface concentrations. Panel (d) shows the simulated (blue line) and observed (black line) seasonal cycle of dissolved phosphate concentration at the surface calculated as monthly means averaged over the Southern Ocean.

Southern Ocean uptake of PO4 in the sunlit surface by autotrophic phytoplankton is known to be limited by the availability of light and the availability of the micronutrient iron. Figure 25 features a comparison of simulated iron concentrations with observations. Even though the spatial and temporal coverage of iron measurements is still sparse, the emerging pattern is one where the simulated biotic iron drawdown at the surface appears to be too strong. Surface iron concentrations are biased low, just like the PO4 concentrations, and they appear so throughout an annual cycle. Such deficient model behavior can be caused by insufficient throttling of phytoplankton growth by both iron and light limitation. Looking closer at seasonal model–data misfits, however, suggests that a deficient iron limitation is more likely to be the cause: the dependency of growth (and associated micro- and macronutrient drawdown at the surface) is known to be a highly non-linear function of environmental drivers. We find that the bias in surface PO4 concentrations is almost constant over the course of a seasonal cycle (Fig. 24d), even though the photosynthetically available radiation varies dramatically from season to season in the respective latitudes (also because radiation experienced by phytoplankton cells dispersed in surface waters is a function of the seasonally varying surface mixed layer depth). By chance, this could be the result of non-linear forcing modulating a deficient non-linear formulation of PO4 limitation such that the model bias stays constant over a wide range of environmental conditions (here seasons). But this is unlikely. Looking into the seasonal bias of simulated iron concentrations (Fig. 25d), we find that it varies substantially from season to season compared to the respective PO4 variability (Fig. 24d) – just as is expected when a deficient non-linear model formulation is exposed to substantial variations in driving environmental conditions.

Figure 25Comparison of simulated with observed dissolved iron concentrations. Panel (a) refers to simulated surface iron concentrations (reference simulation, averaged over nominal years 1993 to 1998). The color bar denotes iron concentrations in units of nM Fe. Panel (b) shows simulated (blue line) and observed (black line) concentrations averaged horizontally in the Southern Ocean (i.e., south of 40 S), along with their respective spatial standard deviations (grey and blue horizontal bars). Panel (c) refers to observed surface iron concentrations (compiled by Tagliabue et al.2012, and Mawji et al.2014). (The color bar matches panel a.) Panel (d) depicts the simulated (blue line) and observed (black line) seasonal cycles of dissolved iron concentration at the surface calculated as monthly means averaged over the Southern Ocean.

The distribution of PO4 and dissolved oxygen in the interior (cf. Figs. 26 and 27) is consistent with an overly vivid biological pump that is not throttled enough by iron limitation: a prolonged overestimated biological production drives an overly high export of organic matter to depth (cf. spurious PO4 maximum around 400 m depth south of 60 S). Once in the thermocline, this PO4 is exported equatorwards and leaves the Southern Ocean (e.g., Marinov et al.2006). Becoming part of the meridional overturning circulation, this excess PO4 will return to the Southern Ocean at the surface. Because of insufficient iron limitation, the excess will be used up and exported to depth into, e.g., the Subantarctic Mode Water which leaves the Southern Ocean before it reaches higher latitudes. By this spurious mechanism, the Southern Ocean loses PO4 to the rest of the world's ocean.

Figure 26Zonally averaged meridional section of phosphate concentration in units of mmol P m−3. Panel (a) refers to simulated concentrations averaged over the nominal 1993–1998 period of the reference simulation; (b) refers to observed climatological concentrations (WOA09); (c) refers to the difference between simulated and observed concentrations, with blue colors denoting modeled concentrations that are biased low. The thick black, red, cyan, blue and green contours refer to mean densities (Sallée et al.2013b, their Table 2) of subtropical water, mode water, intermediate water, circumpolar deep water and bottom water, respectively.


Figure 27Zonally averaged meridional section of dissolved oxygen concentration in units of mmol O2 m−3. Panel (a) refers to simulated concentrations averaged over the nominal 1993–1998 period of the reference simulation; (b) refers to observed climatological concentrations (WOA09); (c) refers to the difference between simulated and observed concentrations, with blue colors denoting modeled concentrations that are biased low. The thick black, red, cyan, blue and green contours refer to mean densities (Sallée et al.2013b, their Table 2) of subtropical water, mode water, intermediate water, circumpolar deep water and bottom water, respectively.


The simulated patterns of dissolved oxygen bias in Fig. 27 are consistent with an overly efficient biological pump which strips the surface waters of PO4 on their way from the world's ocean to higher latitudes, exports them to depth where they are remineralized and imprints a spurious oxygen deficit. Superimposed onto this, as a rough calculation suggests, is a negative bias of the order of 5 mmol O2 m−3 caused by the temperature bias via the temperature dependence of oxygen solubility. Note that a more comprehensive analysis is beyond the scope of this model description which is intended to describe the configuration rather than evaluate it.

Figure 28 summarizes the model's fidelity in reproducing observed patterns of biogeochemical: ranked against the very similar, albeit coarser, non-eddying reference configuration used in Dietze et al. (2017) and UVic (here, referring to the reference version used in Löptien and Dietze2019, which is based on Keller et al.2012), MOMSO appears to perform in a very similar manner.

Figure 28Taylor diagram of simulated phosphate concentration referenced to Southern Ocean WOA09 annual mean climatology (Garcia et al.2010). The units for the standard deviation and root mean square deviations are µmol P m−3. “Coarse” refers to the very similar albeit coarser, non-eddying reference configuration used in Dietze et al. (2017) and “UVic” refers to the non-eddying Earth system model used, e.g., in climate engineering assessments (cf. reference version in Löptien and Dietze2019; Keller et al.2012).


3.8 Computational cost

Initial development, testing and spin-up of the physical configuration was carried out on two 32-core workstations, based on a 6320 AMD Opteron (Abu Dhabi) (8-core CPU, 2.80 GHz, 16 MB L3 cache, DDR3 1600) interconnected with a QDR InfiniBand. As of 2013, the system cost  EUR 10 000. The fully coupled ocean-circulation biogeochemical model has been run at a supercomputing center (project shk00027 at the Norddeutscher Verbund für Hoch- und Höchstleistungsrechnen, HLRN). Figure 29 documents excellent parallel performance up to 500 tasks – a setting which optimized our throughput on the machine/queueing system. On 500 Intel Xeon CPU E5-2670 cores (Sandy Bridge), it takes 10 h to simulate 1 year. Referring to pricing published by HLRN in 2018, this amounts to  EUR 250 per integrated model year. We use compiler settings that optimize for maximum computational efficiency. They violate computational reproducibility (cf. Intel Fortran compiler option -fp-model and summary by Corden and Kreitzer2018). A description of the associated computational uncertainty is provided in Sect. 3.8 of Dietze et al. (2014), which uses the same code (applied to a different domain).

4 Research questions

The purpose of this section is to outline research questions for which MOMSO may serve as a tool. MOMSO is (to our knowledge) unique in that it (1) features a high spatial resolution in the Southern Ocean which allows an eddying circulation in combination with (2) a biogeochemistry that is sufficiently equilibrated such that the remaining trend is small enough to allows for studying decadal signals in SO carbon uptake; and (3) it is a “free-forward-in-time-running” model in a sense that (projected or hindcasted) anomalies to air–sea fluxes can be applied easily (as opposed to in data-assimilated or inverse models where this can be difficult). In the following, we sketch – very roughly – two ideas.

4.1 Wind-induced changes in Southern Ocean carbon uptake: parameterizing versus resolving the effects of eddies

Winds have been changing in the past decades over the Southern Ocean and they are anticipated to do so in the future. Lovenduski et al. (2013) estimate a linear increase at a rate of 14 % in 50 years throughout the period 1958 to 2007. Dietze et al. (2017) find in a coarse-resolution model that such a trend in winds is associated with a decreasing trend in Southern Ocean natural carbon uptake between 4.2 and 5.4 PgC 1000−1 yr−2 – depending on the choice of the underlying eddy parameterization. An open question is how these estimates compare against a model configuration that is actually resolving the eddying circulation rather than just parameterizing it. The main problem preventing an answer to this question so far has been the computational expense associated with running the carbon dynamics of a realistic eddying coupled ocean-circulation biogeochemical model into equilibrium.

MOMSO is identical to the coarse-resolution configurations used in Dietze et al. (2017), with the only difference being the higher spatial resolution and eddying circulation. Figures 13, 17, 30 and 31 show that the forcing of MOMSO with wind anomalies corresponding to a linear increase at a rate of 14 % in 50 years, as estimated by Lovenduski et al. (2013), drives trends in sea surface temperatures, sea surface salinities and the Eulerian overturning which are clearly distinguishable from the remaining model drift. In line with these clearly detectable trends in physical properties in response to increasing winds, Fig. 32 reveals that the associated trend in the oceanic uptake of natural carbon is also clearly detectable: the remaining drift in the MOMSO reference simulation is 1.0 PgC 1000−1 yr−2, allowing for a clear detection of the effect of increasing winds in simulation WIND which yields a substantially larger decline of 5.0 PgC 1000−1 yr−2.

Putting the simulated Southern Ocean carbon uptake in response to increasing winds into perspective, we find that MOMSO yields a drift-corrected 5 PgC 1000−1 yr−2 1 PgC 1000−1 yr−2) = 4 PgC 1000−1 yr−2. Using the same code framework and exact same biogeochemical module (but much coarser spatial resolution), Dietze et al. (2017) find between 4.2 and 5.4 PgC 1000−1 yr−2 depending on the underlying eddy parameterization. So this is an indication that contemporary eddy parameterizations may provide results very similar to high-resolution eddying configurations.

Further investigations, which are beyond the scope of this sketch of a potential application of MOMSO, must include an in-depth analysis of eddy compensation, which is capable of canceling changes in wind-driven overturning (cf. Marshall and Radko2003; Viebahn and Eden2010; Abernathey et al.2011, 2016; Hallberg and Gnanadesikan2006; Thompson et al.2014; Tamsitt et al.2017).

4.2 Wind-induced anomalies in Antarctic continental shelf bottom water temperatures

The realistic seabed temperatures (shown in Fig. 33), in combination with its biogeochemical module, provide a starting point which – with further development – may be suited to explore the following feedback loop with MOMSO: atmospherically driven changes in ocean-circulation drive additional heat supply which leads to basal melting. The buoyant lens of meltwater may, for one, suppress the AABW formation (Williams et al.2016). This, in turn, affects oceanic carbon sequestration. Second, the meltwater carries bioavailable iron to the Southern Ocean – which affects oceanic primary productivity and the associated export of organic carbon to depth (Grotti et al.2005; Lannuzel et al.2008, 2010; Raiswell et al.2008; Smith et al.2007; Smith and Nelson1986; van der Merwe et al.2009).

Figure 34 shows that the simulated Antarctic continental shelf bottom water (ASBW) changes its temperature in response to a linear wind increase (at rate of 14 % in 50 years, as estimated by Lovenduski et al.2013). The magnitude of ASBW temperature is distinguishable from the persisting model drift in the reference simulation (cf. Fig. 34a, b) and simulated magnitudes are similar to recent observational estimates by Schmidtko et al. (2014) (cf. Fig. 34b, c). The sign, however, is reversed, which suggests that observed changes of ASBW temperatures cannot be attributed to the recent increases in wind strength alone. In any case, Fig. 34 suggests that the equilibration of the model has set in sufficiently so that model responses to decadal forcing variability are detectable.

5 Summary and conclusions

We set out to develop a near-global coupled ocean-circulation biogeochemical model which explicitly resolves effects of mesoscale eddy dynamics on the uptake of carbon in the Southern Ocean – in contrast to relying on respective eddy parameterizations. The setup is dubbed MOMSO, with MOM referring to the underlying framework (GFDL's Modular Ocean Model version MOM4p1) and SO referring to the Southern Ocean. We use the biogeochemical carbon module BLING developed by Galbraith et al. (2010).

Figure 29Computational performance as a function of CPU cores put to work simultaneously. On 500 Intel Xeon CPU E5-2670 cores, it takes 10 h to simulate 1 year.


Figure 30Maximum Southern Ocean Eulerian meridional overturning circulation south of 40 S in units of 106 m3 s−1. The black (red) line refers to the reference (increasing wind) simulation.


Overall, we find in a climatological simulation (REF) an eddying surface circulation that is in reasonable agreement with observations from space. Further, simulated temperatures and sea surface salinity show a close agreement with observations. Sea ice cover is biased high but is still in the range of observed values during particularly cold winters. A remaining caveat is a low bias in the Drake Passage (99 Sv compared to observational estimates ranging from 110 to 170 Sv). This may be associated with locally enhanced viscosity over rough topography at the exit of the Drake Passage (cf. Fig. 2) which has been necessitated by numerical stability. (Whether this is associated with our very low choice of the Smagorinsky isotropic viscosity and/or the lack of biharmonic viscosity/diffusion as suggested by reviewer 1 still warrants investigation.)

Figure 31Simulated linear trend of Southern Ocean Eulerian meridional overturning circulation effected by increasing winds during nominal years 1980 to 2024 in units of 106 m3 s−1 yr−1. Positive values denote increasing overturning.


Figure 32Simulated oceanic uptake of natural carbon south of 40 S. The black (red) line refers to the reference (increasing wind) simulation.


Figure 33Temperature of ASBW at the seabed for depths shallower than 1500 m. Panels (a) and (b) refer to the simulated 1993–1998 climatology and observations reproduced from Schmidtko et al. (2014), respectively.

Figure 34Temporal trend in the temperature of the ASBW at the seabed for depths shallower than 1500 m. Panel (a) refers to the drift that is still persistent in the reference simulation during nominal years 1980–2024. Panel (b) refers to the trend in simulation in the simulation with increasing winds corrected by the drift that still persists in the reference. Panel (c) refers to an observational estimate reproduced from Schmidtko et al. (2014).

The simulated biogeochemistry is also biased. Surface PO4 concentrations are too low in the Southern Ocean, which is indicative of a deficient formulation of the limitation of phytoplankton growth. Seasonally varying biases in simulated surface iron concentrations suggest that the problem is associated with an up-to-date, incomprehensive quantitative understanding of iron dynamics. The model performance with respect to biogeochemistry is, however, similar to what the state of the art is in coarse-resolution models.

MOMSO is a step forward in the field of eddying coupled ocean-circulation biogeochemical carbon modeling in that it allows to investigate the effects of decadal changes in atmospheric boundary conditions on the oceanic carbon uptake with an eddy-resolving model that features a realistic geometry and topography. Previous attempts have been hindered by the computational cost that is associated with running simulations into a semi-equilibrated state which features trends considerably lower than the climate signals (as affected by prescribed anomalies in boundary conditions) under investigation. To this end, MOMSO benefitted from (1) an ever-increasing ease of access to computing power that is associated with Moore's law and (2) the fortunate coincidence that the spun-up coarse-resolution restart from Dietze et al. (2017) was close enough to the equilibrated high-resolution state of MOMSO such that it provided a relatively equilibrated starting point.

In this model description paper, we showcased that the remaining drift in MOMSO's Southern Ocean carbon uptake is substantially lower than changes driven by typical decadal variability of atmospheric variations. We illustrate that the respective sensitivity experiment is suitable for a comparison to coarser-resolution model versions and present a first impression in Fig. 32. Similarly, we find that the drift in simulated temperatures of the ASBW is small enough to allow to link atmospheric decadal variability to oceanic temperature variations at the boundary to ice shelves around Antarctica.

We illustrate that the respective sensitivity experiment is suitable for a comparison to coarser-resolution model versions and present a first impression in Fig. 25.

On a final note, the computational cost associated with integrating MOMSO is substantial. With present-day hardware, full spin-ups are in the five-figure Euro range so that extensive tuning of the biogeochemical model, which still features some substantial misfits to observations, is uneconomical. Starting from the already spun-up state presented here, however, is in the reasonable range of EUR 2500 per decade. The corresponding energy needs (operation only) are of the order of 1000 kWh per decade. Using a carbon intensity of 450 gCO2 kWh−1 (e.g., Moro and Lonza2018), this yields an emission of around 0.5 t of CO2 per simulated decade.

Code and data availability

The circulation model code MOM4p1 is distributed by NOAA's Geophysical Fluid Dynamics Laboratory (, last access: 7 February 2017). We use the original code without applying any changes to it, apart from very minor changes (≈10 lines of code) which levitate sea ice (see Sect. 2.3) and increase local viscosity (see Sect. 2.2 and 2.5). The respective code, initial conditions, forcing, namelists and model output are accessible via, (Dietze2018).

Author contributions

HD and UL were involved in setting up and running the model configuration. All three authors contributed to the interpretation of model results and to the outlining and writing of the paper.

Competing interests

The authors declare that they have no conflict of interest.


Eric Galbraith, contributor to the MOM (, last access: 7 February 2017) community and developer of BLING (, last access: 11 March 2017), shared his biogeochemical model code with us. We are grateful to him and the rest of the MOM community! Discussions with Ivy Frenger are appreciated. The authors acknowledge the North-German Supercomputing Alliance (HLRN) for providing HPC resources that have substantially contributed to the research results reported in this paper. The authors acknowledge computing time on the joint GEOMAR/University Kiel NEC HPC Linux cluster. Initial configuration development and testing have been carried out on the FB2/BM compute clusters ( and, located at the GEOMAR Helmholtz Centre for Ocean Research, west shore campus, Kiel, Germany. We acknowledge help from Kai Grunau for cluster maintenance and IT services. The altimeter products were produced by Ssalto/Duacs and distributed by Aviso, with support from CNES (, last access: 1 September 2019). SMOS ocean surface salinity was distributed in NetCDF format by the Integrated Climate Data Center (ICDC;, last access: 1 September 2019) University of Hamburg, Hamburg, Germany. Some of the data used in this study were collected and made freely available by the international Argo program and the national programs that contribute to it (, last access: 1 September 2019;, last access: 1 September 2019). The Argo program is part of the Global Ocean Observing System. We are grateful to all observational programs pursuing an open data policy. Constructive comments by two knowledgable reviewers helped to improve the manuscript substantially. Finally, we thank our editor, Paul Halloran.

Financial support

This research has been supported by the Helmholtz Gemeinschaft (grant no. ZT-I-0010) and the Deutsche Forschungsgemeinschaft (grant no. SCHN 762/5-1).

The article processing charges for this open-access
publication were covered by a Research
Centre of the Helmholtz Association.

Review statement

This paper was edited by Paul Halloran and reviewed by two anonymous referees.


Abernathey, R. J., Marshal, J., and Ferreira, D.: The dependence of Southern Ocean meridional overturning on wind stress, Ocean Model., 72, 1–16,, 2011. a, b

Abernathey, R. P., Cerovecki, I., Holland, P. R., Newsom, E. Mazloff, M., and Talley, L. D.: Water-mass transformation by sea ice in the upper branch of the Southern Ocean overturning, Nature, 9, 596–601,, 2016. a, b, c

Adcroft, A. and Campin, J.-M.: Rescaled height coordinates for accurate representation of free-surface flows in ocean circulation models, Ocean Model., 7, 269–284,, 2004. a

Anderson, R. F., Ali, S., Bradtmiller, L. I., Nielsen, S. H. H., Fleisher, M. Q., Anderson, B. E., and Burckle, L. H.: Wind-driven upwelling in the Southern Ocean and the deglacial rise in atmospheric CO2, Science, 323, 1443–1448,, 2009. a

Antonov, J. I., Seidov, D., Boyer, T. P., Locarnini, A., Mishonov, A. V., Garcia, H. E., Baranova, O. K., Zweng, M. M., and Johnson, D. R.: World Ocean Atlas 2009, Volume 2: Salinity, edited by: Levitus, S., NOAA Atlas NESDIS 69, U.S. Government Printing Office, Washington DC, 184 pp., 2010. a, b, c, d, e, f

Arakawa, A. and Lamb, V. R.: Computational design of the basic dynamical processes of the UCLA general circulation model, Methods in Computational Physics, edited by: Chang, J., Academic Press, Vol. 17, 173–265,, 1977. a

Barnier, B., Madec, G., Penduff, T., Molines, J.-M., Treguier, A.-M., Le Sommer, J., Beckmann, A., Biastoch, A., Böning, C., Dengg, J., Derval, C., Durand, E., Gulev, S., Remy, E., Tlandier, C., Theetten, S., Maltrud, M., McClean, J., and De Cuevas, B.: Impact of partial steps and momentum advection schemes in a global ocean circulation model at eddy permitting resolution, Ocean Dynam., 56, 543–567,, 2006. a

Bernardello, R., Marinov, I., Palter, J. P., Sarmiento, J. L., Galbraith, E. D., and Slater, R. .D.: Response of the Ocean Natural Carbon Storage to Projected Twenty-First-Century Climate Change, J. Climate, 27, 2033–2053,, 2014. a

Böning, C. W., Dispert, A., Visbeck, M., Rintoul, S. R., and Schwarzkopf, F. U: The response of the Antarctic Circumpolar Current to recent climate change, Nat. Geosci., 1, 864–869,, 2008. a

Boyd, P. W. and Ellwood, M. J.: The biogeochemical cycle of iron in the ocean, Nat. Geosci., 3, 675–682,, 2010. a

Campin, J.-M., Adcroft, A., Hill, C., and Marshall, J.: Conservation of properties in a free-surface model, Ocean Model., 221–244,, 2004. a

Carpenter Jr., R. L., Droegmeier, K. K., Woodward, P. R., and Hane, C. E.: Application of the Piecewise Parabolic Method (PPM) to Meteorological Modeling, Mon. Weather Rev., 118, 586–612,<0586:AOTPPM>2.0.CO;2, 1989. a

Chidichimo, M. P., Donohue, K. A., Watts, D. R., and Tracey, K. L.: Baroclinic transport time series of the Antarctic Circumpolar Current measured in Drake Passage, J. Phys. Oceanogr., 44, 1829–1853,, 2014. a

Chu, P. C. and Fan, C.: An inverse model for calculation of global volume transport from wind and hydrographic data, J. Marine Syst., 65, 376–399,, 2007. a

Colella, P. and Woodward, P. R.: The Piecewise Parabolic Method (PPM) for Gas-Dynamical Simulations, J. Comput. Phys., 54, 174–201,, 1984. a

Condie, S. A. and Dunn, J. R.: Seasonal characteristics of the surface mixed layer in the Australasian region: implications for primary production regimes and biogeography, Mar. Freshwater Res., 57, 569–590,, 2006. a, b, c

Corden, J. and Kreitzer, D.: Consistency of Floating-Point Results using the Intel Compiler or Why doesn't my application always give the same answer? Software Services Group, Intel Corporation, 1–25, available at: (last access: 6 January 2020), 2018. a

Cunningham, S. A., Alderson, S. G., King, B. A., and Brandon, M. A.: Transport and variability of the Antarctic Circumpolar Current in the Drake Passage, J. Geophys. Res., 108, 8084,, 2003. a

Delworth, T. L., Rosati, A., Anderson, W., Adcroft, A. J., Balaji, V., Benson, R., Dixon, K., Griffies, S. M., Lee, H.-C., Pacanowski, R. C., Vecchi, G. A., Wittenberg, A. T., Zeng, F., and Zhang, R.: Simulated Climate and Climate Change in the GFDL CM2.5 High-Resolution Coupled Climate Model, J. Climate, 25, 2755–2781,, 2012. a

DeVries, T., Holzer, M., and Primeau, F.: Recent increase in oceanic carbon uptake driven by weaker upper-ocean overturning, Nature, 542, 215–218,, 2017. a

Dietze, H.: Dietze et al. 2018, available at:, last access: 19 December 2018. a

Dietze, H. and Kriest, I.: 137Cs off Fukushima Dai-ichi, Japan – model based estimates of dilution and fate, Ocean Sci., 8, 319–332,, 2012. a, b

Dietze, H. and Loeptien, U.: Revisiting “nutrient trapping” in global coupled biogeochemical ocean circulation models, Global Biogeochem. Cy., 27, 1–20,, 2013. a

Dietze, H., Löptien, U., and Getzlaff, K.: MOMBA 1.1 – a high-resolution Baltic Sea configuration of GFDL's Modular Ocean Model, Geosci. Model Dev., 7, 1713–1731,, 2014. a, b, c

Dietze, H., Getzlaff, J., and Löptien, U.: Simulating natural carbon sequestration in the Southern Ocean: on uncertainties associated with eddy parameterizations and iron deposition, Biogeosciences, 14, 1561–1576,, 2017. a, b, c, d, e, f, g, h, i, j

Dufour, C. O., Griffies, S. M., de Souza, G. F., Frenger, I., Morrison, A. K., Palter, J. B., Sarmiento, J. L., Galbraith, E. D., Dunne, J. P., Anderson, W. G., and Slater, R. D: Role of mesoscale eddies in cross-frontal transport of heat and biogeochemical tracers in the Southern Ocean, J. Phys. Oceanogr., 45, 3057–3081,, 2015. a

Fratantoni, D. M.: North Atlantic surface circulation during the 1990's observed with satellite-tracked drifters, J. Geophys. Res., 106, 22067–22093,, 2001. a, b

Galbraith, E. D., Gnanadesikan, A., Dunne, J. P., and Hiscock, M. R.: Regional impacts of iron-light colimitation in a global biogeochemical model, Biogeosciences, 7, 1043–1064,, 2010. a, b, c, d, e, f, g

Garcia, H. E., Locarnini, R. A., Boyer, T. P., and Antonov, J. I.: World Ocean Atlas 2009, Volume 4: Nutrients (phosphate, nitrate, silicate), edited by: Levitus, S., NOAA Atlas NESDIS 71, U.S. Government Printing Office, Washington DC, 398 pp., 2010. a, b

Getzlaff, J., Dietze, H., and Oschlies, A.: Simulated effects of southern hemispheric wind changes on the Pacific oxygen minimum zone, Geophys. Res. Lett., 43, 728–734,, 2016. a

Griffies, S. M.: Elements of MOM4p1. GFDL Ocean Group Technical Report No. 6, NOAA/Geophysical Fluid Dynamics Laboratory, 444 pp., Version 16 December 2009. a, b, c

Griffies, S. M. and Hallberg, R. W.: Biharmonic friction with a Smagorinsky-like viscosity for use in large-scale eddy-permitting ocean models, Mon. Weather Rev., 128, 2935–2946,<2935:BFWASL>2.0.CO;2, 2000. a, b

Grotti, M., Soggia, F., Ianni, C., and Frache, R.: Trace metals distribution in coastal sea ice of Terra Nova Bay, Ross Sea, Antarctica, Antarct. Sci., 17, 289–300,, 2005. a

Hall, A. and Visbeck, M.: Synchronous variability in the Southern Hemisphere atmosphere, sea ice, and ocean resulting from the Annular Mode, J. Climate, 15, 3043–3057,<3043:SVITSH>2.0.CO, 2002. a

Hallberg, R. and Gnanadesikan, A.: The role of eddies in determining the structure and response of the wind-driven southern hemisphere overturning: results from the modeling eddies in the southern ocean (MESO) project, J. Phys. Oceanogr., 36, 2232–2252,, 2006. a, b, c

Heinze, C., Meyer, S., Goris, N., Anderson, L., Steinfeldt, R., Chang, N., Le Quéré, C., and Bakker, D. C. E.: The ocean carbon sink – impacts, vulnerabilities and challenges, Earth Syst. Dynam., 6, 327–358,, 2015. a

Hidake, K. and Tsuchiya, M.: On the Antarctic Circumpolar Current, J. Mar. Res., 12, 214–222, 1953. a

Hogg, A. McC., Meredith, M. P., Blundell, J. R., and Wilson, C.: Eddy heat flux in the Southern Ocean: Response to variable wind forcing, J. Climate, 21, 608–620,, 2008. a

Holte, J., Talley, L. D., Gilson, J., and Roemmich, D.: An Argo mixed layer climatology and database, Geophys. Res. Lett., 44, 5618–5626,, 2017. a, b, c, d

Hordoir, R., Axell, L., Höglund, A., Dieterich, C., Fransner, F., Gröger, M., Liu, Y., Pemberton, P., Schimanke, S., Andersson, H., Ljungemyr, P., Nygren, P., Falahat, S., Nord, A., Jönsson, A., Lake, I., Döös, K., Hieronymus, M., Dietze, H., Löptien, U., Kuznetsov, I., Westerlund, A., Tuomi, L., and Haapala, J.: Nemo-Nordic 1.0: a NEMO-based ocean model for the Baltic and North seas – research and operational applications, Geosci. Model Dev., 12, 363–386,, 2019. a

Hunke, E. C. and Dukowicz, J. K.: An Elastic-Viscous-Plastic Model for Sea Ice Dynamics, J. Phys. Oceanogr., 27, 1849–1867,<1849:AEVPMF>2.0.CO;2, 1997. a

Jochum, M., Danabasoglu, G., Holland, M., Kwon, Y.-O., and Large, W. G.: Ocean viscosity and climate, J. Geophys. Res., 113, C06017,, 2008. a

Keller, D. P., Oschlies, A., and Eby, M.: A new marine ecosystem model for the University of Victoria Earth System Climate Model, Geosci. Model Dev., 5, 1195–1220,, 2012. a, b

Köhler, J., Martins, M. S., Serra, N., and Stammer, D.: Quality assessment of spaceborne sea surface salinity observations over the northern North Atlantic, J. Geophys. Res.-Oceans, 120, 94–112,, 2015. a, b

Langlais, C. E., Lenton, A., Matear, R., Monselesan, D., and Legresy, B., Cougnon, E., and Rintoul, S.: Stationary Rossby waves dominate subduction of anthropogenic carbon in the Southern Ocean, Nature, 7, 17076,, 2017. a

Lannuzel, D., Schoemann, V., de Jong, J., Chou, L., Delille, B., Becquevort, S., and Tison, J.: Iron study during a time series in the western Weddell pack ice, Mar. Chem., 108, 85–95,, 2008. a

Lannuzel, D., Schoemann, V., de Jong, J., Pasquer, B., van der Merwe, P., Masson, F., Tison, J.-L., and Bowie, A.: Distribution of dissolved iron in Antarctic sea ice: Spatial, seasonal, and interannual variability, J. Geophys. Res., 115, G03022,, 2010. a

Large, W. G. and Yeager, S. G.: Diurnal to decadal global forcing for ocean and sea-ice models: the data sets and flux climatologies, NCAR Technical Note, TN-460, 2004. a

Large, W. G., McWilliams, J. C., and Doney, S. C.: Oceanic vertical mixing – A review and a model with nonlocal boundary-layer parameterization, Rev. Geophys., 32, 363–403,, 1994. a

Lenton, A. and Matear, R. J.: Role of the southern annular mode 80 (SAM) in Southern Ocean CO2 uptake, Global Biogeochem. Cy., 21, GB2016,, 2007. a

Le Quéreé, C., Raupach, M. R., Canadell, J. G., Marland, G., Bopp, L., Ciais, P., Conway, T. J., Doney, S. C., Feely, R. A., Foster, P., Friedlingstein, P., Gurney, K., Houghton, R. A., House, J. I., Huntingford, C., Levy, P. E., Lomas, M. R., Majkut, J., Metzl, N., Ometto, J. P., Peters, G. P., Prentice, I. C., Randerson, J. T., Running, S. W., Sarmiento, J. L., Schuster, U., Sitch, S., Takahashi, T., Viovy, N., van der Werf, G. R., and Woodward, F. I.: Trends in the sources and sinks of carbon dioxide, Nat. Geosci., 2, 831–836,, 2009. a

Lévy, M., Estublier, A., and Madec, G.: Choice of an Advection Scheme for Biogeochemical Models, Geophys. Res. Lett., 28, 3725–3728,, 2001. a

Limpasuvan, V. and Hartmann, D. L.: Eddies and the annular modes of climate variability, Geophys. Res. Lett., 26, 3133–3136,, 1999. a

Liu, N., Eden, C., Dietze, H., Wu, D., and Lin, X.: Model-based estimates of the heat budget in the East China Sea, J. Geophys. Res.-Oceans, 115, C08026,, 2010. a

Locarnini, R. A., Mishonov, A. V., Antonov, J. I., Boyer, T. P., Garcia, H. E., Baranova, O. K., Zweng, M. M., and Johnson, S. G: World Ocean Atlas 2009, Volume 1: Temperature, edited by: Levitus, S., NOAA Atlas NESDIS 69, U.S. Government Printing Office, Washington DC, 184 pp., 2010. a, b, c, d, e

Löptien, U. and Dietze, H.: Constraining parameters in marine pelagic ecosystem models – is it actually feasible with typical observations of standing stocks?, Ocean Sci., 11, 573–590,, 2015. a

Löptien, U. and Dietze, H.: Effects of parameter indeterminancy in pelagic biogeochemical modules of Earth System Models on projections into a warming future: The scale of the problem, Global Biogeochem. Cy., 31, 1155–1172,, 2017. a

Löptien, U. and Dietze, H.: Reciprocal bias compensation and ensuing uncertainties in model-based climate projections: pelagic biogeochemistry versus ocean mixing, Biogeosciences, 16, 1865–1881,, 2019. a, b, c

Lovenduski, N. S., Gruber, N., and Doney, S. C.: Toward a mechanistic understanding of the decadal trends in the Southern Ocean carbon sink, Global Biogeochem. Cy., 22, GB3016,, 2008. a

Lovenduski, N. S., Long, M. C., Gent, P. R., and Lindsay, K.: Multi-decadal trends in the advection and mixing of natural carbon in the Southern Ocean, Geophys. Res. Lett., 40, 139–142,, 2013. a, b, c, d, e

Marinov, I., Gnanadesikan, A., Toggweiler, J. R., and Sarmiento, J. L.: The Southern Ocean biogeochemical divide, Nature, 441, 964–967,, 2006. a

Marshall, G. J.: Trends in the Southern Annular Mode from observations and reanalyses, J. Climate, 16, 4134–4143,<4134:TITSAM>2.0.CO;2, 2003. a

Marshall, J. and Radko, T.: Residual-Mean Solutions for the Antarctic Circumpolar Current and Its Associated Overturning Circulation, J. Phys. Oceanogr., 33, 11, 2341–2354,<2341:RSFTAC>2.0.CO;2, 2003. a, b

Marshall J. and Speer K.: Closure of the meridional overturning circulation through Southern Ocean upwelling, Nat. Geosci., 5, 171–180,, 2012. a

Mawji, E., Schlitzer, R., Dodas, E. M., Abadie, C., Abouchami, W., Anderson, R. F., Baars, O., Bakker, K., Baskaran, M., Bates, N. R., Bluhm, K., Bowie, A., Bown, J., Boye, M., Boyle, E. A., Branellec, P., Bruland, K. W., Brzezinski, M. A., Bucciarelli, E., Buesseler, K., Butler, E., Cai, P., Cardinal, D., Casciotti, K., Chaves, J., Cheng, H., Chever, F., Church, T. M., Colman, A. S., Conway, T. M., Croot, P. L., Cutter, G. A., de Baar, H. J. W., de Souza, G. F., Dehairs, F., Deng, F., Dieu, H. T., Dulaquais, G., Echegoyen-Sanz, Y., Edwards, R. L., Fahrbach, E., Fitzsimmons, J., Fleisher, M., Frank, M., Friedrich, J., Fripiat, F., Galer, S. J. G., Gamo, T., Solsona, E. G., Gerringa, L. J. A., Godoy, J. M., Gonzalez, S., Grossteffan, E., Hatta, M., Hayes, C. T., Heller, M. I., Henderson, G., Huang, K.-F., Jeandel, C., Jenkins, W. J., John, S., Kenna, T. C., Klunder, M., Kretschmer, S., Kumamoto, Y., Laan, P., Labatut, M., Lacan, F., Lam, P. J., Lannuzel, D., le Moigne, F., Lechtenfeld, O. J., Lohan, M. C., Lu, Y., Masqué, P., McClain, C. R., Measures, C., Middag, R., Moffett, J., Navidad, A., Nishioka, J., Noble, A., Obata, H., Ohnemus, D. C., Owens, S., Planchon, F., Pradoux, C., Puigcorbé, V., Quay, P., Radic, A., Rehkämper, M., Remenyi, T., Rijkenberg, M. J. A., Rintoul, S., Robinson, L. F., Roeske, T., Rosenberg, M., Rutgers van der Loeff, M., Ryabenko, E., Saito, M. A., Roshan, S., Salt, L., Sarthou, G., Schauer, U., Scott, P., Sedwick, P. N., Sha, L., Shiller, A. M., Sigman, D. M., Smethie, W., Smith, G. J., Sohrin, Y., Speich, S., Stichel, T., Stutsman, J., Swift, J. H., Tagliabue, A., Thomas, A., Tsunogai, U., Twining, B. S., van Aken, H. M., van Heuven, S., van Ooijen, J., van Weerlee, E., Venchiarutti, C., Voelker, A. H. L., Wake, B., Warner, M. J., Malcolm, E., Woodward, S., Wu, J., Wyatt, N., Yoshikawa, H., Zheng, X.-Y., Xue, Z., Zieringer, M., and Zimmer, L. A.: The GEOTRACES Intermediate Data Product 2014, Mar. Chem., 177, 1–8,, 2015. a

Mazloff, M. R., Heimbach, P., and Wunsch, C.: An Eddy-Permitting Southern Ocean State Estimate, J. Phys. Oceanogr., 40, 880–899,, 2010. a, b, c, d

McIntosh, P. and McDougall, T.: Isopycnal averaging and the residual mean circulation, J. Phys. Oceanogr., 26, 1655–1660,<1655:IAATRM>2.0.CO;2, 1996. a

Moro, A., and Lonza, L.: Electricity carbon intensity in European Member States: Impacts on GHG emissions of electric vehicles, Transport. Res. D-Tr. E, 64, 5–14,, 2018. a

Munday, D. R., Johnson, H. L., and Marshall, D. P.: Impacts and effects of mesoscale ocean eddies on ocean carbon storage and atmospheric pCO2, Global Biogeochem. Cy., 28, 877–896,, 2014. a

Orr, J. C., Najjar, R. G., Aumont, O., Bopp, L., Bullister, J. L., Danabasoglu, G., Doney, S. C., Dunne, J. P., Dutay, J.-C., Graven, H., Griffies, S. M., John, J. G., Joos, F., Levin, I., Lindsay, K., Matear, R. J., McKinley, G. A., Mouchet, A., Oschlies, A., Romanou, A., Schlitzer, R., Tagliabue, A., Tanhua, T., and Yool, A.: Biogeochemical protocols and diagnostics for the CMIP6 Ocean Model Intercomparison Project (OMIP), Geosci. Model Dev., 10, 2169–2199,, 2017. a

Pickard, G. L. and Emery, W. J.: Descriptive Physical Oceanography, An Introduction, Pergamon Press, 5th Edn., 173–176, 1990. a

Raiswell, R., Benning, L. G., Tranter, M., and Tulaczyk, S.: Bioavailable iron in the Southern Ocean: The significance of the iceberg conveyor belt, Geochem. T., 9, 7,, 2008. a

Rayner, N. A., Parker, D. E., Horton, E. B., Folland, C. K., Alexander, L. V., Rowell, D. P., Kent, E. C., and Kaplan, A.: Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century, J. Geophys. Res., 108, 4407,, 2003. a, b, c, d, e

Ritter, R., Landschützer, P., Gruber, N.,Fay, A. R., Lida, Y., Jones, S., Nakaoka, S., Park, G.-H., Peylin, P., Rödenbeck, C., Rodgers, K. B., Shutler, J. D., and Zeng, J.: Observation-based trends of the Southern Ocean carbon sink, Geophys. Res. Lett., 44, 12339–12348,, 2017. a

Roemmich, D. and Gilson, J.: The 2004–2008 mean and annual cycle of temperature, salinity, and steric height in the global ocean from the Argo Program, Prog. Oceanogr., 82, 81–100,, 2009. a, b, c, d, e

Russell, J., Kamenkovich, I., Bitz, C., Ferrari, R., Gille, S. T., Goodman, P. J., Hallberg, R., Johnson, K., Khazmutdinova, K., Marinov, I., Mazloff, M., Riser, S., Sarmiento, J. L., Speer, K., Talley, L. D., and Wanninkhof, R.: Metrics for the Evaluation of the Southern Ocean in Coupled Climate Models and Earth System Models, J. Geophys. Res.-Oceans, 123, 3120–3143,, 2018. a, b

Saenko, O. A., Fyfe, J. C., and England, M. H.: On the response of the oceanic wind-driven circulation to atmospheric CO2 increase, Climate Dynam., 25, 415–426,, 2005. a

Sallée, J.-B., Shuckburgh, E., Bruneau, N., Meijers, A. J. S., Bracegirdle, T. J., and Wang, Z.: Assessment of Southern Ocean mixed-layer depths in CMIP5 models: Historical bias and forcing response, J. Geophys. Res.-Oceans, 118, 1845–1862,, 2013a. a, b

Sallée, J.-B., Shuckburgh, E., Bruneau, N., Meijers, A. J. S., Bracegirdle, T. J., Wang, Z., and Roy, T.: Assessment of Southern Ocean water mass circulation and characteristics in CMIP5 models: Historical bias and forcing response, J. Geophys. Res.-Oceans, 118, 1830–1844,, 2013b. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o

Schmidtko, S., Heywood, K. J., Thompson, A. F., and Aoki, S.: Multidecadal warming of Antarctic waters, Science, 346, 1227–1231,, 2014. a, b, c

Schröder, M. and Fahrbach, E.: On the structure and the transport of the eastern Weddell Gyre, Deep-Sea Res. Pt. II, 46, 501–527,, 1999. a

Screen, J. A., Gillett, N. G., Stevens, D. P., Marshall, G. J., and Roscoe, H. K.: The role of eddies in the Southern Ocean temperature response to the Southern Annular Mode, J. Climate, 22, 806–818,, 2009. a

Shapiro, R.: Smoothing, filtering, and boundary effects, Rev. Geophys. Space Phys., 8, 359–387,, 1970. a

Simpkins, G. R. and Karpechko, A. Y.: Sensitivity of the southern annular mode to greenhouse gas emission scenarios, Clim. Dynam., 38, 563–572,, 2012. a

Smagorinsky, J.: General circulation experiments with the primitive equations: I. The basic experiment, Mon. Weather Rev., 91, 99–164,<0099:GCEWTP>2.3.CO;2, 1963. a

Smagorinsky, J.: Some historical remarks on the use of nonlinear viscosities, in: Large Eddy Simulation of Complex Engineering and Geophysical Flows, edited by: Galperin, B. and Orszag, S. A., Cambridge University Press, 3–36, 1993. a

Smith, K., Robison, B. H., Helly, J. J., Kaufmann, R. S., Ruhl, H. A., Shaw, T. J., Twining, B. S., and Vernet, M.: Free-drifting icebergs: Hot spots of chemical and biological enrichment in the Weddell Sea, Science, 317, 5837, 478–482,, 2007. a

Smith, W. and Nelson, D.: Importance of ice edge phytoplankton production in the Southern Ocean, BioScience, 36, 251–257,, 1986. a

Spence, P., Saenko, O. A., Eby, M., and Weaver, A. J.: The Southern Ocean Overturning: Parameterized versus Permitted Eddies, J. Phys. Oceanogr., 39, 1634–1651,, 2009. a

Spence, P., Fyfe, J. C., Montenegro, A., and Weaver, A. J.: Southern Ocean Response to Strengthening Winds in an Eddy-Permitting Global Climate Model, J. Climate, 23, 5332–5343,, 2010. a

Stacey, M. W., Pond, S., and Nowak, Z. P.: A numerical model of the circulation in Knight Inlet British Columbia, Canada, J. Phys. Oceanogr., 25, 1037–1062,<1037:ANMOTC>2.0.CO;2, 1995. a

Sweby, P. K.: High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws, SIAM J. Numer. Anal., 21, 995–1011,, 2012. a

Tagliabue, A., Mtshali, T., Aumont, O., Bowie, A. R., Klunder, M. B., Roychoudhury, A. N., and Swart, S.: A global compilation of dissolved iron measurements: focus on distributions and processes in the Southern Ocean, Biogeosciences, 9, 2333–2349,, 2012. a

Takahashi, T., Sweeney, C., Hales, B., Chipman, D. W., Newberger, T., Goddard, J. G., Iannuzzi, R. A., and Sutherland, S. C.: The changing carbon cycle in the Southern Ocean, Oceanography, 25, 26–37,, 2012. a

Tamsitt V., Drake, H. F., Morrison, A. K., Talley, L. D., Dufour, C. O., Gray, A. R., Griffies, S. M., Mazloff, M. R., Sarmiento, J. L., Wang, J., and Weijer, W.: Spiraling pathways of global deep waters to the surface of the Souther Ocean, Nat. Commun., 8, 172,, 2017. a, b, c

Thompson, D. W. J. and Solomon, S.: Interpretation of Recent Southern Hemisphere Climate Change, Science 296, 895–899,, 2002. a, b

Thompson, D. W. J., Solomon, S., Kushner, P. J., England, M. H., Grise, K. M., and Karoly, D. J.: Signatures of the Antarctic ozone hole in Southern Hemisphere surface climates change, Nat. Geosci., 4, 741–749,, 2014.  a, b

Turner, J., Bracegirdle, T. J., Phillips, T., Marshall, G. J., and Hosking, J. S.: An Initial Assessment of Antarctic Sea Ice Extent in the CMIP5 Models, J. Climate, 26, 1473–1484,, 2013. a, b

van der Merwe, P., Lannuzel, D., Nichols, C. A. M., Meiners, K., Heil, P., Norman, L., Thomas, D. N., and Bowie, A. R.: Biogeochemical observations during the winter-spring transition in East Antarctic sea ice: Evidence of iron and exopolysaccharide controls, Mar. Chem., 115, 163–175,, 2009. a

van Heuven, S. M., Hoppema, M., Jones, E. M., and de Baar, H. J.: Rapid invasion of anthropogenic CO2 into the deep circulation of the Weddell Gyre, Philos. T. Roy. Soc. A, 372, 20130056,, 2014. a

Verdy, A. and Mazloff, M. R.: A data assimilating model for estimating Southern Ocean biogeochemistry, J. Geophys. Res.-Oceans, 122, 6968–6988,, 2017. a

Verdy, A., Dutkiewicz, S., Follows, M. J., Marshall, J., and Czaja, A.: Carbon dioxide and oxygen fluxes in the Southern Ocean: Mechanisms of interannual variability, Global Biogeochem. Cy., 21, GB2020,, 2007. a

Viebahn, J. and Eden, C.: Towards the impact of eddies on the response of the Southern Ocean to climate change, Ocean Model., 34, 150–165,, 2010. a, b, c

Williams, G. D., Herraiz-Borreguero, L., Roquet, F., Tamura, T., Ohshima, K. I., Fukamachi, Y., Fraser, A. D., Gao, L., Chen, H., McMahon, C. R., Harcourt, R., and Hindell, M.: The suppression of Antarctic bottom water formation by melting ice shelves in Prydz Bay, Nat. Commun., 7, 12577,, 2016. a

Whithworth, T.: Monitoring the transport of the Antarctic Circumpolar Current at Drake Passage, J. Phys. Oceanogr., 13, 2045–2057,,2045:MTTOTA.2.0.CO;2, 1983. a

Yongliang, D., Hongwei, L., Weidong, Y., and Yijun, H.: The mean properties and variations of the Southerm Hemisphere subpolar gyres estimated by Simple Ocean Data Assimilation (SODA) products, Acta Oceanol. Sin., 35, 8–13,, 2016. a, b

Zickfeld, K., Fyfe, J. C., Saenko, O. A., Eby, M., and Weaver, A. J.: Response of the global carbon cycle to human-induced changes in Southern Hemisphere winds, Geophys. Res. Lett., 34, L12712,, 2007. a

Short summary
We present a new near-global coupled biogeochemical ocean-circulation model configuration of the Southern Ocean. The configuration features both a relatively equilibrated oceanic carbon inventory and an explicit representation of mesoscale eddies. In this paper, we document the model configuration and showcase its potential to tackle research questions such as the Southern Ocean carbon uptake dynamics on decadal timescales.