Sylvia T. Cole

Characterizing the eddy field in the Arctic Ocean halocline

Ice-Tethered Profilers (ITP), deployed in the Arctic Ocean between 2004 and 2013, have provided detailed temperature and salinity measurements of an assortment of halocline eddies. A total of 127 mesoscale eddies have been detected, 95% of which were anticyclones, the majority of which had anomalously cold cores. These cold-core anticyclonic eddies were observed in the Beaufort Gyre region (Canadian water eddies) and the vicinity of the Transpolar Drift Stream (Eurasian water eddies). An Arctic-wide calculation of the first baroclinic Rossby deformation radius Rd has been made using ITP data coupled with climatology; Rd ∼ 13 km in the Canadian water and ∼8 km in the Eurasian water. The observed eddies are found to have scales comparable to Rd. Halocline eddies are in cyclogeostrophic balance and can be described by a Rankine vortex with maximum azimuthal speeds between 0.05 and 0.4 m/s. The relationship between radius and thickness for the eddies is consistent with adjustment to the ambient stratification. Eddies may be divided into four groups, each characterized by distinct core depths and core temperature and salinity properties, suggesting multiple source regions and enabling speculation of varying formation mechanisms.

Horizontal density structure and restratification of the Arctic Ocean surface layer

Ice-tethered profiler (ITP) measurements from the Arctic Ocean’s Canada Basin indicate an ocean surface layer beneath sea ice with significant horizontal density structure on scales of hundreds of kilometers to the order 1 km submesoscale. The observed horizontal gradients in density are dynamically important in that they are associated with restratification of the surface ocean when dense water flows under light water. Such restratification is prevalent in wintertime and competes with convective mixing upon buoyancy forcing (e.g., ice growth and brine rejection) and shear-driven mixing when the ice moves relative to the ocean. Frontal structure and estimates of the balanced Richardson number point to the likelihood of dynamical restratification by isopycnal tilt and submesoscale baroclinic instability. Based on the evidence here, it is likely that submesoscale processes play an important role in setting surface-layer properties and lateral density variability in the Arctic Ocean.

Ekman Veering, Internal Waves, and Turbulence Observed under Arctic Sea Ice

The ice‐ocean system is investigated on inertial to monthly time scales using winter 2009‐10 observations from the first ice-tethered profiler (ITP) equipped with a velocity sensor (ITP-V). Fluctuations in surface winds, ice velocity, and ocean velocity at 7-m depth were correlated. Observed ocean velocity was primarily directed to the right of the ice velocity and spiraled clockwise while decaying with depth through the mixed layer. Inertial and tidal motions of the ice and in the underlying ocean were observed throughout the record. Just below the ice‐ocean interface, direct estimates of the turbulent vertical heat, salt, and momentum fluxes and the turbulent dissipation rate were obtained. Periods of elevated internal wave activity were associated with changes to the turbulent heat and salt fluxes as well as stratification primarily within the mixed layer. Turbulent heat and salt fluxes were correlated particularly when the mixed layer was closest to the freezing temperature. Momentum flux is adequately related to velocity shear using a constant ice‐ocean drag coefficient, mixing length based on the planetary and geometric scales, or Rossby similarity theory. Ekman viscosity described velocity shear over the mixed layer. The ice‐ocean drag coefficient was elevated for certain directions of the ice‐ocean shear, implying an ice topography that was characterized by linear ridges. Mixinglengthwasbestestimatedusingthewavenumber ofthebeginningoftheinertialsubrangeoravariable drag coefficient. Analyses of this and future ITP-V datasets will advance understanding of ice‐ocean interactions and their parameterizations in numerical models.

Observations of Tidal Internal Wave Beams at Kauai Channel, Hawaii

