The geography of numerical mixing in a suite of global ocean models

2021 ◽  
Author(s):  
Ryan Holmes ◽  
Jan Zika ◽  
Stephen Griffies ◽  
Andrew Hogg ◽  
Andrew Kiss ◽  
...  
Keyword(s):  
Ocean Circulation ◽  
General Circulation ◽  
Horizontal Resolution ◽  
Global Ocean ◽  
Vertical Diffusivity ◽  
Mixing Parameters ◽  
Ocean Models ◽  
Numerical Mixing ◽  
Comparable Magnitude ◽  
Ice Model

<p>Numerical mixing, the physically spurious diffusion of tracers due to the numerical discretization of advection, is known to contribute to biases in ocean circulation models. However, quantifying numerical mixing is non-trivial, with most studies utilizing specifically targeted experiments in idealized settings. Here, we present a precise method based on water-mass transformation for quantifying numerical mixing, including its spatial structure, that can be applied to any conserved variable in global general circulation ocean models. The method is applied to a suite of global MOM5 ocean-sea ice model simulations with differing grid spacings and sub-grid scale parameterizations. In all configurations numerical mixing drives across-isotherm heat transport of comparable magnitude to that associated with explicitly-parameterized mixing. Numerical mixing is prominent at warm temperatures in the tropical thermocline, where it is sensitive to the vertical diffusivity and resolution. At colder temperatures, numerical mixing is sensitive to the presence of explicit neutral diffusion, suggesting that much of the numerical mixing in these regions acts as a proxy for neutral diffusion when it is explicitly absent. Comparison of equivalent (with respect to vertical resolution and explicit mixing parameters) 1/4-degree and 1/10-degree horizontal resolution configurations shows only a modest enhancement in numerical mixing at the eddy-permitting 1/4-degree resolution. Our results provide a detailed view of numerical mixing in ocean models and pave the way for future improvements in numerical methods.</p>

On the impact of sea ice in a global ocean circulation model

1997 ◽  
Vol 25 ◽  
pp. 111-115 ◽  
Author(s):  
Achim Stössel
Keyword(s):  
Southern Ocean ◽  
Sea Ice ◽  
Ocean Circulation ◽  
General Circulation ◽  
Thermohaline Circulation ◽  
Circulation Model ◽  
Deep Ocean ◽  
Global Ocean ◽  
Convective Activity ◽  
The Impact

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.

scholarly journals On the discretization of the ice thickness distribution in the NEMO3.6-LIM3 global ocean–sea ice model

2019 ◽  
Vol 12 (8) ◽  
pp. 3745-3758 ◽  
Author(s):  
François Massonnet ◽  
Antoine Barthélemy ◽  
Koffi Worou ◽  
Thierry Fichefet ◽  
Martin Vancoppenolle ◽  
...  
Keyword(s):  
Sea Ice ◽  
General Circulation ◽  
Large Scale ◽  
Thickness Distribution ◽  
Circulation Model ◽  
Grid Cell ◽  
The Arctic ◽  
Global Ocean ◽  
Ice Thickness ◽  
Ice Model

Abstract. The ice thickness distribution (ITD) is one of the core constituents of modern sea ice models. The ITD accounts for the unresolved spatial variability of sea ice thickness within each model grid cell. While there is a general consensus on the added physical realism brought by the ITD, how to discretize it remains an open question. Here, we use the ocean–sea ice general circulation model, Nucleus for European Modelling of the Ocean (NEMO) version 3.6 and Louvain-la-Neuve sea Ice Model (LIM) version 3 (NEMO3.6-LIM3), forced by atmospheric reanalyses to test how the ITD discretization (number of ice thickness categories, positions of the category boundaries) impacts the simulated mean Arctic and Antarctic sea ice states. We find that winter ice volumes in both hemispheres increase with the number of categories and attribute that increase to a net enhancement of basal ice growth rates. The range of simulated mean winter volumes in the various experiments amounts to ∼30 % and ∼10 % of the reference values (run with five categories) in the Arctic and Antarctic, respectively. This suggests that the way the ITD is discretized has a significant influence on the model mean state, all other things being equal. We also find that the existence of a thick category with lower bounds at ∼4 and ∼2 m for the Arctic and Antarctic, respectively, is a prerequisite for allowing the storage of deformed ice and therefore for fostering thermodynamic growth in thinner categories. Our analysis finally suggests that increasing the resolution of the ITD without changing the lower limit of the upper category results in small but not negligible variations of ice volume and extent. Our study proposes for the first time a bi-polar process-based explanation of the origin of mean sea ice state changes when the ITD discretization is modified. The sensitivity experiments conducted in this study, based on one model, emphasize that the choice of category positions, especially of thickest categories, has a primary influence on the simulated mean sea ice states while the number of categories and resolution have only a secondary influence. It is also found that the current default discretization of the NEMO3.6-LIM3 model is sufficient for large-scale present-day climate applications. In all cases, the role of the ITD discretization on the simulated mean sea ice state has to be appreciated relative to other influences (parameter uncertainty, forcing uncertainty, internal climate variability).

