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 2021.

Description: This thesis explores the volatile content of the mantle, subducted oceanic crust, and arc magmas as well as the structure of slow spreading ocean crust and the heterogeneity of Earth’s upper mantle. In Chapter 2, I directly explore the halogen (F and Cl) content of mantle minerals in situ, then use these measurements to assess the halogen content of the upper mantle. In Chapter 3, I investigate the volatile content of Raspas eclogites (SW Ecuador), a proxy for deeply subducted oceanic crust, to evaluate volatile transfer from crustal generation at divergent plate boundaries (e.g., mid-ocean ridges) to recycling of ocean crust at subduction zones. In Chapter 4, I use the H2O content of nominally anhydrous minerals in plutonic arc cumulates to elucidate the H2O content of the melts from which the rocks crystallized. In this way, I assert that primitive arc magmas may contain 4–10 wt.% H2O and through fractional crystallization up to ~20 wt.% H2O, making them far more hydrous than traditional methods (i.e., olivine-hosted melt inclusions) surmise. In Chapter 5, I show that mantle peridotite exposed along the 16ºN region of the Mid-Atlantic Ridge originated in an arc setting and has been remixed into subridge mantle, indicating that the sub-ridge mantle is more heterogeneous and depleted than inferences made from mid-ocean ridge basalts suggest. Chapter 6 surveys the life cycle of oceanic core complexes through zircon geochronology and posits an updated framework for understanding the termination of oceanic core complexes, and more broadly oceanic detachment faults. Together, this contribution highlights the chemical heterogeneity of the mantle, and quantifies the full extent of volatiles hosted by mantle and crustal reservoirs.

