ALBERT

Description: Rupture directivity and its potential frequency dependence is an open issue within the seismological community, especially for small‐to‐moderate events. Here, we provide a statistical overview based on empirical evidence of seismological observations, thanks to the large amount of high‐quality seismic recordings (more than 30,000 waveforms) from Central Italy, which represents an excellent and almost unique natural laboratory of normal faulting earthquakes in the magnitude range between 3.4 and 6.5 within the time frame 2008–2018. In order to detect an anisotropic distribution of ground motion amplitudes due to the rupture directivity, we fit the smoothed Fourier Amplitude Spectra (FAS) cleared of source‐, site‐ and path‐effects. According to our criteria, about 36% of the analyzed events (162 out of 456) are directive and the distribution of rupture direction is aligned with the strikes of the major faults of the Central Apennines. We find that the directivity is a band‐limited phenomenon whose width may extend up to five times the corner frequency. The results of this research provide useful insights to parameterize directivity, to be explicitly implemented in future ground motion modeling and scenario predictions.

Description: Published

Description: e2021JB023498

Description: 5T. Sismologia, geofisica e geologia per l'ingegneria sismica

Description: JCR Journal

Description: We introduce a back-projection method to locate tremor sources using products of cross-correlation envelopes of time series between seismic stations. For a given subset of n stations, we calculate the (n − 1)th-order product of cross-correlation envelopes and we stack the back-projected products over combinations of station subsets.We show that compared to existing correlation methods and for realistic signal and noise characteristics, this way of combining information can significantly reduce the effects of correlated (spurious or irrelevant signals) and uncorrelated noise. Each back-projected product constitutes an individual localized estimate of the source locations, as opposed to a hyperbola for the existing correlation techniques, assuming a uniform velocity in two dimensions. We demonstrate the method with synthetic examples and a real-data example from tremor at Katla Volcano, Iceland, in July 2011. Despite very complex near-surface structure, including strong topography and thick ice cover, the method appears to produce robust estimates of tremor location.

Description: Published

Description: 3156–3164

Description: 3T. Fisica dei terremoti e Sorgente Sismica

Description: 4V. Processi pre-eruttivi

Description: 5V. Processi eruttivi e post-eruttivi

Description: JCR Journal

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: Post-orogenic back-arc magmatism is accompanied by hydrothermal ore deposits and mineralizations derived from mantle and crustal sources. We investigate Zannone Island (ZI), back-arc Tyrrhenian basin, Italy, to define the source(s) of mineralizing hydrothermal fluids and their relationships with the regional petrological-tectonic setting. On ZI, early Miocene thrusting was overprinted by late Miocene post-orogenic extension and related hydrothermal alteration. Since active submarine hydrothermal outflow is reported close to the island, Zannone provides an ideal site to determine the P-T-X evolution of the long-lived hydrothermal system. We combined field work with microstructural analyses on syn-tectonic quartz veins and carbonate mineralizations, X-ray diffraction analysis, microthermometry and element mapping of fluid inclusions (FIs), C, O, and clumped isotopes, and analyses of noble gases (He-Ne-Ar) and CO2 content in FIs. Our results document the evolution of a fluid system of magmatic origin with increasing mixing of meteoric fluids. Magmatic fluids were responsible for quartz veins precipitation at ∼125 to 150 MPa and ∼300°C–350°C. With the onset of extensional faulting, magmatic fluids progressively interacted with carbonate rocks and mixed with meteoric fluids, leading to (a) host rock alteration with associated carbonate and minor ore mineral precipitation, (b) progressive fluid neutralization, (c) cooling of the hydrothermal system (from ∼320°C to ∼86°C), and (d) embrittlement and fracturing of the host rocks. Both quartz and carbonate mineralizations show noble gases values lower than those from the adjacent active volcanic areas and submarine hydrothermal systems, indicating that the fossil-to-active hydrothermal history is associated with the emplacement of multiple magmatic intrusions.

Description: Published

Description: e2022GC010474

Description: 1T. Struttura della Terra

Description: 2T. Deformazione crostale attiva

Description: 7T. Variazioni delle caratteristiche crostali e "precursori"

Description: 9T. Geochimica dei fluidi applicata allo studio e al monitoraggio di aree sismiche

Description: 2TR. Ricostruzione e modellazione della struttura crostale

Description: JCR Journal

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

