Meltwater Intrusions Reveal Mechanisms for Rapid Submarine Melt at a Tidewater Glacier

Submarine melting has been implicated as a driver of glacier retreat and sea level rise, but to date melting has been difficult to observe and quantify. As a result, melt rates have been estimated from parameterizations that are largely unconstrained by observations, particularly at the near-vertical termini of tidewater glaciers. With standard coefficients, these melt parameterizations predict that ambient melting (the melt away from subglacial discharge outlets) is negligible compared to discharge-driven melting for typical tidewater glaciers. Here, we present new data from LeConte Glacier, Alaska, that challenges this paradigm. Using autonomous kayaks, we observe ambient meltwater intrusions that are ubiquitous within 400 m of the terminus, and we provide the first characterization of their properties, structure, and distribution. Our results suggest that ambient melt rates are substantially higher (×100) than standard theory predicts and that ambient melting is a significant part of the total submarine melt flux. We explore modifications to the prevalent melt parameterization to provide a path forward for improved modeling of ocean-glacier interactions.This work was funded by NSF OPP Grants 1503910, 1504191, 1504288, and 1504521 and National Geographic Grant CP4-171R-17. Additionally, this research was supported by the NOAA Climate and Global Change Postdoctoral Fellowship Program, administered by UCAR’s Cooperative Programs for the Advancement of Earth System Science (CPAESS) under award #NA18NWS4620043B. These observations would not be possible without the skilled engineering team who developed the autonomous kayaks—including Jasmine Nahorniak, June Marion, Nick McComb, Anthony Grana, and Corwin Perren—and also the Captain and crew of the M/V Amber Anne. We thank Donald Slater and an anonymous reviewer for valuable feedback that improved this manuscript. Data availability: All of the oceanographic data collected by ship and kayak have been archived with the National Centers for Environmental Information (Accession 0189574, https://accession.nodc.noaa.gov/ 0189574). The glacier data have been archived at the Arctic Data Center (https://doi.org/10.18739/A22G44).Ye

A “cold path” for the Gulf Stream–troposphere connection

The mechanism by which the Gulf Stream sea surface temperature (SST) front anchors a band of precipitation on its warm edge is still a matter of debate, and little is known about how synoptic activity contributes to the mean state. In the present study, the influence of the SST front on precipitation is investigated during the course of a single extratropical cyclone using a regional configuration of the Met Office Unified Model. The comparison of a control run with a simulation in which SST gradients were smoothed brought the following conclusions: a band of precipitation is reproduced for a single extratropical cyclone, and the response to the SST gradient is dominated by a change of convective precipitation in the cold sector of the storm. Several climatological features described by previous studies, such as surface wind convergence on the warm edge or a meridional circulation cell across the SST front, are also reproduced at synoptic time scales in the cold sector. Based on these results, a simple boundary layer model is proposed to explain the convective and dynamical response to the SST gradient in the cold sector. In this model, cold and dry air parcels acquire more buoyancy over a sharp SST gradient and become more convectively unstable. The convection sets a pressure anomaly over the entire depth of the boundary layer that drives wind convergence. This case study offers a new pathway by which the SST gradient can anchor a climatological band of precipitation

Observations and numerical simulations of large-eddy circulation in the ocean surface mixed layer

Author Posting. © American Geophysical Union, 2014. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 41 (2014): 7584–7590, doi:10.1002/2014GL061637.Two near-surface dye releases were mapped on scales of minutes to hours temporally, meters to order 1 km horizontally, and 1–20 m vertically using a scanning, depth-resolving airborne lidar. In both cases, dye evolved into a series of rolls with their major axes approximately aligned with the wind and/or near-surface current. In both cases, roll spacing was also of order 5–10 times the mixed layer depth, considerably larger than the 1–2 aspect ratio expected for Langmuir cells. Numerical large-eddy simulations under similar forcing showed similar features, even without Stokes drift forcing. In one case, inertial shear driven by light winds induced large aspect ratio large-eddy circulation. In the second, a preexisting lateral mixed layer density gradient provided the dominant forcing. In both cases, the growth of the large-eddy structures and the strength of the resulting dispersion were highly dependent on the type of forcing.Support for the 2004 field experiment was provided by the Cecil H. and Ida M. Green Technology Innovation Fund and Coastal Ocean Institute grant 27001545, both through Woods Hole Oceanographic Institution, and by Office of Naval Research grant N00014-01-1-0984. Support for the 2011 field experiments was provided by ONR grants N00014-09-1-0194, N00014-09-1-0175, N00014-11-WX-21010, N00014-12-WX-21031, and N00014-09-1-0460 and NSF grants OCE-0751734 and OCE-0751653. Simulations were supported under grant N00014-09-1-0268.2015-05-0

