scholarly journals Diathermal Heat Transport in a Global Ocean Model

2019 ◽  
Vol 49 (1) ◽  
pp. 141-161 ◽  
Author(s):  
Ryan M. Holmes ◽  
Jan D. Zika ◽  
Matthew H. England

AbstractThe rate at which the ocean moves heat from the tropics toward the poles, and from the surface into the interior, depends on diabatic surface forcing and diffusive mixing. These diabatic processes can be isolated by analyzing heat transport in a temperature coordinate (the diathermal heat transport). This framework is applied to a global ocean sea ice model at two horizontal resolutions (1/4° and 1/10°) to evaluate the partioning of the diathermal heat transport between different mixing processes and their spatial and seasonal structure. The diathermal heat transport peaks around 22°C at 1.6 PW, similar to the peak meridional heat transport. Diffusive mixing transfers this heat from waters above 22°C, where surface forcing warms the tropical ocean, to temperatures below 22°C where midlatitude waters are cooled. In the control 1/4° simulation, half of the parameterized vertical mixing is achieved by background diffusion, to which sensitivity is explored. The remainder is associated with parameterizations for surface boundary layer, shear instability, and tidal mixing. Nearly half of the seasonal cycle in the peak vertical mixing heat flux is associated with shear instability in the tropical Pacific cold tongue, highlighting this region’s global importance. The framework presented also allows for quantification of numerical mixing associated with the model’s advection scheme. Numerical mixing has a substantial seasonal cycle and increases to compensate for reduced explicit vertical mixing. Finally, applied to Argo observations the diathermal framework reveals a heat content seasonal cycle consistent with the simulations. These results highlight the utility of the diathermal framework for understanding the role of diabatic processes in ocean circulation and climate.

Ocean Science ◽  
2017 ◽  
Vol 13 (2) ◽  
pp. 235-257 ◽  
Author(s):  
Reiner Onken

Abstract. The Regional Ocean Modeling System (ROMS) has been employed to explore the sensitivity of the forecast skill of mixed-layer properties to initial conditions, boundary conditions, and vertical mixing parameterisations. The initial and lateral boundary conditions were provided by the Mediterranean Forecasting System (MFS) or by the MERCATOR global ocean circulation model via one-way nesting; the initial conditions were additionally updated through the assimilation of observations. Nowcasts and forecasts from the weather forecast models COSMO-ME and COSMO-IT, partly melded with observations, served as surface boundary conditions. The vertical mixing was parameterised by the GLS (generic length scale) scheme Umlauf and Burchard (2003) in four different set-ups. All ROMS forecasts were validated against the observations which were taken during the REP14-MED survey to the west of Sardinia. Nesting ROMS in MERCATOR and updating the initial conditions through data assimilation provided the best agreement of the predicted mixed-layer properties with the time series from a moored thermistor chain. Further improvement was obtained by the usage of COSMO-ME atmospheric forcing, which was melded with real observations, and by the application of the k-ω vertical mixing scheme with increased vertical eddy diffusivity. The predicted temporal variability of the mixed-layer temperature was reasonably well correlated with the observed variability, while the modelled variability of the mixed-layer depth exhibited only agreement with the observations near the diurnal frequency peak. For the forecasted horizontal variability, reasonable agreement was found with observations from a ScanFish section, but only for the mesoscale wave number band; the observed sub-mesoscale variability was not reproduced by ROMS.


2020 ◽  
Author(s):  
Casimir de Lavergne ◽  
Clément Vic ◽  
Gurvan Madec ◽  
Fabien Roquet ◽  
Amy Waterhouse ◽  
...  

<p>Vertical mixing is often regarded as the Achilles’ heel of ocean models. In particular, few models include a comprehensive and energy-constrained parameterization of mixing by internal ocean tides. Here, we present an energy-conserving mixing scheme which accounts for the local breaking of high-mode internal tides and the distant dissipation of low-mode internal tides. The scheme relies on four static two-dimensional maps of internal tide dissipation, constructed using mode-by-mode Lagrangian tracking of energy beams from sources to sinks. Each map is associated with a distinct dissipative process and a corresponding vertical structure. Applied to an observational climatology of stratification, the scheme produces a global three-dimensional map of dissipation which compares well with available microstructure observations and with upper-ocean finestructure mixing estimates. Implemented in the NEMO global ocean model, the scheme improves the representation of deep water-mass transformation and obviates the need for a constant background diffusivity.</p>


2017 ◽  
Vol 47 (10) ◽  
pp. 2429-2454 ◽  
Author(s):  
Eric D. Skyllingstad ◽  
Jenessa Duncombe ◽  
Roger M. Samelson

