Thesis (S.M.)--Joint Program in Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Dept. of Ocean Engineering; and, the Woods Hole Oceanographic Institution), 2003. / Includes bibliographical references (p. 118-120). / by Christopher George Bohner. / S.M.

Thesis (Ph. D.)--Joint Program in Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Department of Mechanical Engineering; and the Woods Hole Oceanographic Institution), 2015. / This electronic version was submitted by the student author. The certified thesis is available in the Institute Archives and Special Collections. / Cataloged from student-submitted PDF version of thesis. / Includes bibliographical references (pages 181-201). / Real-time cooperation between autonomous vehicles can enable time-critical missions such as tracking and pursuit of a dynamic target or environmental feature, but relies on wireless communications. Underwater, communication over distances beyond about one hundred meters is almost exclusively accomplished through acoustics, which bring challenges such as propagation delays, low data rates, packet loss, and scheduling constraints due to interference and limited bandwidth. These limitations make underwater pursuit missions preeminent applications of networked control. Motivated by such applications, this thesis presents contributions towards multi-vehicle feedback control in the presence of severe communication constraints. The first major area of work considers the formulation and solution of new underwater multi-vehicle tracking and pursuit problems using closed-loop control. We begin with a centralized robust optimization approach for multicast routing and power control which is suitable for integration with vehicle control. Next, we describe field experiments in range-based target pursuit at high tracking bandwidths in a challenging shallow-water environment. Finally, we present a methodology for pursuit of dynamic ocean features such as fronts, which we validate using hindcast ocean model data. The primary innovation is a projection algorithm which carries out linearization of ocean model forecast dynamics and uncertainty directly in vehicle coordinates via a forward model technique. The resulting coupled linear stochastic system is suitable for networked control. The second area of work presents a unified formalism for multi-vehicle control and estimation with measurement, control, and acknowledgment packets all subject to scheduling, delays and packet loss. The modular framework we develop is built around a jump linear system description incorporating receding horizon optimization and buffering at actuators. Integration of these elements enables synthesis of a novel technique for estimation using delayed and lossy control acknowledgments-a desirable and practical capability of fielded systems that has not been considered to date. Simulations and field experiments demonstrate the effectiveness of our approach. / by Brooks Louis-Kiguchi Reed. / Ph.D.

Thesis (Ph. D.)--Joint Program in Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2001. / "September 2001." Vita. Page 223 blank. / Includes bibliographical references. / This thesis studies interactions between mid-ocean ridges and mantle plumes using geophysics, geochemistry, and geodynamical modeling. Chapter 1 investigates the effects of the Marion and Bouvet hotspots on the ultra-slow spreading, highly-segmented Southwest Indian Ridge (SWIR). Gravity data indicate that both Marion and Bouvet impart high-amplitude mantle Bouguer anomaly lows to the ridge axis, and suggest that long-offset transforms may diminish along-axis plume flow. Building upon this observation, Chapter 2 presents a series of 3D numerical models designed to quantify the sensitivity of along-axis plume-driven mantle flow to transform offset length, spreading rate, and mantle viscosity structure. The calculations illustrate that long-offset transforms in ultra-slow spreading environments may significantly curtail plume dispersion. Chapter 3 investigates helium isotope systematics along the western SWIR as well as near a global array of hotspots. The first part of this study reports uniformly low 3He/4He ratios of 6.3-7.3 R/Ra along the SWIR from 9⁰-24⁰E, compared to values of 8 +/- 1 Ra for normal mid-ocean ridge basalt. The favored explanation for these low values is addition of (U+Th) into the mantle source by crustal and/or lithospheric recycling. Although high He/4He values have been observed along the SWIR near Bouvet Island to the west, there is no evidence for elevated 3He/4He ratios along this section of the SWIR. The second part of Chapter 3 investigates the relationship between 3He/4He ratios and geophysical indicators of plume robustness for nine hotspots. / (cont.) A close correlation between a plume's flux and maximum 3He/4He ratio suggests a link between plume upwelling strength and origination in the deep, relatively undegassed mantle. Chapter 4 studies 3D mantle flow and temperature patterns beneath oceanic ridge-ridge-ridge triple junctions (TJs). In non-hotspot-affected TJs with geometry similar to the Rodrigues TJ, temperature and upwelling velocity along the slowest-spreading of the three ridges are predicted to increase within a few hundred kilometers of the TJ, to approach those of the fastest-spreading ridge. Along the slowest-spreading branch in hotspot-affected TJs such as the Azores, a strong component of along-axis flow directed away from the TJ is predicted to advect a hotspot thermal anomaly away from its deep-seated source. / by Jennifer E. Georgen. / Ph.D.