Description: Radioactive isotopes act as nuclear clocks that are utilized to trace and measure rates of chemical, biological, physical, and geological oceanographic processes. This thesis seeks to utilize both artificial (e.g., released from anthropogenic sources) and natural radioisotopes as tracers within the Pacific Ocean basin. Artificial radioisotopes released as a result of the 2011 Fukushima Daiichi nuclear power plants accident have the potential to negatively impact human and environmental health. This study evaluates 137Cs, 90Sr, and 129I concentrations in seawater off the coast of Japan, reconciles the sources of contaminated waters, and assesses the application of 137Cs/90Sr, 129I/137Cs, and 129I/90Sr as oceanic tracers. The analysis of activity ratios suggests a variety of sources, including ongoing sporadic and independent releases of radiocontaminants. Though decreasing, concentrations remain elevated compared to preaccident levels. Future planned releases of stored water from the reactor site may affect the surrounding environment; and thus, continued efforts to understand the distribution and fate of these radionuclides are warranted. Naturally-occurring radioisotopes (e.g., the 238U-234Th series used in this thesis) can give insight into surface export and remineralization of particulate organic carbon (POC) and trace metals (TMs). POC and TMs play a vital role in regulating the biological carbon pump (BCP), which in turn helps to moderate atmospheric CO2 levels by transporting carbon to the deep ocean, where it can be sequestered on timescales of centuries to millennia. Through this thesis we utilize the 238U:234Th disequilibrium method throughout the GEOTRACES GP15 Pacific Meridional Transect in order o provide basin-scale estimates of POC export and remineralization. There is only limited, recent use of this method to constrain TM fluxes, and as such this study also seeks to further develop this method for use in understanding TM cycling through comparative flux studies in the North Pacific.

Description: Funding sources for this thesis include the Gordon and Betty Moore Foundation, the Deerbrook Charitable Trust, NSF OCE award no. 1356630 and no. 1735445, and the NSF GRFP.

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 May 2020.

Description: Highly productive marine microbial communities in the coastal Southern Ocean sustain the broader Antarctic ecosystem and play a key role in Earth’s climate via the biological pump. Regional phytoplankton growth is primarily limited by iron and co-limited by cobalamin (vitamin B12), a trace cobalt-containing organometallic compound only synthesized by some bacteria and archaea. These micronutrients impact primary production and the microbial ecology of the two keystone phytoplankton types: diatoms and Phaeocystis antarctica. This thesis investigates microbe-driven cobalamin cycling in Antarctic seas across multiple spatiotemporal scales. I conducted laboratory culture experiments with complementary proteomics and transcriptomics to investigate the B12-ecophysiology of P. antarctica strain CCMP 1871 morphotypes under iron-B12 co-limitation. We observed colony formation under higher iron treatments, and a facultative use of B12-dependent (MetH) and B12-independent (MetE) methionine synthase isoforms in response to vitamin availability, demonstrating that this strain is not B12-auxotrophic. Through comparative ’omics, we identified a putative MetE protein in P. antarctica abundant under low B12, which is also found in other marine microbes. Across Antarctic seas, community-scale cobalt and B12 uptake rates were measured by 57Co radiotracer incubation experiments and integrated with hydrographic and phytoplankton pigment data. I observed significant correlations between uptake fluxes and environmental variables, providing evidence for predominantly diatom-driven uptake of these micronutrients in warmer, fresher surface waters with notable regional differences. To date, this work is the most comprehensive attempt to elucidate the processes governing the co-cycling of cobalt and B12 in any marine system. At the ecosystem-scale, I developed and tested a hypothesis of micronutrient-driven community dynamics through a trait-based model with cross-feeding interactions. The model demonstrates how the observed seasonal succession of springtime P. antarctica from solitary to colonial cells, bacterioplankton, and summertime diatoms may be explained by the microbial cycling of iron, dissolved organic carbon, and B12. Overall, this dissertation provides new information about the micronutrient-driven ecology of Antarctic marine microbes and adds to our understanding of the interconnections between organismal life cycle, trace metals, and trace organics in marine environments.

Description: My training as a scientist during my time in the MIT–WHOI Joint Program (2014-2020) and the work presented in this dissertation were financially supported by the Academic Programs Office (APO) at the Woods Hole Oceanographic Institution (WHOI) and various funding agencies. My first semester was supported by the WHOI Von Damm Fellowship (2014). Subsequent years and endeavors were supported by awards from the Gordon and Betty Moore Foundation to Professor Michael Follows (Award 3778, M.J.F.) and Simons Collaboration on Computational Biogeochemical Modeling of Marine Ecosystems (Award 549931, M.J.F.); National Science Foundation (NSF) grant to Dr. Stephanie Dutkiewicz (Grant number 1434007, S.D.), and NSF Office of Polar Program (OPP) grant to Dr. Makoto Saito (M.S.) for the CICLOPS research expedition (OPP-1643684, OPP-1643845, and OPP-1644073).

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

Description: Little is known about how Atlantic hurricane activity changes on long timescales. This thesis uses proxy development and proxy-model integration to constrain the spatiotemporal variability in hurricane activity in the Bahama Archipelago over the past millennium. I present annually-resolved archives of storm activity derived from sediment cores from blue holes on three islands in the Bahama Archipelago: South Andros, Long Island, and Middle Caicos. Dramatic differences between these records suggest localized controls on the hurricane patterns observed by each island. Thus, compiling these records together more accurately captures regional variations in hurricane strikes. Integrating our new Bahama Archipelago compilation with compiled paleohurricane records from the U.S. coastline indicates shifting patterns of hurricane activity over the past millennium between the Gulf Coast and the Bahama Archipelago/New England. Finally, I address whether variability in hurricane strikes observed in Bahamian paleohurricane records is related to climate or random variability using hurricane model output. The signal observed in any individual record of paleohurricane activity from the Bahama Archipelago is driven more by random variability in hurricane tracks than by climate. This thesis lays the groundwork for creating high-resolution paleohurricane records from blue holes and using hurricane models to inform our interpretations of these records.

