ALBERT

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

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

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

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

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

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

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

Description: Author Posting. © American Geophysical Union, 2012. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 13 (2012): Q08014, doi:10.1029/2012GC004163.

Description: Mapping and sampling of 18 eruptive units in two study areas along the Galápagos Spreading Center (GSC) provide insight into how magma supply affects mid-ocean ridge (MOR) volcanic eruptions. The two study areas have similar spreading rates (53 versus 55 mm/yr), but differ by 30% in the time-averaged rate of magma supply (0.3 × 106 versus 0.4 × 106 m3/yr/km). Detailed geologic maps of each study area incorporate observations of flow contacts and sediment thickness, in addition to sample petrology, geomagnetic paleointensity, and inferences from high-resolution bathymetry data. At the lower-magma-supply study area, eruptions typically produce irregularly shaped clusters of pillow mounds with total eruptive volumes ranging from 0.09 to 1.3 km3. At the higher-magma-supply study area, lava morphologies characteristic of higher effusion rates are more common, eruptions typically occur along elongated fissures, and eruptive volumes are an order of magnitude smaller (0.002–0.13 km3). At this site, glass MgO contents (2.7–8.4 wt. %) and corresponding liquidus temperatures are lower on average, and more variable, than those at the lower-magma-supply study area (6.2–9.1 wt. % MgO). The differences in eruptive volume, lava temperature, morphology, and inferred eruption rates observed between the two areas along the GSC are similar to those that have previously been related to variable spreading rates on the global MOR system. Importantly, the documentation of multiple sequences of eruptions at each study area, representing hundreds to thousands of years, provides constraints on the variability in eruptive style at a given magma supply and spreading rate.

Description: This work was supported by the National Science Foundation grants OCE08–49813, OCE08–50052, and OCE08– 49711.

Description: 2013-02-25

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 2012

Description: The structure of the oceanic lithosphere results from magmatic and extensional processes taking place at mid-ocean ridges (MORs). The temporal and spatial scales of the variability of these two processes control the degree of heterogeneity of the oceanic lithosphere, represented by two end-member models: the classical Penrose Model exemplified by layered magmatic crust formed along fast-spreading MORs, e.g., East Pacific Rise (EPR); and the recently defined Chapman Model describing heterogeneous mafic and ultramafic lithosphere formed in settings of oceanic detachment faulting common along slow-spreading MORs, e.g., Mid-Atlantic Ridge (MAR). This thesis is using advanced marine geophysical methods (including finite-difference wave propagation modeling, 3D multi-channel seismic reflection imaging, waveform inversion, streamer tomography, and near-bottom magnetics) to study lithospheric accretion processes in MORs characterized by contrasting tectono-magmatic settings: the magmatically dominated EPR axis between 9°30'-10°00’N, and the Kane Oceanic Core Complex (KOCC), a section of MAR lithosphere (23°20’-23°38’N) formed by detachment faulting. At the EPR study area, I found that the axial magma chamber (AMC) melt sill is segmented into four prominent 2-4-km-long sections spaced every ~5- 10 km along the ridge axis characterized by high melt content (〉95%). In contrast, within the intervening sections, the AMC sill has a lower melt content (41-46%). The total magma volume extracted from the AMC sill was estimated of ~46 × 106 m3, with ~24 × 106 m3 left unerupted in the upper crust as dikes after 2005-06 eruption. At the KOCC, I used streamer tomography to constrain the shallow seismic velocity structure. Lithological interpretation of the seismic tomographic models provides insights into the temporal and spatial evolution of the melt supply at the spreading axis as the KOCC formed and evolved. Investigation of a magnetic polarity reversal boundary in crosssection at the northern boundary of KOCC suggests that the boundary (representing both a frozen isotherm and an isochron) dips away from the ridge axis along the Kane transform fault scarp, with a west-dipping angle of ~45° in the shallow (〈1 km) crust and 〈20° in the deeper crust.

Description: This thesis was funded by National Science Foundation grants OCE-9987004, OCE- 0621660 and OCE-0327885, WHOI Academic Program Office and WHOI Deep Ocean Exploration Institute.

Description: Author Posting. © American Geophysical Union, 2012. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 13 (2012): Q04007, doi:10.1029/2011GC003990.

