ALBERT

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Gosselin, K. M., Nelson, R. K., Spivak, A. C., Sylva, S. P., Van Mooy, B. A. S., Aeppli, C., Sharpless, C. M., O’Neil, G. W., Arrington, E. C., Reddy, C. M., & Valentine, D. L. Production of two highly abundant 2-methyl-branched fatty acids by blooms of the globally significant marine cyanobacteria Trichodesmium erythraeum. ACS Omega, 6(35), (2021): 22803–22810, https://doi.org/10.1021/acsomega.1c03196.

Description: The bloom-forming cyanobacteria Trichodesmium contribute up to 30% to the total fixed nitrogen in the global oceans and thereby drive substantial productivity. On an expedition in the Gulf of Mexico, we observed and sampled surface slicks, some of which included dense blooms of Trichodesmium erythraeum. These bloom samples contained abundant and atypical free fatty acids, identified here as 2-methyldecanoic acid and 2-methyldodecanoic acid. The high abundance and unusual branching pattern of these compounds suggest that they may play a specific role in this globally important organism.

Description: This work was funded with grants from the National Science Foundation grants OCE-1333148, OCE-1333162, and OCE-1756254 and the Woods Hole Oceanographic Institution (IR&D). GCxGC analysis made possible by WHOI’s Investment in Science Fund.

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

Description: Marine microbes are key drivers of biogeochemical transformations within the world’s oceans. Although seawater appears uniform at scales that humans often interact with and sample, the world that marine microbes inhabit can be highly heterogeneous, with numerous biological and physical processes giving rise to resource hotspots where nutrient concentrations exceed background levels by orders of magnitude. While the impact of this microscale heterogeneity has been investigated in the laboratory with microbial isolates and theoretical models, microbial ecologists have lacked adequate tools to interrogate microscale processes directly in the natural environment. Within this thesis I introduce three new technologies that enable interrogation of microbial processes at the microscale in natural marine communities. The IFCB-Sorter acquires images and sorts individual phytoplankton cells, directly from seawater, allowing studies exploring connections between the diversity of forms present in the plankton and genetic variability at the single-cell level. The In Situ Chemotaxis Assay (ISCA) is a field-going microfluidic device designed to probe the distribution and role of motility behavior among microbes in aquatic environments. By creating microscale hotspots that simulate naturally occurring ones, the ISCA makes it possible to examine the role of microbial chemotaxis in resource acquisition, phytoplankton-bacteria interactions, and host-symbiont systems. Finally, the Millifluidic In Situ Enrichment (MISE) is an instrument that enables the study of rapid shifts in gene expression that permit microbial communities to exploit chemical hotspots in the ocean. The MISE subjects natural microbial communities to a chemical amendment and preserves their RNA in a minute-scale time series. Leveraging an array of milliliter-volume wells, the MISE allows comparison of community gene expression in response to a chemical stimulus to that of a control, enabling elucidation of the strategies employed by marine microbes to survive and thrive in fluctuating environments. Together, this suite of instruments enables culture-independent examination of microbial life at the microscale and will empower microbial ecologists to develop a more holistic understanding of how interactions at the scale of individual microbes impact processes in marine ecosystems at a global scale.

Description: I’d like to thank the Gordon and Betty Moore Foundation, the National Science Foundation, and NSERC for funding portions of my research.

Description: I’d like to thank the Gordon and Betty Moore Foundation, the National Science Foundation, and NSERC for funding portions of my research.

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 2014

Description: Viral predation on bacteria in the ocean liberates carbon from the particulate fraction, where it is accessible to higher trophic levels, and redirects it to the dissolved fraction, where it supports microbial growth. Although viruses are highly abundant in the ocean little is known about how their interactions with bacteria are structured. This challenge arises because the diversity of both bacteria and viruses is exceedingly high and interactions between them are mediated by specific molecular interactions. This thesis uses heterotrophic bacteria of the genus Vibrio as a model to quantify virus-host interactions in light of host population structure and ecology. The methods developed in this thesis include streamlining of standard bacteriophage protocols, such as the agar overlay, and facilitate higher throughput in the isolation and characterization of novel environmental virus-host systems. Here, 〉1300 newly isolated Vibrio are assayed for infection by viral predators and susceptibility is found to be common, though total concentrations of predators are highly skewed, with most present at low abundance. The largest phylogenetically-resolved host range cross test available to date is conducted, using 260 viruses and 277 bacterial strains, and highly-specific viruses are found to be prevalent, with nearly half infecting only a single host in the panel. Observations of blocks of multiple viruses with nearly identical infection profiles infecting sets of highly-similar hosts suggest that increases in abundance of particular lineages of bacteria may be important in supporting the replication of highly specific viruses. The identification of highly similar virus genomes deriving from different sampling time points also suggests that interactions for some groups of viruses and hosts may be stable and persisting. Genome sequencing reveals that members of the largest broad host-range viral group recovered in the collection have sequence homology to non-tailed viruses, which have been shown to be dominant in the surface oceans but are underrepresented in culture collections. By integrating host population structure with sequencing of over 250 viral genomes it is found that viral groups are genomically cohesive and that closely-related and co-occurring populations of bacteria are subject to distinct regimes of viral predation.

Description: I also gratefully acknowledge the WHOI Ocean Ventures Fund, which provided funding for the sequencing of over 250 viral genomes of the Nahant Collection and thereby contributed immensely to the impact of the thesis work presented here. Work presented in this thesis was also made possible by support from National Science Foundation grant DEB 0821391, National Institute of Environmental Health Sciences grant P30- ES002109, the Moore Foundation and the Broad Institute’s SPARC program.

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 December 1996

Description: Two-thirds of the surface of the Earth is created at mid-ocean ridges where magmas rise from the mantle and cool to form the oceanic crust. The objective of this Thesis is to examine the influence of magma supply and eruptive processes on axial morphology, crustal construction, and the properties of crustal magma chambers at intermediate and fast spreading ridges. Variations in magma supply on time scales of ~100 Kyr generate along-axis changes in crustal thickness and temperature. Magma sill properties and hydrothermal activity are closely linked to spreading events which occur on much shorter time scales (ca. 10-100 yr) than the longer-term variations in magma supply reflected in along-axis changes in ridge morphology. The seismically constrained depths of ridge crest magma sills (〉1-2 km) are considerably deeper than the level of neutral buoyancy (100-400 m). The apparent inverse relationship between magma sill depth and spreading rate suggests that a thermally controlled permeability boundary, such as the solidus horizon, controls the depth at which magma ponds beneath mid-ocean ridges. Recent thermo-mechanical models predict that, at intermediate spreading rates, rift valley and magma sill formation are sensitive to small changes in crustal thickness and mantle temperature. Analysis of gravity at an intermediate spreading ridge shows that small differences in crustal thickness (300-700 m) and mantle temperature (10-15°C) are indeed sufficient to produce major changes in lithospheric strength and axial morphology. A stochastic model for the emplacement of dikes and lava flows with a bimodal distribution of lava flows is required to satisfy geological and geophysical constraints on the construction of the extrusive section. Most dikes are intruded within a narrow zone at the ridge axis. Short flows build up approximately half the extrusive volume. Occasional flows that pond at a considerable distance off-axis build up the remainder of the extrusive section. This Thesis underlines the importance of eruption dynamics in the emplacement of the uppermost volcanic layer of the crust and of the crustal thermal structure in controlling local variations in magma sill depth and ridge morphology.

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 1996