Description: This work was funded by the National Science Foundation Graduate Research Fellowship (to E.J.W.), National Science Foundation grant OCE-1356708 (to J.P.D. and P.J.vH.) and the Dalio Explore 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 2020.

Description: Speleothems, or sedimentary rocks formed in caves, act as valuable archives of past climate change due to their suitability for U-series dating and high-resolution proxy analysis. These records can provide insights into water availability and controls on hydrology prior to the instrumental record. In this thesis, I present three records from newly-analyzed Mexican stalagmites using stable isotope (oxygen and carbon) and trace element to calcium (Mg/Ca and Sr/Ca) ratios as proxies for changing hydroclimate. Chapter 2 presents a precisely dated, mid-Holocene record of high rainfall and limited precipitation variability in the Yucatan Peninsula, Mexico. Chapters 3 and 4 present novel climate records from northeastern Mexico, an understudied region of North America. Both records come from cave sites within the Mexican arid zone, which is simultaneously experiencing increased water scarcity and a rapidly growing population. In Chapter 3, I examine a speleothem from the first millennium of the Common Era, which showed that there is a precipitation dipole between northern and southern Mexico. Chapter 4 highlights, for the first time at decadal resolution, the northeast Mexican response to the 8.2 ka event and the Younger Dryas. These chapters show that the San Luis Potosí region is vulnerable to droughts under multiple climate mean states, and is subject to drying as Atlantic Meridional Overturning Circulation weakens due to anthropogenic climate change. The climate records detailed in this thesis improve our understanding of controls on Mexican hydroclimate and can serve as benchmarks for climate models.

Description: This work was funded by US National Science Foundation (NSF) grants AGS-1702848 (M. Medina-Elizalde), AGS–1502877 (S. Burns), AGS-1804512 and AGS-1806090 (K. Johnson and D. McGee). I was also supported by the NSF Graduate Research Fellowship and the MIT School of Science Dean’s Fellowship. Fieldwork and analysis were funded by the WHOI Ocean Ventures Fund, the MIT EAPS Student Research Fund, and the MIT International Science and Technology Initiatives (MISTI) Mexico program. Initial work for this project was also supported by UC MEXUS-CONACYT Collaborative Grant from the University of California Institute for Mexico and the United States (UC MEXUS CN-16-120). Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author and do not necessarily reflect the views of the National Science Foundation.

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

Description: Marine diatoms are abundant photoautotrophic algae that contribute significantly to photosynthetic carbon fixation and export throughout the oceans. Zinc is an important micronutrient in algal metabolism, with scarce dissolved concentrations in the upper euphotic zone reflecting high biological demand. In this thesis, I investigated the response of marine diatoms to Zn scarcity to characterize metabolic mechanisms used to combat Zn stress. I began by assaying the ability to metabolically substitute cobalt (Co) in place of Zn in four diatom species and found that enhanced abilities to use Co are likely an adaptation to high surface dCo:dZn ratios in the native environment. I next demonstrated that Zn/Co metabolic substitution in diatoms is not universal using culture studies of Chaetoceros neogracile RS19, which has an absolute Zn requirement. Using global proteomic analysis, I then identified and characterized diatom ZCRP-A and ZCRP-B, a putative Zn-chaperone and membrane-tethered Zn acquisition protein, respectively, as two proteins involved in the low-Zn response. I demonstrated that these proteins are widespread in marine phytoplankton and can be deployed as protein biomarkers of Zn stress in the field. I furthermore documented both the detection of ZCRPs in the Southern Ocean and the existence of Zn/Fe co-limitation within the natural phytoplankton population in Terra Nova Bay, demonstrating that Zn co-limitation can indeed occur in the field, even in high macronutrient waters. Lastly, I explored the relative demand of Zn and cadmium (Cd) within the Southern Ocean community using stable 67Zn and 110Cd tracers, documenting a high demand for both metals during the austral 2017-2018 summer season and investigating the cycling of these elements within this important region. Overall, this dissertation provides new information regarding Zn acquisition and homeostasis mechanisms within marine algae and demonstrates that Zn co-limitation in the field is not only possible, but detectable via protein biomarkers.

Description: I am very fortunate to have been financially supported over the course of my PhD by awards granted to the Saito Lab, specifically by National Institutes of Health (NIH) grant R01GM135709, Gordon and Betty Moore Foundation grant 3782, and National Science Foundation (NSF) grants OCE-1657766, OCE-1658030, OCE-1736599, OCE-1850719, and OCE-1924554.

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 Master of Science in Chemical Oceanography at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 2020.

