Publications

Recent evidence from mooring data in the equatorial Atlantic reveals that semiannual and longer time scale ocean current variability is close to being resonant with equatorial basin modes. Here we show that intraseasonal variability, with time scales of tens of days, provides the energy to maintain these resonant basin modes against dissipation. The mechanism is analogous to that by which storm systems in the atmosphere act to maintain the atmospheric jet stream. We demonstrate the mechanism using an idealized model setup that exhibits equatorial deep jets. The results are supported by direct analysis of available mooring data from the equatorial Atlantic Ocean covering a depth range of several thousand meters. The analysis of the mooring data suggests that the same mechanism also helps maintain the seasonal variability.

---------------------------------------------------------------------

Arctic Ocean measurements reveal a near doubling of ocean heat content relative to the freezing temperature in the Beaufort Gyre halocline over the past three decades (1987–2017). This warming is linked to anomalous solar heating of surface waters in the northern Chukchi Sea, a main entryway for halocline waters to join the interior Beaufort Gyre. Summer solar heat absorption by the surface waters has increased fivefold over the same time period, chiefly because of reduced sea ice coverage. It is shown that the solar heating, considered together with subduction rates of surface water in this region, is sufficient to account for the observed halocline warming. Heat absorption at the basin margins and its subsequent accumulation in the ocean interior, therefore, have consequences for Beaufort Gyre sea ice beyond the summer season.

---------------------------------------------------------------------

The Arctic ice cover influences the generation, propagation, and dissipation of internal waves, which in turn may affect vertical mixing in the ocean interior. The Arctic internal wave field and its relationship to the ice cover is investigated using observations from Ice-Tethered Profilers with Velocity and Seaglider sampling during the 2014 Marginal Ice Zone experiment in the Canada Basin. Ice roughness, ice concentration, and wind forcing all influenced the daily to seasonal changes in the internal wave field. Three different ice concentration thresholds appeared to determine the evolution of internal wave spectral energy levels: 1) the initial decrease from 100% ice concentration after which dissipation during the surface reflection was inferred to increase, 2) the transition to 70-80% ice concentration when the local generation of internal waves increased, and 3) the transition to open water that was associated with larger amplitude internal waves. Ice roughness influenced internal wave properties for ice concentrations greater than approximately 70-80%: smoother ice was associated with reduced local internal wave generation. Richardson numbers were rarely supercritical, consistent with weak vertical mixing under all ice concentrations. On decadal timescales, smoother ice may counteract the effects of lower ice concentration on the internal wave field complicating future predictions of internal wave activity and vertical mixing.

The Argo Program has been implemented and sustained for almost two decades, as a global array of about 4000 profiling floats. Argo provides continuous observations of ocean temperature and salinity versus pressure, from the sea surface to 2000 dbar. The successful installation of the Argo array and its innovative data management system arose opportunistically from the combination of great scientific need and technological innovation. Through the data system, Argo provides fundamental physical observations with broad societally-valuable applications, built on the cost-efficient and robust technologies of autonomous profiling floats. Following recent advances in platform and sensor technologies, even greater opportunity exists now than 20 years ago to (i) improve Argo’s global coverage and value beyond the original design, (ii) extend Argo to span the full ocean depth, (iii) add biogeochemical sensors for improved understanding of oceanic cycles of carbon, nutrients, and ecosystems, and (iv) consider experimental sensors that might be included in the future, for example to document the spatial and temporal patterns of ocean mixing. For Core Argo and each of these enhancements, the past, present, and future progression along a path from experimental deployments to regional pilot arrays to global implementation is described. The objective is to create a fully global, top-to-bottom, dynamically complete, and multidisciplinary Argo Program that will integrate seamlessly with satellite and with other in situ elements of the Global Ocean Observing System (Legler et al., 2015). The integrated system will deliver operational reanalysis and forecasting capability, and assessment of the state and variability of the climate system with respect to physical, biogeochemical, and ecosystems parameters. It will enable basic research of unprecedented breadth and magnitude, and a wealth of ocean-education and outreach opportunities.

