The ocean energy cycle is calculated using a new available potential energy (APE) decomposition, which partitions adiabatic buoyancy fluxes from diapycnal mixing, applied to results from the Estimating the Circulation and Climate of the Ocean, Phase II (ECCO2), eddy-permitting ocean state estimate and observed surface buoyancy fluxes from the WHOI OAFlux project. Compared with the traditional Lorenz energy cycle, this framework provides a more accurate estimate of the background potential energy (PE) of the global oceans and the surface generation and interior fluxes of APE. Calculations of the global energy budget using 16 yr of ECCO2 output suggest that the adiabatic portion of the general circulation is maintained by a balance between the mean wind-driven upwelling that increases APE (10.27 TW) and time-fluctuating processes, including mesoscale eddies, which release APE (20.27 TW). The APE generated by surface buoyancy fluxes (0.46 TW) is comparable to the generation by the mean winds. The global rate of irreversible mixing (0.46 TW), which balances surface APE generation, is consistent with previous estimates of the diapycnal fluxes associated with maintaining deep stratification (see Munk and Wunsch) and a global diapycnal diffusivity of O(1 3 10 24 ) m 2 s 21 . However, the net contribution of diapycnal mixing to the total potential energy is negligible, which suggests that mixing, contrary to one current paradigm, does not place a global demand on kinetic energy dissipation. However, there are regions where mixing is significant, for example, between 3000 and 5000 m (in ECCO2), where mixing increases PE by 0.1 TW. The work provides a new framework for separating adiabatic-diabatic fluxes and for monitoring the global rate of diapycnal mixing rate using measurable surface properties such as SST and heat flux.

Iron-rich layers are known to form in the stellar subsurface through a combination of gravitational settling and radiative levitation. Their presence, nature, and detailed structure can affect the excitation process of various stellar pulsation modes and must therefore be modeled carefully in order to better interpret Kepler asteroseismic data. In this paper, we study the interplay between atomic diffusion and fingering convection in A-type stars, as well as its role in the establishment and evolution of iron accumulation layers. To do so, we use a combination of three-dimensional idealized numerical simulations of fingering convection (which neglect radiative transfer and complex opacity effects) and one-dimensional realistic stellar models. Using the three-dimensional simulations, we first validate the mixing prescription for fingering convection recently proposed by Brown et al. (within the scope of the aforementioned approximation) and identify what system parameters (total mass of iron, iron diffusivity, thermal diffusivity, etc.) play a role in the overall evolution of the layer. We then implement the Brown et al. prescription in the Toulouse-Geneva Evolution Code to study the evolution of the iron abundance profile beneath the stellar surface. We find, as first discussed by Théado et al., that when the concurrent settling of helium is ignored, this accumulation rapidly causes an inversion in the mean molecular weight profile, which then drives fingering convection. The latter mixes iron with the surrounding material very efficiently, and the resulting iron layer is very weak. However, taking helium settling into account partially stabilizes the iron profile against fingering convection, and a large iron overabundance can accumulate. The opacity also increases significantly as a result, and in some cases it ultimately triggers dynamical convection. The direct effects of radiative acceleration on the dynamics of fingering convection (especially in the nonlinear regime) remain to be added in the future to improve the quantitative predictions of the model.

Linear theory for steady stratified flow over topography sets the range for topographic wavenumbers over which freely propagating internal waves are generated, and the radiation and breaking of these waves contribute to energy dissipation away from the ocean bottom. However, previous numerical work demonstrated that dissipation rates can be enhanced by flow over large scale topographies with wavenumbers outside of the lee wave radiative range. We conduct idealized 3D numerical simulations of steady stratified flow over 1D topography in a rotating domain and quantify vertical distribution of kinetic energy dissipation. We vary two parameters: the first determines whether the topographic obstacle is within the lee wave radiative range and the second, proportional to the topographic height, measures the degree of flow non-linearity. For certain combinations of topographic width and height, breaking occurs in pulses every inertial period, such that kinetic energy dissipation develops inertial periodicity. In these simulations, kinetic energy dissipation rates are also enhanced in the interior of the domain. In the radiative regime the inertial motions arise due to resonant wave-wave interactions. In the small wavenumber non-radiative regime, instabilities downstream of the obstacle can facilitate the generation and propagation of non-linearly forced inertial motions, especially as topographic height increase. In our simulations, dissipation rates for tall and wide non-radiative topography are comparable to those of radiative topography, even away from the bottom, which is relevant to the ocean where the topographic spectrum is such that wider abyssal hills also tend to be taller.

Uptake of atmospheric carbon by the ocean, especially at high latitudes, plays an important role in offsetting anthropogenic emissions. At the surface of the Southern Ocean south of 30∘S, the ocean carbon uptake, which had been weakening in 1990s, strengthened in the 2000s. However, sparseness of in-situ measurements in the ocean interior make it difficult to compute changes in carbon storage below the surface. Here we develop a machine-learning model, which can estimate concentrations of dissolved inorganic carbon (DIC) in the Southern Ocean up to 4 km depth only using data available at the ocean surface. Our model is fast and computationally inexpensive. We apply it to calculate trends in DIC concentrations over the past three decades and find that DIC decreased in the 1990s and 2000s, but has increased, in particular in the upper ocean since the 2010s. However, the particular circulation dynamics that drove these changes may have differed across zonal sectors of the Southern Ocean. While the near-surface decrease in DIC concentrations would enhance atmospheric CO2 uptake continuing the previously-found trends, weakened connectivity between surface and deep layers and build-up of DIC in deep waters could reduce the ocean’s carbon storage potential.

One of the proposed mechanisms for energy loss in the ocean is through dissipation of internal waves, in particular above rough topography where internal lee waves are generated. Rates of dissipation and diapycnal mixing are often estimated using linear internal wave generation theory and a constant value for mixing efficiency. However, previous oceanographic measurements found that non-linear dynamics may be important close to topography. In order to investigate the role of non-linear interactions, we conduct idealized 3D direct numerical simulations (DNS) of steady flow over 1D topography and vary the topographic height, which correlates to the degree of flow non-linearity. We analyze the spatial distribution of energy transfer rates between internal waves and the non-geostrophic portion of the time-mean flow, and of dissipation and diapycnal mixing rates. In our simulations with taller, more non-linear topographies, energy transfer rates are similar to previously unexplained oceanographic observations near topography: internal waves gain energy from time-mean flow through horizontal straining and lose energy through vertical shearing. In the tall topography simulations, buoyancy fluxes also play a significant role, consistent with observations but contrary to linear wave theory, suggesting that quasigeostrophy-based approximations and linear theory may not hold in some regions above rough topography. Both dissipation and mixing rates increase with topographic height, but their vertical distributions differ between topographic regimes. As such, the vertical profile of mixing efficiency is different for linear and non-linear topographic regimes, which may need to be incorporated into parameterizations of small-scale processes in models and estimates of ocean energy loss.