This thesis aims to investigate the stresses at mid-ocean ridge offsets, and particularly at the particular class of offsets represented by oceanic microplates. Amongthese, the Easter microplate is one of the best surveyed. This thesis first studies the stress field associated with mid-ocean ridges and simple types of ridge offsets, and then uses the stress field observed at Easter to constrain the driving mechanism of microplates. Two-dimensional finite element modelling is used to predict the lithospheric stress indicators, which are then compared with observations. Extensional structures at high angles (> 35 ) to ridge trends are often observed at ridge-transform intersections and non-tranform offsets, but remained unexplained until now. This study proposes that the topographic loading created by the elevation of mid-ocean ridges relative to old seafloor is a source of ridge parallel tensile stresses, and shows they can be explained by the rotation of ridge parallel tensile stresses at locked offsets. The elasto-plastic rheology is used to investigate the evolution of normal faults near mid-ocean ridges. It is shown that variations in the lithospheric strength, caused entirely by variations in the brittle layer thickness, can account for the observed variations in fault character with spreading rate and along-axis position. Plasticity is shown to prevent the achievement of large fault throws in thin brittle layers. Consequently, it may be important at fast spreading ridges. A new dynamic model is proposed for Easter microplate. It mainly consists of: 1) driving forces along the East and West Rifts, resulting from the combination of a regional tensile stress with an increasing ridge strength towards rift tips, 2) mantle basal drag resisting the microplate rotation, and contributing with less than 20% to the total resisting torque, and 3) resisting forces along the northern and southern boundaries. To explain both the earthquake focal mechanism evidence and theexistence of compressional ridges in the Nazca plate, the boundary conditions alongthe northern boundary are required to change with time, from completely locked tolocked in the normal direction only. This study does not invalidate the microplate kinematic model proposed by Schouten et al. (1993), but shows that normal resisting forces along the northern and southern boundaries of Easter microplate must exist in order to explain the stress observations. Also, it suggests that ridge strength variations play an important role in the dyamics of mid-ocean ridge overlap regions.

Hydrothermal circulation at mid-ocean ridges plays an important role in the interaction between oceanic lithosphere and the overlying ocean. Changes in fluid flux within hydrothermal systems may directly impact ocean circulation, temperature, and chemistry, and hence the lives of biological organisms in hydrothermal vent environments. The permeability structure within a hydrothermal environment is an important control on fluid flow throughout the system. Common suggestions for mechanisms that might increase permeability within the system include thermal cracking due to contraction of the hot rock from interaction with cold seawater, fluid pressure of the water moving through the rock, tectonics, or tidal forces. Additional factors such as mineral precipitation can decrease permeability in the system, further complicating the permeability structure. Though there are many factors to consider within a hydrothermal system, few quantitative studies of these cracking mechanisms exist. This study examines the role of thermal cracking near hydrothermal vents via a numerical model created in Matlab. Flow was modeled using the Laplace equation, and the heat transfer equation was used to determine temperature differences in the rock, which lead to thermal cracking. The numerical results were compared with microearthquakes observed by Tolstoy et al. (2008) near a hydrothermal vent field along the East Pacific Rise. The model suggests that thermal cracking does occur, and this cracking occurs within the area of microearthquakes observed. Though thermal cracking is important for increasing permeability within the system, there are no obvious spatial or temporal trends within the earthquake data that support a direct relation between the modeled thermal cracking and observed earthquakes. It is likely that the observed earthquakes are due to a combination of cracking mechanisms, such as cracking due to fluid pressure, tectonics, or tides, in addition to the modeled thermal stresses.

No description available.