There is a growing need for operational oceanographic predictions in both the Arctic and Antarctic polar regions. In the former, this is driven by a declining ice cover accompanied by an increase in maritime traffic and exploitation of marine resources. Oceanographic predictions in the Antarctic are also important, both to support Antarctic operations and also to help elucidate processes governing sea ice and ice shelf stability. However, a significant gap exists in the ocean observing system in polar regions, compared to most areas of the global ocean, hindering the reliability of ocean and
sea ice forecasts. This gap can also be seen from the spread in ocean and sea ice reanalyses for polar regions which provide an estimate of their uncertainty. The reduced reliability of polar predictions may affect the quality of various applications including search and rescue, coupling with numerical weather and seasonal predictions, historical reconstructions (reanalysis), aquaculture and environmental management including environmental emergency response. Here, we outline the status of existing near-real time ocean observational efforts in polar regions, discuss gaps, and explore perspectives for the future. Specific recommendations include a renewed call for open access to data, especially real-time data, as a critical capability for improved sea ice and weather forecasting and other environmental prediction needs. Dedicated efforts are also needed to make use of additional observations made as part of the Year of Polar Prediction (YOPP; 2017–2019) to inform optimal observing system design. To provide a polar extension to the Argo network, it is recommended that a network of ice-borne
sea ice and upper-ocean observing buoys be deployed and supported operationally in ice-covered areas together with autonomous profiling floats and gliders (potentially with ice detection capability) in seasonally ice covered seas. Finally, additional efforts to better measure and parameterize surface exchanges in polar regions are much needed
to improve coupled environmental prediction.

--------------------------------------------------------------------------

The history of observing the ocean is reviewed. The evolution of particular classes of ocean measurements (e.g. shipboard, moorings and drifting floats) are summarized along with some of the discoveries and dynamical understanding they made possible. By the 1970s, isolated and "expedition" observational approaches were evolving into experimental campaigns that covered large ocean areas and addressed multi-scale phenomena using diverse instrumental suites and associated modelling and analysis teams. The Mid-Ocean Dynamics Experiment (MODE) addressed mesoscale "eddies" and their interaction with larger-scale currents using new ocean modeling and experiment design techniques and a suite of developing observational methods.  Following MODE, new instrument networks were established to study processes that dominated ocean behavior in different regions. The Tropical Ocean Global Atmosphere program gathered multi-year time series in the tropical Pacific to understand, and eventually predict, evolution of coupled ocean-atmosphere phenomena like El Niño Southern Oscillation (ENSO). The World Ocean Circulation Experiment (WOCE) sought to quantify ocean transport throughout the global ocean using temperature, salinity and other tracer measurements along with fewer direct velocity measurements with floats and moorings. Western and eastern boundary currents attracted comprehensive measurements and various coastal regions, each with its unique scientific and societally important phenomena, became home to regional observing systems. Today the trend toward networked observing arrays of many instrument types continues to be a productive way to understand and predict large-scale ocean phenomena.

---------------------------------------------------------------------

Biogenic matter characteristics and radiocarbon contents of organic carbon (OC) were examined on sinking particle samples intercepted at three nominal depths of 1000 m, 2000 m, and 3000 m (~50 m above the seafloor) during a 3-year sediment trap program on the New England slope in the Northwest Atlantic. We have sought to characterize the sources of sinking particles in the context of vertical export of biogenic particles from the overlying water column and lateral supply of resuspended sediment particles from adjacent margin sediments. High aluminum (Al) abundances and low OC radiocarbon contents indicated contributions from resuspended sediment which was greatest at 3000 m but also significant at shallower depths. The benthic source (i.e., laterally supplied resuspended sediment) of opal appears negligible based on the absence of a correlation with Al fluxes. In comparison, CaCO3fluxes at 3000 m showed a positive correlation with Al fluxes.  Benthic sources accounted for 42 ~ 63 % of the sinking particle flux based on radiocarbon mass balance and the relationship between Al flux and CaCO3flux. Episodic pulses of Al at 3000 m were significantly correlated with the near-bottom current at a nearby hydrographic mooring site, implying the importance of current variability in lateral particle transport. However, Al fluxes at 1000 m and 2000 m were coherent but differed from those at 3000 m, implying more than one mode of lateral supply of particles in the water column.

---------------------------------------------------------------------

