ALBERT

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): 2910-2925, doi:10.1175/2009JPO4139.1.

Description: The propagation of Rossby waves on a midlatitude β plane is investigated in the presence of density diffusion with the aid of linear hydrostatic theory. The search for wave solutions in a vertically bounded medium subject to horizontal (vertical) diffusion leads to an eigenvalue problem of second (fourth) order. Exact solutions of the problem are obtained for uniform background stratification (N), and approximate solutions are constructed for variable N using the Wentzel–Kramers–Brillouin method. Roots of the eigenvalue relations for free waves are found and discussed. The barotropic wave of adiabatic theory is also a solution of the eigenvalue problem as this is augmented with density diffusion in the horizontal or vertical direction. The barotropic wave is undamped as fluid parcels in the wave move only horizontally and are therefore insensitive to the vortex stretching induced by mixing. On the other hand, density diffusion modifies the properties of baroclinic waves of adiabatic theory. In the presence of horizontal diffusion the baroclinic modes are damped but their vertical structure remains unaltered. The ability of horizontal diffusion to damp baroclinic waves stems from its tendency to counteract the deformation of isopycnal surfaces caused by the passage of these waves. The damping rate increases (i) linearly with horizontal diffusivity and (ii) nonlinearly with horizontal wavenumber and mode number. In the presence of vertical diffusion the baroclinic waves suffer both damping and a change in vertical structure. In the long-wave limit the damping is critical (wave decay rate numerically equal to wave frequency) and increases as the square roots of vertical diffusivity and zonal wavenumber. Density diffusion in the horizontal or vertical direction reduces the amplitude of the phase speed of westward-propagating waves. Observational estimates of eddy diffusivities suggest that horizontal and vertical mixing strongly attenuates baroclinic waves in the ocean but that vertical mixing is too weak to notably modify the vertical structure of the gravest modes.

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

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, 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): 1091-1106, doi:10.1175/2007JPO3805.1.

Description: A model of deep ocean circulation driven by turbulent mixing is produced in a long, rectangular laboratory tank. The salinity difference is substituted for the thermal difference between tropical and polar regions. Freshwater gently flows in at the top of one end, dense water enters at the same rate at the top of the other end, and an overflow in the middle removes the same amount of surface water as is pumped in. Mixing is provided by a rod extending from top to bottom of the tank and traveling back and forth at constant speed with Reynolds numbers 〉500. A stratified upper layer (“thermocline”) deepens from the mixing and spreads across the entire tank. Simultaneously, a turbulent plume (“deep ocean overflow”) from a dense-water source descends through the layer and supplies bottom water, which spreads over the entire tank floor and rises into the upper layer to arrest the upper-layer deepening. Data are taken over a wide range of parameters and compared to scaling theory, energetic considerations, and simple models of turbulently mixed fluid. There is approximate agreement with a simple theory for Reynolds number 〉1000 in experiments with a tank depth less than the thermocline depth. A simple argument shows that mixing and plume potential energy flux rates are equal in magnitude, and it is suggested that the same is approximately true for the ocean.

Description: The research was supported by the Ocean Climate Change Institute of Woods Hole Oceanographic Institution.

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): 1012–1021, doi:10.1175/JPO-D-11-0184.1.

Description: Pacific Water flows across the shallow Chukchi Sea before reaching the Arctic Ocean, where it is a source of heat, freshwater, nutrients, and carbon. A substantial portion of Pacific Water is routed through Barrow Canyon, located in the northeast corner of the Chukchi. Barrow Canyon is a region of complex geometry and forcing where a variety of water masses have been observed to coexist. These factors contribute to a dynamic physical environment, with the potential for significant water mass transformation. The measurements of turbulent kinetic energy dissipation presented here indicate diapycnal mixing is important in the upper canyon. Elevated dissipation rates were observed near the pycnocline, effectively mixing winter and summer water masses, as well as within the bottom boundary layer. The slopes of shear/stratification layers, combined with analysis of rotary spectra, suggest that near-inertial wave activity may be important in modulating dissipation near the bottom. Because the canyon is known to be a hotspot of productivity with an active benthic community, mixing may be an important factor in maintenance of the biological environment.

Description: ELS was supported as a WHOI Postdoctoral Scholar through the WHOI Ocean and Climate Change Institute.

Description: 2012-12-01

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 Atmospheric and Oceanic Technology 30 (2013): 1767–1788, doi:10.1175/JTECH-D-12-00140.1.

Description: Seismic images of oceanic thermohaline finestructure record vertical displacements from internal waves and turbulence over large sections at unprecedented horizontal resolution. Where reflections follow isopycnals, their displacements can be used to estimate levels of turbulence dissipation, by applying the Klymak–Moum slope spectrum method. However, many issues must be considered when using seismic images for estimating turbulence dissipation, especially sources of random and harmonic noise. This study examines the utility of seismic images for estimating turbulence dissipation in the ocean, using synthetic modeling and data from two field surveys, from the South China Sea and the eastern Pacific Ocean, including the first comparison of turbulence estimates from seismic images and from vertical shear. Realistic synthetic models that mimic the spectral characteristics of internal waves and turbulence show that reflector slope spectra accurately reproduce isopycnal slope spectra out to horizontal wavenumbers of 0.04 cpm, corresponding to horizontal wavelengths of 25 m. Using seismic reflector slope spectra requires recognition and suppression of shot-generated harmonic noise and restriction of data to frequency bands with signal-to-noise ratios greater than about 4. Calculation of slope spectra directly from Fourier transforms of the seismic data is necessary to determine the suitability of a particular dataset to turbulence estimation from reflector slope spectra. Turbulence dissipation estimated from seismic reflector displacements compares well to those from 10-m shear determined by coincident expendable current profiler (XCP) data, demonstrating that seismic images can produce reliable estimates of turbulence dissipation in the ocean, provided that random noise is minimal and harmonic noise is removed.

Description: This work was funded by NSF Grants 0452744, 0405654, and 0648620, and ONR/DEPSCoR Grant DODONR40027.

Description: 2014-02-01

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): 909–917, doi:10.1175/2007JPO3535.1.

Description: The classical two-box model of Stommel is extended in two directions: replacing the buoyancy constraint with an energy constraint and including the wind-driven gyre. Stommel postulated a buoyancy constraint for the thermohaline circulation, and his basic idea has evolved into the dominating theory of thermohaline circulation; however, recently, it is argued that the thermohaline circulation is maintained by mechanical energy from wind stress and tides. The major difference between these two types of models is the bifurcation structure: the Stommel-like model has two thermal modes (one stable and another one unstable) and one stable haline mode, whereas the energy-constraint model has one stable thermal mode and two saline modes (one stable and another one unstable). Adding the wind-driven gyre changes the threshold value of thermohaline bifurcation greatly; thus, the inclusion of the wind-driven gyre is a vital step in completely modeling the physical processes related to thermohaline circulation.

Description: YPG was supported by the National Science Foundation of China (NSFC, 40676022), the National Basic Research Program of China (2006CB403605), and the Guangdong Natural Science Foundation (5003672). RXH was supported by the National Oceanic and Atmospheric Administration through CICOR Cooperative Agreement NA17RJ1223 to the Woods Hole Oceanographic Institution.

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.