ALBERT

Description: A parameterization for eddy buoyancy fluxes for use in coarse-grid models is developed and tested against eddy-resolving simulations. The development is based on the assumption that the eddies are adiabatic (except near the surface) and the observation that the flux of buoyancy is affected by barotropic, depth-independent eddies. Like the previous parameterizations of Gent and McWilliams (GM) and Visbeck et al. (VMHS), the horizontal flux of a tracer is proportional to the local large-scale horizontal gradient of the tracer through a transfer coefficient assumed to be given by the product of a typical eddy velocity scale and a typical mixing length. The proposed parameterization differs from GM and VMHS in the selection of the eddy velocity scale, which is based on the kinetic energy balance of baroclinic eddies. The three parameterizations are compared to eddy-resolving computations in a variety of forcing configurations and for several sets of parameters. The VMHS and the energy balance parameterizations perform best in the tests considered here.

Description: It is demonstrated that eddy fluxes of buoyancy at the eastern and western boundaries maintain alongshore buoyancy gradients along the coast. Eddy fluxes arise near the eastern and western boundaries because on both coasts buoyancy gradients normal to the boundary are strong. The eddy fluxes are accompanied by mean vertical flows that take place in narrow boundary layers next to the coast where the geostrophic constraint is broken. These ageostrophic cells have a velocity component normal to the coast that balances the geostrophic mean velocity. It is shown that the dynamics in these thin ageostrophic boundary layers can be replaced by effective boundary conditions for the interior flow, relating the eddy flux of buoyancy at the seaward edge of the boundary layers to the buoyancy gradient along the coast. These effective boundary conditions are applied to a model of the thermocline linearized around a mean stratification and a state of rest. The linear model parameterizes the eddy fluxes of buoyancy as isopycnal diffusion. The linear model produces horizontal gradients of buoyancy along the eastern coast on a vertical scale that depends on both the vertical diffusivity and the eddy diffusivity. The buoyancy field of the linear model agrees very well with the mean state of an eddy-resolving computation. Because the east–west difference in buoyancy is related to the zonally integrated meridional velocity, the linear model successfully predicts the meridional overturning circulation.

Description: The role of the relative geometry of mechanical forcing (wind stress) and buoyancy forcing (prescribed surface temperature) in the maintenance of the main thermocline is explored. In particular, the role of the wind stress curl in enhancing or suppressing the generation of baroclinic eddies is studied in simplified domains. The dependence of key quantities, such as the depth of the thermocline and the maximum heat transport, on the external parameters such as diapycnal mixing and dissipation rate is examined. Qualitatively different regimes are found depending on the relative phase of the wind stress and surface buoyancy distribution. The most efficient arrangement for eddy generation has Ekman pumping (suction) in conjunction with high (low) surface buoyancy. This corresponds to the situation found in the midlatitudes, where the surface Ekman flow carries heat toward the warmer region (i.e., upgradient of the surface temperature). In this case, strong eddy fluxes are generated in order to counteract the upgradient heat transport by the Ekman cell. The result is a thermocline whose depth is independent of the diapycnal diffusivity. However, the competition between these opposing heat fluxes leads to a weak net heat transport, proportional to the diffusivity responsible for the diabatic forcing. This arrangement of wind stress provides a large source of available potential energy on which eddies can grow, so the mechanical energy balance for the eddies is consistent with a substantial eddy heat flux. When the same surface temperature distribution is paired with the opposite wind stress curl, the mean flow produces a sink, rather than a source, of available potential energy and eddies are suppressed. With this arrangement, typical of low latitudes and the subpolar regions, the Ekman overturning cell carries heat downgradient of the surface temperature. Thus, the net heat transport is almost entirely due to the Ekman flow and is independent of the diapycnal diffusivity. At the same time the thermocline is a thin, diffusive boundary layer. Quantitative scalings for the thermocline depth and the poleward heat transport in these two limiting cases are contrasted and successfully compared with eddy-resolving computations.

Description: The effect of the pole-to-pole surface temperature difference on the deep stratification and the strength of the global meridional overturning circulation (MOC) is examined in an eddy-resolving ocean model configured in an idealized domain roughly representing the Atlantic sector. Mesoscale eddies lead to qualitative differences in the mean stratification and the MOC compared to laminar (i.e., eddy free) models. For example, the spreading of fluid across the model’s representation of the Antarctic Circumpolar Current (ACC) no longer relies on the existence of a sill in the ACC. In addition, the deep- and bottom-water masses—roughly representing North Atlantic Deep Water (NADW) and Antarctic Bottom Water (ABW), respectively—are eroded by the eddies so that their zonal and meridional extents are much smaller than in the laminar case. It is found that if the north pole temperature is sufficiently warm, the formation of northern deep water is suppressed and the middepth cell is small and weak while the deep cell is large and vigorous. In contrast, if the north pole temperature is in the range of the southern channel temperatures, the middepth cell is large and strong while the deep cell has a reduced amplitude. This result is consistent with the predictions of the laminar theory of the MOC. In contrast to the laminar theory, realistically strong deep stratification is formed even if the temperature at the northern sinking site is warmer than any temperature found in the channel. Indeed, middepth stratification is actually stronger in the latter case than the former case.

Description: The equilibrium of an idealized flow driven at the surface by wind stress and rapid relaxation to nonuniform buoyancy is analyzed in terms of entropy production, mechanical energy balance, and heat transport. The flow is rapidly rotating, and dissipation is provided by bottom drag. Diabatic forcing is transmitted from the surface by isotropic diffusion of buoyancy. The domain is periodic so that zonal averaging provides a useful decomposition of the flow into mean and eddy components. The statistical equilibrium is characterized by quantities such as the lateral buoyancy flux and the thermocline depth; here, scaling laws are proposed for these quantities in terms of the external parameters. The scaling theory predicts relations between heat transport, thermocline depth, bottom drag, and diapycnal diffusivity, which are confirmed by numerical simulations. The authors find that the depth of the thermocline is independent of the diapycnal mixing to leading order, but depends on the bottom drag. This dependence arises because the mean stratification is due to a balance between the large-scale wind-driven heat transport and the heat transport due to baroclinic eddies. The eddies equilibrate at an amplitude that depends to leading order on the bottom drag. The net poleward heat transport is a residual between the mean and eddy heat transports. The size of this residual is determined by the details of the diapycnal diffusivity. If the diffusivity is uniform (as in laboratory experiments) then the heat transport is linearly proportional to the diffusivity. If a mixed layer is incorporated by greatly increasing the diffusivity in a thin surface layer then the net heat transport is dominated by the model mixed layer.