Description: We analyze bathymetric and gravity anomalies at five plume-ridge systems to constrain crustal and mantle density structure at these prominent oceanic features. Numerical models are then used to explore the physical mechanisms controlling plume-ridge interaction and to place theoretical constraints on the temperature anomalies, dimensions, and fluxes of the Icelandic and Galapagos plumes. In Chapter 1 we analyze bathymetric and gravity anomalies along the hotspot-influenced Galapagos Spreading Center. We find that the Galapagos plume generates along-axis bathymetric and mantle-Bouguer gravity anomalies (MBA) that extend 〉500 km east and west of the Galapagos Islands. The along-axis MBA becomes increasingly negative towards the plume center, reaching a minimum of ~-90 mGal near 91°W, and axial topography shallows by ~1.1 km toward the plume. These variations in MBA and bathymetry are attributed to the combined effects of crustal thickening and anomalously low mantle densities, both of which are due to a mantle temperature anomaly imposed beneath the ridge by the Galapagos plume. Passive mantle flow models predict a temperature anomaly of 50±25°C is sufficient to produce the 2-4 km excess crust required to explain the along-axis anomalies. 70-75% of the along-axis bathymetric and MBA variations are estimated to arise from the crust with the remaining 25-30% generated by the anomalously hot, thus low-density mantle. Along Cocos-plate isochrons, bathymetric and MBA variations increase with increasing isochron age, suggesting the subaxial mantle temperature anomaly was greater in the past when the plume was closer, to the ridge axis. In addition to the Galapagos plume-ridge system, in Chapter 2 we examine alongisochron bathymetric and MBA variations at four other plume-ridge systems associated with the Iceland, Azores, Easter and Tristan hotspots. We show that residual bathymetry (up to 4.7 km) and mantle-Bouguer gravity anomalies (up to -340 mGal) are greatest at on-axis plumes and decreases with increasing ridge-hotspot separation distance, until becoming insignificant at a plume-ridge separation of ~500 km. Along-isochron widths of bathymetric anomalies (up to 2700 km) decrease with increasing paleo-spreading rate, reflecting the extent to which plume material flows along-axis before being swept away by the spreading lithosphere. Scaling arguments suggest an average ridgeward plume flux of -2.2x106 km/my. Assuming that the amplitudes of the MBA and bathymetric anomalies reflect crustal thickness and mantle density variations, passive mantle flow models predict maximum subaxial mantle temperature anomalies to be 150-225°C for ridge-center plumes, which decrease as the ridges migrate away from the plumes. The dynamics of mantle flow and melting at ridge-centered plumes are investigated in Chapters 3 using three-dimensional, variable-viscosity, numerical models. Three buoyancy sources are examined: temperature, melt depletion, and melt retention. The width W to which a plume spreads along a ridge axis depends on plume volume flux Q, full spreading rate U, buoyancy number B = (QΔρg)/(48η0U2), and ambient/plume viscosity contrast ϒ according to W=2.37(Q/U)l/2(Bϒ)0.04. Thermal buoyancy is first order in controlling along-axis plume spreading while latent heat loss due to melting, and depletion and retention buoyancy forces contribute second order effects. Two end-member models of the Iceland-Mid-Atlantic Ridge (MAR) system are examined. The first endmember model has a broad plume source of radius 300 km, temperature anomaly of 75°C, and volume flux of 1.2xl07 km3/my. The second model has a narrower plume source of radius 60 km, temperature anomaly of l70°C, and flux of 2.1 x106 km3/my. The first model predicts successfully the observed crustal thickness, topographic, and MBA variations along the MAR, but the second model requires substantial along-axis melt transport in order to explain the observed along-axis variations in crustal thickness, bathymetry, and gravity. We favor this second model because it predicts a mantle P-wave velocity reduction in the plume of ~2% as consistent with recent seismic observations beneath Iceland. Finally in Chapter 4 we use three-dimensional numerical models to investigate the interaction of plumes and migrating midocean ridges. Scaling laws of axial plume spreading width Ware derived first for stationary ridges and off-axis plumes, which yield results consistent with those obtained from independent studies of Ribe [1996]. Wand the maximum plume-ridge interaction distance Xmax again scale with (Q/U)l/2 as in the case of ridge-centered plumes and increase with ϒ and buoyancy number. In the case of a migrating ridge, Xmax is reduced when a ridge migrates toward the plume due to excess drag of the faster-moving leading plate, and enhanced when a ridge migrates away from the plume due to reduced drag of the slower-moving trailing plate. Thermal erosion of the lithospheric boundary layer by the previously ridge-centered plume further enhances Wand Xmax but to a degree that is secondary to the differential migration rates of the two plates. Model predictions are compared with observed along-isochron bathymetric and MBA variations at the Galapagos plume-ridge system. The anomaly amplitudes and widths, as well as the increase in anomaly amplitude with age are predicted with a plume source temperature anomaly of 80-120°C, radius of 80-100 km, and volume flux of 4.5x106 km3/m.y. Our numerical models also predict crustal production rates of the Galapagos Islands consistent with those estimated independently using the observed island topography. Predictions of the geochemical signature of the plume along the present-day ridge suggest that mixing between the plume and ambient mantle sources is unlikely to occur in the asthenosphere or shallow crust, but most likely deeper in the mantle possibly by entrainment of ambient mantle as the plume ascends through the depleted portion of the mantle from its deep source reservoir.

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 March 1997

Description: The formation of new oceanic crust is the result of a complex geodynamic system in which mantle rises beneath spreading centers and undergoes decompression melting. The melt segregates from the matrix and is focused to the rise axis, where it is eventually intruded and/or erupted to form the oceanic crust. This thesis combines surface observations with laboratory studies and geodynamic modeling to study this crustal-production system. Quantitative modeling of the crustal and mantle contributions to the axial gravity and topography observed at the East Pacific Rise shows that the retained melt fraction in the mantle is small (〈3%) and is focused into a narrow column extending up to 70 km beneath the ridge axis. Consistent with geochemical constraints, the extraction of melt from the mantle therefore appears to be efficiently focus melt toward the ridge axis. A combination of laboratory and numerical studies are used to constrain the pattern of mantle flow beneath highly-segmented ridges. Even when the buoyant component of mantle flow is constrained to be two-dimensional, laboratory studies show that a segmented ridge will drive three-dimensional mantle upwelling. However, using reasonable mantle parameters in numerical models, it is difficult to induce large-amplitude three-dimensional mantle upwelling at the relatively short wavelengths of individual segments (~50 km). Instead, a simple model of three-dimensional melt migration shows that the observed segment-scale variations in crustal thickness can be explained by focusing of melt as it upwells through a more two-dimensional mantle flow field. At the Reykjanes Ridge, the melt appears to accumulate in small crustal magma chambers, before erupting in small batches to form numerous overlapping hummocky lava flows and small volcanoes. This suggests that crustal accretion, particularly at slow-spreading centers, may be a highly discontinuous process. Long-wavelength variations in crustal accretion may be dominated by variations in mantle upwelling while short-wavelength, segment-scale variations are more likely controlled by a complex three-dimensional processes of melt extraction and magma eruption.