Description: Arctic marine and lacustrine systems are experiencing rapid warming due to climate change. These changes are especially important at the interface between sediments and surface waters because they are hotspots for biogeochemical transformations such as redox reactions, nutrient consumption and regeneration, organic matter leaching and degradation, and mineral weathering. Radium isotopes (223Ra, 224Ra, 226Ra, 228Ra) and radon-222, naturally occurring radioactive isotopes produced in sediments, are well-suited as tracers of nutrients, trace metals, and organic matter cycling processes at the sediment-water interface. In this thesis, I have applied radon-222 and the quartet of radium isotopes to study fundamental processes in subarctic lakes and on the Arctic continental shelf. First, radon-222 is used to quantify groundwater discharge into a shallow, tundra lake on the Yukon-Kuskokwim Delta in Alaska in summer of 2017. Radon-derived groundwater fluxes were then paired with methane (CH4) measurements to determine delivery rates of methane into the lake via groundwater. Groundwater CH4 fluxes significantly exceeded diffusive air-water fluxes from the lake to the atmosphere, suggesting that groundwater is an important source of CH4 to Arctic lakes and may drive observed CH4 emissions. Higher CH4 emissions were observed compared to those reported previously in high latitude lakes, like due to higher CH4 concentrations in groundwater. These findings indicate that deltaic lakes across warmer permafrost regions may act as important hotspots for methane release across Arctic landscapes. Then, the quartet of radium isotopes is used to study the impacts of storms and sea ice formation as drivers of sediment-water interaction on the Alaskan Beaufort shelf. The timeseries presented in this study is among the first to document the combined physical and chemical signals of winter water formation in the Beaufort Sea, made possible by repeat occupations of the central Beaufort shelf. Radium measurements are combined with inorganic nitrogen and hydrographic measurements to elucidate the episodic behavior of winter water formation and its ability to drive exchange with bottom sediments during freeze-up.

Description: Financial support for Chapter 2 was funded by National Science Foundation awards OCE-1458305 to M.A.C., 1561437 to S.M.N, J.D.S., and R.M.H and 1624927 to S.M.N., P.J.M. and R.M.H. The work completed for Chapter 3 was funded by the Montrym Fund at the Massachusetts Institute of Technology, the Academic Programs Office at Woods Hole Oceanographic Institution, and the NSF Arctic GEOTRACES (OCE-1458305), Pacific GEOTRACES (OCE-1736277), and Arctic Observing Network programs (OPP-1733564).

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

Description: Robotic swarms are increasingly complex above the waterline due to reliable communication links. However, the limited propagation of similar signals in the ocean has impacted advances in undersea robotics. Underwater vehicles often rely on acoustics for navigation solutions; however, this presents challenges for robotic swarms. Many localization methods rely on precision time synchronization or two-way communication to estimate ranges. The cost of Chip-scale Atomic Clocks (CSACs) and acoustic modems is limiting for large-scale swarms due to the cost-per-vehicle and communications structure. We propose a single vehicle with reliable navigation as a "leader" for a scalable swarm of lower-cost vehicles that receive signals via a single hydrophone. This thesis outlines range estimation methods for sources with known signal content, including frequency and power at its origin. Transmission loss is calculated based on sound absorption in seawater and geometric spreading loss to estimate range through the Signal Absorption-Based Range Estimator (SABRE). SABRE's objective is to address techniques that support low-cost undersea swarming. This thesis's contributions include a novel method for range estimation onboard underwater vehicles that supports relative navigation through Doppler-shift methods for target bearing. This thesis develops the theory, algorithms, and analytical tools for real-world data range estimation.

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

Description: This thesis examines the transition of a vessel from the open ocean, where collisions are rare, to a high risk and heavy traffic area such as a Traffic Separation Scheme (TSS). Previous autonomy approaches generally view path planning and collision avoidance as two separate functions, i.e. a vessel will follow the planned path until conditions are met for collision avoidance algorithms to take over. Here an intermediate phase is proposed with the goal of adjusting the time of arrival to a high vessel density area so that the risk of collision is reduced. A general algorithm that calculates maximum future traffic density for all choices in the speed domain is proposed and implemented as a MOOS-IvP behavior. This behavior gives the vessel awareness of future collision risks and aids the collision avoidance process. This new approach improves the safety of the vessel by reducing the number of risky encounters that will likely require the vessel to maneuver for safety.

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: The early life stages of marine fishes play a critical role in population dynamics, largely due to their high abundance, high mortality, and ease of transport in ocean currents. This dissertation demonstrates the value of combining larval data, collected in the field and the laboratory, with model simulations. In Chapter 2, analyses of field observations of ontogenetic vertical distributions of coral reef fish revealed a diversity of behaviors both between and within families. In Caribbean-wide particle-tracking simulations of representative behaviors, surface-dwelling larvae were generally transported longer distances with greater population connectivity amongst habitat patches, while the evenly-distributed vertical behavior and downward ontogenetic vertical migration were similar to one another and led to greater retention near natal sites. However, hydrodynamics and habitat vailability created some local patterns that contradicted the overall expectation. Chapter 3 presents evidence of tuna spawning inside a large no-take marine protected area, the Phoenix Islands Protected Area (PIPA). Despite variation in temperature and chlorophyll, the larval tuna distributions were similar amongst years, with skipjack (Katsuwonus pelamis) and Thunnus spp. tunas observed in all three years. Backtracking simulations indicated that spawning occurred inside PIPA in all 3 study years, demonstrating that PIPA is protecting viable tuna spawning habitat. In Chapter 4, several lines of larval evidence support the classification of the Slope Sea as a major spawning ground for Atlantic bluefin tuna with conditions suitable for larval growth. The abundance of bluefin tuna larvae observed in the Slope Sea aligns with typical observations on the other two spawning grounds. Age and growth analyses of bluefin tuna larvae collected in the Slope Sea and the Gulf of Mexico in 2016 did not show a growth rate difference between regions, but did suggest that Slope Sea larvae are larger at the onset of exogenous feeding. Collected larvae were backtracked to locations north of Cape Hatteras and forward tracked to show that they would have been retained within the Slope Sea until the onset of swimming. As a whole, this thesis presents valuable contributions to the study of larval fishes and the attendant implications for marine resource management.

