Hydrographic data collected from research cruises, bottom-anchored moorings, driftingIce-Tethered Profilers, and satellite altimetry in the Beaufort Gyre region of the Arctic Ocean document an increase of more than 6,400 km 3 of liquid freshwater content from 2003 to 2018: a 40% growth relative to the climatology of the 1970s. This fresh water accumulation is shown to result from persistent anticyclonic atmospheric wind forcing accompanied by sea ice melt, a wind-forced redirection of Mackenzie River discharge from predominantly eastward to westward flow, and a contribution of low salinity waters of Pacific Ocean origin via Bering Strait. Despite significant uncertainties in the different observations, this study has demonstrated the synergistic value of having multiple diverse datasets to obtain a more comprehensive understanding of Beaufort Gyre freshwater content variability. For example, Beaufort Gyre Observational System (BGOS) surveys clearly show the interannual increase in freshwater content, but without satellite or Ice-Tethered Profiler measurements, it is not possible to resolve the seasonal cycle of freshwater content, which in fact is larger than the year-to-year variability, or the more subtle interannual variations. Plain Language AbstractThe Beaufort Gyre centered in the Canada Basin of the Arctic Ocean is the major reservoir of fresh water in the Arctic. The primary focus of this study is on quantifying variability and trends in liquid (water) and solid (sea ice) freshwater content in this region. The Beaufort Gyre Exploration Program was initiated in 2003 to synthesize results of historical data analysis, design and conduct long-term observations, and to provide information for numerical modeling under the umbrella of the FAMOS (Forum for Arctic Observing and Modeling Synthesis) project. The data collected from research cruises, moorings, Ice-Tethered Profiler observations, and satellite altimetry document an increase of more than 6,400 km 3 of liquid freshwater content from 2003 to 2018, a 40% growth relative to the climatology of the 1970s. This fresh water volume is comparable to the fresh water volume released to the sub-arctic seas during the Great Salinity Anomaly episode of the 1970s. Thus, since the 2000s, the stage has been set for another possible release of fresh water to lower latitudes with accompanying climate impacts, including changes to sea ice conditions, ocean circulation, and ecosystems of the Sub-Arctic similar to the influence of the Great Salinity Anomaly observed in the 1970s.

Signatures of ocean eddies, fronts, and filaments are commonly observed within marginal ice zones (MIZs) from satellite images of sea ice concentration, and in situ observations via ice‐tethered profilers or underice gliders. However, localized and intermittent sea ice heating and advection by ocean eddies are currently not accounted for in climate models and may contribute to their biases and errors in sea ice forecasts. Here, we explore mechanical sea ice interactions with underlying submesoscale ocean turbulence. We demonstrate that the release of potential energy stored in meltwater fronts can lead to energetic submesoscale motions along MIZs with spatial scales O(10 km) and Rossby numbers O(1). In low‐wind conditions, cyclonic eddies and filaments efficiently trap the sea ice and advect it over warmer surface ocean waters where it can effectively melt. The horizontal eddy diffusivity of sea ice mass and heat across the MIZ can reach O(200 m2 s−1). Submesoscale ocean variability also induces large vertical velocities (order 10 m d−1) that can bring relatively warm subsurface waters into the mixed layer. The ocean‐sea ice heat fluxes are localized over cyclonic eddies and filaments reaching about 100 W m−2. We speculate that these submesoscale‐driven intermittent fluxes of heat and sea ice can contribute to the seasonal evolution of MIZs. With the continuing global warming and sea ice thickness reduction in the Arctic Ocean, submesoscale sea ice‐ocean processes are expected to become increasingly prominent.

Recently, the Beaufort Gyre has accumulated over 20,000 km3 of freshwater in response to strong anticyclonic atmospheric winds that have prevailed over the gyre for almost two decades. Here we explore key physical processes affecting the accumulation and release of freshwater within an idealized eddy‐resolving model of the Beaufort Gyre. We demonstrate that a realistic halocline can be achieved when its deepening tendency due to Ekman pumping is counteracted by the cumulative action of mesoscale eddies. Based on this balance, we derive analytical scalings for the depth of the halocline and its spin‐up time scale and emphasize their explicit dependence on eddy dynamics. Our study further suggests that the Beaufort Gyre is currently in a state of high sensitivity to atmospheric winds. However, an intensification of surface stress would inevitably lead to a saturation of the freshwater content—a constraint inherently set by the intricacies of the mesoscale eddy dynamics.

In the Amundsen Sea, modified Circumpolar Deep Water (mCDW) intrudes into ice shelf cavities, causing high ice shelf melting near the ice sheet grounding lines, accelerating ice flow, and controlling the pace of future Antarctic contributions to global sea level. The pathways of mCDW towards grounding lines are crucial as they directly control the heat reaching the ice. A realistic representation of mCDW circulation, however, remains challenging due to the sparsity of in-situ observations and the difficulty of ocean models to reproduce the available observations. In this study, we use an unprecedentedly high-resolution (200 m horizontal and 10 m vertical grid spacing) ocean model that resolves shelf-sea and sub-ice-shelf environments in qualitative agreement with existing observations during austral summer conditions. We demonstrate that the waters reaching the Pine Island and Thwaites grounding lines follow specific, topographically-constrained routes, all passing through a relatively small area located around 104°W and 74.3°S. The temporal and spatial variabilities of ice shelf melt rates are dominantly controlled by the sub-ice shelf ocean current. Our findings highlight the importance of accurate and high-resolution ocean bathymetry and subglacial topography for determining mCDW pathways and ice shelf melt rates.

