Thesis (Ph. D.)--Joint Program in Oceanography/Applied Ocean Science and Engineering (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. 191-202). / Melting and crystallization processes on the Earth and Moon are explored in this thesis, and the topics of melt generation, transport, and crystallization are discussed in three distinct geologic environments: the Moon's mantle, the Greenland ice sheet, and the Earth's crust. Experiments have been conducted to determine the conditions of origin for two high-titanium magmas from the Moon. The lunar experiments (Chapter 2) were designed to explore the effects of variable oxygen fugacity (fo₂) on the high pressure and high temperature crystallization of olivine and orthopyroxene in high-Ti magmas. The results of these experiments showed that the source regions for the high-Ti lunar magmas are distributed both laterally and vertically within the lunar mantle, and that it is critical to estimate the pre-eruptive oxygen fugacity in order to determine true depth of origin for these magmas within the lunar mantle. Chapter 3 models the behavior of water flow through the Greenland ice sheet driven by hydrofracture of water through ice. The results show that melt water in the ablation zone of Greenland has almost immediate access to the base of the ice sheet in areas with up two kilometers of ice. Chapter 4 is an experimental study of two hydrous high-silica mantle melts from the Mt. Shasta, CA region. Crystallization is simulated at H₂O saturated conditions at all crustal depths, and a new geobarometer-hygrometer based on amphibole magnesium number is calibrated. In Chapter 5 I use the new barometer to study a suite of mafic enclaves from the Mt. Shasta region, and apply it to amphiboles in these enclaves. Evidence for pre-eruptive H₂O contents of up to 14 wt% is presented, and bulk chemical analyses of the inclusions are used to show that extensive magma mixing has occurred at all crustal depths up to 35 km beneath Mt. Shasta. / by Michael James Krawczynski. / Ph.D.

Thesis: Ph. D., Joint Program in Oceanography (Massachusetts Institute of Technology, Department of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2014. / 122 / Cataloged from PDF version of thesis. / Includes bibliographical references (pages 149-158). / Understanding the fate of primary production in the ocean is a challenging task because once produced, organic material is oxidized over timescales which range from minutes, to millions of years. This timescale diversity is matched by an equal heterogeneity in both the local physical and chemical environment. In this thesis we explore the relationship between the distinct reservoirs of organic carbon in the ocean and their underlying complexity. First, we show how the heterogeneity of portions of the carbon cycle can be packaged in terms of age structured models and their accompanying age and rate distributions. We further relate the moments of the rate distributions to bulk reservoir properties like average age and flux. Explicit relationships are then derived for the specific case of a single turnover time and a lognormal distribution. We apply these ideas to the problem of dissolved organic carbon (DOC) cycling in the ocean. Current models of bulk concentration and isotope data suggest a microbially sourced DOC reservoir consisting of two components. A nearly homogeneous background component with a long turnover time (> 6000 years) is joined by a component of fast turnover time (~ 1 year) and equal concentration in the surface ocean. We confirm the presence of isotopically enriched, modern DOC co-cycling with an isotopically depleted older fraction in the upper ocean. However, our results show that up to 30% of the deep DOC reservoir is modern and supported by a 1 Pg per year carbon flux, ten times higher than inferred from bulk isotope measurements. Isotopically depleted material turns over at an apparent time scale of 30, 000 years, far slower than indicated by bulk isotope measurements. These results are consistent with global DOC measurements and explain both the fluctuations in deep DOC concentration and the anomalous radiocarbon values of DOC in the Southern Ocean. Finally, the thesis explores methods for determining the validity of diffusion limitation as the mechanism behind the power-law slowdown in organic remineralization in sediment. We find that diffusion limitation connects the decay behavior of organic material to the correlations found between mineral surface area and organic matter content in sediments. / by Christopher L. Follett. / 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 (pages 173-186). / This thesis investigates the long-term evolution of rift-bounding normal faults. To first order, the observed diversity of extensional tectonic styles reflects differences in the maximum offset that can be accommodated on individual faults during their life span. My main objective is to develop a theoretical framework that explains these differences in terms of a few key mechanical and geological controls. I start by laying out the energy cost associated with slip on a normal fault, which consists of (1) overcoming the frictional resistance on the fault, (2) bending the faulted layer and (3) sustaining the growth of topography. In Chapter 2, I propose that flexural rotation of the active fault plane enables faults to evolve along a path of minimal energy, thereby enhancing their life span. Flexural rotation occurs more rapidly in thinner faulted layers, and can potentially explain the wide range of normal fault dips documented with focal mechanisms. In Chapter 3, I show that surface processes can enhance the life span of continental normal faults by reducing the energy cost associated with topography buildup. In Chapter 4, I focus on lithospheric bending induced by fault growth, which is well described by elasto-plastic flexure models. I demonstrate that numerical models that treat the lithosphere as a visco-plastic solid can properly predict fault evolution only when the rate-dependent viscous flexural wavelength of the lithosphere is accommodated within the numerical domain. In Chapter 5, I consider the interplay of faulting and crustal emplacement at a slow mid-ocean ridge. I show that a depth-variable rate of magma emplacement can reconcile the formation of long-lived detachment faults, which requires a moderate melt supply, and the exhumation of large volumes of lower crustal material. Finally, in Chapter 6 I investigate the three-dimensional interactions between normal faults in a lithosphere of varying thickness. I suggest that large along-axis gradients in lithospheric thickness can prevent the growth of continuous faults along-axis, and instead decouple the modes of faulting at the segment center and at the segment end. / by Jean-Arthur Louis Olive. / Ph. D.

Thesis (S.M.)--Joint Program in Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Dept. of Electrical Engineering and Computer Science; and the Woods Hole Oceanographic Institution), 2005. / Includes bibliographical references (leaves 120-121). / This thesis develops and utilizes a method for analyzing data from the North Pacific Acoustic Laboratory's (NPAL) Basin Acoustic Seamount Scattering Experiment (BASSEX). BASSEX was designed to provide data to support the development of analytical techniques and methods which improve the understanding of sound propagation around underwater seamounts. The depth-dependent sound velocity profile of typical ocean waveguides force sound to travel in convergence zones about a minimum sound speed depth. This ducted nature of the ocean makes modeling the acoustic field around seamounts particularly challenging, compared to an isovelocity medium. The conical shape of seamounts also adds to the complexity of the scatter field. It is important to the U.S. Navy to understand how sound is diffracted around this type of topographic feature. Underwater seamounts can be used to conceal submarines by absorbing and scattering the sound they emit. BASSEX measurements have characterized the size and shape of the forward scatter field around the Kermit-Roosevelt Seamount in the Pacific Ocean. Kermit-Roosevelt is a large, conical seamount which shoals close to the minimum sound speed depth, making it ideal for study. Acoustic sources, including M-sequence and linear frequency-modulated sources, were stationed around the seamount at megameter ranges. A hydrophone array was towed around the seamount to locations which allowed measurement of the perturbation zone. Results from the method developed in this thesis show that the size and shape of the perturbation zone measured coincides with theoretical and experimental results derived in previous work. / by Joseph J. Sikora, III. / S.M.