Description: Here we demonstrate with a study of the Lucky Strike hydrothermal field that image mosaicing over large seafloor areas is feasible with new image processing techniques, and that repeated surveys allow temporal studies of active processes. Lucky Strike mosaics, generated from 〉56,000 images acquired in 1996, 2006, 2008 and 2009, reveal the distribution and types of diffuse outflow throughout the field, and their association with high-temperature vents. In detail, the zones of outflow are largely controlled by faults, and we suggest that the spatial clustering of active zones likely reflects the geometry of the underlying plumbing system. Imagery also provides constraints on temporal variability at two time-scales. First, based upon changes in individual outflow features identified in mosaics acquired in different years, we document a general decline of diffuse outflow throughout the vent field over time-scales up to 13 years. Second, the image mosaics reveal broad patches of seafloor that we interpret as fossil outflow zones, owing to their association with extinct chimneys and hydrothermal deposits. These areas encompass the entire region of present-day hydrothermal activity, suggesting that the plumbing system has persisted over long periods of time, loosely constrained to hundreds to thousands of years. The coupling of mosaic interpretation and available field measurements allow us to independently estimate the heat flux of the Lucky Strike system at ~200 to 1000 MW, with 75% to 〉90% of this flux taken up by diffuse hydrothermal outflow. Based on these heat flux estimates, we propose that the temporal decline of the system at short and long time scales may be explained by the progressive cooling of the AMC, without replenishment. The results at Lucky Strike demonstrate that repeated image surveys can be routinely performed to characterize and study the temporal variability of a broad range of vent sites hosting active processes (e.g., cold seeps, hydrothermal fields, gas outflows, etc.), allowing a better understanding of fluid flow dynamics from the sub-seafloor, and a quantification of fluxes.

Description: This project was funded by CNRS/IFREMER through the 2006, 2008, 2009 and 2010 cruises within the MoMAR program (France), by ANR (France) Mothseim Project NT05-3 42213 to J. Escartín, and by grant CTM2010-15216/MAR from the Spanish Ministry of Science to R. Garcia and J. Escartín. T. Barreyre was supported by University Paris Diderot (Paris 7– France) and Institut de Physique du Globe de Paris (IPGP, France). E. Mittelstaedt was supported by the International Research Fellowship Program of the U.S. National Science Foundation (OISE-0757920).

Description: 2012-10-19

Description: Author Posting. © American Geophysical Union, 2011. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 12 (2011): Q05008, doi:10.1029/2010GC003439.

Description: Variations in topography and seismic structure are observed along the Juan de Fuca (JdF) Ridge axis and in the vicinity of pseudofaults on the ridge flanks left by former episodes of ridge propagation. Here we analyze gravity data coregistered with multichannel seismic data from the JdF Ridge and flanks in order to better understand the origin of crustal structure variations in this area. The data were collected along the ridge axis and along three ridge-perpendicular transects at the Endeavor, Northern Symmetric, and Cleft segments. Negative Mantle Bouguer anomalies of −21 to −28 mGal are observed at the axis of the three segments. Thicker crust at the Endeavor and Cleft segments is inferred from seismic data and can account for the small differences in axial gravity anomalies (3–7 mGal). Additional low densities/elevated temperatures within and/or below the axial crust are required to explain the remaining axial MBA low at all segments. Gravity models indicate that the region of low densities is wider beneath the Cleft segment. Gravity models for pseudofaults crossed along the three transects support the presence of thinner and denser crust within the pseudofault zones that we attribute to iron-enriched crust. On the young crust side of the pseudofaults, a 10–20 km wide zone of thicker crust is found. Reflection events interpreted as subcrustal sills underlie the zones of thicker crust and are the presumed source for the iron enrichment.

Description: This work was supported by the National Science Foundation grants OCE‐0648303 to Lamont‐Doherty Earth Observatory, OCE‐0648923 to 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 September 1998

