ALBERT

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