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

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

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

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)

Fair winds and following seas remotely: modifying perceptions of fieldwork as a requirement in marine science to aid in diversifying the discipline

Pursuing an academic career in marine science requires a range of skills that can be applied across different contexts, including experimental or computational proficiency, policy engagement, teaching, and seagoing fieldwork. The tendency to advertise careers in marine science with imagery of research expeditions results in the perception that it is a requirement for a career in marine science, an indicator of competitiveness in this discipline. Historically, those participating in remote fieldwork over extended periods of time were perceived as “adventurous explorers, with a strong bias towards western, able-bodied men” (Nash et al., 2019). Use of imagery reinforcing such notions for marine scientists fails to recognize that this perception can be discouraging to individuals from other backgrounds who may be excluded from the discipline by a range of real and perceived participatory barriers. Such exclusionary factors include: caring responsibilities, physical mobility, challenging social environments, isolating and physically uncomfortable working environments, mental health challenges, and access to opportunity (Giles et al., 2020). Such barriers disproportionately affect diverse, underrepresented, and marginalized groups, who may therefore struggle to identify with marine science as a potential discipline in which to pursue a successful career. Current work toward achieving net zero targets within ocean research emphasizes the use of autonomous vehicles as alternatives to ocean-going ships (Storey, 2023), and the proposed concept of digital twinning would incorporate similar remote technology coupled with simulations and shore-based decision-making. The concept of digital twinning refers to the use of responsive autonomous platforms that can both collect data and be operated in response to that data, which could provide a non-field-based approach to delivering marine science while also potentially expanding the opportunities available for individuals not able or interested in working in the field. In distinguishing digital twinning from current approaches such as data assimilating models, Kritzinger et al. (2018) note the importance of a two-way data flow between the physical environment and its virtual representation, called a “digital twin,” which, for example, may lead to changes in deployment strategy or data collection by researchers. Because these twins can be controlled and simulated anywhere with access to sufficient computing power, shore-based individuals can interact with a virtual version of the physical environment without being physically present at sea. The technology to support a fully realized digital twin of the ocean is still under development, but its use would require a broader range of skills and roles in the discipline, many of which are not accurately conveyed by the prevailing marketing of field-based disciplines (see Mol and Atchinson, 2019, regarding geosciences). In order to fully integrate this new approach into marine science, employment of individuals with experience and training across a wide range of disciplines from software engineering to traditional field sampling is essential while also presenting the potential for making marine science more inclusive. Individuals for whom working at sea is not possible and/or desirable would be able to make equally valid contributions to such research projects via digital routes, without facing the many barriers fieldwork may present. This study explores the expectations of marine scientists, from both early and more established career stages, around the importance of field experience as a precursor or requirement for a successful marine science career, and also examines the advantages and disadvantages of using digital twinning as a complement to traditional field-based marine science

Temporal Variability of Diapycnal Mixing in Shag Rocks Passage

Diapycnal mixing rates in the oceans have been shown to have a great deal of spatial variability, but the temporal variability has been little studied. Here we present results from a method developed to calculate diapycnal diffusivity from moored Acoustic Doppler Current Profiler (ADCP) velocity shear profiles. An 18-month time series of diffusivity is presented from data taken by a LongRanger ADCP moored at 2400 m depth, 600 m above the sea floor, in Shag Rocks Passage, a deep passage in the North Scotia Ridge (Southern Ocean). The Polar Front is constrained to pass through this passage, and the strong currents and complex topography are expected to result in enhanced mixing. The spatial distribution of diffusivity in Shag Rocks Passage deduced from lowered ADCP shear is consistent with published values for similar regions, with diffusivity possibly as large as 90 × 10-4 m2 s-1 near the sea floor, decreasing to the expected background level of ~ 0.1 × 10-4 m2 s-1 in areas away from topography. The moored ADCP profiles spanned a depth range of 2400 to 1800 m; thus the moored time series was obtained from a region of moderately enhanced diffusivity. The diffusivity time series has 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. Using the observed stratification and an axisymmetric seamount, of similar dimensions to the one west of the mooring, in a model of baroclinic topographically-trapped waves, produces periods of 3.8 and 2.6 days, in agreement with the signals observed. The diffusivity is anti-correlated with the rotary coefficient (indicating that stronger mixing occurs during times of upward energy propagation), which suggests that mixing occurs due to the breaking of internal waves generated at topography

A comparison of five surface mixed layer models with a year of observations in the North Atlantic

Five upper ocean mixed layer models driven by ERA-Interim surface forcing are compared with a year of hydrographic observations of the upper 1000 m, taken at the Porcupine Abyssal Plain observatory site using profiling gliders. All the models reproduce sea surface temperature (SST) fairly well, with annual mean warm biases of 0.11  ° C (PWP model), 0.24  ° C (GLS), 0.31  ° C (TKE), 0.91  ° C (KPP) and 0.36  ° C (OSMOSIS). The main exception is that the KPP model has summer SSTs which are higher than the observations by nearly 3 ° . Mixed layer salinity (MLS) is not reproduced well by the models and the biases are large enough to produce a non-trivial density bias in the Eastern North Atlantic Central Water which forms in this region in winter. All the models develop mixed layers which are too deep in winter, with average winter mixed layer depth (MLD) biases between 160 and 228 m. The high variability in winter MLD is reproduced more successfully by model estimates of the depth of active mixing and/or boundary layer depth than by model MLD based on water column properties. After the spring restratification event, biases in MLD are small and do not appear to be related to the preceding winter biases. There is a very clear relationship between MLD and local wind stress in all models and in the observations during spring and summer, with increased wind speeds leading to deepening mixed layers, but this relationship is not present during autumn and winter. We hypothesize that the deepening of the MLD in autumn is so strongly driven by the annual cycle in surface heat flux that the winds are less significant in the autumn. The surface heat flux drives a diurnal cycle in MLD and SST from March onwards, though this effect is much more significant in the models than in the observations. We are unable to identify one model as definitely better than the others. The only clear differences between the models are KPP’s inability to accurately reproduce summer SSTs, and the OSMOSIS model’s more accurate reproduction of MLS