The refreezing of melt ponds on Arctic sea ice

The presence of melt ponds on the surface of Arctic sea ice significantly reduces its albedo, inducing a positive feedback leading to sea ice thinning. While the role of melt ponds in enhancing the summer melt of sea ice is well known, their impact on suppressing winter freezing of sea ice has, hitherto, received less attention. Melt ponds freeze by forming an ice lid at the upper surface, which insulates them from the atmosphere and traps pond water between the underlying sea ice and the ice lid. The pond water is a store of latent heat, which is released during refreezing. Until a pond freezes completely, there can be minimal ice growth at the base of the underlying sea ice. In this work, we present a model of the refreezing of a melt pond that includes the heat and salt balances in the ice lid, trapped pond, and underlying sea ice. The model uses a two-stream radiation model to account for radiative scattering at phase boundaries. Simulations and related sensitivity studies suggest that trapped pond water may survive for over a month. We focus on the role that pond salinity has on delaying the refreezing process and retarding basal sea ice growth. We estimate that for a typical sea ice pond coverage in autumn, excluding the impact of trapped ponds in models overestimates ice growth by up to 265 million km3, an overestimate of 26%

Toward Regional Characterizations of the Oceanic Internal Wavefield

Many major oceanographic internal wave observational programs of the last 4 decades are reanalyzed in order to characterize variability of the deep ocean internal wavefield. The observations are discussed in the context of the universal spectral model proposed by Garrett and Munk. The Garrett and Munk model is a good description of wintertime conditions at Site-D on the continental rise north of the Gulf Stream. Elsewhere and at other times, significant deviations in terms of amplitude, separability of the 2-D vertical wavenumber - frequency spectrum, and departure from the model's functional form are noted. Subtle geographic patterns are apparent in deviations from the high frequency and high vertical wavenumber power laws of the Garrett and Munk spectrum. Moreover, such deviations tend to co-vary: whiter frequency spectra are partnered with redder vertical wavenumber spectra. Attempts are made to interpret the variability in terms of the interplay between generation, propagation and nonlinearity using a statistical radiative balance equation. This process frames major questions for future research with the insight that such integrative studies could constrain both observationally and theoretically based interpretations

Langmuir turbulence and surface heating in the ocean surface boundary layer

This study uses large-eddy simulation to investigate the structure of the ocean surface boundary layer (OSBL) in the presence of Langmuir turbulence and stabilizing surface heat fluxes. The OSBL consists of a weakly stratified layer, despite a surface heat flux, above a stratified thermocline. The weakly stratified (mixed) layer is maintained by a combination of a turbulent heat flux produced by the wave-driven Stokes drift and downgradient turbulent diffusion. The scaling of turbulence statistics, such as dissipation and vertical velocity variance, is only affected by the surface heat flux through changes in the mixed layer depth. Diagnostic models are proposed for the equilibrium boundary layer and mixed layer depths in the presence of surface heating. The models are a function of the initial mixed layer depth before heating is imposed and the Langmuir stability length. In the presence of radiative heating, the models are extended to account for the depth profile of the heating

Impact of melt ponds on Arctic sea ice simulations from 1990 to 2007