Abstract To observe the across-ridge structure of internal tides, density and velocity were measured using SeaSoar and a Doppler sonar over the upper 400–600 m of the ocean extending 152 km on each side of the Hawaiian Ridge at Kauai Channel. Eighteen sections were completed in about 18 days with sampling intentionally detuned from the lunar semidiurnal (M2) tide so that averaging over all sections was equivalent to phase averaging the M2 tide. Velocity and displacement variance and several covariances involving velocity and displacement showed one M2 internal wave beam on each side of the ridge and reflection of the beams off of the surface. Theoretical ray slopes aligned with the observed beams and originated from the sides of the ridge. Energy flux was in agreement with internal wave generation at the ridge. Inferred turbulent dissipation was elevated relative to open ocean values near tidal beams. Energy flux was larger than total dissipation almost everywhere across the ridge. Internal wave energy fl...

Eddy stirring and horizontal diffusivity from Argo float observations: Geographic and depth variability

Stirring along isopycnals is a significant factor in determining the distribution of tracers within the ocean. Salinity anomalies on density surfaces from Argo float profiles are used to investigate horizontal stirring and estimate eddy mixing lengths. Eddy mixing length and velocity fluctuations from the ECCO2 global state estimate are used to estimate horizontal diffusivity at a 300 km scale in the upper 2000 m with near-global coverage. Diffusivity varies by over two orders of magnitude with latitude, longitude, and depth. In all basins, diffusivity is elevated in zonal bands corresponding to strong current regions, including western boundary current extension regions, the Antarctic Circumpolar Current, and equatorial current systems. The estimated mixing lengths and diffusivities provide an observationally based data set that can be used to test and constrain predictions and parameterizations of eddy stirring.

The Beaufort Gyre intensification and stabilization: A model-observation synthesis

A model-observation synthesis is conducted to investigate changes in the upper ocean circulation and stratification in the Canada Basin [CB] of the Arctic Ocean. Results show that the Beaufort Gyre [BG] has been generally intensifying during 1992–2015 in conjunction with changes in sea ice and the upper ocean including increasing sea surface height [SSH], sea ice and ocean speed, Ekman transport convergence and downwelling, and freshwater content, decreasing ice thickness and upper ocean salinity, shoaling summer halocline and mixed layer, and deepening winter halocline and mixed layer. Increasing Ekman transport convergence draws more water from surrounding areas into the CB, thus lowering SSH in those areas and raising SSH in the CB. The rate of change in the CB began to decrease in 2008 and the BG circulation appears to be stabilizing, if not relaxing slightly. This is reflected in the general plateauing of SSH, the intensity of the sea ice and ocean circulation, and various measures of the CB thermohaline stratification. The BG intensification and subsequent stabilization appear to have been strongly controlled by atmospheric changes in the CB characterized by generally increasing anticyclonic wind circulation and sea level pressure (SLP) before 2008 and falling wind strength and SLP to below-average levels in some years after 2008. Changes in SLP are highly correlated with changes in ocean surface stress curl and downwelling. Since 2008, the magnitude of the stress curl and downwelling in much of the CB has declined, contributing to BG stabilization. The general leveling-off of sea ice thickness also contributes to the stabilization by limiting melt water input to the CB that increases freshwater content. Temperatures in the Near Surface Temperature Maximum layer trended upward slightly over 1992–2015, which is closely correlated with decreasing sea ice thickness. Upper ocean heat content increased over the study period mainly due to strong temperature increases in the summer Pacific Water layer.

Evolution of a Canada Basin ice-ocean boundary layer and mixed layer across a developing thermodynamically forced marginal ice zone