Description: National Science Foundation Graduate Research Fellowship Program (to C.M.H.), the Woods Hole Oceanographic Institution Ocean Life Institute (Grant 22569.01 to J. Llopiz and C.M.H.), the Adelaide and Charles Link Foundation, the Phoenix Islands Protected Area Trust, the J. Seward Johnson Endowment in support of the Woods Hole Oceanographic Institution's Marine Policy Center, 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: 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 Doctor of Philosophy in Mechanical and Oceanographic Engineering at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution February 2021.

Description: The satellite ocean color remote sensing paradigm developed by government space agencies enables the assessment of ocean color products on global scales at kilometer resolutions. A similar paradigm has not yet been developed for regional scales at sub-meter resolutions, but it is essential for specific ocean color applications (e.g., mapping algal biomass in the marginal ice zone). While many aspects of the satellite ocean color remote sensing paradigm are applicable to sub-meter scales, steps within the paradigm must be adapted to the optical character of the ocean at these scales and the opto-electronics of the available sensing instruments. This dissertation adapts the three steps of the satellite ocean color remote sensing paradigm that benefit the most from reassessment at sub-meter scales, namely the correction for surface-reflected light, the design and selection of the opto-electronics and the post-processing of over-sampled regions. First, I identify which surface-reflected light removal algorithm and view angle combination are optimal at sub-meter scales, using data collected during a field deployment to the Martha’s Vineyard Coastal Observatory. I find that of the three most widely used glint correction algorithms, a spectral optimization based approach applied to measurements with a 40∘ view angle best recovers the remotesensing reflectance and chlorophyll concentration despite centimeter scale variability in the surface-reflected light. Second, I develop a simulation framework to assess the impact of higher optical and electronics noise on ocean color product retrieval from unique ocean color scenarios. I demonstrate the framework’s power as a design tool by identifying hardware limitations, and developing potential solutions, for estimating algal biomass from high dynamic range sensing in the marginal ice zone. Third, I investigate a spectral super-resolution technique for application to spatially over-sampled oceanic regions. I determine that this technique more accurately represents spectral frequencies beyond the Nyquist and that it can be trained to be invariant to noise sources characteristic of ocean color remote sensing on images with similar statistics as the training dataset. Overall, the developed and critically assessed sub-meter ocean color remote sensing paradigm enables researchers to collect high fidelity sub-meter data from imaging spectrometers in unique ocean color scenarios.

Description: Ryan O’Shea was supported by the Department of Defense (DoD) through the National Defense Science & Engineering Graduate Fellowship (NDSEG) Program. This research was funded by Woods Hole Oceanographic Institution’s Edwin W. Hiam Ocean Science and Technology Award Fund, its Ocean Venture Funds, its Academic Programs Office, and the National Aeronautics and Space Administration via grant number CCE NNX17AI72G to Dr. Samuel Laney. The raw data for Figures 3-3 and 3-4 were provided through Australian Antarctic Science grants 2678 and 4390.

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 Master of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution February 2022.

Description: Diving operations are inherently complex due to navigation and communication limitations. Until recently, fixed-beacon acoustic localization techniques have served as the primary means of improving diver navigation. However, modern artificial intelligence and acoustic modem technologies have enabled accurate relative navigation methods between a diver and an autonomous vehicle. Human-robot collaboration takes advantage of each member’s strengths to create the most effective team. This concept proves especially advantageous within the ocean domain, where humans are naturally deficient navigators. Yet humans serve as the team’s creative spirit, offering the critical thinking and flexibility needed to succeed in an unpredictable and dynamic environment. Recent underwater human-robot cooperative navigation systems typically rely on autonomous surface vehicles (ASVs), specially designed underwater vehicles, or stereo cameras. This thesis proposes a diver navigation method exhibiting significantly improved accuracy over dead reckoning without relying on a surface presence, cameras, or fixed acoustic beacons. Specifically, we develop and evaluate the communication architecture and autonomous behaviors required to guide a diver to a target location using subsurface humanautonomous underwater vehicle (AUV) teaming with no requirement for ocean current data or exact diver speeds. By depending on acoustic communication and commercial AUV navigation capabilities, our method has increased accessibility, applicability, and robustness over former techniques. We utilize the Woods Hole Oceanographic Institution (WHOI) Micromodem 2’s twoway-travel-time (TWTT) capability to enable range-only single-beacon navigation between two kayaks serving as proxies for the diver and Remote Environmental Monitoring Units (REMUS) 100 AUV. During processing, a nonlinear least-squares (NLS) method, called incremental smoothing and mapping 2 (iSAM2), utilizes odometry and range measurements to provide real-time diver position estimates given unknown ocean currents. Field experiments demonstrate an average online endpoint error of 4.53 meters after transits four hundred meters long. Additionally, simulations test our method’s performance in more challenging situations than those experienced in the field. Overall, this research progresses the interoperability of divers and AUVs.

