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

2020 ◽  
Vol 13 (1) ◽  
pp. 71-97
Author(s):  
Heiner Dietze ◽  
Ulrike Löptien ◽  
Julia Getzlaff

Abstract. 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.

2019 ◽  
Author(s):  
Heiner Dietze ◽  
Ulrike Löptien ◽  
Julia Getzlaff

Abstract. 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 GFDL's Modular Ocean Model coupled to the Biology Light Iron Nutrients and Gasses (BLING) ecosystem model of Galbraith et al. (2010). The configuration is cutting-edge in that it features both a relatively equilibrated oceanic carbon inventory and a realistic representation of eddy kinetic energy – a combination that has, to-date, been precluded by prohibitive computational cost. Results from a simulation with climatological forcing and a sensitivity experiment with increasing winds suggest that the configuration is suited to explore Southern Ocean Carbon uptake dynamics on decadal timescales. Further, the fidelity of simulated bottom water temperatures off and on the Antarctic Shelf suggest that the configuration may be used to provide boundary conditions to ice-sheet models. The configuration is dubbed MOMSO a Modular Ocean Model Southern Ocean configuration.


1997 ◽  
Vol 25 ◽  
pp. 111-115 ◽  
Author(s):  
Achim Stössel

This paper investigates the long-term impact of sea ice on global climate using a global sea-ice–ocean general circulation model (OGCM). The sea-ice component involves state-of-the-art dynamics; the ocean component consists of a 3.5° × 3.5° × 11 layer primitive-equation model. Depending on the physical description of sea ice, significant changes are detected in the convective activity, in the hydrographic properties and in the thermohaline circulation of the ocean model. Most of these changes originate in the Southern Ocean, emphasizing the crucial role of sea ice in this marginally stably stratified region of the world's oceans. Specifically, if the effect of brine release is neglected, the deep layers of the Southern Ocean warm up considerably; this is associated with a weakening of the Southern Hemisphere overturning cell. The removal of the commonly used “salinity enhancement” leads to a similar effect. The deep-ocean salinity is almost unaffected in both experiments. Introducing explicit new-ice thickness growth in partially ice-covered gridcells leads to a substantial increase in convective activity, especially in the Southern Ocean, with a concomitant significant cooling and salinification of the deep ocean. Possible mechanisms for the resulting interactions between sea-ice processes and deep-ocean characteristics are suggested.


2006 ◽  
Vol 3 (3) ◽  
pp. 637-669 ◽  
Author(s):  
S. Natale ◽  
R. Sorgente ◽  
S. Gaberšek ◽  
A. Ribotti ◽  
A. Olita

Abstract. Ocean forecasts over the Central Mediterranean, produced by a near real time regional scale system, have been evaluated in order to assess their predictability. The ocean circulation model has been forced at the surface by a medium, high or very high resolution atmospheric forcing. The simulated ocean parameters have been compared with satellite data and they were found to be generally in good agreement. High and very high resolution atmospheric forcings have been able to form noticeable, although short-lived, surface current structures, due to their ability to detect transient atmospheric disturbances. The existence of the current structures has not been directly assessed due to lack of measurements. The ocean model in the slave mode was not able to develop dynamics different from the driving coarse resolution model which provides the boundary conditions.


2020 ◽  
Vol 13 (11) ◽  
pp. 5465-5483
Author(s):  
Clément Bricaud ◽  
Julien Le Sommer ◽  
Gurvan Madec ◽  
Christophe Calone ◽  
Julie Deshayes ◽  
...  

