Thesis (S.M.)--Joint Program in Oceanography/Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2012. / Cataloged from PDF version of thesis. / Includes bibliographical references (p. 40-41). / The purpose of this study is to understand the interactions of tropical cyclones with ocean eddies. In particular we examine the influence of a cold-core eddy on the cold wake formed during the passage of Typhoon Fanapi (2010). The three-dimensional version of the numerical Price-Weller-Pinkel (PWP) vertical mixing model has previously been used to simulate and study the cold wakes of Atlantic hurricanes. The model has not been used in comparison with observations of typhoons in the Western Pacific Ocean. In 2010 several typhoons were studied during the Impact of Typhoons on the Ocean in the Pacific (ITOP) field campaign and Fanapi was particularly well observed. We use these observations and the 3DPWP to understand the ocean cold wake generated by Fanapi. The cold wake of Fanapi was advected by a cyclonic eddy that was south of the typhoon track. The 3DPWP model outputs with and without an eddy are compared with observations made during the field campaign. These observations are compared to model outputs with eddies in a series of positions right and left of the storm track in order to study effects of mesoscale eddies on ocean vertical mixing in the cold wake of typhoons. / by Alexander Reid Miltenberger. / S.M.

Thesis: Ph. D., Joint Program in Oceanography/Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2019 / Cataloged from PDF version of thesis. / Includes bibliographical references (pages 113-121). / Submesoscale ocean dynamics and instabilities, with characteristic scales 0.1-10 kin, can play a critical role in setting the ocean's surface boundary layer thickness and associated density stratification. Submesoscale instabilities contribute to lateral stirring and tracer dispersal. These dynamics are investigated in the Bay of Bengal, motivated by the upper ocean's potentially coupled interactions with Monsoon winds and convection. The region's excess precipitation and runoff generates strong salinity gradients that typically set density fronts and stratification in the upper 50 m. Since we cannot synoptically measure currents containing fast-evolving and oscillating components across the submesoscale range, we instead analyze passive tracer distributions (spice = density-compensated temperature (T) and salinity (S) anomalies), identifying signatures of flows and testing dynamical theories. / The analysis is based on over 9000 vertical profiles of T and S measured along ~4800 km of ship tracks in the Bay of Bengal during ASIRI and MISO-BOB expeditions in 2013, 2015, and 2018. Observations in the surface mixed layer reveal ~1 km scale-selective correlation of surface T and S, with compensation reducing cross-front density gradients by ~50%. Using a process study ocean model, we show this is caused by submesoscale instabilities slumping fronts, plus surface cooling over the resultant enhanced salinity stratification, potentially thwarting the forward cascade of energy. In the stratified interior, we present a spectral analysis of horizontal spice variance statistics from wavenumber k ~0.01 cpkm to ~1 cpkm. At scales <10 km, stratified layers that are closer to the surface exhibit redder passive tracer spectra (power spectra k⁻³, gradient spectra k⁻¹) than predicted by quasi-geostrophic or frontogenetic theories. / Complimentary observations reveal spice patterns with multiple, parallel, ~10 m thin layers, crossing isopycnals with O(10⁻⁴) slopes, coherent over at least 30-80 kin, with coincident layers of stratification anomalies. Comparison with shear measurements, and a numerical process study, suggest that both submesoscale sheared eddies, and thin near-inertial waves, form such layers. Fast formation timescales and large aspect ratios suggest they enhance horizontal mixing by shear dispersion, reducing variance at ~1-10 km scales. / by Gualtiero Victor Rudi Spiro Jaeger. / Ph. D. / Ph.D. Joint Program in Oceanography/Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution)