Description: The upper ocean crust contains a comprehensive record of the shallow geological processes active along the world's mid-ocean ridge system. This thesis examines the magnetic and seismic structure of the upper crust at two contrasting ridges-the fast spreading East Pacific Rise (EPR) and the slow spreading Mid-Atlantic Ridge (MAR)-to build a more complete understanding about the roles of volcanic emplacement, tectonic disruption and hydrothermal alteration in the near-ridge environment. A technique that inverts potential field measurements. directly from an uneven observation track is developed and applied to near-bottom magnetic data from the spreading segments north of the Kane transform on the MAR. It is concluded that the central anomaly magnetization high marks the locus of focused volcanic emplacement. A cyclic faulting model is proposed to explain the oscillatory magnetization pattern associated with discrete blocks of crust being transported out of the rift valley between intensely altered fault zones. Seismic waveform and amplitude analyses of the magma sill along the EPR reveal it to be a thin (〈100 m) body of partial melt. These characteristics have important implications for melt availability and transport within the cycle of eruption and replenishment. A genetic algorithm-based seismic waveform inversion technique is developed and applied to on- and near-axis multichannel data from 17°20'S on the EPR and the spreading segment south of the Oceanographer transform (MAR) to map and compare for the first time the detailed velocity structure of the upper crust at two different spreading rates. Combined with conventionally processed seismic profiles, our results show that, while final extrusive thickness is comparable at all spreading ridges (300-500 m), the style of thickening may vary. While a thin (≤100 m) extrusive carapace quadruples in thickness within 1-4 km of the EPR crest, the extrusive section at the MAR achieves its final thickness within the inner valley. Both show evidence for a narrow zone of volcanic emplacement. Vigorous hydrothermalism at the EPR may produce a more rapid increase in basement velocities relative to the MAR. Rapid modification of the extrusive/dike transition at both ridges indicates that hydrothermalism is enhanced in this interval. Along-axis transport of lavas may thicken the extrusive pile at slow spreading segment ends, strengthening the magnetic highs generated by lava chemistry.

Description: ONR graduate fellow. The magnetics portion of this thesis was also supported by NSF grants OCE-9204141, OCE-9200905 (M. A. Tivey & H. Schouten) and NERC GR3/7702 (R. C. Searle). Research for Chapter 4 was partially funded by NSF grant OCE-9402933 (R. S. Detrick). The final two science chapters were supported by NSF grants OCE-9012707, OCE-9300450 (R. S. Detrick), OCE-9401717 (G. M. Kent & R. S. Detrick), OCE-9400623 (M. A. Tivey) and the Education office.

Description: Author Posting. © American Geophysical Union, 2011. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 38 (2011): L14307, doi:10.1029/2011GL047798.

Description: The heat output and thermal regime of fast and intermediate spreading centers are strongly controlled by boundary layer processes between the hydrothermal system and the underlying crustal magma chamber (AMC), which remain to be fully understood. Here, we model the interactions between a shallow two-dimensional cellular hydrothermal system at temperatures 〈700°C, and a deeper AMC at temperatures up to 1200°C. We show that hydrothermal cooling can freeze the AMC in years to decades, unless melt injections occur on commensurate timescales. Moreover, the differential cooling between upflow and downflow zones can segment the AMC into mush and melt regions that alternate on sub-kilometric length scales. These predictions are consistent with along-axis variations in AMC roof depth observed in ophiolites and oceanic settings. In this respect, fine-scale geophysical investigations of the structure of AMCs may help constrain hydrothermal recharge locations associated with active hydrothermal sites.

Description: Author Posting. © American Geophysical Union, 2011. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 116 (2011): B10102, doi:10.1029/2011JB008259.

Description: Crustal thickness variations at the ultraslow spreading 10–16°E region of the Southwest Indian Ridge are used to constrain melt migration processes. In the study area, ridge morphology correlates with the obliquity of the ridge axis with respect to the spreading direction. A long oblique “supersegment”, nearly devoid of magmatism, is flanked at either end by robust magmatic centers (Joseph Mayes Seamount and Narrowgate segment) of much lesser obliquity. Plate-driven mantle flow and temperature structure are calculated in 3-D based on the observed ridge segmentation. Melt extraction is assumed to occur in three steps: (1) vertical migration out of the melting region, (2) focusing along an inclined permeability barrier, and (3) extraction when the melt enters a region shallower than ∼35 km within 5 km of the ridge axis. No crust is predicted in our model along the oblique supersegment. The formation of Joseph Mayes Seamount is consistent with an on-axis melt anomaly induced by the local orthogonal spreading. The crustal thickness anomaly at Narrowgate results from melt extracted at a tectonic damage zone as it travels along the axis toward regions of lesser obliquity. Orthogonal spreading enhances the Narrowgate crustal thickness anomaly but is not necessary for it. The lack of a residual mantle Bouguer gravity high along the oblique supersegment can be explained by deep serpentization of the upper mantle permissible by the thermal structure of this ridge segment. Buoyancy-driven upwelling and/or mantle heterogeneities are not required to explain the extreme focusing of melt in the study area.