scholarly journals Simulating Southern Hemisphere extra-tropical climate variability with an idealised coupled atmosphere-ocean model

2012 ◽  
Vol 5 (5) ◽  
pp. 1161-1175 ◽  
Author(s):  
H. Kurzke ◽  
M. V. Kurgansky ◽  
K. Dethloff ◽  
D. Handorf ◽  
S. Erxleben ◽  
...  
Keyword(s):  
Climate Variability ◽  
Southern Hemisphere ◽  
Ocean Circulation ◽  
Coupled Model ◽  
Circulation Model ◽  
Horizontal Resolution ◽  
Global Ocean ◽  
Simplified Model ◽  
Cmip5 Models ◽  
Decadal Climate Variability

Abstract. A quasi-geostrophic model of Southern Hemisphere's wintertime atmospheric circulation with horizontal resolution T21 has been coupled to a global ocean circulation model with a resolution of 2° × 2° and simplified physics. This simplified coupled model reproduces qualitatively some features of the first and the second EOF of atmospheric 833 hPa geopotential height in accordance with NCEP data. The variability patterns of the simplified coupled model have been compared with variability patterns simulated by four complex state-of-the-art coupled CMIP5 models. The first EOF of the simplified model is too zonal and does not reproduce the right position of the centre of action over the Pacific Ocean and its extension to the tropics. The agreement in the second EOF between the simplified and the CMIP5 models is better. The total variance of the simplified model is weaker than the observational variance and those of the CMIP5 models. The transport properties of the Southern Ocean circulation are in qualitative accord with observations. The simplified model exhibits skill in reproducing essential features of decadal and multi-decadal climate variability in the extratropical Southern Hemisphere. Notably, 800 yr long coupled model simulations reveal sea surface temperature fluctuations on the timescale of several decades in the Antarctic Circumpolar Current region.

An Adjoint Analysis of the Meridional Overturning Circulation in a Hybrid Coupled Model

2006 ◽  
Vol 19 (15) ◽  
pp. 3751-3767 ◽  
Author(s):  
Véronique Bugnion ◽  
Chris Hill ◽  
Peter H. Stone
Keyword(s):  
Boundary Conditions ◽  
Wind Stress ◽  
Ocean Circulation ◽  
General Circulation ◽  
Meridional Overturning Circulation ◽  
Global Ocean ◽  
Model Parameters ◽  
Overturning Circulation ◽  
Near Surface ◽  
Meridional Overturning

Abstract Multicentury sensitivities in a realistic geometry global ocean general circulation model are analyzed using an adjoint technique. This paper takes advantage of the adjoint model’s ability to generate maps of the sensitivity of a diagnostic (i.e., the meridional overturning’s strength) to all model parameters. This property of adjoints is used to review several theories, which have been elaborated to explain the strength of the North Atlantic’s meridional overturning. This paper demonstrates the profound impact of boundary conditions in permitting or suppressing mechanisms within a realistic model of the contemporary ocean circulation. For example, the so-called Drake Passage Effect in which wind stress in the Southern Ocean acts as the main driver of the overturning’s strength, is shown to be an artifact of boundary conditions that restore the ocean’s surface temperature and salinity toward prescribed climatologies. Advective transports from the Indian and Pacific basins play an important role in setting the strength of the overturning circulation under “mixed” boundary conditions, in which a flux of freshwater is specified at the ocean’s surface. The most “realistic” regime couples an atmospheric energy and moisture balance model to the ocean. In this configuration, inspection of the global maps of sensitivity to wind stress and diapycnal mixing suggests a significant role for near-surface Ekman processes in the Tropics. Buoyancy also plays an important role in setting the overturning’s strength, through direct thermal forcing near the sites of convection, or through the advection of salinity anomalies in the Atlantic basin.

scholarly journals A multi-year timeseries of observation-based 3D horizontal and vertical quasi-geostrophic global ocean currents

2020 ◽  
Author(s):  
Bruno Buongiorno Nardelli
Keyword(s):  
Ocean Circulation ◽  
General Circulation ◽  
Numerical Models ◽  
Ocean Currents ◽  
Ocean Dynamics ◽  
Global Ocean ◽  
Diagnostic Model ◽  
Geostrophic Balance ◽  
Consiglio Nazionale Delle Ricerche ◽  
Marine Environment Monitoring

