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 February 2009

Description: Observations and inverse models suggest that small-scale turbulent mixing is enhanced in the Southern Ocean in regions above rough topography. The enhancement extends 1 km above the topography suggesting that mixing is supported by breaking of gravity waves radiated from the ocean bottom. In other regions, gravity wave radiation by bottom topography has been primarily associated with the barotropic tide. In this study, we explore the alternative hypothesis that the enhanced mixing in the Southern Ocean is sustained by internal waves generated by geostrophic motions flowing over bottom topography. Weakly-nonlinear theory is used to describe the internal wave generation and the feedback of the waves on the zonally averaged flow. A major finding is that the waves generated at the ocean bottom at finite inverse Froude numbers drive vigorous inertial oscillations. The wave radiation and dissipation at equilibrium is therefore the result of both geostrophic flow and inertial oscillations and differs substantially from the classical lee wave problem. The theoretical predictions are tested versus two-dimensional and three-dimensional high resolution numerical simulations with parameters representative of the Drake Passage region. Theory and fully nonlinear numerical simulations are used to estimate internal wave radiation from LADCP, CTD and topography data from two regions in the Southern Ocean: Drake Passage and the Southeast Pacific. The results show that radiation and dissipation of internal waves generated by geostrophic motions reproduce the magnitude and distribution of dissipation measured in the region.

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): 789-801, doi:10.1175/2009JPO4039.1.

Description: The issue of internal wave–mesoscale eddy interactions is revisited. Previous observational work identified the mesoscale eddy field as a possible source of internal wave energy. Characterization of the coupling as a viscous process provides a smaller horizontal transfer coefficient than previously obtained, with vh 50 m2 s−1 in contrast to νh 200–400 m2 s−1, and a vertical transfer coefficient bounded away from zero, with νυ + (f2/N2)Kh 2.5 ± 0.3 × 10−3 m2 s−1 in contrast to νυ + (f2/N2)Kh = 0 ± 2 × 10−2 m2 s−1. Current meter data from the Local Dynamics Experiment of the PolyMode field program indicate mesoscale eddy–internal wave coupling through horizontal interactions (i) is a significant sink of eddy energy and (ii) plays an O(1) role in the energy budget of the internal wave field.

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): 2556-2574, doi:10.1175/2008JPO3666.1.

Description: Vertical profiles of horizontal velocity obtained during the Mid-Ocean Dynamics Experiment (MODE) provided the first published estimates of the high vertical wavenumber structure of horizontal velocity. The data were interpreted as being representative of the background internal wave field, and thus, despite some evidence of excess downward energy propagation associated with coherent near-inertial features that was interpreted in terms of atmospheric generation, these data provided the basis for a revision to the Garrett and Munk spectral model. These data are reinterpreted through the lens of 30 years of research. Rather than representing the background wave field, atmospheric generation, or even near-inertial wave trapping, the coherent high wavenumber features are characteristic of internal wave capture in a mesoscale strain field. Wave capture represents a generalization of critical layer events for flows lacking the spatial symmetry inherent in a parallel shear flow or isolated vortex.

Description: Salary support for this analysis was provided by Woods Hole Oceanographic Institution bridge support funds.

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 2000

