Thesis: Ph. D., Joint Program in Physical Oceanography (Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2017. / Cataloged from PDF version of thesis. / Includes bibliographical references (pages 163-174). / The Gulf Stream and Deep Western Boundary Current (DWBC) shape the distribution of heat and carbon in the North Atlantic, with consequences for global climate. This thesis employs a combination of theory, observations and models to probe the dynamics of these two western boundary currents. First, to diagnose the dynamical balance of the Gulf Stream, a depth-averaged vorticity budget framework is developed. This framework is applied to observations and a state estimate in the subtropical North Atlantic. Budget terms indicate a primary balance of vorticity between wind stress forcing and dissipation, and that the Gulf Stream has a significant inertial component. The next chapter weighs in on an ongoing debate over how the deep ocean is filled with water from high latitude sources. Measurements of the DWBC at Line W, on the continental slope southeast of New England, reveal water mass changes that are consistent with changes in the Labrador Sea, one of the sources of deep water thousands of kilometers upstream. Coherent patterns of change are also found along the path of the DWBC. These changes are consistent with an advective-diffusive model, which is used to quantify transit time distributions between the Labrador Sea and Line W. Advection and stirring are both found to play leading order roles in the propagation of water mass anomalies in the DWBC. The final study brings the two currents together in a quasi-geostrophic process model, focusing on the interaction between the Gulf Stream's northern recirculation gyre and the continental slope along which the DWBC travels. We demonstrate that the continental slope restricts the extent of the recirculation gyre and alters its forcing mechanisms. The recirculation gyre can also merge with the DWBC at depth, and its adjustment is associated with eddy fluxes that stir the DWBC with the interior. This thesis provides a quantitative description of the structure of the overturning circulation in the western North Atlantic, which is an important step towards understanding its role in the climate system. / by Isabela Astiz Le Bras. / Ph. D.

Thesis (Ph. D.)--Joint Program in Physical Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2011. / Cataloged from PDF version of thesis. / Includes bibliographical references (p. 221-227). / Eastern oceanic boundary currents are subject to hydrodynamic instability, generate small scale features that are visible in satellite images and may radiate westward into the interior, where they can be modified by the large-scale circulations. This thesis studies the stability of an eastern boundary current with and without the large-scale flow influence in an idealized framework represented by barotropic quasi-geostrophic dynamics. The linear stability analysis of a meridional current with a continuous velocity profile shows that meridional eastern and western boundary currents support a limited number of radiating modes with long meridional and zonal wavelengths and small growth rates. However, the linearly stable, long radiating modes of an eastern boundary current can become nonlinearly unstable by resonating with short trapped unstable modes. This phenomenon is clearly demonstrated in the weakly nonlinear simulations. Results suggest that linearly stable longwave modes deserve more attention when the radiating instability of a meridional boundary current is considered. A large-scale flow affects the short trapped unstable mode and long radiating mode through different mechanisms. The large-scale flow modifies the structure of the boundary current to stabilize or destabilize the unstable modes, leading to a meridionally localized maximum in the perturbation kinetic energy field. The shortwave mode is accelerated or decelerated by the meridional velocity adjustment of the large-scale flow to have an elongated or a squeezed meridional structure, which is confirmed both in a linear WKB analysis and in nonlinear simulations. The squeezed or elongated unstable mode detunes the nonlinear resonance with the longwave modes, which then become less energetic. These two modes show different meridional structures in kinetic energy field because of the different mechanisms. In spite of the model simplicity, these results can potentially explain the formation of the zonal jets observed in altimeter data, and indicate the influence of the large-scale wind-driven circulation on eastern boundary upwelling systems in the real ocean. Studies with more realistic configurations remain future challenges. / by Jinbo Wang. / Ph.D.

Thesis: Ph. D., Joint Program in Physical Oceanography (Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2014. / Cataloged from PDF version of thesis. / Includes bibliographical references (pages 175-180). / This thesis explores the buoyancy-driven circulation in the Red Sea, using a combination of observations, as well as numerical modeling and analytical method. The first part of the thesis investigates the formation mechanism and spreading of Red Sea Overflow Water (RSOW) in the Red Sea. The preconditions required for open-ocean convection, which is suggested to be the formation mechanism of RSOW, are examined. The RSOW is identified and tracked as a layer with minimum potential vorticity and maximum chlorofluorocarbon-12. The pathway of the RSOW is also explored using numerical simulation. If diffusivity is not considered, the production rate of the RSOW is estimated to be 0.63 Sv using Walin's method. By comparing this 0.63 Sv to the actual RSOW transport at the Strait of Bab el Mandeb, it is implied that the vertical diffusivity is about 3.4 x10-5 m 2 s-1. The second part of the thesis studies buoyancy-forced circulation in an idealized Red Sea. Buoyancy-loss driven circulation in marginal seas is usually dominated by cyclonic boundary currents on f-plane, as suggested by previous observations and numerical modeling. This thesis suggests that by including [beta]-effect and buoyancy loss that increases linearly with latitude, the resultant mean Red Sea circulation consists of an anticyclonic gyre in the south and a cyclonic gyre in the north. In mid-basin, the northward surface flow crosses from the western boundary to the eastern boundary. The observational support is also reviewed. The mechanism that controls the crossover of boundary currents is further explored using an ad hoc analytical model based on PV dynamics. This ad hoc analytical model successfully predicts the crossover latitude of boundary currents. It suggests that the competition between advection of planetary vorticity and buoyancy-loss related term determines the crossover latitude. The third part of the thesis investigates three mechanisms that might account for eddy generation in the Red Sea, by conducting a series of numerical experiments. The three mechanisms are: i) baroclinic instability; ii) meridional structure of surface buoyancy losses; iii) cross-basin wind fields. / by Ping Zhai. / Ph. D.

Thesis (S.M.)--Joint Program in Physical Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2008. / Includes bibliographical references (p. 57-58). / The mid-to-deep Arctic Ocean is generally characterized by a cyclonic circulation, contained along shelves and ridges. Here we analyze the general Arctic circulation using an idealized numerical model consisting of a circular basin with two channels acting as inflow and outflow. We analyze the circulation (direction, strength and sensitivity) for wind forcing with and without bathymetry (ridges), and with and without stratification. We find that the circulation is modified drastically by both bathymetry and wind direction, where an altered wind field can change both the direction of the horizontal basin circulation as well as the strength of the inflow and outflow. The idealized circulations imply that the Arctic circulation, and the associated export of freshwater, can easily switch states in a changing climate. / by Peter Joseph Sugimura. / S.M.

Thesis (Ph. D.)--Joint Program in Physical Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2008. / Includes bibliographical references (p. 207-221). / This thesis studies the dynamics of a rotating compressible gas sphere, driven by internal convection, as a model for the dynamics on the giant planets. We develop a new general circulation model for the Jovian atmosphere, based on the MITgcm dynamical core augmenting the nonhydrostatic model. The grid extends deep into the planet's interior allowing the model to compute the dynamics of a whole sphere of gas rather than a spherical shell (including the strong variations in gravity and the equation of state). Different from most previous 3D convection models, this model is anelastic rather than Boussinesq and thereby incorporates the full density variation of the planet. We show that the density gradients caused by convection drive the system away from an isentropic and therefore barotropic state as previously assumed, leading to significant baroclinic shear. This shear is concentrated mainly in the upper levels and associated with baroclinic compressibility effects. The interior flow organizes in large cyclonically rotating columnar eddies parallel to the rotation axis, which drive upgradient angular momentum eddy fluxes, generating the observed equatorial superrotation. Heat fluxes align with the axis of rotation, contributing to the observed flat meridional emission. We show the transition from weak convection cases with symmetric spiraling columnar modes similar to those found in previous analytic linear theory, to more turbulent cases which exhibit similar, though less regular and solely cyclonic, convection columns which manifest on the surface in the form of waves embedded within the superrotation. We develop a mechanical understanding of this system and scaling laws by studying simpler configurations and the dependence on physical properties such as the rotation period, bottom boundary location and forcing structure. These columnar cyclonic structures propagate eastward, driven by dynamics similar to that of a Rossby wave except that the restoring planetary vorticity gradient is in the opposite direction, due to the spherical geometry in the interior. / (cont.) We further study these interior dynamics using a simplified barotropic annulus model, which shows that the planetary vorticity radial variation causes the eddy angular momentum flux divergence, which drives the superrotating equatorial flow. In addition we study the interaction of the interior dynamics with a stable exterior weather layer, using a quasigeostrophic two layer channel model on a beta plane, where the columnar interior is therefore represented by a negative beta effect. We find that baroclinic instability of even a weak shear can drive strong, stable multiple zonal jets. For this model we find an analytic nonlinear solution, truncated to one growing mode, that exhibits a multiple jet meridional structure, driven by the nonlinear interaction between the eddies. Finally, given the density field from our 3D convection model we derive the high order gravitational spectra of Jupiter, which is a measurable quantity for the upcoming JUNO mission to Jupiter. / by Yohai Kaspi. / Ph.D.

Thesis (Ph. D.)--Joint Program in Physical Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2010. / Cataloged from PDF version of thesis. / Includes bibliographical references (p. 131-149). / Aspects of the circulation and convection in the Irminger Sea are investigated using a variety of in-situ, satellite, and atmospheric reanalysis products. Westerly Greenland tip jet events are intense, small-scale wind phenomena located east of Cape Farewell, and are important to circulation and convection in the Irminger Sea. A climatology of such events was used to investigate their evolution and mechanism of generation. The air parcels constituting the tip jet are shown to have a continental origin, and to exhibit a characteristic deflection and acceleration around southern Greenland. The events are almost invariably accompanied both by a notable coherence of the lower-level tip jet with an overlying upper-level jet stream, and by a surface cyclone located in the lee (east) of Greenland. It is argued that the tip jet arises from the interplay of the synopticscale flow evolution and the perturbing effects of Greenland's topography upon the flow. The Irminger Gyre is a narrow, cyclonic recirculation confined to the southwest Irminger Sea. While the gyre's existence has been previously documented, relatively little is known about its specific features or variability. The mean strength of the gyre's circulation between 1991 and 2007 was 6.8 ± 1.8 Sv. It intensified at a rate of 4.3 Sv per decade over the observed period despite declining atmospheric forcing. Examination of the temporal evolution of the LSW layer thickness across the Irminger Basin suggests that local convection formed LSW during the early 1990s within the Irminger Gyre. In contrast, LSW appeared outside of the gyre in the eastern part of the Irminger Sea with a time lag of 2-3 years, consistent with transit from a remote source in the Labrador Sea. In the winter of 2007-08 deep convection returned to both the Labrador and Irminger seas following years of shallow overturning. The transition to a convective state took place abruptly, without going through a preconditioning phase, which is contrary to general expectations. Changes in the hemispheric air temperature, tracks of storms, flux of freshwater to the Labrador Sea, and distribution of pack ice all conspired to enhance the air-sea heat flux, resulting in the deep overturning. / Kjetil Våge. / Ph.D.

The Monterey Bay Aquarium Research Institute (MBARI) established the Free Ocean Carbon Enrichment (FOCE) experiment to study the long-term effects of decreased ocean pH levels by developing in-situ platforms [1]. Deep FOCE (dpFOCE) was the first platform, which was deployed in 950 meters of water in Monterey Bay. After the conclusion of dpFOCE, MBARI developed an open source shallow water FOCE (swFOCE) platform located at around 250 meter of water to facilitate worldwide shallow water experiments on FOCE [1][2]. A shallow water platform can be more ubiquitous than a deep-water platform as shallow water instruments are less expensive (as it does not have to be designed to withstand the pressure at deep ocean depths) and more easily deployed (they can be deployed right along the coast). The swFOCE experiment is an open source platform, and MBARI has made the plans available online to anyone interested in studying shallow water carbon enrichment. There is a gateway node what is connected to four sensor nodes within the swFOCE. In order to test the sensor node individually, an idea of designing an Oceanographic Instrument Simulator is purposed. The Oceanographic instrument simulator (OIS), described in this paper provides the means for MBARI engineers to test the swFOCE platform without attaching the numerous and expensive oceanographic instruments. The Oceanographic Instrument Simulator simulates the various scientific instruments that could be deployed in an actual experiment. The Oceanographic Instrument Simulator (OIS) system includes the designed circuit board, Arduino Due and an SD Card shield. The designed circuit board will be connected to a computer through a USB cable, and be connected to MBARI’s swFOCE sensor node through a serial connection. When a query is given from the sensor node, the Arduino Due will parse the data given from the sensor node, search through the pre-installed data in the SD card and return the appropriate data back to the sensor node. A user can also manually set up the input current through a computer terminal window to control the simulated signals from the PCB.

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

Thesis (Ph. D.)--Joint Program in Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 1986. / Vita. / Includes bibliographical references (leaves 208-215). / A study of the remineralization of organic carbon was conducted in the organic-rich sediments of Buzzards Bay, MA. Major processes affecting the carbon chemistry in sediments are reflected by changes in the stable carbon isotope ratios of dissolved inorganic carbon (XCO2) in sediment pore water. Six cores were collected seasonally over a period of two years. The following species were measured in the pore waters: JC02, &'3 C-2CO2, P04, JH 2S, Alk, DOC, and Ca. Measurements of pore water collected seasonally show large gradients with depth, which are larger in summer than in winter. The 613C (PDB) of IC02 varies from 1.3 o/oo in the bottom water to approximately -10 o/oo at 30 cm. During all seasons, there was a trend towards more negative values with depth in the upper 8 cm due to the remineralization of organic matter. There was a trend toward more positive values below 8 cm, most likely due to biological irrigation of sediments with bottom water. Below 16-20 cm, a negative gradient was re-established which indicates a return to remineralization as the main process affecting pore water chemistry. Using the XC02 depth profile, it was estimated that 67-85 gC/m 2 are oxidized annually and 5 gC/m 2-yr are buried. The amount of carbon oxidized represented remineralization occurring within the sediments. This estimate indicated that approximately 20% of the annual primary productivity reached the sediments. The calculated remineralization rates varied seasonally with the high of 7.5 x 10-' mol/L-sec observed in August 84 and the low (0.6 x 10-) in December 83. The calculated remineralization rates were dependent on the amount of irrigation in the sediments; if the irrigation parameter is known to ±20%, then the remineralization rates are known to this certainty also. The amount of irrigation in the sediments was estimated using the results of a seasonal study of 2 2 2Rn/ 22 'Ra disequilibria at the same study site (Martin, 1985). Estimates of the annual remineralization in the sediments using solid-phase data indicated that the solid-phase profiles were not at steady-state concentrations. The isotopic signature of 2C02 was used as an indicator of the processes affecting IC02 in pore water. During every month, the oxidation of organic carbon to C02 provided over half of the carbon added to the IC02 pool. However, in every month, the 6'"C of XCO2 added to the pore water in the surface sediments was greater than -15 o/oo, significantly greater than the 613C of solid-phase organic carbon in the sediments (-20.6 o/oo). The 613C of IC02 added to the pore water in the sediments deeper than 7 cm was between -20 and -21 o/oo, similar to the organic carbon in the sediments. Possible explanations of the 1 3C-enrichment observed in the surface sediments were: a) significant dissolution of CaC0 3 (613C = +1.7 o/oo) b) the addition of significant amounts of carbonate ion from bottom water to pore water c) an isotopic difference between the carbon oxidized in the sediments and that remaining in the sediments. The effect of CaC0 3 dissolution was quantified using measured dissolved Ca profiles and was not large enough to explain the observed isotopic enrichment. An additional source of 13C-enriched carbon was bottom water carbonate ion. In every month studied, there was a net flux of 2C0 2 from pore water to bottom water. The flux of pore water 2C02 to bottom water ranged from a minimum of 10 x 10-12 mol/cm 2 -sec in December 83 to a maximum of 50 x 10-12 mol/cm2-sec in August 84. However, because the pH of bottom water was about 8 while that of the pore water was less than or equal to 7, the relative proportion of the different species of inorganic carbon (H2CO, HCO-, C0~) was very different in bottom water and pore water. Thus, while there was a net flux of IC0 2 from pore water to bottom water, there was a flux of carbonate ion from bottom water to pore water. Because bottom water JC02 was more 13C-enriched than pore water JC0 2, the transfer of bottom water carbonate ion to pore water was a source of 13C-enriched carbon to the pore water. If the &'3C of CO2 added to the pore water from the oxidation of organic carbon was -20.6 o/oo, then the flux of C3% from bottom water to pore water must have been 10-30% of the total flux of 2C02 from pore water to bottom water. This is consistent with the amount calculated from the observed gradient in carbonate ion. Laboratory experiments were conducted to determine whether the 613C of C02 produced from the oxidation of organic carbon (613C-OCOX) was different from the 613C of organic carbon in the sediments (613C-SOC). In the laboratory experiments, mud from the sampling site was incubated at a constant temperature. Three depths were studied (0-3, 10-15, and 20-25 cm). For the first study (IEl), sediment was stirred to homogenize it before packing into centrifuge tubes for incubation. For the second study (IE2), sediment was introduced directly into glass incubation tubes by subcoring. The second procedure greatly reduced disturbance to the sediment. Rates of C02 production were calculated from the concentrations of 2C02 measured over up to 46 days. In both studies, the values of Re in the deeper intervals were about 10% of the surface values. This was consistent with the field results, although the rates decreased more rapidly in the field. In all cases, the remineralization rates during the beginning of IEl were much greater than those at the beginning of IE2. The sediment for IEl was collected in February 84. The measured value of Rc in the surface sediment of the laboratory experiment (24 x 10- mol/L-sec) was much greater than the value of Rc observed in the field in another winter month, December 83 (.62 x 10~9). The sediment for IE2 was collected in August 85. The measured values of Re in the surface sediment (6.6-12 x 10~9 mol/L-sec) were consistent with the field values from August 84 (7.5 x 10-9). The XC02 results indicated that IE2 reproduced field conditions more accurately than IEl did. The isotopic results from the experiments strongly suggested that 613C-OCOX in the surface sediments (-17.8 o/oo ± 1.9 o/oo) was greater than 6'3C-SOC (-20.6 ± 0.2 o/oo). The magnitude of the observed fractionation was small enough that the observed values of 613C-C02 in the pore waters could be explained by fractionated oxidation coupled with the diffusion of carbonate ion from bottom water to pore water. The observed fractionation was most likely due to the multiple sources of organic carbon to coastal sediments. A study of the natural levels of radiocarbon in these sediments indicated that the carbon preserved in the sediments is approximately 30% terrestrial while the rest is from phytoplankton. / by Ann P. McNichol. / Ph.D.

Thesis: Ph. D., Joint Program in Oceanography/Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2015. / Cataloged from PDF version of thesis. / Includes bibliographical references. / This thesis investigates lithospheric dynamics of Earth's subduction zones and Martian tectonic provinces on multiple time scales ranging from short-term earthquake deformation to long-term tectonic loading. In Chapter 2, I use geodetic observations to constrain the postseismic viscoelastic deformation following the 1960 M9.5 Valdivia, Chile earthquake and quantify its stress loading on the rupture zone of the 2010 M8.8 Maule, Chile earthquake. Results of analysis reveal that the post-1960 viscoelastic process might have contributed to the triggering of the 2010 earthquake. Chapter 3 presents numerical experiments to investigate elastoplastic deformation and faulting in the overriding plates of subduction zones caused by the movement of subducted seamounts. Numerical simulations show that a group of normal faults would first appear on the seaward side of a subducted seamount, followed by a group of thrust faults on the landward side of the seamount. In Chapter 4, I use the most recent Martian gravity and topography data to constrain spatial variations in lithospheric flexural deformation for various tectonic regions on Mars. The effective lithospheric thickness is estimated to be relatively small for the plain regions in the southern highland, but relatively large for the impact basins in the northern lowland as well as for volcanic montes in the Tharis province. The regional variations in the estimated effective lithospheric thickness might reflect both spatial and temporal changes in the thermal state of Mars. / by Min Ding. / Ph. D.