Modes of crustal accretion and their implications for hydrothermal circulation

Hydrothermal convection at mid-ocean ridges links the ocean's long-term chemical evolution to solid earth processes, forms hydrothermal ore deposits, and sustains the unique chemosynthetic vent fauna. Yet the depth extent of hydrothermal cooling and the inseparably connected question of how the lower crust accretes remain poorly constrained. Here based on coupled models of crustal accretion and hydrothermal circulation, we provide new insights into which modes of lower crust formation and hydrothermal cooling are thermally viable and most consistent with observations at fast-spreading ridges. We integrate numerical models with observations of melt lens depth, thermal structure, and melt fraction. Models matching all these observations always require a deep crustal-scale hydrothermal flow component and less than 50% of the lower crust crystallizing in situ

Faulting and off-axis submarine massive sulfide accumulation at slow spreading mid-ocean ridges: A numerical modeling perspective

The potential of mining seafloor massive sulfide deposits for metals such as Cu, Zn, and Au is currently debated. One key challenge is to predict where the largest deposits worth mining might form, which in turn requires understanding the pattern of subseafloor hydrothermal mass and energy transport. Numerical models of heat and fluid flow are applied to illustrate the important role of fault zone properties (permeability and width) in controlling mass accumulation at hydrothermal vents at slow spreading ridges. We combine modeled mass-flow rates, vent temperatures, and vent field dimensions with the known fluid chemistry at the fault-controlled Logatchev 1 hydrothermal field of the Mid-Atlantic Ridge. We predict that the 135 kilotons of SMS at this site (estimated by other studies) can have accumulated with a minimum depositional efficiency of 5% in the known duration of hydrothermal venting (58,200 year age of the deposit). In general, the most productive faults must provide an efficient fluid pathway while at the same time limit cooling due to mixing with entrained cold seawater. This balance is best met by faults that are just wide and permeable enough to control a hydrothermal plume rising through the oceanic crust. Model runs with increased basal heat input, mimicking a heat flow contribution from along-axis, lead to higher mass fluxes and vent temperatures, capable of significantly higher SMS accumulation rates. Nonsteady state conditions, such as the influence of a cooling magmatic intrusion beneath the fault zone, also can temporarily increase the mass flux while sustaining high vent temperatures

Tectonic evolution and extension at the Møre Margin – Offshore mid-Norway

Highlights • New and reprocessed seismic data improved structural mapping at the Møre Margin. • Time-structure and thickness maps of the Cretaceous units have been constructed. • Stratigraphy reconstruction of a transect reveals 188 km extension. • Average stretching factor is 2.2–3.6 depending on assumed initial crustal thickness. Abstract Lithospheric stretching is the key process in forming extensional sedimentary basins at passive rifted margins. This study explores the stretching factors, resulting extension, and structural evolution of the Møre segment on the Mid-Norwegian continental margin. Based on the interpretation of new and reprocessed high-quality seismic, we present updated structural maps of the Møre margin that show very thick post-rift sediments in the central Møre Basin and extensive sill intrusion into the Cretaceous sediments. A major shift in subsidence and deposition occurred during mid-Cretaceous. One transect across the Møre continental margin from the Slørebotn Subbasin to the continent-ocean boundary is reconstructed using the basin modelling software TecMod. We test different initial crustal configurations and rifting events and compare our structural reconstruction results to stretching factors derived both from crustal thinning and the classical backstripping/decompaction approach. Seismic interpretation in combination with structural reconstruction modelling does not support the lower crustal bodies as exhumed and serpentinised mantle. Our extension estimate along this transect is ~ 188 ± 28 km for initial crustal thickness varying between 30 and 40 km

Pre-breakup magmatism on the Vøring margin: Insight from new sub-basalt imaging and results from Ocean Drilling program hole 642E

