Satish C. Singh

Discovery of a magma chamber and faults beneath a Mid-Atlantic Ridge hydrothermal field

Crust at slow-spreading ridges is formed by a combination of magmatic and tectonic processes, with magmatic accretion possibly involving short-lived crustal magma chambers. The reflections of seismic waves from crustal magma chambers have been observed beneath intermediate and fast-spreading centres, but it has been difficult to image such magma chambers beneath slow-spreading centres, owing to rough seafloor topography and associated seafloor scattering. In the absence of any images of magma chambers or of subsurface near-axis faults, it has been difficult to characterize the interplay of magmatic and tectonic processes in crustal accretion and hydrothermal circulation at slow-spreading ridges. Here we report the presence of a crustal magma chamber beneath the slow-spreading Lucky Strike segment of the Mid-Atlantic Ridge. The reflection from the top of the magma chamber, centred beneath the Lucky Strike volcano and hydrothermal field, is approximately 3 km beneath the sea floor, 3–4 km wide and extends up to 7 km along-axis. We suggest that this magma chamber provides the heat for the active hydrothermal vent field above it. We also observe axial valley bounding faults that seem to penetrate down to the magma chamber depth as well as a set of inward-dipping faults cutting through the volcanic edifice, suggesting continuous interactions between tectonic and magmatic processes.

Melt to mush variations in crustal magma properties along the ridge crest at the southern East Pacific Rise

The determination of along-axis variations in melt properties within the crustal axial magma chamber beneath fast spreading axes is important for understanding melt delivery from the mantle, eruption history along the ridge crest, and the process of crustal accretion. Seismic reflection images have shown the molten sill to be continuous along the ridge crest for many tens of kilometres with varying widths (250–4,500 m), but variations in its seismic properties and thickness have remained elusive, despite several attempts to constrain these properties. Here we report that the melt sill along the southern East Pacific Rise, which is about 50 m thick, undergoes abrupt changes in its internal properties, ranging from pure melt to mush. The 60-km-long ridge-crest segment near 14° 00′ S is characterized by three 2–4-km sections containing pure melt embedded within a magma chamber rich in mush. These small pure melt pockets may represent a fresh supply of magma from the mantle, capable of erupting and forming the upper crust. Conversely, the 80–90% of the magma chamber which is mushy is unlikely to erupt and may influence the lower-crustal accretion.

Elastic properties of hydrate‐bearing sediments using effective medium theory

Accurate and detailed models of the seismic velocity structure of gas hydrate-bearing sediments may be determined by careful analysis of controlled source seismic data. However, interpretation of these velocities in terms of hydrate saturation of the pore space has hitherto relied on semiempirical formulas and/or simple effective medium theory. We develop a rigorous theoretical scheme to relate the seismic properties of a clay-rich hydrate-bearing sediment to its porosity, mineralogy, microstructural features and hydrate saturation. We consider separately the two possible end-members for the distribution of hydrate in the pore space: (1) hydrates are unconnected and located in the pore voids without appreciable grain contact and (2) connected hydrates are forming cement binding around the grains. The scheme is transversely isotropic, to allow for anisotropy due to alignment of clay platelets, and is based on a combination of a self-consistent approximation, a differential effective medium theory, and a method of smoothing for crystalline aggregates. We have applied the scheme to lithological and seismic velocity data from Ocean Drilling Program Site 995 on the Blake Ridge (southeastern U.S. continental margin) to make estimates of the hydrate saturation. It was found that the hydrates are probably unconnected, and their volume concentration varies between approx. 0% at 100 m below the seabed and approx. 9% at 400 m depth, just above the “bottom simulating reflector”, if the clay platelet orientation distribution resembles the function we have used.

Seismic velocity structure at a gas hydrate reflector, offshore western Colombia, from full waveform inversion