Description: This work was supported by grants OCE‐ 0623188 and OCE‐0937277 from the National Science Foundation.

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 2001

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

Description: Funding was provided by the National Science Foundation through grants OCE- 9811924 and OCE-9907630, and by a National Defense Science and Engineering Graduate Fellowship.

Description: Author Posting. © American Geophysical Union 2003. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 4 (2003): 8515, doi:10.1029/2003GC000609.

Description: Complete multibeam bathymetric coverage of the western Galápagos Spreading Center (GSC) between 90.5°W and 98°W reveals the fine-scale morphology, segmentation and influence of the Galápagos hot spot on this intermediate spreading ridge. The western GSC comprises three morphologically defined provinces: A Western Province, located farthest from the Galápagos hot spot west of 95°30′W, is characterized by an axial deep, rift valley morphology with individual, overlapping, E-W striking segments separated by non-transform offsets; A Middle Province, between the propagating rift tips at 93°15′W and 95°30′W, with transitional axial morphology strikes ∼276°; An Eastern Province, closest to the Galápagos hot spot between the ∼90°50′W Galápagos Transform and 93°15′W, with an axial high morphology generally less than 1800 m deep, strikes ∼280°. At a finer scale, the axial region consists of 32 individual segments defined on the basis of smaller, mainly 〈2 km, offsets. These offsets mainly step left in the Western and Middle Provinces, and right in the Eastern Province. Glass compositions indicate that the GSC is segmented magmatically into 8 broad regions, with Mg # generally decreasing to the west within each region. Striking differences in bathymetric and lava fractionation patterns between the propagating rifts with tips at 93°15′W and 95°30′W reflect lower overall magma supply and larger offset distance at the latter. The structure of the Eastern Province is complicated by the intersection of a series of volcanic lineaments that appear to radiate away from a point located on the northern edge of the Galápagos platform, close to the southern limit of the Galápagos Fracture Zone. Where these lineaments intersect the GSC, the ridge axis is displaced to the south through a series of overlapping spreading centers (OSCs); abandoned OSC limbs lie even farther south. We propose that southward displacement of the axis is promoted during intermittent times of increased plume activity, when lithospheric zones of weakness become volcanically active. Following cessation of the increased plume activity, the axis straightens by decapitating southernmost OSC limbs during short-lived propagation events. This process contributes to the number of right stepping offsets in the Eastern Province.

Description: This work was supported by NSF grants OCE98- 18632 to the University of Hawai’i and OCE98-19117 to the Woods Hole Oceanographic Institution; support was provided to M. B. by a CIW/DTM Postdoctoral Fellowship

Description: Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 9 (2008): Q04015, doi:10.1029/2007GC001611.

Description: Near-bottom magnetic data collected along the crest of the East Pacific Rise between 9°55′ and 9°25′N identify the Central Anomaly Magnetization High (CAMH), a geomagnetic anomaly modulated by crustal accretionary processes over timescales of ∼104 years. A significant decrease in CAMH amplitude is observed along-axis from north to south, with the steepest gradient between 9°42′ and 9°36′N. The source of this variation is neither a systematic change in geochemistry nor varying paleointensity at the time of lava eruption. Instead, magnetic moment models show that it can be accounted for by an observed ∼50% decrease in seismic Layer 2A thickness along-axis. Layer 2A is assumed to be the extrusive volcanic layer, and we propose that this composes most of the magnetic source layer along the ridge axis. The 9°37′N overlapping spreading center (OSC) is located at the southern end of the steep CAMH gradient, and the 9°42′–9°36′N ridge segment is interpreted to be a transition zone in crustal accretion processes, with robust magmatism north of 9°42′N and relatively low magmatism at present south of 9°36′N. The 9°37′N OSC is also the only bathymetric discontinuity associated with a shift in the CAMH peak, which deviates ∼0.7 km to the west of the axial summit trough, indicating southward migration of the OSC. CAMH boundaries (defined from the maximum gradients) lie within or overlie the neovolcanic zone (NVZ) boundaries throughout our survey area, implying a systematic relationship between recent volcanic activity and CAMH source. Maximum flow distances and minimum lava dip angles are inferred on the basis of the lateral distance between the NVZ and CAMH boundaries. Lava dip angles average ∼14° toward the ridge axis, which agrees well with previous observations and offers a new method for estimating lava dip angles along fast spreading ridges where volcanic sequences are not exposed.