AbstractGeneration of ocean surface boundary layer turbulence and coherent roll structures is examined in the context of wind-driven and geostrophic shear associated with horizontal density gradients using a large-eddy simulation model. Numerical experiments over a range of surface wind forcing and horizontal density gradient strengths, combined with linear stability analysis, indicate that the dominant instability mechanism supporting coherent roll development in these simulations is a mixed instability combining shear instability of the ageostrophic, wind-driven flow with symmetric instability of the frontal geostrophic shear. Disruption of geostrophic balance by vertical mixing induces an inertially rotating ageostrophic current, not forced directly by the wind, that initially strengthens the stratification, damps the instabilities, and reduces vertical mixing, but instability and mixing return when the inertial buoyancy advection reverses. The resulting rolls and instabilities are not aligned with the frontal zone, with an oblique orientation controlled by the Ekman-like instability. Mean turbulence is enhanced when the winds are destabilizing relative to the frontal orientation, but mean Ekman buoyancy advection is found to be relatively unimportant in these simulations. Instead, the mean turbulent kinetic energy balance is dominated by mechanical shear production that is enhanced when the wind-driven shear augments the geostrophic shear, while the resulting vertical mixing nearly eliminates any effective surface buoyancy flux from near-surface, cold-to-warm, Ekman buoyancy advection.


2016 ◽  
Author(s):  
Reiner Onken

Abstract. The Regional Ocean Modeling System (ROMS) has been employed to explore the sensitivity of the forecast skill of mixed-layer properties to the initial conditions, boundary conditions, and vertical mixing parameterisations. The initial and lateral boundary conditions were provided by the Mediterranean Forecasting System (MFS) or by the MERCATOR global ocean circulation model via one-way nesting; the initial conditions were additionally updated by the assimilation of observations. Nowcasts and forecasts from the weather forecast models COSMO-ME and COSMO-IT, partly melded with observations, served as surface boundary conditions. The vertical mixing was parameterised by the GLS (Generic Length Scale) scheme (Umlauf et al. 2003) in four different setups. All ROMS forecasts were validated against observations which were taken during the REP14-MED oceanographic survey to the west of Sardinia. Nesting ROMS in MERCATOR and updating the initial conditions by data assimilation provided the best agreement of the predicted mixed-layer temperature and the mixed-layer depth with time series from a moored thermistor chain. Further improvement was obtained by the usage of COSMO-ME atmospheric forcing which was melded with real observations, and by the application of the k − ε vertical mixing scheme with increased vertical eddy diffusivity. The predicted temporal variability of the mixed-layer temperature was reasonably well correlated with the observed variability in the frequency range above one cycle per day, while the modelled variability of the mixed-layer depth exhibited only agreement with the observations near the diurnal frequency peak. For the forecasted horizontal variability, reasonable agreement was found with observations from a ScanFish section, but only for the mesoscale wavenumber band; the observed sub-mesoscale variability was not reproduced by ROMS.


2013 ◽  
Vol 43 (4) ◽  
pp. 698-705 ◽  
Author(s):  
Julian J. Schanze ◽  
Raymond W. Schmitt

Abstract Owing to the larger thermal expansion coefficient at higher temperatures, more buoyancy is put into the ocean by heating than is removed by cooling at low temperatures. The authors show that, even with globally balanced thermal and haline surface forcing at the ocean surface, there is a negative density flux and hence a positive buoyancy flux. As shown by McDougall and Garrett, this must be compensated by interior densification on mixing due to the nonlinearity of the equation of state (cabbeling). Three issues that arise from this are addressed: the estimation of the annual input of density forcing, the effects of the seasonal cycle, and the total cabbeling potential of the ocean upon complete mixing. The annual expansion through surface density forcing in a steady-state ocean driven by balanced evaporation–precipitation–runoff (E–P–R) and net radiative budget at the surface Qnet is estimated as 74 000 m3 s−1 (0.07 Sv; 1 Sv ≡ 106 m3 s−1), which would be equivalent to a sea level rise of 6.3 mm yr−1. This is equivalent to approximately 3 times the estimated rate of sea level rise or 450% of the average Mississippi River discharge. When seasonal variations are included, this density forcing increases by 35% relative to the time-mean case to 101 000 m3 s−1 (0.1 Sv). Likely bounds are established on these numbers by using different Qnet and E–P–R datasets and the estimates are found to be robust to a factor of ~2. These values compare well with the cabbeling-induced contraction inferred from independent thermal dissipation rate estimates. The potential sea level decrease upon complete vertical mixing of the ocean is estimated as 230 mm. When horizontal mixing is included, the sea level drop is estimated as 300 mm.


