David Marshall

Abstract:

Ocean ventilation translates atmospheric forcing into the ocean interior. The Southern Ocean is an important ventilation site for heat and carbon and is likely to influence the outcome of anthropogenic climate change. We conduct an extensive backwards-in-time trajectory experiment to identify spatial and temporal patterns of ventilation. Temporally, almost all ventilation occurs between August and November. Spatially, “hotspots” of ventilation account for 60% of open-ocean ventilation on a 30 years timescale; the remaining 40% ventilates in a circumpolar pattern. The densest waters ventilate on the Antarctic shelf, primarily near the Antarctic Peninsula (40%) and the west Ross sea (20%); the remaining 40% is distributed across East Antarctica. Shelf-ventilated waters experience significant densification outside of the mixed layer.

Abstract:

There is an increasing interest in mesoscale eddy parameterizations that are scale-aware, normally interpreted to mean that a parameterization does not require parameter recalibration as the model resolution changes. Here we examine whether Gent–McWilliams (GM) based version of GEOMETRIC, a mesoscale eddy parameterization that is constrained by a parameterized eddy energy budget, is scale-aware in its energetics. It is generally known that GM-based schemes severely damp out explicit eddies, so the parameterized component would be expected to dominate across resolutions, and we might expect a negative answer to the question of energetic scale-awareness. A consideration of why GM-based schemes damp out explicit eddies leads a suggestion for what we term a splitting procedure: a definition of a “large-scale” field is sought, and the eddy-induced velocity from the GM-scheme is computed from and acts only on the large-scale field, allowing explicit and parameterized components to co-exist. Within the context of an idealized re-entrant channel model of the Southern Ocean, evidence is provided that the GM-based version of GEOMETRIC is scale-aware in the energetics as long as we employ a splitting procedure. The splitting procedure also leads to an improved representation of mean states without detrimental effects on the explicit eddy motions.

Abstract:

Estimates of the Weddell Gyre transport vary widely between climate simulations. Here, we investigate if inter-model variability can originate from differences in the horizontal resolution of the ocean model. We run an idealized model of the Weddell Gyre at eddy-parameterized, eddy-permitting, and eddy-rich resolutions and find that the gyre is strongly sensitive to horizontal resolution. The gyre transport is largest at eddy-permitting resolutions (45 Sv with a noisy bathymetry) and smallest at eddy-parameterized resolutions (12 Sv). The eddy-permitting simulations have the largest horizontal density gradients and the weakest stratification over the gyre basin. The large horizontal density gradients induce a significant thermal wind transport and increase the mean available potential energy for mesoscale eddies. The distribution of eddy kinetic energy indicates that explicit eddies in simulations intensify the bottom circulation of the gyre via non-linear dynamics. If climate models adopt horizontal resolutions that the Weddell Gyre is most sensitive to, then simulations of the Weddell Gyre could become more disparate.

Abstract:

Although the surface eddy kinetic energy (EKE) has been well studied using satellite altimeter and surface drifter observations, our knowledge of EKE in the ocean interior is much more limited due to the sparsity of subsurface current measurements. Here we develop a new approach for estimating EKE over the full depth of the global ocean by combining 20 years of satellite altimeter and Argo float data to infer the vertical profile of eddies. The inferred eddy profiles are surface-intensified at low latitudes and deep-reaching at mid- and high latitudes. They compare favorably to the first empirical orthogonal function obtained from current meter velocities. The global-integrated EKE estimated from the inferred profiles is about 3.1 × 10¹⁸ J, which is close to that estimated from the surface mode (3.0 × 10¹⁸ J) but about 30% smaller than that estimated from the traditional flat bottom modes (4.6 × 10¹⁸ J).