Search

Oceanic mesoscale motions including eddies, meanders, fronts, and filaments comprise a dominant fraction of oceanic kinetic energy and contribute to the redistribution of tracers in the ocean such as heat, salt, and nutrients. This reservoir of mesoscale energy is regulated by the conversion of potential energy and transfers of kinetic energy across spatial scales. Whether and under what circumstances mesoscale turbulence precipitates forward or inverse cascades, and the rates of these cascades, remain difficult to directly observe and quantify despite their impacts on physical and biological processes. Here we use global observations to investigate the seasonality of surface kinetic energy and upper-ocean potential energy. We apply spatial filters to along-track satellite measurements of sea surface height to diagnose surface eddy kinetic energy across 60–300-km scales. A geographic and scale-dependent seasonal cycle appears throughout much of the midlatitudes, with eddy kinetic energy at scales less than 60 km peaking 1–4 months before that at 60–300-km scales. Spatial patterns in this lag align with geographic regions where an Argo-derived estimate of the conversion of potential to kinetic energy is seasonally varying. In midlatitudes, the conversion rate peaks 0–2 months prior to kinetic energy at scales less than 60 km. The consistent geographic patterns between the seasonality of potential energy conversion and kinetic energy across spatial scale provide observational evidence for the inverse cascade and demonstrate that some component of it is seasonally modulated. Implications for mesoscale parameterizations and numerical modeling are discussed. Significance Statement This study investigates the seasonality of upper-ocean potential and kinetic energy in the context of an inverse cascade, consisting of energy transfers to and through the mesoscale. Observations show a scale-dependent cycle in kinetic energy that coincides with temporal variability in mixed layer potential energy and progresses seasonally from smaller to larger scales. This pattern appears dominant over large regions of the ocean. Results are relevant to ocean and climate models, where a large fraction of ocean energy is often parameterized. A customizable code repository and dataset are provided to enable comparisons of model-based resolved and unresolved kinetic energy to observational equivalents. Implications result for a range of processes including mixed layer stratification and vertical structure of ocean currents.

Pacific Summer Water eddies and intrusions transport heat and salt from boundary regions into the western Arctic basin. Here we examine concurrent effects of lateral stirring and vertical mixing using microstructure data collected within a Pacific Summer Water intrusion with a length scale of ∼20 km. This intrusion was characterized by complex thermohaline structure in which warm Pacific Summer Water interleaved in alternating layers of m thickness with cooler water, due to lateral stirring and intrusive processes. Along interfaces between warm/salty and cold/freshwater masses, the density ratio was favorable to double-diffusive processes. The rate of dissipation of turbulent kinetic energy (ε) was elevated along the interleaving surfaces, with values up to 3 × 10−8 W kg−1 compared to background ε of less than 10−9 W kg−1. Based on the distribution of ε as a function of density ratio Rρ, we conclude that double-diffusive convection is largely responsible for the elevated ε observed over the survey. The lateral processes that created the layered thermohaline structure resulted in vertical thermohaline gradients susceptible to double-diffusive convection, resulting in upward vertical heat fluxes. Bulk vertical heat fluxes above the intrusion are estimated in the range of 0.2–1 W m−2, with the localized flux above the uppermost warm layer elevated to 2–10 W m−2. Lateral fluxes are much larger, estimated between 1000 and 5000 W m−2, and set an overall decay rate for the intrusion of 1–5 years.

Arctic Ocean properties and processes are highly relevant to the regional and global coupled climate system, yet still scarcely observed, especially in winter. Team OCEAN conducted a full year of physical oceanography observations as part of the Multidisciplinary drifting Observatory for the Study of the Arctic Climate (MOSAiC), a drift with the Arctic sea ice from October 2019 to September 2020. An international team designed and implemented the program to characterize the Arctic Ocean system in unprecedented detail, from the seafloor to the air-sea ice-ocean interface, from sub-mesoscales to pan-Arctic. The oceanographic measurements were coordinated with the other teams to explore the ocean physics and linkages to the climate and ecosystem. This paper introduces the major components of the physical oceanography program and complements the other team overviews of the MOSAiC observational program. Team OCEAN’s sampling strategy was designed around hydrographic ship-, ice- and autonomous platform-based measurements to improve the understanding of regional circulation and mixing processes. Measurements were carried out both routinely, with a regular schedule, and in response to storms or opening leads. Here we present along-drift time series of hydrographic properties, allowing insights into the seasonal and regional evolution of the water column from winter in the Laptev Sea to early summer in Fram Strait: freshening of the surface, deepening of the mixed layer, increase in temperature and salinity of the Atlantic Water. We also highlight the presence of Canada Basin deep water intrusions and a surface meltwater layer in leads. MOSAiC most likely was the most comprehensive program ever conducted over the ice-covered Arctic Ocean. While data analysis and interpretation are ongoing, the acquired datasets will support a wide range of physical oceanography and multi-disciplinary research. They will provide a significant foundation for assessing and advancing modeling capabilities in the Arctic Ocean.

For ice concentrations less than 85%, internal ice stresses in the sea ice pack are small and sea ice is said to be in free drift. The sea ice drift is then the result of a balance between Coriolis acceleration and stresses from the ocean and atmosphere. We investigate sea ice drift using data from individual drifting buoys as well as Arctic-wide gridded fields of wind, sea ice, and ocean velocity. We perform probabilistic inverse modeling of the momentum balance of free-drifting sea ice, implemented to retrieve the Nansen number, scaled Rossby number, and stress turning angles. Since this problem involves a nonlinear, underconstrained system, we used a Monte Carlo guided search scheme—the Neighborhood Algorithm—to seek optimal parameter values for multiple observation points. We retrieve optimal drag coefficients of CA=1.2×10−3 and CO=2.4×10−3 from 10-day averaged Arctic-wide data from July 2014 that agree with the AIDJEX standard, with clear temporal and spatial variations. Inverting daily averaged buoy data give parameters that, while more accurately resolved, suggest that the forward model oversimplifies the physical system at these spatial and temporal scales. Our results show the importance of the correct representation of geostrophic currents. Both atmospheric and oceanic drag coefficients are found to decrease with shorter temporal averaging period, informing the selection of drag coefficient for short timescale climate models.