Thesis (S.M.)--Joint Program in Physical Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences and the Woods Hole Oceanographic Institution), 2000. / Includes bibliographical references (leaves 78-82). / A simple quasi-two dimensional dynamical model of Thermohaline circulation (THC) is developed, assuming that the mixing only occurs near western and eastern boundary layers. When the surface density is prescribed, the climatically important quantities, such as the strength of overturning and meridional heat transport, are related to the zonal integral over the vigorously mixing regions and scaled as (KvΔx)2/3. The numerical results suggest that the density difference between eastern and western boundaries play an important role in the meridional overturning. The eastern boundary is characterized by the upwelling on top of downwelling. The western boundary layer is featured by the universal upwelling. The inefficiency of diffusion heat transport accounts for the narrowness of sinking region and shallowness of overturning cell in one-hemisphere. The experiments with other surface boundary conditions are also explored. The circulation patterns obtained are similar under various surface temperature distributions, suggesting these are very robust features of THC. The role of boundary mixing is further explored in global ocean. The 2 1/2 dynamical model is extended to two-hemisphere ocean. Additional dynamics such as Rayleigh friction and abyssal water properties are taken into account. A set of complicated governing equations are derived and numerically solved to obtain steady state solution. The basic circulation features are revealed in our dynamical model. An equtorially asymmetric meridional circulation is observed due to small perturbation at the surface temperature in the high latitude. The density differences between eastern and western boundaries are distinct in both hemispheres. This is achieved during the spin-up process. Although the dynamical model results agree well with OGCM results in one-hemisphere, several important dynamics are lacking and exposed in two-hemisphere experiments. We need to consider horizontal advection terms which will effectively advect positive density anomalies across the equator and form the deep water for the entire system. / by Hua Ru. / S.M.

Thesis (Ph. D.)--Joint Program in Physical Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 1988. / Includes bibliographical references (leaves 163-172). / A freshwater plume often forms when a river or an estuary discharges water onto the continental shelf. Freshwater plumes are ubiquitous features of the coastal ocean and usually leave a striking signature in the coastal hydrography. The present study combines both hydrographic data and idealized numerical simulations to examine how ambient currents and winds influence the transport and mixing of plume waters. The first portion of the thesis considers the alongshore transport of freshwater using idealized numerical simulations. In the absence of any ambient current, the downstream coastal current only carries a fraction of the discharged fresh water; the remaining fraction recirculates in a continually growing "bulge" of fresh water in the vicinity of the river mouth. The fraction of fresh water transported in the coastal current is dependent on the source conditions at the river mouth. The presence of an ambient current augments the transport in the plume so that its freshwater transport matches the freshwater source. For any ambient current in the same direction as the geostrophic coastal current, the plume will evolve to a steady-state width. A key result is that an external forcing agent is required in order for the entire freshwater volume discharged by a river to be transported as a coastal current. The next section of the thesis addresses the wind-induced advection of a river plume, using hydrographic data collected in the western Gulf of Maine. The observations suggest that the plume's cross-shore structure varies markedly as a function of fluctuations in alongshore wind forcing. Consistent with Ekman dynamics, upwelling favorable winds spread the plume offshore, at times widening it to over 50 km in offshore extent, while downwelling favorable winds narrow the plume width to a few Rossby radii. Near-surface current meters show significant correlations between cross-shore currents and alongshore wind stress, consistent with Ekman theory. Estimates of the terms in the alongshore momentum equation calculated from moored current meter arrays also indicate an approximate Ekman balance within the plume. A significant correlation between alongshore currents and alongshore wind stress suggests that interfacial drag may be important. The final section of the thesis is an investigation of the advection and mixing of a surface-trapped river plume in the presence of an upwelling favorable wind stress, using a three-dimensional model in a simple, rectangular domain. Model simulations demonstrate that the plume thins and is advected offshore by the cross shore Ekman transport. The thinned plume is susceptible to significant mixing due to the vertically sheared horizontal currents. The first order plume response is explained by Ekman dynamics and a Richardson number mixing criterion. / by Derek Allen Fong. / Ph.D.

Thesis (S.M.)--Joint Program in Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), February 2001. / Includes bibliographical references (p. 116-118). / I study the exchange between a boundary current and flanking horizontal recirculations in a 'sliced-cylinder' rotating tank laboratory experiment. Two flow configurations are investigated: a single recirculation and a double, figure-8, recirculation. The latter case involves a hyperbolic point, while the former does not. I investigate the stirring and mixing under both steady and unsteady forcing. I quantify the mixing in each case using effective diffusivity, Keff, and a corollary effective length, Leff, as derived by Nakamura (1995, 1996). This approach involves diagnosing the geometric complexity of a tracer field. Geometric complexity is indicative of advective stirring. Because stirring creates high gradients, flows with high advective stirring also have high diffusion, and stronger overall mixing. I calculate effective length from images of dye in the tank and find much higher values of Leff in the unsteady hyperbolic cases than in the other cases. Slight unsteadiness in flows involving hyperbolic points gives rise to a chaotic advection mechanism known as 'lobe dynamics'. These lobes carry fluid in and out of the recirculations, acting as extremely effective stirring mechanisms. I demonstrate the existence of these exchange lobes in the unsteady hyperbolic (figure-8) flow. The velocity field in the tank is calculated utilizing particle image velocimetry (PIV) techniques and a time series U(t) demonstrates the (forced) unsteadiness in the flow. Images of dye in the tank show exchange lobes forming at this same forcing period, and carrying fluid in and out of the recirculation. Based on the results of these experiments, I am able to confirm that, at least in this controlled environment, basic geometry has a profound effect on the mixing effectiveness of a recirculation. I demonstrate radically increased stirring and mixing in the unsteady hyperbolic flow as compared to steady flows and flows without hyperbolic points. Recirculations are ubiquitous in the world ocean; they occur on a variety of scales, in many different configurations, and at all depths. Some of these configurations involve hyperbolic points, while others do not. Chaotic advection via lobe exchange may be an important component of the mixing at multiple locations in the ocean where hyperbolic recirculation geometries exist. / by Heather E. Deese. / S.M.

