Monthly mean wind stress and Sverdrup transport in the North Atlantic: A comparison of the Hellerman-Rosenstein and Isemer-Hasse climatologies

The monthly mean wind stress climatology of Hellerman and Rosenstein (HR) is compared with the climatology of Isemer and Hasse (IH), which represents a version of the Bunker atlas (BU) for the North Atlantic based on revised parameterizations. The drag coefficients adopted by IH are 21% smaller than the values of BU and HR, and the calculation of wind speed from marine estimates of Beaufort force (Bft) is based on a revised Beaufort equivalent scale similar to the scientific scale recommended by WMO. The latter choice significantly increases wind speed below Bft 8, and effectively counteracts the reduction of the drag coefficients. Comparing the IH stresses with HR reveals substantially enhanced magnitudes in the trade wind region throughout the year. At 15°N the mean easterly stress increases from about 0.9 (HR) to about 1.2 dyn cm−1 (IH). Annual mean differences are smaller in the region of the westerlies. In winter, the effect due to the reduced drag coefficient dominates and leads to smaller stress values in IH; during summer season the revision of the Beaufort equivalents is more effective and leads to increased stresses. Implications of the different wind stress climatologies for forcing the large-scale ocean circulation are discussed by means of the Sverdrup transport streamfunction (ψs): Throughout the subtropical gyre a significant intensification of ψs takes place with IH. At 27°N, differences of more than 10 Sv (1 Sv ≡ 106 m3 s−1) are found near the western boundary. Differences in the seasonality of ψs are more pronounced in near-equatorial regions where IH increase the amplitude of the annual cycle by about 50%. An eddy-resolving model of the North Atlantic circulation is used to examine the effect of the different wind stresses on the seasonal cycle of the Florida Current. The transport predicted by the numerical model is in much better agreement with observations when the circulation is forced by IH than by HR, regarding both the annual mean (29.1 Sv vs 23.2 Sv) and the seasonal range (6.3 Sv vs 3.4 Sv)

Acute sensitivity of global ocean circulation and heat content to eddy energy dissipation time-scale

The global ocean overturning circulation, critically dependent on the global density stratification, plays a central role in regulating climate evolution. While it is well-known that the global stratification profile exhibits a strong dependence to Southern Ocean dynamics and in particular to wind and buoyancy forcing, we demonstrate here that the stratification is also acutely sensitive to the mesoscale eddy energy dissipation time-scale. Within the context of a global ocean circulation model with an energy constrained mesoscale eddy parameterization, it is shown that modest variations in the eddy energy dissipation time-scale lead to significant variations in key metrics relating to ocean circulation, namely the Antarctic Circumpolar Current transport, Atlantic Meridional Overturning Circulation strength, and global ocean heat content, over long time-scales. The results highlight a need to constrain uncertainties associated with eddy energy dissipation for climate model projections over centennial time-scales, but also for paleoclimate simulations over millennial time-scales.Comment: 13 pages, 4 figures; preprint version; accepted by Geophysical Research Letters, post-print version to be made available at a later poin

Stability of the thermohaline circulation examined with a one-dimensional fluid loop

The Stommel box model elegantly demonstrates that the oceanic response to mixed boundary conditions, combining a temperature relaxation with a fixed salt flux forcing, is nonlinear owing to the so-called salt advection feedback. This nonlinearity produces a parameter range of bi-stability associated with hysteresis effects characterised by a fast thermally-driven mode and a slow salinity-driven mode. Here we investigate whether a similar dynamical behaviour can be found in the thermohaline loop model, a one-dimensional analogue of the box model. A semi-analytical method to compute possible steady states of the loop model is presented, followed by a linear stability analysis carried out for a large range of loop configurations. While the salt advection feedback is found as in the box model, a major difference is obtained for the fast mode: an oscillatory instability is observed near the turning point of the fast mode branch, such that the range of bi-stability is systematically reduced, or even removed, in some cases. The oscillatory instability originates from a salinity anomaly that grows exponentially as it turns around the loop, a situation that may occur only when the salinity torque is directed against the loop flow. Factors such as mixing intensity, the relative strength of thermal and haline forcings, the nonlinearity of the equation of state or the loop geometry can strongly affect the stability properties of the loop

Stability and accuracy of Runge-Kutta based split-explicit time-stepping algorithms for free-surface ocean models

International audienceBecause of the Boussinesq assumption employed in the vast majority of oceanic models, the acoustic waves are filtered and the fast dynamics corresponds to the external gravity-wave propagation which is much faster than other (internal) processes. The fast and slow dynamics are traditionally split into separate subproblems where the fast motions are nearly independent of depth. It is thus natural to model these motions with a two-dimensional (barotropic) system of equations while the slow processes are modeled with a three-dimensional (baroclinic) system. However such splitting is inexact, the barotropic mode is not strictly depth-independent meaning that the separation of slow and fast modes is non-orthogonal, even in the linear case. A consequence is that there are some fast components contained in the slow motions which induce instabilities controlled by time filtering of the fast mode. In this talk we present an analysis of the stability and accuracy of the barotropic–baroclinic mode splitting in the case where the baroclinic mode is integrated using a Runge-Kutta scheme and the barotropic mode is integrated explicitly (i.e. the so-called split-explicit approach). By referring to the theoretical framework developed by Demange et al. (2019), the analysis is based on an eigenvector decomposition using the true (depth-dependent) barotropic mode. We investigate several strategies to achieve stable integrations whose performance is assessed first on a theoretical ground and then by idealized linear and nonlinear numerical experiments