Description: The Stanley Watson Fellowship (WHOI) provided financial support during my first year of graduate school. The Academic Programs Office Ocean Venture Fund (WHOI) provided seed funding which initiated Chapters 3 and 4, and ultimately led to two funded NSF proposals. These resources are vital to JP students, and I am incredibly grateful for them. Primary support was provided by the National Science Foundation grants to Veronique Le Roux (EAR P&G #1524311, #1839128, #1855302) and Henry Dick (MG&G #1637130, #1657983).

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 2021.

Description: Oceanic fronts at the mesoscale and submesoscale are associated with enhanced vertical motion, which strengthens their role in global biogeochemical cycling as hotspots of primary production and subduction of carbon from the surface to the interior. Using process study models, theory, and field observations of biogeochemical tracers, this thesis improves understanding of submesoscale vertical tracer fluxes and their influence on carbon cycling. Unlike buoyancy, vertical transport of biogeochemical tracers can occur both due to the movement of isopycnals and due to motion along sloping isopycnals. We decompose the vertical velocity below the mixed layer into two components in a Lagrangian frame: vertical velocity along sloping isopycnal surfaces and the adiabatic vertical velocity of isopycnal surfaces and demonstrate that vertical motion along isopycnal surfaces is particularly important at submesoscales (1-10 km). The vertical flux of nutrient, and consequently the new production of phytoplankton depends not just on the vertical velocity but on the relative time scales of vertical transport and nutrient uptake. Vertical nutrient flux is maximum when the biological timescale of phytoplankton growth matches the vertical velocity frequency. Export of organic matter from the surface and the interior requires water parcels to cross the mixed layer base. Using Lagrangian analysis, we study the dynamics of this process and demonstrate that geostrophic and ageostrophic frontogenesis drive subduction along density surfaces across the mixed layer base. Along-front variability is an important factor in subduction. Both the physical and biological modeling studies described above are used to interpret observations from three research cruises in the Western Mediterranean. We sample intrusions of high chlorophyll and particulate organic carbon below the euphotic zone that are advected downward by 100 meters on timescales of days to weeks. We characterize the community composition in these subsurface intrusions at a lateral resolution of 1–10 km. We observe systematic changes in community composition due to the changing light environment and differential decay of the phytoplankton communities in low-light environments, along with mixing. We conclude that advective fluxes could make a contribution to carbon export in subtropical gyres that is equal to the sinking flux.

Description: The work in this dissertation was funded by a NDSEG fellowship, Martin Fellowship, Grassle fellowship, Montrym grant, WHOI Academic Programs Office, and Office of Naval Research CALYPSO DRI grant N00014-16-1-3130.

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science in Mechanical Engineering at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 2021.

Description: Operations in the Arctic Ocean are increasingly important due to the changing environment and the resulting global implications. These changes range from the availability of new global trade routes, accessibility of newly available resources in the area, and national security interests of the United States in the region. It’s necessary to build a greater understanding of the undersea environment and how it’s changing since these environmental changes have a direct impact on adjusting future operations in the region and looming global changes as less Arctic ice is present. The recent presence of the Beaufort Lens is changing the acoustic propagation paths throughout the Arctic region. Here a network of buoys were employed to communicate with an Autonomous Undersea Vehicle (AUV) while it operated under the ice throughout the Beaufort Lens with the goal of achieving near GPS quality navigation. The acoustic communications paths were compared using a vertical array throughout the Beaufort Lens. This beam forming was compared to the prediction from BELLHOP. As well, since acoustic communications are affected by multi-path, attenuation and interference from other sources it was interesting to note that bottom bounce was sometimes a reliable acoustic path.

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 2021.

Description: Chromium (Cr) isotopes have shown great potential as a paleo-redox proxy to trace the redox conditions of ancient oceans and atmosphere. However, its cycling in modern environments is poorly constrained. In my thesis, I attempt to fill in the gap of our understanding of chromium cycling in the modern ocean, with a focus on the redox processes in global oxygen deficient zones (ODZs). Firstly, we developed a method to analyze Cr isotopes of different Cr redox species. Tests on processing conditions demonstrated its robustness in obtaining accurate Cr isotope data. It is applicable to both frozen and fresh samples. This method allows us to investigate the redox cycling of Cr that is hard to unravel by existing total Cr methods. Secondly, in the Eastern Tropical North Pacific (ETNP), Eastern Tropical South Pacific (ETSP) and Arabian Sea ODZs, their total dissolved Cr profiles show preferential reduction of isotopically light Cr(VI) to Cr(III), which is scavenged and exported to deeper oceans. Applying our new method to ETNP and ETSP ODZ seawater samples, we observed Cr(VI) reduction in both ODZs with a similar fractionation factor. This indicates similar mechanisms may be controlling Cr(VI) reduction in the two ODZs. Cr(III) maximum coincides with Fe(II) and secondary nitrite maximums in the upper core of both ODZs. Shipboard incubations with spiked Fe(II) showed fast Cr(VI) reduction occurring in the ETNP ODZ. But neither Fe(II) nor microbes were reducing Cr(VI) directly. Thirdly, we calculated the isotope effects of Cr scavenging in the ETNP and ETSP ODZs. Thetwo ODZs show a similar isotope partitioning during Cr scavenging. And spatial variability is observed in the ETNP ODZ. Our calculated scavenged Cr isotope ratio is lighter than that of the total dissolved Cr from the same depth. It is also comparable to that of reducing or anoxic sediments, which implies that Cr isotopes can be used as an archive for local redox conditions.

Description: This research was supported by an anonymous MIT Fellowship, Praceis Presidential Fellowship, Frederick A. Middleton Fellowship, the US National Science Foundation (NSF Award No. OCE-1459287, OCE-1736996, OCE-1924050 and DEB-1542240) and the Center for Microbial Oceanography: Research and Education (C-MORE, NSF-OIA Award No. EF-0424599).

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 2022.

Description: Detection, classification, localization, and tracking (DCLT) of unmanned underwater vehicles (UUVs) in the presence of shipping traffic is a critical task for passive acoustic harbor security systems. In general, vessels can be tracked by their unique acoustic signature due to machinery vibration and cavitation noise. However, cavitation noise of UUVs is considerably quieter than ships and boats, making detection significantly more challenging. In this thesis, I demonstrated that it is possible to passively track a UUV from its highfrequency motor noise using a stationary array in shallow-water experiments with passing boats. First, causes of high frequency tones were determined through direct measurements of two UUVs at a range of speeds. From this analysis, common and dominant features of noise were established: strong tones at the motor’s pulse-width modulated frequency and its harmonics. From the unique acoustic signature of the motor, I derived a high-precision, remote sensing method for estimating propeller rotation rate. In shallow-water UUV field experiments, I demonstrated that detecting a UUV from motor noise, in comparison to broadband noise from the vehicle, reduces false alarms from 45% to 8.4% for 90% true detections. Beamforming on the motor noise, in comparison to broadband noise, improved the bearing accuracy by a factor of 3.2×. Because the signal is also high-frequency, the Doppler effect on motor noise is observable and I demonstrate that range rate can be measured. Furthermore, measuring motor noise was a superior method to the “detection of envelope modulation on noise” algorithm for estimating the propeller rotation rate. Extrapolating multiple measurements from the motor signature is significant because Bearing-Doppler-RPM measurements outperform traditional bearing-Doppler target motion analysis. In the unscented Kalman filter implementation, the tracking solution accuracy for bearing, bearing rate, range, and range rate improved by a factor 2.2×, 15.8×, 3.1×, and 6.2× respectively. These findings are significant for improving UUV localization and tracking, and for informing the next-generation of quiet UUV propulsion systems.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Chemical Oceanography at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution June 2021.

Description: The existence of a marine phosphorus (P) redox cycle was recently confirmed when phosphonates, organophosphorus compounds with P in the (III) oxidation state, were found in high molecular weight dissolved organic matter. Although some features of the P redox cycle have come to light since the discovery of phosphonates, many aspects of phosphonate production, cycling and fate remain unknown. To address these gaps in our understanding, we studied phosphonate cycling in the Eastern Mediterranean Sea, a chronically P-limited basin, using 33P and enzymatic assays. We showed that phosphonate production was low but consumption was high, suggesting that phosphonate production and consumption may be spatially or temporally decoupled. We also explored phosphonate production in the model marine cyanobacterium Prochlorococcus SB. Using 31P NMR, we found Prochlorococcus SB allocates ~50% of its cellular P to phosphonates. Allocation of P to phosphonates was conserved under P-limitation, and further investigation revealed phosphonates were associated with proteins. The discovery of phosphonoproteins in Prochlorococcus SB opens new perspectives on the biochemical function of phosphonates and their role in P-cycling. Finally, we developed a new P-targeted method to characterize marine organophosphorus compounds using liquid chromatography coupled to electrospray ionization and inductively coupled plasma mass spectrometry.

Description: This work was supported by the Simons Foundation under grant numbers POP49476 and 721227 [D. Repeta], the Gordon and Betty Moore Foundation under the grant number 6000 [D. Repeta] and the National Science Foundation OCE under the grant number 1634080 [D. Repeta].

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Geophysics at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 2021.

Description: This thesis focuses on interpreting geophysical and geochemical observables in terms of the thermomechanical state of the lithosphere. In Chapter 1, I correlate lower crustal rheology with seismic wave speed. Compositional variation is required to explain half of the total variability in predicted lower crustal stress, implying that constraining regional lithology is important for lower crustal geodynamics. In Chapter 2, I utilize thermobarometry, diffusion models, and thermodynamic modelling to constrain the ultra-high formation conditions and cooling rates of the Gore Mountain Garnet Amphibolite in order to understand the rheology of the lower crust during orogenic collapse. In Chapter 3, I interpret geophysical data along a 74 Myr transect in the Atlantic to the temporal variability and relationship of crustal thickness and normal faults. In Chapter 4, I constrain the error present in the forward-calculation of seismic wave speed from ultramafic bulk composition. I also present a database and toolbox to interpret seismic wave speeds in terms of temperature and composition. Finally, in Chapter 5 I apply the methodology from Chapter 4 to interpret a new seismic tomographic model in terms of temperature, density, and composition in order to show that the shallow lithospheric roots are density unstable.

Description: Funding for this research was provided by an MIT Presidential Fellowship, MIT Student Research Funds, the National Science Foundation Division of Earth Sciences (EAR) and Ocean Sciences (OCE) grants EAR-16-24109, EAR-17-22932, EAR-17-22935, OCE-14-58201, and SCEC Awards 16106 and 17202., SCEC, Geological Society of America Graduate Student Research Fellowship, WHOI Ocean Venture Fund, and 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 2021.

Description: Observations of hydrographic and dynamical properties on the Middle Atlantic Bight shelf document strong variability at time scales spanning events that last a few days to century long trends. This thesis studies individual processes which impact shelf temperature and velocity structure, and quantifies the mean velocity conditions at the shelf break. Chapter 2 uses model output to study the dynamics that lead to the breakdown of summertime thermal stratification, and how the processes which reduce stratification vary from year to year. In summer, the atmosphere heats the surface of the ocean, leading to strong thermal stratification with warm water overlying cool water. During fall, strong storm events with downwelling-favorable winds are found to be the primary process by which stratification is reduced. The timing of these events and the associated destratification varies from year to year. In Chapter 3, the velocity structure of the New Jersey shelf break is examine, with a focus on the Shelfbreak Jet. Using 25 years of velocity measurements, mean velocity sections of the Shelfbreak Jet are created in both Eulerian and stream coordinate frameworks. The jet exhibits strong seasonal variability, with maximum velocities observed in spring and minimum velocities in summer. Evidence is found that Warm Core Rings, originating from the Gulf Stream and passing through the Slope Sea adjacent to the New Jersey shelf, tend to shift the Shelfbreak Jet onshore of its mean position or entirely shutdown the Shelfbreak Jet’s flow. At interannual timescales, variability in the Shelfbreak Jet velocity is correlated with the temperature on the New Jersey Shelf, with temperature lagging by about 2 months. Chapter 4 focuses on the impact of Warm Core Rings on the velocity and temperature structure on the New Jersey shelf. Warm Core Rings that have higher azimuthal velocities and whose cores approach closer to the shelf are found to exert greater influence on the shelf’s along-shelf velocities, with the fastest and closest rings reversing the direction of flow at the shelf break. Warm Core Rings are also observed to exert long-lasting impacts on the shelf temperature, with faster rings cooling the shelf and slower rings warming the shelf. Seasonal changes in thermal stratification strongly affect how rings alter the shelf temperature. Rings in summer tend to cool the shelf, and rings throughout the rest of the year generally warm the shelf.

Description: This research was funded under WHOI Academic Programs Endowed Funds, NSF #OCE-1634094, and NSF #OCE-1924041.

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 2021.

Description: Small pelagic fishes, also termed forage fishes, represent a critical link between secondary production and myriad top predators in marine ecosystems, including the Northeast US shelf. In this dissertation, I analyze the drivers of forage fish distribution throughout the Northeast US shelf and the drivers of the abundance of the ecologically important northern sand lance. Chapter 2 examines the basic ecology of northern sand lance and uses these insights to identify mechanistic drivers of their abundance. I then explore different scenarios of these drivers to project sand lance abundance through the end of the 21st century, which appears precarious for adult sand lance unless current trajectories change. Chapter 3 analyzes the environmental drivers of the distribution of the six dominant, offshore forage fish species (northern sand lance, Atlantic herring, alewife, blueback herring, Atlantic mackerel, and Atlantic butterfish) on the Northeast US shelf to elucidate the role of environmental covariates in shelf occupancy by these taxa. The results of this chapter indicate shelf occupancy of butterfish and Atlantic mackerel are increasing through time while occupancy of sand lance is decreasing with time. The occurrence of most of these species is also moving deeper and northward with time. Chapter 4 assesses the source-sink dynamics of three sand lance hotspots through Lagrangian particle tracking models simulating larval sand lance transport. Connectivity varies among these hotspots with Georges Bank and Stellwagen Bank having notable retention while the Great South Channel relies on larvae from other hotspots. Retention on Stellwagen Bank and Georges Bank are linked to strong wind events during the larval period of sand lance. Collectively, this dissertation improves our understanding of the dynamics driving variability in the Northeast US shelf forage fish complex, particularly for northern sand lance.

Description: The research within this dissertation was funded by a National Science Foundation Graduate Research Fellowship (awarded to JJS), Woods Hole Sea Grant (NA18OAR4170104, Project No. R/O-57), the Bureau of Ocean Energy Management (IA agreement M17PG0019), and the National Marine Sanctuary 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 2021.

Description: An emerging paradigm posits that the abyssal overturning circulation is driven by bottom-enhanced mixing, which results in vigorous upwelling in the bottom boundary layer (BBL) along the sloping seafloor and downwelling in the stratified mixing layer (SML) above; their residual is the overturning circulation. This boundary-controlled circulation fundamentally alters abyssal tracer distributions, with implications for global climate. Chapter 1 describes how a basin-scale overturning circulation arises from the coupling between the ocean interior and mixing-driven boundary layers over rough topography, such as the sloping flanks of mid-ocean ridges. BBL upwelling is well predicted by boundary layer theory, whereas the compensation by SML downwelling is weakened by the upward increase of the basin-wide stratification, which supports a finite net overturning. These simulated watermass transformations are comparable to best-estimate diagnostics but are sustained by a crude parameterization of boundary layer restratification processes. In Chapter 2, I run a realistic simulation of a fracture zone canyon in the Brazil Basin to decipher the non-linear dynamics of abyssal mixing layers and their interactions with rough topography. Using a hierarchy of progressively idealized simulations, I identify three physical processes that set the stratification of abyssal mixing layers (in addition to the weak buoyancy-driven cross-slope circulation): submesoscale baroclinic eddies on the ridge flanks, enhanced up-canyon flow due to inhibition of the cross-canyon thermal wind, and homogenization of canyon troughs below the level of blocking sills. Combined, these processes maintain a sufficiently large near-boundary stratification for mixing to drive globally significant BBL upwelling. In Chapter 3, simulated Tracer Release Experiments illustrate how passive tracers are mixed, stirred, and advected in abyssal mixing layers. Exact diagnostics reveal that while a tracer’s diapycnal motion is directly proportional to the mean divergence of mixing rates, its diapycnal spreading depends on both the mean mixing rate and an additional non-linear stretching term. These simulations suggest that the theorized boundary-layer control on the abyssal circulation is falsifiable: downwelling in the SML has already been confirmed by the Brazil Basin Tracer Release Experiment, while an upcoming experiment in the Rockall Trough will confirm or deny the existence of upwelling in the BBL.

Description: This material is based upon work supported by the National Science Foundation Graduate Research Fellowship under Grant No. 174530. I also acknowledge funding support from National Science Foundation Awards OCE-1536515 and OCE-1736109. This work was partially supported by MIT’s Rosenblith Presidential Fellowship.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Biological Oceanography at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 2021.

Description: Sound is utilized by marine animal taxa for many ecologically important functions, and these taxa are vulnerable to adverse effects of anthropogenic noise on hearing and behavior. However, little is known about marine invertebrates’ responses to anthropogenic noise, and the ambient environmental sounds (“soundscapes”) they detect and respond to. Most acoustic studies report sound pressure (detected by mammals and some fish), but few report particle motion, the back-and-forth vibratory component of sound detected by marine invertebrates. I investigated invertebrate use of and response to sounds in two facets: 1) behavioral responses of longfin squid, Doryteuthis pealeii to anthropogenic noise, and 2) particle motion of coral reef soundscapes in the U.S. Virgin Islands. In laboratory-based experiments I exposed D. pealeii to construction noise originally recorded from an offshore wind farm. I found significant increases in squids’ alarm responses and in failed prey capture attempts during noise. Conversely, noise exposure had no significant effects on reproductive behaviors of groups of D. pealeii, indicating high motivation of these squid to reproduce during this stressor. Collectively, these experiments revealed the importance of considering behavioral context in studies and regulatory decisions regarding invertebrates’ susceptibility to anthropogenic noise impacts. In studying coral reef soundscapes, I reported particle motion trends over several months for coral reefs varying in habitat quality, including coral cover and fish abundance. I found acoustic properties over which particle motion closely scaled with pressure, and others over which it did not. I compared soundscape data with particle motion hearing thresholds, and found that invertebrates may only detect high amplitude and low frequency transient sound cues on reefs, such as those produced by fishes. My research bring new insights on natural and anthropogenic sound cues detectable by marine invertebrates, and how and when invertebrates will be vulnerable to anthropogenic noise pollution.