A comprehensive set of autonomous, ice-ocean measurements were collected across the Canada Basin to study the summer evolution of the ice-ocean boundary layer (IOBL) and ocean mixed layer (OML). Evaluation of local heat and freshwater balances and associated turbulent forcing reveals that melt ponds (MPs) strongly influence the summer IOBL-OML evolution. Areal expansion of MPs in mid-June start the upper ocean evolution resulting in significant increases to ocean absorbed radiative flux (19 W m−2 in this study). Buoyancy provided by MP drainage shoals and freshens the IOBL resulting in a 39 MJ m−2 increase in heat storage in just 19 days (52% of the summer total). Following MP drainage, a near-surface fresh layer deepens through shear-forced mixing to form the summer mixed layer (sML). In late summer, basal melt increases due to stronger turbulent mixing in the thin sML and the expansion of open water areas due in part to wind-forced divergence of the sea ice. Thermal heterogeneities in the marginal ice zone (MIZ) upper ocean led to large ocean-to-ice heat fluxes (100–200 W m−2) and enhanced basal ice melt (3–6 cm d−1), well away from the ice edge. Calculation of the upper ocean heat budget shows that local radiative heat input accounted for at least 89% of the observed latent heat losses and heat storage (partitioned 0.77/0.23). These results suggest that the extensive area of deteriorating sea ice observed away from the ice edge during the 2014 season, termed the “thermodynamically forced MIZ,” was driven primarily by local shortwave radiative forcing.

Evolution of the eddy field in the Arctic Ocean's Canada Basin, 2005–2015

The eddy field across the Arctic Oceans Canada Basin is analyzed using Ice-Tethered Profiler (ITP) and moored measurements of temperature, salinity, and velocity spanning 2005 to 2015. ITPs encountered 243 eddies, 98% of which were anticyclones, with approximately 70% of these having anomalously cold cores. The spatially and temporally varying eddy field is analyzed accounting for sampling biases in the unevenly distributed ITP data and caveats in detection methods. The highest concentration of eddies was found in the western and southern portions of the basin, close to topographic margins and boundaries of the Beaufort Gyre. The number of lower halocline eddies approximately doubled from 2005–2012 to 2013–2014. The increased eddy density suggests more active baroclinic instability of the Beaufort Gyre that releases available potential energy to balance the wind energy input; this may stabilize the Gyre spin-up and associated freshwater increase.

Observational Inferences of Lateral Eddy Diffusivity in the Halocline of the Beaufort Gyre

Using Ekman pumping rates mediated by sea ice in the Arctic Ocean’s Beaufort Gyre (BG), the magnitude of lateral eddy diffusivities required to balance downward pumping is inferred. In this limit—that of vanishing residual-mean circulation—eddy-induced upwelling exactly balances downward pumping. The implied eddy diffusivity varies spatially and decays with depth, with values of 50–400 m2/s. Eddy diffusivity estimated using mixing length theory applied to BG mooring data exhibits a similar decay with depth and range of values from 100 m2/s to more than 600 m2/s. We conclude that eddy diffusivities in the BG are likely large enough to balance downward Ekman pumping, arresting the deepening of the gyre and suggesting that eddies play a zero-order role in buoyancy and freshwater budgets of the BG.

Processing of velocity observations from Ice-Tethered Profilers

Ocean velocity in the Arctic is estimated using a three-axis acoustic travel-time current meter onboard an IceTethered Profiler (ITP). The ITP is moored on an ice floe, with the current meter measuring velocity past the sensor as it profiles vertically through the water column along an inclined wire and translates horizontally with the motion of the overlying sea ice. Accurate estimates of ocean velocity rely on simultaneous estimates of instrument pitch, roll, and heading, the velocity at which the sensor crawls along the wire, and the velocity of the overlying sea ice. An inertial measurement unit is used to derive instrument pitch, roll, and heading from a three-axis flux-gate compass, accelerometer, and angle rate gyro. Instrument speed along the wire is derived from pressure as measured by the CTD. Hourly GPS fixes of the surface package are used to estimate ice velocity. Pre-deployment calibrations are preformed to determine magnetometer, accelerometer, pitch, roll, and heading biases. Post-deployment calibration parameters are determined to account for remaining biases. Post-deployment parameters are chosen so that there is no obvious bias related to wire angle, velocity estimates when the profiler is stationary versus moving, or upward versus downward profiles. These procedures are described as implemented for the second-generation velocity sensor deployed in the Arctic Ocean on five Ice-Tethered Profilers during 2013-2015.