Abstract. Estimates of 3D ocean circulation are needed to improve our understanding of ocean dynamics and to assess its impact on marine ecosystems and Earth climate. Here we present the OMEGA3D product, an observation-based timeseries of (quasi) global 3D ocean currents covering the 1993–2018 period, developed by the Italian Consiglio Nazionale delle Ricerche within the European Copernicus Marine Environment Monitoring Service (CMEMS). This dataset was obtained by applying a diabatic quasi-geostrophic (QG) diagnostic model to CMEMS data-driven ARMOR3D weekly reconstruction of temperature and salinity and ERA-Interim fluxes. Outside the equatorial band, vertical velocities were retrieved in the upper 1500 m, at nominal ¼° resolution, and successively used to compute the horizontal ageostrophic components. Root mean square differences between OMEGA3D total horizontal velocities and totally independent drifter observations at two different depths (15 m and 1000 m) decrease with respect to corresponding estimates obtained from zero-order geostrophic balance, meaning that estimated vertical velocities can also be deemed reliable. OMEGA3D horizontal velocities are also closer to drifter observations than velocities provided by a set of re-analyses spanning a comparable time period, but based on data assimilation in ocean general circulation numerical models. The full OMEGA3D product (released on 31st of March 2020) is available upon free registration at https://doi.org/10.25423/cmcc/multiobs_glo_phy_w_rep_015_007. The reduced subset used here for validation and review purposes is openly available at https://doi.org/10.5281/zenodo.3696885 (Buongiorno Nardelli, 2020).

scholarly journals Global ocean storage of anthropogenic carbon

2012 ◽  
Vol 9 (7) ◽  
pp. 8931-8988 ◽  
Author(s):  
S. Khatiwala ◽  
T. Tanhua ◽  
S. Mikaloff Fletcher ◽  
M. Gerber ◽  
S. C. Doney ◽  
...  
Keyword(s):  
Ocean Circulation ◽  
Global Scale ◽  
Global Ocean ◽  
Observational Methods ◽  
Robust Estimates ◽  
Anthropogenic Carbon ◽  
Ocean Models ◽  
Biogeochemical Parameters ◽  
Range Of Values ◽  
Made In

Abstract. The global ocean is a significant sink for anthropogenic carbon (Cant), absorbing roughly a third of human CO2 emitted over the industrial period. Robust estimates of the magnitude and variability of the storage and distribution of Cant in the ocean are therefore important for understanding the human impact on climate. In this synthesis we review observational and model-based estimates of the storage and transport of Cant in the ocean. We pay particular attention to the uncertainties and potential biases inherent in different inference schemes. On a global scale, three data based estimates of the distribution and inventory of Cant are now available. While the inventories are found to agree within their uncertainty, there are considerable differences in the spatial distribution. We also present a review of the progress made in the application of inverse and data-assimilation techniques which combine ocean interior estimates of Cant with numerical ocean circulation models. Such methods are especially useful for estimating the air-sea flux and interior transport of Cant, quantities that are otherwise difficult to observe directly. However, the results are found to be highly dependent on modeled circulation, with the spread due to different ocean models at least as large as that from the different observational methods used to estimate Cant. Our review also highlights the importance of repeat measurements of hydrographic and biogeochemical parameters to estimate the storage of Cant on decadal timescales in the presence of the variability in circulation that is neglected by other approaches. Data-based Cant estimates provide important constraints on ocean forward models, which exhibit both broad similarities and regional errors relative to the observational fields. A compilation of inventories of Cant gives us a "best" estimate of the global ocean inventory of anthropogenic carbon in 2010 of 155 Pg C with an uncertainty of ±20%. This estimate includes a broad range of values suggesting that a combination of approaches is necessary in order to achieve a robust quantification of the ocean sink of anthropogenic CO2.

scholarly journals Ocean-sea-ice coupling in a global ocean general circulation model

1997 ◽  
Vol 25 ◽  
pp. 116-120 ◽  
Author(s):  
S. Legutke ◽  
E. Maier-Reimkr ◽  
A. Stössel ◽  
A. Hellbach
Keyword(s):  
Sea Ice ◽  
General Circulation Model ◽  
Ocean Circulation ◽  
General Circulation ◽  
Circulation Model ◽  
Ice Cover ◽  
Ocean General Circulation Model ◽  
Global Ocean ◽  
Ocean General Circulation ◽  
The Impact

A global ocean general circulation model has been coupled with a dynamic thermodynamic sea-ice model. This model has been spun-up in a 1000 year integration using daily atmosphere model data. Main water masses and currents are reproduced as well as the seasonal characteristics of the ice cover of the Northern and Southern Hemispheres. Model results for the Southern Ocean, however, show the ice cover as too thin, and there are large permanent polynyas in the Weddell and Ross Seas. These polynyas are due to a large upward oceanic heat flux caused by haline rejection during the freezing of sea ice. Sensitivity studies were performed to test several ways of treating the sea-surface salinity and the rejected brine. The impact on the ice cover, water-mass characteristics, and ocean circulation are described.