Description: During my first three years in the Joint Program, I was supported by an National Science Foundation Graduate Student Fellowship. Other support has been derived from National Science Foundation grants OCE-9296017, OCE-9224738, OCE-9215544, and EAR grant 93-07400.

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 December 1996

Description: The objective of this Thesis was to interpret the structural development of slowspreading ridge segments by: 1) delineating the nature, magnitude, and relative importance of primary tectonic and volcanic processes that control crustal morphology, 2) investigating the spatial and temporal variability of these processes, and 3) examining how rheological variations in the lithosphere control its structural configuration. To that end, this Thesis provides detailed documentation of faults and volcanoes (seamounts) at the Mid-Atlantic Ridge from 25°25'N to 27°10'N and extending from zero-age crust at the ridge axis to -29 Ma crust on the ridge flank. This information was used to analyze the evolution of ocean crust from initial formation in the rift valley to degradation by aging processes on the ridge flank. Accumulation of sediments affects the seafloor morphological expression of ocean crustal structure, and sediment thicknesses were also mapped to facilitate study of the morphological record of crustal accretion and tectonism. In addition, deformation conditions in the lithosphere were analyzed by study of microstructure and geothermometry of abyssal peridotite mylonites recovered from fault zones at slow-spreading ridges.

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 1995

Description: The global mid-ocean ridge system is one of the most striking geological features on the surface of the Earth. In this system, the East Pacific Rise (EPR) is the fastest spreading ridge and is thus considered as the most active magmatically among the plate boundaries. In January and February of 1988, an extensive survey by the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution was conducted along the EPR between 9°05' and 9°55'N to study the crustal structure of the axial region. This thesis, the result of that cruise, comprises four main topics: (1) characterization of normal faulting from Sea Beam bathymetric data, (2) application of mechanical models to explore the hypothesis that buoyancy arising from crustal magma chambers and gravitational spreading of the upper crust are the principal processes leading to the initiation and development of normal faults, (3) investigation of seafloor magnetization anomalies to constrain upper crustal structure, and (4) analysis of gravity anomalies to examine possible correlations between observed variations in seafloor manifestations of volcanism and deformation and underlying structure. Thus, each topic focuses on different levels of the mid-ocean ridge. Together with the results of seismic and other observations, the findings are woven into a better understanding of the tectonic processes and structure of fastspreading mid-ocean ridges. First, to understand the characteristics of normal faults at fast-spreading ridges, we utilized swaths of Sea Beam bathymetry and estimated the distribution and geometry of normal fault zones using the slope of the seafloor as the criterion for a faulted surface. In our survey area, nonnal fault activity begins 2-8 km off-axis and continues at least to 30-40 km from the axis, as indicated by an increase in the total and average throws of normal fault zones versus distance from the axis. There appears to be no significant difference in the plan-view area of inward- and outward-facing nonnal fault zones. The distance from the rise axis to the nearest large-offset fault zone (throw 〉 20 m) on either side of the axis is approximately symmetric to the north of 9°23'N, but the midpoint between nearest largeoffset fault zones is offset 2-3 km to the west of the bathymetric axis to the south of 9°23'N. The continued growth of nonnal fault zones suggests that significant extensional stress persists to greater distances from the axis than previously thought and that the rise axis possesses a finite strength. The argument that the rise axis has finite strength is consistent with recent evidence for solidified axial dikes along magmatically active portions of the EPR from near-bottom seismic refraction experiments, which suggests that, while eruption of magma at the rise axis weakens the axis, the persistence of such weak zones is short-lived and the emplacement zones at any given time are localized along the axis. We examined how the presence of a low-density, low-strength magma chamber within the crust and gravitational spreading of a mechanically strong upper crust over an underlying substrate contribute to the fonnation of faults at a fast spreading mid-ocean ridge by comparing the predicted stress field with the observed pattern of normal fault zones. We employed boundary element methods to incorporate buoyancy and gravitational spreading as body forces in an elastic medium, and we detennined stress and strain fields for a variety of rise axis conditions and a range of possible sets of material properties for different parts of the mid-ocean ridge. Our results show that the strength of the rise axis is one of the most crucial factors governing the near-axis stress field. If the rise axis is mechanically weak, the maximum extensional stress from buoyancy occurs at shallow depth off the rise axis. A weak rise axis may result from recent magmatism such as the intrusion of dikes into the upper crust. On the other hand, if the rise axis is mechanically strong, which may result after solidification and cooling of the dike zone, the maximum surface extensional stress occurs on the rise axis. However, the reduction in size of a magma chamber that would accompany cessation of dike injection would lead to less buoyancy and thus a lower likelihood of stress levels sufficient for faulting. For a given set of material strengths and a given magnitude of buoyancy force, the flexural rigidity of the upper crust plays an important role in detennining if a zone of extension will develop off axis and, if so, the position and horizontal extent of that zone. A thin or mechanically weak upper crust is more likely to develop a zone of extension than one that is thick or mechanically strong. The stress field resulting from gravitational spreading is similarly affected by the strength of the rise axis. While buoyancy can explain a consistent distance at which normal faults initiate off-axis, gravitational spreading can account for continued activity on normal faults to a greater distance from the axis than can buoyancy. The existence of a magma lens can play an important role in reducing the magnitude of the stress field for a weak rise axis, as the crust above the magma lens can slide and thus relieve the thickness-averaged extensional stress. Next, we inverted surface ship measurements of the scalar magnetic field along the EPR between 9°10' and 9°50'N. We examined whether the axial magnetization high, which increases in amplitude to the south in our area, can best be explained by variations in the thickness or in the magnetization intensity of the source layer. The variation in axial magnetization is too large to be explained solely by the variations in the depth to the top of the axial magma chamber indicated by reflection seismology. For a magnetic source layer that is 500 or 750 m thick, the observed along-axis variations in FeO and Ti02 explain only 36 and 60%, respectively, of the total variance of axial magnetization anomalies. Therefore, a combination of variations in magnetic layer thickness and in intensity of magnetization (by variations in the FeO and Ti02 contents of the source rock or by other mechanisms) is needed to explain the along-axis variation of axial magnetization. In addition to the increase in amplitude to the south, the axial magnetization high exhibits at least three marked changes in magnitude and offsets in its along-axis linearity ('magnetic devals') (at 9°25', 9°37', and 9°45'N) which appear to be related to boundaries or offsets between the segments of the axial summit caldera (ASC). Because the amplitudes of the axial magnetization anomalies are highest at the midpoints of the ASC segments, we speculate that midpoints of the ASC segments are the loci of more frequent lava eruptions, and the seafloor basalts at the midpoints are thus younger and more magnetic, than at the segment ends. The magnetization shows distinct short-wavelength (~ 5 km) banding to the north of 9°25'N over a region that does not appear to have been affected by an overlapping spreading center. Among the possible explanations for these off-axis magnetization anomalies are short geomagnetic reversal events within the Brunhes epoch, variations in the paleointensity of the Earth's field, variations in the magnetization intensity of the source rock due to variability in the magmatic supply, and variations in the degree of hydrothermal alteration at the rise axis. On the basis of comparisons of forward models and observations, short geomagnetic reversal events appear to be the most likely explanation of these anomalies. The analysis of sea-surface gravity field measurements shows an axial residual mantle Bouguer gravity anomaly too large to be explained by the anomalous temperature of the mantle or by changes in the thickness of the crust. The broad axial residual gravity low is interpreted as a signal arising largely from the upper mantle, presumably by presence of partial melt along the rise axis. A northward increase in the width of the low implies a greater melt fraction in the region to the north than to the south, especially on the Pacific plate side. The residual gravity anomaly also shows several short-wavelength local lows along the axis (e.g., 9°21', 9°32', and 9°42'N) which correlate with along-axis variations in axial magnetization and tomographic images of mid-crustal seismic velocities. Along axis the local lows have an amplitude of 1.5-3 mGal and appear at a nearly regular spacing (10- 15 km). Across the axis, however, the local lows show a greater variation (3-5 mGal), suggesting that there is an additional gravity anomaly signal arising from a low-density structure that is approximately continuous along the axis. The anomalous masses producing the local lows are interpreted as zones of relatively high melt concentration, formed within the crust by recent replenishment of magma from the upper mantle, that are surrounded by a region of lesser melt concentration corresponding to the low-velocity volume imaged by seismic tomography. If the zone of high melt concentration are modeled as circular rods of radius 1 km, along-axis length 10 km, and center of mass 2.25 km below the seafloor, density contrasts of 200-350 kg/m3 are needed to match the observed anomalies. For larger anomalous mass volumes, the density contrasts would be lower. The findings of this study support the hypothesis that the axis of the EPR can be divided into segments 10-15 km in length, with each segment defmed by the locus and timing of most recent emplacement of magma in the axial crust. The segments in the study area appear to be in different phases of a magmatic cycle, but the period of such a magmatic cycle is not known. By this view, the discrete emplacement of magma bodies gives rise to along-axis variations in crustal structure manifested as short-wavelength residual gravity anomalies and magnetic devals. Another consequence of a rise axis at which magma is emplaced at discrete locations is that the mechanical strength of the axial upper crust varies with position along the axis and over time. During active magmatism, the rise axis acts as a weak zone and the buoyancy of the axial magma chamber and surrounding low-velocity volume can lead to initiation of off-axis normal faulting. However, for a long segment of the rise bounded by transform faults, the axis will have sections with a solidified rucial injection zone as well as sections undergoing active magmatism, and thus the rise overall may appear to have finite strength. If such a finitestrength ridge axis is subject to significant extensional stress as a result of gravitational spreading, mantle convective tractions, or differential cooling, then continued normal fault activity would extend over a broad region to distances of at least several tens of kilometers from the spreading axis.

