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 Colson, B. C., & Michel, A. P. M. Flow-through quantification of microplastics using impedance spectroscopy. ACS Sensors, 6(1), (2021): 238–244, doi:10.1021/acssensors.0c02223.

Description: Understanding the sources, impacts, and fate of microplastics in the environment is critical for assessing the potential risks of these anthropogenic particles. However, our ability to quantify and identify microplastics in aquatic ecosystems is limited by the lack of rapid techniques that do not require visual sorting or preprocessing. Here, we demonstrate the use of impedance spectroscopy for high-throughput flow-through microplastic quantification, with the goal of rapid measurement of microplastic concentration and size. Impedance spectroscopy characterizes the electrical properties of individual particles directly in the flow of water, allowing for simultaneous sizing and material identification. To demonstrate the technique, spike and recovery experiments were conducted in tap water with 212–1000 μm polyethylene beads in six size ranges and a variety of similarly sized biological materials. Microplastics were reliably detected, sized, and differentiated from biological materials via their electrical properties at an average flow rate of 103 ± 8 mL/min. The recovery rate was ≥90% for microplastics in the 300–1000 μm size range, and the false positive rate for the misidentification of the biological material as plastic was 1%. Impedance spectroscopy allowed for the identification of microplastics directly in water without visual sorting or filtration, demonstrating its use for flow-through sensing.

Description: The authors thank the Richard Saltonstall Charitable Foundation and the National Academies Keck Futures Initiative (NAKFI DBS13) for their funding support.

Description: Submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics and Astronautics at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution June 2021.

Description: Global temperature rise and increased atmospheric carbon dioxide (CO2) levels have affected the health of the world’s ocean and water ecosystems, impacting the balances of natural carbon cycling and causing ocean acidification. Additionally, as global temperatures rise, thawing permafrost has stimulated increased release of methane (CH4), a gas with a shorter lifetime in the atmosphere but with even more heat trapping ability than CO2. In situ analysis of dissolved gas content in surface waters is currently performed with large, expensive instruments, such as spectrometers, which are coupled with gas equilibration systems, which extract dissolved gas from water and feed it to the sensor. Accurate, low cost, and portable sensors are needed to measure the dissolved CH4 and CO2 concentration in water systems to quantify their release and understand their relationship to the global carbon budget. At the same time, while greenhouse gases are well established threats to water ecosystems, the ubiquity and potential consequences of microplastics in aqueous environments are just beginning to be recognized by the environmental research community. Microplastics (MPs) are small particles of polymer debris, commonly defined as being between 1 μm and 1000 μm. Despite the pervasiveness of MPs, our ability to characterize MPs in the environment is limited by the lack of technologies for rapidly and accurately identifying and quantifying MPs. This thesis is concerned with the engineering challenges prompted by the need for high quality and quantity environmental data to better study and the impact, cycling, and prevalence of these pollutants in aqueous environments. Three distinct investigations are presented here. First, the design of the Low-Cost Gas Extraction and Measurement System (LC-GEMS) for dissolved CO2 is presented. At just under $600 dollar to build, the LC-GEMS is an ultra-portable, toolbox-sized instrument for dissolved gas sensing in near-surface waters. The LCGEMS was characterized in the lab and demonstrated linear relationships with dissolved CO2 as well as temperature. Lab calibrations and subsequent field testing in the Little Sippewissett Marsh, in Falmouth, Massachusetts showed that the LCGEMS captures both diurnal and minute-time scale trends in dissolved CO2. Second, this thesis presents the novel design of three simple and low-cost planar nanophotonic and plasmonic structures as optical transducers for measuring dissolved CH4. Through simulations, the sensitivity of the structures are evaluated and found to exhibit superior performance in the reflectance intensity readout mode to that of the standard surface-plasmon-polariton-mode Spreeta sensor. A practical, small, and low-cost implementation of this chip with a simple intensity-based measurement scheme is proposed. This design is novel in the space of dissolved gas monitoring because it shows potential to measure directly in the water phase while being robust and low-cost to implement. Finally, this thesis presents a literature review and perspective to motivate the development of field-deployable microplastic sensing techniques. A framework for field-deployable microplastic sensing is presented and seeks to inform the MP community of the potential in both traditional MP analysis techniques and unconventional methods for creating rapid and automated MP sensors. The field-deployabilty framework addresses a full scope of practical/technological trade-offs to be considered for portable MP detection.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 2019.