Description: The thesis develops and demonstrates methods of classifying ocean processes using an underwater moving platform such as an Autonomous Underwater Vehicle (AUV). The "mingled spectrum principle" is established which concisely relates observations from a moving platform to the frequency-wavenumber spectrum of the ocean process. It clearly reveals the role of the AUV speed in mingling temporal and spatial information. For classifying different processes, an AUV is not only able to jointly utilize the time-space information, but also at a tunable proportion by adjusting its cruise speed. In this respect, AUVs are advantageous compared with traditional oceanographic platforms. Based on the mingled spectrum principle, a parametric tool for designing an AUVbased spectral classifier is developed. An AUV's controllable speed tunes the separability between the mingled spectra of different processes. This property is the key to optimizing the classifier's performance. As a case study, AUV-based classification is applied to distinguish ocean convection from internal waves. The mingled spectrum templates are derived from the MIT Ocean Convection Model and the Garrett-Munk internal wave spectrum model. To allow for mismatch between modeled templates and real measurements, the AUVbased classifier is designed to be robust to parameter uncertainties. By simulation tests on the classifier, it is demonstrated that at a higher AUV speed, convection's distinct spatial feature is highlighted to the advantage of classification. Experimental data are used to test the AUV-based classifier. An AUV-borne flow measurement system is designed and built, using an Acoustic Doppler Velocimeter (ADV). The system is calibrated in a high-precision tow tank. In February 1998, the AUV acquired field data of flow velocity in the Labrador Sea Convection Experiment. The Earth-referenced vertical flow velocity is extracted from the raw measurements. The classification test result detects convection's occurrence, a finding supported by more traditional oceanographic analyses and observations. The thesis work provides an important foundation for future work in autonomous detection and sampling of oceanographic processes.

Description: This thesis research has been funded by the Office of Naval Research (ONR) under Grants NOOOl4-95-1-1316, NOO0l4-97-1-0470, and by the MIT Sea Grant College Program under Grant NA46RG0434.

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): 2605–2623, doi:10.1175/2010JPO4132.1.

Description: Steady scale-invariant solutions of a kinetic equation describing the statistics of oceanic internal gravity waves based on wave turbulence theory are investigated. It is shown in the nonrotating scale-invariant limit that the collision integral in the kinetic equation diverges for almost all spectral power-law exponents. These divergences come from resonant interactions with the smallest horizontal wavenumbers and/or the largest horizontal wavenumbers with extreme scale separations. A small domain is identified in which the scale-invariant collision integral converges and numerically find a convergent power-law solution. This numerical solution is close to the Garrett–Munk spectrum. Power-law exponents that potentially permit a balance between the infrared and ultraviolet divergences are investigated. The balanced exponents are generalizations of an exact solution of the scale-invariant kinetic equation, the Pelinovsky–Raevsky spectrum. A small but finite Coriolis parameter representing the effects of rotation is introduced into the kinetic equation to determine solutions over the divergent part of the domain using rigorous asymptotic arguments. This gives rise to the induced diffusion regime. The derivation of the kinetic equation is based on an assumption of weak nonlinearity. Dominance of the nonlocal interactions puts the self-consistency of the kinetic equation at risk. However, these weakly nonlinear stationary states are consistent with much of the observational evidence.

Description: This research is supported by NSF CMG Grants 0417724, 0417732 and 0417466. YL is also supported by NSF DMS Grant 0807871 and ONR Award N00014-09-1-0515.

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): 1524–1547, doi:10.1175/JPO-D-11-0117.1.

Description: Evidence is presented for the transfer of energy from low-frequency inertial–diurnal internal waves to high-frequency waves in the band between 6 cpd and the buoyancy frequency. This transfer links the most energetic waves in the spectrum, those receiving energy directly from the winds, barotropic tides, and parametric subharmonic instability, with those most directly involved in the breaking process. Transfer estimates are based on month-long records of ocean velocity and temperature obtained continuously over 80–800 m from the research platform (R/P) Floating Instrument Platform (FLIP) in the Hawaii Ocean Mixing Experiment (HOME) Nearfield (2002) and Farfield (2001) experiments, in Hawaiian waters. Triple correlations between low-frequency vertical shears and high-frequency Reynolds stresses, uiw∂Ui/∂z, are used to estimate energy transfers. These are supported by bispectral analysis, which show significant energy transfers to pairs of waves with nearly identical frequency. Wavenumber bispectra indicate that the vertical scales of the high-frequency waves are unequal, with one wave of comparable scale to that of the low-frequency parent and the other of much longer scale. The scales of the high-frequency waves contrast with the classical pictures of induced diffusion and elastic scattering interactions and violates the scale-separation assumption of eikonal models of interaction. The possibility that the observed waves are Doppler shifted from intrinsic frequencies near f or N is explored. Peak transfer rates in the Nearfield, an energetic tidal conversion site, are on the order of 2 × 10−7 W kg−1 and are of similar magnitude to estimates of turbulent dissipation that were made near the ridge during HOME. Transfer rates in the Farfield are found to be about half the Nearfield values.