Seismic reflection profiles across many continental margins have imaged bottom simulating reflectors (BSRs), which have been interpreted as being formed at the base of a methane hydrate stability field. Such reflectors might arise either from an impedance contrast between high-velocity, partially hydrated sediments and water-saturated sediments or from a contrast with partially gas-saturated sediments. These alternatives may be hard to distinguish by conventional amplitude-versus-offset or waveform modeling approaches. Here we investigate the origin of a high amplitude BSR in the accretionary wedge offshore of western Colombia by seismic waveform inversion. The inversion procedure consists of three steps: firstly, determination of root-mean-square velocities and hence estimates of the interval velocities between major reflectors by a global grid search for maximum normalized energy along elliptical trajectories in the intercept time-slowness domain; secondly, determination of accurate interval velocities between these reflectors by a Monte Carlo search for maximum energy; and thirdly, a waveform fit in the frequency-slowness domain, using differential reflectivity seismograms and a conjugate-gradient optimization algorithm to minimize the sample-by-sample waveform misfit between data and synthetic. At two locations, near a structural high, we find a ∼30-m thick low-velocity zone beneath the BSR, with the properties of a partially gas-saturated zone, while at a third location, where the BSR amplitude is lower, we find no evidence for anomalously low velocities. The preferential development of the BSR in structures that would tend to intercept fluid flow or migrating gas and the presence of free gas beneath the BSR indicate a mechanism of BSR formation in which free methane gas migrates upward into the hydrate stability field or is carried there in advecting pore water.

Evidence from three-dimensional seismic reflectivity images for enhanced melt supply beneath mid-ocean-ridge discontinuities

Quantifying the melt distribution and crustal structure across ridge-axis discontinuities is essential for understanding the relationship between magmatic, tectonic and petrologic segmentation of mid-ocean-ridge spreading centres. The geometry and continuity of magma bodies beneath features such as overlapping spreading centres can strongly influence the composition of erupted lavas and may give insight into the underlying pattern of mantle flow. Here we present three-dimensional images of seismic reflectivity beneath a mid-ocean ridge to investigate the nature of melt distribution across a ridge-axis discontinuity. Reflectivity slices through the 9° 03′ N overlapping spreading centre on East Pacific Rise suggest that it has a robust magma supply, with melt bodies underlying both limbs and ponding of melt beneath large areas of the overlap basin. The geometry of melt distribution beneath this offset is inconsistent with large-scale, crustal redistribution of melt away from centres of upwelling. The complex distribution of melt seems instead to be caused by a combination of vertical melt transport from the underlying mantle and subsequent focusing of melt beneath a magma freezing boundary in the mid-crust.

Seismic velocity studies of a gas hydrate bottom‐simulating reflector on the northern Cascadia continental margin: Amplitude modeling and full waveform inversion

On the northern Cascadia subduction margin, the multichannel seismic amplitude-versus-offset (AVO) behavior of a bottom-simulating reflector (BSR) suggests a P wave velocity change from high-velocity hydrate-bearing sediment to lower velocity sediment containing a small amount of free gas. The observed nonlinear AVO behavior, constant or slightly decreasing amplitudes at near-to-mid offsets and a large-amplitude increase at far offsets, can be reproduced in the models if an S wave velocity enhancement is assumed as expected from hydrate cementation. The AVO behavior of the hydrate BSR is found not to be as useful as was earlier thought for determining the amount of free gas below the BSR. This is because Poissons ratio change below the BSR due to gas is likely very small in high-porosity unconsolidated sediments. The uncertainty in the models is large, as there is no reliable S wave velocity information for the sediments containing hydrate or gas. AVO modeling alone is not sufficient to distinguish different velocity models across the BSR. Our interpretation of the BSR amplitude behavior is that a P wave velocity increase above the BSR is the main cause of the BSR reflection amplitude increase at large incidence angles. Caution must be taken in applying AVO analysis, as little is known about S wave velocities. However, subtle differences in BSR amplitude behavior and reflection waveform can provide constraints through very careful full waveform inversion. A well-defined reference velocity-depth profile is also required to represent water-saturated sediment unaffected by either hydrate or free gas. Using full waveform velocity inversion, a high-resolution velocity model for the hydrate BSR has been derived. The best fit model for the seismic data near the Ocean Drilling Program (ODP) sites 889/890 consists of a high-velocity zone above the BSR and a thin low-velocity layer below, in agreement with the ODP downhole velocity data.