Description: This thesis addresses the dynamics of estuarine networks, based on hydrographic observations in Newark Bay, a sub-estuarine network connected to the Hudson River estuary through New York Harbor. Estuarine networks differ from simple estuaries in that they may have multiple connections to the ocean, multiple freshwater sources, and often contain complex junctions between estuarine segments. The Newark Bay estuarine network is connected to the sea through two tidal straits, and is fed by multiple internal and external sources of fresh water. The estuarine network is also naturally divided into a series of reaches, each of which is characterized by a different cross-sectional geometry. This thesis focuses on the hydrographic variability and varying exchange flow within the Newark Bay estuarine network. Shipboard hydrographic measurements reveal the time-dependent formation of salinity fronts between reaches of the estuary. Each front is generated by a different mechanism; however, all are generated by tidal flow through channel junctions during ebb tide, and are advected landward during flood tide. Mooring-based measurements confirm that these fronts form during nearly every tidal cycle, and that the fronts are associated with substantial changes in local salinity on tidal timescales. The effect of tidal processes, such as frontal advection, on the exchange flow is investigated by applying the isohaline total exchange flow (TEF) framework to mooring-based observations in multiple reaches of the estuarine network. This reveals that over half of the exchange flow is driven by tidal processes at all sites within the estuary. Both the TEF-based salt balance and the standard Eulerian salt balance indicate that tidal processes are also responsible for at least half of the landward salt flux at most sites within the estuary; TEF and Eulerian salt balances are nearly identical. Tidal processes within the estuary are in large part associated with fronts. The large influence of tidal processes on the exchange flow in Newark Bay is thus likely due to the prevalence of channel junctions within the estuarine network.

Description: The studies contained in this thesis were largely funded as part of a National Science Foundation Coastal SEES project (Grant OCE-1325136), which was developed to investigate the effects of anthropogenic modifications on the physical processes in estuaries. Additional funding was provided by the J. Seward Johnson Fund at Woods Hole Oceanographic Institution, and by Hudson River Foundation Graduate Fellowship GF/01/17.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Oceanography/Applied Ocean Science and Engineering at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution June 2019.

Description: Submesoscale ocean dynamics and instabilities, with characteristic scales 0.1–10 km, can play a critical role in setting the ocean’s surface boundary layer thickness and associated density stratification. Submesoscale instabilities contribute to lateral stirring and tracer dispersal. These dynamics are investigated in the Bay of Bengal, motivated by the upper ocean’s potentially coupled interactions with Monsoon winds and convection. The region’s excess precipitation and runoff generates strong salinity gradients that typically set density fronts and stratification in the upper 50 m. Since we cannot synoptically measure currents containing fast-evolving and oscillating components across the submesoscale range, we instead analyze passive tracer distributions (spice ⌘ density-compensated temperature (T) and salinity (S) anomalies), identifying signatures of flows and testing dynamical theories. The analysis is based on over 9000 vertical profiles of T and S measured along ⇠4800 km of ship tracks in the Bay of Bengal during ASIRI and MISO-BOB expeditions in 2013, 2015, and 2018. Observations in the surface mixed layer reveal ⇠1 km scale-selective correlation of surface T and S, with compensation reducing cross-front density gradients by ⇠50%. Using a process study ocean model, we show this is caused by submesoscale instabilities slumping fronts, plus surface cooling over the resultant enhanced salinity stratification, potentially thwarting the forward cascade of energy. In the stratified interior, we present a spectral analysis of horizontal spice variance statistics from wavenumber k ⇠0.01 cpkm to ⇠1 cpkm. At scales 〈10 km, stratified layers that are closer to the surface exhibit redder passive tracer spectra (power spectra k−3, gradient spectra k−1) than predicted by quasi-geostrophic or frontogenetic theories. Complimentary observations reveal spice patterns with multiple, parallel, ⇠10 m thin layers, crossing isopycnals with O(10−4) slopes, coherent over at least 30–80 km, with coincident layers of stratification anomalies. Comparison with shear measurements, and a numerical process study, suggest that both submesoscale sheared eddies, and thin near-inertial waves, form such layers. Fast formation timescales and large aspect ratios suggest they enhance horizontal mixing by shear dispersion, reducing variance at ⇠1–10 km scales.

