Aspects of southern ocean transport and mixing

Understanding and quantifying the circulation of the oceans and the driving mechanisms thereof is an important step in developing models which can accurately predict future climate change. In particular, model studies have shown that the spatial variability of diapycnal diffusivity, which represents the rate at which deep water returns to shallower depths by means of turbulent diapycnal mixing, is a critical factor controlling the strength and structure of the circulation. Efforts are therefore ongoing to measure diffusivity as extensively as possible, but temporal variability in diffusivity has not been widely addressed. Results from three Southern Ocean studies are presented in this thesis. Firstly, a high resolution hydrographic survey carried out on the northern flank of the Kerguelen Plateau identifies a complex meandering current system carrying a total eastward volume transport of 174 ± 22 Sv, mostly associated with the blended Subtropical Front/Subantarctic Front. Significant water mass transformation across isopycnals is not required to balance the budgets in this region. Secondly, results are presented which cast doubt on the advisability of using density profiles acquired using Conductivity-Temperature-Depth instruments to estimate diapycnal diffusivity (an attractive proposition due to low cost and widespread data availability) in areas of weak stratification such as the Southern Ocean, because the noise characteristics of the data result in inaccurate diffusivity estimates. Finally, a method is developed for estimating diffusivity from profiles of velocity shear acquired by moored acoustic Doppler current profilers. An 18-month time series of diffusivity estimates is derived with a median of 3.3 × 10−4 m2 s−1 and a range of 0.5 × 10−4 m2 s−1 to 57 × 10−4 m2 s−1. There is no significant signal at annual or semiannual periods, but there is evidence of signals at periods of approximately fourteen days (likely due to the spring-neaps tidal cycle), and at periods of 3.8 and 2.6 days most likely due to topographically-trapped waves propagating around the local seamount. More widespread application of this method would allow for an assessment of natural climate variability in diapycnal diffusivity

The vertical structure of upper ocean variability at the Porcupine Abyssal Plain during 2012-2013

This study presents the characterization of variability in temperature, salinity and oxygen concentration, including the vertical structure of the variability, in the upper 1000m of the ocean over a full year in the northeast Atlantic. Continuously profiling ocean gliders with vertical resolution between 0.5-1m provide more information on temporal variability throughout the water column than time series from moorings with sensors at a limited number of fixed depths. The heat, salt and dissolved oxygen content are quantified at each depth. While the near surface heat content is consistent with the net surface heat flux, heat content of the deeper layers is driven by gyre-scale water mass changes. Below ~150m, heat and salt content display intraseasonal variability which has not been resolved by previous studies. A mode-1 baroclinic internal tide is detected as a peak in the power spectra of water mass properties. The depth of minimum variability is at ~415m for both temperature and salinity, but this is a depth of high variability for oxygen concentration. The deep variability is dominated by the intermittent appearance of Mediterranean Water, which shows evidence of filamentation. Susceptibility to salt fingering occurs throughout much of the water column for much of the year. Between about 700-900m, the water column is susceptible to diffusive layering, particularly when Mediterranean Water is present. This unique ability to resolve both high vertical and temporal resolution highlights the importance of intraseasonal variability in upper ocean heat and salt content, variations that may be aliased by traditional observing techniques

Internal tide energy flux over a ridge measured by a co-located ocean glider and moored Acoustic Doppler Current Profiler

Internal tide energy flux is an important diagnostic for the study of energy pathways in the ocean, from large-scale input by the surface tide, to small-scale dissipation by turbulent mixing. Accurate calculation of energy flux requires repeated full-depth measurements of both potential density (ρ) and horizontal current velocity (u) over at least a tidal cycle and over several weeks to resolve the internal spring-neap cycle. Typically, these observations are made using full-depth oceanographic moorings that are vulnerable to being ‘fished-out’ by commercial trawlers when deployed on continental shelves and slopes. Here we test an alternative approach to minimise these risks, with u measured by a low-frequency ADCP moored near the seabed and ρ measured by an autonomous ocean glider holding station by the ADCP. The method is used to measure the semidiurnal internal tide radiating from the Wyville Thompson Ridge in the North Atlantic. The observed energy flux (4.2±0.2 kW m-1) compares favourably with historic observations and a previous numerical model study. Error in the energy flux calculation due to imperfect co-location of the glider and ADCP is estimated by sub-sampling potential density in an idealised internal tide field along pseudorandomly distributed glider paths. The error is considered acceptable (<10%) if all the glider data is contained within a ‘watch circle’ with a diameter smaller than 1/8 the mode-1 horizontal wavelength of the internal tide. Energy flux is biased low because the glider samples density with a broad range of phase shifts, resulting in underestimation of vertical isopycnal displacement and available potential energy. The negative bias increases with increasing watch circle diameter. If watch circle diameter is larger than 1/8 the mode-1 horizontal wavelength, the negative bias is more than 3% and all realisations within the 95% confidence interval are underestimates. Over the Wyville Thompson Ridge, where the semidiurnal mode-1 horizontal wavelength is ≈100 km and all the glider dives are within a 5 km diameter watch circle, the observed energy flux is estimated to have a negative bias of only 0.4% and an error of less than 3% at the 95% confidence limit. With typical glider performance, we expect energy flux error due to imperfect co-location to be <10% in most mid-latitude shelf slope regions