Description: My graduate work was funded in part by the US Department of Interior, Bureau of Ocean Energy Management Environmental Studies Program through Interagency Agreement Number M17PG00029 with the U.S. Department of Commerce, National Oceanic and Atmospheric Administration (funding to Aran Mooney and Jenni Stanley). My work was also supported by the NSF Biological Oceanography award OCE-1536782 (funding to Aran Mooney). I received tuition and stipend support from the National Science Foundation Graduate Research Fellowship Program [Grant No. 2388357]. The Academic Program Office at the Woods Hole Oceanographic Institution provided tuition and stipend support as well as travel support. The MIT Student Assistance Fund, the Aquatic Noise 2019 Organizing Committee, and the Acoustical Society of America also provided travel support.

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science in Mechanical Engineering at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 2021.

Description: A novel performance metric to improve underwater digital acoustic communication, called Multipath Penalty (MPP), is proposed as an alternative to traditional signal-to-noise ratio (SNR) methods in the context of the Arctic Beaufort Sea. MPP and SNR are compared alongside a third performance metric, Minimum Achievable Error (MAE), which replicates the operation of a channel estimate-based decision feedback equalizer in an acoustic modem. The three metrics are then tested in a hardware-in-the-loop Virtual Ocean simulator for an autonomous undersea vehicle (AUV) communicating with a collaborator. Using field data of modem statistics obtained duringICEX20 and expanded data supplied by the simulator, calibration of the three metrics to modem packet success is evaluated, resulting in a proposed recalibration for MAE. The AUV’s ability to communicate when adaptively choosing its depth is analyzed above and below the Beaufort Lens, and settings for MPP’s engineering variables are obtained. The results show MPP generally improves reception and demodulation of acoustic transmissions over SNR by approximately 5% within an operational range of 8 km, while achieving similar results to the more robust metric MAE. MPP is an improved utility for underwater digital acoustic communication in both marine autonomy and as a tactical decision aid.