Description: Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Physical Oceanography at the Massachusetts Institute of Technology and the Woods Hole Oceanographic Institution September 2019.

Description: The Galápagos Cold Pool (GCP) is a region of anomalously cold sea surface temperature (SST) just west of the Galápagos Archipelago. Modeling studies have shown that the GCP is maintained by wind- and current-driven upwelling. The Galápagos Archipelago lies on the equator, in the path of the Pacific Equatorial Undercurrent (EUC) as it flows eastward across the Pacific at the depth of the thermocline. It is hypothesized that the EUC upwells into the GCP as it reaches the topographical barrier of the Galápagos Archipelago. The path of the EUC in the vicinity of the archipelago is not well understood. The ‘Repeat Observations by Gliders in the Equatorial Region’ (ROGER) program deployed a fleet of Spray autonomous underwater gliders in the region just west of the Galápagos Archipelago from 2013 – 2016 with the goal of continuously occupying three transects that form a closed area, with the archipelago as the eastern boundary. Gliders obtained subsurface measurements of temperature, salinity, and velocity with unprecedented temporal and spatial resolution. These measurements are used to observe the path of the EUC as it bifurcates into a north and south branch around the Galápagos Archipelago. Net horizontal transport into the volume defined by the closed area formed by the glider transects is used to estimate an average vertical velocity profile in the region of the GCP, indicating upwelling in the upper 300 m. The bifurcation latitude of the EUC, estimated to be approximately 0.4∘S from volume transport as a function of salinity, is coincident with the meridional center of the archipelago, suggesting the bifurcation latitude is topographically controlled. Ertel potential vorticity and a Bernoulli function are qualitatively conserved, supporting an inertial model of the EUC. Average spectral variance from Argo profiling float observations is used to show that tropical instability waves propagate with frequency and wavelength consistent with linearized, equatorial 𝛽-plane model results and may impact the GCP, according to their vertical structure.

Description: Funding for my thesis research was provided by the National Science Foundation (grants OCE-1232971 and OCE-1233282), the NASA Earth and Space Science Fellowship Program (grant 80NSSC17K0443), the J. Seward Johnson Fund, and the Karen L. Von Damm Fellowship.

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 2016

Description: Eddies in the ocean move westwards. Those shed by western boundary currents must then interact with continental shelf-slope topography at the western boundary. The presence of other eddies and mean lows complicates this simple picture, yet satellite images show that mesoscale eddies translating near the shelfbreak routinely affect the continental shelves of the Mid-Atlantic Bight, the Gulf of Mexico etc. The consequent cross-shelfbreak transports are currently of unknown importance to shelf budgets of heat, salt and volume. Thus motivated, this thesis uses idealized continuously stratified numerical experiments to explore eddy-slope interactions under four questions: 1. Can the continental slope prevent an eddy from reaching the shelfbreak? 2. What is the structure of the eddy-driven offshore low? 3. How is the continental shelf affected by an eddy at the shelfbreak? 4. Given surface observations, can one estimate the volume of water transported across the shelfbreak? The experiments show that the efficiency of Rossby wave radiation from the eddy controls whether it can cross isobaths: by radiating energy the eddy becomes shallow enough to move into shallower water. For wide continental slopes, relative to an eddy diameter, a slope can prevent an anticyclone from reaching the shelfbreak by shutting down such radiation. For narrow continental slopes, the interaction repeatedly produces dipoles, whose cyclonic halves contain shelf-slope water stacked over eddy water. The formation of such cyclones is explained. Then, the structure of shelf lows forced by the eddy are studied: their vertical structures are rationalized and scalings derived for their cross-isobath scales; for example, the extent to which the eddy influences the shelf. A recipe for estimating cross-isobath transports based on eddy surface properties is put forward. Finally, the findings are tested against observations in the Middle Atlantic Bight of the northeastern United States.

