ALBERT

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

Description: Ocean surface transport is at the core of many environmental disasters, including the spread of marine plastic pollution, the Deepwater Horizon oil spill and the Fukushima nuclear contamination. Understanding and predicting flow transport, however, remains a scientific challenge, because it operates on multiple length- and time-scales that are set by the underlying dynamics. Building on the recent emergence of Lagrangian methods, this thesis investigates the present-day abilities to describe and understand the organization of flow transport at the ocean surface, including the abilities to detect the underlying key structures, the regions of stirring and regions of coherence within the flow. Over the past four years, the field of dynamical system theory has adapted several algorithms from unsupervised machine learning for the detection of Lagrangian Coherent Structures (LCS). The robustness and applicability of these tools is yet to be proven, especially for geophysical flows. An updated, parameter-free spectral clustering approach is developed and a noise-based cluster coherence metric is proposed to evaluate the resulting clusters. The method is tested against benchmarks flows of dynamical system theory: the quasi-periodic Bickley jet, the Duffng oscillator and a modified, asymmetric Duffing oscillator. The applicability of this newly developed spectral clustering method, along with several common LCS approaches, such as the Finite-Time Lyapunov Exponent, is tested in several field studies. The focus is on the ability to predict these LCS in submesoscale ocean surface flows, given all the uncertainties of the modeled and observed velocity fields, as well as the sparsity of Lagrangian data. This includes the design and execution of field experiments targeting LCS from predictive models and their subsequent Lagrangian analysis. These experiments took place in Scott Reef, an atoll system in Western Australia, and off the coast of Martha's Vineyard, Massachusetts, two case studies with tidally-driven channel flows. The FTLE and spectral clustering analyses were particularly helpful in describing key transient flow features and how they were impacted by tidal forcing and vertical velocities. This could not have been identified from the Eulerian perspective, showing the utility of the Lagrangian approach in understanding the organization of transport.

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, 1982

Description: Oceanic fluctuations are dependent on geographical location. Near intense currents, the eddy field is highly energetic and has broad meridional extent. It is likely that the energy arises from instabilities of the intense current. However, the meridional extent of the linearly most unstable modes of such intense jets is much narrower than the observed region of energetic fluctuations. It is proposed here that weaker instabilities, in the linear sense, which are very weakly trapped to the current, may be the dominant waves in the far field. As a preliminary problem, the (barotropic) instability of parallel shear flow on the beta plane is discussed. An infinite zonal flow with a continuous cross-stream velocity gradient is approximated with segments of uniform flow, joined together by segments of uniform potential vorticity. This simplification allows an exact dispersion relation to be found. There are two classes of linearly unstable solutions. One type is trapped to the source of energy and has large growth rates. The second type are weaker instabilities of the shear flow which excite Rossby waves in the far field: the influence of these weaker instabilities extends far beyond that of the most unstable waves. The central focus of the thesis i: the linear stability of thin, twolayer, zonal jets on the beta plane, with both horizontal and vertical shear. The method used for the parallel shear flow is extended to the two-layer flow. Each layer of the jet has uniform velocity in the center, bordered by shear zones with zero potential vorticity gradient. The velocity in each layer outside the jet is constant in latitude. Separate linearly unstable modes arise from horizontal and vertical shear. The energy source for the vertical shear modes is nearly all potential while the source for the horizontal shear modes is both kinetic and potential. The most unstable waves are tightly trapped to the jet, within two or three deformation radii for small but nonzero beta. Rossby waves and baroclinically unstable waves (in the presence of vertical shear) exist outside the jet because of a nonzero potential vorticity gradient there. Weakly growing jet instabilities can force these waves when their phase speeds and wavelengths match. In particular, westward jets and any jets with vertical shear exterior to the jet can radiate in this sense. The radiating modes influence a large region, their decay scales inversely proportional to the growth rate. Two types of radiating instability are found: (1) a subset of the main unstable modes near marginal stability and (2) modes which appear to be destabilized neutral modes. Westward jets have more vigorously unstable radiating modes. Applications of the model are made to the eddy field south of the Gulf Stream, using data from the POLYMODE settings along 55°W and farther into the gyre at MODE. The energy decay scale and the variation of vertical structure with latitude in different frequency bands can be roughly explained by the model. The lower frequency disturbances decay more slowly and become more surface intensified in the far field. These disturbances are identified with the weak, radiating instabilities of the model. The higher frequency disturbances are more trapped and retain their vertical structure as they decay, and are identified with the trapped, strongly unstable modes of the jet.

Description: This work was supported by a grant from the National Science Foundation, Office of Atmospheric Science.

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 2017

