ALBERT

Description: The transport of bacteria in turbulent river-like environments is addressed, where bacterial populations are frequently encountered attached to solids. This transport mode is investigated by studying the transient settling of heavy particles in turbulent channel flows featuring sediment beds. A numerical method is used to fully resolve turbulence and finite-size particles, which enables the assessment of the complex interplay between flow structures, suspended solids and river sediment.

Description: Soot formation in a laminar and a turbulent diffusion flame as well as the formation of iron oxide nanoparticles in an iron pentacarbonyl doped near stoichiometric H?/O?/Ar-premixed flame have been investigated numerically. The particular modus operandi of the simulated test cases is described. In the case of the turbulent flame the impact of the turbulent closure of the single particle source terms on particle size distribution is shown.

Description: The first in situ measurements of seawater density that referred to a geographical position at sea and time of the year were carried out by Count Luigi Ferdinando Marsili between 1679 and 1680 in the Adriatic Sea, Aegean Sea, Marmara Sea, and the Bosporus. Not only was this the first investigation with documented oceanographic measurements carried out at stations, but themeasurements were described in such an accurateway that the authorswere able to reconstruct the observations in modern units. These first measurements concern the ‘‘specific gravity’’ of seawaters (i.e., the ratio between fluid densities). The data reported in the historical oceanographic treatise Osservazioni intorno al Bosforo Tracio (Marsili) allowed the reconstruction of the seawater density at different geographic locations between 1679 and 1680. Marsili’s experimental methodology included the collection of surface and deep water samples, the analysis of the samples with a hydrostatic ampoule, and the use of a reference water to standardize the measurements.Acomparison of reconstructed densities with present-day values shows an agreement within 10%–20% uncertainty, owing to various aspects of the measurement methodology that are difficult to reconstruct from the documentary evidence. Marsili also measured the current speed and the depth of the current inversion in the Bosporus, which are consistent with the present-day knowledge. The experimental data collected in the Bosporus enabledMarsili to enunciate a theory on the cause of the two-layer flow at the strait, demonstrated by his laboratory experiment and later confirmed by many analytical and numerical studies.

Description: American Meteorological Society.

Description: Published

Description: 845 - 860

Description: 4A. Oceanografia e clima

Description: JCR Journal

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution February 2003

Description: Oceanic observations indicate that abyssal mixing is localized in regions of rough topography. How locally mixed fluid interacts with the ambient fluid is an open question. Laboratory experiments explore the interaction of mechanically induced boundary mixing and an interior body of linearly stratified rotating fluid. A single oscillating bar produces a small region of turbulence along the wall at middepth. Mixed fluid quickly reaches a steady state height set by a turbulent-buoyant balance, independent of rotation. Initially, the bar is exposed on three sides. Mixed fluid intrudes directly into the interior rather than forming a boundary current. The circulation patterns suggest a model of unmixed fluid being laterally entrained into the turbulent zone. In accord with the model, observed outflux is constant, independent of stratification and restricted by rotation. Later the bar is laterally confines between two walls, which form a channel opening into the basin. A small percentage of mixed fluid enters a boundary current, which exits the channel. The bulk forms a cyclonic circulation in front of the bar, which blocks the channel and restricts horizontal entrainment. In the confined case, the volume flux of mixed fluid decays with time.

Description: This work was supported by the Ocean Ventures Fund, the Westcott Fund and the WHOI Education Office. Financial support was also provided by the National Science Foundation through grant OCE-9616949.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 2000

Description: This thesis investigates the physical mechanisms of air-water gas transfer through direct measurements of turbulence at the air--water interface. To enable this study, a new approach to the particle image velocimetry (PIV) technique is developed in order to quantify free-surface flows. Two aspects of this work are innovative. First, the use of a three-dimensional laser light cone and optical filtering of the camera allow for the motion of fluorescent flow tracers at the water surface to be isolated and measured. Validation experiments indicate that this measurement reflects the fluid motion within the upper few hundred microns. A key benefit to this approach is the ability to deal with deforming surfaces, provided the amplitudes are not prohibitively large. This feature was used in this thesis to explore the surface flow induced by mechanically generated waves. Second, a new hybrid PIV image processing algorithm was developed that provides high accuracy velocity estimation with improved computational efficiency. This algorithm combines the concepts of dynamic Fourier-domain cross-correlation with a localized direct multiplicative correlation. In order to explore relationships between free-surface hydrodynamics and air-water gas transfer, an oscillating grid-stirred tank was constructed. By its design, this tank can be managed for chemical cleanliness, offers an unobstructed free surface, and is suited for turbulent mixing and air--water gas-exchange studies. A series of acoustic Doppler velocimeter, PIV, and infrared imaging experiments are presented that characterize the flow for the grid forcing conditions studied. Results indicate that the flows are stationary and reasonably repeatable. In addition, the flows exhibit near-isotropic turbulence and are quasi-homogeneous in horizontal planes. Secondary circulations are revealed and investigated. Finally, PIV measurements of free-surface turbulence are performed with concurrent measurements of gas transfer in the grid tank for a range of turbulent mixing and surface conditions. Surface turbulence, vorticity, and divergence are all affected by the presence of a surface film, with significant effects realized for relatively small surface pressures. Results show that while a relationship between surface turbulence and the gas-transfer velocity is an obvious improvement over that found using an estimate of the bulk flow turbulence, this relationship is dependent on the flow regime. This is revealed through additional surface wave studies. However, the data from both the wave experiments and the grid turbulence experiments can be reconciled by a single relationship between the gas-transfer velocity and the 1/2-power of the surface divergence, which agrees with previous conceptual models. These results (1) further our understanding of interfacial transport processes, (2) demonstrate the important role of surface divergence in air-water gas exchange, and (3) relate, in a physically meaningful way, the interactions between surface renewal, surfactants, and gas transfer.

Description: I was supported as an Office of Naval Research Graduate Fellow. This assistance was the impetus for my pursuit of a doctoral degree and is gratefully acknowledged. Very special thanks also to the WHOI Ocean Ventures Fund Program, Mr. F. Thomas Westcott, and the WHOI Education Department for generous financial support over the past several years.

Description: Author Posting. © American Meteorological Society, 2011. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 24 (2011): 4844–4858, doi:10.1175/2011JCLI4130.1.

Description: The factors that determine the heat transport and overturning circulation in marginal seas subject to wind forcing and heat loss to the atmosphere are explored using a combination of a high-resolution ocean circulation model and a simple conceptual model. The study is motivated by the exchange between the subpolar North Atlantic Ocean and the Nordic Seas, a region that is of central importance to the oceanic thermohaline circulation. It is shown that mesoscale eddies formed in the marginal sea play a major role in determining the mean meridional heat transport and meridional overturning circulation across the sill. The balance between the oceanic eddy heat flux and atmospheric cooling, as characterized by a nondimensional number, is shown to be the primary factor in determining the properties of the exchange. Results from a series of eddy-resolving primitive equation model calculations for the meridional heat transport, overturning circulation, density of convective waters, and density of exported waters compare well with predictions from the conceptual model over a wide range of parameter space. Scaling and model results indicate that wind effects are small and the mean exchange is primarily buoyancy forced. These results imply that one must accurately resolve or parameterize eddy fluxes in order to properly represent the mean exchange between the North Atlantic and the Nordic Seas, and thus between the Nordic Seas and the atmosphere, in climate models.

Description: This study was supported by the National Science Foundation under Grants OCE-0726339 and OCE-0850416.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 1996

Description: This thesis addresses the question of how a highly energetic eddy field could be generated in the interior of the ocean away from the swift boundary currents. The energy radiation due to the temporal growth of non-trapped (radiating) disturbances in such a boundary current is thought to be one of the main sources for the described variability. The problem of stability of an energetic current, such as the Gulf Stream, is formulated. The study then focuses on the ability of the current to support radiating instabilities capable of significant penetration into the far-field and their development with time. The conventional model of the Gulf Stream as a zonal current is extended to allow the jet axis to make an angle to a latitude circle. The linear stability of such a nonzonal flow, uniform in the along-jet direction on a beta-plane, is first studied. The stability computations are performed for piece-wise constant and continuous velocity profiles. New stability properties of nonzonal jets are discussed. In particular, the destabilizing effect of the meridional tilt of the jet axis is demonstrated. The radiating properties of nonzonal currents are found to be very different from those of zonal currents. In particular, purely zonal flows do not support radiating instabilities, whereas flows with a meridional component are capable of radiating long and slowly growing waves. The nonlinear terms are then included in the consideration and the effects of the nonlinear interactions on the radiating properties of the solution are studied in detail. For these purposes, the efficient numerical code for solving equation for the QG potential vorticity with open boundary conditions of Orlanski's type is constructed. The results show that even fast growing linear solutions, which are trapped during the linear stage of developement, can radiate energy in the nonlinear regime if the basic current is nonzonal. The radiation starts as soon as the initial fast exponential growth significantly slows. The initial trapping of those solutions is caused by their fast temporal growth. The new mechanism for radiation is related to the nonzonality of a current.

Description: This work was supported by NSF Grant OCE 9301845.

Description: Author Posting. © American Meteorological Society, 2013. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 26 (2013): 2453–2466, doi:10.1175/JCLI-D-12-00023.1.

Description: The North Atlantic Oscillation (NAO) is one of the most important modes of variability in the global climate system and is characterized by a meridional dipole in the sea level pressure field, with centers of action near Iceland and the Azores. It has a profound influence on the weather, climate, ecosystems, and economies of Europe, Greenland, eastern North America, and North Africa. It has been proposed that around 1980, there was an eastward secular shift in the NAO’s northern center of action that impacted sea ice export through Fram Strait. Independently, it has also been suggested that the location of its southern center of action is tied to the phase of the NAO. Both of these attributes of the NAO have been linked to anthropogenic climate change. Here the authors use both the one-point correlation map technique as well as empirical orthogonal function (EOF) analysis to show that the meridional dipole that is often seen in the sea level pressure field over the North Atlantic is not purely the result of the NAO (as traditionally defined) but rather arises through an interplay among the NAO and two other leading modes of variability in the North Atlantic region: the East Atlantic (EA) and the Scandinavian (SCA) patterns. This interplay has resulted in multidecadal mobility in the two centers of action of the meridional dipole since the late nineteenth century. In particular, an eastward movement of the dipole has occurred during the 1930s to 1950s as well as more recently. This mobility is not seen in the leading EOF of the sea level pressure field in the region.

Description: GWKM was supported by the Natural Sciences and Engineering Research Council of Canada. IAR was supported in part by NE/C003365/1. RSP was supported by Grant OCE-0959381 from the U.S. National Science Foundation.

Description: 2013-10-15

Description: 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 Journal of Physical Oceanography 45 (2015): 2006–2024, doi:10.1175/JPO-D-14-0234.1.

Description: The effects of wind-driven whitecapping on the evolution of the ocean surface boundary layer are examined using an idealized one-dimensional Reynolds-averaged Navier–Stokes numerical model. Whitecapping is parameterized as a flux of turbulent kinetic energy through the sea surface and through an adjustment of the turbulent length scale. Simulations begin with a two-layer configuration and use a wind that ramps to a steady stress. This study finds that the boundary layer begins to thicken sooner in simulations with whitecapping than without because whitecapping introduces energy to the base of the boundary layer sooner than shear production does. Even in the presence of whitecapping, shear production becomes important for several hours, but then inertial oscillations cause shear production and whitecapping to alternate as the dominant energy sources for mixing. Details of these results are sensitive to initial and forcing conditions, particularly to the turbulent length scale imposed by breaking waves and the transfer velocity of energy from waves to turbulence. After 1–2 days of steady wind, the boundary layer in whitecapping simulations has thickened more than the boundary layer in simulations without whitecapping by about 10%–50%, depending on the forcing and initial conditions.

Description: We thank Skidmore College for financial and infrastructure support, and Skidmore and the National Science Foundation for funding travel to meetings where early versions of this work were presented. We also thank the National Science Foundation, Oregon State University, Jonathan Nash, and Joe Jurisa for funding and hosting a workshop on River Plume Mixing in October, 2013, where ideas and context for this paper were developed.

Description: 2016-02-01

Description: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 1555-1566, doi:10.1175/JPO-D-17-0231.1.

Description: A primary challenge in modeling flow over shallow coral reefs is accurately characterizing the bottom drag. Previous studies over continental shelves and sandy beaches suggest surface gravity waves should enhance the drag on the circulation over coral reefs. The influence of surface gravity waves on drag over four platform reefs in the Red Sea is examined using observations from 6-month deployments of current and pressure sensors burst sampling at 1Hz for 4–5min. Depth-average current fluctuations U0 within each burst are dominated by wave orbital velocities uw that account for 80%–90%of the burst variance and have a magnitude of order 10 cm s21, similar to the lower-frequency depth-average current Uavg. Previous studies have shown that the cross-reef bottom stress balances the pressure gradient over these reefs. A bottom stress estimate that neglects the waves (rCdaUavgjUavgj, where r is water density and Cda is a drag coefficient) balances the observed pressure gradient when uw is smaller than Uavg but underestimates the pressure gradient when uw is larger than Uavg (by a factor of 3–5 when uw 5 2Uavg), indicating the neglected waves enhance the bottom stress. In contrast, a bottom stress estimate that includes the waves [rCda(Uavg 1 U0)jUavg 1 U0j)] balances the observed pressure gradient independent of the relative size of uw and Uavg, indicating that this estimate accounts for the wave enhancement of the bottom stress. A parameterization proposed by Wright and Thompson provides a reasonable estimate of the total bottom stress (including the waves) given the burst-averaged current and the wave orbital velocity.

Description: The Red Sea field program was supported by Awards USA 00002 and KSA 00011 made by KAUST. S. Lentz was supported for the analysis by NSF Award OCE-1558343.

Description: 2019-01-13

Description: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Atmospheric and Oceanic Technology 35 (2018): 893-910, doi:10.1175/JTECH-D-17-0102.1.

Description: Rotary sidescan sonars are widely used to image the seabed given their high temporal and spatial resolution. This high resolution is necessary to resolve bedform dynamics and evolution; however, sidescan sonars do not directly measure bathymetry, limiting their utility. When sidescan sonars are mounted close to the seabed, bedforms may create acoustical “shadows” that render previous methods that invert the backscatter intensity to estimate bathymetry and are based on the assumption of a fully insonified seabed ineffective. This is especially true in coastal regions, where bedforms are common features whose large height relative to the water depth may significantly influence the surrounding flow. A method is described that utilizes sonar shadows to estimate bedform height and asymmetry. The method accounts for the periodic structure of bedform fields and the projection of the shadows onto adjacent bedforms. It is validated with bathymetric observations of wave-orbital ripples, with wavelengths ranging from 0.3 to 0.8 m, and tidally reversing megaripples, with wavelengths from 5 to 8 m. In both cases, bathymetric-measuring sonars were deployed in addition to a rotary sidescan sonar to provide a ground truth; however, the bathymetric sonars typically measure different and smaller areas than the rotary sidescan sonar. The shadow-based method and bathymetric-measuring sonar data produce estimates of bedform height that agree by 34.0% ± 27.2% for wave-orbital ripples and 16.6% ± 14.7% for megaripples. Errors for estimates of asymmetry are 1.9% ± 2.1% for wave-orbital ripples and 11.2% ± 9.6% for megaripples.

Description: This project is partially supported by the National Science Foundation through a Graduate Research Fellowship and a Massachusetts Institute of Technology Energy Initiative Fellowship. Additionally, funding used in developing the method was obtained from NSF Grants OCE-1634481 and OCE-1635151. Field work was funded under ONR Grants N00014-06-10329 and N00014-13-1-0364.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution August, 1977

Description: Stimulated by new evidence from both "in situ" oceanic observations and results from numerical modelling, a laboratory study of quasigeostrophic flow and turbulence in a rotating homogeneous fluid has been undertaken. Two dimensional turbulence driven by a uniform distribution of sources and sinks which oscillate in time, can be fairly well reproduced in this context. Inertial time scales are about ten times smaller than Ekman spinup time, and typical Reynolds numbers read 2000. The observations emphasize the spectral tendency of the energy containing eddies. The case of no topography is first discussed. In steadily forced turbulence, it is observed that the energy containing scale is significantly larger than the forcing scale. In the decaying stage the red cascade is observed and rates of interaction are measured. Theoretical arguments for both behaviors are presented; the former concerning the forced turbulence case is believed to be new. The forcing is next applied over various large scale topographies, modelling the geophysical beta effect. The polar beta plane geometry preserves the above spectral characteristics but at the same time introduces anisotropy into the flow pattern. A broad westward mean flow develops in the north and is surrounded by a belt of cyclones lying on its southward side. The calculated second-order Eulerian mean flows induced by steadily and uniformly forced Rossby waves in a long zonal channel, exhibit much of the same momentum distribution in the inertial regime. In contrast, the "sliced cylinder" geometry which possesses no closed geostrophic contours drastically modifies the above picture. Both mean flow production and a large scale tendency for the eddies are inhibited. The geographical distribution of the eddy intensities and scales is now wildly inhomogeneous. The second aspect of this work is a study of the interaction of Rossby waves with mean flows. A zonally traveling, forced wave is generated near the southern boundary of a polar beta plane. Due to energy radiation in the free interior and (or) potential vorticity mixing by the finite amplitude waves, a westward zonal flow develops. The effect of the mean flow upon the forced steady waves is to weaken the anticyclones and intensify the cyclones. Pressure time series reveal a growth of harmonics and general spectral broadening as the waves travel freely inwards, suggesting active nonlinear interactions. An experimental test of Rhines' (1977) potential vorticity mixing theory is also presented at free latitudes. The decay period when the driving is suppressed shows that a net transfer from the waves to the mean flow kinetic energy occurs. Connection with hydrodynamic stability theory is discussed. Interaction of Rossby waves with an externally generated westward mean flow allows one to make a controlled study of the critical layer problem. For small amplitude waves, the mean flow is accelerated in the entire region between the forcing and the critical latitude which acts as a wall for mean wave momentum. In nonlinear runs the steady profile of the westward flow indicates that an accelerating force is acting everywhere, revealing the increasing transmission of wave momentum through the critical layer. At the same time, pressure measurements near the critical point show considerable fine structure developing over a long time scale. The third part deals with steady isolated source-sink flows in the sliced cylinder geometry. The response of the fluid to a meridionally oriented steady dipole extends exclusively westward of the forcing. The viscously balanced solutions are discussed and relevance to oceanic abyssal circulation is emphasized. With strong driving, the combination of a cyclone to the north and an anticyclone to the south is absolutely stable although the reverse configuration is not. A connection with a certain class of free, steady, isolated, inertial solutions developed recently by Stern (1976) is made.

Description: The DGRST . (FRACE) and the Joint Program in Oceanography, Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution offered a fellowship for the first two years. The National Science Foundation under Grant OCE75-2l 674 and the Office of Naval Research under Contract N00014-74-C0262-NR-083-004 supported this study for the final two years.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution January 1977

Description: Millimeter scale fluctuations in refractive index recorded with a freely sinking shadowgraph system are correlated with finestructure profiles of temperature, salinity and density and compared to models of ocean turbulence. Images with vertically aligned periodic structure, called bands, are identified as salt fingers, while others with chaotic structure are turbulent. Images are found on interfaces that are 1-10m thick and have gradients at least several times the mean. From 6 profiles in the Mediterranean Outflow region of the eastern North Atlantic between 1.0 and 1.9 km depth, 398 interfaces have been identified and a significant fraction (about 1/3) of these have detectable images. High contrast images, including bands, are most often found below warm, saline intrusions and within stepped structure where there is a regular sequence of homogeneous mixed layers separated by interfaces. As the interfacial salinity gradient increases in the sense that allows salt finger convection, the fraction of interfaces with images increases. The horizontal spacing of bands (~5 mm) is consistent with calculated salt finger diameter. The calculated and observed length of ocean salt fingers (10-20 cm) is a small fraction of the interface thickness. High levels of small scale variability in the shadowgraphs is reflected in high levels of variance in the finestructure band of the temperature spectra. The temperature gradient spectra have a slope of -1, indicative of turbulence affected by buoyancy forces, and there is a relative peak at a wavelength near the observed salt finger length. The high contrast images are found at interfaces within the enhanced mean salinity gradient below saline intrusions. For very strong salinity gradients there is a solitary interface with intense images, but for weaker mean gradients the convection takes the form of stepped structure. The steps may evolve from the solitary interface as the salinity gradient is run down by salt finger convection. This study identifies parts of the ocean where salt finger convection is prevalent and includes the first comprehensive description of salt fingers in the ocean. Existing models of salt fingers are evaluated in light of ocean observations, and models of ocean turbulence are compared to measurements.