Description: The research project was funded by National Science Foundation under grants OCE-9819261 and OCE- 0096468.

Description: Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 9 (2008): Q08O10, doi:10.1029/2008GC001965.

Description: We use 2-D numerical models to explore the thermal and mechanical effects of magma intrusion on fault initiation and growth at slow and intermediate spreading ridges. Magma intrusion is simulated by widening a vertical column of model elements located within the lithosphere at a rate equal to a fraction, M, of the total spreading rate (i.e., M = 1 for fully magmatic spreading). Heat is added in proportion to the rate of intrusion to simulate the thermal effects of magma crystallization and the injection of hot magma into the crust. We examine a range of intrusion rates and axial thermal structures by varying M, spreading rate, and the efficiency of crustal cooling by conduction and hydrothermal circulation. Fault development proceeds in a sequential manner, with deformation focused on a single active normal fault whose location alternates between the two sides of the ridge axis. Fault spacing and heave are primarily sensitive to M and secondarily sensitive to axial lithosphere thickness and the rate that the lithosphere thickens with distance from the axis. Contrary to what is often cited in the literature, but consistent with prior results of mechanical modeling, we find that thicker axial lithosphere tends to reduce fault spacing and heave. In addition, fault spacing and heave are predicted to increase with decreasing rates of off-axis lithospheric thickening. The combination of low M, particularly when M approaches 0.5, as well as a reduced rate of off-axis lithospheric thickening produces long-lived, large-offset faults, similar to oceanic core complexes. Such long-lived faults produce a highly asymmetric axial thermal structure, with thinner lithosphere on the side with the active fault. This across-axis variation in thermal structure may tend to stabilize the active fault for longer periods of time and could concentrate hydrothermal circulation in the footwall of oceanic core complexes.

Description: Funding for this research was provided by NSF grants OCE-0327018 (M.D.B.), OCE-0548672 (M.D.B.), OCE- 0327051 (G.I.), and OCE-03-51234 (G.I.).

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 June 2010

Description: Mid-ocean ridge basalts (MORBs) exhibit a wide range of CO2 concentrations, reflecting saturation to supersaturation (and rarely undersaturation) relative to their emplacement depths. In this study, we explore the mechanisms of CO2 degassing and the implications this has for estimating the advance rates and durations of seafloor eruptions. We present dissolved volatile concentrations (mainly of CO2 and H2O) and vesicle size distributions (VSDs) for a unique suite of MORB glasses collected at the East Pacific Rise, ~9° 50′ N. These MORB glasses were collected at ~200 m intervals along an across-axis track over a single flow pathway within the recently emplaced 2005-06 eruption boundaries; systematic sample collection provides one of the first opportunities to characterize intra-flow geochemical and physical evolution during a single eruption at a fast-spreading ridge. Compared to measurements of MORB volatiles globally, dissolved H2O concentrations are relatively uniform (0.10 - 0.16 weight percent), whereas dissolved CO2 contents exhibit a range of concentrations (154 - 278 ppm) and decrease with distance from the EPR axis (i.e., eruptive vent). Ion microprobe analyses of dissolved volatiles within the MORB glasses suggest that the magma erupted supersaturated (pressure equilibrium with 920 - 1224 mbsf) and in near-equilibrium with the melt lens of the axial magma chamber (~1250 - 1500 mbsf), and degassed to near equilibrium (299 - 447 mbsf) with seafloor depths over the length of the flow. The decrease in CO2 concentrations spans nearly the full range of dissolved CO2 contents observed at the EPR and shows that the varying degrees of volatile saturation that have been observed in other MORB sample suites may be explained by degassing during emplacement. Vesicularity (0.1 - 1.2%) increases with decreasing dissolved CO2 concentrations. We use vesicle size distributions (VSDs)—vesicle sizes and number densities—to quantify the physical evolution of the CO2 degassing process. VSDs suggest that diffusion of CO2 into preexisting vesicles, and not nucleation of new vesicles, is the dominant mechanism of increasing CO2 in the vapor phase. We also use VSDs, along with estimates of vesicle growth rates, to constrain emplacement time of the 2005-06 eruption to 〈~24 hours and to resolve variations in advance rate with downflow distance.