Description: The work in this thesis was supported by the National Science Foundation under grants OCE-8615797, OCE-8615892, and OCE-9000177.

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 1999

Description: Changes in morphology of the Marquesas Fracture Zone are correlated with small changes in Pacific-Farallon relative motion. The simple flexural signal of a locked fracture zone may be obscured by tectonic effects, and there is no evidence for the release of shear stress on the fracture zone by vertical slip after leaving the active transform. One such small change in plate motion is documented in the South ern Austral Island region of the South Pacific. A twelve degree clock wise change in Pacific-Farallon relative motion occurred around fifty million years ago. This Eocene change in spreading direction and rate is locally constrained with observations of magnetic anomalies and spreading fabric orientation. At the southeastern end of the Cook-Austral Island chain, multiple episodes of volcanism have left a diverse population of seamounts. Volume estimates from geophysical data and modeling show that one-half to two-thirds of the volcanic material is over thirty million years old, while the remainder is less than five million years old. Seismic and bathymetric data imply the presence of abyssal basalt flows in the flexural moat of the Austral Islands, probably associated with Austral Islands volcanism, which may contribute a significant amount of material to the archipelagic apron.

Description: The research presented in Chapter 2 was supported by National Science Founda tion grants OCE-9012949 and OCE-9012529. Chapters 3, 4 and 5 were supported by National Science Foundation grant OCE-9415930. A National Science Foundation graduate fellowship supported my first three years of graduate study.

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 May 1988

Description: Heat flow, tectonic subsidence and crustal thickness distributions in the Ligurian Basin are best explained by asymmetric lithospheric thinning mechanisms. Over 150 heat flow measurements are made on several transects between Nice, France and Calvi, Corsica on continental slope and rise settings. Thermal gradient determinations are improved using an optimization technique. Piston core data and surface sediment 3.5 kHz reflectivity patterns help constrain thermal conductivity obtained from over 100 in situ stations. Plio-Quaternary stratigraphy is revised using new seismic reflection profiles: a boundary fault system associated with postrift margin uplift, a Pleistocene-age Var Fan construction, and recent diapirism of Messinian salt are indicated. After assessing local thermal disturbances (mass-wasting, microtopography, and salt refraction), positive heat flow corrections are made for multi-lithologic sedimentation histories and glacial paleotemperatures. Using boundary-layer cooling models, equilibrium heat flow estimates support geologic evidence for Oligocene and early Miocene rifting. Heat flow maxima correlate well with two "oceanic" sub-basins, suggesting that the southeastern trough near Corsica is ~5 Myr younger, consistent with the southeastern progression of volcanism and back arc rifting in the Western Mediterranean. Tectonic subsidence-crustal thickness trends indicate lithospheric stretching, with heat flow supporting asymmetric sub-crustal lithospheric thinning during the conjugate margin formation.

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 May 1988