Description: I acknowledge high-performance computing support from Yellowstone (ark:/85065/d7wd3xhc) provided by NCAR’s Computational and Information Systems Laboratory, sponsored by the National Science Foundation. The research presented here was funded by National Science Foundation grants OCE-1059632 and OCE-1433953. Funding support from the Academic Programs Oice, WHOI is also gratefully acknowledged.

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 February 2013

Description: The interpretation of echoes collected by active remote-sensing systems, such as sonar and radar, is often ambiguous due to the complexities in the scattering processes involving the scatterers, the environment, and the sensing system. This thesis addresses this challenge using a combination of laboratory and fi eld experiments, theoretical modeling, and numerical simulations in the context of acoustic scattering by marine organisms. The unifying themes of the thesis are 1) quantitative characterization of the spectral, temporal, and statistical features derived from echoes collected using both broadband and narrowband signals, and 2) the interpretation of echoes by establishing explicit links between echo features and the sources of scattering through physics principles. This physics-based approach is distinct from the subjective descriptions and empirical methods employed in most conventional fisheries acoustic studies. The fi rst part focuses on understanding the dominant backscattering mechanisms of live squid as a function of orientation. The study provides the first broadband backscattering laboratory data set from live squid at all angles of orientation, and conclusively con firms the fluidlike, weakly-scattering material properties of squid through a series of detailed comparisons between data and predictions given by models derived based on the distorted-wave Born approximation. In the second part, an exact analytical narrowband model and a numerical broadband model are developed based on physics principles to describe the probability density function of the amplitudes of echo envelopes (echo pdf) of arbitrary aggregations of scatterers. The narrowband echo pdf model signi cantly outperforms the conventional mixture models in analyzing simulated mixed assemblages. When applied to analyze fish echoes collected in the ocean, the numerical density of sh estimated using the broadband echo pdf model is comparable to the density estimated using echo integration methods. These results demonstrate the power of the physics-based approach and give a rst-order assessment of the performance of echo statistics methods in echo interpretation. The new data, models, and approaches provided here are important for advancing the eld of active acoustic observation of the ocean.

Description: Taiwan Merit Scholarship (NSC-095-SAF-I-564-021-TMS), Office of Naval Research (ONR; grants N00014-10-1-0127, N00014-08-1-1162, N00014-07-1-1034), National Science Foundation (NSF; grant OCE-0928801), Naval Oceanographic Offi ce (grant N62306007-D9002), WHOI Ocean Life Institute, and the WHOI Academic Programs O ffice funds.

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