Description: This work would not be possible without the extraordinary support of the United States Navy, which provided funding for this research, my degree, and my livelihood as an active duty submarine officer.

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science in Mechanical Engineering at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 2021.

Description: There are many significant challenges for unmanned autonomous platforms at sea including predicting the likely scenarios for the ocean environment, quantifying regional uncertainties, and updating forecasts of the evolving dynamics using their observations. Due to the operational constraints such as onboard power, memory, bandwidth, and space limitations, efficient adaptive reduced order models (ROMs) are needed for onboard predictions. In the first part, several reduced order modeling schemes for regional ocean forecasting onboard autonomous platforms at sea are described, investigated, and evaluated. We find that Dynamic Mode Decomposition (DMD), a data-driven dimensionality reduction algorithm, can be used for accurate predictions for short periods in ocean environments. We evaluate DMD methods for ocean PE simulations by comparing and testing several schemes including domain splitting, adjusting training size, and utilizing 3D inputs. Three new approaches that combine uncertainty with DMD are also investigated and found to produce practical and accurate results, especially if we employ either an ensemble of DMD forecasts or the DMD of an ensemble of forecasts. We also demonstrate some results from projecting / compressing high-fidelity forecasts using schemes such as POD projection and K-SVD for sparse representation due to showing promise for distributing forecasts efficiently to remote vehicles. In the second part, we combine DMD methods with the GMM-DO filter to produce DMD forecasts with Bayesian data assimilation that can quickly and efficiently be computed onboard an autonomous platform. We compare the accuracy of our results to traditional DMD forecasts and DMD with Ensemble Kalman Filter (EnKF) forecast results and show that in Root Mean Square Error (RMSE) sense as well as error field sense, that the DMD with GMM-DO errors are smaller and the errors grow slower in time than the other mentioned schemes. We also showcase the DMD of the ensemble method with GMM-DO. We conclude that due to its accurate and computationally efficient results, it could be readily applied onboard autonomous platforms. Overall, our contributions developed and integrated stochastic DMD forecasts and efficient Bayesian GMM-DO updates of the DMD state and parameters, learning from the limited gappy observation data sets.

Description: I wish to thank the U.S. Navy’s Civilian Institution Program along with the MITWHOI Joint Program for providing the funding and resources that made this research and continuing my education possible.

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 September 2021.

Description: The Beaufort Gyre region of the Arctic Ocean is strongly stratified at the base of the wintertime mixed layer, which impedes the vertical transport of heat, energy, and other tracers. Ice-Tethered Profiler observations during 2004-2018 were used to characterize and investigate the seasonal and interannual variability of the strength, depth, density, and thickness of this highly stratified layer at the base of the mixed layer. This includes investigating the remnant stratification maximum, which formed when the summer mixed layer shoaled. Seasonally, the stratification maximum was never in a steady state. It was largest in October (4.8 × 10−3 rad2/sec2) and decreased during all winter months (to 2.3 × 10−3rad2/sec2 in June), indicating that surface forcing and interior vertical mixing were never in equilibrium during the year. Interannually, the period from 2011-2018 had a higher stratification maximum than then the period from 2005-2010 regardless of the season. The remnant stratification maximum was consistently weaker than the winter stratification maximum from which it formed. The initial evolution of the remnant stratification maximum is used to estimate an effective vertical diffusivity of order 10−6m2/s. No significant geographic variability was found, in part due to high temporal and small scale variability of the stratification maximum layer. Implications for heat transport through to the sea ice cover are discussed.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Biological Oceanography at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution February 2022.

Description: Protists are taxonomically and metabolically diverse drivers of energy and nutrient flow in the marine environment, with recent research suggesting significant roles in global carbon cycling throughout the water column. Top-down controls on planktonic protists include grazing and parasitism, processes that both contribute to nutrient transfer and biogeochemical cycling in the global ocean. Recent global surveys of eukaryotic small subunit ribosomal RNA molecular signatures have highlighted the fact that parasites belonging to the marine alveolate order Syndiniales are both abundant and ubiquitous in coastal and open ocean environments, suggesting a major role for this taxon in marine food webs. Two coastal sites, Saanich Inlet (Vancouver Island, BC) and Salt Pond (Falmouth, MA, USA) were selected as model ecosystems to examine the impacts of Syndinian parasitism on protist communities. Data presented in this thesis combines high-resolution sampling, water chemistry (including nutrients) analyses, molecular marker gene analyses, fluorescence in situ hybridization, and modeling to address key knowledge gaps regarding syndinian ecology. Information is presented on previously undescribed putative host taxa, the prevalence of syndinian parasites and infections on different hosts in coastal waters, and a framework for modeling host-parasite interactions based on field observations.

Description: Research was supported by the WHOI Ocean Venture Fund, the National Science Foundation Biological Oceanography OCE-1851012, and the National Science Foundation Graduate Research Fellowship under Grant No. 1745302.

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 2022.