Continuous observations in the Canada Basin from March through December 2014 are 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 velocity in combination with turbulent momentum flux measurements from four clusters of drifting instrument systems are complemented by direct drill-hole measurements and IceBridge over flights in March, as well as remote sensing imagery about the instrument clusters. Observations captured near-complete seasonal melting of the ice cover. Spatially, 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%. In July and August, indirect indicators of ice-ocean drag and roughness suggested increased momentum transfer in low ice concentration conditions. The ice-ocean drag coefficient parameterization was found to be invalid in September with low ice concentrations and small ice floe sizes. The differing and temporally-variable nearby ice topography affected mixed layer currents, including the ocean’s response to storm events. The kinetic energy budget of the ice showed that wind work and ocean work on the ice were the dominant terms controlling ice motion throughout the entire melt season. These seasonal changes in the air-ice-ocean system show that ice topography, ice concentration, and the shallow summer mixed layer all impact the ocean’s response to forcing events.

---------------------------------------------------------------------

A moored array spanning the continental slope southeast of Cape Cod sampled the equatorward-flowing Deep Western Boundary Current (DWBC) for a 10-year period: May 2004 - May 2014.  Daily profiles of subinertial velocity, temperature, salinity and neutral density are constructed for each mooring site and cross-line DWBC transport time series are derived for specified water mass layers.  Time-averaged transports based on daily estimates of the flow and density fields in stream coordinates are contrasted with those derived from the Eulerian-mean flow field, modes of DWBC transport variability are investigated through compositing, and comparisons are made to transport estimates for other latitudes.   Integrating the daily velocity estimates over the neutral density range of 27.8 - 28.125 kg/m³ (encompassing Labrador Sea and Overflow Water layers), a mean equatorward DWBC transport of 22.8 x 10⁶ m³/s ± 1.9 x 10⁶ m³/s is obtained.  Notably, a statistically-significant trend of decreasing equatorward transport is observed in several of the DWBC components as well as the current as a whole.  The largest linear change (a 4% decrease per year) is seen in the layer of Labrador Sea Waters renewed in the early 1990s whose transport fell from 9.0 x 10⁶ m³/s at the beginning of the field program to 5.8 x 10⁶ m³/s at its end.  The corresponding linear fit to the combined Labrador Sea and Overflow Water DWBC transport decreases from 26.4 x 10⁶ m³/s to 19.1 x 10⁶ m³/s.  In contrast, no long-term trend is observed in upper-ocean Slope Water transport.  These trends are discussed in the context of decadal observations of the North Atlantic circulation, and subpolar air-sea interaction/water mass transformation.

---------------------------------------------------------------------

Observations of the Deep Western Boundary Current (DWBC) at Line W on the western North Atlantic continental slope southeast of Cape Cod from 1995 to 2014 reveal water mass changes that are consistent with changes in source water properties upstream in the Labrador Sea. This is most evident in the cold, dense, and deep class of  Labrador Sea Water (dLSW) that was created and progressively replenished and deepened by recurring winter convection during the severe winters of 1987-1994. The arrival of this record cold, fresh and low potential vorticity anomaly at Line W lags it formation in the Labrador Sea by 3-7 years. Complementary observations along path of the DWBC provide further evidence that this anomaly is advected along the boundary and indicate that stirring between the boundary and the interior intensifies south of the Flemish Cap.  Finally, the consistency of the data with realistic advective and mixing time scales is assessed using the Waugh and Hall [2005] model framework. The data are found to be best represented by a mean transit time of 5 years from the Labrador Sea to Line W, with a  leading order role for both advection by the DWBC and mixing between the boundary flow and interior waters.

---------------------------------------------------------------------