Description: Support from ONR contract N00014-66-C0241 NR. 083-004 is also acknowledged.

Description: Author Posting. © American Meteorological Society, 2009. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 39 (2009): 1077–1096, doi:10.1175/2008JPO4044.1.

Description: Observations of turbulent kinetic energy (TKE) dynamics in the ocean surface boundary layer are presented here and compared with results from previous observational, numerical, and analytic studies. As in previous studies, the dissipation rate of TKE is found to be higher in the wavy ocean surface boundary layer than it would be in a flow past a rigid boundary with similar stress and buoyancy forcing. Estimates of the terms in the turbulent kinetic energy equation indicate that, unlike in a flow past a rigid boundary, the dissipation rates cannot be balanced by local production terms, suggesting that the transport of TKE is important in the ocean surface boundary layer. A simple analytic model containing parameterizations of production, dissipation, and transport reproduces key features of the vertical profile of TKE, including enhancement near the surface. The effective turbulent diffusion coefficient for heat is larger than would be expected in a rigid-boundary boundary layer. This diffusion coefficient is predicted reasonably well by a model that contains the effects of shear production, buoyancy forcing, and transport of TKE (thought to be related to wave breaking). Neglect of buoyancy forcing or wave breaking in the parameterization results in poor predictions of turbulent diffusivity. Langmuir turbulence was detected concurrently with a fraction of the turbulence quantities reported here, but these times did not stand out as having significant differences from observations when Langmuir turbulence was not detected.

Description: The Office of Naval Research funded this work as a part of CBLAST-Low.

Description: Author Posting. © American Meteorological Society, 2008. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 38 (2008): 880–895, doi:10.1175/2007JPO3750.1.

Description: The oceanic response to overflows is explored using a two-layer isopycnal model. Overflows enter the open ocean as dense gravity currents that flow along and down the continental slope. While descending the slope, overflows typically double their volume transport by entraining upper oceanic water. The upper oceanic layer must balance this loss of mass, and the resulting convergent flow produces significant vortex stretching. Overflows thus represent an intense and localized mass and vorticity forcing for the upper ocean. In this study, simulations show that the upper ocean responds to the overflow-induced forcing by establishing topographic β plumes that are aligned more or less along isobaths and that have a transport that is typically a few times larger than that of the overflows. For the topographic β plume driven by the Mediterranean overflow, the occurrence of eddies near Cape St. Vincent, Portugal, allows the topographic β plume to flow across isobaths. The modeled topographic β-plume circulation forms two transatlantic zonal jets that are analogous to the Azores Current and the Azores Countercurrent. In other cases (e.g., the Denmark Strait overflow), the same kind of circulation remains trapped along the western boundary and hence would not be readily detected.

Description: SK’s support during the time of his Ph.D. research in the MIT/WHOI Joint Program was provided by the National Science Foundation through Grant OCE04-24741. JP and JY have also received support from the Climate Process Team on Gravity Current Entrainment, NSF Grant OCE-0611530.

Description: Author Posting. © American Meteorological Society, 2007. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 37 (2007): 1764-1777, doi:10.1175/jpo3098.1.

Description: The vertical structure of the dissipation of turbulence kinetic energy was observed in the nearshore region (3.2-m mean water depth) with a tripod of three acoustic Doppler current meters off a sandy ocean beach. Surface and bottom boundary layer dissipation scaling concepts overlap in this region. No depth-limited wave breaking occurred at the tripod, but wind-induced whitecapping wave breaking did occur. Dissipation is maximum near the surface and minimum at middepth, with a secondary maximum near the bed. The observed dissipation does not follow a surfzone scaling, nor does it follow a “log layer” surface or bottom boundary layer scaling. At the upper two current meters, dissipation follows a modified deep-water breaking-wave scaling. Vertical shear in the mean currents is negligible and shear production magnitude is much less than dissipation, implying that the vertical diffusion of turbulence is important. The increased near-bed secondary dissipation maximum results from a decrease in the turbulent length scale.

Description: Funding was provided by NSF and ONR.

Description: Author Posting. © American Meteorological Society 2006. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 19 (2006): 5366–5387, doi:10.1175/JCLI3892.1.

Description: The oceanic Ekman transport and pumping are among the most important parameters in studying the ocean general circulation and its variability. Upwelling due to the Ekman transport divergence has been identified as a leading mechanism for the seasonal to interannual variability of the upper-ocean heat content in many parts of the World Ocean, especially along coasts and the equator. Meanwhile, the Ekman pumping is the primary mechanism that drives basin-scale circulations in subtropical and subpolar oceans. In those ice-free oceans, the Ekman transport and pumping rate are calculated using the surface wind stress. In the ice-covered Arctic Ocean, the surface momentum flux comes from both air–water and ice–water stresses. The data required to compute these stresses are now available from satellite and buoy observations. But no basin-scale calculation of the Ekman transport in the Arctic Ocean has been done to date. In this study, a suite of satellite and buoy observations of ice motion, ice concentration, surface wind, etc., will be used to calculate the daily Ekman transport over the whole Arctic Ocean from 1978 to 2003 on a 25-km resolution. The seasonal variability and its relationship to the surface forcing fields will be examined. Meanwhile, the contribution of the Ekman transport to the seasonal fluxes of heat and salt to the Arctic Ocean mixed layer will be discussed. It was found that the greatest seasonal variations of Ekman transports of heat and salt occur in the southern Beaufort Sea in the fall and early winter when a strong anticyclonic wind and ice motion are present. The Ekman pumping velocity in the interior Beaufort Sea reaches as high as 10 cm day−1 in November while coastal upwelling is even stronger. The contributions of the Ekman transport to the heat and salt flux in the mixed layer are also considerable in the region.

Description: This study has been supported by NASA Cryospheric Science Program (Grant NNG04GP34G) and by the NSF Office of Polar Program (Grant OPP0424074).

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 1999

Description: Observations of diapycnal mixing rates are examined and related to diapycnal advection for both double-diffusive and turbulent regimes. The role of double-diffusive mixing at the site of the North Atlantic Tracer Release Experiment is considered. The strength of salt-finger mixing is analyzed in terms of the stability parameters for shear and double-diffusive convection, and a nondimensional ratio of the thermal and energy dissipation rates. While the model for turbulence describes most dissipation occurring in high shear, dissipation in low shear is better described by the salt-finger model, and a method for estimating the salt-finger enhancement of the diapycnal haline diffusivity over the thermal diffusivity is proposed. Best agreement between tracer-inferred mixing rates and microstructure based estimates is achieved when the salt-finger enhancement of haline flux is taken into account. The role of turbulence occurring above rough bathymetry in the abyssal Brazil Basin is also considered. The mixing levels along sloping bathymetry exceed the levels observed on ridge crests and canyon floors. Additionally, mixing levels modulate in phase with the spring-neap tidal cycle. A model of the dissipation rate is derived and used to specify the turbulent mixing rate and constrain the diapycnal advection in an inverse model for the steady circulation. The inverse model solution reveals the presence of a secondary circulation with zonal character. These results suggest that mixing in abyssal canyons plays an important role in the mass budget of Antarctic Bottom Water.

Description: This work was supported by contracts N00014-92-1323 and N00014-97-10087 of the Office of Naval Research and grant OCE94-15589 of the National Science Foundation.

Description: Author Posting. © American Meteorological Society, 2011. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 41 (2011): 1874–1893, doi:10.1175/2011JPO4604.1.

Description: A two-dimensional cross-shelf model of the New England continental shelf and slope is used to investigate the mean cross-shelf and vertical circulation at the shelf break and their seasonal variation. The model temperature and salinity fields are nudged toward climatology. Annual and seasonal mean wind stresses are applied on the surface in separate equilibrium simulations. The along-shelf pressure gradient force associated with the along-shelf sea level tilt is tuned to match the modeled and observed depth-averaged along-shelf velocity. Steady-state model solutions show strong seasonal variation in along-shelf and cross-shelf velocity, with the strongest along-shelf jet and interior onshore flow in winter, consistent with observations. Along-shelf sea level tilt associated with the tuned along-shelf pressure gradient increases shoreward because of decreasing water depth. The along-shelf sea level tilt varies seasonally with the wind and is the strongest in winter and weakest in summer. A persistent upwelling is generated at the shelf break with a maximum strength of 2 m day−1 at 50-m depth in winter. The modeled shelfbreak upwelling differs from the traditional view in that most of the upwelled water is from the upper continental slope instead of from the shelf in the form of a detached bottom boundary layer.

Description: WGZ was supported by the Woods Hole Oceanographic Institution postdoctoral scholarship program. GGGandDJMwere supported byONRGrant N-00014- 06-1-0739.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 1992

Description: Oceanic profiles of temperature, salinity, horizontal velocity, rate of dissipation of turbulent kinetic energy (ε) and rate of dissipation of thermal variance (χ) are used to examine the parameterization of turbulent mixing in the ocean due to internal waves. Turbulent mixing is quantified through eddy diffusivity parameterizations of the mass (Kρ; Osborn, 1980) and heat fluxes (Kτ; Osborn and Cox, 1972) in turbulent production/dissipation balances. Turbulence in the ocean is generally held to result from the occurrence of shear instability in regions where the Richardson number is locally supercritical (i.e. Ri ≤ 1/4), permitting the growth of small-scale waves which break and result in turbulent mixing. The occurrence of shear instability results from the local intensification of the shear in the internal wave field. The energy dissipated in such events is provided by the energy flux to higher wavenumber due to nonlinear wave/wave interactions on scales of 10's to 100's of meters. In turn, the strength of the wave/wave interactions depends generally on the energy content of the internal wave field, which can vary considerably over even larger scales due to the presence of topography or background flows. The magnitude of turbulent mixing is linked to internal wave dynamics by equating the turbulent dissipation with the energy flux through the vertical wavenumber spectrum under the priviso that the model spectrum which forms the basis for the analysis is statistically stationary with respect to the nonlinear interactions. Dynamical models (McComas and Muller, 1981; Henyey et al., 1986) indicate that the Garrett and Munk (GM; Munk, 1981) spectrum is stationary. Observations from the far field of a seamount in a region of negligible large-scale flow were examined to address the issue of the buoyancy scaling of ε. These data exhibited large variations in background stratification with depth, but the internal wave characteristics were not substantially differentiable from the GM prescription. The magnitude of ε and its functional dependence upon internal wave energy levels (E) and buoyancy frequency (N) was best described by the dynamical model ofHenyey et al. (1986) (ε ~ E2N2). The Richardson number scaling model of Kunze et al. (1990) produced consistent estimates. A second dynamical model, McComas and Muller (1981), predicted an appropriate (E,N) scaling, but overestimated the observed dissipation rates by a factor of five. Two kinematical dissipation parameterizations (Garmett and Holloway (1984) and Munk (1981)) predicted buoyancy scalings of N3/2 which were inconsistent with the observed scaling. Data from an upper-ocean front, a warm core ring and a region of steep topography were analyzed in order to examine the parameter dependence of E in internal wave fields which exhibited potentially nonstationary characteristics. Evidence was provided which implied the internal wave field in an upper ocean front was interacting with and modified by the background flow. Inhomogeneity and anisotropy of the internal wave field were noted in that data set. The model of Gregg (1989), which in turn was based upon the model of Henyey et al., effectively collapsed the observed diffusivity estimates from the front. The warm core ring profiles were noted to be anisotropic, dominated by near-inertial frequencies and to have a peaked vertical wavenumber shear spectrum. The data from a region of steep topography were noted to have a peaked vertical wavenumber spectrum and were characterized by higher than GM frequency motions. For the latter two data sets, application of a frequency based correction to the Henyey et al. model (Henyey, 1991) reduced more than an order of magnitude scatter in the parameterized estimates of E to less than a factor of four. Of the possible non-equilibrium conditions in the internal wave field, the (E,N) scaled dissipation rates were most sensitive to deviations in wave field frequency content. On the basis of a number of theoretical Richardson number probability distributions (Ri = N2/S2, where S2 is the sum of the squared vertical derivatives of horizontal velocity), the nominal dissipation scaling of the Kunze et al. model was determined to be E2N3. This scaling is altered to the observed ε ~ E2N2 scaling by a statistical dependence between N2 and S2 which reduces the occurrence of supercritical Ri values. This statistical dependence is hypothesized to be an effect of the turbulent momentum and buoyancy fluxes on the internal wave shear and strain profiles caused by shear instability. The statistical dependence between N2 and S2 exhibited a buoyancy scaling which was interpreted as resulting from the decreasing ratio between the time scale of the shear instability mechanism [T- 2π/N] and the adiabatic time scale [T - 2π/(Nf)1/2] of the internal wave field (f is the Coriolis parameter). This phenomenology is interpreted in light of saturated spectral theories which suggest that the magnitude and shape of the vertical wavenumber spectrum is controlled by instability mechanisms at large wavenumber ( ≥ .1 cpm). We argue that saturated spectral theories are valid only in the limit where a separation exists between the two time scales, i.e. for large N, low internal wave frequency content, and small f. These results have immediate implications for oceanic mixing driven by internal wave motions. First, background diffusivities are small: at GM energy levels, Kρ - .03x10-4 m2/s (Kρ = .25ε/N2). Secondly, since Kρ is independent of N at constant E, some process or collection of processes must be responsible for heightened E values in the abyss if internal waves cause the 0(1-10x10-4 m2/s) diffusivities generally inferred from deep ocean hydrographic data. We view internal wave reflection and/or internal wave generation associated with topographic features to be likely candidates.

Description: Author Posting. © American Meteorological Society, 2012. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 42 (2012): 2143–2152, doi:10.1175/JPO-D-12-027.1.

Description: Direct measurements of turbulence levels in the Drake Passage region of the Southern Ocean show a marked enhancement over the Phoenix Ridge. At this site, the Antarctic Circumpolar Current (ACC) is constricted in its flow between the southern tip of South America and the northern tip of the Antarctic Peninsula. Observed turbulent kinetic energy dissipation rates are enhanced in the regions corresponding to the ACC frontal zones where strong flow reaches the bottom. In these areas, turbulent dissipation levels reach 10−8 W kg−1 at abyssal and middepths. The mixing enhancement in the frontal regions is sufficient to elevate the diapycnal turbulent diffusivity acting in the deep water above the axis of the ridge to 1 × 10−4 m2 s−1. This level is an order of magnitude larger than the mixing levels observed upstream in the ACC above smoother bathymetry. Outside of the frontal regions, dissipation rates are O(10−10) W kg−1, comparable to the background levels of turbulence found throughout most mid- and low-latitude regions of the global ocean.

Description: This work was supported by the U.S. National Science Foundation and by the Natural Environment Research Council of the United Kingdom.

Description: 2013-06-01

Description: Author Posting. © American Meteorological Society, 2012. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 42 (2012): 2206–2228, doi:10.1175/JPO-D-11-0191.1.

Description: This study investigates the anisotropic properties of the eddy-induced material transport in the near-surface North Atlantic from two independent datasets, one simulated from the sea surface height altimetry and one derived from real-ocean surface drifters, and systematically examines the interactions between the mean- and eddy-induced material transport in the region. The Lagrangian particle dispersion, which is widely used to characterize the eddy-induced tracer fluxes, is quantified by constructing the “spreading ellipses.” The analysis consistently demonstrates that this dispersion is spatially inhomogeneous and strongly anisotropic. The spreading is larger and more anisotropic in the subtropical than in the subpolar gyre, and the largest ellipses occur in the Gulf Stream vicinity. Even at times longer than half a year, the spreading exhibits significant nondiffusive behavior in some parts of the domain. The eddies in this study are defined as deviations from the long-term time-mean. The contributions from the climatological annual cycle, interannual, and subannual (shorter than one year) variability are investigated, and the latter is shown to have the strongest effect on the anisotropy of particle spreading. The influence of the mean advection on the eddy-induced particle spreading is investigated using the “eddy-following-full-trajectories” technique and is found to be significant. The role of the Ekman advection is, however, secondary. The pronounced anisotropy of particle dispersion is expected to have important implications for distributing oceanic tracers, and for parameterizing eddy-induced tracer transfer in non-eddy-resolving models.

Description: IR was supported by Grant NSF-OCE-0725796. IK would like to acknowledge support by the National Science foundation Grant OCE-0842834.

Description: 2013-06-01

Description: Author Posting. © American Meteorological Society, 2017. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 47 (2017): 633-647, doi:10.1175/JPO-D-16-0089.1.

Description: 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.

Description: DGE was supported by a Natural Environment Research Council studentship award at the University of Southampton. JMT’s contribution was supported by the U.S. National Science Foundation (Grant OCE-1332667). GF’s contribution was supported by the U.S. National Science Foundation through Grant OCE-0961713 and by the U.S. National Oceanic and Atmospheric Administration through Grant NA10OAR4310135. The contributions of JDZ and AJGN were supported by the NERC Grant ‘‘Climate scale analysis of air and water masses’’ (NE/ K012932/1). ACNG gratefully acknowledges support from the Leverhulme Trust, the Royal Society, and the Wolfson Foundation. LY was supported by NASA Ocean Vector Wind Science Team (OVWST) activities under Grant NNA10AO86G.

Description: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 435-453, doi:10.1175/JPO-D-17-0122.1.

Description: Observations of surface waves, currents, and turbulence at the Columbia River mouth are used to investigate the source and vertical structure of turbulence in the surface boundary layer. Turbulent velocity data collected on board freely drifting Surface Wave Instrument Float with Tracking (SWIFT) buoys are corrected for platform motions to estimate turbulent kinetic energy (TKE) and TKE dissipation rates. Both of these quantities are correlated with wave steepness, which has previously been shown to determine wave breaking within the same dataset. Estimates of the turbulent length scale increase linearly with distance from the free surface, and roughness lengths estimated from velocity statistics scale with significant wave height. The vertical decay of turbulence is consistent with a balance between vertical diffusion and dissipation. Below a critical depth, a power-law scaling commonly applied in the literature works well to fit the data. Above this depth, an exponential scaling fits the data well. These results, which are in a surface-following reference frame, are reconciled with results from the literature in a fixed reference frame. A mapping between free-surface and mean-surface reference coordinates suggests 30% of the TKE is dissipated above the mean sea surface.

Description: Funding for this project was provided by the Office of Naval Research as part of the RIVET-II DRI, and for the DARLA group.

Description: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 607-623, doi:10.1175/JPO-D-17-0189.1.

Description: The roles of straining and dissipation in controlling stratification are derived analytically using a vertical salinity variance method. Stratification is produced by converting horizontal variance to vertical variance via straining, that is, differential advection of horizontal salinity gradients, and stratification is destroyed by the dissipation of vertical variance through turbulent mixing. A numerical model is applied to the Changjiang estuary in order to demonstrate the salinity variance balance and how it reveals the factors controlling stratification. The variance analysis reveals that dissipation reaches its maximum during spring tide in the Changjiang estuary, leading to the lowest stratification. Stratification increases from spring tide to neap tide because of the increasing excess of straining over dissipation. Throughout the spring–neap tidal cycle, straining is almost always larger than dissipation, indicating a net excess of production of vertical variance relative to dissipation. This excess is balanced on average by advection, which exports vertical variance out of the estuarine region into the plume. During neap tide, tidal straining shows a general tendency of destratification during the flood tide and restratification during ebb, consistent with the one-dimensional theory of tidal straining. During spring tide, however, positive straining occurs during flood because of the strong baroclinicity induced by the intensified horizontal salinity gradient. These results indicate that the salinity variance method provides a valuable approach for examining the spatial and temporal variability of stratification in estuaries and coastal environments.

Description: X. Li was supported by the China Scholarship Council. W. R. Geyer was supported by NSF Grants OCE 1736539 and OCE 1634480. J. Zhu was supported by the National Natural Science Foundation of China (41476077 and 41676083). H. Wu was supported by the National Natural Science Foundation of China (41576088 and 41776101).