Thesis (S.M.)--Joint Program in Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 1999. / Includes bibliographical references (leaves 16-21). / Today, deep waters produced in the North Atlantic are exported through the western South Atlantic. Antarctic intermediate water (AAIW) also enters the Atlantic in this region. Circumpolar deep water (CDW) fills the depths below AAIW and above and below northern source waters. A depth transect of cores from 1567-3909 m water depth in the western South Atlantic are ideally located to monitor inter-ocean exchange of deep water, and variations in the relative strength of northern versus southern source water production. Last glacial maximum (LGM) Cd/Ca and 813C data indicate a nutrient-depleted intermediate-depth water mass. In the mid-depth western South Atlantic, a simple conversion of LGM 813C data suggests significantly less nutrient enrichment than LGM Cd/Ca ratios, but Cd/Ca and 613C data can be reconciled when plotted in CdW/ 13C space. Paired LGM Cd/Ca and S13C data from mid-depth cores suggest increasingly nutrient rich waters below 2000 m, but do not require an increase in Southern Ocean water contribution relative to today. Cd/Ca data suggest no glacial-interglacial change in the hydrography of the deepest waters of the region. To maintain relatively low Cd/Ca ratios (low nutrients) in the deepest western South Atlantic waters, and in CDW in general, during the LGM requires an increased supply of nutrient-depleted glacial North Atlantic intermediate water (GNAIW) and/or nutrient-depleted glacial Subantarctic surface waters to CDW to balance reduced NADW contribution to CDW. LGM Cd/Ca and 513C data suggest strong GNAIW influence in the western South Atlantic which in turn implies export of GNAIW from the Atlantic, and entrainment of GNAIW into the Antarctic Circumpolar current. / by Michael Horowitz. / S.M.

Thesis (Ph. D.)--Joint Program in Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), February 2000. / Includes bibliographic references. / The bottom boundary layer is an important dynamical region of shallow water flows. In this thesis, the problem of turbulent mixing in the coastal bottom boundary layer is investigated with a unique set of field measurements of velocity and sound speed that span a significant fraction of the boundary layer obtained over a six-week long period in the late summer of 1996 on the New England shelf. The energetics of the turbulent fluctuations are investigated by testing simplified budgets for turbulent kinetic energy and scalar variance. The turbulent kinetic energy budget is locally balanced while the scalar variance budget is not, probably due to turbulent diffusion. The direct effects of stratification are consistently significant only in the outer part of the boundary layer, where the flux Richardson number is approximately equal to a critical value of 0.2. Turbulence closure is investigated in terms of non-dimensional profiles of velocity and sound speed. Close to the bottom, the results are consistent with Monin-Obukhov similarity theory, while in the outer part of the boundary layer other scales including the height of the boundary layer are important for setting the turbulent length scale. / by William J. Shaw. / Ph.D.