Shipboard velocity and water property data from 18 transects across the North Atlantic Deep Western Boundary Current (DWBC) near 40 °N are examined to study the evolution of the Denmark Strait Overflow Water (DSOW) component of the DWBC and mixing between DSOW and the interior. The examined transects along Line W – which stretches from the continental shelf south of New England to Bermuda – were made between 1994 and 2014. The shipboard data comprise measurements at regular stations of velocity from lowered acoustic Doppler current profilers, CTD profiles and trace gas chlorofluorocarbon (CFC) concentrations from bottle samples at discrete depths. Comparison of the Line W velocity sections with concurrent sea surface height maps from satellite altimetry indicates that large cyclones in the deep ocean accompany intermittent quasi-stationary meander troughs in the Gulf Stream path at Line W. A composite of 5 velocity sections along Line W suggests that a typical cyclone reaches swirl speeds of greater than 30 cm s⁻¹ at 3400-m depth and has a radius (distance between the center and the maximum velocity) of 75 km. Tracer data suggest that these cyclones affect not only the deep velocity structure along Line W, but also provide a mechanism for water exchange between the DWBC's DSOW and the interior. Vigorous exchange is corroborated by a mismatch in the CFC-11:CFC-12 and CFC-113:CFC-12 ratio ages calculated for DSOW at Line W. During the most recent 5-year period (2010–2014), a decrease in DSOW density has been driven by warming (increasing by almost 0.1 °C) as salinity has increased only slightly (by 0.003, which is close to the 0.002 uncertainty of the measurements). The abyssal ocean offshore of the DWBC and Gulf Stream and deeper than 3000-m depth has freshened at a rate of 6×10⁻⁴ yr⁻¹ since at least 2003. Density here remains nearly unchanged over this period, due to temperature compensation, though a linear cooling trend in the abyssal ocean (to compensate the freshening) is not statistically significant.

------------------------------------------------------------------------

Direct measurements of oceanic turbulent parameters were taken upstream of and across Drake Passage, in the region of the Subantarctic and Polar Fronts. Values of turbulent kinetic energy dissipation rate ε estimated by microstructure are up to two orders of magnitude lower than previously published estimates in the upper 1000 m. Turbulence levels in Drake Passage are systematically higher than values upstream, regardless of season. The dissipation of thermal variance χ is enhanced at middepth throughout the surveys, with the highest values found in northern Drake Passage, where water mass variability is the most pronounced. Using the density ratio, evidence for double-diffusive instability is presented. Subject to double-diffusive physics, the estimates of diffusivity using the Osborn–Cox method are larger than ensemble statistics based on ε and the buoyancy frequency.

-------------------------------------------------

Seasonal variability of the tropical Atlantic circulation is dominated by the annual cycle, but semiannual variability is also pronounced, despite weak forcing at that period. This study uses multiyear, full-depth velocity measurements from the central equatorial Atlantic to analyze the vertical structure of annual and semiannual variations of zonal velocity. A baroclinic modal decomposition finds that the annual cycle is dominated by the fourth mode and the semiannual cycle is dominated by the second mode. Similar local behavior is found in a high-resolution general circulation model. This simulation reveals that the annual and semiannual cycles of the respective dominant baroclinic modes are associated with characteristic basinwide structures. Using an idealized, linear, reduced-gravity model to simulate the dynamics of individual baroclinic modes, it is shown that the observed circulation variability can be explained by resonant equatorial basin modes. Corollary simulations of the reduced-gravity model with varying basin geometry (i.e., square basin vs realistic coastlines) or forcing (i.e., spatially uniform vs spatially variable wind) show a structural robustness of the simulated basin modes. A main focus of this study is the seasonal variability of the Equatorial Undercurrent (EUC) as identified in recent observational studies. Main characteristics of the observed EUC including seasonal variability of transport, core depth, and maximum core velocity can be explained by the linear superposition of the dominant equatorial basin modes as obtained from the reduced-gravity model.

--------------------------------------------------

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⁻² in this study). Buoyancy provided by MP drainage shoals and freshens the IOBL resulting in a 39 MJ m⁻² 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⁻²) and enhanced basal ice melt (3–6 cm d⁻¹), 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.

----------------------------------------------------