Description: A portion of the northeast Pacific ocean was chosen within which to evaluate and use altimetric data from the U.S. Navy Geodetic Satellite GEOSAT. The zero-order accuracy of the major GEOSAT geophysical data record (GDR) channels was verified, and occasional gaps in the altimeter coverage were noted. GEOSAT'S 17-day repeat orbit allowed use of collinear-track processing to create profiles of the difference between the sea surface height along a given satellite repeat, and the mean sea surface height along that repeat's groundtrack. Detrending of sea surface bias and tilt on each repeat reduced orbit and other long wavelength errors in the difference profiles. The corrections provided on the GEOSAT GDR were examined for their effects on the difference profiles of three test arcs. It was found that only the ocean tide, electromagnetic bias, and inverted barometer corrections varied enough over the arc lengths (~4400 km) to have any noticeable effect on the difference profiles. Only the ocean tide correction was accurate enough to warrant using it to adjust the sea surface heights. The recommended processing of GEOSAT data for the area included making the ocean tide correction, three-point block averaging successive sea surface heights, and forming the mean height profiles from 18 repeat cycles (to reduce aliasing of the M2 tidal component). A set of difference proftles for one GEOSAT arc indicated that a reasonable estimate of GEOSAT's system precision was -4.5 cm (RMS). The mid wavelength range (100-500 km) of these profiles was found to be the only range in which oceanic mesoscale features could be separated from altimeter errors. Mean alongtrack wavenumber spectra of oceanic variability for two GEOSAT arcs were compared with a SEASAT-derived regional spectrum of Fu (1983). Agreement was good, with GEOSAT showing less system noise at short wavelengths, and greater oceanic variability at long wavelengths. The GEOSAT spectra were fit well by a k-1.5 slope at wavelengths from 100 to 1000 km. Sea surface temporal variability as a function of location was compared with the SEASAT results of Cheney et al. (1983). Qualitative agreement was excellent and quantitative differences were largely accounted for. GEOSAT picked up the variability of the major current systems of the northeast Pacific, including the Alaskan, Californian, and North Equatorial currents. Error bounds on GEOSAT-derived oceanic variability showed that the effects of uncorrected electromagnetic bias, inverted barometer, and wet troposphere were significant. Further work in the areas of error modelling, orbit determination, and geoid calculation were called for.

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

Description: This thesis examines the nature of eddy-mean flow interactions in western boundary current jets and recirculation gyre dynamics from both theoretical and observational perspectives. It includes theoretical studies of eddy-mean flow interactions in idealized configurations relevant to western boundary current jet systems, namely (i) a study of the mechanism by which eddies generated from a localized forcing drive mean recirculation gyres through the process of nonlinear rectification; and (ii) a study of the role of eddies in the downstream evolution of a baroclinic jet subject to mixed instabilities. It also includes an observational analysis to characterize eddy-mean flow interactions in the Kuroshio Extension using data from the downstream location of maximum eddy kinetic energy in the jet. New insights are presented into a rectification mechanism by which eddies drive the recirculation gyres observed in western boundary current systems. Via this mechanism, eddies drive the recirculations by an up-gradient eddy potential vorticity flux inside a localized region of eddy activity. The effectiveness of the process depends on the properties of the energy radiation from the region, which in turn depends on the population of waves excited. In the zonally-evolving western boundary current jet, eddies also act to stabilize the unstable jet through down-gradient potential vorticity fluxes. In this configuration, the role of eddies depends critically on their downstream location relative to where the unstable time-mean jet first becomes stabilized by the eddy activity. The zonal advection of eddy activity from upstream of this location is fundamental to the mechanism permitting the eddies to drive the mean flows. Observational results are presented that provide the first clear evidence of a northern recirculation gyre in the Kuroshio Extension, as well as support for the hypothesis that the recirculations are, at least partially, eddy-driven. Support for the idealized studies’ relevance to the oceanic regime is provided both by indications that various model simplifications are appropriate to the observed system, as well as by demonstrated consistencies between model predictions and observational results in the downstream development of time-mean and eddy properties.

Description: Funding was for this research and my education was provided by the MIT Presidential Fellowship and NSF grants OCE-0220161 and OCE-0825550. The financial assistance of the Houghton Fund, the MIT Student Assistance Fund, and WHOI Academic Programs is also gratefully acknowledged.

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 2008