Description: The ventilation of intermediate waters in the Labrador Sea has important implications for the strength of the Atlantic Meridional Overturning Circulation. Boundary current-interior interactions regulate the exchange of properties between the slope and the basin, which in turn regulates the magnitude of interior convection and the export of ventilated waters from the subpolar gyre. This thesis characterizes theWest Greenland Boundary Current System near Cape Farewell across a range of spatio-temporal scales. The boundary current system is composed of three velocity cores: (1) the West Greenland Coastal Current (WGCC), transporting Greenland and Arctic meltwaters on the shelf; (2) the West Greenland Current (WGC), which advects warm, saline Atlantic-origin water at depth, meltwaters at the surface, and newly-ventilated Labrador Sea Water (LSW); and (3) the Deep Western Boundary Current, which carries dense overflow waters ventilated in the Nordic Seas. The seasonal presence of the LSW and Atlantic-origin water are dictated by air-sea buoyancy forcing, while the seasonality of the WGCC is governed by remote wind forcing and the propagation of coastally trapped waves from East Greenland. Using mooring data and hydrographic surveys, we demonstrate mid-depth intensified cyclones generated at Denmark Strait are found offshore of the WGC and enhance the overflow water transport at synoptic timescales. Using mooring, hydrographic, and satellite data, we demonstrate that the WGC undergoes extensive meandering due to baroclinic instability that is enhanced in winter due to LSW formation adjacent to the current. This leads to the production of small-scale, anticyclonic eddies that can account for the entirety of wintertime heat loss within the Labrador Sea. The meanders are shown to trigger the formation of Irminger Rings downstream. Using mooring, hydrographic, atmospheric, and Lagrangian data, and a mixing model, we find that strong atmospheric storms known as forward tip jets cause upwelling at the shelfbreak that triggers offshore export of freshwater. This freshwater flux can explain the observed lack of ventilation in the eastern Labrador Sea. Together, this thesis documents previously unobserved interannual, seasonal, and synoptic-scale variability and dynamics within the West Greenland boundary current system that must be accounted for in future modeling.

Description: The work in this dissertation was funded by the National Science Foundation grants OCE-1259618 and OCE-1756361.

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 2021.

Description: Subduction zones host the greatest earthquakes on earth and pose great threat to human society. The largest slip in megathrust earthquakes often occurs in the 10–50 km depth range, yet seismic imaging of the material properties in this region has proven difficult. This thesis focuses on developing methods to utilize high frequency (2–12 Hz) seismic waves scattered from the megathrust plate interface to constrain its fine-scale velocity structures and to investigate the relationship between velocity structures and megathrust slip behaviors. Chapter 2 investigates the locking condition of the subducted Gorda plate by simulating afterslip that would be expected as a result of the stress changes from offshore strike-slip earthquakes. Chapter 3 develops array analysis methods to identify P-to-S and S-to-P seismic converted phases that convert at the subducted Gorda plate interface from local earthquakes and uses them to constrain the geometry and material properties of the plate boundary fault of the subducted Gorda plate between 5–20 km depth. Chapters 4 and 5 use a dense nodal array and numerical modeling methods to study the seismic guided waves that propagate along the thin low velocity layer at the boundary of the subducted Gorda plate. Taken together, our results indicate that material properties of the subduction plateboundary fault is highly heterogeneous and the plate-boundary fault is potentially contained in a low velocity layer with significant porosity and fluid content at seismogenic depths.

Description: Funding for this research was provided by National Science Foundation Division of Earth Sciences (EAR) award #1520690 and the WHOI Academic Programs Office.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Oceanography at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution February 2021.

Description: Anthropogenic emissions of greenhouse gases are driving rapid changes in ocean conditions. Shallow-water coral reefs are experiencing the brunt of these changes, including intensifying marine heatwaves (MHWs) and rapid ocean acidification (OA). Consequently, coral reefs are in broad-scale decline, threatening the livelihoods of hundreds of millions of people. Ensuring survival of coral reefs in the 21st century will thus require a new management approach that incorporates robust understanding of reef-scale climate change, the mechanisms by which these changes impact corals, and their potential for adaptation. In this thesis, I extract information from within coral skeletons to 1) Quantify the climate changes occurring on coral reefs and the effects on coral growth, 2) Identify differences in the sensitivity of coral reefs to these changes, and 3) Evaluate the adaptation potential of the keystone reef-building coral, Porites. First, I develop a mechanistic Porites growth model and reveal the physicochemical link between OA and skeletal formation. I show that the thickening (densification) of coral skeletal framework is most vulnerable to OA and that, under 21st century climate model projections, OA will reduce Porites skeletal density globally, with greatest impact in the Coral Triangle. Second, I develop an improved metric of thermal stress, and use a skeletal bleaching proxy to quantify coral responses to intensifying heatwaves in the central equatorial Pacific (CEP) since 1982. My work reveals a long history of bleaching in the CEP, and reef-specific differences in thermal tolerance linked to past heatwave exposure implying that, over time, reef communities have adapted to tolerate their unique thermal regimes. Third, I refine the Sr-U paleo-thermometer to enable monthly-resolved sea surface temperatures (SST) generation using laser ablation ICPMS. I show that laser Sr-U accurately captures CEP SST, including the frequency and amplitude of MHWs. Finally, I apply laser Sr-U to reconstruct the past 100 years of SST at Jarvis Island in the CEP, and evaluate my proxy record of bleaching severity in this context. I determine that Porites coral populations on Jarvis Island have not yet adapted to the pace of anthropogenic climate change.

Description: This research was supported by US National Science Foundation Awards OCE-1220529, ANT-1246387, OCE-1737311, CE-1601365, OCE-1805618, OCE-1537338, OCE-2016133, and from the Woods Hole Oceanographic Institution through the Ocean Life Institute, the Ocean Ventures Fund, the Grassle Fellowship Fund, and the MIT-WHOI Academic Programs Office. Additional funding was provided by the Taiwan MOST Grant 104-2628-M-001-007-MY3, the Robertson Foundation, the Leverhulme Trust in UK, the Atlantic Donor Advised Fund, The Prince Albert 2 of Monaco Foundation, the Akiko Shiraki Dynner Fund, the New England Aquarium, the Martin Family Society Fellowship for Sustainability, the Gates Millenium Scholarship, the Arthur Vining Davis Foundation, the NOAA Coral Reef Conservation Program, and from the Woods Hole Oceanographic Institution through Investment in Science Fund, the Early Career Award, and the Access to the Sea Award.

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science in Physical Oceanography at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution February 2021.