Description: The age of the Nazca plate where it enters the Peru and northern Chile trenches varies from 30 Ma in the north to 45 Ma in the south as its dip beneath the South American continent steepens from 13° to 30°. If the elastic thickness Te of oceanic lithosphere depends only on its age, and therefore thermal state, we would expect that Te determined from fitting the flexure of the lithosphere over the outer rise as revealed in the depth and geoid anomalies would increase from the Peru Trench in the north to the northern Chile Trench further south. We find that just the opposite is true: the lithosphere appears stiffer outboard of the Peru Trench than it does further south, and the isotherm controlling the elastic/ductile transition must be 600°C or greater if the thermal structure of the plate is that predicted by the standard thermal plate model. Because the decrease in plate stiffness to the south is correlated with a decrease in the minimum radius of curvature of the flexed plate over the outer rise and outer trench wall, we interpret our result in terms of inelastic yielding of the oceanic lithosphere when bent to high strains. The fact that the more highly bent segment of subducting lithosphere also dips at a steeper angle at greater depth beneath the continent might suggest that the amount of inelastic weakening of the lithosphere could be predicted from seismic images of the down going slab, but we find little support for this correlation worldwide. The forces and moments controlling the shallow deformation of the plate seaward of the trench do not appear to be linked to upper mantle processes which impose the dip at greater depth. Finally, we consider the possibility that the elastic thickness of the lithosphere would be reduced for trenches that are highly arcuate in map view, again due to inelastic yielding. If such a relationship exists, the effect for oceanic lithosphere is much smaller than what is documented for continental plates where they underthrust highly arcuate fold belts.

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 1984

Description: The distribution and feeding behavior of bacterivorous micro flagellates (2-20 μm protozoa) and their ingestion by copepods were examined in an attempt to assess the importance of these protozoa as a trophic link between planktonic bacteria and zooplankton. The abundance of microflagellates relative to other picoplankton (0.2-2.0 μm) and nanoplankton (2-20 μm) populations in water samples in the North Atlantic and in Lake Ontario and on macroaggregates in the North Atlantic was determined using direct microscopical and culture estimation techniques. Seasonal, vertical and geographical changes in the density of microflagellates were generally not greater than one order of magnitude. Microscopical counts of heterotrophic nanoplankton (presumably microflagellates) typically ranged from a few hundred to a few thousand m1-1 for a variety of planktonic environments. They constituted approximately 1/3 to 1/2 of the nanoplankton in the euphotic zone and dominated the nanoplankton in the aphotic zone. Most Probable Number (MPN) estimation of the density of bacterivorous protozoa indicated that microflagellates were, on average, an order of magnitude more abundant than bacterivorous ciliates and amoebae. MPN and direct microscopical counts of microflagellates differed by as much as 104. This discrepancy was smaller in eutrophic environments (e.g. Continental Shelf and Lake Ontario) and on macroscopic detrital aggregates. All microbial populations enumerated were highly concentrated on macroscopic detrital aggregates relative to their abundance in the water surrounding the aggregates. Enrichment factors (the ratio of abundance of a population on a macroaggregate to its abundance in the surrounding water) increased along a eutrophic-to-oligotrophic gradient because of the combined effects of an increased abundance of microorganisms on macroaggregates in oligotrophic environments and a decreased abundance in the surrounding water in these same environments. Average enrichment factors for direct microscopical counts of heterotrophic nanoplankton (range = 17-114) were not as large as enrichment factors observed for MPN estimates of the number of bacterivorous microflagellates (range = 273-18400). Microflagellates numerically dominated the bacterivorous protozoa cultured from macroaggregates by one to two orders of magnitude, but ciliates and amoebae were also highly enriched on macroaggregates. Microenvironments are therefore a potentially important aspect for the ecology of planktonic microorganisms. Observations on the microbial colonization of mucus sloughed by ctenophores and discarded appendicularian houses suggest that these materials may be important sources of macroaggregates. Batch and continuous culture experiments were conducted with clonal cultures of microflagellates to test their ability to grow on various types and densities of bacteria. The doubling time of Monas sp. 1 ranged from 43 hr (when fed the cyanobacterium Synechococcus Strain WH 8101) to 6.9 hr (when fed the heterotrophic bacterium Serratia marinorubra). Cell yields (i.e. the conversion of bacterial biomass into protozoan biomass) of Monas sp. 1 fed two species of heterotrophic bacteria were greater than yields for the microflagellate fed two species chroococcoid cyanobacteria (range = 7-68%). Cell yields of two other species of microflagellates (Monas sp. 2 and Cryptobia maris) were 48% and 61%, respectively, on the bacterium Pseudomonas halodurans. Microflagellates grew in continuous culture at concentrations of bacteria which were lower than bacterial densities required for the growth of ciliates as shown by other investigations. Therefore, microflagellates appear to be well-adapted for grazing bacterioplankton. Microflagellates were also investigated for their ability to graze bacteria attached to particles. Bodo nanorensis and Rhynchomonas nasuta both showed a marked ability to graze attached bacteria and a limited ability to graze unattached cells. These results suggest that microflagellates may also be important consumers of bacteria attached to particles in the plankton and may explain the highly elevated densities of microflagellates on macroaggregates. Grazing experiments performed with the copepod Acartia tonsa indicated that heterotrophic microflagellates were ingested by the copepods at rates comparable to the ingestion of phytoplankton of similar size. The presence of heterotrophic microflagellates did not depress filtration rates of the copepods, and one species (Cryptobia maris) appeared to be selectively grazed. Survival of A. tonsa on a diet of heterotrophic microflagellates was similar to survival on a diet of phytoplankton and was significantly longer than survival of starved Controls or copepods fed only bacteria. Due to their ability to grow at in-situ densities of planktonic bacteria, their relatively high cell yields, and their acceptability as food for zooplankton, it is concluded that bacterivorous microflagellates may constitute an important trophic link between bacteria and zooplankton. This link may provide a mechanism whereby organic material and energy from the detrital food chain can be returned to the classical phytoplankton-copepod-fish food chain.

Description: This research was supported by National Science Foundation grants OCE80-2444l and OCE82-l4928 and Ocean Industry Program grant 4473 awarded to Dr. Laurence P. Madin, NSF Doctoral Dissertation grant OCE8l-l299l, the Woods Hole Oceanographic Institution Education Program and the Wood Hole Oceanographic Institution Biology Department.

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 2009

Description: Trichodesmium is a colonial, N2-fixing cyanobacterium found in tropical oceans. Species of Trichodesmium are genetically similar but several species exist together in the same waters. In order to coexist, Trichodesmium spp. may occupy different niche spaces through differential utilization of resources such as nutrients and light, and through responses to physical characteristics such as temperature and turbulence. To investigate niche differentiation in Trichodesmium, I characterized cultured strains of Trichodesmium, identified and enumerated Trichodesmium clades in the field, and investigated P stress and N2 fixation in field populations. Species of Trichodesmium grouped into two clades based on sequences from 16S rDNA, the internal transcribed spacer (ITS), and the heterocyst differentiation gene hetR. Clade I contained Trichodesmium erythraeum and Trichodesmium contortum, and clade II contained Trichodesmium thiebautii, Trichodesmium tenue, Trichodesmium hildebrandtii, and Trichodesmium pelagicum. Each clade was morphologically diverse, but species within each clade had similar pigmentation. I developed a quantitative polymerase chain reaction (qPCR) method to distinguish between these two clades. In field populations of the Atlantic and Pacific Oceans, the qPCR method revealed that clade II Trichodesmium spp. were more prominent than clade I in the open ocean. Concentrations of Trichodesmium did not correlate with nutrient concentrations, but clade I had wider temperature and depth distributions than clade II. Temperature and light are physical characteristics that may define niche spaces for species of Trichodesmium. Clade I and II concentrations correlated with each other in the Pacific but not in the Atlantic, indicating that the two clades were limited by the same factors in the Pacific while different factors were limiting the abundance of the two clades in the Atlantic. Trichodesmium populations in the North Atlantic were more P stressed and had higher N2 fixation rates than populations in the western Pacific. While nutrient concentrations didn’t directly correlate with Trichodesmium concentrations, the contrasting nutrient regimes found in the Atlantic and Pacific Oceans might influence distributions of the two clades differently. Unraveling the differences among species of Trichodesmium begins to explain their coexistence and enables us to understand factors controlling global N2 fixation.