Description: The subtidal circulation of the southeast Greenland shelf is described using a set of highresolution hydrographic and velocity transects occupied in summer 2004. The main feature present is the East Greenland Coastal Current (EGCC), a low-salinity, highvelocity jet with a wedge-shaped hydrographic structure characteristic of other surface buoyancy-driven currents. The EGCC was observed along the entire Greenland shelf south of Denmark Strait, while the transect north of the strait showed only a weak shelf flow. This observation, combined with evidence from chemical tracer measurements that imply the EGCC contains a significant Pacific Water signal, suggests that the EGCC is an inner branch of the polar-origin East Greenland Current (EGC). A set of idealized laboratory experiments on the interaction of a buoyant current with a submarine canyon also supported this hypothesis, showing that for the observed range of oceanic parameters, a buoyant current such as the EGC could exhibit both flow across the canyon mouth or into the canyon itself, setting the stage for EGCC formation. Repeat sections occupied at Cape Farewell between 1997 and 2004 show that the alongshelf wind stress can also have a strong influence on the structure and strength of the EGCC and EGC on timescales of 2-3 days. Accounting for the wind-induced effects, the volume transport of the combined EGC/EGCC system is found to be roughly constant (~2 Sv) over the study domain, from 68°N to Cape Farewell near 60°N. The corresponding freshwater transport increases by roughly 60% over this distance (59 to 96 mSv, referenced to a salinity of 34.8). This trend is explained by constructing a simple freshwater budget of the EGCC/EGC system that accounts for meltwater runoff, melting sea-ice and icebergs, and net precipitation minus evaporation. Variability on interannual timescales is examined by calculating the Pacific Water content in the EGC/EGCC from 1984-2004 in the vicinity of Denmark Strait. The PW content is found to correlate significantly with the Arctic Oscillation index, lagged by 9 years, suggesting that the Arctic Ocean circulation patterns bring varying amounts of Pacific Water to the North Atlantic via the EGC/EGCC.

Description: Funding for the cruise and analysis was provided by National Science Foundation grant OCE-0450658, which along with NSF grant OCE- 0095427 provided funds for my tuition and stipend as well.

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 August 1980

Description: Observational evidence of seasonal variability below the main thermocline in the eastern North Atlantic is described, and a theoretical model of oceanic response to seasonally varying windstress forcing is constructed to assist in the interpretation of the observations. The observations are historical conductivity-temperature-depth data from the Bay of Biscay region (2° to 20°W, 42° to 52°N), a series of eleven cruises over the three years 1972 through 1974, spaced approximately three months apart. The analysis of the observations utilizes a new technique for identifying the adiabatically leveled density field corresponding to the observed density field. The distribution of salinity anomaly along the leveled surfaces is examined, as are the vertical displacements of observed density surfaces from the leveled reference surfaces, and the available potential energy. Seasonal variations in salinity anomaly and vertical displacement occur as westward propagating disturbances with zonal wavelength 390 (±50) km, phase 71 (±30) days from 1 January, and maximum amplitudes of ±30 ppm and ±20 db respectively. The leveled density field varies seasonally with an amplitude corresponding to a thermocline displacement of ±15 db. The observations are consistent with the predictions of a model in which an ocean of variable stratification with a surface mixed layer and an eastern boundary is forced by seasonal changes in a sinusoidal windstress pattern, when windstress parameters calculated from the observations of Bunker and Worthington (1976) are applied.

Description: This work was supported by the Office of Naval Research under contract N00014~76-C-197, NR 083-400.

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 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: The series of observations described in this report were planned with the double purpose of measuring the evaporation and transport of water vapor from the ocean into an unstable atmosphere, and of studying the diffusion processes operating in air of this stability class. Measured values of the evaporation from ocean surfaces were conspicuously absent from the meteorological literature until Craig and Montgomery (1949) published values for hydrostatically stable air. The present set of measurements extends our knowledge to include evaporation into a hydrostatically unstable air mass. In addition to evaporation values at the surface, net transports of water vapor at many levels up to 2000 meters have been measured.