Cooling of the lower oceanic crust

Thermal models of mid-ocean ridges that balance the influx of heat from magmatic sources with the removal of heat by conduction and hydrothermal circulation allow quantification of cooling of young oceanic crust. These models reproduce key observations relating to crustal accretion and hydrothermal cooling at fast-spreading ridges. The rate of cooling is constrained both by the bathymetry of ridge axes and by olivine compositions from ophiolite gabbros. Successful models involve extensive hydrothermal cooling of the lower crust within 20 km of the ridge, with ~50–70 kW of hydrothermal cooling for every 1 m of ridge axis at crustal ages of <0.1–0.4 Ma. These timates can be used to refine global models of geochemical and thermal fluxes close to spreading ridges.

Seismic evidence for a hydrothermal layer above the solid roof of the axial magma chamber at the southern East Pacific Rise

A full-waveform inversion of two-ship, wide-aperture, seismic reflection data from a ridge-crest seismic line at the southern East Pacific Rise indicates that the axial magma chamber here is about 50 m thick, is embedded within a solid roof, and has a solid floor. The 50--60-m-thick roof is overlain by a 150--200-m-thick low-velocity zone that may correspond to a fracture zone that hosts the hydrothermal circulation, and the roof itself may be the transition zone separating the magma chamber from circulating fluids. Furthermore, enhanced hydrothermal activity at the sea floor seems to be associated with a fresh supply of magma in the crust from the mantle. The presence of the solid floor indicates that at least the upper gabbros of the oceanic lower crust are formed by cooling and crystallization of melt in magma chambers.

Velocity structure of a bottom simulating reflector offshore Peru: Results from full waveform inversion

Abstract Much of our knowledge of the worldwide distribution of submarine gas hydrates comes from seismic observations of Bottom Simulating Reflectors (BSRs). Full waveform inversion has proven to be a reliable technique for studying the fine structure of BSRs using the compressional wave velocity. We applied a non-linear full waveform inversion technique to a BSR at a location offshore Peru. We first determined the large-scale features of seismic velocity variations using a statistical inversion technique to maximise coherent energy along travel-time curves. These velocities were used for a starting velocity model for the full waveform inversion, which yielded a detailed velocity/depth model in the vicinity of the BSR. We found that the data are best fit by a model in which the BSR consists of a thin, low-velocity layer. The compressional wave velocity drops from 2.15 km/s down to an average of 1.70 km/s in an 18 m thick interval, with a minimum velocity of 1.62 km/s in a 6 m interval. The resulting compressional wave velocity was used to estimate gas content in the sediments. Our results suggest that the low velocity layer is a 6–18 m thick zone containing a few percent of free gas in the pore space. The presence of the BSR coincides with a region of vertical uplift. Therefore, we suggest that gas at this BSR is formed by a dissociation of hydrates at the base of the hydrate stability zone due to uplift and subsequently a decrease in pressure.

Seismic reflection images of the Moho underlying melt sills at the East Pacific Rise

The determination of melt distribution in the crust and the nature of the crust–mantle boundary (the ‘Moho’) is fundamental to the understanding of crustal accretion processes at oceanic spreading centres. Upper-crustal magma chambers have been imaged beneath fast- and intermediate-spreading centres but it has been difficult to image structures beneath these magma sills. Using three-dimensional seismic reflection images, here we report the presence of Moho reflections beneath a crustal magma chamber at the 9° 03′ N overlapping spreading centre, East Pacific Rise. Our observations highlight the formation of the Moho at zero-aged crust. Over a distance of less than 7 km along the ridge crest, a rapid increase in two-way travel time of seismic waves between the magma chamber and Moho reflections is observed, which we suggest is due to a melt anomaly in the lower crust. The amplitude versus offset variation of reflections from the magma chamber shows a coincident region of higher melt fraction overlying this anomalous region, supporting the conclusion of additional melt at depth.