At fast spreading ridges such as the East Pacific Rise (EPR) volcanic eruptions are predicted to occur on a decadal timescale. Due to the limited ability to observe submarine eruptions, little is known about the magmatic processes occurring in the underlying magma chamber leading up to a volcanic event, including differentiation and magma mixing. The recent 2005-06 eruption at 9˚46’N-9˚56’N along the EPR provides a unique opportunity to gain a better understanding of rates of magma transport and magma replenishment associated with a typical eruption. This study examines the geochemistry of phenocrysts from the 2005-06 eruption in order to determine if they are in equilibrium with their host melt, or if magma mixing occurred prior to eruption. A diffusion model is used to model those crystals which are out of equilibrium with their erupted host to determine timescales of magma mixing. The major and trace element contents of plagioclase and olivine phenocrysts provides evidence for melts both more evolved (> 3.5 wt % MgO) and more primitive (< 8.8 wt % MgO) than those found within the host lava (7.7-8.3 wt % MgO; Goss et al., 2010). Glomerocrysts and resorbed crystals in equilibrium with evolved melts (3.5-6.5 wt % MgO) suggests an origin in a roof mush zone, and were disrupted and entrained into their host melt within days of eruption. Modelling of the zoning profiles of phenocrysts suggest the 2005-06 eruption was likely triggered by an influx of hotter, more primitive melt (~ 9.0 wt % MgO) which was injected into the melt lens a few weeks to months prior to the eruption. With decreasing time before eruption, there is an overall increase in the number of crystals with modelled timescales representing mixing events in the magma chamber. This increase in modelled timescales appears to correlate with the increase in seismic activity recorded prior to the eruption (Tolstoy et al., 2006). This suggests magma mixing events within the underlying magma chamber may be linked to seismic activity at fast spreading ridges. / Graduate / 0372 / 0996

Seafloor hydrothermal systems have been intensively studied for the past few decades; however, the location of recharge zones and details of fluid circulation patterns are still largely uncertain. To better understand the effects of anhydrite precipitation on hydrothermal flow paths, we conduct 2-D numerical simulations of hydrothermal circulation at a mid-ocean ridge using a NaCl-H2O numerical code. The simulations focus on East Pacific Rise hydrothermal system at 950N due to availability of key observational data to constrain the models. Seismicity data that is available suggests that fluid flow is primarily along axis and that recharge is focused into a small zone near a 4th order discontinuity in the ridge axis. Simulations are carried out in an open-top square box 1500 m on a side maintained at a surface pressure of 25 MPa, and nominal seawater temperature of 10 C. The sides of the box are assumed to be impermeable and insulated. A constant temperature distribution is maintained along the bottom of the box consisting of a 1000 m long central-heated region maintained at 450 C to represent the axial magma chamber and ensure P-T conditions for phase separation; a linearly decreasing temperature profile from 450 to 300 C is maintained along the 250 m long segments adjacent to the heated region to delineate the recharge zone. We constructed a homogeneous model with a uniform cell size of 25 m with a permeability of 10-13 m2 and a similar model with a 200 m thick layer 2A region with a permeability of 10-12 m2. For the homogeneous model the simulations were run for 100 years to approximate steady state conditions and the model with layer 2A was run for 50 years. Assuming that anhydrite precipitation resulted from the decrease in solubility with increasing temperature as downwelling fluid gets heated, the rate of porosity decrease and sealing time was calculated at 50 and 100 years. The results showed that sealing occurred most rapidly at the bottom of the recharge areas near the base of the high-temperature plumes, where complete sealing occurred after ~55-625 years for an initial porosity of 0.1. The simulations also suggested that sealing would occur more slowly at the margins of the ascending plumes, with times ranging between ~ 80 and 5000 years. The sealing times in the deep recharge zone determined in these simulations are considerably greater than estimated from 1D analytical calculations, suggesting that with a 2D model, focused recharge at the EPR 950N site may occur, at least on a decadal time scale. More detailed analyses are needed to determine whether such focused recharge can be maintained for longer times. / Master of Science