Description: The United States Navy funded my graduate education. The Office of Naval Research also partially supported this work under grant N00014-18-1-2832.

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

Description: Marine microbes require copper (Cu) for a variety of key enzymes and can therefore experience limitation when concentrations are low. However, when Cu concentrations are too high, it becomes toxic causing decreased cell growth and even cell death. Laboratory culture experiments have shown that a diverse array of microbes produce organic ligands that complex Cu (CuL) and buffer the free ion concentration, which is the most bioavailable fraction. In this way, the microbes impose a control on the speciation of Cu, decreasing the toxic effects of Cu and making seawater conditions favorable for growth. Studies have shown that CuL complexes produced in laboratory cultures have similar complexation strengths to those found in seawater samples, which suggests a biological source of CuLs in seawater where dissolved Cu is almost entirely bound by organic ligands. However, information about individual CuL complexes is lacking which limits our understanding of the sources, sinks, and cycling of dissolved Cu. In order to fill this gap in knowledge, molecular level information about CuL complexes produced in culture and found in seawater must be obtained. To investigate this, liquid chromatography (LC) was coupled to two mass spectrometers (MS), an inductively coupled plasma (ICP) MS and an electrospray ionization (ESI) MS. By using data supplied by both techniques, the molecular charateristics of CuLs were determined laboratory cultures of the marine diatom Phaeodactylum tricornutum and the cyanobacterium Synechococcus, as well as investigating the distribution of CuLs in natural seawater samples along a line from 56°N to 20°S, along 152°W through the north and central Pacific Ocean. The CuLs identified in laboratory cultures had molecular formulae and fragmentation patterns characteristic of linear tetrapyrroles, a group of organic compounds commonly found in biological systems. This identification was further supported by absorbance and nuclear magnetic resonance spectroscopy. The distribution of CuLs in the Pacific Ocean showed a highly dynamic and complex mixture of ligands, closely tied to biological cycles.

Description: Funding for this work was provided by the National Science Foundation Graduate Research Fellowship Program (NSF award 1122374) for providing three years of funding. Thank you to the National Science Foundation Chemical Oceanography Program (NSF award OCE-1736280 and OCE-2045223) and the Simons Collaboration on Ocean Processes and Ecology (award P49476). A portion of this work was performed at the National High Magnetic Field Laboratory in Tallahassee, Florida, which is supported through NSF DMR 11-57490, and the State of Florida.

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

Description: Removal of particulate organic carbon (POC) from sunlit surface waters into the deep ocean represents a climatically important sink of atmospheric carbon dioxide (CO2), linking the biogeochemical cycling of POC to CO2-driven climate change. As POC is not well preserved in the sediment record, other proxies, including the chemistry of barium (Ba) in the ocean and through the sedimentary record, offer an avenue to investigate oceanic carbon export through Earth’s history. This thesis seeks to constrain the controls on the formation, cycling, and isotopic signature of the main particulate phase of marine barium, the mineral barite (BaSO4) through its inception in the water column, during deposition, and ultimately into the rock record. To that end, I characterize the depth, spatial region, and general controls on particulate Ba formation in the South Pacific Ocean through shipboard experimentation and find that particulate Ba forms mainly in the surface of the Polar Frontal Zone in the presence of large particles and microbial activity. Next, I characterize the effect of ion exchange on BaSO4, a process previously unstudied under marine conditions, in a laboratory setting. Ion exchange occurs rapidly between dissolved Ba and BaSO4 and imparts a characteristic net offset between the Ba isotope composition of the dissolved and solid phase, which arises through a combination of Ba isotope fractionation during both precipitation and dissolution. Finally, I investigate the role of ion exchange in marine settings using co-located pore fluids and sedimented BaSO4. Modeling constrained by data from natural samples produce results that are consistent with the laboratory study, suggesting that this mode of isotopic fractionation impacts Ba isotopes in the environment and must be accounted for when applying Ba based climate proxies.