The halocline of the Beaufort Gyre varies significantly on interannual to decadal time scales, affecting the freshwater content (FWC) of the Arctic Ocean. This study explores the role of eddies in the Ekmandriven gyre variability. Following the transformed Eulerian-mean paradigm, the authors develop a theory that links the FWC variability to the stability of the large-scale gyre, defined as the inverse of its equilibration time. The theory, verified with eddy-resolving numerical simulations, demonstrates that the gyre stability is explicitly controlled by the mesoscale eddy diffusivity. An accurate representation of the halocline dynamics requires the eddy diffusivity of 300 6 200 m 2 s 21 , which is lower than what is used in most low-resolution climate models. In particular, on interannual and longer time scales the eddy fluxes and the Ekman pumping provide equally important contributions to the FWC variability. However, only large-scale Ekman pumping patterns can significantly alter the FWC, with spatially localized perturbations being an order of magnitude less efficient. Lastly, the authors introduce a novel FWC tendency diagnostic-the Gyre Index-that can be conveniently calculated using observations located only along the gyre boundaries. Its strong predictive capabilities, assessed in the eddy-resolving model forced by stochastic winds, suggest that the Gyre Index would be of use in interpreting FWC evolution in observations as well as in numerical models.

The Beaufort Gyre freshwater content has increased since the 1990s, potentially stabilizing in recent years. The mechanisms proposed to explain the stabilization involve either mesoscale eddy activity that opposes Ekman pumping or the reduction of Ekman pumping due to reduced sea ice-ocean surface stress. However, the relative importance of these mechanisms is unclear. Here, we present observational estimates of the Beaufort Gyre mechanical energy budget and show that energy dissipation and freshwater content stabilization by eddies increased in the late-2000s. The loss of sea ice and acceleration of ocean currents after 2007 resulted in enhanced mechanical energy input but without corresponding increases in potential energy storage. To balance the energy surplus, eddy dissipation and its role in gyre stabilization must have increased after 2007. Our results imply that declining Arctic sea ice will lead to an increasingly energetic Beaufort Gyre with eddies playing a greater role in its stabilization.

The Western Arctic Ocean is a host to major ocean circulation systems, many of which generate eddies that can transport water masses and corresponding tracers over long distances from their formation sites. However, comprehensive observations of critical eddy characteristics are currently not available and are limited to spatially and temporally sparse in situ observations. Here we use high-resolution spaceborne synthetic aperture radar measurements to detect eddies from their surface imprints in ice-free sea surface roughness, and in sea ice patterns throughout marginal ice zones. We provide the first estimate of eddy characteristics extending over the seasonally ice-free and marginal ice zone regions of the Western Arctic Ocean, including their locations, diameters, and monthly distribution. Using available synthetic aperture radar data, we identified over 4,000 open ocean eddies, as well as over 3,500 eddies in marginal ice zones from June to October in 2007, 2011, and 2016. Eddies range in size between 0.5 and 100 km and are frequently found over the shelf and near continental slopes but also present in the deep Canada Basin and over the Chukchi Plateau. We find that cyclonic eddies are twice more frequent compared to anticyclonic eddies at the surface, distinct from the dominating anticyclonic eddies observed at depth by in situ moorings and ice-tethered profilers. Our study supports the notion that eddies are ubiquitous in the Western Arctic Ocean even in the presence of sea ice and emphasizes the need for improved ocean observations and modeling at eddy scales. Plain Language SummaryOcean eddies play an important role in the transport of heat, salt, and pollutants over long distances from their formation sites. However, their observations in the Arctic Ocean are complex due to severe weather and sea ice cover. Here we present results of high-resolution satellite observations over the ice-free ocean and in the marginal ice zones. Detailed eddy characteristics are for the first time presented for the Western Arctic Ocean. These results provide observational evidence that eddies are ubiquitous in this Arctic region even in the presence of sea ice and emphasize the need for improved ocean observations and modeling at eddy scales. Key Points: • First account of eddies in the Western Arctic Ocean is presented based on satellite observations over open ocean and marginal ice zones • Eddies range in size between 0.5 and 100 km and are ubiquitous over deep and shelf regions of the Western Arctic Ocean • Cyclonic eddies are twice more frequent compared to anticyclones

Mesoscale eddies shape the Beaufort Gyre response to Ekman pumping, but their transient dynamics are poorly understood. Climate models commonly use the Gent-McWilliams (GM) parameterization, taking the eddy streamfunction c* to be proportional to an isopycnal slope s and an eddy diffusivity K. This local-in-time parameterization leads to exponential equilibration of currents. Here, an idealized, eddy-resolving Beaufort Gyre model is used to demonstrate that c* carries a finite memory of past ocean states, violating a key GM assumption. As a consequence, an equilibrating gyre follows a spiral sink trajectory implying the existence of a damped mode of variability-the eddy memory (EM) mode. The EM mode manifests during the spinup as a 15% overshoot in isopycnal slope (2000 km 3 freshwater content overshoot) and cannot be explained by the GM parameterization. An improved parameterization is developed, such that c* is proportional to an effective isopycnal slope s*, carrying a finite memory g of past slopes. Introducing eddy memory explains the model results and brings to light an oscillation with a period 2p ffiffiffiffiffiffiffiffiffi T E g p ' 50 yr, where the eddy diffusion time scale T E ; 10 yr and g ' 6 yr are diagnosed from the eddy-resolving model. The EM mode increases the Ekman-driven gyre variance by g/T E ' 50% 6 15%, a fraction that stays relatively constant despite both time scales decreasing with increased mean forcing. This study suggests that the EM mode is a general property of rotating turbulent flows and highlights the need for better observational constraints on transient eddy field characteristics.