Thesis (Ph. D.)--Joint Program in Oceanography/Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2008. / Includes bibliographical references (p. 110-119). / This study presents observations of turbulence dynamics made during the low winds portion of the Coupled Boundary Layers and Air-Sea Transfer experiment (CBLAST-Low). Observations were made of turbulent fluxes, turbulent kinetic energy, and the length scales of flux-carrying and energy-containing eddies in the ocean surface boundary layer. A new technique was developed to separate wave and turbulent motions spectrally, using ideas for turbulence spectra that were developed in the study of the bottom boundary layer of the atmosphere. The observations of turbulent fluxes allowed the closing of heat and momentum budgets across the air-sea interface. The observations also show that flux-carrying eddies are similar in size to those expected in rigid-boundary turbulence, but that energy-containing eddies are smaller than those in rigid-boundary turbulence. This suggests that the relationship between turbulent kinetic energy, depth, and turbulent diffusivity are different in the ocean surface boundary layer than in rigid-boundary turbulence. The observations confirm previous speculation that surface wave breaking provides a surface source of turbulent kinetic energy that is transported to depth where it dissipates. A model that includes the effects of shear production, wave breaking and dissipation is able to reproduce the enhancement of turbulent kinetic energy near the wavy ocean surface. However, because of the different length scale relations in the ocean surface boundary layer, the empirical constants in the energy model are different from the values that are used to model rigid-boundary turbulence. The ocean surface boundary layer is observed to have small but finite temperature gradients that are related to the boundary fluxes of heat and momentum, as assumed by closure models. However, the turbulent diffusivity of heat in the surface boundary layer is larger than predicted by rigid-boundary closure models. Including the combined effects of wave breaking, stress, and buoyancy forcing allows a closure model to predict the turbulent diffusivity for heat in the ocean surface boundary layer. / by Gregory Peter Gerbi. / Ph.D.

Dans les océans profonds linéairement stratifiés, la déstabilisation des ondes de gravité internes est importante car elle contribue probablement au mélange turbulent et à la circulation thermohaline.À l'aide de simulations numériques directes, nous créons un faisceau d'onde interne progressive. Cette situation est équivalente à une onde produite par l'oscillation de la marée sur une topographie sous-marine. Nous retrouvons les résultats expérimentaux obtenus par cite{Bourget13} : le faisceau se déstabilise en un mode petite échelle. Nous regardons l'effet d'un écoulement horizontal moyen sur cette instabilité en prenant soin d'abaisser la fréquence de forçage afin de compenser l'effet doppler et de conserver localement la même onde. Un cas limite apparaît lorsque le forçage devient stationnaire, ce qui équivaut à une onde de sillage issue d'un écoulement constant au dessus d'une topographie.Les écoulements à petite vitesse voient une instabilité petite échelle similaire au cas marée alors que les écoulement intermédiaires restent stables. Les écoulements plus rapides (jusqu'au cas sillage) voient, par contre, une instabilité bien plus grande échelle que celle dans le cas marée. Cette sélection d'échelle est robuste aux variations du nombre de Froude, de Reynolds, de la taille du faisceau ou de l'angle de l'onde.Nous montrons que ces instabilités peuvent être décrites comme des triades résonantes et que les différentes échelles correspondent à différentes branches triadiques. Nous confirmons la présence de cas stables pour des vitesses intermédiaires en calculant les modes propres comme des modes de Floquet à l'aide d'un algorithme d'Arnoldi--Krylov, et en montrant qu'ils sont associés à des taux de croissance négatifs.Le cas sillage est instable et nous le stabilisons par une méthode deselective frequency damping cite{Akervik06} afin d'obtenir un écoulement de base stationnaire autour duquel nous calculons les perturbations optimales qui maximisent l'énergie totale à différents horizons temporels. Pour des horizons courts, la perturbation optimale est petite échelle alors que pour des horizons longs, elle est grande échelle et converge vers la solution non-linéaire obtenue précédemment. Les horizons courts voient une instabilité triadique petite échelle advectée par l'écoulement et les horizons longs développent une instabilité d'une branche triadique grande échelle capable de se maintenir dans le faisceau malgré l'écoulement.Nous interprétons cette sélection de mode par le biais de la théorie des instabilités absolue ou convective. Dans le cas de l'onde de sillage l'instabilité grande échelle est absolue alors que la petite échelle est convective (et domine la croissance transitoire puisque son taux de croissance local est supérieur). Les rôles s'inversent dans le cas marée et l'instabilité petit échelle devient absolue alors que la grande échelle est convective. Nous confirmons cette hypothèse en calculant la réponse impulsionnelle d'une onde plane monochromatique dans un domaine 2Dpériodique. L'évolution spatio-temporelle d'une perturbation localisée en temps et en espace montre la formation de trois paquets d'onde, chacun étant associé à une branche triadique que nous identifions par une extension de la théorie triadique prenant en compte un désaccordage cite{McEwan77} et permettant de calculer la vitesse de groupe des sommets des paquets. En calculant ensuite le taux de croissance absolu le long de rayons à x/t et z/t constant, nous validons notre hypothèse. / Internal gravity waves that exist in a continuously stratified fluid are particularly important in the ocean. They transport energy and are thought to generate turbulent mixing, which contribute to the deep ocean circulation.We generate an internal wave beam that propagates in a continuously stratified fluid with direct numerical simulations. This situation is equivalent to a tidal wave, where the tidal flow oscillates over a topography and generates a wave. Experimental results obtained by cite{Bourget13} are recovered, ie. the beam destabilizes into a small scale mode. We consider the effect of an horizontal mean flow on the instability and lower the forcing frequency in order to compensate for the doppler effect and to keep locally the same wave. A limit case appears when the forcing becomes stationary. This case is equivalent to a lee wave appearing when a stratified fluid flows over a topography.For small mean flow, small scale instabilities develop as in the tidal case. The beam then stabilizes at intermediate mean flows and destabilizes again for increasing flow speed. At this second threshold, down to the lee wave case, the instability is of much larger scale than for the tidal case. Varying the Reynolds number, the Froude number, the wave angle or the beam size doesn't affect the instability scale selection : a small scale instability in the tidal regime, and large scale instability in the lee regime.We show that the instability mechanism may be interpreted using the triadic instability. Scale selection corresponds to different branches of triadic resonance. We confirm the presence of a stability region for intermediate value of the mean advection velocity by computing the linear eigenmode as Floquet mode with an Arnoldi-Krylov technique and show that the leading eigenmode has a negative growth rate.In the lee wave, case the flow is unstable and a selective frequency damping method cite{Akervik06} is used to compute a steady base flow. We then implement a linear direct-adjoint method to compute the optimal perturbations that maximizes the total energy at different time horizons. At short time horizon, the optimal perturbation is small scale while at large time the perturbation switches to a large scale solution and converges to the large scale mode observed through the nonlinear simulations. Short time transients correspond to the small scale triadic instability advected by the flow whereas the long time large scale instability corresponds to large scale branch of the triadic instability that is able to sustain the flow.We propose an interpretation of the selection of these different instabilities in term of absolute and convective instability. In the case of the lee wave, the large scale instability is absolute whereas the small scale instability is convective (and dominates the short time transient growth because it has a larger local growth rate). When the mean flow is varied, the properties of small scale and large scale instabilities exchange: in the tidal case the short scale instability is absolute and the large scale convective. This conjecture is confirmed by computing the impulse response around a plane monochromatic internal gravity wave in an extended two dimensional periodic domain. The spatio temporal evolution of a perturbation localized in space and time points out the formation of three different wave packets corresponding to different branches of triadic instability. Using the triadic theory with finite detuning cite{McEwan77},we derive the group velocity at the maximum growth rate of the three different branches of triadic instability and find a good agreement with the velocity of the three wave paquet maxima in the impulse response. Analyzing the impulse response along rays, i.e. at x/t and z/tconstant, we compute the absolute growth rate along all possible rays and validate our conjecture.