Description: This work was supported by the National Science Foundation and the Office of Naval Research.

Description: 2013-03-01

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 January 18, 1988

Description: We have contrived a regional model Φ(K, ω, η, φ, λ) for the distribution of low frequency variability energy in horizontal wavenumber, frequency, vertical mode and geography. We assume horizontal isotropy, Φ(K, ω, η, φ, λ) = 2πKψ(k, l, ω, η, φ, λ), with K designating the amplitude of total horizontal wavenumber. The parameters of Φ(K, ω, η, φ, λ) can be derived from observations: (i) satellite altimetry measurements yield the surface eddy kinetic energy wavenumber and frequency spectra and the geographic distribution of surface eddy kinetic energy magnitude, (ii) XBT measurements yield the temperature wavenumber spectra, (iii) current meter and thermistor measurements yield the frequency spectra of kinetic energy and temperature, (iv) tomographic measurements yield the frequency spectra of range— and depth—averaged temperature, and (v) the combination of satellite altimetry and current meter measurements yields the vertical partitioning of kinetic energy among dynamical modes. We assume the form of the geography—independent part of our model Φ(K, ω, η) ∝Kpωq. The observed kinetic energy and temperature wavenumber spectra suggest p = 3/2 at K 〈 K0 and p = —2 at K 〉 K0 for the barotropic mode, and p = —1/2 at K 〈 K0 and p = —3 at K 〉 K0 for the baroclinic mods, where K0 is the transitional wavenumber of the wavenumber spectra. The observed frequency spectra of temperature and kinetic energy suggest that q = —1/2 for ω 〈 ω0 and q = —2 for ω 〉 ω0, where ω0 is the transitional frequency of the frequency spectra. The combination of satellite altimetry and current meter measurements suggests the vertical structure of the low frequency variability is governed by the first few modes. The geography—dependent part of our model is the energy magnitude. Although we have shown analytically that the tomographic measurements behave as a low—pass filter, it is impossible to identify this filtering effect in the real data due to the strong geographic variability of the energy magnitude and the vertical gradient of the mean temperature. The model wavenumber spectrum is appropriate only where the statistical properties are relatively homogeneous in space.

Description: My first year in the Joint Program was supported by the National Science Foundation under grant OCE 92-16628, then were supported by the National Science Foundation under grant OCE 95-29545.

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): 17–28, doi:10.1175/JPO-D-11-0108.1.

Description: Observational evidence is presented for transfer of energy from the internal tide to near-inertial motions near 29°N in the Pacific Ocean. The transfer is accomplished via parametric subharmonic instability (PSI), which involves interaction between a primary wave (the internal tide in this case) and two smaller-scale waves of nearly half the frequency. The internal tide at this location is a complex superposition of a low-mode waves propagating north from Hawaii and higher-mode waves generated at local seamounts, making application of PSI theory challenging. Nevertheless, a statistically significant phase locking is documented between the internal tide and upward- and downward-propagating near-inertial waves. The phase between those three waves is consistent with that expected from PSI theory. Calculated energy transfer rates from the tide to near-inertial motions are modest, consistent with local dissipation rate estimates. The conclusion is that while PSI does befall the tide near a critical latitude of 29°N, it does not do so catastrophically.

Description: This work was sponsored by NSF OCE 04-25283.

Description: 2013-07-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 April, 1977