Abstract. Ocean biogeochemical models are key tools for both scientific and operational applications. Nevertheless the cost of these models is often expensive because of the large number of biogeochemical tracers. This has motivated the development of multi-grid approaches where ocean dynamics and tracer transport are computed on grids of different spatial resolution. However, existing multi-grid approaches to tracer transport in ocean modelling do not allow the computation of ocean dynamics and tracer transport simultaneously. This paper describes a new multi-grid approach developed for accelerating the computation of passive tracer transport in the Nucleus for European Modelling of the Ocean (NEMO) ocean circulation model. In practice, passive tracer transport is computed at runtime on a grid with coarser spatial resolution than the hydrodynamics, which reduces the CPU cost of computing the evolution of tracers. We describe the multi-grid algorithm, its practical implementation in the NEMO ocean model, and discuss its performance on the basis of a series of sensitivity experiments with global ocean model configurations. Our experiments confirm that the spatial resolution of hydrodynamical fields can be coarsened by a factor of 3 in both horizontal directions without significantly affecting the resolved passive tracer fields. Overall, the proposed algorithm yields a reduction by a factor of 7 of the overhead associated with running a full biogeochemical model like PISCES (with 24 passive tracers). Propositions for further reducing this cost without affecting the resolved solution are discussed.


2021 ◽  
pp. 50-66
Author(s):  
V. N. Stepanov ◽  
◽  
Yu. D. Resnyanskii ◽  
B. S. Strukov ◽  
A. A. Zelen’ko ◽  
...  

The quality of simulation of model fields is analyzed depending on the assimilation of various types of data using the PDAF software product assimilating synthetic data into the NEMO global ocean model. Several numerical experiments are performed to simulate the ocean–sea ice system. Initially, free model was run with different values of the coefficients of horizontal turbulent viscosity and diffusion, but with the same atmospheric forcing. The model output obtained with higher values of these coefficients was used to determine the first guess fields in subsequent experiments with data assimilation, while the model results with lower values of the coefficients were assumed to be true states, and a part of these results was used as synthetic observations. The results are analyzed that are assimilation of various types of observational data using the Kalman filter included through the PDAF to the NEMO model with real bottom topography. It is shown that a degree of improving model fields in the process of data assimilation is highly dependent on the structure of data at the input of the assimilation procedure.


2018 ◽  
Vol 48 (5) ◽  
pp. 1139-1150 ◽  
Author(s):  
Lachlan Stoney ◽  
Kevin J. E. Walsh ◽  
Steven Thomas ◽  
Paul Spence ◽  
Alexander V. Babanin

Abstract A parameterization of turbulent mixing from unbroken surface waves is included in a 16-yr simulation within a high-resolution ocean circulation model (MOM5). This “surface wave mixing” (SWM) derives from the wave orbital motion and is parameterized as an additional term in a k-epsilon model. We show that SWM leads to significant changes in sea surface temperatures but smaller changes in ocean heat content, and show the extent to which these changes can reduce pre-existing model biases with respect to observed data. Specifically, SWM leads to a widespread improvement in sea surface temperature in both hemispheres in summer and winter, while for ocean heat content the improvements are less clear. In addition, we show that introducing SWM can lead to an accumulation of wave-induced ocean heat content between years. While it has been well established that secular positive trends exist in global wave heights, we find that such trends are relatively unimportant in driving the accumulation of wave-induced ocean heat content. Rather, in response to the new source of mixing, the simulated ocean climate evolves toward a new equilibrium with greater total ocean heat content.


2016 ◽  
Vol 46 (5) ◽  
pp. 1399-1419 ◽  
Author(s):  
Maarten C. Buijsman ◽  
Joseph K. Ansong ◽  
Brian K. Arbic ◽  
James G. Richman ◽  
Jay F. Shriver ◽  
...  

AbstractThe effects of a parameterized linear internal wave drag on the semidiurnal barotropic and baroclinic energetics of a realistically forced, three-dimensional global ocean model are analyzed. Although the main purpose of the parameterization is to improve the surface tides, it also influences the internal tides. The relatively coarse resolution of the model of ~8 km only permits the generation and propagation of the first three vertical modes. Hence, this wave drag parameterization represents the energy conversion to and the subsequent breaking of the unresolved high modes. The total tidal energy input and the spatial distribution of the barotropic energy loss agree with the Ocean Topography Experiment (TOPEX)/Poseidon (TPXO) tidal inversion model. The wave drag overestimates the high-mode conversion at ocean ridges as measured against regional high-resolution models. The wave drag also damps the low-mode internal tides as they propagate away from their generation sites. Hence, it can be considered a scattering parameterization, causing more than 50% of the deep-water dissipation of the internal tides. In the near field, most of the baroclinic dissipation is attributed to viscous and numerical dissipation. The far-field decay of the simulated internal tides is in agreement with satellite altimetry and falls within the broad range of Argo-inferred dissipation rates. In the simulation, about 12% of the semidiurnal internal tide energy generated in deep water reaches the continental margins.