The formation of the ocean’s anthropogenic carbon reservoir

The shallow overturning circulation of the oceans transports heat from the tropics to the mid-latitudes. This overturning also influences the uptake and storage of anthropogenic carbon (Cant). We demonstrate this by quantifying the relative importance of ocean thermodynamics, circulation and biogeochemistry in a global biochemistry and circulation model. Almost 2/3 of the Cant ocean uptake enters via gas exchange in waters that are lighter than the base of the ventilated thermocline. However, almost 2/3 of the excess Cant is stored below the thermocline. Our analysis shows that subtropical waters are a dominant component in the formation of subpolar waters and that these water masses essentially form a common Cant reservoir. This new method developed and presented here is intrinsically Lagrangian, as it by construction only considers the velocity or transport of waters across isopycnals. More generally, our approach provides an integral framework for linking ocean thermodynamics with biogeochemistry

Water masses as a unifying framework for understanding the Southern Ocean Carbon Cycle

International audienceThe scientific motivation for this study is to understand the processes in the ocean interior controlling carbon transfer across 30° S. To address this, we have developed a unified framework for understanding the interplay between physical drivers such as buoyancy fluxes and ocean mixing, and carbon-specific processes such as biology, gas exchange and carbon mixing. Given the importance of density in determining the ocean interior structure and circulation, the framework is one that is organized by density and water masses, and it makes combined use of Eulerian and Lagrangian diagnostics. This is achieved through application to a global ice-ocean circulation model and an ocean biogeochemistry model, with both components being part of the widely-used IPSL coupled ocean/atmosphere/carbon cycle model. Our main new result is the dominance of the overturning circulation (identified by water masses) in setting the vertical distribution of carbon transport from the Southern Ocean towards the global ocean. A net contrast emerges between the role of Subantarctic Mode Water (SAMW), associated with large northward transport and ingassing, and Antarctic Intermediate Water (AAIW), associated with a much smaller export and outgassing. The differences in their export rate reflects differences in their water mass formation processes. For SAMW, two-thirds of the surface waters are provided as a result of the densification of thermocline water (TW), and upon densification this water carries with it a substantial diapycnal flux of dissolved inorganic carbon (DIC). For AAIW, principal formatin processes include buoyancy forcing and mixing, with these serving to lighten CDW. An additional important formation pathway of AAIW is through the effect of interior processing (mixing, including cabelling) that serve to densify SAMW. A quantitative evaluation of the contribution of mixing, biology and gas exchange to the DIC evolution per water mass reveals that mixing and, secondarily, gas exchange, effectively nearly balance biology on annual scales (while the latter process can be dominant at seasonal scale). The distribution of DIC in the northward flowing water at 30° S is thus primarily set by the DIC values of the water masses that are involved in the formation processes

Sliding or stumbling on the staircase: numerics of ocean circulation along piecewise-constant coastlines

Coastlines in most ocean general circulation models are piecewise constant. Accurate representation of boundary currents along staircase-like coastlines is a long-standing issue in ocean modelling. Pioneering work by Adcroft and Marshall (1998) revealed that artificial indentation of model coastlines, obtained by rotating the numerical mesh within an idealized square basin, generates a \textit{spurious form drag} that slows down the circulation. Here, we revisit this problem and show how this spurious drag may be eliminated. First, we find that \textit{physical} convergence (i.e. the main characteristics of the flow are insensitive to the increase of the mesh resolution) allows simulations to become independent of the mesh orientation. An advection scheme with a wider stencil also reduces sensitivity to mesh orientation from coarser resolution. Second, we show that indented coastlines behave as straight and slippery shores when a true mirror boundary condition on the flow is imposed. This finding applies to both symmetric and rotational-divergence formulations of the stress tensor, and to both flux and vector-invariant forms of the equations. Finally, we demonstrate that the detachment of a vortex flowing past an outgoing corner of the coastline is faithfully simulated with exclusive implementation of impermeability conditions. These results provide guidance for a better numerical treatment of coastlines (and isobaths) in ocean general circulation models

Wave effects on coastal upwelling and water level

Traditional atmosphere, ocean and wave models are run independently of each other. This means that the energy and momentum fluxes do not fully account for the impact of the oceanic wave field at the air-sea interface. In this study, the Stokes drift impact on mass and tracer advection, the Stokes-Coriolis forcing, and the sea-state-dependent momentum and energy fluxes are introduced into an ocean circulation model and tested for a domain covering the Baltic Sea and the North Sea. Sensitivity experiments are designed to investigate the influence on the simulation of storms and Baltic Sea upwelling. Inclusion of wave effects improves the model performance compared with the stand-alone circulation model in terms of sea level height, temperature and circulation. The direct sea-state-dependent momentum and turbulent kinetic energy fluxes prove to be of higher importance than the Stokes drift related effects investigated in this study (i.e., Stokes-Coriolis forcing and Stokes drift advection on tracers and on mass). The latter affects the mass and tracer advection but largely balances the influence of the Stokes-Coriolis forcing. The upwelling frequency changes by >10% along the Swedish coast when wave effects are included. In general, the strong (weak) upwelling probability is reduced (increased) when adding the wave effects. From the results, we conclude that inclusion of wave effects can be important for regional, high-resolution ocean models even on short time scales, suggesting that they should be introduced in operational ocean circulation models. However, care should be taken when introducing the Stokes-Coriolis forcing as it should be balanced by the Stokes drift in mass and tracer advection