Description: 2018-09-08

Description: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 1831-1848, doi:10.1175/JPO-D-18-0068.1.

Description: We present a simplified theory using reduced-gravity equations for North Atlantic Deep Water (NADW) and its variation driven by high-latitude deep-water formation. The theory approximates layer thickness on the eastern boundary with domain-averaged layer thickness and, in tandem with a mass conservation argument, retains fundamental physics for cross-equatorial flows on interannual and longer forcing time scales. Layer thickness anomalies are driven by a time-dependent northern boundary condition that imposes a southward volume flux representative of a variable source of NADW and damped by diapycnal mixing throughout the basin. Moreover, an outflowing southern boundary condition imposes a southward volume flux that generally differs from the volume flux at the northern boundary, giving rise to temporal storage of NADW within the Atlantic basin. Closed form analytic solutions for the amplitude and phase are provided when the variable source of NADW is sinusoidal. We provide a nondimensional analysis that demonstrates that solution behavior is primarily controlled by two parameters that characterize the meridional extent of the southern basin and the width of the basin relative to the equatorial deformation radius. Similar scaling applied to the time-lagged equations of Johnson and Marshall provides a clear connection to their results. Numerical simulations of reduced-gravity equations agree with analytic predictions in linear, turbulent, and diabatic regimes. The theory introduces a simple analytic framework for studying idealized buoyancy- and wind-driven cross-equatorial flows on interannual and longer time scales.

Description: This research was supported by the National Science Foundation under Grant OCE- 1634468.

Description: 2019-02-15

Description: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 31 (2018): 4847-4863, doi:10.1175/JCLI-D-17-0802.1.

Description: The sensitivity of sea ice to the temperature of inflowing Atlantic water across the Greenland–Scotland Ridge is investigated using an eddy-resolving configuration of the Massachusetts Institute of Technology General Circulation Model with idealized topography. During the last glacial period, when climate on Greenland is known to have been extremely unstable, sea ice is thought to have covered the Nordic seas. The dramatic excursions in climate during this period, seen as large abrupt warming events on Greenland and known as Dansgaard–Oeschger (DO) events, are proposed to have been caused by a rapid retreat of Nordic seas sea ice. Here, we show that a full sea ice cover and Arctic-like stratification can exist in the Nordic seas given a sufficiently cold Atlantic inflow and corresponding low transport of heat across the Greenland–Scotland Ridge. Once sea ice is established, continued sea ice formation and melt efficiently freshens the surface ocean and makes the deeper layers more saline. This creates a strong salinity stratification in the Nordic seas, similar to today’s Arctic Ocean, with a cold fresh surface layer protecting the overlying sea ice from the warm Atlantic water below. There is a nonlinear response in Nordic seas sea ice to Atlantic water temperature with simulated large abrupt changes in sea ice given small changes in inflowing temperature. This suggests that the DO events were more likely to have occurred during periods of reduced warm Atlantic water inflow to the Nordic seas.

Description: The research was supported by the Centre for Climate Dynamics at the Bjerknes Centre for Climate Research. The research leading to these results is part of the ice2ice project funded by the European Research Council under the European Community Seventh Framework Programme (FP7/2007-2013)/ERC Grant Agreement 610055.

Description: Part I: It is to be expected that in the future, measurements of pressure fluctuations beneath the ocean surface will provide basic data for solution of many practical problems on the state of the sea. As a result of the wartime impetus, several types of underwater pressure recorders (of similar instrumentation principles) were developed both in this country and abroad. Essentially, the instrument consists of an underwater unit which electrically transmits pressure impulses near the sea bottom to a clockwork recorder installed on the shore. The underwater pressure unit is adjusted for pressure fluctuations resulting from surface waves within the spectrum band of periods set up by winds acting on the sea surface. The resulting records may be scaled for height and period of the pressure fluctuations over known time intervals. An accessory wave analyzer has been constructed for rapid periodogram analysis of the pressure records. This investigation is concerned with a comparative study of observed sea surface waves and recorded sea bottom pressure fluctuations. It was undertaken for the purpose of evaluating sea surface wave heights from sea bottom pressure recordings in the vicnity of Woods Hole and Bermuda. Part II: The results presented here comprise data on the state of the sea surface in the vicinity of Woods Hole (off Cuttyhunk Island, Massachusetts) July 1946 to May 1947, and off Bermuda, British West Indies, February 1947 to May 1947, obtained from analyses of automatic wave recordings. At both locations measuring elements of the wave recording instruments were located on the sea bottom (approximate depth 75 feet at Cuttyhunk and 120 feet at Bermuda) and electrically connected to shore recorders several miles away. The instrument operational schedules were continuous twenty minute recordings every two hours at Bermuda and every six hours at Cuttyhunk.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution February 1978

Description: The design and operation of a unique flow measuring instrument for bottom boundary layer studies in the marine environment is documented. The effectiveness of the instrument in acquiring data with which models of near bottom flows in the ocean can be tested is demonstrated by the results of a field experiment in Vineyard Sound. The instrument uses four sensors which measure the mean and fluctuating parts of the three components of the velocity vector at four heights above the sea bed. The sensors employ the acoustic travel time difference technique, and are designed to minimize sensor-induced flow disturbances. BASS, an acronym for Benthic Acoustic Stress Sensor, has a resolution of .033 cm/sec per least bit, a range of ±62 cm/sec, noise of .07 cm/sec in 10 sec, and an estimated accuracy of ±.5 cm/sec, referred to an in situ zero point. A complete set of velocity measurements is made every .750 seconds, each measurement being the vector component averaged over 15 cm. The data is internally recorded on digital cassette tape. Eight hours of continuous data can be recorded. BASS was deployed in a tidal flow in Vineyard Sound at a depth of 10 m where a time series of u, v, and w velocities at 26 cm, 46 cm, 96 cm, and 210 cm above the bottom was recorded. The mean velocity was determined by fitting each 6 hour series with a sixth order polynomial and the deviations from the polynomial, the fluctuating velocity components, were correlated to produce Reynolds stress profiles. The stress series shows very few negative stress events while the dominant positive events have an average duration of 5 seconds and exceed 30 dynes/cm2. Zero offset was removed from the mean by assuming a log profile at maximum ebb. Deviations from a log profile developed when the current dropped below 40% of maximum, i.e., when the flow could no longer be considered steady. A break in the Reynolds stress profile at 1 m suggested a larger length scale than the 1 cm bottom roughness was present in the flow. A value of u* was determined by using the quadratic drag law (u* = 1.56 cm/sec), the log profile method (u* = 1.60 cm/sec), and the eddy correlation method (u* = 1.91 cm/sec). Integral length scales of 5 m cross-stream, and 2.5 m vertically were identified by correlation calculations. Two length scales were present in the downstream direction, 5 m within 1 meter of the wall and 8 m further from the wall.

Description: Support from the National Science Foundation is acknowledged.

Description: Author Posting. © American Meteorological Society, 2009. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 39 (2009): 1035-1049, doi:10.1175/2008JPO3920.1.

Description: Seasonal variability of near-inertial horizontal kinetic energy is examined using observations from a series of McLane Moored Profiler moorings located at 39°N, 69°W in the western North Atlantic Ocean in combination with a one-dimensional, depth-integrated kinetic energy model. The time-mean kinetic energy and shear vertical wavenumber spectra of the high-frequency motions at the mooring site are in reasonable agreement with the Garrett–Munk internal wave description. Time series of depth-dependent and depth-integrated near-inertial kinetic energy are calculated from available mooring data after filtering to isolate near-inertial-frequency motions. These data document a pronounced seasonal cycle featuring a wintertime maximum in the depth-integrated near-inertial kinetic energy deriving chiefly from the variability in the upper 500 m of the water column. The seasonal signal in the near-inertial kinetic energy is most prominent for motions with vertical wavelengths greater than 100 m but observable wintertime enhancement is seen down to wavelengths of the order of 10 m. Rotary vertical wavenumber spectra exhibit a dominance of clockwise-with-depth energy, indicative of downward energy propagation and implying a surface energy source. A simple depth-integrated near-inertial kinetic energy model consisting of a wind forcing term and a dissipation term captures the order of magnitude of the observed near-inertial kinetic energy as well as its seasonal cycle.

Description: Funding to initiate the McLane Moored Profiler observations at Line W were provided by grants from the G. Unger Vetlesen Foundation and the Comer Charitable Fund to the Woods Hole Oceanographic Institution’s Ocean and Climate Change Institute. Ongoing moored observations at Line W are supported by the National Science Foundation (NSF Grant OCE-0241354).

Description: Author Posting. © American Meteorological Society, 2007. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 37 (2007): 1496-1511, doi:10.1175/jpo3071.1.

Description: Measurements collected in the York River estuary, Virginia, demonstrate the important impact that tidal and lateral asymmetries in turbulent mixing have on the tidally averaged residual circulation. A reduction in turbulent mixing during the ebb phase of the tide caused by tidal straining of the axial density gradient results in increased vertical velocity shear throughout the water column during the ebb tide. In the absence of significant lateral differences in turbulent mixing, the enhanced ebb-directed transport caused by tidal straining is balanced by a reduction in the net seaward-directed barotropic pressure gradient, resulting in laterally uniform two-layer residual flow. However, the channel–shoal morphology of many drowned river valley estuaries often leads to lateral gradients in turbulent mixing. Tidal straining may then lead to tidal asymmetries in turbulent mixing near the deeper channel while the neighboring shoals remain relatively well mixed. As a result, the largest lateral asymmetries in turbulent mixing occur at the end of the ebb tide when the channel is significantly more stratified than the shoals. The reduced friction at the end of ebb delays the onset of the flood tide, increasing the duration of ebb in the channel. Conversely, over the shoal regions where stratification is more inhibited by tidal mixing, there is greater friction and the transition from ebb to flood occurs more rapidly. The resulting residual circulation is seaward over the channel and landward over the shoal. The shoal–channel segregation of this barotropically induced estuarine residual flow is opposite to that typically associated with baroclinic estuarine circulation over channel–shoal bathymetry.

Description: Support for this research was provided by the National Science Foundation Division of Ocean Sciences Grant OCE- 9984941.

Description: Author Posting. © American Meteorological Society, 2007. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography. 37 (2007): 2363-2386, doi:10.1175/jpo3118.1.

Description: Intrinsic low-frequency variability is studied in the idealized, quasigeostrophic, midlatitude, wind-driven ocean gyres operating at large Reynolds number. A robust decadal variability mode driven by the transient mesoscale eddies is found and analyzed. The variability is a turbulent phenomenon, which is driven by the competition between the eddy rectification process and the potential vorticity anomalies induced by changes of the intergyre transport

Description: Funding for Pavel Berloff was provided by NSF Grants OCE-0091836 and OCE- 0344094, by the U.K. Royal Society Fellowship, and by the Newton Trust Award, A. M. Hogg was supported by an Australian Research Council Postdoctoral Fellowship (DP0449851) during this work, and William K. Dewar was supported by NSF Grants OCE-0424227 and OCE-0550139.

Description: Author Posting. © American Meteorological Society, 2011. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 41 (2011): 166-185, doi:10.1175/2010JPO4470.1.

Description: Field observations of turbulent kinetic energy (TKE), dissipation rate ε, and turbulent length scale demonstrate the impact of both density stratification and nonlocal turbulent production on turbulent momentum flux. The data were collected in a highly stratified salt wedge estuary using the Mobile Array for Sensing Turbulence (MAST). Estimates of the dominant length scale of turbulent motions obtained from the vertical velocity spectra provide field confirmation of the theoretical limitation imposed by either the distance to the boundary or the Ozmidov scale, whichever is smaller. Under boundary-limited conditions, anisotropy generally increases with increasing shear and decreased distance to the boundary. Under Ozmidov-limited conditions, anisotropy increases rapidly when the gradient Richardson number exceeds 0.25. Both boundary-limited and Ozmidov-limited conditions demonstrate significant deviations from a local production–dissipation balance that are largely consistent with simple scaling relationships for the vertical divergence in TKE flux. Both the impact of stratification and deviation from equilibrium turbulence observed in the data are largely consistent with commonly used turbulence closure models that employ “nonequilibrium” stability functions. The data compare most favorably with the nonequilibrium version of the L. H. Kantha and C. A. Clayson stability functions. Not only is this approach more consistent with the observed critical gradient Richardson number of 0.25, but it also accounts for the large deviations from equilibrium turbulence in a manner consistent with the observations.

Description: The funding for this research was obtained from ONR Grant N00014-06-1-0292 and NSF Grants and OCE-08-25226 and OCE-08-24871.

Description: Author Posting. © American Meteorological Society, 2012. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 42 (2012): 855–868, doi:10.1175/JPO-D-10-05010.1.

Description: Data from the Hudson River estuary demonstrate that the tidal variations in vertical salinity stratification are not consistent with the patterns associated with along-channel tidal straining. These observations result from three additional processes not accounted for in the traditional tidal straining model: 1) along-channel and 2) lateral advection of horizontal gradients in the vertical salinity gradient and 3) tidal asymmetries in the strength of vertical mixing. As a result, cross-sectionally averaged values of the vertical salinity gradient are shown to increase during the flood tide and decrease during the ebb. Only over a limited portion of the cross section does the observed stratification increase during the ebb and decrease during the flood. These observations highlight the three-dimensional nature of estuarine flows and demonstrate that lateral circulation provides an alternate mechanism that allows for the exchange of materials between surface and bottom waters, even when direct turbulent mixing through the pycnocline is prohibited by strong stratification.

Description: The funding for this research was obtained from NSF Grant OCE-08-25226.

Description: 2012-11-01

Description: Author Posting. © American Meteorological Society, 2013. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 43 (2013): 259–282, doi:10.1175/JPO-D-11-0194.1.

Description: This study reports on observations of turbulent dissipation and internal wave-scale flow properties in a standing meander of the Antarctic Circumpolar Current (ACC) north of the Kerguelen Plateau. The authors characterize the intensity and spatial distribution of the observed turbulent dissipation and the derived turbulent mixing, and consider underpinning mechanisms in the context of the internal wave field and the processes governing the waves’ generation and evolution. The turbulent dissipation rate and the derived diapycnal diffusivity are highly variable with systematic depth dependence. The dissipation rate is generally enhanced in the upper 1000–1500 m of the water column, and both the dissipation rate and diapycnal diffusivity are enhanced in some places near the seafloor, commonly in regions of rough topography and in the vicinity of strong bottom flows associated with the ACC jets. Turbulent dissipation is high in regions where internal wave energy is high, consistent with the idea that interior dissipation is related to a breaking internal wave field. Elevated turbulence occurs in association with downward-propagating near-inertial waves within 1–2 km of the surface, as well as with upward-propagating, relatively high-frequency waves within 1–2 km of the seafloor. While an interpretation of these near-bottom waves as lee waves generated by ACC jets flowing over small-scale topographic roughness is supported by the qualitative match between the spatial patterns in predicted lee wave radiation and observed near-bottom dissipation, the observed dissipation is found to be only a small percentage of the energy flux predicted by theory. The mismatch suggests an alternative fate to local dissipation for a significant fraction of the radiated energy.

Description: SW acknowledges the support of the Grantham Institute for Climate Change, Imperial College London. ACNG acknowledges the support of a NERC Advanced Research Fellowship (Grant NE/C517633/1). KLP acknowledges support from Woods Hole Oceanographic Institution bridge support funds.

Description: 2013-08-01

Description: Author Posting. © American Meteorological Society, 2013. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 43 (2013): 222–230, doi:10.1175/JPO-D-12-099.1.

Description: Observations with fine horizontal resolution are used to identify the horizontal scales of variability over the Middle Atlantic Bight (MAB) shelf break and continental rise. Spray gliders collected observations along two alongshelf transects over the continental rise in March–April 2006 and along 16 cross-shelf transects over the shelf break and continental rise during July–October 2007. Horizontal resolution varied from 1 km or finer over the shelf to 6 km in deep water. These observations allow horizontal thermohaline variability offshore of the MAB shelf break to be examined for the first time. Structure functions of temperature and salinity, the mean square difference between observations separated by specified distances, reveal the horizontal spatial scales in the region. Exponential (e-folding) scales of temperature and salinity increase from 8–13 km near the shelf break to about 30 km over the continental rise. Just offshore of the shelf break, alongshelf structure functions exhibit periodicity with a 40–50-km wavelength that matches the wavelength of shelfbreak frontal meanders. Farther offshore, alongshelf structure functions suggest a dominant wavelength of 175–250 km, but these scales are only marginally resolved by the available observations. Examination of structure functions of along-isopycnal salinity (i.e., spice) suggests that interleaving of shelf and slope water masses contributes most of the horizontal variability near the MAB shelf break, but heaving of isopycnals is the primary source of horizontal variability over the continental rise.

Description: Glider observations in March–April 2006 were supported by the National Science Foundation through Grant OCE-0220769. Glider observations in July–October 2007 were supported by a grant from Raytheon. RET was supported by the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Cooperative Institute for the North Atlantic Region. GGG was supported by the National Science Foundation under Grant OCE-1129125.

Description: 2013-07-01

Description: Author Posting. © American Meteorological Society, 2014. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 44 (2014): 427–444, doi:10.1175/JPO-D-13-070.1.

Description: Between 25 September 2007 and 28 September 2009, a heavily instrumented mooring was deployed in the Labrador Sea, offshore of the location where warm-core, anticyclonic Irminger rings are formed. The 2-year time series offers insight into the vertical and horizontal structure of newly formed Irminger rings and their heat and salt transport into the interior basin. In 2 years, 12 Irminger rings passed by the mooring. Of these, 11 had distinct properties, while 1 anticyclone likely passed the mooring twice. Eddy radii (11–35 km) were estimated using the dynamic height signal of the anticyclones (8–18 cm) together with the observed velocities. The anticyclones show a seasonal cycle in core properties when observed (1.9°C in temperature and 0.07 in salinity at middepth) that has not been described before. The temperature and salinity are highest in fall and lowest in spring. Cold, fresh caps, suggested to be an important source of freshwater, were seen in spring but were almost nonexistent in fall. The heat and freshwater contributions by the Irminger rings show a large spread (from 12 to 108 MJ m−2 and from −0.5 to −4.7 cm, respectively) for two reasons. First, the large range of radii leads to large differences in transported volume. Second, the seasonal cycle leads to changes in heat and salt content per unit volume. This implies that estimates of heat and freshwater transport by eddies should take the distribution of eddy properties into account in order to accurately assess their contribution to the restratification.

Description: This work was supported by the U.S. National Science Foundation and the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Devonshire Foundation.

Description: 2014-08-01

Description: Author Posting. © American Meteorological Society, 2016. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 46 (2016): 327–348, doi:10.1175/JPO-D-15-0112.1.

Description: Potential vorticity structure in two segments of the North Atlantic’s western boundary current is examined using concurrent, high-resolution measurements of hydrography and velocity from gliders. Spray gliders occupied 40 transects across the Loop Current in the Gulf of Mexico and 11 transects across the Gulf Stream downstream of Cape Hatteras. Cross-stream distributions of the Ertel potential vorticity and its components are calculated for each transect under the assumptions that all flow is in the direction of measured vertically averaged currents and that the flow is geostrophic. Mean cross-stream distributions of hydrographic properties, potential vorticity, and alongstream velocity are calculated for both the Loop Current and the detached Gulf Stream in both depth and density coordinates. Differences between these mean transects highlight the downstream changes in western boundary current structure. As the current increases its transport downstream, upper-layer potential vorticity is generally reduced because of the combined effects of increased anticyclonic relative vorticity, reduced stratification, and increased cross-stream density gradients. The only exception is within the 20-km-wide cyclonic flank of the Gulf Stream, where intense cyclonic relative vorticity results in more positive potential vorticity than in the Loop Current. Cross-stream gradients of mean potential vorticity satisfy necessary conditions for both barotropic and baroclinic instability within the western boundary current. Instances of very low or negative potential vorticity, which predispose the flow to various overturning instabilities, are observed in individual transects across both the Loop Current and the Gulf Stream.

Description: Glider operations in the Gulf Stream were supported by the National Science Foundation under Grant OCE-0220769. Glider operations in the Gulf of Mexico were supported by BP. R.E.T. was supported by the Penzance Endowed Fund in Support of Assistant Scientists and the Independent Research and Development Program at WHOI.

