The barotropic vorticity (BV) balance is fundamental when interpreting the ocean gyre circulation. Here we propose an intercomparison of vorticity equations for the depth‐integrated flow applied to ...ocean models. We review four distinct variants of the BV balances, each giving access to diagnostic equations for the depth‐integrated ocean circulation, either meridional, across geostrophic contours or its divergence. We then formulate those balances in the Vorticity Balances in NEMO (VoBiN) diagnostic package aimed at the NEMO ocean platform and more generally C‐grid ocean models. We show that spatial discretization of the equations of motion have profound implications for those vorticity balances. Finally, we diagnose the main balances of a global ocean climate simulation. In all vorticity balances, topographic torques arise from interactions of the flow with slanting topography. We identify significant spurious topographic torques related to the model's C‐grid discretizations, and we suggest ways to address them. In the depth‐integrated and BV balances, bottom vortex stretching and bottom pressure torque drive the flow interaction with topography, respectively. Contrary to Sverdrup theory, the wind stress curl, although dominant in the interior Subtropics, becomes a minor player anywhere significant bottom velocities prevail. The geostrophic contour vorticity balance highlights the limits of barotropic models of the ocean circulation through the so‐called JEBAR term. Finally, the transport divergence vorticity balance stresses the limitations of Ekman plus geostrophic dynamics for the mass balance closure in ocean models. This framework should encourage ocean modellers to diagnose more routinely momentum and vorticity equations.
Plain Language Summary
Ocean gyre theories involve the key role played by the wind variations across latitudes to force an interior flow. However, recent work has put forward the role played by bottom topography as a guide, or an obstacle, to the gyre circulation. The general framework employed in those theories is the so‐called barotropic vorticity equation involving the balancing of the spin induced by the Earth rotation in an ocean in motion. This work proposes a review on oceanic vorticity balances that synthesizes the informations that they provide about the ocean circulation. We then apply this framework to an ocean model used for climate projections. Our analysis confirms the key role played by interactions with topography in driving the gyre circulation. The ocean bottom topography is a geological constant, so that it constitutes a long‐lasting constraint for the circulation. Finally, we also stress the large effect of the ocean circulation formulation in the computer programs that make up an ocean model and should be designed carefully.
Key Points
The large‐scale ocean circulation includes four distinct vorticity balances for the depth‐integrated flow
Those vorticity balances are discretized in the VoBiN diagnostic module designed for C‐grid ocean models such as the NEMO platform
The barotropic vorticity equation of a global NEMO climate simulation is dominated by Sverdrup and topographic balances
The Strait of Gibraltar is a narrow and shallow channel that controls the Mediterranean Sea thermohaline and biogeochemical balances. Strong tidal currents significantly modulate exchanges across ...this strait and induce an intense vertical mixing, impacting both the Mediterranean Sea and the Atlantic Ocean on a climatic scale. However, the turbulent processes controlling the tidal mixing location, timing, and magnitude remain unclear. To fill this gap, we investigate tidal mixing at the Strait of Gibraltar in yearly twin tidal and non-tidal simulations from a regional configuration of the three-dimensional numerical model MITgcm, using a high spatial resolution around the Strait of Gibraltar (1/200°, 100 vertical levels). More specifically, we investigate the model turbulence closure scheme, based on a turbulent kinetic energy budget, and illustrate that vertical buoyancy fluxes should be preferred to diapycnal diffusivities as mixing indicators. In agreement with previous literature, we find that tides strongly intensify vertical mixing and motions within the Strait of Gibraltar. We then demonstrate that tidal mixing relies on two main ingredients: a sustained vertical shear of horizontal velocities and a local weakening of stratification. In the Mediterranean layer, the former drives diapycnal mixing near the seafloor and the latter in shallower areas above the prominent sills of the strait. We also evidence the frequent but irregular occurrences of static instabilities in the vicinity of these sills. In the Atlantic layer, both vertical shear and stratification are involved in diapycnal mixing that develops along the trail of the eastward internal bore released at the Camarinal sill. At high frequency, the local weakening of stratification results from convergence and divergence patterns in the vicinity of the Camarinal and Espartel sills, feeding recirculation cells between the Atlantic and Mediterranean layers. In addition, we highlight that diapycnal mixing mainly develops during the westward tidal phase in the Mediterranean layer and the eastward tidal phase in the Atlantic layer. We conclude by proposing a revised conceptual view of tidal mixing at the Strait of Gibraltar, where tidally-induced recirculation cells play an instrumental role in transforming the exchanged water masses. Overall, this study emphasizes the relevance of a realistic representation of both tides and abrupt topography to simulate the exchanges through the Strait of Gibraltar and argues for the use of a specific tidal mixing parameterization otherwise.
•Tidal mixing mechanisms are investigated in a highresolution simulation of the Strait of Gibraltar.•Tidal convergence and divergence patterns, that feed recirculation cells, play an instrumental role in enhanced mixing at the Strait of Gibraltar.•Tides drive diapycnal mixing by enhancing vertical shear and weakening vertical stratification.•Accurate representation of both tides and abrupt topography is key for representing exchanges through the Strait of Gibraltar.
The interannual variability of the Atlantic cold tongue (ACT) is studied by means of a mixed-layer heat budget analysis. A method to classify extreme cold and warm ACT events is proposed and applied ...to ten various analysis and reanalysis products. This classification allows 5 cold and 5 warm ACT events to be selected over the period 1982–2007. Cold (warm) ACT events are defined by the presence of negative (positive) sea surface temperature (SST) anomalies at the center of the equatorial Atlantic in late boreal spring, preceded by negative (positive) zonal wind stress anomalies in the western equatorial Atlantic. An ocean general circulation model capable of reconstructing the interannual variability of the ACT correctly is used to demonstrate that cold ACT events develop rapidly from May to June mainly due to intense cooling by vertical mixing and horizontal advection. The simulated cooling at the center of the basin is the result of the combined effects of non-local and local processes. The non-local process is an upwelling associated with an eastward-propagating Kelvin wave, which makes the mixed-layer more shallow and preconditions the upper layers to be cooled by an intense heat loss at the base of the mixed-layer, which is amplified by a stronger local injection of energy from the atmosphere. The early cooling by vertical mixing in March is also shown to be a good predictor of June cooling. In July, horizontal advection starts to warm the mixed-layer abnormally and damps SST anomalies. The advection anomalies, which result from changes in the horizontal temperature gradient, are associated in some cases with the propagation of Rossby waves along the equator. During warm ACT events, processes are reversed, generating positive SST anomalies: a downwelling Kelvin wave triggers stratification anomalies and mixed-layer depth anomalies, amplified by a weaker injection of energy from the atmosphere in May–June. In July, warm ACT events are abnormally cooled due to negative horizontal advection anomalies resulting from processes similar to those that occur during cold ACT events. This additional cooling process extends the period of cooling of the ACT, reducing SST anomalies.
The formation of the Atlantic cold tongue (ACT) is the dominant seasonal sea surface temperature signal in the eastern equatorial Atlantic (EEA). A comprehensive analysis of variability in its ...spatial extent, temperature, and onset is presented. Then, the physical mechanisms which initiate ACT onset, as well as the feedbacks from the ACT to the maritime boundary layer, and how the ACT influences the onset of the West African monsoon (WAM) are discussed. We argue that in the EEA, the air‐sea coupling between the ACT and WAM occurs in two phases. From March to mid‐June, the ACT results from the intensification of the southeastern trades associated with the St. Helena anticyclone. Steering of surface winds by the basin shape of the EEA imparts optimal wind stress for generating the maximum upwelling south of the equator. During the second phase (mid‐June–August), wind speeds north of the equator increase as a result of the northward progression of the intensifying trades and as a result of significant surface heat flux gradients produced by the differential cooling between the ACT and the tropical waters circulating in the Gulf of Guinea (GG). It is anticipated that the atmospheric divergence induced at low levels north of the equator reduces convection over the GG and that increased northward winds shift convection over land. Correlations between the ACT and the WAM onset dates over the last 26 years (1982–2007) measure as much as 0.8. This suggests that the ACT plays a key role in the WAM onset.
Key Points
Coupling between the Atlantic cold tongue and the West African monsoon
Variability of the Atlantic cold tongue
Interaction between the SST, winds, and precipitation
Estimating the mixed‐layer heat budget is a key issue for understanding the cold tongue development in the eastern equatorial Atlantic. A high‐resolution ocean regional model is used to diagnose the ...mixed‐layer heat budget online during the EGEE‐3 experiment from May to August 2006. The heat budget shows the major role of the horizontal advection and turbulent mixing in the mixed‐layer temperature balance in the cold tongue. The surface net heat flux and entrainment processes play a minor role. The equatorial cooling is mainly induced by low‐frequency advection, which is balanced by high‐frequency zonal and meridional advections. The high‐frequency advections are organized in patterns along the northern edge of the cold tongue, where they are associated with strong sea surface temperature gradients and well‐developed tropical instability waves in the western Atlantic. Special attention is paid to the wind energy flux, which controls horizontal advection and turbulent mixing. We suggest that the wind energy flux drives the vertical velocity, which in turn adjusts the mixed‐layer depth, its stratification, and the vertical shear of the horizontal current. Although vertical advection is not essential in providing cold water in the Atlantic cold tongue, it is shown that the vertical velocity plays a central role in preconditioning the mixed layer and maximizes the turbulent mixing.
Key Points
Horizontal advection and turbulent mixing drive the cold tongue development
Low/high‐frequency advection cools/warms the Equatorial Atlantic
Surface wind‐energy flux controls the equatorial mixed‐layer dynamics
Estimating the processes that control the north equatorial sea surface temperature (SST)-front on the northern edge of the cold tongue in the tropical Atlantic is a key issue for understanding the ...dynamics of the oceanic equatorial Atlantic and the West African Monsoon. Diagnosis of the frontogenetic forcings on a realistic high-resolution simulation was used to identify the processes involved in the formation and evolution of the equatorial SST-front. The turbulent forcing associated with the mixed-layer turbulent heat flux was found to be systematically frontolytic while the dynamic forcing associated with currents was found to be frontogenetic for the equatorial SST-front. Nevertheless, the low-frequency component of the turbulent forcing was frontogenetic and initiated the SST-front which was then amplified and maintained by the leading dynamic forcing. This forcing was mainly driven by the meridional convergence of the northern South Equatorial Current (nSEC) and the Guinea Current, which points out the essential role played by the circulation in the equatorial SST-front evolution. The quasi-biweekly variability of the equatorial SST-front and its forcings were found to be more strongly coupled to the wind energy flux (WEF) than to the surface wind stress. In fact the WEF controlled the convergence/divergence of the nSEC and Guinea Current and thus the meridional component of the leading dynamic forcing. The WEF explains the equatorial SST-front development better than the wind does because it is a coupled ocean-atmosphere process.
A new one‐dimensional (1‐D) parameterization of penetrative convection has been developed in order to have a better representation of the vertical mixing in ocean general circulation models. Our ...approach is inspired from atmospheric parameterizations of shallow convection which assumes that in the convective boundary layer, the subgrid‐scale fluxes result from two different mixing scales: small eddies, which are represented by an Eddy‐Diffusivity (ED) contribution, and large eddies associated with thermals, which are represented by a mass‐flux contribution. In the present work, the local (small eddies) and nonlocal (large eddies) contributions are unified into an Eddy‐Diffusivity‐Mass‐Flux (EDMF) parameterization which treats simultaneously the whole vertical mixing. EDMF is implemented in the community ocean model NEMO and tested in its 1‐D column version. Deepening of dense water in analytic cases, successfully reproduced in LES simulations, is more realistic with EDMF than with standard diffusion parameterizations. Also the convective events observed in the western Mediterranean at the Lion station and in the North Pacific Ocean at the PAPA station are more realistic in terms of sequencing and amplitude with EDMF.
Plain Language Summary
A new representation of oceanic convection has been developed in order to have a better representation of vertical mixing in ocean models. Our proposition is to represent oceanic convection consistently to atmospheric convection. We want to represent the convective plumes in oceans in the same way that the cumulus clouds are represented in the atmosphere. In this unified approach, the oceanic vertical mixing is viewed as a combination of large eddies associated with strong nonbuoyant downdrafts and small eddies which induce local turbulence. This new paradigm of oceanic mixing leads to more realistic simulations of the hydrological properties of water masses. In the future, it is expected to obtain more reliable climate projections.
Key Points
An Eddy‐Diffusivity Mass‐Flux (EDMF) parameterization is proposed to represent oceanic convection consistently to atmospheric convection
The EDMF parameterization unifies diffusion and convection processes in ocean models
The EDMF parameterization represents the penetrative convection and the associated counter‐gradient heat flux in the stratified thermocline
A generalized ω‐equation is used to identify the contributions from different processes that force upward motions in the Atlantic Marine ITCZ (AMI) from a numerical mesoscale simulation of June 2010. ...This ω‐equation separates the diabatic heating contributions, which lie at the core of the Weak Temperature Gradient (WTG) framework, from the dynamical terms. Three layers of atmosphere are found with different balance. In the Marine Atmospheric Boundary‐Layer (MABL), the upward motions in the AMI are induced by the frontogenesis and buoyancy components, which are regulated by the ageostrophic adjustment due to the presence of thermal‐wind imbalance. The balance of these three processes well captures the variability of the vertical velocity and the associated precipitation, meaning that boundary‐layer processes play a central role in the AMI dynamics. In the layer 600–2,000 m, a zone of strong vertical wind‐shear just above the MABL, the upward motions are induced by the ageostrophic adjustment and radiative components, which are counteracted by evaporation of convective precipitation. Above 2,000 m the ascending motions are driven by the deep convection heating, as expected by the WTG framework, and more surprisingly by the ageostrophic adjustment term within the Tropical Easterly Jet. Thanks to the use of the ω‐equation, these results extend the current WTG framework to the boundary layer, where it is not expected to hold. In the free troposphere, the WTG framework only accounts for half of the AMI ascent, the other half being forced by the dynamical terms.
Plain Language Summary
Upward motions in the Atlantic Marine ITCZ are generally explained as resulting from latent heat release associated with deep convection in the troposphere. This study shows that this approach is challenged by other processes on the vertical. In the lowest atmospheric layers, upward motions are driven by frontogenetic processes induced by wind convergence and by differential heating associated with the turbulent heat fluxes. In the middle troposphere, upward motions are controlled by ageostrophic adjustment, characterizing an unbalanced system, and radiative heating. In the deep troposphere, ascending motions are controlled by heating induced by the convective latent heat release, as expected by classic theories, but more surprisingly also by ageostrophic adjustment within the Tropical Easterly Jet. These results shows that dynamic forcings have an equivalent role to diabatic forcings in the production of upward motions within the ITCZ.
Key Points
Upward motions in the Atlantic Marine ITCZ are not resulting from latent heat release in the deep troposphere only
Buoyancy fluxes, frontogenesis, and ageostrophic circulation are also key factors that force the ascent from the surface to the top
The surface buoyancy flux is found to strongly influence the vertical velocity in the marine atmospheric boundary‐layer, unlike the sea surface temperature
The correct representation of the Maud Rise open‐ocean polynya in the Weddell Sea remains a challenge for ocean models. Here we reproduce the most recent polynya openings in 2016–2017 using a ...regional configuration, and assess their dependencies on vertical convective mixing schemes and freshwater forcing, both separately and in combination. We test three vertical convective mixing schemes: the enhanced vertical diffusion (EVD), the Eddy‐Diffusivity Mass‐Flux (EDMF) parameterization, and a modified version of EDMF accounting for thermobaric effects. Using simulations for the period 2007–2017, we find that the modified EDMF reproduces the observed climatological evolution of the mixed layer depth better than the original EDMF and the EVD, but a polynya fails to open due to excessive freshwater forcing. We thus use the modified EDMF to perform sensitivity experiments with reduced precipitation during 2012–2017. The imposed freshwater forcing strongly affects the number of years with polynyas. The simulation with the best representation of the 2016–2017 polynyas is analyzed to evaluate the triggering mechanisms. The 2016 polynya was induced by the action of thermobaric instabilities on a weak ambient stratification. This opening preconditioned the water column for 2017, which produced a stronger polynya. By examining the impacts of the different convective mixing schemes, we show that the modified EDMF generates more realistic patterns of deep convection. Our results highlight the importance of surface freshwater forcing and thermobaricity in governing deep convection around Maud Rise, and the need to represent thermobaric instabilities to realistically model Maud Rise polynyas.
Plain Language Summary
We investigate the impacts of representing numerical vertical mixing and surface freshwater forcing in a regional ocean model on polynyas (large openings in the pack ice) at Maud Rise, Southern Ocean. Maud Rise is prone to hosting polynyas, often associated with deep convection, which is a local vertical mixing process homogenizing the water column between surface and depths of several hundred meters. Numerical models often use simplistic strategies to represent this process, but improved parameterizations have recently become available. In this work, we test the impact of the representation of convective mixing in a particularly sensitive region. The last Maud Rise polynyas were observed in 2016 and 2017. Our regional simulation is capable of reproducing these polynyas, which has long been a challenge for ocean‐sea ice models. We show that the 2016 polynya resulted from the action of a vertical instability at depth acting on weak ambient stratification. This event preconditioned the stronger 2017 polynya and deep convection. We conclude that representing convective plumes as a sub‐grid scale process in models leads to a more realistic representation of open‐ocean polynyas and associated convection events.
Key Points
The Eddy‐Diffusivity Mass‐Flux (EDMF) parameterization is tested in a regional simulation of the ocean around Maud Rise
Thermobaric effects on convective plumes are enabled by modifying the EDMF parameterization
Simulations of Maud Rise polynyas are highly sensitive to freshwater forcing and mixing schemes
A one-dimensional Atmospheric Boundary Layer (ABL1D) model is coupled with the NEMO ocean model and implemented over the Iberian–Biscay–Ireland (IBI) area at 1/36° resolution to investigate the ...damping effect of the current and the thermal feedback on the kinetic energy (KE) at the mesoscale. This type of coupling between an ocean model and an ABL1D is a newly proposed approach as an alternative of intermediate complexity between bulk forcing and full coupling with an atmosphere model. In ABL1D, the prognostic tracers are nudged toward large-scale variables and the wind is guided by a low-frequency geostrophic wind provided from the ERA-Interim reanalyses. First, the ABL1D is successfully validated against satellite observations regarding the wind, and the dynamic coupling coefficient (linking the near surface wind and wind-stress to the of the surface currents) are consistent with the literature, over the period 2016–2017. Our results show that the thermal feedback has a negligible impact on kinetic energy (KE) and does not influence the strength of the current feedback in the region. Given the ABL1D physics, this further indicates that the changes in the vertical wind structure caused by CFB are primarily governed by local mechanical mechanisms associated with surface wind-stress condition, rather than by thermodynamic or non-local processes within the planetary boundary layer. The induced KE reduction by the current feedback amounts to 14% at the surface and propagates down to 2000 m, indicating that it can modify the vertical distribution of KE throughout the water column. KE reductions in the surface boundary layer (0 – 300 m) and in the interior (300 – 2000 m) are attributed to a reduction of the surface wind work by 4%, and of the pressure work by 7%, respectively. The Ekman pumping anomalies induced by the current feedback tend to attenuate eddy activity and horizontal pressure gradients at depth, illustrating the potential of the current feedback to induce a geostrophic adjustment on the water column. These results illustrate the relevance of the proposed ABL1D coupling approach for reproducing the wind-current coupling (a.k.a. current feedback effect) which cannot be taken into account straightforwardly with simple bulk forcing.
•A one-dimensional Atmospheric Boundary Layer (ABL1D) model is coupled with the NEMO ocean model over the Iberian–Biscay–Ireland area.•Our model is validated against satellite data and used to study current and thermal feedback effects on kinetic energy (KE).•The SST-wind coupling has a negligible impact on KE and on the current feedback (CFB) strength over the study area.•The wind structure is influenced by CFB mainly through local mechanical processes, not thermodynamic or non-local ones.•The CFB reduces surface KE by 14% and extends down to 2000m, altering KE distribution throughout the water column.•Our results highlights the relevance of ABL1D coupling to capture wind-current interaction compared to bulk forcing methods.