Description: National Science Foundation (NSF) Biocomplexity Program Grant (OCE-0323332); the Center for Microbial Oceanography Research and Education (C-MORE), an NSF Science and Technology Center (EF-0424599); the Woods Hole Oceanographic Institution (WHOI) Ocean Life Institute (OLI) grant to J. Waterbury, 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 September 1978

Description: This thesis consists of three papers examining problems related to the crustal structure, isostasy and subsidence history of aseismic ridges and mid-plate island chains. Analysis of gravity and bathymetry data across the Ninetyeast and eastern Walvis Ridges indicates these features are locally compensated by an over thickening of the oceanic crust. Maximum crustal thicknesses are 15-30 km. The western Walvis Ridge is also compensated by crustal thickening; however, the isostasy of this part of the ridge is best explained by a plate model of compensation with elastic plate thicknesses of 5-8 km. These results are consistent with the formation of the Ninetyeast and Walvis Ridges near spreading centers on young lithosphere with flexural rigidities at least an order of magnitude less than those typically determined from flexural studies in older parts of the ocean basins. As the lithosphere cools and thickens, its rigidity increases, explaining the differences in isostasy between aseismic ridges and mid-plate island chains. The long-term subsidence of aseismic ridges and island/ seamount chains can also be explained entirely by lithospheric cooling. Aseismic ridges form near ridge crests and subside at nearly the same rate as normal oceanic crust Mid-plate island chains subside at slower rates because they are built on older crust. However, some island chains have subsided faster than expected based on the age of the surrounding sea floor, probably because of lithospheric thinning over midplate hot spots, like Hawaii. This lithospheric thinning model has major implications both for lithospheric and mantle convection studies as well as the origin of continental rift systems.

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 January 1972

Description: Free-air and simple Bouguer anomaly maps of the Venezuelan continental margin (from 60°W to 72°W and from 7°N to 13°N) are presented. The major features of the free-air map are: the large lows associated with the deep sedimentary basins, -200 mgal in the Eastern Venezuela basin and -164 mgal in the Maracaibo basin; the high of greater than 300 mgal over the Venezuelan Andes; and a belt of highs associated with the offshore islands extending from Blanquilla to Curacao and then over the Guajira peninsula, where they terminate. The Bouguer anomaly map shows a large low (-196 mgal) over the Eastern Venezuela basin and relative minimums over the coastal mountains. A minimum associated with the Venezuelan Andes is shifted to the northwest of the topographic axis and lies over the flank of the Andes and part of the Maracaibo basin. Using the gravity data, structural sections were constructed for a series of profiles across the Venezuelan Andes and Caribbean mountains. They show that there is no light crustal root under the Andes, the relative mass excess is as much as 600 kg/cm2, and that there is an excess of low density material under the Maracaibo basin. This appears to be caused by a combination of a southeastward dipping shear zone in the lithosphere under the basin-mountain boundary and a component of compressive stress perpendicular to this zone, both of which have resulted in the uplift of the crust under the Andes, and downwarp under the basin. The apparent flexural rigidity of the lithosphere under the Maracaibo basin is 0.6 x 1023 newton-m, a normal value for lithosphere deformations of Miocene age. The Caribbean mountains have a light crustal root which has been formed by the sliding of blocks of crustal material from the north over the rocks to the south, and perhaps by the underthrusting of oceanic crust under the continental crust. This underthrusting may have been a result of the formation of a downgoing slab of lithosphere along the Venezuelan continental margin during the late Cretaceous. The downgoing slab may have existed until mid-Eocene time. The gravity minimum over the Eastern Venezuela basin is due to the downwarping of lighter crustal material into the higher density mantle. This may be a result of compression from the north along a north-south direction causing plastic downbuckling of the lithosphere. The present deformation along the northern boundary appears to be due to differences in relative motion between the North and South American plates. Because the Caribbean mountains are partially isostatically compensated, while the Venezuelan Andes are above isostatic equilibrium, this suggests that the relative motion of the Caribbean plate with respect to the South American plate is eastward. The compressive stress across the boundary in the region of the Venezuelan Andes is probably greater than the compressive stress across the Caribbean mountains.

Description: This investigation was supported by a grant from the National Science Foundation GA-12204, and by contract N00014-66-C-024l with the Office of Naval Research.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Science at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution January 1981

Description: The results of a detailed geophysical survey are used in conjunction with all available information in a study of the tectonic development of the Cayman Trough and the Greater Antilles Ridge. This development is connected with the relative motions of the North and South Americas and the eastern Pacific plates. Thus, the pre-Tertiary history of the region is one of simple convergence. This contrasts with the complex tectonism of primary translation, with secondary convergence and divergence during the Tertiary. The ancestral Greater Antillean Arc suffered fracturing during collision with the Bahamas stable platform in the Late Cretaceous. Oblique convergence re-established itself across the remnant fragments of the ancestral arc in the Tertiary, producing a sheared welt partially decoupled from both the North American and Caribbean plates. Pronounced temporal and structural heterogeneity occurs within this Plate Boundary Zone. Along its northern margin secondary convergence with the North American plate formed the massive subduction complex of the Cuchillas Uplift and the Sierra Septentrional. Convergence between the Plate Boundary Zone and the Caribbean plate resulted in the triple virgation of the fold belts extending westward from the Los Muertos Trough to Oriente Province (Cuba), the Cayman Trough and the Nicaraguan Rise. Tectonism along these fold belts youngs southwestward preserving the stratigraphy of the Caribbean Basin at the time of their formation during the early, middle, and late Tertiary. The Caribbean/North American Plate boundary occurred along the zones of major strain accomodation within the Plate Boundary Zone. The Cayman Trough was produced during a period of divergence between the Nicaraguan Rise and the North American plates during the Miocene. Since the Pliocene, the shear boundary within the Cayman Trough occurs along the Oriente Deep proceeding via the Windward Passage Deep and the Valle del Cibao to the Puerto Rico Trench. Convergence and shear predominate the present tectonic framework of the Plate Boundary Zone.

Description: Cruise #97 of the R. V. ATLANTIS II was sponsored by the National Science Foundation (OCE78-20/ 11336.00). Further support was received from the Ocean Industry Program and the Educational Program of 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 June, 1979