Description: With the advance in knowledge of oceanic circulation there now exists a demand for additional identifying properties which will serve to trace the origin and movements of water masses in the sea, and to check earlier conclusions based entirely on temperature and salinity distribution. Of all the known identifying properties (except temperature and salinity) oxygen appears to be the most useful, not only because of the ease with which it can be accurately measured at sea, but also because of the large amount of data available on its distribution in the open ocean.

Description: With the opening of the Woods Hole Oceanographic Institution in July 1931, there was inaugurated a program of investigations in the deep waters of the western North Atlantic for which there had long been a great need. In contrast to most deep-sea investigations, which have had to be planned as single expeditions, the Institution was able to initiate a general program which could be carried out gradually in order to take advantage of knowledge gained during the course of the work. Suffcient funds having been provided for the continuous operation of its research vessel "Atlantis," work could be planned for all seasons of the year. Although these investigations have not been in progress for long and new data are continually being brought in by the "Atlantis," there are several reasons that make it seem desirable at this time to publish a preliminary report based on the completed temperature and salinity observations. In the first place, the problem of oceanic circulation is such that we cannot hope for a satisfactory solution for a long time to come. Moreover, it would be unwise to allow too much data to accumulate, because several years may pass before we can arrive at more important conclusions. Secondly, both the chemical and biological programs undertaken at the same time, require as a background the general scheme of circulation in the western North Atlantic as well as the distribution of temperature and salinity. It is, in fact, the necessity of taking into consideration the movements of the sea water which ties together the whole subject of oceanography. Therefore, it is the duty of those interested in ocean circulation to make available their findings as soon as possible for investigators of other problems in the same area. The "Atlantis" temperature and salinity observations discussed in these pages were planned with two main purposes in view. The first objective was an intensive study of seasonal changes along sections running from the southwestern corner of Nova Scotia to Bermuda and from Bermuda to the mouth of Chesapeake Bay.! This, of course, included an examination of fluctuations in the Gulf Stream, as well as of the variations in the water masses on each side of it. Second, there has been planned and partly carried out, a more general survey of the western North Atlantic, where accurate, deep stations have been sadly lacking.

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: During the past four years a deliberate effort has been made at the Woods Hole Oceanographic Institution to devise methods of kinematic observation generally suited to the needs of oceanographers. One result of this work, the electromagnetic method, has been brought from the experimental stage to one of useful maturity. Many of the theoretical potentialities of the method are still to be explored and developed. Nevertheless it seems likely that this remaining work may be done more soundly if present developments of the theory and instrumentation are made available for use and evaluation by, others. These studies in methods of kinematic observation have been supported mainly under the provisions of Bureau of Ships Contract NObs-2083, and Office of Naval Research Contract N6onr-277-1. This support and the assistance of the Naval Ordnance Laboratory, the Hydrographic Office (Oceanographic Division), the United States Coast Guard, and the David Taylor Model Basin of the United States Navy is gratefully acknowledged.

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.