Description: 2016-07-01

Description: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 905-923, doi:10.1175/JPO-D-17-0133.1.

Description: Observations of turbulent kinetic energy, dissipation, and turbulent stress were collected in the middle reaches of Chesapeake Bay and were used to assess second-moment closure predictions of turbulence generated beneath breaking waves. Dissipation scaling indicates that the turbulent flow structure observed during a 10-day wind event was dominated by a three-layer response that consisted of 1) a wave transport layer, 2) a surface log layer, and 3) a tidal, bottom boundary layer limited by stable stratification. Below the wave transport layer, turbulent mixing was limited by stable stratification. Within the wave transport layer, where dissipation was balanced by a divergence in the vertical turbulent kinetic energy flux, the eddy viscosity was significantly underestimated by second-moment turbulence closure models, suggesting that breaking waves homogenized the mixed surface layer to a greater extent than the simple model of TKE diffusing away from a source at the surface. While the turbulent transport of TKE occurred largely downgradient, the intermittent downward sweeps of momentum generated by breaking waves occurred largely independent of the mean shear. The underprediction of stress in the wave transport layer by second-moment closures was likely due to the inability of the eddy viscosity model to capture the nonlocal turbulent transport of the momentum flux beneath breaking waves. Finally, the authors hypothesize that large-scale coherent turbulent eddies played a significant role in transporting momentum generated near the surface to depth.

Description: This work was supported by National Science Foundation Grants OCE-1061609 and OCE-1339032.

Description: 2018-10-19

Description: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 1815-1830, doi:10.1175/JPO-D-17-0275.1.

Description: Recent progress in direct numerical simulations (DNSs) of stratified turbulent flows has led to increasing attention to the validity of the constancy of the dissipation flux coefficient Γ in the Osborn’s eddy diffusivity model. Motivated by lack of observational estimates of Γ, particularly under weakly stratified deep-ocean conditions, this study estimates Γ using deep microstructure profiles collected in various regions of the North Pacific and Southern Oceans. It is shown that Γ is not constant but varies significantly with the Ozmidov/Thorpe scale ratio ROT in a fashion similar to that obtained by previous DNS studies. Efficient mixing events with Γ ~ O(1) and ROT ~ O(0.1) tend to be frequently observed in the deep ocean (i.e., weak stratification), while moderate mixing events with Γ ~ O(0.1) and ROT ~ O(1) tend to be observed in the upper ocean (i.e., strong stratification). The observed negative relationship between Γ and ROT is consistent with a simple scaling that can be derived from classical turbulence theories. In contrast, the observed results exhibit no definite relationships between Γ and the buoyancy Reynolds number Reb, although Reb has long been thought to be another key parameter that controls Γ.

Description: This study was supported by MEXT KAKENHI Grant JP15H05824 and JSPS KAKENHI Grant JP15H02131.

Description: 2019-02-15

Description: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018):1941-1950, doi:10.1175/JPO-D-17-0194.1.

Description: Subglacial discharges have been observed to generate buoyant plumes along the ice face of Greenland tidewater glaciers. These plumes have been traditionally modeled using classical plume theory, and their characteristic parameters (e.g., velocity) are employed in the widely used three-equation melt parameterization. However, the applicability of plume theory for three-dimensional turbulent wall plumes is questionable because of the complex near-wall plume dynamics. In this study, corrections to the classical plume theory are introduced to account for the presence of a wall. In particular, the drag and entrainment coefficients are quantified for a three-dimensional turbulent wall plume using data from direct numerical simulations. The drag coefficient is found to be an order of magnitude larger than that for a boundary layer flow over a flat plate at a similar Reynolds number. This result suggests a significant increase in the melting estimates by the current parameterization. However, the volume flux in a wall plume is found to be one-half that of a conical plume that has 2 times the buoyancy flux. This finding suggests that the total entrainment (per unit area) of ambient water is the same and that the plume scalar characteristics (i.e., temperature and salinity) can be predicted reasonably well using classical plume theory.

Description: This work was supported by the Linné FLOW Centre at KTH and the Academy of Finland Center of Excellence Programme Grant 307331 (author Ezhova) and by VR Swedish Research Council GrantVR2014-5001 (author Brandt). Support to author Cenedese was given by NSF Project OCE-1434041.

Description: 2019-02-27

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 1996

Description: Nonlinear quasigeostrophic flows in two layers over a topographic slope are considered. The evolution depends on the size of two parameters which indicate the degree of nonlinearity at depth. The first measures the importance of relative vorticity advection and the second of stretching vorticity. Two types of isolated vortex are used to examine the parameter dependence. An initially barotropic vortex remains barotropic only when the first parameter is large, otherwise topographic waves dominate at depth. An Initially surface-trapped vortex larger than deformation scale is baroclinically unstable when the second is large, but is stabilized by the slope otherwise. Both parameters are also relevant to cascading geostrophic turbulence. If the stretching parameter is large, a "barotropic cascade" occurs at the deformation radius (Rhines, 1977) and the cascade "arrests" when the relative vorticity parameter is order unity. If small, layer coupling is hindered and the cascade is arrested at the deformation scale, with the flow dominated by isotropic surface vortices. In both cases, the distinction between vortices and waves is transparent when viewing potential vorticity. It is more difficult to identify waves and vortices from the streamfunction fields, because the waves are present in both layers.

Description: Funding for this research was provided by Office of Naval Research Coastal Science Code, grants N00014-92-J-1643 and N00014-92-J-1528.

Description: Author Posting. © American Meteorological Society, 2009. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 39 (2009): 3162-3175, doi:10.1175/2009JPO4239.1.

Description: This study analyzes anisotropic properties of the material transport by eddies and eddy-driven zonal jets in a general circulation model of the North Atlantic through the analysis of Lagrangian particle trajectories. Spreading rates—defined here as half the rate of change in the particle dispersion—in the zonal direction systematically exceed the meridional rates by an order of magnitude. Area-averaged values for the upper-ocean zonal and meridional spreading rates are approximately 8100 and 1400 m2 s−1, respectively, and in the deep ocean they are 2400 and 200 m2 s−1. The results demonstrate that this anisotropy is mainly due to the action of the transient eddies and not to the shear dispersion associated with the time-mean jets. This property is consistent with the fact that eddies in this study have zonally elongated shapes. With the exception of the upper-ocean subpolar gyre, eddies also cause the superdiffusive zonal spreading, significant variations in the spreading rate in the vertical and meridional directions, and the difference between the westward and eastward spreading.

Description: Funding for IK was provided by NSF Grants OCE 0346178, 0749722, and 0842834. Funding for PB was provided by NSF Grants OCE 0344094 and OCE 0725796 and by the research grant from the Newton Trust of the University of Cambridge. For JP the acknowledgement is to NSF OCE-0451086.

Description: Author Posting. © American Meteorological Society, 2009. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 39 (2009): 1361-1379, doi:10.1175/2008JPO4096.1.

Description: Multiple zonal jets are observed in satellite data–based estimates of oceanic velocities, float measurements, and high-resolution numerical simulations of the ocean circulation. This study makes a step toward understanding the dynamics of these jets in the real ocean by analyzing the vertical structure and dynamical balances within multiple zonal jets simulated in an eddy-resolving primitive equation model of the North Atlantic. In particular, the authors focus on the role of eddy flux convergences (“eddy forcing”) in supporting the buoyancy and relative/potential vorticity (PV) anomalies associated with the jets. The results suggest a central role of baroclinic eddies in the barotropic and baroclinic dynamics of the jets, and significant differences in the effects of eddy forcing between the subtropical and subpolar gyres. Additionally, diabatic potential vorticity sources and sinks, associated with vertical diffusion, are shown to play an important role in supporting the potential vorticity anomalies. The resulting potential vorticity profile does not resemble a “PV staircase”—a distinct meridional structure observed in some idealized studies of geostrophic turbulence.

Description: Funding for IK was provided by NSF Grants OCE 0346178 and 0749722. Funding for PB was provided by NSF Grants OCE 0344094 and OCE 0725796 and by the research grant from the Newton Trust of the University of Cambridge. For JP the acknowledgement is to NSF OCE-0451086.

Description: Author Posting. © American Meteorological Society, 2011. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 41 (2011): 2307–2327, doi:10.1175/JPO-D-10-05004.1.

Description: Results from a high-resolution (~2 km) numerical simulation of the Irminger Basin during summer 2003 are presented. The focus is on the East Greenland Spill Jet, a recently discovered component of the circulation in the basin. The simulation compares well with observations of surface fields, the Denmark Strait overflow (DSO), and the hydrographic structure of typical sections in the basin. The model reveals new aspects of the circulation on scales of O(0.1–10) days and O(1–100) km. The model Spill Jet results from the cascade of dense waters over the East Greenland shelf. Spilling can occur in various locations southwest of the strait, and it is present throughout the simulation but exhibits large variations on periods of O(0.1–10) days. The Spill Jet sometimes cannot be distinguished in the velocity field from surface eddies or from the DSO. The vorticity structure of the jet confirms its unstable nature with peak relative and tilting vorticity terms reaching twice the planetary vorticity term. The average model Spill Jet transport is 4.9 ±1.7 Sv (1 Sv ≡ 106 m3 s−1) equatorward, about 2½ times larger than has been previously reported from a single ship transect in August 2001. Kinematic analysis of the model results suggests two different types of spilling events. In the first case (type I), a local perturbation results in dense waters descending over the shelf break into the Irminger Basin. In the second case (type II), surface cyclones associated with DSO deep domes initiate the spilling process. During summer 2003, more than half of the largest Spill Jet transport values are of type II.

Description: The research is supported by the National Science Foundation Grants OCE-0726393 and OCI-0904640 (MGM and TWNH) and OCE-0726640 (RSP).

Description: 2012-06-01

Description: Author Posting. © American Meteorological Society, 2011. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 41 (2011): 2223–2241, doi:10.1175/2011JPO4344.1.

Description: Results are presented from an observational study of stratified, turbulent flow in the bottom boundary layer on the outer southeast Florida shelf. Measurements of momentum and heat fluxes were made using an array of acoustic Doppler velocimeters and fast-response temperature sensors in the bottom 3 m over a rough reef slope. Direct estimates of flux Richardson number Rf confirm previous laboratory, numerical, and observational work, which find mixing efficiency not to be a constant but rather to vary with Frt, Reb, and Rig. These results depart from previous observations in that the highest levels of mixing efficiency occur for Frt 〈 1, suggesting that efficient mixing can also happen in regions of buoyancy-controlled turbulence. Generally, the authors find that turbulence in the reef bottom boundary layer is highly variable in time and modified by near-bed flow, shear, and stratification driven by shoaling internal waves.

Description: Funding was provided by grants from the National Oceanic and Atmospheric Administration’s National Undersea Research Program, National Science Foundation Grants OCE-0622967 and OCE- 0824972 to SGM, and the Singapore Stanford Program. Kristen Davis was supported by a National Defense Science and Engineering Graduate Fellowship and an ARCS Foundation Fellowship.

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 1993

Description: A turbulent mixing experiment was conducted to observe the dynamics and the energetics of layer formation along with the region of layer formation in the Reynolds number (Re) and the overall Richardson number (Rio) space. A salt stratified fluid was mixed uniformly throughout its depth with a vertical rod that moved horizontally at a constant speed. The evolution of density was measured with a conductivity probe. As the instability theory of Phillips (1972) and Posmentier (1977) shows, an initially uniform density profile turns into a series of steps when Rio is larger than a critical value Ric, which forms a stability boundary. For fixed Re, as Rio decreases to Ric, the steps get weaker; the density difference across the interface and the difference of density gradient between layers and interfaces become small. Ric increases as Re increases with a functional relation log Ric ≈ Re/900. The steps evolve over time, with small steps forming first, and larger steps appearing later through merging and decay of the interfaces. After some time the interior seems to reach an equilibrium state and the evolution of the interior steps stops. The length scale of the equilibrium step, ls, is a linear function of U /Ni, where U is the speed of the rod and Ni is the buoyancy frequency of the initial profile. The functional relationship is ls = 2.6U / Ni + l.Ocm. For Rio 〈 Ric, the mixing efficiency, Rf, monotonically decreases to the end of a run. However, for Rio 〉 Ric, the evolution of Rf is closely related to the evolution of the density field. Rf changes rapidly during the initiation of the steps. For Rio » Ric, R1 increases initially, while for Rio ≥ Ric, Rf ecreases initially. When the interior reaches an equilibrium state, Rf becomes uniform. Posmentier (1977) theorized that when steps reach an equilibrium state, a density flux is independent of the density gradient. The present experiments show a uniform density flux in the layered interior irrespective of the density structure, and this strongly supports the theory of Posmentier. The density flux generated in the bottom boundary mixed layer goes through the interior all the way to the top boundary mixed layer without changing the interior density structure. Thus, turbulence can transport scalar properties further than the characteristic length scale of active eddies without changing a density structure. When the fluid becomes two mixed layers, the relation between Rf and Ril was found for Ril 〉 1. Here, Ril is the local Richardson number based on the thickness of the interface. R, does decrease as Ril increases, which is the most crucial assumption of the instability theory.

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 1995

Description: This thesis describes the steps necessary to develop an acoustic vorticity meter for the atmosphere. The analysis is based on Benthic Acoustic Stress Sensor (BASS) technology that is currently used for similar acoustic measurements in the ocean. Compared to sonic anemometer measurements, the BASS measurements of velocity are not only made in a different fluid but in a different way. Due to these differences, the physical make up of BASS needed to be altered, and the validity of the measurement technique had to be explored. The alterations to the BASS hardware occurred for several reasons. Because attenuation of sound is much higher in air than in water for the same frequencies, it was necessary to change the transducers. The generally faster and unidirectional mean flows that are present in the air encourage open measurement volumes which the BASS vorticity meters do not have. The difference in group speed of sound is different for water and air, and this forced a change to the timing and burst generation board of the BASS vorticity meter. The measurement technique used by the BASS instrumentation is validated by the error analysis in the text. Because the BASS instrumentation actually provides a time difference, the equation used by the BASS instrumentation to compute velocity was assumed throughout the error analysis. The error analysis shows that the combination of BASS measurement techniques with a temperature sensor will provide errors that are less than 2% of the velocity. The types of measurements that an atmospheric vorticity meter would provide to a researcher are described in the text to show the meter's potential. If deployed on a buoy, a vorticity meter could measure shearing of the wind close to the surface of the waves. If deployed at heights much greater than its path lengths, an atmospheric vorticity meter could provide three-dimensional vorticity measurements which would provide a unique measurement of a fundamental characteristic of turbulent flows.

Description: 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 Journal of Physical Oceanography 45 (2015): 2381–2406, doi:10.1175/JPO-D-14-0086.1.

Description: While near-inertial waves are known to be generated by atmospheric storms, recent observations in the Kuroshio Front find intense near-inertial internal-wave shear along sloping isopycnals, even during calm weather. Recent literature suggests that spontaneous generation of near-inertial waves by frontal instabilities could represent a major sink for the subinertial quasigeostrophic circulation. An unforced three-dimensional 1-km-resolution model, initialized with the observed cross-Kuroshio structure, is used to explore this mechanism. After several weeks, the model exhibits growth of 10–100-km-scale frontal meanders, accompanied by O(10) mW m−2 spontaneous generation of near-inertial waves associated with readjustment of submesoscale fronts forced out of balance by mesoscale confluent flows. These waves have properties resembling those in the observations. However, they are reabsorbed into the model Kuroshio Front with no more than 15% dissipating or radiating away. Thus, spontaneous generation of near-inertial waves represents a redistribution of quasigeostrophic energy rather than a significant sink.

Description: “The Study of Kuroshio Ecosystem Dynamics for Sustainable Fisheries (SKED)” supported by MEXT, MIT-Hayashi Seed Fund, ONR (Awards N000140910196 and N000141210101), NSF (Award OCE 0928617, 0928138) for support.

Description: 2016-03-01

Description: 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 Journal of Physical Oceanography 45 (2015): 2497–2521, doi:10.1175/JPO-D-14-0128.1.

Description: Oceanic density overturns are commonly used to parameterize the dissipation rate of turbulent kinetic energy. This method assumes a linear scaling between the Thorpe length scale LT and the Ozmidov length scale LO. Historic evidence supporting LT ~ LO has been shown for relatively weak shear-driven turbulence of the thermocline; however, little support for the method exists in regions of turbulence driven by the convective collapse of topographically influenced overturns that are large by open-ocean standards. This study presents a direct comparison of LT and LO, using vertical profiles of temperature and microstructure shear collected in the Luzon Strait—a site characterized by topographically influenced overturns up to O(100) m in scale. The comparison is also done for open-ocean sites in the Brazil basin and North Atlantic where overturns are generally smaller and due to different processes. A key result is that LT/LO increases with overturn size in a fashion similar to that observed in numerical studies of Kelvin–Helmholtz (K–H) instabilities for all sites but is most clear in data from the Luzon Strait. Resultant bias in parameterized dissipation is mitigated by ensemble averaging; however, a positive bias appears when instantaneous observations are depth and time integrated. For a series of profiles taken during a spring tidal period in the Luzon Strait, the integrated value is nearly an order of magnitude larger than that based on the microstructure observations. Physical arguments supporting LT ~ LO are revisited, and conceptual regimes explaining the relationship between LT/LO and a nondimensional overturn size are proposed. In a companion paper, Scotti obtains similar conclusions from energetics arguments and simulations.

Description: B.D.M. and S.K.V. gratefully acknowledge the support of the Office of Naval Research under Grants N00014-12-1-0279, N00014-12-1-0282, and N00014-12-1-0938 (Program Manager: Dr. Terri Paluszkiewicz). S.K.V. also acknowledges support of the National Science Foundation under Grant OCE-1151838. L.S.L. acknowledges support for BBTRE by the National Science Foundation by Contract OCE94-15589 and NATRE and IWISE by the Office of Naval Research by Contracts N00014-92-1323 and N00014-10-10315. J.N.M. was supported through Grant 1256620 from the National Science Foundation and the Office of Naval Research (IWISE Project).

Description: 2016-04-01

Description: 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 Journal of Physical Oceanography 45 (2015): 2621–2639, doi:10.1175/JPO-D-14-0239.1.

Description: Measurements made as part of a large-scale experiment to examine wind-driven circulation and mixing in Chesapeake Bay demonstrate that circulations consistent with Langmuir circulation play an important role in surface boundary layer dynamics. Under conditions when the turbulent Langmuir number Lat is low (〈0.5), the surface mixed layer is characterized by 1) elevated vertical turbulent kinetic energy; 2) decreased anisotropy; 3) negative vertical velocity skewness indicative of strong/narrow downwelling and weak/broad upwelling; and 4) strong negative correlations between low-frequency vertical velocity and the velocity in the direction of wave propagation. These characteristics appear to be primarily the result of the vortex force associated with the surface wave field, but convection driven by a destabilizing heat flux is observed and appears to contribute significantly to the observed negative vertical velocity skewness. Conditions that favor convection usually also have strong Langmuir forcing, and these two processes probably both contribute to the surface mixed layer turbulence. Conditions in which traditional stress-driven turbulence is important are limited in this dataset. Unlike other shallow coastal systems where full water column Langmuir circulation has been observed, the salinity stratification in Chesapeake Bay is nearly always strong enough to prevent full-depth circulation from developing.

Description: The funding for this research was provided by the National Science Foundation Grants OCE-1339032 and OCE-1338518.

Description: 2016-04-01

Description: Author Posting. © American Meteorological Society, 2016. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 46 (2016): 3155-3163, doi:10.1175/JPO-D-16-0123.1.