Description: A major goal in the study of plate tectonics is the acquisition of a knowledge of the history of relative motion among the rigid plates of the earth's lithosphere. The three papers of this thesis contribute to this effort and demonstrate that studies of the stability and evolution of triple junctions and of the finite rotations of systems of three plates can yield significantly more accurate tectonic histories than can studies of the relative motions between two plates alone. Topographic and magnetic investigation of the Southwest Indian Ridge and reconstruction of the plate system of the Indian Ocean shows that both Africa and Antarctica are rigid plates and their pole of relative rotation has remained fixed near 8°N, 42°W since the Eocene. A detailed survey of the Indian Ocean triple junction reveals that the Indian Ocean plate motions have remained constant since 10 Ma. The stability conditions of the junction show that the general morphology of the Southwest Indian Ridge results from the evolution of the Indian Ocean triple junction. A method is presented for determining the finite rotations best reconstructing the past relative positions of three plates around a triple junction. The method is illustrated by reconstructions of the plates around the Labrador Sea triple junction at the times of anomalies 24 (56 Ma) and 21 (50 Ma). The region of uncertainty of the Greenland-North America finite pole is mapped for each reconstruction, and it demonstrates that consideration of the three plate system yields more well-constrained results than does a treatment of the two plates alone.

Description: This work was supported by the Office of Naval Research contract N00014-75-C-0291 with the Massachusetts Institute of Technology.

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 April 1978

Description: The structure of the membrane-free nucleoid of Escherichia coli and of unfolded chromosomal DNA was investigated by sedimentation on neutral sucrose gradients after irradiation with 60Co gamma-rays and ultraviolet light (2S4nm). Irradiation both in vivo and in vitro was used as a molecular probe of the constraints on DNA~packaging in the bacterial chromosome. The extremely gentle lysis and unfolding procedures which were developed yielded undamaged, replicating genomes, thus permitting direct measurement of the formation and repair of DNA double-strand breaks at biologically-significant doses of ionizing radiation. In vitro UV-irradiation of nucleoids resulted in an increase in the observed rate of sedimentation due to the formation of an unknown photo-product. In contrast, UV-irradiation of wild-type cells in vivo showed evidence of the formation of incision breaks which resulted in the relaxation of supercoiling in the nucleoid. Strand breakage was also observed following in vivo UV-irradiation of a uvrB-5 strain, but at a lower rate and also accompanied by considerable unfolding of the chromosome. Such lesions may have been the result of direct photochemical reactions in the nucleoid, or enzyme activity associated with a uvr-independent mode of repair. The number of domains of supercoiling was estimated at 170 per genome equivalent of DNA based on measurements of relaxation caused by single-strand break formation in in vivo- and in vitro-gamma-irradiated folded chromosomes. Similar estimates based on the target size of RNA molecules responsible for maintaining the compact packaging of the nucleoid predicted negligible unfolding due to the formation of RNA single-strand breaks at doses up-to 10 Krad, and were born out by experimental measurements. Unfolding of the nucleoid in vitro by limit-digestion with RNase or by heating at 70° resulted in DNA complexes with sedimentation coefficients of 1030±59S and 625±15S respectively. The difference in these rates was apparently due to more complete deproteinization and thus less mass in the heated material. These structures are believed to represent intact, replicating genomes in the form of complex-theta structures containing 2-3 genome equivalents of DNA. The rate of formation of double-strand breaks was determined from molecular weight measurements of thermally unfolded chromosomal DNA gamma-irradiated in vitro. Break formation was linear with dose up to 10 Krad, resulting in 0.27 double-strand breaks per kilorad per genome equivalent of DNA and requiring 1080 eV/double-strand break. The influence of possible non-linear DNA conformations of these calculations is discussed. Repair of ionizing radiation damage to folded chromosomes was observed within 2-3 hours of post-irradiation incubation in growth medium. A model based on recombinational repair is proposed to explain the formation of 2200-2300S material during early stages of incubation and subsequent changes in the gradient profiles. Such behavior is not observed for post-irradiation incubation of wild-type cells in buffer or for a recA-13 strain incubated in growth medium. Association of unrepaired DNA with plasma membrane is proposed to explain the formation of a peak of rapidly sedimenting material (〉〉3100S) during the later stages of repair. Direct evidence of repair of double-strand breaks during post-irradiation incubation in growth medium was obtained from gradient profiles of DNA from RNAse-digested chromosomes. The sedimentation coefficient of broken molecules was restored to the value of unirradiated DNA after 2-3 hours of incubation, and the fraction of the DNA repaired in this fashion was equal to the fraction of cells which survived at the same dose. An average of 2.7 double-strand breaks per genome per lethal event was observed, suggesting that 1-2 double-strand breaks per genome are repairable in this strain of E. coli.

Description: Submitted in the partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution August, 1996

Description: Two-thirds of the Earth's surface is oceanic crust formed by magmatic and tectonic processes along mid-ocean ridges. Slow-spreading ridges, such as the Mid-Atlantic Ridge, are discontinuous and composed of ridge segments. Segments are thus fundamental units of magmatic accretion and tectonic deformation that control the evolution of the crust. The objective of this Thesis is to constrain the tectonic processes that occur at the scale of slowspreading segments, to identify the factors controlling segment propagation, and to provide constraints on lithospheric strength with laboratory deformation experiments. In chapter 2, bathymetry and gravity from various areas along the global mid-ocean ridge system are analyzed to quantify systematic variations at the scale of individual segments. There is a marked asymmetry in bathymetry and gravity in the vicinity of segment offsets. We develop a model of faulting to explain these observations. Low-angle faults appear to accommodate tectonic extension at the inside corners of ridge-offset intersections, and result in substantially uplifted terrain with thin crust with respect to that at the outside corners or centers of segments. Results from Chapter 3 indicate that the crust magmatically emplaced on axis is not maintained off-axis. This transition is revealed by both statistical and spectral analyses of bathymetry and gravity. Tectonic extension varies along the length of a segment, resulting in thinning and uplift of the crust at ridge-offset inside corners, and a decorrelation between bathymetry and gravity patterns. Tectonic deformation substantially reshapes the oceanic crust that is magmatically emplaced on-axis, and strongly controls the crustal structure and seafloor morphology off-axis. Satellite gravity data over the Atlantic shown in Chapter 4 reveal a complex history of ridge segmentation, and provides constraints on the processes driving the propagation of segments. The pattern of segmentation is controlled mainly by the geometry of the ridge axis, and secondarily by hot spots. Segments migrate primarily down regional gradients associated with hot spot swells. However, the lack of correlation between gradients and propagation rate, and the propagation up gradient of some offsets, suggest that additional factors control propagation (e.g., variations in lithospheric strength). Most non-transform offsets are short-lived and migrating, while transform offsets are long-lived and stable. Both the propagation of segments (Chapter 4) tectonism along a segment (Chapters 2 and 3) are controlled by the lithospheric rheology. In Chapter 5 I present results from laboratory deformation experiments on serpentinite. These experiments demonstrate that serpentinites are considerably weaker than peridotites or gabbros, display a non-dilatant style of brittle deformation, and strain is accommodated by shear cracking. Serpentinites may weaken the lithosphere, enhance strain localization along faults, and control the style of faulting.