Thesis (Ph. D.)--Joint Program in Physical Oceanography (Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2013. / Cataloged from PDF version of thesis. / Includes bibliographical references (p. 187-195). / Large-scale thermal forcing and freshwater fluxes play an essential role in setting temperature and salinity in the ocean. A number of recent estimates of the global oceanic freshwater balance as well as the global oceanic surface net heat flux are used to investigate the effects of heat- and freshwater forcing at the ocean surface. Such forcing induces changes in both density and density-compensated temperature and salinity changes ('spice'). The ratio of the relative contributions of haline and thermal forcing in the mixed layer is maintained by large-scale surface fluxes, leading to important consequences for mixing in the ocean interior. In a stratified ocean, mixing processes can be either along lines of constant density (isopycnal) or across those lines (diapycnal). The contribution of these processes to the total mixing rate in the ocean can be estimated from the large-scale forcing by evaluating the production of thermal variance, salinity variance and temperature-salinity covariance. Here, I use new estimates of surface fluxes to evaluate these terms and combine them to generate estimates of the production of density and spice variance under the assumption of a linear equation of state. As a consequence, it is possible to estimate the relative importance of isopycnal and diapycnal mixing in the ocean. While isopycnal and diapycnal processes occur on very different length scales, I find that the surface-driven production of density and spice variance requires an approximate equipartition between isopycnal and diapycnal mixing in the ocean interior. In addition, consideration of the full nonlinear equation of state reveals that surface fluxes require an apparent buoyancy gain (expansion) of the ocean, which allows an estimate of the amount of contraction on mixing due to cabbeling in the ocean interior. / by Julian J. Schanze. / Ph.D.