The equatorial deep jets (EDJs) are a ubiquitous feature of the equatorial oceans; in the Atlantic Ocean, they are the dominant mode of interannual variability of the zonal flow at intermediate depth. On the basis of more than 10 years of moored observations of zonal velocity at 23°W, the vertically propagating EDJs are best described as superimposed oscillations of the 13th to the 23rd baroclinic modes with a dominant oscillation period for all modes of 1650 days. This period is close to the resonance period of the respective gravest equatorial basin mode for the dominant vertical modes 16 and 17. It is argued that since the equatorial basin mode is composed of linear equatorial waves, a linear reduced-gravity model can be employed for each baroclinic mode, driven by spatially homogeneous zonal forcing oscillating with the EDJ period. The fit of the model solutions to observations at 23°W yields a basinwide reconstruction of the EDJs and the associated vertical structure of their forcing. From the resulting vertical profile of mean power input and vertical energy flux on the equator, it follows that the EDJs are locally maintained over a considerable depth range, from 500 to 2500 m, with the maximum power input and vertical energy flux at 1300 m. The strong dissipation closely ties the apparent vertical propagation of energy to the vertical distribution of power input and, together with the EDJs’ prevailing downward phase propagation, requires the phase of the forcing of the EDJs to propagate downward.

-----------------------------------------------------------------------

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.

---------------------------------------------------------------------

Understanding the physical and biogeochemical processes that control CO₂ and dissolved oxygen (DO) dynamics in the Arctic Ocean (AO) is crucial for predicting future air-sea CO₂ fluxes and ocean acidification. Past studies have primarily been conducted on the AO continental shelves during low-ice periods and we lack information on gas dynamics in the deep AO basins where ice typically inhibits contact with the atmosphere. To study these gas dynamics, in situ time-series data have been collected in the Canada Basin during late summer to autumn of 2012. Partial pressure of CO₂ (pCO₂), DO concentration, temperature, salinity, and chlorophyll-a fluorescence (Chl-a) were measured in the upper ocean in a range of sea ice states by two drifting instrument systems. Although the two systems were on average only 222 km apart, they experienced considerably different ice cover and external forcings during the 40–50 day periods when data were collected. The pCO₂ levels at both locations were well below atmospheric saturation whereas DO was almost always slightly supersaturated. Modeling results suggest that air-sea gas exchange, net community production (NCP), and horizontal gradients were the main sources of pCO₂ and DO variability in the sparsely ice-covered AO. In areas more densely covered by sea ice, horizontal gradients were the dominant source of variability, with no significant NCP in the surface mixed layer. If the AO reaches equilibrium with atmospheric CO₂ as ice cover continues to decrease, aragonite saturation will drop from a present mean of 1.00 ± 0.02 to 0.86 ± 0.01.

---------------------------------------------------------------------------

The Arctic Ocean thermohaline stratification frequently exhibits a staircase structure overlying the Atlantic Water Layer that can be attributed to the diffusive form of double-diffusive convection. The staircase consists of multiple layers of O(1) m in thickness separated by sharp interfaces, across which temperature and salinity change abruptly. Through a detailed analysis of Ice-Tethered Profiler measurements from 2004 to 2013, the double-diffusive staircase structure is characterized across the entire Arctic Ocean. We demonstrate how the large-scale Arctic Ocean circulation influences the small-scale staircase properties. These staircase properties (layer thicknesses and temperature and salinity jumps across interfaces) are examined in relation to a bulk vertical density ratio spanning the staircase stratification. We show that the Lomonosov Ridge serves as an approximate boundary between regions of low density ratio (approximately 3–4) on the Eurasian side and higher density ratio (approximately 6–7) on the Canadian side. We find that the Eurasian Basin staircase is characterized by fewer, thinner layers than that in the Canadian Basin, although the margins of all basins are characterized by relatively thin layers and the absence of a well-defined staircase. A double-diffusive 4/3 flux law parametrization is used to estimate vertical heat fluxes in the Canadian Basin to be O(0.1) W m⁻². It is shown that the 4/3 flux law may not be an appropriate representation of heat fluxes through the Eurasian Basin staircase. Here molecular heat fluxes are estimated to be between O(0.01) and O(0.1) W m⁻². However, many uncertainties remain about the exact nature of these fluxes.

-------------------------------------------------------------------