Description: A fellowship from Caixa de Pensions "La Caixa" in Barcelona provided me with all the required financial support to come to WHOI. The work presented in this thesis was also supported by the National Science Foundation grants OCE-90l2576, OCE-930078, OCE-9313812, and Office of Naval Research grant N00014-9l-J-1433.

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 1976

Description: This thesis is a collection and analysis of seafloor magnetic anomalies, bathymetry, and the paleomagnetism of DSDP sediments and basalt in the West Philippine Basin, in an attempt to resolve questions about its origin as a marginal basin. Our results suggest that this basin was formed in an Eocene pulse of rapid spreading (v1/2 = 41-44 mm/yr) in a direction (N 21°E) significantly different from later pulses which opened the more eastern basins of the Philippine Sea. The Central Basin Fault appears to be intimately associated with this spreading by nature of its structure and trend, and it may be a remanent of a former ridge system. Our preliminary calculation of paleopole positions also suggests that there was a large amount (60°) of clockwise rotation between this basin and the magnetic pole. This is consistent with rotations of the Pacific plate with respect to the magnetic pole and current directions of Philippine- Pacific'relative rotations. Basement depths of 6 km in the West philippine Basin imply that its crustal and/or lithospheric structure is different from Pacific structure of the same age.

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 January, 1977

Description: This thes is is a collection of papers on the paleomagnetics of samples from several Deep Sea Drilling Project (DSDP) sites in the Indian Ocean. These papers present the basic paleomagnetic data, discuss the statistical methods for analyzing such data from DSDP cores, and examine the implications of the paleolatitudes for the origin of the Ninetyeast Ridge and the northward motion of India. Rarely do DSDP paleolatitudes approach the reliability of good continental pole positions. However, the reliability of such paleolatitudes can be markedly improved by using comparisons with paleolatitudes of different ages from the same site, paleolatitudes of similar ages from different sites on the same plate, estimates of paleolatitude from the skewness of marine magnetic anomalies, and continental paleopole. positions. Using such comparisons, a new paleomagnetic pole of upper Cretaceous age has been defined for the Pacific plate. A middle Cretaceous pole has been defined for the Wharton Basin plate, and it suggests that there may have been left lateral motion between Australia and the Wharton Basin. Paleolatitudes from the Ninetyeast Ridge are consistent with the pole position for the Deccan Traps. These data indicate that India and the Ninetyeast Ridge moved northwards with respect to the South Pole at 14.9 ± 4.5 cm/yr from 70 to 40 mybp and at 5.2 ± .8 cm/yr from 40 mybp until the present. However, when this paleomotion is compared to the Australian paleomagnetic data (by removing the relative motion components), a major inconsistency appears between 40 and 50 mybp. The Australian data indicate that India should be 13° further north than the positions implied by the Ninetyeast Ridge data. Basal paleolatitudes on the Ninetyeast Ridge indicate that its volcanic source was approximately fixed in latitude near 50°S, supporting the hypothesis that the ridge is the trace of the Kerguelen hotspot on the northward moving Indian plate. There is considerable geologic evidence in favor of such an hypothesis, and there is none to contradict it.

Description: National Science Foundation (Grant DES-74-22552).

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, 1975

Description: The Mid-Atlantic Ridge is one of the most well known and yet poorly understood spreading centers in the world. A detailed investigation of the Mid-Atlantic Ridge crest near 37°N (FAMOUS) was conducted using a deeply towed instrument package. The objective was to study the detailed structure and spreading history of the Mid-Atlantic Ridge median valley, to explore the roles of volcanism and faulting in the evolution of oceanic crust, and to study the morphologic expression and structural history of the zone of crustal accretion. In addition, microearthquake surveys were conducted using arrays of free-floating hydrophones. The most recent expression of the accreting plate boundary in the Famous Rift is an alternating series of linear central volcanoes and depressions 1.5 km wide which lie within the inner floor. This lineament is marked by a sharp maximum in crustal magnetization only 2-3 km wide. Magnetic studies indicate that over 90% of the extrusive volcanism occurs within the rift inner floor, a zone 1 to 12 km wide, while volcanism is extremely rare in the rift mountains. Volcanoes created in the inner floor are transported out on, block faults, becoming a lasting part of the topography. Magnetic anomaly transition widths vary from 1 km to 8 km with time and appear to reflect a bi-stable median valley structure. The valley has either a wide inner floor and narrow terraces, in which case the volcanic zone is wide and magnetic anomalies are poorly recorded (wide transition widths); or it has a narrow inner floor and wide terraces, the volcanic zone is then narrow and anomalies are clearly recorded (narrow transition widths). The median valley of any ridge segment varies between these two structures with time. At present the. Famous Rift has a narrow inner floor and volcanic zone (1-3 km) while the south Famous Rift is at the opposite end of the cycle with a wide inner floor and volcanic zone (10-12 km). Over 95% of the large scale (〉2 km) relief of the median valley is accounted for by normal faults dipping toward the valley axis. Normal faulting along fault planes dipping away from the valley begins just past the outer walls of the valley. Outward facing normal faulting accounts for most of the decay of median valley relief in the rift mountains while crustal tilting accounts for less than 20%. The pattern of normal faulting creates a broad, undulating horst and graben relief. Volcanic features contribute little to the large scale relief, but contribute to the short wavelength (〈2km) roughness of the topography. Spreading in the Famous area is highly asymmetric with rates twice as high to the east as to the west. At 1.7 m.y.b.p. the sense of asymmetry reverses in direction with spreading faster to the west, resulting in a gross symmetry when averaged through time. The change in spreading asymmetry occurred in less than 0.15 m.y. Structural studies indicate that the asymmetric spreading is accomplished through asymmetric crustal extension as well as asymmetric crustal accretion. Spreading in the Famous area is 17° oblique. Even on a fine scale there is no indication of readjustment to an orthogonal plate boundary system. Spreading has been stably oblique for at least 6 m.y., even through a change in spreading direction. Magnetic studies reveal that the deep DSDP hole at site 332 was drilled into a magnetic polarity transition, and may have sampled rocks which recorded the earth i s field behavior during a reversal. The presence of negative polarity crust within the Brunhes normal epoch in the inner floor has been determined, and may be due to old crust left behind or recording of a geomagnetic field event. Crustal magnetization decays to lie of its initial value in less than 0.6 m.y. The rapid decay may be facillitated by very intense crustal fracturing observed in the inner floor. Microearthquake, magnetic and structural studies indicate that both the spreading and transform plate boundaries are very narrow (1-2 km) and well-defined for short periods, but migrate over zones 10-20 km wide through time.