Description: Estimating turbulence in the marine-atmospheric boundary layer is critical to many industrial, commercial and scientific fields, but of particular importance to the wind energy industry. Contributing to both the efficiency of energy extraction and the life-cycle cost of the turbine itself, turbulence in the atmospheric boundary layer is estimated within the wind energy industry as Turbulence Intensity (TI) and more recently by Turbulent Kinetic Energy (TKE). Traditional in-situ methods to measure turbulence are extremely difficult to deploy in the marine environment, resulting in a recent movement to and dependence on remote sensing methods. One type of remote sensing instrument, Doppler lidars, have shown to reliably estimate the wind speed and atmospheric turbulence while being cost effective and easily deployable, and hence are being increasingly utilized as a standard for wind energy assessments. In this thesis, the ability of lidars to measure turbulence up to a height of 200 m above mean sea level in the marine-atmospheric boundary layer was tested using a 7-month data set spanning winter to early summer. Lidar-based TI and TKE were estimated by three methods using observations from a highly validated lidar system and compared under both convective and stable atmospheric stability conditions. Convective periods were found to have higher turbulence at all the heights compared to stable conditions, while mean wind speed and shear were higher during stable conditions. The study period was characterized by generally low turbulent conditions with high turbulence events occurring at timescales of a few days. Mean vertical profiles of TKE were non-uniformly distributed in height during low turbulent conditions. During highly turbulent events, TKE increased more strongly with height. The definition of TI– following the industry or meteorology conventions – had no real effect on the results, and differences between cup or sonic anemometers and lidar TI values were small except at low wind speeds. All the three lidar-based TKE methods tested corresponded closely to independent estimates, and differences between the methods were small relative to the temporal variability of TKE observed at the offshore site.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Physical Oceanography at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution June 2021.

Description: The apparent global increase in harmful algal blooms (HABs) includes Pseudo-nitzschia blooms in the Gulf of Maine, where shellfishery closures can cost millions of dollars. Temperatures in the gulf are warming, which can affect the severity of some HABs. Yet Pseudo-nitzschia in the region are understudied. Pseudo-nitzschia bloom dynamics, P. australis introduction, and potential future changes thereof were investigated in the Gulf of Maine. Data from ship surveys and moorings were used, as well as hydrodynamic, climate, and Lagrangian particle tracking models. Pseudonitzschia bloom toxicity was driven primarily by species composition, not environmental factors. P. australis was introduced to the region in 2016 via a coastal current from the Scotian Shelf. Climate change might intensify Pseudo-nitzschia blooms, shift bloom timing 1–2 weeks earlier in the spring or 4–6 weeks later in the fall, or lengthen the growing season by 3 weeks. It might also affect species composition and connectivity within the gulf. This work has implications for the monitoring of current and future blooms in the Gulf of Maine and for our understanding of HAB introduction to the region. It can also be used to develop predictive models for Pseudo-nitzschia, which could be applied to other HABs.

Description: This research was funded by the National Science Foundation (Grants OCE-1314642 and OCE-1840381), the National Institute of Environmental Health Sciences (Grants 1P01ES021923-01 and P01 ES028938-01), the Woods Hole Center for Oceans and Human Health, WHOI Academic Programs Funds, the Vannevar Bush Faculty Fellowship, and the National Oceanic and Atmospheric Administration’s HAB Event Response Program (2012 and 2016).

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 2021.

Description: The Arctic Ocean is a vital component of Earth’s climate system experiencing dramatic environmental changes. The changes are reflected in its underwater ambient soundscape as its origin and propagation are primarily dependent on properties of the ice cover and water column. The first component of this work examines the effect on ambient noise characteristics due to changes to the Beaufort Sea sound speed profile (SSP) and ice cover. Specifically, the emergence of a warm water intrusion near 70 m depth has altered the historical Arctic SSP while the ice cover has become thinner and younger due to the rise in average global temperature. Hypothesized shifts to the ambient soundscape and surface noise generation due to these changes are verified by comparing the measured noise data during two experiments to modeled results. These changes include a broadside notch in noise vertical directionality as well as a shift from uniform surface noise generation to discrete generation at specific ranges. Motivated by our data analyses, the second component presents several tools to facilitate ambient noise characterization and generation monitoring. One is a convolutional neural network (CNN) approach to noise range estimation. Its robustness to SSP and bottom depth mismatch is compared with conventional matched field processing. We further explore how the CNN approach achieves its performance by examining its intermediate outputs. Another tool is a frequency domain, transient event detection algorithm that leverages image processing and hierarchical clustering to identify and categorize noise transients in data spectrograms. The spectral content retained by this method enables insight into the generation mechanism of the detected events by the ice cover. Lastly, we present the deployment of a seismo-acoustic system to localize transient events. Two forward approaches that utilize time-difference-ofarrival are described and compared with a more conventional, inverse technique. The examination of this system’s performance prompts recommendations for future deployments. With our ambient noise analysis and algorithm development, we hope these contributions provide a stronger foundation for continued study of the Arctic ambient soundscape as the region continues to grow in significance.

Description: Office of Naval Research (ONR) via the University of California - San Diego (UCSD) under award number N00014-16-1-2129. Defense Advanced Research Projects Agency (DARPA) via Applied Physical Sciences Corp. (APS) under award number HR0011-18-C-0008. Office of Naval Research (ONR) under award number N00014-17-1-2474. Office of Naval Research (ONR) under award number N00014-19-1-2741. National Science Foundation (NSF) under grant number 2389237.

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

Description: Nearly 1.5 million people inhabit barrier islands along the U.S. Atlantic and Gulf Coasts and coastal groundwater dynamics influence the availability of freshwater, ecosystem health, pollutant transport, and flooding in these densely populated communities. However, groundwater dynamics, including the aquifer head distribution and subsurface salinity structure, in coastal aquifers are affected by multiple environmental forcings, such as waves, tides, storm surges, and precipitation that act on a variety of spatial and temporal scales, making coastal groundwater dynamics complex and difficult to predict. Here, measurements of groundwater heads, salinities, and temperatures collected for 3 years across a 550-m-wide barrier island are used in conjunction with observations of ocean tides, surge, waves, sound level, and rainfall to characterize the dynamics of the surface aquifer. Infiltration from surge, tides, and waves during storms caused up to 2 m increases in the groundwater level under the dune. The head gradients owing to these storm-induced groundwater bulges suggest flows become inland directed on the ocean-side of the island during storms. An upper saline plume (20-30 PSU) was observed above fresher (10 PSU) water up to 30 m inland of the dune face, which was the maximum wave runup location. Differences in inland propagation between tidal- and storm-induced groundwater head fluctuations are explained using analytical theories for intermediate depth aquifers. Additionally, a separate analytical water-table evolution model driven with estimated ocean shoreline water levels (based on the 36-hr-averaged offshore tide, surge, and wave height) and measured precipitation is validated by citizen-science flood reports and predicts the maximum water-table height within 0.1 m of the observed levels across the barrier island.

Description: Funding for this research was provided by the U.S. Coastal Research Program, the National Science Foundation, a National Science Foundation Graduate Research Fellowship, the Woods Hole Oceanographic ISP program, and National Security Science & Engineering and Vannevar Bush Faculty Fellowships.

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 2021.