Thesis (Ph. D.)--Joint Program in Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2006. / Includes bibliographical references. / Between 90-25° E on the ultraslow-spreading Southwest Indian Ridge lie two sharply contrasting supersegments. One 630 km long supersegment erupts N-MORB that is progressively enriched in incompatible element concentrations from east to west. The second 400 km long supersegment contains three separate volcanic centers erupting E-MORB and connected by long amagmatic accretionary segments, where mantle is emplaced directly to the seafloor with only scattered N-MORB and E-MORB erupted. Rather than a major break in mantle composition at the discontinuity between the supersegments, this sharp contrast in geometry, physiography, and chemistry reflects "source" versus "process" dominated generation of basalt. Robust along-axis correlation of ridge characteristics (i.e. morphology, upwelling rate, lithospheric thickness), basalt chemistry, and crustal thickness (estimated from gravity) provides a unique opportunity to compare the influence of spreading geometry and rate on MORB generation. What had not been well established until now is the importance of melting processes rather than source at spreading rates < 20 mm/yr. / (cont.) Along the orthogonally spreading supersegment (14 mm/yr) moderate degrees of partial melting effectively sample the bulk mantle source, while on the obliquely spreading supersegment (7-14 mm/yr) suppression of mantle melting to low degrees means that the bulk source is not uniformly sampled, and thus "process" rather than "source" dominates melt chemistry. / by Jared Jeffrey Standish. / Ph.D.

Thesis (S.M.)--Joint Program in Marine Geology and Geophysics (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2008. / Includes bibliographical references (p. 30-37). / We evaluate the applicability of plagioclase and gabbro flow laws by comparing predicted and observed deformation mechanisms in gabbroic shear zones. Gabbros and layered gabbro mylonites were collected from the Southwest Indian Ridge (SWIR), ODP Hole 735B. Deformation temperatures are constrained by two-pyroxene thermometry, stress is estimated from grain size, and deformation mechanisms are analyzed by microstructure and the presence or absence of a lattice preferred orientation (LPO). Our analyses indicate that mylonite layers deformed at a strain rate in the range of 1012 to 101 s-1, while coarse-grained gabbro deformed at a strain rate of approximately 10-14 to 1013 s-1. Plagioclase in pure plagioclase mylonite layers exhibit strong LPOs indicating they deform by dislocation creep. Plagioclase grain size in mixed plagioclase-pyroxene mylonite layers is finer than in pure plagioclase layers, and depends on the size and proportion of pyroxenes. Progressive mixing of pyroxene and plagioclase within gabbro mylonite layers is accompanied by weakening of the LPO indicating that phase mixing promotes a transition to diffusion creep processes that involve grain boundary sliding. Our results indicate that experimental flow laws are accurate at geologic strain rates, although the strain rate for diffusion creep of fine-grained gabbro may be underestimated. At the conditions estimated for the SWIR crust, our calculations suggest that strain localization leads to a factor of two to four decrease in lower crustal viscosity. Even so, the viscosity of lower gabbroic crust is predicted to be similar to that of dry upper mantle. / by Luc Mehl. / S.M.

Thesis (S.M.)--Joint Program in Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2008. / Includes bibliographical references (p. 45-47). / Biogeochemical cycling of phosphate is a key component in the overall production rate of coastal ecosystems. Mineral phases in the near-shore sediments play a significant role in the return of phosphate remineralized in the upper sediments to the water column. Sequential Extraction (SEDEX) of the solid-phase associated P04-3 yielded reservoir profiles of phosphate at three sites off of the Massachusetts coast. These extractions found Fe-associated P04 to be the dominant phase associated with rapid porewater-solid P exchange. Additionally, a seasonal enrichment/depletion pattern of phosphate fluxes relative to total carbon was observed from the sediments. These observations established the behavior of phosphate in coastal sediments as interconnected with the ongoing Fe-cycling in the sediments as well. / by Richard Hunter Oates, Jr. / S.M.

Thesis (S.M.)--Joint Program in Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences and the Woods Hole Oceanographic Institution), 2002. / Includes bibliographical references (p. 22-28). / We have investigated the near liquidus phase relations of a primitive absarokite from the Mascota region in Western Mexico. Sample M. 102 was chosen because it has high MgO contents, a high Mg# and Fo90 olivine phenocrysts, indicating it is primitive mantle melt. Highpressure experiments on a synthetic analogue of the absarokite composition with a H20 content of either -1.7 wt% or -5.1 wt% were carried out in a piston cylinder apparatus. The composition with -1.7 wt% H20 is multiply saturated with olivine and orthopyroxene as liquidus phases at 1.6 GPa and 14000C. At the same pressure clinopyroxene appears 300C below the liquidus. With a H20 content of -5.1 wt% composition M.102 is multiply saturated with olivine and orthopyroxene on the liquidus at 1.7 GPa and 13000 C. Assuming batch melting, we suggest that absarokite M. 102 segregated from a depleted lherzolite or harzburgite residue at depth -50 km depth in the mantle wedge. Unlike most lavas in the region, the absarokite has not ponded and fractionated at the crust mantle interface (-35-40 km), and the temperatures of multiple saturation indicate that the mantle wedge beneath the Jalisco block is hotter than previously thought. The low degree batch melting of an original metasomatised harzburgite source, can produce the observed trace element abundances. The liquidus phase relations are not consistent with the presence of non-peridotitic veins at the depth of last equilibration. Therefore, we propose that the Mascota absarokites segregated at an apparent melt fraction of less than 5% from a depleted peridotitic source. They initially formed by a small degree of melting of a metasomatised original source at greater depth. / by Marc Hesse. / S.M.