The extent and thickness of the Arctic sea ice cover has decreased dramatically in the past few decades with minima in sea ice extent in September 2007 and 2011 and climate models did not predict this decline. One of the processes poorly represented in sea ice models is the formation and evolution of melt ponds. Melt ponds form on Arctic sea ice during the melting season and their presence affects the heat and mass balances of the ice cover, mainly by decreasing the value of the surface albedo by up to 20%. We have developed a melt pond model suitable for forecasting the presence of melt ponds based on sea ice conditions. This model has been incorporated into the Los Alamos CICE sea ice model, the sea ice component of several IPCC climate models. Simulations for the period 1990 to 2007 are in good agreement with observed ice concentration. In comparison to simulations without ponds, the September ice volume is nearly 40% lower. Sensitivity studies within the range of uncertainty reveal that, of the parameters pertinent to the present melt pond parameterization and for our prescribed atmospheric and oceanic forcing, variations of optical properties and the amount of snowfall have the strongest impact on sea ice extent and volume. We conclude that melt ponds will play an increasingly important role in the melting of the Arctic ice cover and their incorporation in the sea ice component of Global Circulation Models is essential for accurate future sea ice forecasts

A new expression for the form stress term in the vertically Lagrangian mean framework for the effect of surface waves on the upper ocean circulation

There is an ongoing discussion in the community concerning the wave-averaged momentum equations in the hybrid vertically Lagrangian and horizontally Eulerian (VL) framework and, in particular, the form stress term (representing the residual effect of pressure perturbations) which is thought to restrict the handling of higher order waves in terms of a perturbation expansion. The present study shows that the traditional pressure-based form stress term can be transformed into a set of terms that do not contain any pressure quantities but do contain the time derivative of a wave-induced velocity. This wave-induced velocity is referred to as the pseudomomentum in the VL framework, as it is analogous to the generalized pseudomomentum in Andrews and McIntyre. This enables the second expression for the wave-averaged momentum equations in the VL framework (this time for the development of the total transport velocity minus the VL pseudomomentum) to be derived together with the vortex force. The velocity-based expression of the form stress term also contains the residual effect of the turbulent viscosity, which is useful for understanding the dissipation of wave energy leading to transfer of momentum from waves to circulation. It is found that the concept of the virtual wave stress of Longuet-Higgins is applicable to quite general situations: it does not matter whether there is wind forcing or not, the waves can have slow variations, and the viscosity coefficient can vary in the vertical. These results provide a basis for revisiting the surface boundary condition used in numerical circulation models

The LatMix summer campaign : submesoscale stirring in the upper ocean

Author Posting. © American Meteorological Society, 2015. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Bulletin of the American Meteorological Society 96 (2015): 1257–1279, doi:10.1175/BAMS-D-14-00015.1.Lateral stirring is a basic oceanographic phenomenon affecting the distribution of physical, chemical, and biological fields. Eddy stirring at scales on the order of 100 km (the mesoscale) is fairly well understood and explicitly represented in modern eddy-resolving numerical models of global ocean circulation. The same cannot be said for smaller-scale stirring processes. Here, the authors describe a major oceanographic field experiment aimed at observing and understanding the processes responsible for stirring at scales of 0.1–10 km. Stirring processes of varying intensity were studied in the Sargasso Sea eddy field approximately 250 km southeast of Cape Hatteras. Lateral variability of water-mass properties, the distribution of microscale turbulence, and the evolution of several patches of inert dye were studied with an array of shipboard, autonomous, and airborne instruments. Observations were made at two sites, characterized by weak and moderate background mesoscale straining, to contrast different regimes of lateral stirring. Analyses to date suggest that, in both cases, the lateral dispersion of natural and deliberately released tracers was O(1) m2 s–1 as found elsewhere, which is faster than might be expected from traditional shear dispersion by persistent mesoscale flow and linear internal waves. These findings point to the possible importance of kilometer-scale stirring by submesoscale eddies and nonlinear internal-wave processes or the need to modify the traditional shear-dispersion paradigm to include higher-order effects. A unique aspect of the Scalable Lateral Mixing and Coherent Turbulence (LatMix) field experiment is the combination of direct measurements of dye dispersion with the concurrent multiscale hydrographic and turbulence observations, enabling evaluation of the underlying mechanisms responsible for the observed dispersion at a new level.The bulk of this work was funded under the Scalable Lateral Mixing and Coherent Turbulence Departmental Research Initiative and the Physical Oceanography Program. The dye experiments were supported jointly by the Office of Naval Research and the National Science Foundation Physical Oceanography Program (Grants OCE-0751653 and OCE-0751734).2016-02-0