Description: This thesis investigates the evolution of the oceanic lithosphere in a broad sense from formation to subduction, in a focused case at the ridge, and in a focused case proximal to subduction. In general, alteration of the oceanic lithosphere begins at the ridge through focused and diffuse hydrothermal flow, continues off axis through low temperature circulation, and may occur approaching subduction zones as bending related faulting provides fluid pathways. In Chapter 2 I use a dataset of thousands of microearthquakes recorded at the Rainbow massif on the Mid-Atlantic Ridge to characterize the processes which are responsible for the long-term, high-temperature, hydrothermal discharge found hosted in this oceanic core complex. I find that the detachment fault responsible for the uplift of the massif is inactive and that the axial valleys show no evidence for faulting or active magma intrusion. I conclude that the continuous, low-magnitude seismicity located in diffuse pattern in a region with seismic velocities indicating ultramafic host rock suggests that serpentinization may play a role in microearthquake generation but the seismic network was not capable of providing robust focal mechanism solutions to constrain the source characteristics. In Chapter 3 I find that the Juan de Fuca plate, which represents the young/hot end-member of oceanic plates, is lightly hydrated at upper crustal levels except in regions affected by propagator wakes where hydration of lower crust and upper mantle is evident. I conclude that at the subduction zone the plate is nearly dry at upper mantle levels with the majority of water contained in the crust. Finally, in Chapter 4 I examine samples of cretaceous age serpentinite sampled just before subduction at the Puerto Rico Trench. I show that these upper mantle rocks were completely serpentinized under static conditions at the Mid-Atlantic Ridge. Further, they subsequently underwent 100 Ma of seafloor weathering wherein the alteration products of serpentinization themselves continue to be altered. I conclude that complete hydration of the upper mantle is not the end point in the evolution of oceanic lithosphere as it spreads from the axis to subduction.

Description: Funding was provided by the National Science Foundation through grants OCE-1029305 and OCE-0961680, the Deep Ocean Exploration Institute - Ocean Ridge Initiative, and by the WHOI Academic Programs office

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 2017

Description: The daily heating of the ocean by the sun can create a stably stratified near-surface layer when the winds are slight and solar insolation is strong. This type of shallow stable layer is called a Diurnal Warm Layer (DWL). This thesis examines the physics and dynamics of DWLs from observations of the subtropical North Atlantic Ocean associated with the Salinity Processes in the Upper ocean Regional Study (SPURS-I). Momentum transferred from the atmosphere to the ocean through wind stress becomes trapped within the DWL, generating shear across the layer. During SPURS-I, strong diurnal shear across the DWL was coincident with enhanced turbulent kinetic energy (TKE) dissipation (𝜖, 𝜖 〉 10−5 W/kg) observed from glider microstructure profiles of the near-surface. However, a scale analysis demonstrated that surface forcing, including diurnal shear, could not be the sole mechanism for the enhanced TKE dissipation. High-frequency internal waves (𝜔 ≫ 𝑓) were observed in the upper ocean during the daytime within the DWL. Internal waves are able to transfer energy from the deep ocean into the DWL through the unstratified remnant mixed layer, which is the intervening layer between the DWL and seasonal thermocline. As the strength of the stratification of the DWL increases, so does the shear caused by the tunneling internal waves. The analysis demonstrates that internal waves can generate strong enough shear to cause a shear-induced instability, and are a plausible source of the observed enhanced TKE dissipation. Vertically-varying horizontal transport across the upper ocean occurs because a diurnal current exists within the DWL, but not in the unstratified remnant mixed layer below. Therefore, when a DWL is present, the water within DWL is horizontally transported a different distance than the water below. Coupled with nocturnal convection that mixes the DWL with the unstratified layer at night, this cycle is a mechanism for submesoscale (1-10 km) lateral diffusion across the upper ocean. Estimates of a horizontal diffusion coefficient are similar in magnitude to current estimates of submesoscale diffusion based on observations, and are likely an important source of horizontal diffusion in the upper ocean.

Description: Supported by the Department of Defense (DoD) through the National Defense Science & Engineering Graduate Fellowship (NDSEG) Program and the National Science Foundation under Grant No. OCE-1129646. The collection and analysis of data from the SPURS-I central mooring were supported under National Aeronautics and Space Administration (NASA) Grant No. NNX11AE84G and NNX14AH38G.

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 2016

Description: Since the Last Glacial Maximum (LGM, ~ 20,000 years ago) air temperatures warmed, sea level rose roughly 130 meters, and atmospheric concentrations of carbon dioxide increased. This thesis combines global models and paleoceanographic observations to constrain the ocean’s role in storing and transporting heat, salt, and other tracers during this time, with implications for understanding how the modern ocean works and how it might change in the future. • By combining a kinematic ocean model with “upstream” and “downstream” deglacial oxygen isotope time series from benthic and planktonic foraminifera, I show that the data are in agreement with the modern circulation, quantify their power to infer circulation changes, and propose new data locations. • An ocean general circulation model (the MITgcm) constrained to fit LGM sea surface temperature proxy observations reveals colder ocean temperatures, greater sea ice extent, and changes in ocean mixed layer depth, and suggests that some features in the data are not robust. • A sensitivity analysis in the MITgcm demonstrates that changes in winds or in ocean turbulent transport can explain the hypothesis that the boundary between deep Atlantic waters originating from Northern and Southern Hemispheres was shallower at the LGM than it is today.

Description: Support for this work came from an MIT Presidential Fellowship, an NSF Graduate Research Fellowship, and grants NASA NNX12AJ93G – Gravity data for ocean circulation and climate studies, NSF OCE-0961713 – Collaborative Research: The Physics and Statistics of Global Sea Level Change, NSF OCE-1060735 – Collaborative Research: Beyond the Instrumental Record - the Ocean Circulation at the last Glacial maximum and the deglacial sequence, and NASA NNX08AR33G – Application of Satellite Altimetry Gravity Winds and in Situ Data to Problems of the Ocean Circulation.