2021 ◽  
Vol 2 (1) ◽  
Author(s):  
Carine G. van der Boog ◽  
Henk A. Dijkstra ◽  
Julie D. Pietrzak ◽  
Caroline A. Katsman

AbstractDouble-diffusive processes enhance diapycnal mixing of heat and salt in the open ocean. However, observationally based evidence of the effects of double-diffusive mixing on the global ocean circulation is lacking. Here we analyze the occurrence of double-diffusive thermohaline staircases in a dataset containing over 480,000 temperature and salinity profiles from Argo floats and Ice-Tethered Profilers. We show that about 14% of all profiles contains thermohaline staircases that appear clustered in specific regions, with one hitherto unknown cluster overlying the westward flowing waters of the Tasman Leakage. We estimate the combined contribution of double-diffusive fluxes in all thermohaline staircases to the global ocean’s mechanical energy budget as 7.5 GW [0.1 GW; 32.8 GW]. This is small compared to the estimated energy required to maintain the observed ocean stratification of roughly 2 TW. Nevertheless, we suggest that the regional effects, for example near Australia, could be pronounced.


2017 ◽  
Vol 837 ◽  
pp. 341-380 ◽  
Author(s):  
Peter P. Sullivan ◽  
James C. McWilliams

The evolution of upper ocean currents involves a set of complex, poorly understood interactions between submesoscale turbulence (e.g. density fronts and filaments and coherent vortices) and smaller-scale boundary-layer turbulence. Here we simulate the lifecycle of a cold (dense) filament undergoing frontogenesis in the presence of turbulence generated by surface stress and/or buoyancy loss. This phenomenon is examined in large-eddy simulations with resolved turbulent motions in large horizontal domains using${\sim}10^{10}$grid points. Steady winds are oriented in directions perpendicular or parallel to the filament axis. Due to turbulent vertical momentum mixing, cold filaments generate a potent two-celled secondary circulation in the cross-filament plane that is frontogenetic, sharpens the cross-filament buoyancy and horizontal velocity gradients and blocks Ekman buoyancy flux across the cold filament core towards the warm filament edge. Within less than a day, the frontogenesis is arrested at a small width,${\approx}100~\text{m}$, primarily by an enhancement of the turbulence through a small submesoscale, horizontal shear instability of the sharpened filament, followed by a subsequent slow decay of the filament by further turbulent mixing. The boundary-layer turbulence is inhomogeneous and non-stationary in relation to the evolving submesoscale currents and density stratification. The occurrence of frontogenesis and arrest are qualitatively similar with varying stress direction or with convective cooling, but the detailed evolution and flow structure differ among the cases. Thus submesoscale filament frontogenesis caused by boundary-layer turbulence, frontal arrest by frontal instability and frontal decay by forward energy cascade, and turbulent mixing are generic processes in the upper ocean.


2016 ◽  
Vol 46 (4) ◽  
pp. 1097-1115 ◽  
Author(s):  
Zhan Su ◽  
Andrew P. Ingersoll ◽  
Andrew L. Stewart ◽  
Andrew F. Thompson

AbstractThe energetics of thermobaricity- and cabbeling-powered deep convection occurring in oceans with cold freshwater overlying warm salty water are investigated here. These quasi-two-layer profiles are widely observed in wintertime polar oceans. The key diagnostic is the ocean convective available potential energy (OCAPE), a concept introduced in a companion piece to this paper (Part I). For an isolated ocean column, OCAPE arises from thermobaricity and is the maximum potential energy (PE) that can be converted into kinetic energy (KE) under adiabatic vertical parcel rearrangements. This study explores the KE budget of convection using two-dimensional numerical simulations and analytical estimates. The authors find that OCAPE is a principal source for KE. However, the complete conversion of OCAPE to KE is inhibited by diabatic processes. Further, this study finds that diabatic processes produce three other distinct contributions to the KE budget: (i) a sink of KE due to the reduction of stratification by vertical mixing, which raises water column’s center of mass and thus acts to convert KE to PE; (ii) a source of KE due to cabbeling-induced shrinking of the water column’s volume when water masses with different temperatures are mixed, which lowers the water column’s center of mass and thus acts to convert PE into KE; and (iii) a reduced production of KE due to diabatic energy conversion of the KE convertible part of the PE to the KE inconvertible part of the PE. Under some simplifying assumptions, the authors also propose a theory to estimate the maximum depth of convection from an energetic perspective. This study provides a potential basis for improving the convection parameterization in ocean models.


2021 ◽  
Author(s):  
Ryan Holmes ◽  
Jan Zika ◽  
Stephen Griffies ◽  
Andrew Hogg ◽  
Andrew Kiss ◽  
...  

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


Sign in / Sign up

Export Citation Format

Share Document