Description: Autonomous underwater vehicles (AUVs) are an increasingly capable robotic platform, with embedded acoustic sensing to facilitate navigation, communication, and collaboration. The global positioning system (GPS), ubiquitous for air- and terrestrial-based drones, cannot position a submerged AUV. Current methods for acoustic underwater navigation employ a deterministic sound speed to convert recorded travel time into range. In acoustically complex propagation environments, however, accurate navigation is predicated on how the sound speed structure affects propagation. The Arctic’s Beaufort Gyre provides an excellent case study for this relationship via the Beaufort Lens, a recently observed influx of warm Pacific water that forms a widespread yet variable sound speed lens throughout the gyre. At short ranges, the lens intensifies multipath propagation and creates a dramatic shadow zone, deteriorating acoustic communication and navigation performance. The Arctic also poses the additional operational challenge of an ice-covered, GPSdenied environment. This dissertation demonstrates a framework for a physics-based, model-aided, real-time conversion of recorded travel time into range—the first of its kind—which was essential to the successful AUV deployment and recovery in the Beaufort Sea, in March 2020. There are three nominal steps. First, we investigate the spatio-temporal variability of the Beaufort Lens. Second, we design a human-in-the-loop graphical decision-making framework to encode desired sound speed profile information into a lightweight, digital acoustic message for onboard navigation and communication. Lastly, we embed a stochastic, ray-based prediction of the group velocity as a function of extrapolated source and receiver locations. This framework is further validated by transmissions among GPS-aided modem buoys and improved upon to rival GPS accuracy and surpass GPS precision. The Arctic is one of the most sensitive regions to climate change, and as warmer surface temperatures and shrinking sea ice extent continue to deviate from historical conditions, the region will become more accessible and navigable. Underwater robotic platforms to monitor these environmental changes, along with the inevitable rise in human traffic related to trade, fishing, tourism, and military activity, are paramount to coupling national security with international climate security.

Description: Office of Naval Research (N00014-14-1-0214) — GOATS’14 Adaptive and Collaborative Exploitation of 3-Dimensional Environmental Acoustics in Distributed Undersea Networks Draper Laboratory Incorporated (SC001-0000001039) — Positioning System for Deep Ocean Navigation (POSYDON) Office of Naval Research (N00014-16-1-2129) — MURI: The Information Content of Ocean Noise: Theory and Experiment Office of Naval Research (N00014-17-1-2474) — Environmentally Adaptive Acoustic Communication and Navigation in the New Arctic Office of Naval Research (N00014-19-1-2716) — TFO: Assessing Realism and Uncertainties in Navy Decision Aids Department of Defense, Office of Naval Research — National Defense, Science, and Engineering Graduate Fellowship

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics and Astronautics at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution June 2021.

Description: Global temperature rise and increased atmospheric carbon dioxide (CO2) levels have affected the health of the world’s ocean and water ecosystems, impacting the balances of natural carbon cycling and causing ocean acidification. Additionally, as global temperatures rise, thawing permafrost has stimulated increased release of methane (CH4), a gas with a shorter lifetime in the atmosphere but with even more heat trapping ability than CO2. In situ analysis of dissolved gas content in surface waters is currently performed with large, expensive instruments, such as spectrometers, which are coupled with gas equilibration systems, which extract dissolved gas from water and feed it to the sensor. Accurate, low cost, and portable sensors are needed to measure the dissolved CH4 and CO2 concentration in water systems to quantify their release and understand their relationship to the global carbon budget. At the same time, while greenhouse gases are well established threats to water ecosystems, the ubiquity and potential consequences of microplastics in aqueous environments are just beginning to be recognized by the environmental research community. Microplastics (MPs) are small particles of polymer debris, commonly defined as being between 1 μm and 1000 μm. Despite the pervasiveness of MPs, our ability to characterize MPs in the environment is limited by the lack of technologies for rapidly and accurately identifying and quantifying MPs. This thesis is concerned with the engineering challenges prompted by the need for high quality and quantity environmental data to better study and the impact, cycling, and prevalence of these pollutants in aqueous environments. Three distinct investigations are presented here. First, the design of the Low-Cost Gas Extraction and Measurement System (LC-GEMS) for dissolved CO2 is presented. At just under $600 dollar to build, the LC-GEMS is an ultra-portable, toolbox-sized instrument for dissolved gas sensing in near-surface waters. The LCGEMS was characterized in the lab and demonstrated linear relationships with dissolved CO2 as well as temperature. Lab calibrations and subsequent field testing in the Little Sippewissett Marsh, in Falmouth, Massachusetts showed that the LCGEMS captures both diurnal and minute-time scale trends in dissolved CO2. Second, this thesis presents the novel design of three simple and low-cost planar nanophotonic and plasmonic structures as optical transducers for measuring dissolved CH4. Through simulations, the sensitivity of the structures are evaluated and found to exhibit superior performance in the reflectance intensity readout mode to that of the standard surface-plasmon-polariton-mode Spreeta sensor. A practical, small, and low-cost implementation of this chip with a simple intensity-based measurement scheme is proposed. This design is novel in the space of dissolved gas monitoring because it shows potential to measure directly in the water phase while being robust and low-cost to implement. Finally, this thesis presents a literature review and perspective to motivate the development of field-deployable microplastic sensing techniques. A framework for field-deployable microplastic sensing is presented and seeks to inform the MP community of the potential in both traditional MP analysis techniques and unconventional methods for creating rapid and automated MP sensors. The field-deployabilty framework addresses a full scope of practical/technological trade-offs to be considered for portable MP detection.

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 September 2021.

Description: As the western boundary current of the North Atlantic, the Gulf Stream is a well established area of interest for the United States Navy, predominately due to its proximity to the continental shelf and the associated challenges of acoustic propagation across large property gradients. Autonomous underwater gliders conduct routine, high-resolution surveys along the U.S. East Coast, including within the Gulf Stream. These observations are assimilated into the operational Navy Coastal Ocean Model (NCOM). An investigation of the forecast-to-nowcast changes in the model for 2017 demonstrates the impact of the observations on the model. The magnitude of model change as a function of distance from nearest new observation reveals relatively large impact of glider observations within a radius of 𝒪(100) km. Glider observations are associated with larger local impact than Argo data, likely due to glider sampling focusing on large spatial gradients. Due to the advective nature of the Gulf Stream system, the impact of glider observations in the model is anisotropic with larger impacts extending downstream from observation locations. Forecast-to-nowcast changes in modeled temperature, salinity, and density result in improved agreement between observed and modeled ocean structure within the upper 200 m over the 24 hours between successive model runs.