This thesis discusses the results from the first multi-source and multi-streamer three-dimensional multichannel seismic experiment conducted over a mid-ocean ridge environment. Prestack time migration was applied to the dataset resulting in the most detailed reflection images of a spreading center and its flanks to date. The key products from this work are maps of crustal velocities, crustal thickness, and Moho transition zone (MTZ) reflection character for a section of the fast-spreading East Pacific Rise (EPR) from 9°37.5’N to 9°57’N, excluding the area from 9°40’N to 9°42’N where no data were collected. Moho reflections were imaged within ~92% of the study area. The derived average crustal thickness and average crustal velocity for the investigated ~880 km2 area are 5920±320 m and 6320±290 m/s, respectively. The average crustal thickness varies little from Pacific to Cocos plate suggesting mostly uniform crustal production in the last ~180 Ka. Detailed analysis of the crustal thickness and MTZ reflection character shows that the third-order segmentation is governed by melt extraction processes within the uppermost mantle while the fourth-order ridge segmentation arises from mid- to upper-crustal processes. This analysis also suggests that both the mechanism of lower-crustal accretion and the volume of melt delivered to the crust vary along the investigated section of the EPR. More efficient mantle melt extraction is inferred at latitudes from 9°42’N to 9°51.5’N, with greater proportion of the lower crust accreted from the AML than for the rest of the study area. Larger volume of melt is delivered to the crust from 9°37.5’N to 9°40’N than to the investigated crust further north. At some locations, the Moho reflections are for the first time unambiguously imaged below the AML away from any ridge discontinuity suggesting that the Moho is formed at zero age at least at some sections of the spreading centers. The first study of the melt content of mid-crustal off-axis magma lenses (OAML), done using amplitude variation with offset technique calibrated for a magmatic plumbing system, shows that these magma bodies contain 0 to 20% melt. This suggests that OAMLs likely contribute little to the overall crustal formation.

We present 2-D numerical simulations of two phase flow in seafloor hydrothermal systems using the finite control volume numerical scheme FISHES. The FISHES code solves the coupled non-linear equations for mass, momentum, energy, and salt conservation in a NaCl-H2O fluid to model the seafloor hydrothermal processes. These simulations use homogeneous box geometries at a fixed seafloor pressure of 25 MPa with constant bottom temperature boundary conditions that represent a sub-axial magma chamber to explore the effects of permeability, maximum bottom temperature and system depth on the evolution of vent fluid temperature and salinity, and heat output. We also study the temporal and spatial variability in hydrothermal circulation. The two-phase simulation results show that permeability plays an important role in plume structure and heat output of hydrothermal systems, but it has little effect on vent fluid temperature and salinity, given the same bottom temperature. For some permeability values, multiple plumes can vent at the seafloor above the simulated magma chamber. Temporal variability of vent fluid temperature and salinity and the complexity of phase separation suggest that pressure and temperature conditions at the top of the axial magma chamber cannot be easily inferred from vent fluid temperature and salinity alone. Vapor and brine derived fluids can vent at the seafloor simultaneously, even from neighboring locations that are fed by the same plume. / Master of Science

Ever since the discovery of seafloor hydrothermal systems at mid ocean ridges, scientists have been trying to understand the complex dynamic processes by which thermal energy is transported advectively by chemically reactive aqueous fluids from Earth's interior to the surface. Hydrothermal systems are generally assumed to consist of a heat source and a fluid circulation system. Understanding the interconnected physical, chemical, biological, and geological processes at oceanic spreading centers is important because these processes affect the global energy and biogeochemical budgets of the Earth system. Despite two decades of focused study of hydrothermal systems, several key questions remain concerning the behavior and evolution of hydrothermal vent systems. Among these are: (a) the partitioning of heat transport between focused and diffuse flow, and (b) the spatial extent and distribution of hydrothermal recharge. These are the main topics of investigation in this thesis. To address these issues, I first use a single-pass modeling approach using a variety of observational data in a simple parametric scale analysis of a hydrothermal vent field to determine fundamental parameters associated with the circulation and magmatic heat transfer for a number of seafloor hydrothermal systems for which the constraining data are available. To investigate the partitioning of heat flux between focused high temperature and diffuse flow I extend the one-limb single pass model to incorporate two single-pass limbs to represent deep and shallow circulation pathways. As a result, I find that 90% of the heat output is from high temperature fluid circulating in the deep limb even though much of the heat loss appears at the seafloor as low-temperature diffuse flow. Next, I use the parametric description of hydrothermal circulation to investigate hydrothermal recharge at the East Pacific Rise 9°50′ N hydrothermal site. Using a 1-D model of recharge through an area of 10⁵ m² elucidated by microseismicity in the oceanic crust I find that anhydrite precipitation is likely to result in rapid sealing of pore space in the recharge zone. This would lead to rapid decay of hydrothermal venting, which is contrary to observations. Then I consider two-dimensional numerical models of hydrothermal circulation in a porous box heated from below. The preliminary results of these models suggests that the anhydrite precipitation zone will be more diffuse, but additional work is needed to test whether anhydrite precipitation will seal the pore space. / Master of Science