Thesis (Ph. D.)--Joint Program in Oceanography/Applied Ocean Science and Engineering (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. 199-205). / In a stratified rotating fluid, frictionally driven circulations couple with the buoyancy field over sloping topography. Analytical and numerical methods are used to quantify the impact of this coupling on the vertical circulation, spindown of geostrophic flows, and the formation of a shelfbreak jet. Over a stratified. slope, linear spindown of a geostrophic along-isobath flow induces cross-isobath Ekman flows. Ekman advection of buoyancy weakens the vertical circulation and slows spindown. Upslope (downslope) Ekman flows tend to inject (remove) potential vorticity into (from) the ocean. Momentum advection and nonlinear buoyancy advection are examined in setting asymmetries in the vertical circulation and the vertical relative vorticity field. During nonlinear homogeneous spindown over a flat bottom, momentum advection weakens Ekman pumping and strengthens Ekman suction, while cyclonic vorticity decays faster than anticyclonic vorticity. During nonlinear stratified spindown over a slope, nonlinear advection of buoyancy enhances the asymmetry in Ekman pumping and suction, whereas anticyclonic vorticity can decay faster than cyclonic vorticity outside of the boundary layers. During the adjustment of a spatially uniform geostrophic current over a shelfbreak, coupling between the Ekman flow and the buoyancy field generates Ekman pumping near the shelfbreak, which leads to the formation of a jet. Scalings are presented for the upwelling strength, the length scale over which it occurs, and the timescale for jet formation. The results are applied to the Middle Atlantic Bight shelfbreak. / by Jessica A. Benthuysen. / 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), 2007. / Includes bibliographical references. / The origin of symmetric alternating magnetic polarity stripes on the seafloor is investigated in two marine environments; along the ridge axis of the fast spreading East Pacific Rise (EPR) (90 25'-90 55'N) and at Kane Megamullion (KMM) (230 40'N), near the intersection of the slow-spreading Mid Atlantic Ridge with Kane Transform Fault. Marine magnetic anomalies and magnetic properties of seafloor samples are combined to characterize the magnetic source layer in both locations. The EPR study suggests that along-axis variations in the observed axial magnetic anomaly result from changing source layer thickness alone, consistent with observed changes in seismic Layer 2a. The extrusive basalts of the upper crust therefore constitute the magnetic source layer along the ridge axis and long term crustal accretion patterns are reflected in the appearance of the axial anomaly. At KMM the C2r.2r/C2An. In (- 2.581 Ma) polarity reversal boundary cuts through lower crust (gabbro) and upper mantle (serpentinized peridotites) rocks exposed by a detachment fault on the seafloor, indicating that these lithologies can systematically record a magnetic signal. Both lithologies have stable remanent magnetization, capable of contributing to the magnetic source layer. The geometry of the polarity boundary changes from the northern to the central regions of KMM and is believed to be related to changing lithology. In the northern region, interpreted to be a gabbro pluton, the boundary dips away from the ridge axis and is consistent with a rotated conductively cooled isotherm. In the central region the gabbros have been removed and the polarity boundary, which resides in serpentinized peridotite, dips towards the ridge axis and is thought to represent an alteration front. The linear appearance of the polarity boundary across both regions indicates that the two lithologies acquired their magnetic remanence during approximately the same time interval. Seismic events caused by detachment faulting at Kane and Atlantis Transform Faults are investigated using hydroacoustic waves (T-phases) recorded by a hydrophone array. Observations and ray trace models of event propagation show bathymetric blockage along propagation paths, but suggest current models of T-phase excitation and propagation need to be improved to explain observed characteristics of T-phase data. / by Clare Margaret Williams. / Ph.D.

Thesis (S. M.)--Joint Program in Oceanography/Applied Ocean Science and Engineering (Massachusetts Institute of Technology, Dept. of Earth, Atmospheric, and Planetary Sciences; and the Woods Hole Oceanographic Institution), 2009. / 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 (p. 41-43). / Mechanisms for the tidal component of salt flux in the Hudson River estuary are investigated using a 3D numerical model. Variations with river discharge, fortnightly tidal forcing, and along channel variability are explored. Four river discharge conditions were considered: 1200 m3 s-1, 600 m3 s-1, 300 m3 s-1, 150 m3 s-1. Tide-induced residual salt flux was found to be variable along the channel, with locations of counter-gradient flux during both neap and spring tide. The magnitude of tidal salt flux scales with the river flow and has no clear dependence on the spring-neap tidal forcing. The diffusive fraction, ?, has a value of -0.25 to 0.46 in the lower estuary, increasing to -0.23 to 1 near the head of salt. The phase lag between tidal salinity and velocity is analyzed at three cross-sections with: large positive, negative, and weak tidal salt flux. The composite Froude number, G2, is calculated along the channel and indicates nearly ubiquitous supercritical flow for maximum flood and ebb during both neap and spring tides. Subcritical flow occurs during slack water and at geographically locked locations during neap floods. Application of two-layer, frictional hydraulic theory reveals how variations in channel width and depth generate tidal asymmetries in cross-sectional salinity, the key ingredient of tidal salt flux. / by Patricia Ann Engel. / S.M.

Thesis (Ph. D.)--Joint Program in Oceanography/Applied Ocean Science and Engineering (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. / Large earthquake rupture and triggering mechanisms that drive seismicity in subduction zones are investigated in this thesis using a combination of earthquake observations, statistical and physical modeling. A comparison of the rupture characteristics of M 7.5 earthquakes with fore-arc geological structure suggests that long-lived frictional heterogeneities (asperities) are primary controls on the rupture extent of large earthquakes. To determine when and where stress is accumulating on the megathrust that could cause one of these asperities to rupture, this thesis develops a new method to invert earthquake catalogs to detect space-time variations in stressing rate. This algorithm is based on observations that strain transients due to aseismic processes such as fluid flow, slow slip, and afters lip trigger seismicity, often in the form of earthquake swarms. These swarms are modeled with two common approaches for investigating time-dependent driving mechanisms in earthquake catalogs: the stochastic Epidemic Type Aftershock Sequence model [Ogata, 1988] and the physically-based rate-state friction model [Dieterich, 1994]. These approaches are combined into a single model that accounts for both aftershock activity and variations in background seismicity rate due to aseismic processes, which is then implemented in a data assimilation algorithm to invert catalogs for space-time variations in stressing rate. The technique is evaluated with a synthetic test and applied to catalogs from the Salton Trough in southern California and the Hokkaido corner in northeastern Japan. The results demonstrate that the algorithm can successfully identify aseismic transients in a multi-decade earthquake catalog, and may also ultimately be useful for mapping spatial variations in frictional conditions on the plate interface. / by Andrea Lesley Llenos. / 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), 2017. / Cataloged from PDF version of thesis. / Includes bibliographical references (pages 104-115). / This thesis investigates the evolution of the oceanic lithosphere in a broad sense from formation to subduction, in a focused case at the ridge, and in a focused case proximal to subduction. In general, alteration of the oceanic lithosphere begins at the ridge through focused and diffuse hydrothermal flow, continues off axis through low temperature circulation, and may occur approaching subduction zones as bending related faulting provides fluid pathways. In Chapter 2 1 use a dataset of thousands of microearthquakes recorded at the Rainbow massif on the Mid-Atlantic Ridge to characterize the processes which are responsible for the long-term, high-temperature, hydrothermal discharge found hosted in this oceanic core complex. I find that the detachment fault responsible for the uplift of the massif is inactive and that the axial valleys show no evidence for faulting or active magma intrusion. I conclude that the continuous, low-magnitude seismicity located in diffuse pattern in a region with seismic velocities indicating ultramafic host rock suggests that serpentinization may play a role in microearthquake generation but the seismic network was not capable of providing robust focal mechanism solutions to constrain the source characteristics. In Chapter 3 I find that the Juan de Fuca plate, which represents the young/hot end-member of oceanic plates, is lightly hydrated at upper crustal levels except in regions affected by propagator wakes where hydration of lower crust and upper mantle is evident. I conclude that at the subduction zone the plate is nearly dry at upper mantle levels with the majority of water contained in the crust. Finally, in Chapter 4 I examine samples of cretaceous age serpentinite sampled just before subduction at the Puerto Rico Trench. I show that these upper mantle rocks were completely serpentinized under static conditions at the Mid-Atlantic Ridge. Further, they subsequently underwent 100 Ma of seafloor weathering wherein the alteration products of serpentinization themselves continue to be altered. I conclude that complete hydration of the upper mantle is not the end point in the evolution of oceanic lithosphere as it spreads from the axis to subduction. / by Gregory Horning. / 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. / Page 250 blank. / Includes bibliographical references. / Half of the ocean crust is formed at spreading centers with total opening rates less than 40 km/Myr. The objective of this Thesis is to investigate temporal variations in active ridge processes and crustal aging at slow-spreading centers by comparing axial crustal structure with that on conjugate flanks of the slow-spreading Mid-Atlantic Ridge (MAR) (full rate, 20 km/Myr) and the ultra-slow spreading Southwest Indian Ridge (SWIR) (full rate, 14 km/Myr). Seismic refraction data collected along the rift valley and flanking rift mountains of the OH-1 segment (35ʻN) at the MAR show that the entire crustal section is constructed within a zone that is less than 5 km wide. Shallow-level hydrothermal circulation within the axial valley is suggested by the rift mountain seismic profiles, which show that the upper crust is 20% thinner and 16% faster along strike than zero-age crust. These effects probably result from fissure sealing within the extrusive crust. Deeper crustal velocities remain relatively constant at the segment midpoint within the first 2 Myr, but are reduced near the segment offsets presumably by faulting and fracturing associated with uplift out of the rift valley. / (cont.) A temporal variation in axial melt supply is suggested by a 15% difference in along-strike crustal thickness between the rift valley and rift mountains, with relatively less melt supplied today than 2 Ma. Crustal accretion at the SWIR appears to occur in a similar manner as at the MAR, although gravity and seismic data indicate that the average crustal thickness is 2-4 km less at theultra-slow spreading SWIR. A 25 Myr record on both flanks of the ridge shows that seafloor spreading has been highly asymmetric through time, with 35% faster crustal accretion on the Antarctic (south) plate. A small-offset non-transform discontinuity between two ridge segments is just as stable as two neighboring transform discontinuities, although a single mantle Bouguer gravity anomaly centered over the non-transform offset indicates that this boundary does not significantly perturb underlying mantle flow. Off-axis magnetic anomalies are recorded with high fidelity despite the very low spreading rates and the absence of a basaltic upper crust in one area. The lower crust may be the dominant off-axis carrier of the magnetic signal, contrary to traditionalmodels of crustal magnetic structure. Morphological and gravity data show evidence of asymmetric crustal accretion across the SWIR ridge axis, with slightly warmer mantletemperatures beneath the slower-spreading African (north) plate. / by Allegra Hosford. / Ph.D.