Description: Funding for this work was provided by the National Science Foundation (OCE-2023456 & OCE-1827401), the Woods Hole Oceanographic Institution Ocean Ventures Fund, a National Science Foundation Graduate Research Fellowship (2017250048), and Woods Hole Oceanographic Institution.

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: The Gulf Stream, the western boundary current of the subtropical North Atlantic, plays a key role in the Earth’s climate system with its poleward volume and heat transports being major components of the upper limb of the Atlantic Meridional Overturning Circulation. Extensive observations collected using Spray autonomous underwater gliders from 2004 through 2020 fill a 1500-km-long gap in longer-term sustained subsurface measurements of the Gulf Stream. The gliders provide concurrent, high-resolution measurements of Gulf Stream hydrography and velocity over more than 15 degrees of latitude between Florida and New England. These observations are used to characterize the along-stream evolution of Gulf Stream volume transport; its long-known poleward increase is shown to result primarily from entrainment of subthermocline waters. Antarctic Intermediate Water, which makes up the deepest waters within the Gulf Stream in the Florida Strait, is eroded through both vertical mixing and lateral stirring as it flows downstream. Satellite-based observations of sea surface height coincident with the glider observations are used to evaluate the efficacy of inferring Gulf Stream transport from remotely sensed measurements. The detailed analyses of Gulf Stream transport and water property evolution herein provide targets for regional and global circulation models to replicate.

Description: We gratefully acknowledge funding from the National Science Foundation (OCE-0220769, OCE-1633911, OCE-1923362, OCE-1558521), NOAA’s Global Ocean Monitoring and Observing Program (NA14OAR4320158, NA19OAR4320074), the Office of Naval Research (N000141713040), WHOI’s Oceans and Climate Change Institute, Eastman Chemical Company, 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 Doctor of Philosophy in Chemical Oceanography at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution February 2022.

Description: Coastal ecosystems provide key services that benefit human wellbeing yet are undergoing rapid degradation due to natural and anthropogenic pressures. This thesis seeks to understand how disturbances impact salt marsh and estuarine ecosystem functioning in order to refine their role in coastal ecosystem service delivery and predict future resilience. Salt marsh survival relative to sealevel rise increasingly relies on the accumulation and preservation of soil organic carbon (SOC). Firstly, I characterized SOC development and turnover in a New England salt marsh and found that salt marsh soils typically store marsh grass-derived compounds that are reworked over centuries-to-millennia. Next, I assessed how two common marsh disturbances – natural ponding and anthropogenic mosquito ditching – affect salt marsh carbon cycling and storage. Salt marsh ponds deepen through soil erosion and decomposition of long-buried marsh peat. Further, the SOC lost during pond development is not fully recouped once drained ponds are revegetated and virtually indistinguishable from the surrounding marsh. Mosquito ditches, which were installed in ~ 90% of New England salt marshes during the Great Depression, did not significantly alter marsh carbon storage. In Buzzards Bay, Massachusetts, a US National Estuary, we tested relationships among measures of estuarine water quality, recreational activity, and local socioeconomic conditions to understand how the benefits of cultural ecosystem services are affected by shifts in water quality associated with global change and anthropogenic activity. Over a 24-year period, water quality degradation coinciding with increases in Chlorophyll a is associated with declines in fishery abundance and cultural ecosystem service values ($0.08 – 0.67 million USD). In combination, incorporation of both anthropogenic and natural disturbances to coastal ecosystem functioning and service delivery can produce improved estimates of ecosystem service valuation for effective resource decision-making under future climate scenarios.

Description: Funding for this work was provided by John D. and Catherine T. MacArthur Foundation (Grant no. 14-106159-000-CFP), National Science Foundation (OCE1233678), National Oceanic and Atmospheric Administration, National Oceanic and Atmospheric Administration – National Estuaries Research Reserve Collaborative (NA14OAR4170104 and NA- 14NOS4190145), Woods Hole Sea Grant (NA14OAR4170104), MIT Sea Grant (subaward number 5710004045), Ocean Ventures Fund, the Marine Policy Center Johnson Endowment, and Woods Hole Oceanographic Institution.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Marine Geology & Geophysics at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution May 2022.

Description: It is a scientifically accepted fact that the Earth’s climate is presently undergoing significant changes with the potential for immense negative impacts on human society. As evidence of these impacts become clear and common, it becomes ever more important to constrain the nature, magnitude, and speed of changes to Earth systems. A fundamentally important tool to this understanding is the Earth’s past, recorded in the geologic record. There, lie examples of climate change under various forcings: important data for understanding the fundamental dynamics of climate change on our planet. However, when a climate signal is written in the geologic record, it is coded into the language of proxies and distorted by time. This thesis endeavors to decode that record using a variety of computational methods on a number of challenging proxies, to draw more information from the climate past than has previously been possible. First, machine learning and computer vision are used to decipher the primary, centimeter-scale textures of carbonate deposits in Searles Valley and Mono Lake, California. This work is able to connect facies in the tufa at Searles, grown during the Last Glacial Period, and those forming presently at Mono Lake. Next, the tracks of icebergs purged during Heinrich Events are simulated using the MIT General Circulation Model. This work, running multiple experiments exploring different aspects internal and external to the icebergs, reveals wind and sediment partitioning as centrally important to the spatial extent of Heinrich Layers. Each of these works considers a traditional geologic archive – a carbonate facies, a marine sediment layer – and uses computational methods to approach that archive from a different perspective. By applying these new methods, more information can be gleaned from the geologic record, building a richer narrative of the Earth’s climate history. The final chapter of this thesis discusses effective teaching and strategies for building communities to support teaching practice in Earth Science departments.