Description: The New York Bight consists of the waters lying between Cape May, New Jersey, and Montauk Point, Long Island. A portion of the general southwesterly current known as the Coastal Drift lies in the seaward part of the Bight. Inshore from the Coastal Drift is an area of complex hydrography where the combined outflows of the Hudson River and other rivers enter the sea. In the region where the New Jersey and Long Island coastlines converge, an area 25 nautical miles on each side has been studied at all seasons of the year. This area extends from Sandy Hook southward to a point off Seaside Heights, eastward to 73°15' W longitude, north to the Long Island shore, and westward to Rockaway Inlet. The depth of water in the area averages about 90 feet, except in the innermost part of the Hudson Canyon which runs roughly northwest-southeast across most of the survey area. In the Canyon, depths in excess of 240. feet are found within the limits of the area studied. The hydrographic conditions in the area are in essence similar to those off the mouths of other large rivers. The combined flows of the Hudson and other rivers entering the surveyed area discharge enough fresh water annually to replace about one-half of the total volume of water under the 600 square miles of sea surface extensively surveyed. The salinity within the area is nearly as high as that of adjacent coastal water, however, and the actual quantity of river water within the area at any time rarely exceeds one percent of the total volume of water. Quantitative evaluation of these factors has led to the conclusion that there is an active circulation within the area which rapidly disperses the introduced river effluent. Many surveys of coastal and estuarine waters have been made. Outstanding among these are the survey of the River Tees, (1931, 1935), of the Tamar Estuary, (Hartley and Spooner, 1938; Milne, 1938), and of Alberni Inlet, (Tully, 1949). The general principles of estuarine circulations may be summarized as follows: In order to remove the added river water there must be a non-tidal drift of mixed water in a net seaward direction. When river flow remains constant, a steady state distribution of fresh and salt water throughout the estuary is attained, and at such times the net transport of river water seaward through any complete cross section of the estuary exactly equals the contribution of fresh water from the river during the same interval of time. As the mixture containing the river water moves seaward it gets progressively more saline, as additional sea water is entrained. In order to provide this sea water there must be a counter drift having a net flow in a landward direction. Superimposed on these necessary parts of the circulation are tidal and wind currents. The velocities of the tidal currents are commonly much greater than the velocity of the non-tidal drift, making the latter difficult to measure directly. It can be inferred, however, from the distribution of river water, as derived from the salinity distribution. Using the river water in this way we have evaluated the exchanges of the waters within the New York Bight. Tully (1949) has analyzed the circulation in Alberni Inlet by similar methods. Tidal current measurements made by the Coast and Geodetic Survey at various locations in the northwestern corner of the surveyed area are summarized by Marmer (1935). At Scotland Lightship, which is the location of the stations at the western end of Section A in Figure 1, the total excursion which results from the flood or ebb tidal currents is less than two miles. The currents at Ambrose Lightship, about five miles to the eastward, produce displacements only about half as great. The tidal displacements throughout the rest of the area are presumed to be less than these. The pattern of distribution of properties will be displaced, therefore, a distance less than ±1 mile at various stages of the tide. This distance is small in comparison to the size of the area surveyed, especially when considering the fact that distances between stations ranged from 5 to 8 miles. It was unnecessary, therefore, to attempt to take comparable stations at similar stages of the tide. Other considerations, beside its interesting hydrography, have contributed to the choice of this area for study. Because it is adjacent to centers of dense population and heavy industrial concentration, the New York Bight serves the conflicting purposes of waste disposal and recreation. Sewer effluents and industrial wastes enter the area by way of the rivers. Sewage sludges are barged out and dumped within the region studied. During the period covered by our surveys, The National Lead Company commenced operations to barge and discharge at sea the waste from its titanium plant at Sayreville, New Jersey. Since iron was a major constituent of this waste, analyses for iron in the water were made at each station, and the results have been valuable in checking the rate of the circulation which was computed from the distribution of river effluent. The New York Bight is also used extensively for recreational purposes. Because the area is readily and cheaply accessible by public transportation it must serve the recreational demands of a large part of the population of metropolitan New York. Sport fishing, bathing and boating are the principal recreational activities. Small but valuable commercial fisheries for shellfish and fin-fish also exist. The purpose of this study was to investigate the hydrographic processes in the New York Bight since they have an important bearing on the general problems of coastal oceanography and a knowledge of them should lead to a more successful evaluation and utilization of the area for the diverse purposes it must serve.

Description: This paper is a discussion of possible discrepancies in computations of ocean currents (based on horizontal variations of dynamic topography calculated from arbitrary deep lying reference surfaces), because of time variations of temperature and salinity at fixed depths in the sea (illustrated for a 24-hour period at "Atlantis" Station 2639). The results contained herein, while based chiefly on information from the western North Atlantic, are of general applicability, since time variations of the same order of magnitude have been observed over extensive areas of the Atlantic ocean. In selecting material for analysis of dynamic situations in the region concerned, consideration has been given only to those favorably located stations from which the structural features could most conveniently be obtained for illustrating the points in question.