Search

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.

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.

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. Mixing length was best estimated using the wavenumber of the beginning of the inertial subrange or a variable drag coefficient. Analyses of this and future ITP-V datasets will advance understanding of ice–ocean interactions and their parameterizations in numerical models.

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.

[1] Observations of the spatial distribution and persistence of thermohaline structure are presented, and show how advection and diffusion affect a passive tracer. More than two years of underwater glider observations in the central subtropical North Pacific showed thermohaline variability over horizontal scales from 5 to 1300 km. Thermohaline fluctuations along isopycnals (spice fluctuations) were elevated in layers throughout the water column with the largest fluctuations near the surface and subtropical frontal regions. Fluctuations were uncorrelated between the layers but stirred by the same velocity field. Spice variance had local extrema in the vertical because of differences in source water properties and the influence of neighboring water masses. Spice variance spanned about three orders of magnitude along deeper isopycnals with larger variance where different water masses met and where velocity and vorticity variance were elevated. Horizontal wave number spectra of spice had slopes of −2 everywhere in the upper 1000 m. Submesoscale spice fluctuations had slopes in physical space near the ratio of the Coriolis parameter to the buoyancy frequency (f/N), consistent with predictions of quasi-geostrophic theory. In the mixed layer, thermohaline structure had a significant annual cycle with smaller interannual differences. Thermohaline fluctuations left behind during restratification and isolated from the mixed layer decayed with time because of diffusion along isopycnals. Horizontal diffusivity estimates in the remnant mixed layer were 0.4 m2 s−1 at 15–28 km wavelengths and 0.9 m2 s−1 at 35–45 km wavelengths.

[1] We discuss the seasonal evolution of upper-ocean thermohaline structure at small horizontal scales. The upper 350 m of a 1000 km long section in the subtropical North Pacific was observed in winter, spring, and summer with 3–14 km horizontal resolution. Four vertical regions had distinct density and salinity structure: the mixed layer, remnant mixed layer, high-stratification layer, and permanent thermocline. The remnant mixed layer consists of water from the winter mixed layer left over after restratification. The remnant mixed layer was most similar to the mixed layer in winter and spring, and most similar to the high-stratification layer below in summer. The high-stratification layer had elevated stratification that varied seasonally. The permanent thermocline varied little seasonally and was horizontally and vertically uniform in comparison. In all seasons, density ratios showed that mixed-layer θ-S differences tended to compensate in density with the strongest tendency toward compensation in winter. Density ratios were temperature dominated in the remnant mixed layer consistent with salt-fingering. Salinity anomalies were largest at the surface and decayed with depth in all seasons. Spectra of isopycnal depth and θ-S anomalies along isopycnals are compared between the three seasons and four vertical layers. Isopycnal depth variance at 30–46 km wavelengths decreased from winter to spring to summer by a factor of 2–10 in stratified regions. By treating salinity anomalies as a tracer, the effective isopycnal diffusivity in the remnant mixed layer was estimated to be 1.4 m2 s−1 over 30–46 km wavelengths.