Estimating oceanic primary production using vertical irradiance and chlorophyll profiles from ocean gliders in the North Atlantic

An autonomous underwater vehicle (Seaglider) has been used to estimate marine primary production (PP) using a combination of irradiance and fluorescence vertical profiles. This method provides estimates for depth-resolved and temporally evolving PP on fine spatial scales in the absence of ship-based calibrations. We describe techniques to correct for known issues associated with long autonomous deployments such as sensor calibration drift and fluorescence quenching. Comparisons were made between the Seaglider, stable isotope (13C), and satellite estimates of PP. The Seaglider-based PP estimates were comparable to both satellite estimates and stable isotope measurements

Open-ocean submesoscale motions: a full seasonal cycle of mixed layer instabilities from gliders

The importance of submesoscale instabilities, particularly mixed-layer baroclinic instability and symmetric instability, on upper ocean mixing and energetics is well documented in regions of strong, persistent fronts such as the Kuroshio and the Gulf Stream. Less attention has been devoted to studying submesoscale flows in the open ocean, far from long-term mean geostrophic fronts, characteristic of a large proportion of the global ocean. We present a year-long, submesoscale-resolving time series of near-surface buoyancy gradients, potential vorticity and instability characteristics, collected by ocean gliders, that provides insight into open-ocean submesoscale dynamics over a full annual cycle. The gliders continuously sampled a 225 km2 region in the subtropical northeast Atlantic, measuring temperature, salinity and pressure along 292 short (\~{}20 km) hydrographic sections.Glider observations show a seasonal cycle in near-surface stratification. Throughout the fall (September through November), the mixed layer deepens, predominantly through gravitational instability, indicating that surface cooling dominates submesoscale restratification processes. During winter (December through March), mixed layer depths are more variable, and estimates of the balanced Richardson number, which measures the relative importance of lateral and vertical buoyancy gradients, depict conditions favorable to symmetric instability. The importance of mixed layer instabilities on the restratification of the mixed layer, as compared with surface heating and cooling, shows that submesoscale processes can reverse the sign of an equivalent heat flux up to 25{\%} of the time during winter. These results demonstrate that the open-ocean mixed layer hosts various forced and unforced instabilities, which become more prevalent during winter, and emphasize that accurate parameterizations of submesoscale processes are needed throughout the ocean

Seasonal extrema of sea surface temperature in CMIP6 models

CMIP6 model sea surface temperature (SST) seasonal extrema averaged over 1981–2010 are assessed against the World Ocean Atlas (WOA18) observational climatology. We propose a mask to identify and exclude regions of large differences between three commonly used climatologies (WOA18, WOCE-Argo Global Hydrographic climatology (WAGHC) and the Hadley Centre Sea Ice and Sea Surface Temperature data set (HadISST)). The biases in SST seasonal extrema are largely consistent with the annual mean SST biases. However, the amplitude and spatial pattern of SST bias vary seasonally in the 20 CMIP6 models assessed. Large seasonal variations in the SST bias occur in eastern boundary upwelling regions, polar regions, the North Pacific and the eastern equatorial Atlantic. These results demonstrate the importance of evaluating model performance not simply against annual mean properties. Models with greater vertical resolution in their ocean component typically demonstrate better representation of SST extrema, particularly seasonal maximum SST. No significant relationship of SST seasonal extrema with horizontal ocean model resolution is found

Wind-induced variability of warm water on the Southern Bellingshausen Sea continental shelf

The Bellingshausen Sea hosts heat transport onto the continental shelf, potentially enhancing ice shelf basal melt. Here we use the GLORYS12V1 1993-2018 reanalysis to identify physical processes that set seasonal and interannual variability of water mass properties in the Eltanin and Latady bays on the southern Bellingshausen Sea continental shelf. Annual means of potential temperature from 300 m to the seabed reveal interannual variability and allow separation into warm and cold regimes. The Amundsen Sea Low is more intense and extends further east during the warm regime than the cold regime. In the warm regime, a wind-induced reduction of sea ice concentration near the coast increases surface heat loss, convection, and formation of cold dense water in winter, associated with a decrease in heat content of the southern Bellingshausen Sea over time and a net northward heat transport. In contrast, in the cold regime, increased sea ice concentration reduces surface heat loss and thus formation of cold, dense water. Combined with an increase in heat content over time and a net southward heat transport, this results in a warming of the southern Bellingshausen Sea. This suggests that variability in the deep water temperature in the southern Bellingshausen Sea is primarily due to local surface heat fluxes above the shelf. The variability of surface heat fluxes is related to the variability of the Amundsen Sea Low and its influence on sea ice extent and local formation of cold, dense water in winter