Description: Idealized laboratory experiments have been conducted in a two-layer stratified fluid to investigate the leading-order dynamics that control submarine melting and meltwater export near a vertical ice–ocean interface as a function of subglacial discharge. In summer, the discharge of surface runoff at the base of a glacier (subglacial discharge) generates strong buoyant plumes that rise along the glacier front entraining ambient water along the way. The entrainment enhances the heat transport toward the glacier front and hence the submarine melt rate increases with the subglacial discharge rate. In the laboratory, the effect of subglacial discharge is simulated by introducing freshwater at freezing temperature from a point source at the base of an ice block representing the glacier. The circulation pattern observed both with and without subglacial discharge resembles those observed in previous observational and numerical studies. Buoyant plumes rise vertically until they find either their neutrally buoyant level or the free surface. Hence, the meltwater can deposit within the interior of the water column and not entirely at the free surface, as confirmed by field observations. The heat budget in the tank, calculated following a new framework, gives estimates of submarine melt rate that increase with the subglacial discharge and are in agreement with the directly measured submarine melting. This laboratory study provides the first direct measurements of submarine melt rates for different subglacial discharges, and the results are consistent with the predictions of previous theoretical and numerical studies.

Description: Support to C. C. was given by the NSF project OCE- 1130008 and OCE-1434041. M. G. received support from the ‘‘Gori’’ Fellowship.

Description: 2017-04-07

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution June 2017

Description: This thesis investigates the formation and subsequent motion of oceanic lithospheric plates through geophysical and petrological methods. Ocean crust and lithosphere forms at mid-ocean ridges as the underlying asthenosphere rises, melts, and flows away from the ridge axis. In Chapters 2 and 3, I present the results from partial melting experiments of mantle peridotite that were conducted in order to examine the mantle melting point, or solidus, beneath a mid-ocean ridge. Chapter 2 determines the peridotite solidus at a single pressure of 1.5 GPa and concludes that the oceanic mantle potential temperature must be ~60ºC hotter than current estimates. Chapter 3 goes further to provide a more accurate parameterization of the anhydrous mantle solidus from experiments over a range of pressures. This chapter concludes that the range of potential temperatures of the mantle beneath mid-ocean ridges and plumes is smaller than currently estimated. Once formed, the oceanic plate moves atop the underlying asthenosphere away from the ridge axis. Chapter 4 uses seafloor magnetotelluric data to investigate the mechanism responsible for plate motion at the lithosphere-asthenosphere boundary. The resulting two dimensional conductivity model shows a simple layered structure. By applying petrological constraints, I conclude that the upper asthenosphere does not contain substantial melt, which suggests that either a thermal or hydration mechanism supports plate motion. Oceanic plate motion has dramatically changed the surface of the Earth over time, and evidence for ancient plate motion is obvious from detailed studies of the longer lived continental lithosphere. In Chapter 5, I investigate past plate motion by inverting magnetotelluric data collected over eastern Zambia. The conductivity model probes the Zambian lithosphere and reveals an ancient subduction zone previously suspected from surface studies. This chapter elucidates the complex lithospheric structure of eastern Zambia and the geometry of the tectonic elements in the region, which collided as a result of past oceanic plate motion. Combined, the chapters of this thesis provide critical constraints on ocean plate dynamics.

Description: Funding for this research was provided by the National Science Foundation Division of Earth Sciences (EAR) grant number 1010432, Division of Ocean Sciences (OCE) grant numbers 1459649 and 0928663, WHOI Deep Ocean Exploration Institute, and the WHOI Academic Programs Office.

Description: Author Posting. © American Meteorological Society, 2017. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 47 (2017): 1061-1075, doi:10.1175/JPO-D-16-0248.1.

Description: A major challenge in modeling the circulation over coral reefs is uncertainty in the drag coefficient because existing estimates span two orders of magnitude. Current and pressure measurements from five coral reefs are used to estimate drag coefficients based on depth-average flow, assuming a balance between the cross-reef pressure gradient and the bottom stress. At two sites wind stress is a significant term in the cross-reef momentum balance and is included in estimating the drag coefficient. For the five coral reef sites and a previous laboratory study, estimated drag coefficients increase as the water depth decreases consistent with open channel flow theory. For example, for a typical coral reef hydrodynamic roughness of 5 cm, observational estimates, and the theory indicate that the drag coefficient decreases from 0.4 in 20 cm of water to 0.005 in 10 m of water. Synthesis of results from the new field observations with estimates from previous field and laboratory studies indicate that coral reef drag coefficients range from 0.2 to 0.005 and hydrodynamic roughnesses generally range from 2 to 8 cm. While coral reef drag coefficients depend on factors such as physical roughness and surface waves, a substantial fraction of the scatter in estimates of coral reef drag coefficients is due to variations in water depth.

Description: The Red Sea field program was supported by Awards USA 00002 and KSA 00011 made by KAUST to S. Lentz and J. Churchill. The Palau field program was funded by NSF Award OCE-1220529.

Description: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 1367-1373, doi:10.1175/JPO-D-17-0185.1.

Description: An earlier study indicates that the side melting of icebergs subject to vertically homogeneous horizontal velocities is controlled by two distinct regimes, which depend on the melt plume behavior and produce a nonlinear dependence of side melt rate on velocity. Here, we extend this study to consider ice blocks melting in a two-layer vertically sheared flow in a laboratory setting. It is found that the use of the vertically averaged flow speed in current melt parameterizations gives an underestimate of the submarine side melt rate, in part because of the nonlinearity of the dependence of the side melt rate on flow speed but also because vertical shear in the horizontal velocity profile fundamentally changes the flow splitting around the ice block and consequently the velocity felt by the ice surface. An observational record of 90 icebergs in a Greenland fjord suggests that this effect could produce an average underestimate of iceberg side melt rates of 21%.

Description: A. F. was supported by NA14OAR4320106 from the National Oceanic and Atmospheric Administration, U.S. Department of Commerce. C. C. was supported by NSF OCE-1658079 and F. S. was supported by NSF OCE-1657601 and NSF PLR-1743693.

Description: 2018-12-12

Description: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 31 (2018): 8059-8079, doi:10.1175/JCLI-D-17-0769.1.

Description: We use the method of least squares with Lagrange multipliers to fit an ocean general circulation model to the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) estimate of near sea surface temperature (NSST) at the Last Glacial Maximum (LGM; circa 23–19 thousand years ago). Compared to a modern simulation, the resulting global, last-glacial ocean state estimate, which fits the MARGO data within uncertainties in a free-running coupled ocean–sea ice simulation, has global-mean NSSTs that are 2°C lower and greater sea ice extent in all seasons in both the Northern and Southern Hemispheres. Increased brine rejection by sea ice formation in the Southern Ocean contributes to a stronger abyssal stratification set principally by salinity, qualitatively consistent with pore fluid measurements. The upper cell of the glacial Atlantic overturning circulation is deeper and stronger. Dye release experiments show similar distributions of Southern Ocean source waters in the glacial and modern western Atlantic, suggesting that LGM NSST data do not require a major reorganization of abyssal water masses. Outstanding challenges in reconstructing LGM ocean conditions include reducing effects from model biases and finding computationally efficient ways to incorporate abyssal tracers in global circulation inversions. Progress will be aided by the development of coupled ocean–atmosphere–ice inverse models, by improving high-latitude model processes that connect the upper and abyssal oceans, and by the collection of additional paleoclimate observations.

Description: DEA was supported by a NSF Graduate Research Fellowship and NSF Grant OCE-1060735. OM acknowledges support from the NSF. GF was supported by NASA Award 1553749 and Simons Foundation Award 549931.

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 33(4), (2020): 1535-1545, doi:10.1175/JCLI-D-19-0547.1.

Description: In a transient warming scenario, the North Atlantic is influenced by a complex pattern of surface buoyancy flux changes that ultimately weaken the Atlantic meridional overturning circulation (AMOC). Here we study the AMOC response in the CMIP5 experiment, using the near-geostrophic balance of the AMOC on interannual time scales to identify the role of temperature and salinity changes in altering the circulation. The thermal wind relationship is used to quantify changes in the zonal density gradients that control the strength of the flow. At 40°N, where the overturning cell is at its strongest, weakening of the AMOC is largely driven by warming between 1000- and 2000-m depth along the western margin. Despite significant subpolar surface freshening, salinity changes are small in the deep branch of the circulation. This is likely due to the influence of anomalously salty water in the subpolar intermediate layers, which is carried northward from the subtropics in the upper limb of the AMOC. In the upper 1000 m at 40°N, salty anomalies due to increased evaporation largely cancel the buoyancy increase due to warming. Therefore, in CMIP5, temperature dynamics are responsible for AMOC weakening, while freshwater forcing instead acts to strengthen the circulation in the net. These results indicate that past modeling studies of AMOC weakening, which rely on freshwater hosing in the subpolar gyre, may not be directly applicable to a more complex warming scenario.

Description: We acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups (listed in Table 1 of this paper) for producing and making available their model output. We also thank John Marshall for helpful discussions on the driving mechanisms of the AMOC, and three anonymous reviewers whose comments greatly improved the manuscript. This work was supported by NASA Headquarters under the NASA Earth and Space Science Fellowship Program Award 80NSSC17K0372, and by National Science Foundation Award OCE-1433132.

Description: 2020-07-20

Description: Author Posting. © American Meteorological Society, 2010. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 40 (2010): 2381-2400, doi:10.1175/2010JPO4403.1.

Description: Langmuir circulation (LC) is a turbulent upper-ocean process driven by wind and surface waves that contributes significantly to the transport of momentum, heat, and mass in the oceanic surface layer. The authors have previously performed a direct comparison of large-eddy simulations and observations of the upper-ocean response to a wind event with rapid mixed layer deepening. The evolution of simulated crosswind velocity variance and spatial scales, as well as mixed layer deepening, was only consistent with observations if LC effects are included in the model. Based on an analysis of these validated simulations, in this study the fundamental differences in mixing between purely shear-driven turbulence and turbulence with LC are identified. In the former case, turbulent kinetic energy (TKE) production due to shear instabilities is largest near the surface, gradually decreasing to zero near the base of the mixed layer. This stands in contrast to the LC case in which at middepth range TKE production can be dominated by Stokes drift shear. Furthermore, the Eulerian mean vertical shear peaks near the base of the mixed layer so that TKE production by mean shear flow is elevated there. LC transports horizontal momentum efficiently downward leading to an along-wind velocity jet below LC downwelling regions at the base of the mixed layer. Locally enhanced vertical shear instabilities as a result of this jet efficiently erode the thermocline. In turn, enhanced breaking internal waves inject cold deep water into the mixed layer, where LC currents transport temperature perturbation advectively. Thus, LC and locally generated shear instabilities work intimately together to facilitate strongly the mixed layer deepening process.

Description: This research was supported by the Office of Naval Research through Grants N00014-09-M-0112 (TK) and N00014-06-1-0178 (AP, JT). Author TK also received support from a Woods Hole Oceanographic Institution Cooperative Institute for Climate and Ocean Research Postdoctoral Scholarship.

Description: Author Posting. © American Meteorological Society, 2013. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 43 (2013): 1841–1861, doi:10.1175/JPO-D-12-0231.1.

Description: In this idealized numerical modeling study, the composition of residual sediment fluxes in energetic (e.g., weakly or periodically stratified) tidal estuaries is investigated by means of one-dimensional water column models, with some focus on the sediment availability. Scaling of the underlying dynamic equations shows dependence of the results on the Simpson number (relative strength of horizontal density gradient) and the Rouse number (relative settling velocity) as well as impacts of the Unsteadiness number (relative tidal frequency). Here, the parameter space given by the Simpson and Rouse numbers is mainly investigated. A simple analytical model based on the assumption of stationarity shows that for small Simpson and Rouse numbers sediment flux is down estuary and vice versa for large Simpson and Rouse numbers. A fully dynamic water column model coupled to a second-moment turbulence closure model allows to decompose the sediment flux profiles into contributions from the transport flux (product of subtidal velocity and sediment concentration profiles) and the fluctuation flux profiles (tidal covariance between current velocity and sediment concentration). Three different types of bottom sediment pools are distinguished to vary the sediment availability, by defining a time scale for complete sediment erosion. For short erosion times scales, the transport sediment flux may dominate, but for larger erosion time scales the fluctuation sediment flux largely dominates the tidal sediment flux. When quarter-diurnal components are added to the tidal forcing, up-estuary sediment fluxes are strongly increased for stronger and shorter flood tides and vice versa. The theoretical results are compared to field observations in a tidally energetic inlet.

Description: Project funding was provided by the German Research Foundation (DFG) in the framework of the Project ECOWS (Role of Estuarine Circulation for Transport of Suspended Particulate Matter in the Wadden Sea, BU 1199/11) and by the German Federal Ministry of Research and Education in the framework of the Project PACE [The future of the Wadden Sea sediment fluxes: still keeping pace with sea level rise? (FKZ 03F0634A)].

Description: 2014-03-01

Description: Author Posting. © American Meteorological Society, 2014. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 44 (2014): 1466–1492, doi:10.1175/JPO-D-12-0154.1.

Description: Simultaneous full-depth microstructure measurements of turbulence and finestructure measurements of velocity and density are analyzed to investigate the relationship between turbulence and the internal wave field in the Antarctic Circumpolar Current. These data reveal a systematic near-bottom overprediction of the turbulent kinetic energy dissipation rate by finescale parameterization methods in select locations. Sites of near-bottom overprediction are typically characterized by large near-bottom flow speeds and elevated topographic roughness. Further, lower-than-average shear-to-strain ratios indicative of a less near-inertial wave field, rotary spectra suggesting a predominance of upward internal wave energy propagation, and enhanced narrowband variance at vertical wavelengths on the order of 100 m are found at these locations. Finally, finescale overprediction is typically associated with elevated Froude numbers based on the near-bottom shear of the background flow, and a background flow with a systematic backing tendency. Agreement of microstructure- and finestructure-based estimates within the expected uncertainty of the parameterization away from these special sites, the reproducibility of the overprediction signal across various parameterization implementations, and an absence of indications of atypical instrument noise at sites of parameterization overprediction, all suggest that physics not encapsulated by the parameterization play a role in the fate of bottom-generated waves at these locations. Several plausible underpinning mechanisms based on the limited available evidence are discussed that offer guidance for future studies.

Description: The SOFine project is funded by the United Kingdom’s Natural Environmental Research Council (NERC) (Grant NE/G001510/1). SW acknowledges the support of anARCDiscovery Early CareerResearchAward (Grant DE120102927), as well as the Grantham Institute for Climate Change, Imperial College London, and the ARC Centre of Excellence for Climate System Science (Grant CE110001028). ACNG acknowledges the support of a NERC Advanced Research Fellowship (Grant NE/C517633/1).KLP acknowledges support fromWoods Hole Oceanographic Institution bridge support funds.

Description: 2014-11-01

Description: Author Posting. © American Meteorological Society, 2014. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 44 (2014): 2593–2616, doi:10.1175/JPO-D-13-0120.1.

Description: The first direct estimate of the rate at which geostrophic turbulence mixes tracers across the Antarctic Circumpolar Current is presented. The estimate is computed from the spreading of a tracer released upstream of Drake Passage as part of the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES). The meridional eddy diffusivity, a measure of the rate at which the area of the tracer spreads along an isopycnal across the Antarctic Circumpolar Current, is 710 ± 260 m2 s−1 at 1500-m depth. The estimate is based on an extrapolation of the tracer-based diffusivity using output from numerical tracers released in a one-twentieth of a degree model simulation of the circulation and turbulence in the Drake Passage region. The model is shown to reproduce the observed spreading rate of the DIMES tracer and suggests that the meridional eddy diffusivity is weak in the upper kilometer of the water column with values below 500 m2 s−1 and peaks at the steering level, near 2 km, where the eddy phase speed is equal to the mean flow speed. These vertical variations are not captured by ocean models presently used for climate studies, but they significantly affect the ventilation of different water masses.

Description: NSF support through Awards OCE-1233832, OCE-1232962, and OCE-1048926 is gratefully acknowledged.

Description: 2015-04-01

Description: Author Posting. © American Meteorological Society, 2016. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 46 (2016): 1309-1321, doi:10.1175/JPO-D-15-0068.1.

Description: 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.

Description: This work was supported by grants from the U.S. National Science Foundation.

Description: 2016-10-05

Description: Author Posting. © American Meteorological Society, 2016. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 46 (2016): 1769-1783, doi:10.1175/JPO-D-15-0193.1.

Description: High-resolution observations of velocity, salinity, and turbulence quantities were collected in a salt wedge estuary to quantify the efficiency of stratified mixing in a high-energy environment. During the ebb tide, a midwater column layer of strong shear and stratification developed, exhibiting near-critical gradient Richardson numbers and turbulent kinetic energy (TKE) dissipation rates greater than 10−4 m2 s−3, based on inertial subrange spectra. Collocated estimates of scalar variance dissipation from microconductivity sensors were used to estimate buoyancy flux and the flux Richardson number Rif. The majority of the samples were outside the boundary layer, based on the ratio of Ozmidov and boundary length scales, and had a mean Rif = 0.23 ± 0.01 (dissipation flux coefficient Γ = 0.30 ± 0.02) and a median gradient Richardson number Rig = 0.25. The boundary-influenced subset of the data had decreased efficiency, with Rif = 0.17 ± 0.02 (Γ = 0.20 ± 0.03) and median Rig = 0.16. The relationship between Rif and Rig was consistent with a turbulent Prandtl number of 1. Acoustic backscatter imagery revealed coherent braids in the mixing layer during the early ebb and a transition to more homogeneous turbulence in the midebb. A temporal trend in efficiency was also visible, with higher efficiency in the early ebb and lower efficiency in the late ebb when the bottom boundary layer had greater influence on the flow. These findings show that mixing efficiency of turbulence in a continuously forced, energetic, free shear layer can be significantly greater than the broadly cited upper bound from Osborn of 0.15–0.17.

Description: Holleman was supported by the Devonshire Scholars program. The field study and the coauthors’ contributions were supported by NSF Grant OCE 0926427.

Description: 2016-11-24

Description: © The Author(s), 2016. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Journal of Atmospheric and Oceanic Technology 33 (2016): 873-890, doi:10.1175/JTECH-D-15-0109.1.

Description: Direct covariance flux (DCF) measurements taken from floating platforms are contaminated by wave-induced platform motions that need to be removed before computation of the turbulent fluxes. Several correction algorithms have been developed and successfully applied in earlier studies from research vessels and, most recently, by the use of moored buoys. The validation of those correction algorithms has so far been limited to short-duration comparisons against other floating platforms. Although these comparisons show in general a good agreement, there is still a lack of a rigorous validation of the method, required to understand the strengths and weaknesses of the existing motion-correction algorithms. This paper attempts to provide such a validation by a comparison of flux estimates from two DCF systems, one mounted on a moored buoy and one on the Air–Sea Interaction Tower (ASIT) at the Martha’s Vineyard Coastal Observatory, Massachusetts. The ASIT was specifically designed to minimize flow distortion over a wide range of wind directions from the open ocean for flux measurements. The flow measurements from the buoy system are corrected for wave-induced platform motions before computation of the turbulent heat and momentum fluxes. Flux estimates and cospectra of the corrected buoy data are found to be in very good agreement with those obtained from the ASIT. The comparison is also used to optimize the filter constants used in the motion-correction algorithm. The quantitative agreement between the buoy data and the ASIT demonstrates that the DCF method is applicable for turbulence measurements from small moving platforms, such as buoys.

Description: This work was funded by the National Science Foundation Grant OCE04-24536 as part of the CLIVAR Mode Water Dynamic Experiment (CLIMODE).

Description: © The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Kwon, Y., Seo, H., Ummenhofer, C. C., & Joyce, T. M. Impact of multidecadal variability in Atlantic SST on winter atmospheric blocking. Journal of Climate, 33(3), (2020): 867-892, doi: 10.1175/JCLI-D-19-0324.1.

Description: Recent studies have suggested that coherent multidecadal variability exists between North Atlantic atmospheric blocking frequency and the Atlantic multidecadal variability (AMV). However, the role of AMV in modulating blocking variability on multidecadal times scales is not fully understood. This study examines this issue primarily using the NOAA Twentieth Century Reanalysis for 1901–2010. The second mode of the empirical orthogonal function for winter (December–March) atmospheric blocking variability in the North Atlantic exhibits oppositely signed anomalies of blocking frequency over Greenland and the Azores. Furthermore, its principal component time series shows a dominant multidecadal variability lagging AMV by several years. Composite analyses show that this lag is due to the slow evolution of the AMV sea surface temperature (SST) anomalies, which is likely driven by the ocean circulation. Following the warm phase of AMV, the warm SST anomalies emerge in the western subpolar gyre over 3–7 years. The ocean–atmosphere interaction over these 3–7-yr periods is characterized by the damping of the warm SST anomalies by the surface heat flux anomalies, which in turn reduce the overall meridional gradient of the air temperature and thus weaken the meridional transient eddy heat flux in the lower troposphere. The anomalous transient eddy forcing then shifts the eddy-driven jet equatorward, resulting in enhanced Rossby wave breaking and blocking on the northern flank of the jet over Greenland. The opposite is true with the AMV cold phases but with much shorter lags, as the evolution of SST anomalies differs in the warm and cold phases.