2019 ◽  
Author(s):  
Svetlana N. Losa ◽  
Stephanie Dutkiewicz ◽  
Martin Losch ◽  
Julia Oelker ◽  
Mariana A. Soppa ◽  
...  

Abstract. This study highlights recent advances and challenges of applying coupled physical-biogeochemical modeling for investigating the distribution of the key phytoplankton groups in the Southern Ocean, an area of strong interest for understanding biogeochemical cycling and ecosystem functioning under present climate change. Our simulations of the phenology of various Phytoplankton Functional Types (PFTs) are based on a version of the Darwin biogeochemical model coupled to the Massachusetts Institute of Technology (MIT) general circulation model (Darwin-MITgcm). The ecological module version was adapted for the Southern Ocean by: 1) improving coccolithophores abundance relative to the original model by introducing a high affinity for nutrients and an ability to escape grazing control for coccolithophores; 2) including two different (small vs. large) size classes of diatoms; and 3) accounting for two distinct life stages for Phaeocystis (single cell vs. colonial). This new model configuration describes best the competition and co-occurrence of the PFTs in the Southern Ocean. It improves significantly relative to an older version the agreement of the simulated abundance of the coccolithophores and diatoms with in situ scanning electron microscopy observations in the Subantarctic Zone as well as with in situ diatoms and haptophytes (including coccolithophores and Phaeocystis) chlorophyll a concentrations within the Patagonian Shelf and along the Western Antarctic Peninsula obtained by diagnostic pigment analysis. The modeled Southern Ocean PFT dominance also agrees well with satellite-based PFT information.


2016 ◽  
Author(s):  
Sergey Danilov ◽  
Dmitry Sidorenko ◽  
Qiang Wang ◽  
Thomas Jung

Abstract. Version 2 of the unstructured-mesh sea ice – ocean circulation model FESOM is presented. It builds upon FESOM1.4 (Wang et al., 2014, Geosci. Mod. Dev., 7, 663–693) but differs by its dynamical core (finite volumes instead of finite elements) and is formulated using the Arbitrary Lagrangian Eulerian (ALE) vertical coordinate, which increases model flexibility. The model inherits the framework and sea ice model from the previous version, which minimizes the efforts needed from a user to switch from one version to the other. The ocean states simulated with FESOM1.4 and FESOM2.0 driven by CORE-II forcing are compared on a mesh used for CORE-II intercomparison project. Additionally the performance on an eddy-permitting mesh with uniform resolution is discussed. The new version improves numerical efficiency of FESOM in terms of CPU time by at least three times while retaining its fidelity in simulating sea ice and ocean. From this it is argued that FESOM2.0 provides a major step forward in establishing unstructured-mesh models as valuable tools in climate research.


2007 ◽  
Vol 14 (6) ◽  
pp. 777-788
Author(s):  
A. D. Terwisscha van Scheltinga ◽  
H. A. Dijkstra

Abstract. We propose an efficient method for estimating a time-mean state of an ocean model subject to given observations using implicit time-stepping. The new method uses (i) an implicit implementation of the 4D-Var method to fit the model trajectory to the observations, and (ii) a pre-processor which applies a multi-channel singular spectrum analysis to enhance the signal-to-noise ratio of the observational data and to filter out the high frequency variability. This approach enables one to estimate the time-mean model state using larger time-steps than is possible with an explicit model. The performance of the method is presented for two test cases within a barotropic quasi-geostrophic nonlinear model of the wind-driven double-gyre ocean circulation. The method turns out to be accurate and, in comparison with the time-mean state computed with an explicit version of the model, relatively cheap in computational cost.


Sign in / Sign up

Export Citation Format

Share Document