Description: This thesis work was funded by the MIT EAPS Rasmussen and Whiteman Fellowships, NSF Project Number NSF-EAR-1903544, and the WHOI Academic Programs Office.

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

Description: This thesis presents an Autonomous Underwater Glider (AUG) architecture that is intended for basin-scale unattended survey of Arctic sea-ice. The distinguishing challenge for AUG operations in the Arctic environment is the presence of year-round sea-ice cover which prevents vehicle surfacing for localization updates and shore-side communication. Due to the high cost of operating support vessels in the Arctic, the proposed AUG architecture minimizes external infrastructure requirements to brief and infrequent satellite updates on the order of once per day. This is possible by employing onboard acoustic sensing for sea-ice observation and navigation, along with intelligent management of onboard resources. To enable unattended survey of Arctic sea-ice with an AUG, this thesis proposes a hierarchical acoustics-based sea-ice characterization scheme to perform science data collection and assess environment risk, a multi-factor terrain-aided navigation method that leverages bathymetric features and active ocean current sensing to limit localization error, and a set of energy-optimal propulsive and hotel policies that react to evolving environmental conditions to improve AUG endurance. These methods are evaluated with respect to laboratory experiments and preliminary field data, and future Arctic sea-ice survey mission concepts are discussed.

Description: Support for this research was provided through the National Science Foundation Navigating the New Arctic Grant #1839063 and the NASA PSTAR Grant #NNX16AL08G. Additionally, this research was supported by the Walter A. Rosenblith Presidential Fellowship.

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science in Oceanography/Applied Ocean Science and Engineering at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 2020.

Description: Underwater Vehicles generally have control fins located only near their aft end, for making controllable changes in directions. This design allows for stability of control; but, the turns are typically large in comparison to the vehicle’s body length. Some bony fish, such as tuna, however, have deployable fins located towards the front of their body, in addition to their other fins. Their deployable fins allow them to modulate their hydrodynamic behavior in response to their environment. Tunas keep these fins retracted during steady cruising, and then deploy them during rapid maneuvers. However, the details of these hydrodynamic effects are not well understood. To investigate this phenomena, using a REMUS 100 as a model, a pair of vertical fins was added at different hull positions, to investigate the effects of fin location on the horizontal plane hydrodynamics, through: stability parameters, nonlinear simulation, and towing tank experiments. Depending on the added fin location, the vehicle stability changed, thereby affecting the maneuverability. As fins were placed forward on the vehicle, maneuverability increased, with effects tapering off at 0.2 BL ahead of the vehicle's center of buoyancy. This investigation explored how rigid underwater vehicles could benefit from added fins, without drastically changing the design of current vehicles.

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

Description: Sea ice thickness has long been an under-measured quantity, even in the satellite era. The snow surface elevation, which is far easier to measure, cannot be directly converted into sea ice thickness estimates without knowledge or assumption of what proportion of the snow surface consists of snow and ice. We do not fully understand how snow is distributed upon sea ice, in particular around areas with surface deformation. Here, we show that deep learning methods can be used to directly predict snow depth, as well as sea ice thickness, from measurements of surface topography obtained from laser altimetry. We also show that snow surfaces can be texturally distinguished, and that texturally-similar segments have similar snow depths. This can be used to predict snow depth at both local (sub-kilometer) and satellite (25 km) scales with much lower error and bias, and with greater ability to distinguish inter-annual and regional variability than current methods using linear regressions. We find that sea ice thickness can be estimated to ∼20% error at the kilometer scale. The success of deep learning methods to predict snow depth and sea ice thickness suggests that such methods may be also applied to temporally/spatially larger datasets like ICESat-2.

Description: This research was funded by National Aeronautics and Space Administration grant numbers NNX15AC69G and 80NSSC20K0972, the US National Science Foundation grant numbers ANT-1341513, ANT-1341606, ANT-1142075 and ANT-1341717, and the WHOI Academic Programs Office.

Description: At Mt. Etna volcano, flank eruptions are often accompanied by seismic swarms with damaging earthquakes; the most recent case-history is the 2018 flank eruption, associated with a destructive earthquake (Mw 5.0). In this paper, we analyze the earthquake and eruptive catalogs from 1800 to 2018, to produce quantitative estimates of the earthquake rate under the influence of flank eruptions. We quantify that 30% of the flank eruption onsets precede a major (i.e., damaging, Ix ≥ V–VI EMS) earthquake by 30 days or less; 18% of the major earthquakes follow a flank eruption onset in 30 days or less. Thus, we show that the probability of major earthquakes increases 5–10 times after the onset of flank eruptions and this effect lasts for 30–45 days. This is also observed after the end of the eruptions. Results indicate different relationships depending on the location of the volcano-tectonic systems considered individually (eruptive fissures, seismogenic faults). For instance, we describe a 10–20 times increased probability of earthquakes for 65–70 days after eruptions on the northeastern flank, and of new flank eruptions for 45–70 days after earthquakes of the Pernicana fault.

Description: Published

Description: e2022JB024145

Description: 6V. Pericolosità vulcanica e contributi alla stima del rischio

Description: JCR Journal