Description: We gratefully acknowledge support from the NSF Climate and Large-scale Dynamics Program (AGS-1355339) to Y-OK, HS, CCU, and TMJ, the NASA Physical Oceanography Program (NNX13AM59G) to Y-OK, HS, and TMJ, NOAA CPO Climate Variability and Predictability Program (NA13OAR4310139) and DOE CESD Regional and Global Model Analysis Program (DE-SC0019492) to Y-OK, and NSF Physical Oceanography Program (OCE-1419235) to HS. We are very grateful to the three anonymous reviewers and editor Dr. Mingfang Ting, for their thorough and insightful suggestions. The NOAA 20CR dataset was downloaded from the NOAA Earth System Research Laboratory Physical Science Division webpage (https://www.esrl.noaa.gov/psd/data/20thC_Rean/). Support for the 20CR Project version 2c dataset is provided by the U.S. Department of Energy, Office of Science Biological and Environmental Research (BER), and by the National Oceanic and Atmospheric Administration Climate Program Office. The HadISST dataset was downloaded from the U.K. Met Office Hadley Centre webpage (https://www.metoffice.gov.uk/hadobs/hadisst/). The ERA-20C dataset was downloaded from the ECMWF webpage (https://apps.ecmwf.int/datasets/data/era20c-daily/). The ERSST5 dataset was provided by the NOAA Earth System Research Laboratory Physical Science Division (https://www.esrl.noaa.gov/psd/data/gridded/data.noaa.ersst.v5.html).

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(2), (2020): 415-437, doi:10.1175/JPO-D-19-0019.1.

Description: Results are presented from two dye release experiments conducted in the seasonal thermocline of the Sargasso Sea, one in a region of low horizontal strain rate (~10−6 s−1), the second in a region of intermediate horizontal strain rate (~10−5 s−1). Both experiments lasted ~6 days, covering spatial scales of 1–10 and 1–50 km for the low and intermediate strain rate regimes, respectively. Diapycnal diffusivities estimated from the two experiments were κz = (2–5) × 10−6 m2 s−1, while isopycnal diffusivities were κH = (0.2–3) m2 s−1, with the range in κH being less a reflection of site-to-site variability, and more due to uncertainties in the background strain rate acting on the patch combined with uncertain time dependence. The Site I (low strain) experiment exhibited minimal stretching, elongating to approximately 10 km over 6 days while maintaining a width of ~5 km, and with a notable vertical tilt in the meridional direction. By contrast, the Site II (intermediate strain) experiment exhibited significant stretching, elongating to more than 50 km in length and advecting more than 150 km while still maintaining a width of order 3–5 km. Early surveys from both experiments showed patchy distributions indicative of small-scale stirring at scales of order a few hundred meters. Later surveys show relatively smooth, coherent distributions with only occasional patchiness, suggestive of a diffusive rather than stirring process at the scales of the now larger patches. Together the two experiments provide important clues as to the rates and underlying processes driving diapycnal and isopycnal mixing at these scales.

Description: Results are presented from two dye release experiments conducted in the seasonal thermocline of the Sargasso Sea, one in a region of low horizontal strain rate (~10−6 s−1), the second in a region of intermediate horizontal strain rate (~10−5 s−1). Both experiments lasted ~6 days, covering spatial scales of 1–10 and 1–50 km for the low and intermediate strain rate regimes, respectively. Diapycnal diffusivities estimated from the two experiments were κz = (2–5) × 10−6 m2 s−1, while isopycnal diffusivities were κH = (0.2–3) m2 s−1, with the range in κH being less a reflection of site-to-site variability, and more due to uncertainties in the background strain rate acting on the patch combined with uncertain time dependence. The Site I (low strain) experiment exhibited minimal stretching, elongating to approximately 10 km over 6 days while maintaining a width of ~5 km, and with a notable vertical tilt in the meridional direction. By contrast, the Site II (intermediate strain) experiment exhibited significant stretching, elongating to more than 50 km in length and advecting more than 150 km while still maintaining a width of order 3–5 km. Early surveys from both experiments showed patchy distributions indicative of small-scale stirring at scales of order a few hundred meters. Later surveys show relatively smooth, coherent distributions with only occasional patchiness, suggestive of a diffusive rather than stirring process at the scales of the now larger patches. Together the two experiments provide important clues as to the rates and underlying processes driving diapycnal and isopycnal mixing at these scales.

Description: 2020-08-06

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of the Atmospheric and Oceanic Technology 37(5), (2020): 789-806, doi:10.1175/JTECH-D-18-0244.1.

Description: Realistic ocean state prediction and its validation rely on the availability of high quality in situ observations. To detect data errors, adequate quality check procedures must be designed. This paper presents procedures that take advantage of the ever-growing observation databases that provide climatological knowledge of the ocean variability in the neighborhood of an observation location. Local validity intervals are used to estimate binarily whether the observed values are considered as good or erroneous. Whereas a classical approach estimates validity bounds from first- and second-order moments of the climatological parameter distribution, that is, mean and variance, this work proposes to infer them directly from minimum and maximum observed values. Such an approach avoids any assumption of the parameter distribution such as unimodality, symmetry around the mean, peakedness, or homogeneous distribution tail height relative to distribution peak. To reach adequate statistical robustness, an extensive manual quality control of the reference dataset is critical. Once the data have been quality checked, the local minima and maxima reference fields are derived and the method is compared with the classical mean/variance-based approach. Performance is assessed in terms of statistics of good and bad detections. It is shown that the present size of the reference datasets allows the parameter estimates to reach a satisfactory robustness level to always make the method more efficient than the classical one. As expected, insufficient robustness persists in areas with an especially low number of samples and high variability.

Description: This study has been conducted using EU Copernicus Marine Service Information and was supported by the European Union within the EU Copernicus Marine Service In Situ phase-I and phase-II contracts led by Ifremer. The publication was also supported by SOERE CTDO2 in France. The Argo data were collected and made freely available by the International Argo Program and the national programs that contribute to it (see http://www.argo.ucsd.edu, http://argo.jcommops.org). The Argo Program is part of the Global Ocean Observing System (http://doi.org/10.17882/42182). The marine mammal data were collected and made freely available by the International MEOP Consortium and the national programs that contribute to it (see http://www.meop.net; https://doi.org/10.17882/45461). Aleix Gelabert and Dídac Costa were the skippers of the OPOO, sponsored by the Intergovernmental Oceanographic Commission (UNESCO) and Pharmaton. The BWR is a periodic oceanic race organized by the Fundació Navegació Oceànica de Barcelona (FNOB). Reviewer D. Briand provided some useful comments on the final version of the draft paper before submission.

Description: 2020-11-04

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(5), (2020): 1245-1263, doi:10.1175/JPO-D-19-0213.1.

Description: We use laboratory experiments and theoretical modeling to investigate the surface expression of a subglacial discharge plume, as occurs at many fjords around Greenland. The experiments consider a fountain that is released vertically into a homogeneous fluid, adjacent either to a vertical or a sloping wall, that then spreads horizontally at the free surface before sinking back to the bottom. We present a model that separates the fountain into two separate regions: a vertical fountain and a horizontal, negatively buoyant jet. The model is compared to laboratory experiments that are conducted over a range of volume fluxes, density differences, and ambient fluid depths. It is shown that the nondimensionalized length, width, and aspect ratio of the surface expression are dependent on the Froude number, calculated at the start of the negatively buoyant jet. The model is applied to observations of the surface expression from a Greenland subglacial discharge plume. In the case where the discharge plume reaches the surface with negative buoyancy the model can be used to estimate the discharge properties at the base of the glacier.

Description: We gratefully acknowledge technical assistance from Anders Jensen and thank anonymous reviewers for improving the clarity of the manuscript. CM thanks the Weston Howard Jr. Scholarship for funding. Support to CC was given by NSF project OCE-1434041 and OCE-1658079.

Description: 2020-10-27

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(10), (2022): 2325–2341, https://doi.org/10.1175/jpo-d-21-0015.1.

Description: The ocean surface boundary layer is a gateway of energy transfer into the ocean. Wind-driven shear and meteorologically forced convection inject turbulent kinetic energy into the surface boundary layer, mixing the upper ocean and transforming its density structure. In the absence of direct observations or the capability to resolve subgrid-scale 3D turbulence in operational ocean models, the oceanography community relies on surface boundary layer similarity scalings (BLS) of shear and convective turbulence to represent this mixing. Despite their importance, near-surface mixing processes (and ubiquitous BLS representations of these processes) have been undersampled in high-energy forcing regimes such as the Southern Ocean. With the maturing of autonomous sampling platforms, there is now an opportunity to collect high-resolution spatial and temporal measurements in the full range of forcing conditions. Here, we characterize near-surface turbulence under strong wind forcing using the first long-duration glider microstructure survey of the Southern Ocean. We leverage these data to show that the measured turbulence is significantly higher than standard shear-convective BLS in the shallower parts of the surface boundary layer and lower than standard shear-convective BLS in the deeper parts of the surface boundary layer; the latter of which is not easily explained by present wave-effect literature. Consistent with the CBLAST (Coupled Boundary Layers and Air Sea Transfer) low winds experiment, this bias has the largest magnitude and spread in the shallowest 10% of the actively mixing layer under low-wind and breaking wave conditions, when relatively low levels of turbulent kinetic energy (TKE) in surface regime are easily biased by wave events.

Description: This paper is VIMS Contribution 4103. Computational resources were provided by the VIMS Ocean-Atmosphere and Climate Change Research Fund. AUSSOM was supported by the OCE Division of the National Science Foundation (1558639).

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science in Physical Oceanography at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution February 2021.

Description: Estimating turbulence in the marine-atmospheric boundary layer is critical to many industrial, commercial and scientific fields, but of particular importance to the wind energy industry. Contributing to both the efficiency of energy extraction and the life-cycle cost of the turbine itself, turbulence in the atmospheric boundary layer is estimated within the wind energy industry as Turbulence Intensity (TI) and more recently by Turbulent Kinetic Energy (TKE). Traditional in-situ methods to measure turbulence are extremely difficult to deploy in the marine environment, resulting in a recent movement to and dependence on remote sensing methods. One type of remote sensing instrument, Doppler lidars, have shown to reliably estimate the wind speed and atmospheric turbulence while being cost effective and easily deployable, and hence are being increasingly utilized as a standard for wind energy assessments. In this thesis, the ability of lidars to measure turbulence up to a height of 200 m above mean sea level in the marine-atmospheric boundary layer was tested using a 7-month data set spanning winter to early summer. Lidar-based TI and TKE were estimated by three methods using observations from a highly validated lidar system and compared under both convective and stable atmospheric stability conditions. Convective periods were found to have higher turbulence at all the heights compared to stable conditions, while mean wind speed and shear were higher during stable conditions. The study period was characterized by generally low turbulent conditions with high turbulence events occurring at timescales of a few days. Mean vertical profiles of TKE were non-uniformly distributed in height during low turbulent conditions. During highly turbulent events, TKE increased more strongly with height. The definition of TI– following the industry or meteorology conventions – had no real effect on the results, and differences between cup or sonic anemometers and lidar TI values were small except at low wind speeds. All the three lidar-based TKE methods tested corresponded closely to independent estimates, and differences between the methods were small relative to the temporal variability of TKE observed at the offshore site.

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 51(1), (2021): 19-35, https://doi.org/10.1175/JPO-D-19-0233.1.

Description: In the Beaufort Sea in September of 2015, concurrent mooring and microstructure observations were used to assess dissipation rates in the vicinity of 72°35′N, 145°1′W. Microstructure measurements from a free-falling profiler survey showed very low [O(10−10) W kg−1] turbulent kinetic energy dissipation rates ε. A finescale parameterization based on both shear and strain measurements was applied to estimate the ratio of shear to strain Rω and ε at the mooring location, and a strain-based parameterization was applied to the microstructure survey (which occurred approximately 100 km away from the mooring site) for direct comparison with microstructure results. The finescale parameterization worked well, with discrepancies ranging from a factor of 1–2.5 depending on depth. The largest discrepancies occurred at depths with high shear. Mean Rω was 17, and Rω showed high variability with values ranging from 3 to 50 over 8 days. Observed ε was slightly elevated (factor of 2–3 compared with a later survey of 11 profiles taken over 3 h) from 25 to 125 m following a wind event which occurred at the beginning of the mooring deployment, reaching a maximum of ε= 6 × 10−10 W kg−1 at 30-m depth. Velocity signals associated with near-inertial waves (NIWs) were observed at depths greater than 200 m, where the Atlantic Water mass represents a reservoir of oceanic heat. However, no evidence of elevated ε or heat fluxes was observed in association with NIWs at these depths in either the microstructure survey or the finescale parameterization estimates.

Description: This work was supported by NSF Grants PLR 14-56705 and PLR-1303791 and by NSF Graduate Research Fellowship Grant DGE-1650112.

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 51(1),(2021): 3-17, https://doi.org/10.1175/JPO-D-20-0064.1.

Description: The strong El Niño of 2014–16 was observed west of the Galápagos Islands through sustained deployment of underwater gliders. Three years of observations began in October 2013 and ended in October 2016, with observations at longitudes 93° and 95°W between latitudes 2°N and 2°S. In total, there were over 3000 glider-days of data, covering over 50 000 km with over 12 000 profiles. Coverage was superior closer to the Galápagos on 93°W, where gliders were equipped with sensors to measure velocity as well as temperature, salinity, and pressure. The repeated glider transects are analyzed to produce highly resolved mean sections and maps of observed variables as functions of time, latitude, and depth. The mean sections reveal the structure of the Equatorial Undercurrent (EUC), the South Equatorial Current, and the equatorial front. The mean fields are used to calculate potential vorticity Q and Richardson number Ri. Gradients in the mean are strong enough to make the sign of Q opposite to that of planetary vorticity and to have Ri near unity, suggestive of mixing. Temporal variability is dominated by the 2014–16 El Niño, with the arrival of depressed isopycnals documented in 2014 and 2015. Increases in eastward velocity advect anomalously salty water and are uncorrelated with warm temperatures and deep isopycnals. Thus, vertical advection is important to changes in heat, and horizontal advection is relevant to changes in salt. Implications of this work include possibilities for future research, model assessment and improvement, and sustained observations across the equatorial Pacific.

Description: We gratefully acknowledge the support of the National Science Foundation (OCE-1232971, OCE-1233282) and the Ocean Observing and Monitoring Division of the National Oceanographic and Atmospheric Administration (NA13OAR4830216).

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of the Atmospheric and Oceanic Technology 38(1), (2021): 3-16, https://doi.org/10.1175/JTECH-D-20-0110.1.

Description: Airborne expendable bathythermographs (AXBTs) are air-launched, single-use temperature–depth probes that telemeter temperature observations as VHF-modulated frequencies. This study describes the AXBT Real-Time Editing System (ARES), which is composed of two components: the ARES Data Acquisition System, which receives telemetered temperature–depth profiles with no external hardware other than a VHF radio receiver, and the ARES Profile Editing System, which quality controls AXBT temperature–depth profiles. The ARES Data Acquisition System performs fast Fourier transforms on windowed segments of the demodulated signal transmitted from the AXBT. For each segment, temperature is determined from peak frequency and depth from elapsed time since profile start. Valid signals are distinguished from noise by comparing peak signal levels and signal-to-noise ratios to predetermined thresholds. When evaluated using 387 profiles, the ARES Data Acquisition System produced temperature–depth profiles nearly identical to those generated using a Sippican MK-21 processor, while reducing the amount of noise from VHF interference included in those profiles. The ARES Profile Editor applies a series of automated checks to identify and correct common profile discrepancies before displaying the profile on an editing interface that provides simple user controls to make additional corrections. When evaluated against 1177 tropical Atlantic and Pacific AXBT profiles, the ARES automated quality control system successfully corrected 87% of the profiles without any required manual intervention. Necessary future work includes improvements to the automated quality control algorithm and algorithm evaluation against a broader dataset of temperature–depth profiles from around the world across all seasons.

Description: This work was sponsored by the Office of Naval Research (Grants N000141812819 and N0001420WX00345) and the U.S. Navy’s Civilian Institution Office with the MIT–WHOI Joint Program.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Spingys, C. P., Garabato, A. C. N., Legg, S., Polzin, K. L., Abrahamsen, E. P., Buckingham, C. E., Forryan, A., & Frajka-Williams, E. E. Mixing and transformation in a deep western boundary current: a case study. Journal of Physical Oceanography, 51(4), (2021): 1205-1222, https://doi.org/10.1175/JPO-D-20-0132.1

Description: Water-mass transformation by turbulent mixing is a key part of the deep-ocean overturning, as it drives the upwelling of dense waters formed at high latitudes. Here, we quantify this transformation and its underpinning processes in a small Southern Ocean basin: the Orkney Deep. Observations reveal a focusing of the transport in density space as a deep western boundary current (DWBC) flows through the region, associated with lightening and densification of the current’s denser and lighter layers, respectively. These transformations are driven by vigorous turbulent mixing. Comparing this transformation with measurements of the rate of turbulent kinetic energy dissipation indicates that, within the DWBC, turbulence operates with a high mixing efficiency, characterized by a dissipation ratio of 0.6 to 1 that exceeds the common value of 0.2. This result is corroborated by estimates of the dissipation ratio from microstructure observations. The causes of the transformation are unraveled through a decomposition into contributions dependent on the gradients in density space of the: dianeutral mixing rate, isoneutral area, and stratification. The transformation is found to be primarily driven by strong turbulence acting on an abrupt transition from the weakly stratified bottom boundary layer to well-stratified off-boundary waters. The reduced boundary layer stratification is generated by a downslope Ekman flow associated with the DWBC’s flow along sloping topography, and is further regulated by submesoscale instabilities acting to restratify near-boundary waters. Our results provide observational evidence endorsing the importance of near-boundary mixing processes to deep-ocean overturning, and highlight the role of DWBCs as hot spots of dianeutral upwelling.

Description: CS, ACNG, AF, and EFW were supported by the U.K. Natural Environment Research Council (NERC) Grant NE/K013181/1. ACNG was supported by the Royal Society and Wolfson Foundation. EPA and CEB were supported by NERC Grant NE/K012843/1. CEB was funded by an MSCA grant (No. 798319) from the European Union’s Horizon 2020 program. EPA was supported by NERC Grant NE/N018095/1. SL and KP were supported by U.S. National Science Foundation Grants OCE-1536453 and OCE-1536779. SL acknowledges support of Award NA18OAR4320123 from the National Oceanic and Atmospheric Administration, U.S. Department of Commerce. The statements, findings, conclusions, and recommendations are those of the authors, and do not necessarily reflect the views of the National Oceanic and Atmospheric Administration, or the U.S. Department of Commerce.

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 51(3), (2021): 955–973, https://doi.org/10.1175/JPO-D-20-0240.1.

Description: Fresh Arctic waters flowing into the Atlantic are thought to have two primary fates. They may be mixed into the deep ocean as part of the overturning circulation, or flow alongside regions of deep water formation without impacting overturning. Climate models suggest that as increasing amounts of freshwater enter the Atlantic, the overturning circulation will be disrupted, yet we lack an understanding of how much freshwater is mixed into the overturning circulation’s deep limb in the present day. To constrain these freshwater pathways, we build steady-state volume, salt, and heat budgets east of Greenland that are initialized with observations and closed using inverse methods. Freshwater sources are split into oceanic Polar Waters from the Arctic and surface freshwater fluxes, which include net precipitation, runoff, and ice melt, to examine how they imprint the circulation differently. We find that 65 mSv (1 Sv ≡ 106 m3 s−1) of the total 110 mSv of surface freshwater fluxes that enter our domain participate in the overturning circulation, as do 0.6 Sv of the total 1.2 Sv of Polar Waters that flow through Fram Strait. Based on these results, we hypothesize that the overturning circulation is more sensitive to future changes in Arctic freshwater outflow and precipitation, while Greenland runoff and iceberg melt are more likely to stay along the coast of Greenland.

Description: We gratefully acknowledge the U.S. National Science Foundation: this work was supported by Grants OCE-1258823, OCE-1756272, OCE-1948335, and OCE-2038481. L.H.S. thanks the U.S. Norway Fulbright Foundation for the Norwegian Arctic Chair Grant 2019-20 that made the visit to Scripps Institution of Oceanography possible. N.P.H. acknowledges support by the U.K. Natural Environment Research Council (NERC) National Capability program CLASS (NE/R015953/1), and Grants U.K.-OSNAP (NE/K010875/1, NE/K010875/2) and U.K.-OSNAP Decade (NE/T00858X/1). We acknowledge the World Climate Research Programme, which, through its Working Group on Coupled Modelling, coordinated and promoted CMIP6.

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(4), (2022): 597–616, https://doi.org/10.1175/jpo-d-21-0121.1.

Description: We provide a first-principles analysis of the energy fluxes in the oceanic internal wave field. The resulting formula is remarkably similar to the renowned phenomenological formula for the turbulent dissipation rate in the ocean, which is known as the finescale parameterization. The prediction is based on the wave turbulence theory of internal gravity waves and on a new methodology devised for the computation of the associated energy fluxes. In the standard spectral representation of the wave energy density, in the two-dimensional vertical wavenumber–frequency (m–ω) domain, the energy fluxes associated with the steady state are found to be directed downscale in both coordinates, closely matching the finescale parameterization formula in functional form and in magnitude. These energy transfers are composed of a “local” and a “scale-separated” contributions; while the former is quantified numerically, the latter is dominated by the induced diffusion process and is amenable to analytical treatment. Contrary to previous results indicating an inverse energy cascade from high frequency to low, at odds with observations, our analysis of all nonzero coefficients of the diffusion tensor predicts a direct energy cascade. Moreover, by the same analysis fundamental spectra that had been deemed “no-flux” solutions are reinstated to the status of “constant-downscale-flux” solutions. This is consequential for an understanding of energy fluxes, sources, and sinks that fits in the observational paradigm of the finescale parameterization, solving at once two long-standing paradoxes that had earned the name of “oceanic ultraviolet catastrophe.”

Description: The authors gratefully acknowledge support from the ONR Grant N00014-17-1-2852. YL gratefully acknowledges support from NSF DMS Award 2009418.

Description: 2022-09-25

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of the Atmospheric and Oceanic Technology 39(2), (2022): 223–235, https://doi.org/10.1175/JTECH-D-21-0110.1.

Description: Previous work with simulations of oceanographic high-frequency (HF) radars has identified possible improvements when using maximum likelihood estimation (MLE) for direction of arrival; however, methods for determining the number of emitters (here defined as spatially distinct patches of the ocean surface) have not realized these improvements. Here we describe and evaluate the use of the likelihood ratio (LR) for emitter detection, demonstrating its application to oceanographic HF radar data. The combined detection–estimation methods MLE-LR are compared with multiple signal classification method (MUSIC) and MUSIC parameters for SeaSonde HF radars, along with a method developed for 8-channel systems known as MUSIC-Highest. Results show that the use of MLE-LR produces similar accuracy, in terms of the RMS difference and correlation coefficients squared, as previous methods. We demonstrate that improved accuracy can be obtained for both methods, at the cost of fewer velocity observations and decreased spatial coverage. For SeaSondes, accuracy improvements are obtained with less commonly used parameter sets. The MLE-LR is shown to be able to resolve simultaneous closely spaced emitters, which has the potential to improve observations obtained by HF radars operating in complex current environments.

Description: This work was supported by the National Science Foundation (NSF) under Grant OCE-1658475. Computing resources were provided by the UCSB Center for Scientific Computing through an NSF MRSEC (DMR-1720256) and NSF CNS-1725797.

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(1),(2022): 75–97, https://doi.org/10.1175/JPO-D-21-0099.1.

Description: Mesoscale eddies contain the bulk of the ocean’s kinetic energy (KE), but fundamental questions remain on the cross-scale KE transfers linking eddy generation and dissipation. The role of submesoscale flows represents the key point of discussion, with contrasting views of submesoscales as either a source or a sink of mesoscale KE. Here, the first observational assessment of the annual cycle of the KE transfer between mesoscale and submesoscale motions is performed in the upper layers of a typical open-ocean region. Although these diagnostics have marginal statistical significance and should be regarded cautiously, they are physically plausible and can provide a valuable benchmark for model evaluation. The cross-scale KE transfer exhibits two distinct stages, whereby submesoscales energize mesoscales in winter and drain mesoscales in spring. Despite this seasonal reversal, an inverse KE cascade operates throughout the year across much of the mesoscale range. Our results are not incompatible with recent modeling investigations that place the headwaters of the inverse KE cascade at the submesoscale, and that rationalize the seasonality of mesoscale KE as an inverse cascade-mediated response to the generation of submesoscales in winter. However, our findings may challenge those investigations by suggesting that, in spring, a downscale KE transfer could dampen the inverse KE cascade. An exploratory appraisal of the dynamics governing mesoscale–submesoscale KE exchanges suggests that the upscale KE transfer in winter is underpinned by mixed layer baroclinic instabilities, and that the downscale KE transfer in spring is associated with frontogenesis. Current submesoscale-permitting ocean models may substantially understate this downscale KE transfer, due to the models’ muted representation of frontogenesis.

Description: The OSMOSIS experiment was funded by the U.K. Natural Environment Research Council (NERC) through Grants NE/1019999/1 and NE/101993X/1. ACNG acknowledges the support of the Royal Society and the Wolfson Foundation, and XY that of a China Scholarship Council PhD studentship.

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(8), (2022): 1677-1691, https://doi.org/10.1175/jpo-d-21-0269.1.

Description: Oceanic mesoscale motions including eddies, meanders, fronts, and filaments comprise a dominant fraction of oceanic kinetic energy and contribute to the redistribution of tracers in the ocean such as heat, salt, and nutrients. This reservoir of mesoscale energy is regulated by the conversion of potential energy and transfers of kinetic energy across spatial scales. Whether and under what circumstances mesoscale turbulence precipitates forward or inverse cascades, and the rates of these cascades, remain difficult to directly observe and quantify despite their impacts on physical and biological processes. Here we use global observations to investigate the seasonality of surface kinetic energy and upper-ocean potential energy. We apply spatial filters to along-track satellite measurements of sea surface height to diagnose surface eddy kinetic energy across 60–300-km scales. A geographic and scale-dependent seasonal cycle appears throughout much of the midlatitudes, with eddy kinetic energy at scales less than 60 km peaking 1–4 months before that at 60–300-km scales. Spatial patterns in this lag align with geographic regions where an Argo-derived estimate of the conversion of potential to kinetic energy is seasonally varying. In midlatitudes, the conversion rate peaks 0–2 months prior to kinetic energy at scales less than 60 km. The consistent geographic patterns between the seasonality of potential energy conversion and kinetic energy across spatial scale provide observational evidence for the inverse cascade and demonstrate that some component of it is seasonally modulated. Implications for mesoscale parameterizations and numerical modeling are discussed.

Description: This work was generously funded by NSF Grants OCE-1912302, OCE-1912125 (Drushka), and OCE-1912325 (Abernathey) as part of the Ocean Energy and Eddy Transport Climate Process Team.

Description: This paper contains twelve charts of the North Atlantic Ocean west of the 60th meridian which show the average temperature of the surface water for each month of the year. Two additional charts show respectively the maximum and minimum values reached by the average temperature for any month. The next chart presents the yearly range, that is, the difference between the average temperature of the warmest and coldest months. The last chart of the series shows the depth of the virtually isothermal water during the summer and winter months.

Description: Author Posting. © American Meteorological Society, 2016. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 46 (2016): 1717-1734, doi:10.1175/JPO-D-15-0124.1.

Description: The contribution of warm-core anticyclones shed by the Irminger Current off West Greenland, known as Irminger rings, to the restratification of the upper layers of the Labrador Sea is investigated in the 1/12° Family of Linked Atlantic Models Experiment (FLAME) model. The model output, covering the 1990–2004 period, shows strong similarities to observations of the Irminger Current as well as ring observations at a mooring located offshore of the eddy formation region in 2007–09. An analysis of fluxes in the model shows that while the majority of heat exchange with the interior indeed occurs at the site of the Irminger Current instability, the contribution of the coherent Irminger rings is modest (18%). Heat is provided to the convective region mainly through noncoherent anomalies and enhanced local mixing by the rings facilitating further exchange between the boundary and interior. The time variability of the eddy kinetic energy and the boundary to interior heat flux in the model are strongly correlated to the density gradient between the dense convective region and the more buoyant boundary current. In FLAME, the density variations of the boundary current are larger than those of the convective region, thereby largely controlling changes in lateral fluxes. Synchronous long-term trends in temperature in the boundary and the interior over the 15-yr simulation suggest that the heat flux relative to the temperature of the interior is largely steady on these time scales.

Description: The authors were supported in this work by the U.S. National Science Foundation.

Description: Author Posting. © American Meteorological Society, 2016. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 46 (2016): 1823-1837, doi:10.1175/JPO-D-15-0165.1.

Description: Measurements just beneath the ocean surface demonstrate that the primary mechanism by which energy from breaking waves is transmitted into the water column is through the work done by the covariance of turbulent pressure and velocity fluctuations. The convergence in the vertical transport of turbulent kinetic energy (TKE) balances the dissipation rate of TKE at first order and is nearly an order of magnitude greater than the sum of the integrated Eulerian and Stokes shear production. The measured TKE transport is consistent with a simple conceptual model that assumes roughly half of the surface flux of TKE by wave breaking is transmitted to depths greater than the significant wave height. During conditions when breaking waves are inferred, the direction of momentum flux is more aligned with the direction of wave propagation than with the wind direction. Both the energy and momentum fluxes occur at frequencies much lower than the wave band, consistent with the time scales associated with wave breaking. The largest instantaneous values of momentum flux are associated with strong downward vertical velocity perturbations, in contrast to the pressure work, which is associated with strong drops in pressure and upward vertical velocity perturbations.

Description: Funding for this research was provided by the National Science Foundation Grants OCE-1339032 and OCE-1338518

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution, August 1969.

Description: A series of laboratory experiments were performed to ascertain the extent and manner of turbulence modification induced by low concentration suspended loads of near neutral buoyancy. Hot-wire measurements of the fluid velocity field of free surface flows were obtained in a specially designed flume recirculating a dielectric liquid and 0.5 mm" diameter spherical plastic particles. In the fixed Reynolds number flow (16,800) the data obtained at six different concentration levels, ranging from 0 to 3.5% by volume, indicate that the presence of particles produces substantial turbulence changes. Even at this low level the mean velocity profile shows an increasing gradient near the bed and sharp deviation from a logarithmic profile. The rms level of each of the velocity components u', v' and w' increases, indicating a general rise in turbulence intensity. The Reynolds stress ρ u'v' increases, and its maximum value shifts away from the bed. The overall scale of turbulence appears to remain unchanged. The data indicate that offhand neglect of suspended particle presence is an oversimplification. There is a similarity between these data and those obtained under adverse pressure gradients. Some effort is made to clarify the altered turbulence production meechanism, and some future experimental work is proposed.

Description: This work was supported by the National Science Foundation under grant No. GA-1582.

Description: Author Posting. © American Meteorological Society, 2014. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 44 (2014): 1306–1328, doi:10.1175/JPO-D-12-0191.1.

Description: 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.

Description: Support for this study and the overall ITP program was provided by the National Science Foundation and Woods Hole Oceanographic Institution. Support for S. Cole was partially though the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Devonshire Foundation.

Description: 2014-11-01

Description: Author Posting. © American Meteorological Society, 2016. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 29 (2016): 1545-1571, doi:10.1175/JCLI-D-15-0509.1.

Description: Three sediment records of sea surface temperature (SST) are analyzed that originate from distant locations in the North Atlantic, have centennial-to-multicentennial resolution, are based on the same reconstruction method and chronological assumptions, and span the past 15 000 yr. Using recursive least squares techniques, an estimate of the time-dependent North Atlantic SST field over the last 15 kyr is sought that is consistent with both the SST records and a surface ocean circulation model, given estimates of their respective error (co)variances. Under the authors’ assumptions about data and model errors, it is found that the 10°C mixed layer isotherm, which approximately traces the modern Subpolar Front, would have moved by ~15° of latitude southward (northward) in the eastern North Atlantic at the onset (termination) of the Younger Dryas cold interval (YD), a result significant at the level of two standard deviations in the isotherm position. In contrast, meridional movements of the isotherm in the Newfoundland basin are estimated to be small and not significant. Thus, the isotherm would have pivoted twice around a region southeast of the Grand Banks, with a southwest–northeast orientation during the warm intervals of the Bølling–Allerød and the Holocene and a more zonal orientation and southerly position during the cold interval of the YD. This study provides an assessment of the significance of similar previous inferences and illustrates the potential of recursive least squares in paleoceanography.

Description: OM acknowledges support from the U.S. National Science Foundation. CW acknowledges support from the European Research Council ERC Grant ACCLIMATE 339108.

Description: 2016-08-19

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 51(10), (2021): 3235–3252, https://doi.org/10.1175/JPO-D-20-0288.1.

Description: Recent mooring measurements from the Overturning in the Subpolar North Atlantic Program have revealed abundant cyclonic eddies at both sides of Cape Farewell, the southern tip of Greenland. In this study, we present further observational evidence, from both Eulerian and Lagrangian perspectives, of deep cyclonic eddies with intense rotation (ζ/f 〉 1) around southern Greenland and into the Labrador Sea. Most of the observed cyclones exhibit strongest rotation below the surface at 700–1000 dbar, where maximum azimuthal velocities are ~30 cm s−1 at radii of ~10 km, with rotational periods of 2–3 days. The cyclonic rotation can extend to the deep overflow water layer (below 1800 dbar), albeit with weaker azimuthal velocities (~10 cm s−1) and longer rotational periods of about one week. Within the middepth rotation cores, the cyclones are in near solid-body rotation and have the potential to trap and transport water. The first high-resolution hydrographic transect across such a cyclone indicates that it is characterized by a local (both vertically and horizontally) potential vorticity maximum in its middepth core and cold, fresh anomalies in the deep overflow water layer, suggesting its source as the Denmark Strait outflow. Additionally, the propagation and evolution of the cyclonic eddies are illustrated with deep Lagrangian floats, including their detachments from the boundary currents to the basin interior. Taken together, the combined Eulerian and Lagrangian observations have provided new insights on the boundary current variability and boundary–interior exchange over a geographically large scale near southern Greenland, calling for further investigations on the (sub)mesoscale dynamics in the region.

Description: OOI mooring data are based upon work supported by the National Science Foundation under Cooperative Agreement 1743430. S. Zou, A. Bower, and H. Furey gratefully acknowledge the support from the Physical Oceanography Program of the U.S. National Science Foundation Grant OCE-1756361. R.S. Pickart acknowledges support from National Science Foundation Grants OCE-1259618 and OCE-1756361. N. P. Holliday and L. Houpert were supported by NERC programs U.K. OSNAP (NE/K010875) and U.K. OSNAP-Decade (NE/T00858X/1).

Description: Author Posting. © American Meteorological Society , 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Zippel, S. F., Farrar, J. T., Zappa, C. J., Miller, U., St Laurent, L., Ijichi, T., Weller, R. A., McRaven, L., Nylund, S., & Le Bel, D. Moored turbulence measurements using pulse-coherent doppler sonar. Journal of Atmospheric and Oceanic Technology, 38(9), (2021): 1621–1639, https://doi.org/10.1175/JTECH-D-21-0005.1.

Description: Upper-ocean turbulence is central to the exchanges of heat, momentum, and gases across the air–sea interface and therefore plays a large role in weather and climate. Current understanding of upper-ocean mixing is lacking, often leading models to misrepresent mixed layer depths and sea surface temperature. In part, progress has been limited by the difficulty of measuring turbulence from fixed moorings that can simultaneously measure surface fluxes and upper-ocean stratification over long time periods. Here we introduce a direct wavenumber method for measuring turbulent kinetic energy (TKE) dissipation rates ϵ from long-enduring moorings using pulse-coherent ADCPs. We discuss optimal programming of the ADCPs, a robust mechanical design for use on a mooring to maximize data return, and data processing techniques including phase-ambiguity unwrapping, spectral analysis, and a correction for instrument response. The method was used in the Salinity Processes Upper-Ocean Regional Study (SPURS) to collect two year-long datasets. We find that the mooring-derived TKE dissipation rates compare favorably to estimates made nearby from a microstructure shear probe mounted to a glider during its two separate 2-week missions for O(10−8) ≤ ϵ ≤ O(10−5) m2 s−3. Periods of disagreement between turbulence estimates from the two platforms coincide with differences in vertical temperature profiles, which may indicate that barrier layers can substantially modulate upper-ocean turbulence over horizontal scales of 1–10 km. We also find that dissipation estimates from two different moorings at 12.5 and at 7 m are in agreement with the surface buoyancy flux during periods of strong nighttime convection, consistent with classic boundary layer theory.

Description: This work was funded by NASA as part of the Salinity Processes in the Upper Ocean Regional Study (SPURS), supporting field work for SPURS-1 (NASA Grant NNX11AE84G), for SPURS-2 (NASA Grant NNX15AG20G), and for analysis (NASA Grant 80NSSC18K1494). Funding for early iterations of this project associated with the VOCALS project and Stratus 9 mooring was provided by NSF (Awards 0745508 and 0745442). Additional funding was provided by ONR Grant N000141812431 and NSF Award 1756839. The Stratus Ocean Reference Station is funded by the Global Ocean Monitoring and Observing Program of the National Oceanic and Atmospheric Administration (CPO FundRef Number 100007298), through the Cooperative Institute for the North Atlantic Region (CINAR) under Cooperative Agreement NA14OAR4320158. Microstructure measurements made from the glider were supported by NSF (Award 1129646).

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(3), (2022): 363–382, https://doi.org/10.1175/jpo-d-21-0084.1.

Description: Meltwater from Greenland is an important freshwater source for the North Atlantic Ocean, released into the ocean at the head of fjords in the form of runoff, submarine melt, and icebergs. The meltwater release gives rise to complex in-fjord transformations that result in its dilution through mixing with other water masses. The transformed waters, which contain the meltwater, are exported from the fjords as a new water mass Glacially Modified Water (GMW). Here we use summer hydrographic data collected from 2013 to 2019 in Upernavik, a major glacial fjord in northwest Greenland, to describe the water masses that flow into the fjord from the shelf and the exported GMWs. Using an optimum multi-parameter technique across multiple years we then show that GMW is composed of 57.8% ± 8.1% Atlantic Water (AW), 41.0% ± 8.3% Polar Water (PW), 1.0% ± 0.1% subglacial discharge, and 0.2% ± 0.2% submarine meltwater. We show that the GMW fractional composition cannot be described by buoyant plume theory alone since it includes lateral mixing within the upper layers of the fjord not accounted for by buoyant plume dynamics. Consistent with its composition, we find that changes in GMW properties reflect changes in the AW and PW source waters. Using the obtained dilution ratios, this study suggests that the exchange across the fjord mouth during summer is on the order of 50 mSv (1 Sv ≡ 106 m3 s−1) (compared to a freshwater input of 0.5 mSv). This study provides a first-order parameterization for the exchange at the mouth of glacial fjords for large-scale ocean models.

Description: This work was partially supported by the Centre for Climate Dynamics (SKD) at the Bjerknes Centre for Climate Research. The authors thank NASA and the OMG consortium for making observational data freely available, and acknowledge M. Morlighem for good support in the early stages of this project. MM and LHS and would also like to thank Ø. Paasche, the ACER project, and the U.S. Norway Fulbright Foundation for the Norwegian Arctic Chair Grant 2019–20 that made the visit to Scripps Institution of Oceanography possible. FS acknowledges support from the DOE Office of Science Grant DE-SC0020073, Heising-Simons Foundation and from NSF and OCE-1756272. DAS acknowledges support from U.K. NERC Grants NE/P011365/1, NE/T011920/1, and NERC Independent Research Fellowship NE/T011920/1. MW was supported by an appointment to the NASA Postdoctoral Program at the Jet Propulsion Laboratory, California Institute of Technology, administered by the Universities Space Research Association under contract with NASA. CSA would like to acknowledge Geocenter Denmark for support to the project “Upernavik Glacier.”

Description: Author Posting. © American Meteorological Society, 2021. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 34(22), (2021): 8971–8987, https://doi.org/10.1175/JCLI-D-20-0610.1.

Description: The impact of increasing Greenland freshwater discharge on the subpolar North Atlantic (SPNA) remains unknown as there are uncertainties associated with the time scales of the Greenland freshwater anomaly (GFWA) in the SPNA. Results from numerical simulations tracking GFWA and an analytical approach are employed to estimate the response time, suggesting that a decadal time scale (13 years) is required for the SPNA to adjust for increasing GFWA. Analytical solutions obtained for a long-lasting increase of freshwater discharge show a non-steady-state response of the SPNA with increasing content of the GFWA. In contrast, solutions for a short-lived pulse of freshwater demonstrate different responses of the SPNA with a rapid increase of freshwater in the domain followed by an exponential decay after the pulse has passed. The derived theoretical relation between time scales shows that residence time scales are time dependent for a non-steady-state case and asymptote the response time scale with time. The residence time of the GFWA deduced from Lagrangian experiments is close to and smaller than the response time, in agreement with the theory. The Lagrangian analysis shows dependence of the residence time on the entrance route of the GFWA and on the depth. The fraction of the GFWA exported through Davis Strait has limited impact on the interior basins, whereas the fraction entering the SPNA from the southwest Greenland shelf spreads into the interior regions. In both cases, the residence time of the GFWA increases with depth demonstrating long persistence of the freshwater anomaly in the subsurface layers.

Description: D. S. Dukhovskoy and E. P. Chassignet were funded by the DOE (Award DE-SC0014378) and HYCOM NOPP (Award N00014-19-1-2674). The HYCOM-CICE simulations were supported by a grant of computer time from the DoD High-Performance Computing Modernization Program at NRL SSC. G. Platov was funded by the RSF N19-17-00154. P. G. Myers was funded by an NSERC Discovery Grant (Grant RGPIN 04357). A. Proshutinsky was funded by FAMOS project (NSF Grant NSF 14-584).

Description: © The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Fine, E., MacKinnon, J., Alford, M., Middleton, L., Taylor, J., Mickett, J., Cole, S., Couto, N., Boyer, A., & Peacock, T. Double diffusion, shear instabilities, and heat impacts of a pacific summer water intrusion in the Beaufort Sea. Journal of Physical Oceanography, 52(2), (2022): 189–203, https://doi.org/10.1175/jpo-d-21-0074.1.

Description: Pacific Summer Water eddies and intrusions transport heat and salt from boundary regions into the western Arctic basin. Here we examine concurrent effects of lateral stirring and vertical mixing using microstructure data collected within a Pacific Summer Water intrusion with a length scale of ∼20 km. This intrusion was characterized by complex thermohaline structure in which warm Pacific Summer Water interleaved in alternating layers of O(1) m thickness with cooler water, due to lateral stirring and intrusive processes. Along interfaces between warm/salty and cold/freshwater masses, the density ratio was favorable to double-diffusive processes. The rate of dissipation of turbulent kinetic energy (ε) was elevated along the interleaving surfaces, with values up to 3 × 10−8 W kg−1 compared to background ε of less than 10−9 W kg−1. Based on the distribution of ε as a function of density ratio Rρ, we conclude that double-diffusive convection is largely responsible for the elevated ε observed over the survey. The lateral processes that created the layered thermohaline structure resulted in vertical thermohaline gradients susceptible to double-diffusive convection, resulting in upward vertical heat fluxes. Bulk vertical heat fluxes above the intrusion are estimated in the range of 0.2–1 W m−2, with the localized flux above the uppermost warm layer elevated to 2–10 W m−2. Lateral fluxes are much larger, estimated between 1000 and 5000 W m−2, and set an overall decay rate for the intrusion of 1–5 years.

Description: This work was supported by ONR Grant N00014-16-1-2378 and NSF Grants PLR 14-56705 and PLR-1303791, NSF Graduate Research Fellowship Grant DGE-1650112, as well as by the Postdoctoral Scholar Program at Woods Hole Oceanographic Institution, with funding provided by the Weston Howland Jr. Postdoctoral Scholarship.

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of the Atmospheric and Oceanic Technology 39(4), (2022): 491–502, https://doi.org/10.1175/jtech-d-21-0046.1.

Description: The Air-Launched Autonomous Micro Observer (ALAMO) is a versatile profiling float that can be launched from an aircraft to make temperature and salinity observations of the upper ocean for over a year with high temporal sampling. Similar in dimensions and weight to an airborne expendable bathythermograph (AXBT), but with the same capability as Argo profiling floats, ALAMOs can be deployed from an A-sized (sonobuoy) launch tube, the stern ramp of a cargo plane, or the door of a small aircraft. Unlike an AXBT, however, the ALAMO float directly measures pressure, can incorporate additional sensors, and is capable of performing hundreds of ocean profiles compared to the single temperature profile provided by an AXBT. Upon deployment, the float parachutes to the ocean, releases the air-deployment package, and immediately begins profiling. Ocean profile data along with position and engineering information are transmitted via the Iridium satellite network, automatically processed, and then distributed by the Global Telecommunications System for use by the operational forecasting community. The ALAMO profiling mission can be modified using the two-way Iridium communications to change the profiling frequency and depth. Example observations are included to demonstrate the ALAMO’s utility.

Description: This work was supported by the National Oceanographic and Atmospheric Administration under Grants NA13OAR4830233 (as part of CINAR Sandy Supplemental funding from the Disaster Relief Appropriations Act of 2013) and NA14OAR4320158 and by Office of Naval Research under Grants N0001416WX01384, N0001416WX01262, and N000141512293. ALAMO floats are commercially available from MRV Systems, LLC (https://www.mrvsys.com).

Description: © The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Semper, S., Vage, K., Pickart, R. S., Valdimarsson, H., Torres, D. J., & Jonsson, S. The emergence of the North Icelandic Jet and its evolution from northeast Iceland to Denmark Strait. Journal of Physical Oceanography, 49(10), (2019): 2499-2521, doi:10.1175/JPO-D-19-0088.1.

Description: The North Icelandic Jet (NIJ) is an important source of dense water to the overflow plume passing through Denmark Strait. The properties, structure, and transport of the NIJ are investigated for the first time along its entire pathway following the continental slope north of Iceland, using 13 hydrographic/velocity surveys of high spatial resolution conducted between 2004 and 2018. The comprehensive dataset reveals that the current originates northeast of Iceland and increases in volume transport by roughly 0.4 Sv (1 Sv ≡ 106 m3 s−1) per 100 km until 300 km upstream of Denmark Strait, at which point the highest transport is reached. The bulk of the NIJ transport is confined to a small area in Θ–S space centered near −0.29° ± 0.16°C in Conservative Temperature and 35.075 ± 0.006 g kg−1 in Absolute Salinity. While the hydrographic properties of this transport mode are not significantly modified along the NIJ’s pathway, the transport estimates vary considerably between and within the surveys. Neither a clear seasonal signal nor a consistent link to atmospheric forcing was found, but barotropic and/or baroclinic instability is likely active in the current. The NIJ displays a double-core structure in roughly 50% of the occupations, with the two cores centered at the 600- and 800-m isobaths, respectively. The transport of overflow water 300 km upstream of Denmark Strait exceeds 1.8 ± 0.3 Sv, which is substantially larger than estimates from a year-long mooring array and hydrographic/velocity surveys closer to the strait, where the NIJ merges with the separated East Greenland Current. This implies a more substantial contribution of the NIJ to the Denmark Strait overflow plume than previously envisaged.

Description: Six different research vessels were involved in the collection of the data used in this study: RRS James Clark Ross, R/V Knorr, R/V Bjarni Sæmundsson, R/V Håkon Mosby, NRV Alliance, and R/V Kristine Bonnevie. We thank the captain and crew of each of these vessels for their hard work as well as the many watch standers who have sailed on the cruises and helped collect the measurements. We also thank Frank Bahr for processing the VMADCP data collected on NRV Alliance and Magnús Danielsen for the processing of the hydrographic data collected on R/V Bjarni Sæmundsson. We acknowledge Leah Trafford McRaven for assistance with Fig. 1 and two anonymous reviewers for their helpful comments, which improved the manuscript. Funding for the project was provided by the Bergen Research Foundation Grant BFS2016REK01 (K. Våge and S. Semper), the Norwegian Research Council under Grant Agreement 231647 (K. Våge), and the U.S. National Science Foundation Grants OCE-1259618 and OCE-1756361 (R. S. Pickart and D. J. Torres), as well as OCE-1558742 (R. S. Pickart). The dataset is available on PANGAEA under https://doi.pangaea.de/10.1594/PANGAEA.903535.

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of [publisher] for personal use, not for redistribution. The definitive version was published in Schlundt, M., Farrar, J. T., Bigorre, S. P., Plueddemann, A. J., & Weller, R. A. (2020). Accuracy of wind observations from open-ocean buoys: correction for flow distortion. Journal of Atmospheric and Oceanic Technology, 37(4), 687-703, doi:10.1175/JTECH-D-19-0132.1.

Description: The comparison of equivalent neutral winds obtained from (i) four WHOI buoys in the subtropics and (ii) scatterometer estimates at those locations reveals a root-mean-square (RMS) difference of 0.56–0.76 m s−1. To investigate this RMS difference, different buoy wind error sources were examined. These buoys are particularly well suited to examine two important sources of buoy wind errors because 1) redundant anemometers and a comparison with numerical flow simulations allow us to quantitatively assess flow distortion errors, and 2) 1-min sampling at the buoys allows us to examine the sensitivity of buoy temporal sampling/averaging in the buoy–scatterometer comparisons. The interanemometer difference varies as a function of wind direction relative to the buoy wind vane and is consistent with the effects of flow distortion expected based on numerical flow simulations. Comparison between the anemometers and scatterometer winds supports the interpretation that the interanemometer disagreement, which can be up to 5% of the wind speed, is due to flow distortion. These insights motivate an empirical correction to the individual anemometer records and subsequent comparison with scatterometer estimates show good agreement.

Description: We gratefully acknowledge the help of three anonymous reviewers, whose input greatly improved the paper. In particular, one reviewer pointed out a mistake in our initial interpretation of scatterometer stability, which was corrected in the final manuscript. JTF and MS were supported by NASA Grant NNX14AM71G (International Ocean Vector Winds Science Team). The SPURS observations were supported by NASA (Grants NNX11AE84G, NNX15AG20G, and 80NSSC18K1494). The Stratus, NTAS, and WHOTS ocean reference stations (ORS) are long-term surface moorings deployed as part of the OceanSITES (http://www.oceansites.org) component of the Global Ocean Observing System, and are supported by NOAA’s Climate Program Office’s Ocean Observing and Monitoring Division, as are RAW, AJP, and SPB through the Cooperative Institute for the North Atlantic Region (CINAR) under Cooperative Agreement NA14OAR4320158 with NOAA Climate Program Office (CPO) (FundRef No. 100007298). The technical staff of the UOP Group at WHOI and the crews of NOAA and UNOLS vessels have been essential to the successful long-term maintenance of the ORS.

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(8), (2020): 2251-2270, doi:10.1175/JPO-D-19-0303.1.

Description: The Gulf Stream affects global climate by transporting water and heat poleward. The current’s volume transport increases markedly along the U.S. East Coast. An extensive observing program using autonomous underwater gliders provides finescale, subsurface observations of hydrography and velocity spanning more than 15° of latitude along the path of the Gulf Stream, thereby filling a 1500-km-long gap between long-term transport measurements in the Florida Strait and downstream of Cape Hatteras. Here, the glider-based observations are combined with shipboard measurements along Line W near 68°W to provide a detailed picture of the along-stream transport increase. To account for the influences of Gulf Stream curvature and adjacent circulation (e.g., corotating eddies) on transport estimates, upper- and lower-bound transports are constructed for each cross–Gulf Stream transect. The upper-bound estimate for time-averaged volume transport above 1000 m is 32.9 ± 1.2 Sv (1 Sv ≡ 106 m3 s−1) in the Florida Strait, 57.3 ± 1.9 Sv at Cape Hatteras, and 75.6 ± 4.7 Sv at Line W. Corresponding lower-bound estimates are 32.3 ± 1.1 Sv in the Florida Strait, 54.5 ± 1.7 Sv at Cape Hatteras, and 69.9 ± 4.2 Sv at Line W. Using the temperature and salinity observations from gliders and Line W, waters are divided into seven classes to investigate the properties of waters that are transported by and entrained into the Gulf Stream. Most of the increase in overall Gulf Stream volume transport above 1000 m stems from the entrainment of subthermocline waters, including upper Labrador Sea Water and Eighteen Degree Water.

Description: We gratefully acknowledge funding from the Office of Naval Research (N000141713040), the National Science Foundation (OCE-0220769, OCE-1633911, OCE-1923362), NOAA’s Global Ocean Monitoring and Observing Program (NA14OAR4320158, NA19OAR4320074), WHOI’s Oceans and Climate Change Institute, Eastman Chemical Company, and the W. Van Alan Clark, Jr. Chair for Excellence in Oceanography at WHOI (awarded to Breck Owens).

Description: Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 50(8),(2020): 2315-2321, doi:10.1175/JPO-D-19-0327.1.

Description: Low-frequency currents and eddies transport sediment, pathogens, larvae, and heat along the coast and between the shoreline and deeper water. Here, low-frequency currents (between 0.1 and 4.0 mHz) observed in shallow surfzone waters for 120 days during a wide range of wave conditions are compared with theories for generation by instabilities of alongshore currents, by ocean-wave-induced sea surface modulations, and by a nonlinear transfer of energy from breaking waves to low-frequency motions via a two-dimensional inverse energy cascade. For these data, the low-frequency currents are not strongly correlated with shear of the alongshore current, with the strength of the alongshore current, or with wave-group statistics. In contrast, on many occasions, the low-frequency currents are consistent with an inverse energy cascade from breaking waves. The energy of the low-frequency surfzone currents increases with the directional spread of the wave field, consistent with vorticity injection by short-crested breaking waves, and structure functions increase with spatial lags, consistent with a cascade of energy from few-meter-scale vortices to larger-scale motions. These results include the first field evidence for the inverse energy cascade in the surfzone and suggest that breaking waves and nonlinear energy transfers should be considered when estimating nearshore transport processes across and along the coast.

Description: Funding was provided by a Vannevar Bush Faculty Fellowship [from OUSD(R&E)] and NSF.

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Cusack, J. M., Brearley, J. A., Garabato, A. C. N., Smeed, D. A., Polzin, K. L., Velzeboer, N., & Shakespeare, C. J. Observed eddy-internal wave interactions in the Southern Ocean. Journal of Physical Oceanography, 50(10), (2020): 3042-3062, doi:10.1175/JPO-D-20-0001.1.

Description: The physical mechanisms that remove energy from the Southern Ocean’s vigorous mesoscale eddy field are not well understood. One proposed mechanism is direct energy transfer to the internal wave field in the ocean interior, via eddy-induced straining and shearing of preexisting internal waves. The magnitude, vertical structure, and temporal variability of the rate of energy transfer between eddies and internal waves is quantified from a 14-month deployment of a mooring cluster in the Scotia Sea. Velocity and buoyancy observations are decomposed into wave and eddy components, and the energy transfer is estimated using the Reynolds-averaged energy equation. We find that eddies gain energy from the internal wave field at a rate of −2.2 ± 0.6 mW m−2, integrated from the bottom to 566 m below the surface. This result can be decomposed into a positive (eddy to wave) component, equal to 0.2 ± 0.1 mW m−2, driven by horizontal straining of internal waves, and a negative (wave to eddy) component, equal to −2.5 ± 0.6 mW m−2, driven by vertical shearing of the wave spectrum. Temporal variability of the transfer rate is much greater than the mean value. Close to topography, large energy transfers are associated with low-frequency buoyancy fluxes, the underpinning physics of which do not conform to linear wave dynamics and are thereby in need of further research. Our work suggests that eddy–internal wave interactions may play a significant role in the energy balance of the Southern Ocean mesoscale eddy and internal wave fields.

Description: Funding for DIMES was provided by U.K. Natural Environment Research Council (NERC) Grants NE/E007058/1 and NE/E005667/1. JMC acknowledges the support of a NERC PhD studentship, and ACNG that of the Royal Society and the Wolfson Foundation. NV acknowledges support from the ARC Centre of Excellence for Climate Extremes (CLEX) Honours Scholarship and the ANU PBSA Partnership - Spotless Scholarship. CJS acknowledges support from an ARC Discovery Early Career Researcher Award DE180100087 and an Australian National University Futures Scheme award. Numerical simulations were conducted on the National Computational Infrastructure (NCI) facility, Canberra, Australia. This study has been conducted using E.U. Copernicus Marine Service Information. We thank two anonymous reviewers for their comments which helped to improve the manuscript significantly. Codes and output files are available online at the project repository (https://github.com/jessecusack/DIMES_eddy_wave_interactions).

Description: Author Posting. © American Meteorological Society, 2019. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 32(13), (2019): 3883-3898, doi:10.1175/JCLI-D-18-0735.1.

Description: While it has generally been understood that the production of Labrador Sea Water (LSW) impacts the Atlantic meridional overturning circulation (MOC), this relationship has not been explored extensively or validated against observations. To explore this relationship, a suite of global ocean–sea ice models forced by the same interannually varying atmospheric dataset, varying in resolution from non-eddy-permitting to eddy-permitting (1°–1/4°), is analyzed to investigate the local and downstream relationships between LSW formation and the MOC on interannual to decadal time scales. While all models display a strong relationship between changes in the LSW volume and the MOC in the Labrador Sea, this relationship degrades considerably downstream of the Labrador Sea. In particular, there is no consistent pattern among the models in the North Atlantic subtropical basin over interannual to decadal time scales. Furthermore, the strong response of the MOC in the Labrador Sea to LSW volume changes in that basin may be biased by the overproduction of LSW in many models compared to observations. This analysis shows that changes in LSW volume in the Labrador Sea cannot be clearly and consistently linked to a coherent MOC response across latitudes over interannual to decadal time scales in ocean hindcast simulations of the last half century. Similarly, no coherent relationships are identified between the MOC and the Labrador Sea mixed layer depth or the density of newly formed LSW across latitudes or across models over interannual to decadal time scales.

Description: FL and MSL are thankful for the financial support from the National Science Foundation (NSF) Physical Oceanography Program (NSF-OCE-12-59102, NSF-OCE-12-59103). The NCAR contribution was supported by the National Oceanic and Atmospheric Administration (NOAA) Climate Program Office (CPO) under Climate Variability and Predictability Program (CVP) Grant NA13OAR4310138 and by the NSF Collaborative Research EaSM2 Grant OCE-1243015. NCAR is sponsored by the NSF. NPH is supported by NERC programs U.K. OSNAP (NE/K010875) and ACSIS (National Capability, NE/N018044/1). Y-OK is supported by NOAA CPO CVP (NA17OAR4310111) and NSF EaSM2 grant (OCE-1242989). AR is supported by NASA-ROSES Modeling, Analysis and Prediction 2016 NNX16AC93G-MAP. RZ is supported by NOAA/OAR. Argo data were collected and made freely available by the International Argo Program and the national programs that contribute to it (http://www.argo.ucsd.edu, http://argo.jcommops.org). The Argo Program is part of the Global Ocean Observing System (http://doi.org/10.17882/42182). Data from the RAPID-MOCHA-WBTS array funded by NERC, NSF and NOAA are freely available from www.rapid.ac.uk/rapidmoc. We thank Stephen Griffies for providing access to the GFDL-MOM025 COREII simulation output and Matthew Harrison and Xiaoqin Yan for their comments on the manuscript. We also thank the anonymous reviewers for their valuable comments.

Description: 2020-06-11