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

Feedbacks of sedimentation on crustal heat flow: new insights from the Vøring Basin, Norwegian Sea

Basement heat flow is one of the key unknowns in sedimentary basin analysis. Its quantification is challenging not in the least due to the various feedback mechanisms between the basin and lithosphere processes. This study explores two main feedbacks, sediment blanketing and thinning of sediments during lithospheric stretching, in a series of synthetic models and a reconstruction case study from the Norwegian Sea. Three types of basin models are used: (1) a newly developed one-dimensional (1D) forward model, (2) a decompaction/backstripping approach and (3) the commercial basin modelling software TECMOD2D for automated forward basin reconstructions. The blanketing effect of sedimentation is reviewed and systematically studied in a suite of 1D model runs. We find that even for moderate sedimentation rates (0.5 mm year−1), basement heat flow is depressed by ∼25% with respect to the case without sedimentation; for high sedimentation rates (1.5 mm year−1), basement heat flow is depressed by ∼50%. We have further compared different methods for computing sedimentation rates from the presently observed stratigraphy. Here, we find that decompaction/backstripping-based methods may systematically underestimate sedimentation rates and total subsidence. The reason for this is that sediments are thinned during lithosphere extension in forward basin models while there are not in backstripping/decompaction approaches. The importance of sediment blanketing and differences in modelling approaches is illustrated in a reconstruction case study from the Norwegian Sea. The thermal and structural evolution of a transect across the Vøring Basin has been reconstructed using the backstripping/decompaction approach and TECMOD2D. Computed total subsidence curves differ by up to ∼3 km and differences in computed basement heat flows reach up to 50%. These findings show that strong feedbacks exist between basin and lithosphere processes and that resolving them require integrated lithosphere-scale basin models

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