Author Posting. © American Geophysical Union, 2019. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research-Oceans 124(2), (2019):863-881, doi:10.1029/2018JC014604.

Understanding and predicting sea ice dynamics and ice-ocean feedback processes requires accurate descriptions of momentum fluxes across the ice-ocean interface. In this study, we present observatio...

As the Arctic Ocean transitions into a seasonally ice-covered ocean, the processes that govern the sea ice melt cycle will undergo a fundamental shift. The summer melt cycle passes through several stages, such as snowmelt, meltpond formation and drainage, and basal melting; normally with surface melt processes occurring well before basal melt. Monitoring of atmospheric, sea ice and oceanographic properties from autonomous buoys deployed across the Beaufort Sea, combined with targeted satellite imagery, reveal a fundamental restructuring of these melt processes within the seasonal ice zone – where basal melt occurs before surface melt. We find this seemingly unremarkable change accelerates the basal melting of sea ice by ten-fold, as well as transforming the timing and flux of freshwater and heat into the upper-ocean. These processes play a pivotal role in determining upper-ocean stratification, and when combined with a modification of the under-ice light field, it could impact marine ecosystem function.

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.

The eddy field across the Arctic Ocean's 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.

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.

In ice-covered regions it is challenging to determine constituent budgets – for heat and momentum, but also for biologically and climatically active gases like carbon dioxide and methane. The harsh environment and relative data scarcity make it difficult to characterize even the physical properties of the ocean surface. Here, we sought to evaluate if numerical model output helps us to better estimate the physical forcing that drives the air–sea gas exchange rate (k) in sea ice zones. We used the budget of radioactive 222Rn in the mixed layer to illustrate the effect that sea ice forcing has on gas budgets and air–sea gas exchange. Appropriate constraint of the 222Rn budget requires estimates of sea ice velocity, concentration, mixed-layer depth, and water velocities, as well as their evolution in time and space along the Lagrangian drift track of a mixed-layer water parcel. We used 36, 9 and 2 km horizontal resolution of regional Massachusetts Institute of Technology general circulation model (MITgcm) configuration with fine vertical spacing to evaluate the capability of the model to reproduce these parameters. We then compared the model results to existing field data including satellite, moorings and ice-tethered profilers. We found that mode sea ice coverage agrees with satellite-derived observation 88 to 98 % of the time when averaged over the Beaufort Gyre, and model sea ice speeds have 82 % correlation with observations. The model demonstrated the capacity to capture the broad trends in the mixed layer, although with a significant bias. Model water velocities showed only 29 % correlation with point-wise in situ data. This correlation remained low in all three model resolution simulations and we argued that is largely due to the quality of the input atmospheric forcing. Overall, we found that even the coarse-resolution model can make a modest contribution to gas exchange parameterization, by resolving the time variation of parameters that drive the 222Rn budget, including rate of mixed-layer change and sea ice forcings.

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.

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 m[superscript 2]/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 m[superscript 2]/s to more than 600 m[superscript 2]/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., National Science Foundation (U.S.) (Award 1603557)

The article of record as published may be found at http://dx.doi.org/10.1525/elementa.241/ The interplay between sea ice concentration, sea ice roughness, ocean stratification, and momentum transfer to the ice and ocean is subject to seasonal and decadal variations that are crucial to understanding the present and future air-ice-ocean system in the Arctic. In this study, continuous observations in the Canada Basin from March through December 2014 were used to investigate spatial differences and temporal changes in under-ice roughness and momentum transfer as the ice cover evolved seasonally. Observations of wind, ice, and ocean properties from four clusters of drifting instrument systems were complemented by direct drill-hole measurements and instrumented overhead flights by NASA operation IceBridge in March, as well as satellite remote sensing imagery about the instrument clusters. Spatially, directly estimated ice-ocean drag coefficients varied by a factor of three with rougher ice associated with smaller multi-year ice floe sizes embedded within the first-year-ice/multi-year-ice conglomerate. Temporal differences in the ice-ocean drag coefficient of 20�30% were observed prior to the mixed layer shoaling in summer and were associated with ice concentrations falling below 100%. The ice-ocean drag coefficient parameterization was found to be invalid in September with low ice concentrations and small ice floe sizes. Maximum momentum transfer to the ice occurred for moderate ice concentrations, and transfer to the ocean for the lowest ice concentrations and shallowest stratification. Wind work and ocean work on the ice were the dominant terms in the kinetic energy budget of the ice throughout the melt season, consistent with free drift conditions. Overall, ice topography, ice concentration, and the shallow summer mixed layer all influenced mixed layer currents and the transfer of momentum within the air-ice-ocean system. The observed changes in momentum transfer show that care must be taken to determine appropriate parameterizations of momentum transfer, and imply that the future Arctic system could become increasingly seasonal. Office of Naval Research The deployment of the ITP-Vs and the subsequent analysis effort was supported by the Office of Naval Research under grants N00014-12-10140 and N00014-12-10799.

Ocean velocity in the Arctic is estimated using a three-axis acoustic travel-time current meter onboard an Ice-Tethered 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.