scholarly journals Preformed and regenerated phosphate in ocean general circulation models: can right total concentrations be wrong?

2012 ◽  
Vol 9 (5) ◽  
pp. 1797-1807 ◽  
Author(s):  
O. Duteil ◽  
W. Koeve ◽  
A. Oschlies ◽  
O. Aumont ◽  
D. Bianchi ◽  
...  
Keyword(s):  
Ocean Circulation ◽  
General Circulation ◽  
General Circulation Models ◽  
Global Ocean ◽  
Volume Distribution ◽  
Total Phosphate ◽  
Biogeochemical Models ◽  
Ocean General Circulation Models ◽  
Carbon Pump ◽  
Biological Carbon Pump

Abstract. Phosphate distributions simulated by seven state-of-the-art biogeochemical ocean circulation models are evaluated against observations of global ocean nutrient distributions. The biogeochemical models exhibit different structural complexities, ranging from simple nutrient-restoring to multi-nutrient NPZD type models. We evaluate the simulations using the observed volume distribution of phosphate. The errors in these simulated volume class distributions are significantly larger when preformed phosphate (or regenerated phosphate) rather than total phosphate is considered. Our analysis reveals that models can achieve similarly good fits to observed total phosphate distributions for a~very different partitioning into preformed and regenerated nutrient components. This has implications for the strength and potential climate sensitivity of the simulated biological carbon pump. We suggest complementing the use of total nutrient distributions for assessing model skill by an evaluation of the respective preformed and regenerated nutrient components.

Scattering of internal tide energy to super-tidal frequencies in global HYCOM 

2021 ◽  
Author(s):  
Miguel Solano ◽  
Maarten Buijsman
Keyword(s):  
Internal Tide ◽  
Tidal Energy ◽  
Horizontal Resolution ◽  
Internal Tides ◽  
Numerical Dispersion ◽  
Global Ocean ◽  
Wave Interactions ◽  
Continental Shelves ◽  
Lee Wave ◽  
Ocean Models

<p>Energy decay in realistically forced global ocean models has been mostly studied in the diurnal and semi-diurnal tidal bands and it is unclear how much of the tidal energy in these bands is scattered to higher frequencies. Global ocean models and satellite altimetry have shown that low-mode internal tides can propagate thousands of kilometers from their generation sites before being dissipated in the ocean interior but their pathway to dissipation is obscured due to lee-wave breaking at generation, wave-wave interactions, topographic scattering, shearing instabilities and shoaling on continental shelves. Internal tides from some generation sites, such as the Amazon shelf and the Nicobar and Andaman island chain, have large amounts of energy resulting in a steepening of the internal waves into solitary wave trains due to non-hydrostatic dispersion. In HYCOM, a hydrostatic model, this process is partially simulated by numerical dispersion. However, it is yet unknown how the dissipation of internal tides is affected by the numerical dispersion in hydrostatic models. In this study we use the method of vertical modes and rotary spectra to quantify the scattering of internal tides to higher-frequencies and analyze the dissipation processes in global HYCOM simulations with 4-km horizontal resolution.</p>

scholarly journals A nested Atlantic-Mediterranean Sea general circulation model for operational forecasting

Ocean Science ◽  
2009 ◽  
Vol 5 (4) ◽  
pp. 461-473 ◽  
Author(s):  
P. Oddo ◽  
M. Adani ◽  
N. Pinardi ◽  
C. Fratianni ◽  
M. Tonani ◽  
...  
Keyword(s):  
Mediterranean Sea ◽  
General Circulation Model ◽  
Ocean Circulation ◽  
General Circulation ◽  
Circulation Model ◽  
Ocean Model ◽  
Global Ocean ◽  
Lateral Boundary ◽  
The Mediterranean ◽  
The Impact

Abstract. A new numerical general circulation ocean model for the Mediterranean Sea has been implemented nested within an Atlantic general circulation model within the framework of the Marine Environment and Security for the European Area project (MERSEA, Desaubies, 2006). A 4-year twin experiment was carried out from January 2004 to December 2007 with two different models to evaluate the impact on the Mediterranean Sea circulation of open lateral boundary conditions in the Atlantic Ocean. One model considers a closed lateral boundary in a large Atlantic box and the other is nested in the same box in a global ocean circulation model. Impact was observed comparing the two simulations with independent observations: ARGO for temperature and salinity profiles and tide gauges and along-track satellite observations for the sea surface height. The improvement in the nested Atlantic-Mediterranean model with respect to the closed one is particularly evident in the salinity characteristics of the Modified Atlantic Water and in the Mediterranean sea level seasonal variability.

Sign in / Sign up

Export Citation Format

Share Document