Description: A total of four moorings from POLYMODE Array I and II were analyzed in an investigation of internal wavefield-mean flow interactions. In particular, evidence for wave-mean flow interaction was sought by searching for time correlations between the wavefield vertically-acting Reynolds stress (estimated using the temperature and velocity records), and the mean shear. No significant stress-shear correlations were found at the less energetic moorings, indicating that the magnitude of the eddy viscosity was under 200 cm2/sec, with the sign of the energy transfer uncertain. This is considerably below the 0(4500 cm2/sec) predicted by Müller (1976). An extensive error analysis indicates that the large wave stress predicted by the theory should have been clearly observable under the conditions of measurement. Theoretical computations indicate that the wavefield "basic state" may not be independent of the mean flow as assumed by Müller, but can actually be modified by large-scale vertical shear and still remain in equilibrium. In that case, the wavefield does not exchange momentum with a large-scale vertical shear flow, and, excepting critical layer effects, a small vertical eddy viscosity is to be expected. Using the Garrett-Munk (1975) model internal wave spectrum, estimates were made of the maximum momentum flux (stress) expected to be lost to critical layer absorption. Stress was found to increase almost linearly with the velocity difference across the shear zone, corresponding to a vertical eddy viscosity of -100 cm2 s -1. Stresses indicative of this effect were not observed in the data. The only significantly non-zero stress correlations were found at the more energetic moorings. Associated with the 600 m mean velocity and the shear at the thermocline were a positively correlated stress at 600 m, and a negatively correlated stress at 1000 m. These stress correlations were most clearly observable in the frequency range corresponding to 1 to 8 hour wave periods. The internal wavefield kinetic and potential energy were modulated by the mean flow at both levels, increasing by a factor of two with a factor of ten in the mean flow. The observed stress correlations and energy level changes were found to be inconsistent with ideas of a strictly local eddy viscosity, in which the spectrum of waves is only slightly modified by the shear. When Doppler effects in the temperature equation used to estimate vertical velocity were considered, the observations of stress and energy changes were found to be consistent with generation of short (0.4 to 3 km) internal waves at the level of maximum shear, about 800 m. The intensity of the generated waves increases with the shear, resulting in an effective vertical eddy viscosity (based on the main thermocline shear) of about +100 cm2 s-1 The stresses were not observable at the 1500 m level, indicating that the waves were absorbed within 500 m of vertical travel. The tendency for internal wave currents to be horizontally anisotropic in the presence of a mean current was investigated. Using the Garrett- Munk (1975) model internal wave spectrum, it was found that critical layer absorption cannot induce anisotropies as large as observed. A mechanical noise problem was found to be the cause of large anisotropies measured with Geodyne 850 current meters. It could not be decided, however, whether or not the A.M.F. Vector Averaging Current Meter is able to satisfactorily remove the noise with its averaging scheme.

Description: The research reported here was provided by Office of Naval Research Contract Numer N00014-76-C-0197 NR 083-400.

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 1984

Description: Intermittent, shoreward propagating packets of high frequency first mode internal waves are common on the continental shelf when the water column is stratified and may induce large fluctuations in near bottom velocity. Simple theoretical considerations here lead to an approximate method for estimating those quantities of most interest for the bottom boundary layer interaction problem. Examination of data from the pilot Coastal Ocean Dynamics Experiment (CODE I) shows that near bottom velocity fluctuations in the high frequency internal wave band were dominated by shoreward propagating, intermittent mode 1 internal events. Predictions of CODE I internal wave characteristics using the above approximate method are shown to be good. A boundary layer model is developed, which allows for the nonlinear interaction of surface waves, internal waves, and a steady current over a rough bottom. Modeling results suggest that internal waves will significantly enhance the stress felt by the steady current, and can increase the variability and decrease the reliability of boundary layer measurements by the "log profile" technique, when the waves are present. Theoretical dissipation of internal wave energy in the bottom boundary layer is found to be significantly enhanced in the presence of surface waves and currents, and may be important to the overall internal wave energy balance on the shelf.

Description: My doctoral work was supported for the first three years by an NSF Graduate Fellowship and has been supported since under NSF grant OCE-8014938.