My dissertation research consists of two themes: (a) the analysis of Middle Jurassic - Early Cretaceous marine magnetic anomalies (M-anomalies) in order to construct a comprehensive geomagnetic polarity timescale and (b) the investigation of the upper oceanic crustal architecture using downhole geophysical logs. These themes were chosen to better understand how remotely-sensed geophysical signals elucidate the formation and evolution of oceanic crust. This revised Pacific-wide MGPTS model shows significant improvement in its reliability, exhibits global applicability, and highlights changes in the paleo-Pacific spreading regime. By integrating Atlantic Manomaly analyses with the new MGPTS model and reviewing previous seismic studies, we shed new light on the causes of a ubiquitously distributed ?Atlantic anomaly smooth zone? where little coherency among M5-M15 anomaly sequence is observed. For the second theme, I analyzed the architecture of 15 m.y. old superfast spreading East Pacific Rise crust drilled at Ocean Drilling Program Hole 1256D in the eastern Pacific. An intact upper oceanic crustal section was penetrated at this site to a depth of 1507 mbsf. In situ crustal architecture was mapped from resistivity imagery (electrofacies by Formation MicroScanner) combined with recovered cores and other logs. Highlights of this research are: (1) most of the extrusive section consists of massive flows and fragmented formations including breccias, which has important implications for the magnetic source layer and pathways of hydrothermal alteration; (2) the dike complex is composed of sheeted-dikes dipping away from the paleo-spreading axis consistent with submersible observations at other sites in the eastern Pacific; (3) the crustal construction processess from ridge axis to abyssal plain during 0-50 kyr time are consistent with previous seismic reflection studies based on the integration of our stratigraphy model with lava flow observations from the southern East Pacific Rise.

Mid-ocean ridge hydrothermal systems play an important role in the cycling of energy and mass between the solid earth and oceans. The on-axis low-temperature diffuse fluids (temperature lower than ~100 °C) carry ~90% of the on-axis heat fluxes, but diffuse fluids generation is poorly constrained. This study uses the abundance of the chalcophile elements, which form metal-sulphides in the rock record, to test models for diffuse fluids generation. These include mixing between seawater and high-temperature hydrothermal fluids and conductive cooling of high- temperature hydrothermal fluids. This thesis determined the concentrations of the elements of interest (As, Mo, Ag, Cd, Sn, and Tl) in geological reference materials using standard addition method in ICP-MS. These values were used to calibrate the analysis of samples from Pito Deep to trace the abundance of these elements within the upper oceanic crust. The results show that the Zn, Cu, As, Ag, Cd, Tl, and Pb are generally depleted in sheeted dikes and enriched in the lava unit and/or the transition zone, which is consistent with previous studies on fast-spreading EPR crust at Hole 504B, Hess Deep and Hole 1256D. The enrichment of these elements in the lava unit and/or the transition zone suggests that cooling high-temperature hydrothermal fluids to form diffuse fluids occurred in this iii iv area of the oceanic crust. Molybdenum and Sb are added into all units of the crust by recharging seawater. The concentrations of chalcophile elements in diffuse fluids were calculated by a mass balance. The results of this study favored a diffuse fluids generation model that involves mixing of seawater and high-temperature hydrothermal fluids. Results also show that the observed concentrations of Mo and Sb requires extra input source besides recharging seawater and oceanic crust, possibly particulates in seawater. / Graduate / juliatian2013@gmail.com