Turbulent kinetic energy dissipation rate and associated fluxes in the western tropical Atlantic estimated from ocean glider observations

Ocean gliders enable us to collect the high-resolution microstructure observations necessary to calculate the dissipation rate of turbulent kinetic energy, ε, on timescales of weeks to months: far longer than is normally possible using traditional ship-based platforms. Slocum gliders have previously been used to this end; here, we report the first detailed estimates of ε calculated using the Batchelor spectrum method on observations collected by a FP07 fast thermistor mounted on a Seaglider. We use these same fast thermistor observations to calculate ε following the Thorpe scale method and find very good agreement between the two methods. The Thorpe scale method yields larger values of ε, but the average difference, which is less than an order of magnitude, is smaller than reported elsewhere. The spatio-temporal distribution of ε is comparable for both methods. Maximum values of ε (10−7 W kg−1) are observed in the surface mixed layer; values of approximately 10−9 W kg−1 are observed between approximately 200 and 500 m depth. These two layers are separated by a 100 m thick layer of low ε (10−10 W kg−1), which is co-located with a high-salinity layer of Subtropical Underwater and a peak in the strength of stratification. We calculate the turbulent heat and salt fluxes associated with the observed turbulence. Between 200 and 500 m, ε induces downward fluxes of both properties that, if typical of the annual average, would have a very small influence on the heat and salt content of the overlying salinity-maximum layer. We compare these turbulent fluxes with two estimates of double-diffusive fluxes that occur in regions susceptible to salt fingers, such as the western tropical Atlantic. We find that the double-diffusive fluxes of both heat and salt are larger than the corresponding turbulent fluxes

Stirring across the Antarctic Circumpolar Current's southern boundary at the prime meridian, Weddell Sea

At the southern boundary of the Antarctic Circumpolar Current (ACC), relatively warm ACC waters encounter the colder waters surrounding Antarctica. Strong density gradients across the southern boundary indicate the presence of a frontal jet and are thought to modulate the southward heat transport across the front. In this study, the southern boundary in the Weddell Sea sector at the prime meridian is surveyed for the first time in high resolution over 2 months during an austral summer with underwater gliders occupying a transect across the front on five occasions. The five transects show that the frontal structure (i.e. hydrography, velocities and lateral density gradients) varies temporally. The results demonstrate significant, transient (a few weeks) variability of the southern boundary and its frontal jet in location, strength and width. A mesoscale cold-core eddy is identified to disrupt the southern boundary’s frontal structure and strengthen lateral density gradients across the front. The front's barrier properties are assessed using mixing length scales and potential vorticity to establish the cross-frontal exchange of properties between the ACC and the Weddell Gyre. The results show that stronger lateral density gradients caused by the mesoscale eddy strengthen the barrier-like properties of the front through reduced mixing length scales and pronounced gradients of potential vorticity. In contrast, the barrier-like properties of the southern boundary are reduced when no mesoscale eddy is influencing the density gradients across the front. Using satellite altimetry, we further demonstrate that the barrier properties over the past decade have strengthened as a result of increased meridional gradients of absolute dynamic topography and increased frontal jet speeds in comparison to previous decades. Our results emphasise that locally and rapidly changing barrier properties of the southern boundary are important to quantify the cross-frontal exchange, which is particularly relevant in regions where the southern boundary is located near the Antarctic shelf break (e.g. in the West Antarctic sector)

Seasonality of submesoscale flows in the ocean surface boundary layer

A signature of submesoscale flows in the upper ocean is skewness in the distribution of relative vorticity. Expected to result for high Rossby-number flows, such skewness has implications for mixing, dissipation and stratification within the upper ocean. An array of moorings deployed in the Northeast Atlantic for one year as part of the OSMOSIS experiment reveals that relative vorticity is positively skewed during winter even though the scale of the Rossby number is less than 0.5. Furthermore, this skewness is reduced to zero during spring and autumn. There is also evidence of modest seasonal variations in the gradient Rossby number. The proposed mechanism by which relative vorticity is skewed is that the ratio of lateral to vertical buoyancy gradients, as summarized by the inverse gradient Richardson number, restricts its range during winter but less so at other times of the year. These results support recent observations and model simulations suggesting the upper ocean is host to a seasonal cycle in submesoscale turbulence