Interannual variability in the volumetric water mass distribution within the North Atlantic Subtropical Gyre is described in relation to variability in the Atlantic meridional overturning circulation. The relative roles of diabatic and adiabatic processes in the volume and heat budgets of the subtropical gyre are investigated by projecting data into temperature coordinates as volumes of water using an Argo-based climatology and an ocean state estimate (ECCO version 4). This highlights that variations in the subtropical gyre volume budget are predominantly set by transport divergence in the gyre. A strong correlation between the volume anomaly due to transport divergence and the variability of both thermocline depth and Ekman pumping over the gyre suggests that wind-driven heave drives transport anomalies at the gyre boundaries. This wind-driven heaving contributes significantly to variations in the heat content of the gyre, as do anomalies in the air–sea fluxes. The analysis presented suggests that wind forcing plays an important role in driving interannual variability in the Atlantic meridional overturning circulation and that this variability can be unraveled from spatially distributed hydrographic observations using the framework presented here.

----------------------------------------------------------

Eight Ice-Tethered Profilers were deployed in the Arctic Ocean between 2011 and 2013 to measure vertical distributions of photosynthetically active radiation (PAR) and other bio-optical properties in ice-covered water columns, multiple times a day over periods of up to a year. With the radiometers used on these profilers, PAR could be measured to depths of only ∼20–40 m in the central Arctic in late summer under sea ice ∼1 m thick. At lower latitudes in the Beaufort Gyre, late summer PAR was measurable under ice to depths exceeding 125 m. The maximum depths of measurable PAR followed seasonal trends in insolation, with isolumes shoaling in the fall as solar elevation decreased and deepening in spring and early summer after insolation resumed and sea ice diminished. PAR intensities were often anomalously low above 20 m, likely due to a shading effect associated with local horizontal heterogeneity in light transmittance by the overlying sea ice. A model was developed to parameterize these complex vertical PAR distributions to improve estimates of the water column diffuse attenuation coefficient and other related parameters. Such a model is necessary to separate the effect of surface ice heterogeneity on under-ice PAR profiles from that of the water column itself, so that euphotic zone depth in ice-covered water columns can be computed using canonical metrics such as the 1% light level. Water column diffuse attenuation coefficients derived from such autonomously-collected PAR profile data, using this model, agreed favorably with values determined manually in complementary studies.

--------------------------------------------------------------------

The loss of Arctic sea ice has emerged as a leading signal of global warming. This, together with acknowledged impacts on other components of the Earth system, has led to the term “the new Arctic.” Global coupled climate models predict that ice loss will continue through the twenty-first century, with implications for governance, economics, security, and global weather. A wide range in model projections reflects the complex, highly coupled interactions between the polar atmosphere, ocean, and cryosphere, including teleconnections to lower latitudes. This paper summarizes our present understanding of how heat reaches the ice base from the original sources—inflows of Atlantic and Pacific Water, river discharge, and summer sensible heat and shortwave radiative fluxes at the ocean/ice surface—and speculates on how such processes may change in the new Arctic. The complexity of the coupled Arctic system, and the logistic and technological challenges of working in the Arctic Ocean, require a coordinated interdisciplinary and international program that will not only improve understanding of this critical component of global climate but will also provide opportunities to develop human resources with the skills required to tackle related problems in complex climate systems. We propose a research strategy with components that include 1) improved mapping of the upper- and middepth Arctic Ocean, 2) enhanced quantification of important process, 3) expanded long-term monitoring at key heat-flux locations, and 4) development of numerical capabilities that focus on parameterization of heat-flux mechanisms and their interactions.

----------------------------------------------------------------------

Oceanic internal waves are closely linked to turbulence. Here a relationship between vertical wave number (k_z) spectra of fine-scale vertical kinetic energy (VKE) and turbulent dissipation ε is presented using more than 250 joint profiles from five diverse dynamic regimes, spanning latitudes between the equator and 60°. In the majority of the spectra VKE varies as . Scaling VKE with collapses the off-equatorial spectra to within but underestimates the equatorial spectrum. The simple empirical relationship between VKE and ε fits the data better than a common shear-and-strain fine-scale parameterization, which significantly underestimates ε in the two data sets that are least consistent with the Garrett-Munk (GM) model. The new relationship between fine-scale VKE and dissipation rate can be interpreted as an alternative, single-parameter scaling for turbulent dissipation in terms of fine-scale internal wave vertical velocity that requires no reference to the GM model spectrum.

----------------------------------------------------------