Abstract: On the continental shelf around Greenland, warm-salty Atlantic water at depth fills the deep narrow fjords where Greenland's tidewater glaciers terminate. Changes in the quantity or properties of this water mass starting in the mid 1990s is thought to be largely responsible for increased ocean-driven melting of the Greenland Ice Sheet. Using high-resolution (nominal 0.1-degree) ocean circulation models we cannot accurately resolve small-scale processes on the shelf or within fjords. However, we can assess changes in the flux of heat via Atlantic water onto the continental shelf. To understand the causes of the anomalous heat that has reached the shelf we examine heat content of subtropical gyre water and shifts in the North Atlantic and Atlantic Multidecadal Oscillations.
We compare changes in heat transport in two eddy permitting simulations: a global 0.1 degree (5-7km around Greenland) resolution coupled hindcast (1970-2009) simulation of the Parallel Ocean Program (POP) and a regional 0.08 degree (3-5km around Greenland) resolution coupled HYbrid Coordinate Ocean Model (HYCOM) hindcast (1993-2016) simulation. Both models are coupled to the Los Alamos National Laboratory Community Ice CodE version 4 and forced by atmospheric reanalysis fluxes. In both models we look for processes that could explain the increase in heat; processes that are present in both are likely to be robust causes of warming.

Abstract: Concurrent measurements by surface drifters of different configurations and HF radar reveal substantial differences in estimates of the near-surface seawater velocity. On average, speeds of small ultra-thin (5 cm) drifters are significantly greater than co-located drifters with a traditional shallow drogue design, while velocity measurements from the drogued drifters closely match HF radar velocity estimates. Analysis of directional wave spectra measurements from a nearby buoy reveals that Stokes drift accounts for much of the difference between the velocity measurements from the drogued drifters and the ultra-thin drifters, except during times of wave breaking. Under wave breaking conditions, the difference between the ultra-thin drifter velocity and the drogued drifter velocity is much less than the computed Stokes drift. The results suggest that surface currents measured by more common approaches or simulated in models may underrepresent the velocity at the very surface of the ocean that is important for determining momentum and enthalpy fluxes between the ocean and atmosphere and for estimating transport of material at the ocean surface. However, simply adding an estimate of Stokes drift may also not be an appropriate method for estimating the true surface velocity from models or measurements from drogued drifters or HF radar under all sea conditions.

Abstract: We revisit the challenges and prospects for ocean circulation models following Griffies et al. (2010). Over the past decade, ocean circulation models evolved through improved understanding, numerics, spatial discretization, grid configurations, parameterizations, data assimilation, environmental monitoring, and process-level observations and modeling. Important large scale applications over the last decade are simulations of the Southern Ocean, the Meridional Overturning Circulation and its variability, and regional sea level change. Submesoscale variability is now routinely resolved in process models and permitted in a few global models, and submesoscale effects are parameterized in most global models. The scales where nonhydrostatic effects become important are beginning to be resolved in regional and process models. Coupling to sea ice, ice shelves, and high-resolution atmospheric models has stimulated new ideas and driven improvements in numerics. Observations have provided insight into turbulence and mixing around the globe and its consequences are assessed through perturbed physics models. Relatedly, parameterizations of the mixing and overturning processes in boundary layers and the ocean interior have improved. New diagnostics being used for evaluating models alongside present and novel observations are briefly referenced. The overall goal is summarizing new developments in ocean modeling, including how new and existing observations can be used, what modeling challenges remain, and how simulations can be used to support observations.

Abstract: Zooplankton are important consumers of sinking particles in the ocean's twilight zone. However, the impact of different taxa depends on their feeding mode. In contrast to typical suspension-feeding zooplankton, flux-feeding taxa preferentially consume rapidly sinking particles that would otherwise penetrate into the deep ocean. To quantify the potential impact of two flux-feeding zooplankton taxa [Aulosphaeridae (Rhizaria), and Limacina helicina (euthecosome pteropod)] and the total suspension-feeding zooplankton community, we measured depth-stratified abundances of these organisms during six cruises in the California Current Ecosystem. Using allometric-scaling relationships, we computed the percentage of carbon flux intercepted by flux feeders and suspension feeders. These estimates were compared to direct measurements of carbon flux attenuation (CFA) made using drifting sediment traps and U-238-Th-234 disequilibrium. We found that CFA in the shallow twilight zone typically ranged from 500 to 1000 m mol organic C flux remineralized per 10-m vertical depth bin. This equated to approximately 6-10% of carbon flux remineralized/10 m. The two flux-feeding taxa considered in this study could account for a substantial proportion of this flux near the base of the euphotic zone. The mean flux attenuation attributable to Aulosphaeridae was 0.69%/10 m (median = 0.21%/10 m, interquartile range = 0.04-0.81%) at their depth of maximum abundance (similar to 100 m), which would equate to similar to 10% of total flux attenuation in this depth range. The maximum flux attenuation attributable to Aulosphaeridae reached 4.2%/10 m when these protists were most abundant. L. helicina, meanwhile, could intercept 0.45-1.6% of carbon flux/10 m, which was slightly greater (on average) than the Aulosphaeridae. In contrast, suspension-feeding zooplankton in the mesopelagic (including copepods, euphausiids, appendicularians, and ostracods) had combined clearance rates of 2-81 L m(-3) day(-1) (mean of 19.6 L m(-3) day(-1)). This implies a substantial impact on slowly sinking particles, but a negligible impact on the presumably rapidly sinking fecal pellets that comprised the majority of the material collected in sediment traps. Our results highlight the need for a greater research focus on the many taxa that potentially act as flux feeders in the oceanic twilight zone.

Abstract: Turbulent and radiative exchanges of heat between the ocean and atmosphere (hereafter heat fluxes), ocean surface wind stress, and state variables used to estimate them, are Essential Ocean Variables (EOVs) and Essential Climate Variables (ECVs) influencing weather and climate. This paper describes an observational strategy for producing 3-hourly, 25-km (and an aspirational goal of hourly at 10-km) heat flux and wind stress fields over the global, ice-free ocean with breakthrough 1-day random uncertainty of 15 W m–2 and a bias of less than 5 W m–2. At present this accuracy target is met only for OceanSITES reference station moorings and research vessels (RVs) that follow best practices. To meet these targets globally, in the next decade, satellite-based observations must be optimized for boundary layer measurements of air temperature, humidity, sea surface temperature, and ocean wind stress. In order to tune and validate these satellite measurements, a complementary global in situ flux array, built around an expanded OceanSITES network of time series reference station moorings, is also needed. The array would include 500–1000 measurement platforms, including autonomous surface vehicles, moored and drifting buoys, RVs, the existing OceanSITES network of 22 flux sites, and new OceanSITES expanded in 19 key regions. This array would be globally distributed, with 1–3 measurement platforms in each nominal 10° by 10° box. These improved moisture and temperature profiles and surface data, if assimilated into Numerical Weather Prediction (NWP) models, would lead to better representation of cloud formation processes, improving state variables and surface radiative and turbulent fluxes from these models. The in situ flux array provides globally distributed measurements and metrics for satellite algorithm development, product validation, and for improving satellite-based, NWP and blended flux products. In addition, some of these flux platforms will also measure direct turbulent fluxes, which can be used to improve algorithms for computation of air-sea exchange of heat and momentum in flux products and models. With these improved air-sea fluxes, the ocean’s influence on the atmosphere will be better quantified and lead to improved long-term weather forecasts, seasonal-interannual-decadal climate predictions, and regional climate projections.