Highlights • Sub-basalt imaging improvement on the Vøring Margin • Definition of a new seismic facies unit: the Lower Series Flows • Significant organic carbon content within the melting crustal segment • Apectodinium augustum marker for the PETM is reworked into the Lower Series Flows • The Lower Series Flows, early Eocene in age, predate the Vøring Margin breakup Abstract Improvements in sub-basalt imaging combined with petrological and geochemical observations from the Ocean Drilling Program (ODP) Hole 642E core provide new constraints on the initial breakup processes at the Vøring Margin. New and reprocessed high quality seismic data allow us to identify a new seismic facies unit which we define as the Lower Series Flows. This facies unit is seismically characterized by wavy to continuous subparallel reflections with an internal disrupted and hummocky shape. Drilled lithologies, which we correlate to this facies unit, have been interpreted as subaqueous flows extruding and intruding into wet sediments. Locally, the top boundary of this facies unit is defined as a negative in polarity reflection, and referred as the K-Reflection. This reflection can be correlated with the spatial extent of pyroclastic deposits, emplaced during transitional shallow marine to subaerial volcanic activities during the rift to drift transition. The drilled Lower Series Flows consist of peraluminous, cordierite bearing peperitic basaltic andesitic to dacitic flows interbedded with thick volcano-sedimentary deposits and intruded sills. The peraluminous geochemistry combined with available C (from calcite which fills vesicles and fractures), Sr, Nd, and Pb isotopes data point towards upper crustal rock-mantle magma interactions with a significant contribution of organic carbon rich pelagic sedimentary material during crustal anatexis. From biostratigraphic analyses, Apectodinium augustum was found in the The Lower Series Flows. This species is a marker for the Paleocene – Eocene Thermal Maximum (PETM). However, the absence of very low carbon isotope values (from bulk organic matter), that characterize the PETM, imply that A.augustum was reworked into the early Eocene sediments of this facies unit which predate the breakup time of the Vøring Margin. Finally, a plausible conceptual emplacement model for the Lower Series Flows facies unit is proposed. This model comprises several stages: (1) the emplacement of subaqueous peperitic basaltic andesitic flows intruding and/or extruding wet sediments; (2) a subaerial to shallow marine volcanism and extrusion of dacitic flows; (3) a proto-breakup phase with intense shallow marine to subaerial explosive volcanism responsible for pyroclastic flow deposits which can be correlated with the seismic K-Reflection and (4) the main breakup stage with intense transitional tholeiitic MORB-type volcanism and large subsidence concomitant with the buildup of the Seaward Dipping Reflector wedge

Coupled mechanical and hydrothermal modeling of crustal accretion at intermediate to fast spreading ridges

The genesis of oceanic crust at intermediate to fast spreading ridges occurs by the crystallization of mantle melts accumulated in at least one shallow melt lens situated below the ridge axis. Seismic reflection data suggest that the depth of this melt lens is inversely correlated with spreading rate and thereby magma supply. The heat released in it by crystallization and melt injection is removed by a combination of hydrothermal cooling and diffusion. Due to the different time scales of hydrothermal cooling and crustal accretion, numerical models have so far focused on only one of the two processes. Here we present the results from a coupled model that solves simultaneously for crustal accretion and hydrothermal cooling. Our approach resolves both processes within one 2D finite-element model that self-consistently solves for crustal, mantle, and hydrothermal flow. The formation of new oceanic crust is approximated as a gabbro glacier, in which the entire lower crust crystallizes in one shallow melt lens. We find that the depth of the melt lens and the shape of hot (potentially molten) lower crust are highly dependent on the ridge permeability structure. The predicted depth of the melt lens is primarily controlled by the permeability at the ridge axis, whereas the off-axis permeability determines the width of hot lower crust. A detailed comparison of the modeling results with observed locations of the melt lens at intermediate to fast spreading ridges shows that only a relatively narrow range of crustal permeabilities is consistent with observations. In addition, we find significant deviations between models that resolve or parameterize hydrothermal cooling: the predicted crustal thermal structures show major differences for models that predict the same melt lens location. This illustrates the importance of resolving hydrothermal flow in simulations of crustal accretion

Hybrid shallow on-axis and deep off-axis hydrothermal circulation at fast-spreading ridges

Hydrothermal flow at oceanic spreading centres accounts for about ten per cent of all heat flux in the oceans and controls the thermal structure of young oceanic plates. It also influences ocean and crustal chemistry, provides a basis for chemosynthetic ecosystems, and has formed massive sulphide ore deposits throughout Earth’s history. Despite this, how and under what conditions heat is extracted, in particular from the lower crust, remains largely unclear. Here we present high-resolution, whole-crust, two- and three-dimensional simulations of hydrothermal flow beneath fast-spreading ridges that predict the existence of two interacting flow components, controlled by different physical mechanisms, that merge above the melt lens to feed ridge-centred vent sites. Shallow on-axis flow structures develop owing to the thermodynamic properties of water, whereas deeper off-axis flow is strongly shaped by crustal permeability, particularly the brittle–ductile transition. About 60 per cent of the discharging fluid mass is replenished on-axis by warm (up to 300 degrees Celsius) recharge flow surrounding the hot thermal plumes, and the remaining 40 per cent or so occurs as colder and broader recharge up to several kilometres away from the axis that feeds hot (500–700 degrees Celsius) deep-rooted off-axis flow towards the ridge. Despite its lower contribution to the total mass flux, this deep off-axis flow carries about 70 per cent of the thermal energy released at the ridge axis. This combination of two flow components explains the seismically determined thermal structure of the crust and reconciles previously incompatible models favouring either shallower on-axis or deeper off-axis hydrothermal circulation