Description: This research was funded via the United States Navy’s Civilian Institution Program with the Massachusetts Institute of Technology/Woods Hole Oceanographic Institution Joint Program (MIT/WHOI JP). Glider observations and analyses have been generously supported by the National Science Foundation (OCE-0220769, OCE-1558521, OCE-1633911, OCE-1923362), NOAA’s Global Ocean Monitoring and Observing Program (NA14OAR4320158, NA19OAR4320074), the Office of Naval Research (N000141713040), Eastman Chemical Corporation, WHOI’s Oceans and Climate Change Institute, and the W. Van Alan Clark, Jr. Chair for Excellence in Oceanography at WHOI (awarded to Breck Owens).

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 September 2021.

Description: This thesis analyzes data from two types of unique drifter experiments in order to characterize two aspects of ocean flows that are often difficult to study. First, vertical velocities and their associated transport processes are often challenging to observe in the real ocean since vertical velocities are typically orders of magnitude smaller than horizontal velocities in mesoscale and submesoscale flows. Second, Lagrangian coherent structures (LCS) are features which categorize ocean flows into regimes of distinct behavior. These structures are also difficult to quantify in the real ocean, since sets of gridded trajectories from real ocean data (rather than model fields) are rarely available. The first experiment uses drifters drogued at multiple depths in the Alboran Sea to observe and characterize the ocean’s vertical structure, particularly near a strong front where vertical velocities are expected to be much stronger than other regions of the Ocean. The second experiment uses a roughly gridded pattern of surface drifters in the Gulf of Mexico to study LCSs as quantified by methods from dynamical systems such as finitetime Lyapunov exponents (FTLEs), trajectory arc-length, correlation dimension, dilation, Lagrangian-averaged vorticity deviation (LAVD), and spectral clustering. This thesis includes the first attempt to apply these dynamical systems techniques to real drifters for LCS detection. Overall, these experiments and the methods used in this paper are shown to be promising new techniques for quantifying both the vertical structure of ocean flows and Lagrangian Coherent Structures of flows using real drifter data. Future work may involve modified versions of the experiments, with denser sets of ocean drifters in the horizontal and/or vertical directions.

Description: My Masters studies in the MIT/WHOI Joint Program were funded by the US Navy Civilian Institution Office.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Chemical Oceanography at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 2021.

Description: Phytoplankton are communities of diverse groups of prokaryotic and eukaryotic single-celled organisms responsible for nearly 50% of global primary production. The relative abundance of individual groups changes dynamically in response to environmental perturbations. Recent studies suggest that such changes are primarily driven by the distinct physiological responses employed by each group towards a particular perturbation. Although knowledge of some of these responses has come to light in recent years, many aspects of their metabolisms remain unknown. We attempt to address this gap by studying the metabolism of several phytoplankton groups using metabolomics. Firstly, we developed a method to enhance the analysis of untargeted metabolomics data. Secondly, we constructed two conceptual models describing how metabolism of the raphidophyte Heterosigma akashiwo responds to phosphorus and nitrogen stress. These conceptual models revealed several new stress response mechanisms not previously reported in other phytoplankton. Finally, we compared the metabolic changes of several distinct phytoplankton groups to uncover possible adaptation and acclimations that distinguish them. This analysis revealed several pathways and metabolites that represent the studied groups. The contributions of these pathways and metabolites towards physiology may support the ecological fitness of these organisms.

Description: None of this work would have been possible without a variety of funding sources. I was supported for three years by a National Science Foundation Graduate Research Fellowship and one year with a GEM fellowship. The research was carried out with grants from the MIT Microbiome Center (Award ID #6936800, EBK), the Simons Foundation (Award ID #509034, EBK), the Gordon and Betty Moore Foundation (Award ID #3304 EBK), the National Science Foundation (Award ID #OCE-0619608 to EBK and OCE-1057447 to EBK and MCKS) and the WHOI Ocean Ventures Fund.

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 2022.

Description: Attenuation from fish can reduce the intensity of acoustic signals and significantly decrease detection range for long-range active and passive sensing in the ocean. This makes it important to understand the relevant mechanisms and accurately predict attenuation from fish in underwater acoustic sensing. Formulations for predicting attenuation from fish, however, depend on the accurate characterization of population density and spatial distribution of fish groups along long-range propagation paths, which is difficult to achieve using conventional survey methods. In previous investigations of attenuation from fish, population densities were inferred from reductions in the intensity of long-range acoustic signals caused by diel or seasonal shoaling patterns of fish groups. Here, Ocean Acoustic Waveguide Remote Sensing (OAWRS) is used to instantaneously image massive Norwegian herring shoals that stretch for thousands of square kilometers and simultaneously measure attenuation from these shoals within the active OAWRS transmissions, as well as attenuation to ship-radiated tonals detected by Passive Ocean Acoustic Waveguide Remote Sensing (POAWRS). Reductions in signal intensity are predicted using a normal-mode-based analytical theory derived from first principles for acoustic propagation and scattering through inhomogeneities in an ocean waveguide. The predictions of the waveguide attenuation formulation are in agreement with measured reductions from attenuation, where the position, size, and population density of the fish groups are characterized using OAWRS imagery as well as in situ echosounder measurements of the specific shoals occluding the propagation path. Common heuristic formulations that employ free space scattering assumptions for attenuation from fish groups are not in agreement with measurements here, and waveguide scattering theory is found to be necessary for accurate predictions. It is experimentally and theoretically shown that attenuation can be significant when the sensing frequency is near the resonance frequency of the shoaling fish, where scattering losses from the fish swimbladders and damping from fish flesh is most significant. Negligible attenuation was observed in previous OAWRS and POAWRS surveys because the frequency of the acoustic signals was sufficiently far from the swimbladder resonance peak of the shoaling fish or the packing densities of the fish shoals were not sufficiently high.

Description: This work was supported by: • Office of Naval Research under grant number N00014-17-1-2197. • Office of Naval Research via the Graduate Traineeship Award under grant number N00014-18-1-2085.

Description: Anaerobic oxidation of ammonium (anammox) in oxygen minimum zones (OMZs) is a major pathway of oceanic nitrogen loss. Ammonium released from sinking particles has been suggested to fuel this process. During cruises to the Peruvian OMZ in April–June 2017 we found that anammox rates are strongly correlated with the volume of small particles (128–512 µm), even though anammox bacteria were not directly associated with particles. This suggests that the relationship between anammox rates and particles is related to the ammonium released from particles by remineralization. To investigate this, ammonium release from particles was modelled and theoretical encounters of free-living anammox bacteria with ammonium in the particle boundary layer were calculated. These results indicated that small sinking particles could be responsible for ~75% of ammonium release in anoxic waters and that free-living anammox bacteria frequently encounter ammonium in the vicinity of smaller particles. This indicates a so far underestimated role of abundant, slow-sinking small particles in controlling oceanic nutrient budgets, and furthermore implies that observations of the volume of small particles could be used to estimate N-loss across large areas.