Thesis (M.S.)--Joint Program in Oceanography (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 1987. / Supervised by David G. Aubrey. / Includes bibliographical references (leaves 57-61). / by Barbara Vanston Braatz. / M.S.

Thesis (Ph. D.)--Joint Program in Marine Geology and Geophysics (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 1987. / Includes bibliographical references (v.2, leaves 309-317). / Many features of the oceanic plates cannot be explained by conductive cooling with age. A number of these anomalies require additional convective thermal sources at depths below the plate: mid-plate swells, the evolution of fracture zones, the mean depth and heat flow relationships with age and the observation of small scale (150-250 km) geoid and topography anomalies in the Central Pacific and Indian oceans. Convective models are presented of the formation and evolution of these features. In particular, the effect of a shallow low viscosity layer in the uppermost mantle on mantle flow and its geoid, topography, gravity and heat flow expression is explored. A simple numerical model is employed of convection in a fluid which has a low viscosity layer lying between a rigid bed and a constant viscosity region. Finite element calculations have been used to determine the effects of (1) the viscosity contrast between the two fluid layers, (2) the thickness of the low viscosity zone, (3) the thickness of the conducting lid, and (4) the Rayleigh number of the fluid based on the viscosity of the lower layer. A model simple for mid-plate swells is that they are the surface expression of a convection cell driven by a heat flux from below. The low viscosity zone causes the top boundary layer of the convection cell to thin and, at high viscosity contrasts and Rayleigh numbers, it can cause the boundary layer to go unstable. The low viscosity zone also mitigates the transmission of normal stress to the conducting lid so that the topography and geoid anomalies decrease. The geoid anomaly decreases faster than the topography anomaly, however, so that the depth of compensation can appear to be well within the conducting lid. Because the boundary layer is thinned, the elastic plate thickness also decreases and, since the low viscosity allows the fluid to flow faster in the top layer, the uplift time decreases as well. We have compared the results of this modeling to data at the Hawaii, Bermuda, Cape Verde and Marquesas swells, and have found that it can reproduce their observed anomalies. The viscosity contrasts that are required range from 0.2-0.01, which are in agreement with other estimates of shallow viscosity variation in the upper mantle. Also, the estimated viscosity contrast decreases as the age of the swell increases. This trend is consistent with theoretical estimates of the variation of such a low viscosity zone with age. Fracture zones juxtapose segments of the oceanic plates of different ages and thermal structures. The flow induced by the horizontal temperature gradient at the fracture zone initially downwells immediately adjacent to the fracture zone on the older side, generating cells on either side of the plume. The time scale and characteristic wavelength of this flow depends initially on the viscosity near the largest temperature gradient in the fluid which, in our model, is the viscosity of the low viscosity layer. They therefore depend on both the Rayleigh number and the viscosity contrast between the layers. Eventually the flow extends throughout the box, and the time scales and the characteristic wavelengths of the flow depend on the thickness and viscosity of both layers. When the Rayleigh number based on the viscosity of the top la er, and the depth of both fluid layers, is less than 10 , the geoid anomalies of these flows are dominated by the convective signal. When this Rayleigh number exceeds 106, the geoid anomalies retain a step across the fracture zone out to large ages. We have compared our results to geoid anomalies over the Udintsev fracture zone, and have found that the predicted geoid anomalies, with high effective Rayleigh numbers, agree at longer wavelengths with the observed anomalies and can produce the observed geoid slope-age behaviour. We have also compared the calculated topographic steps to those predicted by the average depth-age relationships observed in the oceans. We have found that only with a low viscosity zone will the flow due to fracture zones not disturb the average depth versus age relationships. We have also applied the model to a numerical study of the effect of a low viscosity zone in the uppermost mantle on the onset and surface expression of convective instabilities in the cooling oceanic plates. We find that the onset and magnitude of the geoid, topography and heat flow anomalies produced by these instabilities are very sensitive to the viscosity contrast and the Rayleigh number, and that the thickness of the low viscosity zone is constrained by the wavelength of the observables. If the Rayleigh number of the low viscosity zone exceeds a critical value then the convection will be confined to the low viscosity zone for a period which depends on the viscosity contrast and the Rayleigh number. The small scale convection will eventually decay into longer wavelength convection which extends throughout the upper mantle, so that the small scale convective signal will eventually be succeeded by a longer wavelength signal. We compare our model to the small scale geoid and topography anomalies observed in the Southeast Pacific. The magnitude (0.50-0.80 m in geoid and 250 m in topography), early onset time (5-10 m.y.) and lifetime (over 40 m.y.) of these anomalies suggest a large viscosity contrast of greater than two orders of magnitude. The trend to longer wavelengths also suggests a high Rayleigh number of near or over 10 and their original 150-250 km wavelength indicates a low viscosity zone of 75- 125 km thickness. We have found that the presence of such small scale convection does not disturb the slope of the depth-age curve but elevates it by up to 250 m, and it is not until the onset of long wavelength convection that the depth-age curves radically depart from a cooling halfspace model. In the Pacific, the depth-age curve is slightly elevated in the region where small scale convection is observed and it does not depart from a halfspace cooling model until an age of 70 m.y.. Models that produce the small scale anomalies predict a departure time between 55 and 65 m.y. / by Elizabeth M. Robinson. / Ph.D.

Thesis (Ph. D.)--Joint Program in Chemical Oceanography (Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2013. / Cataloged from PDF version of thesis. / Includes bibliographical references (p. 30-38). / Thiol substances can form stable complexes with metals (especially copper and mercury) in the surface ocean that can impact cycling and bioavailability of those elements. In this study, I present seven concentration profiles of cysteine and glutathione, two low-molecular weight thiols, from the coastal northwest Atlantic Ocean and the Bermuda Atlantic Time Series (BATS) sampling site in the Sargasso Sea, a first for these regions. These two thiols were found in the upper 200 meters of the ocean at all sites, and the total thiol concentration varied from 0.2 to 3.2 nM. The highest concentration of both thiols was found at the deep chlorophyll maximum in most samples. Thiol concentrations were higher on the continental shelf than in the open ocean. The observed distribution of cysteine and glutathione and thermodynamic stability of copper complexes suggests that Cu(I)-dithiol complexes may be the dominant surface ocean copper and thiol species. Mercury-thiol complexes were also present in thermodynamically modeled seawater, which may provide a vector for mercury uptake in the surface ocean. / by Tristan Kading. / Ph.D.