Description: Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry Geophysics Geosystems 9 (2008): Q08001, doi:10.1029/2007GC001922.

Description: Multichannel seismic and bathymetric data from the Juan de Fuca Ridge (JDFR) provide constraints on axial and ridge flank structure for the past 4–8 Ma within three spreading corridors crossing Cleft, Northern Symmetric, and Endeavour segments. Along-axis data reveal south-to-north gradients in seafloor relief and presence and depth of the crustal magma lens, which indicate a warmer axial regime to the south, both on a regional scale and within individual segments. For young crust, cross-axis lines reveal differences between segments in Moho two-way traveltimes of 200–300 ms which indicate 0.5–1 km thicker crust at Endeavour and Cleft compared to Northern Symmetric. Moho traveltime anomalies extend beyond the 5–15 km wide axial high and coincide with distinct plateaus, 32 and 40 km wide and 200–400 m high, found at both segments. On older crust, Moho traveltimes are similar for all three segments (∼2100 ± 100 ms), indicating little difference in average crustal production prior to ∼0.6 and 0.7 Ma. The presence of broad axis-centered bathymetric plateau with thickened crust at Cleft and Endeavour segments is attributed to recent initiation of ridge axis-centered melt anomalies associated with the Cobb hot spot and the Heckle melt anomaly. Increased melt supply at Cleft segment upon initiation of Axial Volcano and southward propagation of Endeavour segment during the Brunhes point to rapid southward directed along-axis channeling of melt anomalies linked to these hot spots. Preferential southward flow of the Cobb and Heckle melt anomalies and the regional-scale south-to-north gradients in ridge structure along the JDFR may reflect influence of the northwesterly absolute motion of the ridge axis on subaxial melt distribution.

Description: This work was supported by U.S. National Science Foundation grants OCE00-02488 to S.M.C., OCE06-48303 to S.M.C. and M.R.N., OCE-0648923 to J.P.C., and OCE00-02600 to G.M.K. and A.J.H.

Description: Author Posting. © American Geophysical Union, 2004. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 109 (2004): B10310, doi:10.1029/2004JB003066.

Description: Multichannel seismic reflection data are used to infer crustal accretion processes along the intermediate spreading Galapagos Spreading Center. East of 92.5°W, we image a magma lens beneath the ridge axis that is relatively shallow (1.0–2.5 km below the seafloor) and narrow (∼0.5–1.5 km, cross-axis width). We also image a thin seismic layer 2A (0.24–0.42 km) that thickens away from the ridge axis by as much as 150%. West of 92.7°W, the magma lens is deeper (2.5–4.5 km) and wider (0.7–2.4 km), and layer 2A is thicker (0.36–0.66 km) and thickens off axis by 〈40%. The positive correlation between layer 2A thickness and magma lens depth supports the interpretation of layer 2A as the extrusive volcanic layer with thickness controlled by the pressure on the magma lens and its ability to push magma to the surface. Our findings also suggest that narrower magma lenses focus diking close the ridge axis such that lava flowing away from the ridge axis will blanket older flows and thicken the extrusive crust off axis. Flow of lava away from the ridge axis is probably promoted by the slope of the axial bathymetric high, which is largest east of 92.5°W. West of ∼94°W the “transitional” axial morphology lacks a prominent bathymetric high and layer 2A no longer thickens off axis. We detect no magma lens west of 94.7°W where a small axial valley appears. The above changes can be linked to the westward decrease in the magma and heat flux associated with the fading influence of the Galapagos hot spot on the Galapagos Spreading Center.

Description: This project was funded by NSF-OCE- 0002189.